I'm trying to look at the phase and gain margin of an openloop system. Due to presentation/publications, I had to designa function to plot my amplitude and phase response with good font, font dimensions etc. To do so, I'm using the bode function to extract the data from my state space systems. But I have a problem with a phase offset (value at 0Hz). I know we can remove it from the properties editor avalaible in the classical bode plot but is there a solution to correct it when using this syntax [mag,phase,wout] = bode(sys,w) ?