Why I get the wrong peak from bode plot?
61 views (last 30 days)
Show older comments
Bob on 20 Mar 2016
m1 = 2500; % (kg)
m2 = 320; % (kg)
k1 = 80000; % (N/m)
k2 = 500000; % (N/m)
b1 = 350; % (N*s/m)
b2 = 15020; % (N*s/m)
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
[gpeak,fpeak] = getPeakGain(G1)
My resault is wrong according to the zoomed diagram.
According to the picture below I should get:
gpeak = 25 ~ 27
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.
gpeak_dB = 20*log10(gpeak);
and check the value. It should be approx 26.4, which is consistent with the Bode plot.
More Answers (0)
Find more on Time and Frequency Domain Analysis in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!