bode() magnitude

50 views (last 30 days)
aaa
aaa on 2 May 2012
Hi, I Have seen the code [ mag,phase] = bode(sy s) but when i try to see the values of mag, or phase, it returns something like mag(:,:,1). What is the significance of the stuff in the brackets? How can i just get the magnitude values with respect to frequency so I can plot a 2d plot of the magnitude and frequency?
thanks

Accepted Answer

Arkadiy Turevskiy
Arkadiy Turevskiy on 2 May 2012
Magnitude and phase outputs of bode are 3-D arrays, for representing multi-input multi-output systems. The first two dimensions correspond to the number of inputs and outputs, the third dimension corresponds to the number of frequency points. In your case you have a single-input single-output system, so size(mag) is 1 1 n, where n is the number of frequency points.
To get a standard bode plot you can just do
bode(sys)
If you want access to mag and phase as regular vectors, just do
mag=squeeze(mag)
All of this is covered in the doc.
Arkadiy

More Answers (4)

aaa
aaa on 2 May 2012
hey thanks for your answer.
i have a followup question. what i did after. so i have
[mag,phase,wout] = bode(sys);
mag = squeeze(mag);
phase = squeeze(phase);
wout = squeeze(wout);
semilogx(wout,mag)
semilogx(wout,phase)
so when i plot the phase, i get the same result as i would doing bode(sys). The problem is when i plot the magnitude, the values on the y-axis are incorrect. I then added
mag = 20*log10(mag)
semilogx(wout,mag)
but it still appears to be incorrect. Can you provide some insight on maybe where i am going wrong?

Arkadiy Turevskiy
Arkadiy Turevskiy on 2 May 2012
semilogx(wout,20*log10(squeeze(mag)))
gives you the same plot for magnitude as bode.
If it does not look the same, it is because the axes limits do not line up. You can customize that and other plot parameters to your liking as described here.

aaa
aaa on 2 May 2012
i think i found out the problem. wout must be divided by 2pi to get the correct plot.

aaa
aaa on 2 May 2012
SOrry as another question, I noticed that when you do
[mag, phase, wout] = bode(a)
it draws only about 45 points. Is there a way to get more?
Thanks
  1 Comment
Arkadiy Turevskiy
Arkadiy Turevskiy on 3 May 2012
win=logspace(wout(1),wout(45),100);
% replace 100 with desired number of points
[mag,phase] = bode(sys,win);

Sign in to comment.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!