bode() magnitude

22 vues (au cours des 30 derniers jours)
aaa le 2 Mai 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
0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

Réponse acceptée

Arkadiy Turevskiy le 2 Mai 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.
0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

Plus de réponses (4)

aaa le 2 Mai 2012
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 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

Arkadiy Turevskiy le 2 Mai 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.
0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

aaa le 2 Mai 2012
i think i found out the problem. wout must be divided by 2pi to get the correct plot.
0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

aaa le 2 Mai 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 commentaireAfficher -1 commentaires plus anciensMasquer -1 commentaires plus anciens
Arkadiy Turevskiy le 3 Mai 2012
win=logspace(wout(1),wout(45),100);
% replace 100 with desired number of points
[mag,phase] = bode(sys,win);

Connectez-vous pour commenter.

Catégories

En savoir plus sur Plot Customization dans Help Center et File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by