1、廣西大學實驗報告紙【實驗時間】2014年06月15日【實驗地點】(課外)【實驗目的】1、掌握線性系統狀態空間的標準型、解及其模型轉換。【實驗設備與軟件】1、MATLAB數值分析軟件【實驗原理】Matlab提供了非常豐富的線性定常連續系統的狀態空間模型求解(即系統運動軌跡的計算)的功能,主要的函數有 、階躍響應函數step()可用于計算在單位階躍輸入和零初始狀態(條件)下傳遞函數模型 的輸出響應，或狀態空間模型的狀態和輸出響應，其主要調用格式為ste p(sys,t)y,t = ste p( sys,t)y,t,x = ste p( sys,t) 、脈沖激勵下的仿真函數imp ulse()可用于

2、計算在脈沖刺激輸入下傳遞函數模型的輸出 響應,或狀態空間模型的狀態和輸出響應，其主要調用格式為imp ulse(sys,t)y,t = imp ulse(sys,t)y,t,x = impu lse(sys,t) 、任意輸入激勵下的仿真函數lsim()可用于計算在給定的輸入信號序列(輸入信號函數 的采樣值)下傳遞函數模型的輸出響應，其主要調用格式為lsim(sys,u,t,x0)y,t,x = lsim(sys,u,t,x0)【實驗內容、方法、過程與分析】已知線性系統1、利用Matlab求零狀態下的階躍響應(包括狀態和輸出)，生成兩幅圖：第一幅繪制各狀 態響應曲線并標注；第二幅繪制輸出響應曲線

3、。狀態響應曲線：A=-21 19 -20;19 -21 20;40 -40 -40;B=0；1；2;%輸入狀態空間模型各矩陣，若沒有相應值，可賦空矩陣%輸入初始狀態%構造傳遞函數%繪以時間為橫坐標的狀態響應曲線圖C=1 0 2; D=0;X0=0;0;0;sys=ss(A,B,C,D); y,x,t=ste p(sys);plot(t,x);grid;titleC狀態響應曲線')輸出響應曲線：A=-21 19 -20;19 -21 20;40 -40 -40;B=0；1；2;C=1 0 2;D=0;X0=0;0;0n um,de n=ss2tf(A,B,C,D,1); sys=tf(

4、nu m,de n);ste p( sys)gridtitleC輸出響應曲線')$rII2、利用Matlab求零狀態下的沖激響應(包括狀態和輸出)，生成兩幅圖：第一幅繪制各狀 態響應曲線并標注；第二幅繪制輸出響應曲線。狀態響應曲線：A=-21 19 -20;19 -21 20;40 -40 -40;B=0;1;2;C=1 0 2;%輸入狀態空間模型各矩陣，若沒有相應值,可賦空矩陣%輸入初始狀態%構造傳遞函數D=; x0=0;0;0; sys=ss(A,B,C,D);y,x,t= imp ulse(sys);plot(t,x); grid;titleC狀態響應曲線')1- : -

5、 - -II11I1! i1111111“丄_亠-:1-一亠=.4.匸=i2一 2 Ji輸出響應曲線：A=-21 19 -20;19 -21 20;40 -40 -40;B=0；1；2;C=1 0 2;D=0;X0=0;0;0n um,de n=ss2tf(A,B,C,D,1); sys=tf( nu m,de n);impu Ise(sys);grid;titleC輸出響應曲線')83、若控制輸入為，且初始狀態為，求系統的響應，要求a. 在simulink只能夠畫出模型求響應，生成兩幅圖：第一幅繪制各狀態響應曲線并標注； 第二幅繪制輸出響應曲線。b. 編寫.m文件求響應，生成兩幅圖：

6、第一幅繪制各狀態響應曲線并標注；第二幅繪制輸出響應曲線。狀態響應曲線：t=0:0.02:5;if t>=3u=1;elseu=1+ex p(-t).*cos(5*t);endA=-21 19 -20;19 -21 20;40 -40 -40;B=0;1;2;%輸入狀態空間模型各矩陣，若沒有相應值,可賦空矩陣 %輸入初始狀態%就是個條件判斷，只有t=0的時候，U才為“1”%構造傳遞函數C=1 0 2;D=0;X0=0.2;0.2;0.2;u=(t=0); sys=ss(A,B,C,D);y,t,x=lsim(sys,u,t,X0);plot(t,x);grid;titleC狀態響應曲線&#

7、39;)plot(t,y);grid;titleC輸出響應曲線')4、以階躍輸入情況下的，分析各模塊對響應有什么影響。5、求系統的傳遞函數在MATLAB 軟件Comma nd Window窗口中輸入以下程序A=-21 19 -20;19 -21 20;40 -40 -40;B=0;1;2;C=1 0 2;D=0;n um,de n=ss2tf(A,B,C,D,1);prin tsys( nu m,de n)程序運行結果為>> clear» A=-21 ig '20:19 -21 20:40 -40 -40.C=E1 0 2：nuitj derL=s5 2-

8、t-f (& Bj Cj 必 ；pi;int sys (num., den)列3 + S2 5"2 + 3商0 s + 拠QQ6、若采用K增益負反饋，繪制閉環根軌跡圖，并對根軌跡加以描述說明。A=-21 19 -20;19 -21 20;40 -40 -40;B=0；1；2;C=1 0 2;D=0;n um,de n=ss2tf(A,B,C,D,1); rlocus( nu m,de n);gridtitle('K增益負反饋閉環根軌跡圖')采用K增益負反饋，畫出如圖所示的根軌跡圖。由圖可知，共有3條根軌跡，第一條最終趨于原點；第二條收斂在2060之間；第三條最

9、終趨于無窮遠處。7、在Matlab中繪制Bode圖和Nyquist圖，并對圖給予說明。 繪制Bode圖：A=-21 19 -20;19 -21 20;40 -40 -40;B=0；1；2;C=1 0 2;D=0;sys=tf( nu m,de n)bode( nu m,de n)gridtitle('Bode 圖')IQr3S 亠: / 一- - -1 - :二 -1 I - - - - J . I -' 一 r - 二7八一匯出的波特圖如圖所示，由圖可知，對復制響應分析可得，交越頻率在轉折頻率之后， 故復制的變化主要發生在低頻段。對相頻特性進行分析，可知此系統的相頻特性角度均為負 值，并且最后的相角是趨于-9