unwrap problems

[Ron Crummett]

Sorry about this repeated post, I posted on this a while ago and forgot 
to do anything about it.

I am having a problem using unwrap.  I am trying to plot the bode 
diagram and get phase wraparound:

H = tf([330 49.5], [1 20.015 100.3 1.5 0]);
w = logspace(-3, 3, 601);
[mag, ph] = bode(H, w);
semilogx(w, ph)

The plotting command plots a wrapped phase, it jumps from -180 to 180.  
I have tried using the unwrap command but it doesn't seem to do much; 
however, on further review it looks like it unwraps the phase in terms 
of radians, whereas the bode function returns the phase in degrees.

-Wouldn't it be easier if the bode function automatically called the 
unwrap function?  I have edited my bode.m script to do this and now 
everything works fine; is there a reason it hasn't been done?

-Ron



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------