Analysis Of Engineering Systems - Homework 7

David Clark MAE 488 Homework #7

9.2) Input %9.2) sys1=tf([15],[6 2]) sys2=tf([9 0],[8 4]) sys3=tf([6*14 6*7],[10 2]) bode(sys1) bode(sys2) bode(sys3) Output

David Clark MAE 488 Homework #7

David Clark MAE 488 Homework #7

9.13) Input program: %9.13) roots( conv([1 0],[1 6 58]) ); phi1=cos(atan(7/3)); wn1=sqrt(3^2+7^2); %The resonant frequencies are.. wr1a=0 wr1b=wn1*sqrt(1-2*phi1^2) %Part B roots( conv([3 18 174],[2 8 58]) ); phi2a=cos(atan(7/3)); wn2a=sqrt(3^2+7^2); phi2b=cos(atan(5/2)); wn2b=sqrt(2^2+5^2); %The resonant frequencies are.. wr1=wn2a*sqrt(1-2*phi2a^2) wr2=wn2b*sqrt(1-2*phi2b^2) Output wr1a = 0 wr1b = 6.3246 wr1 = 6.3246 wr2 = 4.5826

9.43) Input Program: sys1=tf([4096], [1, 192, 4096]) bode(sys1)

David Clark MAE 488 Homework #7

Output

9.48) Input Program: m1=250; m2=40; c1=1000; k1=1.5e4; k2=1.5e5; c2=0; sys1=tf( [k2*c1, k2*k1], conv([m1, c1, k1],[m2, c1, k1+k2])+conv([0, c1, k1],[0, c1, k1]) ) sys2=tf( [k2*m1, k2*c1, k1*k2], conv([m1, c1, k1],[m2, c1, k1+k2])+conv([0, c1, k1],[0, c1, k1]) ) bode(sys1) bode(sys2)

David Clark MAE 488 Homework #7

Output