实验四线性控制系统的时域响应分析..docxVIP

实验四 线性控制系统的时域响应分析一、实验目的1、熟悉MATLAB有关命令的用法；2、用MATLAB系统命令对给定系统进行时域分析；二、实验内容求连续系统的单位阶跃响应有关命令的用法：命令格式：[y,x,t]=step(num, den) [y,x,t]=step(num, den, t)[y,x,t]=step(A,B,C,D)[y,x,t]=step(A,B,C,D,iu)[y,x,t]=step(A,B,C,D,iu,t)练习4-1. 给定系统的传递函数如下：求该系统的阶跃响应曲线,记录超调量、上升时间、过渡过程时间。num=25; den=[1 4 25]; t=0:0.1:10; [y,x,t]=step(num,den,t); plot(t,y)Step（num,den）M=((max(y)-1)/1)*100; disp([最大超调量M= num2str(M) %])最大超调量M=25.3177%练习4-2. 已知系统的开环传递函数为：求出该系统在单位负反馈下的阶跃响应曲线，记录超调量、上升时间、过渡过程时间。num=20;den=[1 8 36 40 0];[num,den]=cloop(num,den); t=0:0.1:10;[y,x,t]=step(num,den,t);plot(t,y) M=((max(y)-1)/1)*100;disp([最大超调量M= num2str(M) %])最大超调量M=2.5546%step(num,den)练习4-3 已知系统的传递函数为：+——求系统的阶跃响应； s=tf(s);gs=(6.3233*(s+1.4235)*(s+1.4235))/(s*s*(s+1)*(s+5));num=[6.323 18 12.81];den=[1 6 5 0 0];[num,den]=cloop(num,den);t=0:0.1:10;[y,x,t]=step(num,den,t);plot(t,y) 阶跃响应曲线线型用“*”号表示；plot(t,y,’*’) 阶跃响应图应加上横坐标名、纵坐标名和标题名，并加上网格线。title(阶跃响应图);xlabel(t);ylabel(y);set(gca,xgrid,on);set(gca,ygrid,on);练习4-4 求T1、T2、T3系统的阶跃响应；将 T1、T2、T3系统的阶跃响应图画在同一窗口内；num1=2;den1=[1 2 2];num2=[4 2];den2=[1 2 2];num3=1,den3=[2 3 3 1];t=0:0.1:10;[y1,x1,t1]=step(num1,den1,t1);subplot(1,3,1);plot(t1,y1);hold on;[y2,x2,t2]=step(num2,den2,t1);subplot(1,3,2);plot(t2,y2);hold on;[y3,x3,t3]=step(num3,den3,t1);subplot(1,3,3);plot(t3,y3) T1、T2、T3系统的阶跃响应曲线分别用不同的线形和颜色表示；num1=2;den1=[1 2 2];num2=[4 2];den2=[1 2 2];num3=1,den3=[2 3 3 1];t=0:0.1:10;[y1,x1,t1]=step(num1,den1,t1);subplot(1,3,1);plot(t1,y1,r);hold on;[y2,x2,t2]=step(num2,den2,t1);subplot(1,3,2);plot(t2,y2);plot(t1,y1,b);hold on;[y3,x3,t3]=step(num3,den3,t1);subplot(1,3,3);plot(t3,y3);plot(t1,y1,g) 将‘T1、T2、T3’分别标注在对应的曲线上。num1=2;den1=[1 2 2];num2=[4 2];den2=[1 2 2];num3=1,den3=[2 3 3 1];t=0:0.1:10;[y1,x1,t1]=step(num1,den1,t1);subplot(1,3,1);plot(t1,y1,r);title(T1);hold on;[y2,x2,t2]=step(num2,den2,t1);subplot(1,3,2);plot(t2,y2);plot(t1,y1,b);title(T2);hold on;[y3,x3,t3]=step(num3,den3,t1);subplot(1,3,3);plot(t3,y3);plot(t1,y1,g);title(T3); 练习4-5 一个系统的状态空间描述如下：求出 G（S）= Y（S）/U（S）；A=[-1 -1;