bode() magnitude

aaa
2012 5 2
5 个回答
回答已采纳
18 次查看（30 天）

0 个投票

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

请先登录，再回答此问题。

 采纳的回答

1 个投票

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

更多回答（4 个）

aaa
aaa 2012-5-2

0 个投票

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?

0 个投票

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 2012-5-2

0 个投票

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

0 个投票

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 个评论
显示 -1更早的评论 隐藏 -1更早的评论

win=logspace(wout(1),wout(45),100);
% replace 100 with desired number of points
[mag,phase] = bode(sys,win);

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Plot Customization 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by