實(shí)驗(yàn)四線性控制系統(tǒng)的時(shí)域響應(yīng)分析

1、實(shí)驗(yàn)四線性控制系統(tǒng)的時(shí)域響應(yīng)分析一、實(shí)驗(yàn)?zāi)康? 熟悉MATLAB有關(guān)命令的用法；2 用MATLAB系統(tǒng)命令對(duì)給定系統(tǒng)進(jìn)行時(shí)域分析；二、實(shí)驗(yàn)內(nèi)容求連續(xù)系統(tǒng)的單位階躍響應(yīng)有關(guān)命令的用法：命令格式：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)練習(xí)4-1.給定系統(tǒng)的傳遞函數(shù)如下：求該系統(tǒng)的階躍響應(yīng)曲線,記錄超調(diào)量、上升時(shí)間、過(guò)渡過(guò)程時(shí)間。>> num=25;den=1 4 25;step(num,den)練習(xí)4-2.

2、已知系統(tǒng)的開(kāi)環(huán)傳遞函數(shù)為：求出該系統(tǒng)在單位負(fù)反饋下的階躍響應(yīng)曲線，記錄超調(diào)量、上升時(shí)間、過(guò)渡過(guò)程時(shí)間。>> num0 = 20 ;den0=1 8 36 40 0;num,den=cloop(num0,den0);step(num,den)練習(xí)4-3 已知系統(tǒng)的傳遞函數(shù)為：+ 求系統(tǒng)的階躍響應(yīng)； 階躍響應(yīng)曲線線型用“*”號(hào)表示； 躍響應(yīng)圖應(yīng)加上橫坐標(biāo)名、縱坐標(biāo)名和標(biāo)題名，并加上網(wǎng)格線。num1=6.3233*conv(1 1.4235,1 1.4235);den1=1 0;num2=1;den2=conv(1 0,conv(1 1 ,1 5);sys1=tf(num1,den1);

3、sys2=tf(num2,den2);sys0=series(sys1,sys2);figure(1);subplot(2,1,1);step(feedback(sys0,1,1),'*')subplot(2,1,2);step(feedback(sys0,1,-1),'*')練習(xí)4-4 求T1、T2、T3系統(tǒng)的階躍響應(yīng)； 將 T1、T2、T3系統(tǒng)的階躍響應(yīng)圖畫在同一窗口內(nèi)； T1、T2、T3系統(tǒng)的階躍響應(yīng)曲線分別用不同的線形和顏色表示； 將T1、T2、T3分別標(biāo)注在對(duì)應(yīng)的曲線上。num1=2;den1=1 2 2;num2=4 2;den2=1 2 2;num

4、3=1;den3=2 3 3 1;figure(1)hold onstep(num1,den1,'*r')step(num2,den2,'+g')step(num3,den3,'-b')gtext('T1')gtext('T2')gtext('T3')練習(xí)4-5 一個(gè)系統(tǒng)的狀態(tài)空間描述如下： 求出 G（S）= Y（S）/U（S）；繪制該狀態(tài)方程的單位階躍響應(yīng)曲線。>> A=-1 -1;6.5 0;B=1 1;1 0;C=1 0; 0 1 ;D=zeros(2);num1,den1=ss2

5、tf(A,B,C,D,1)num2,den2=ss2tf(A,B,C,D,2)num1 = 0 1.0000 -1.0000 0 1.0000 7.5000den1 = 1.0000 1.0000 6.5000num2 = 0 1.0000 0 0 0 6.5000den2 = 1.0000 1.0000 6.5000step(num1,den1) step(num2,den2)練習(xí)4-6典型二階欠阻尼系統(tǒng)的傳遞函數(shù)為：極點(diǎn)位置：式中：設(shè)a=1, =0.5,1,5 ,求階躍響應(yīng)；wa=1;n=1;for theta=0.5 1 5 subplot(2,2,n) num=wa2+theta2;

6、den=1 2*theta num; step(num,den); n=n+1; str0=strcat('Step Response ¦Ò=',num2str(theta); title(str0);end設(shè)=1 , a=0.5,1,5 ，求階躍響應(yīng)；theta=1;n=1;for wa=0.5 1 5 subplot(2,2,n) num=wa2+theta2; den=1 2*theta num; step(num,den); n=n+1; str0=strcat('Step Response ¦Øa=',num2s

7、tr(wa); title(str0);end設(shè)：求階躍響應(yīng)；kethy=1/sqrt(2);Wn = sqrt(2)/2 sqrt(2) 5/sqrt(2);n=1;for a=Wn Wa=a*sqrt(1-kethy2); theta=a*kethy; subplot(2,2,n) num=Wa2+theta2; den=1 2*theta num; step(num,den); n=n+1; str0=strcat('Step Response ¦Øn =',num2str(a); title(str0);end 設(shè) 求階躍響應(yīng)；delta = pi/6, pi/4, pi/3 ;Wn = sqrt(2);n=1;for a=delta kethy=cos(a); Wa=a*sqrt(1-kethy2); theta=a*kethy; subplot(2,2,n) num=Wa2+theta2; den=1 2*theta num; step(num,den); n=n+1; str0