自动控制18090.docVIP

syms s for zeta=[0:0.2:0.8,1:0.5:2] wn=0.4; wn=sym(num2str(wn)); zet=sym(num2str(zeta)); if zeta==0 figure(1) ezplot(ilaplace(wn^2/s/(s^2+2*2*wn*s+wn^2)),[0 80]); grid on title(\xi=0) elseif zeta==1 figure(2) ezplot(ilaplace(wn^2/s/(s+wn)^2),[0 80]); hold on; else figure(2) ezplot(ilaplace(wn^2/s/(s^2+2*zet*wn*s+wn^2)),[0 80]); hold on; end end end grid on; title(\xi:0.2,0.4,0.6,0.8,1.0,1.5,2.0) axis([0 80 0 2]) gtext(0.4) gtext(1.0) gtext(2.0) 实验3 num=1; den=[4 1]; G=tf(num,den); bode(G,r); 实验3——2 w=[0,logspace(-2,2,200)]; wn=0.8; for zeta=[0.1 0.5 1 1.5 2]; G=tf([wn*wn],[wn^