跃响应的Matlab仿真。um=1;xita=0.6;%阻屁比xita=0.6Wn1=1;den1=[1/Wn1^22*xita/Wn11];%Wn1=1时,闭环传递函数分母多项式系数Wn2=1.5;den2=[1/Wn2^22*xita/Wn21];%Wn2=1.5时,闭环传递函数分母多项式系数Wn3=2.5;den3=[1/Wn3^22*xita/Wn31];Wn4=4;den4=[1/Wn4^22*xita/Wn41];Wn5=6;den5=[1/Wn5^22*xita/Wn51];Wn6=12;den6=[1/Wn6^22*xita/Wn61];t=0:0.05:6;[y1xt]=step(num,den1,t);%step是单位阶跃命令,等号前面的[]是用来把单位阶跃响应的输入时间和对应的输出值返回成数值矩阵。step前面加了等号,不能直接绘单位阶跃响应图。[y2xt]=step(num,den2,t);[y3xt]=step(num,den3,t);[y4xt]=step(num,den4,t);[y5xt]=step(num,den5,t);[y6xt]=step(num,den6,t);figure(1)%生成图框,准备画图plot(t,y1,'b-..',t,y2,'r--',t,y3,'k +',t,y4,'g:',t,y5,'m-.',t,y6,'b');%plot是绘图命令,把前面生成的响应数值绘成图。单引号里边的字母表示图形的颜色,符号表示的线型。默认颜色是黑色,默认的线型是实线。xlabel('t(s)'),ylabel('y(t)');title('固有频率Wn=1,1.5,2.5,4,6,12')legend('Wn1=1','Wn2=1.5','Wn3=2.5','Wn4=4','Wn5=6','Wn6=12')%图形的线型、颜色附加说明。
全文预览
