MATLAB程序.doc

求超调k=10;n1=conv(2.27 1,1.499 1);d1=conv(conv(conv(conv(1 0,1 1),0.25 1),0.11 1),11.992 1);t=0:0.001:50;G1=tf(k*n1,d1);sys=feedback(G1,1);y=step(sys,t);r=1;while y(r)0.98&y(s)1.02;s=s-1;end;Ts=(s-1)*0.001求Bodek0=10;n1=conv(2.27 1,1.499 1);d1=conv(conv(conv(conv(1 0,1 1),0.25 1),0.11 1),11.992 1);mag,phase,w=bode(k0*n1,d1);figure(1);margin(mag,phase,w);grid单位阶跃K0=10;n1=conv(2.27 1,1.499 1);d1=conv(conv(conv(conv(1 0,1 1),0.25 1),0.11 1),11.992 1);t=0:0.001:20;G=tf(K0*n1,d1);sys=feedback(G,1);figure(1);grid;step(sys,t);grid校正装置f1=2.27,1;f2=1.499,1;num=conv(f1,f2);f3=0.11,1;f4=11.992,1;den=conv(f3,f4);bode(num,den)grid校正前斜坡%G0(s) Unit-Ramp Response%num1=10;den1=0.25,1.25,1,10,0,0;t=0:0.1:20;y1,z1,t=step(num1,den1,t);plot(t,y1,-);grid校正后单位斜坡num2=10*conv(2.27,1,1.499,1);den2= conv(conv(conv(conv(1 0,1 1),0.25 1),0.11 1),11.992 1);t=0:0.1:20;y2,z2,t=step(num2,den2,t);plot(t,y2,-);grid未校正时的Bode图，程序如下：k0=10;n1=1;d1=conv(conv(1 0,1 1),0.25 1);mag,phase,w=bode(k0*n1,d1);figure(1);margin(mag,phase,w);grid校正前阶跃K0=10;d1=conv(conv(1 0,1 1),0.25 1);t=0:0.001:20;G=tf(K0,d1);sys=feedback(G,1);figure(1);grid;step(sys,t);grid校正装置的Bode图,程序f1=2.457,1;f2=1.30,1;num=conv(f1,f2);f3=0.102,1;f4=12.99,1;den=conv(f3,f4);bode(num,den)grid校正后的Bode图,程序如下k0=10;n1=conv(2.457 1,1.30 1);d1=conv(conv(conv(conv(1 0,1 1),0.25 1),0.102 1),12.99 1);mag,phase,w=bode(k0*n1,d1);figure(1);margin(mag,phase,w);grid校正后系统阶跃响应曲线的MATLAB程序K0=10;n1=conv(2.457 1,1.30 1)d1=conv(conv(conv(conv(1 0,1 1),0.25 1),0.102 1),12.99 1);t=0:0.001:20;G=tf(K0*n1,d1);sys=feedback(G,1);figure(1);grid;step(sys,t);grid校正后的计算程序如下k=10;n1=conv(2.457 1,1.30 1)d1=conv(conv(conv(conv(1 0,1 1),0.25 1),0.102 1),12.99 1);t=0:0.001:50;G1=tf(k*n1,d1);sys=feedback(G1,1);y=step(sys,t);r=1;while y(r)0.98&y(s)1.02;s=s-1;end;Ts=(s-1)*0.001未校正时在单位斜坡信号作用下的响应曲线%G0(s) Unit-Ramp Response%num1=10;den1=0.25,1.25,1,10,0,0;t=0:0.1:20;y1,z1,t=step(num1,den1,t);plot(t,y1,-);grid校正后在单位斜坡信号作用下num2=10*conv(2.457,1,1.30,1);den2=conv(c