%%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
Best Answer