MATLAB Answers

Why I get the wrong peak from bode plot?

34 views (last 30 days)
Xaris
Xaris on 20 Mar 2016
Edited: Xaris on 21 Mar 2016
%%Parameters
m1 = 2500; % (kg)
m2 = 320; % (kg)
k1 = 80000; % (N/m)
k2 = 500000; % (N/m)
b1 = 350; % (N*s/m)
b2 = 15020; % (N*s/m)
%%Transfer Function
num1 = [(0) (-m1*b2) (-m1*k2) (0) (0)];
den1 = [(m1*m2) (m1*b1+m1*b2+m2*b1) (m1*k1+m1*k2+m2*k1+b1*b2) (b1*k2+k1*b2) (k1*k2)];
G1 = tf(num1,den1);
%%Resonant Frequency of Dynamic Model
bode(G1)
grid on;
[gpeak,fpeak] = getPeakGain(G1)
My resault is wrong according to the zoomed diagram.
gpeak =
20.9719
fpeak =
5.2516
According to the picture below I should get:
gpeak = 25 ~ 27
ZOOMED

  0 Comments

Sign in to comment.

Accepted Answer

Rick Rosson
Rick Rosson on 20 Mar 2016
The function getPeakGain returns the peak gain on an absolute scale, whereas the Bode plot displays the gain on a decibel scale.
Please try:
gpeak_dB = 20*log10(gpeak);
and check the value. It should be approx 26.4, which is consistent with the Bode plot.

  1 Comment

Xaris
Xaris on 21 Mar 2016
Thank you, it works.
I have 2 Questions, I would appreciate it if you could answer me.
1) Can I change the units from Magnitude (dB) to Amplitude (m)?
2) How can I get the 2 peaks from the picture below? Because these commands gives me the only the highest.
[gpeak,fpeak] = getPeakGain(G1)
gpeak_dB = 20*log10(gpeak);

Sign in to comment.

More Answers (0)

Sign in to answer this question.