線性系統(tǒng)時域響應分析WORD

1、文檔可能無法思考全面，請瀏覽后下載! 1 / 10實驗實驗二二 線線性系性系統(tǒng)時統(tǒng)時域響域響應應分析分析一、實驗目的一、實驗目的1熟練掌握 step( )函數(shù)和 impulse( )函數(shù)的使用方法，研究線性系統(tǒng)在單位階躍、單位脈沖及單位斜坡函數(shù)作用下的響應。2通過響應曲線觀測特征參量和對二階系統(tǒng)性能的影響。n二、基礎知識及二、基礎知識及 MATLABMATLAB 函數(shù)函數(shù)（一）基礎知識時域分析法直接在時間域中對系統(tǒng)進行分析，可以提供系統(tǒng)時間響應的全部信息，具有直觀、準確的特點。為了研究控制系統(tǒng)的時域特性，經(jīng)常采用瞬態(tài)響應（如階躍響應、脈沖響應和斜坡響應） 。本次實驗從分析系統(tǒng)的性能指標出發(fā)，給

2、出了在 MATLAB 環(huán)境下獲取系統(tǒng)時域響應和分析系統(tǒng)的動態(tài)性能和穩(wěn)態(tài)性能的方法。用 MATLAB 求系統(tǒng)的瞬態(tài)響應時，將傳遞函數(shù)的分子、分母多項式的系數(shù)分別以 s 的降冪排列寫為兩個數(shù)組 num、den。由于控制系統(tǒng)分子的階次 m 一般小于其分母的階次 n，所以 num 中的數(shù)組元素與分子多項式系數(shù)之間自右向左逐次對齊，不足部分用零補齊，缺項系數(shù)也用零補上。1用 MATLAB 求控制系統(tǒng)的瞬態(tài)響應1）階躍響應求系統(tǒng)階躍響應的指令有： step(num,den) 時間向量 t 的范圍由軟件自動設定，階躍響應曲線隨即繪出step(num,den,t) 時間向量 t 的范圍可以由人工給定（例如t=

3、0:0.1:10）y，x=step(num,den) 返回變量 y 為輸出向量，x 為狀態(tài)向量在 MATLAB 程序中，先定義 num,den 數(shù)組，并調(diào)用上述指令，即可生成單位階躍輸入信號下的階躍響應曲線圖?？紤]下列系統(tǒng)：25425)()(2sssRsC該系統(tǒng)可以表示為兩個數(shù)組，每一個數(shù)組由相應的多項式系數(shù)組成，并且以 s的降冪排列。則 MATLAB 的調(diào)用語句：文檔可能無法思考全面，請瀏覽后下載! 2 / 10 num=0 0 25; %定義分子多項式 den=1 4 25; %定義分母多項式 step(num,den) %調(diào)用階躍響應函數(shù)求取單位階躍響應曲線 grid %畫網(wǎng)格標度線 x

4、label(t/s),ylabel(c(t) %給坐標軸加上說明 title(Unit-step Respinse of G(s)=25/(s2+4s+25) %給圖形加上標題名則該單位階躍響應曲線如圖 2-1 所示：為了在圖形屏幕上書寫文本，可以用 text 命令在圖上的任何位置加標注。例如： text(3.4,-0.06,Y1) 和 text(3.4,1.4,Y2)第一個語句告訴計算機，在坐標點 x=3.4,y=-0.06 上書寫出Y1 。類似地，第二個語句告訴計算機，在坐標點 x=3.4,y=1.4 上書寫出Y2 。若要繪制系統(tǒng) t 在指定時間（0-10s）內(nèi)的響應曲線，則用以下語句：n

5、um=0 0 25; den=1 4 25; t=0:0.1:10; step(num,den,t) 即可得到系統(tǒng)的單位階躍響應曲線在 0-10s 間的部分，如圖 2-2 所示。 2）脈沖響應 求系統(tǒng)脈沖響應的指令有：圖 2-1 二階系統(tǒng)的單位階躍響應 圖 2-2 定義時間范圍的單位階躍響應文檔可能無法思考全面，請瀏覽后下載! 3 / 10 impulse (num,den) 時間向量 t 的范圍由軟件自動設定，階躍響應曲線隨即繪出 impulse (num,den,t) 時間向量 t 的范圍可以由人工給定（例如t=0:0.1:10）y,x=impulse(num,den) 返回變量 y 為輸

6、出向量，x 為狀態(tài)向量y,x,t=impulse(num,den,t) 向量 t 表示脈沖響應進行計算的時間例：試求下列系統(tǒng)的單位脈沖響應： 12 . 01)()()(2sssGsRsC在 MATLAB 中可表示為 num=0 0 1; den=1 0.2 1; impulse(num,den) grid title(Unit-impulse Response of G(s)=1/(s2+0.2s+1)由此得到的單位脈沖響應曲線如圖 2-3 所示： 求脈沖響應的另一種方法應當指出，當初始條件為零時，G (s)的單位脈沖響應與 sG(s)的單位階躍響應相同?？紤]在上例題中求系統(tǒng)的單位脈沖響應，因

7、為對于單位脈沖輸入量，R(s)=1 所以sssssssGsCsRsC112 . 012 . 01)()()()(22因此，可以將 G(s)的單位脈沖響應變換成 sG(s)的單位階躍響應。圖 2-3 二階系統(tǒng)的單位脈沖響應文檔可能無法思考全面，請瀏覽后下載! 4 / 10向 MATLAB 輸入下列 num 和 den，給出階躍響應命令，可以得到系統(tǒng)的單位脈沖響應曲線如圖 2-4 所示。 num=0 1 0; den=1 0.2 1; step(num,den) grid title(Unit-step Response of sG(s)=s/(s2+0.2s+1)3）斜坡響應MATLAB 沒有直

8、接調(diào)用求系統(tǒng)斜坡響應的功能指令。在求取斜坡響應時，通常利用階躍響應的指令?；趩挝浑A躍信號的拉氏變換為 1/s，而單位斜坡信號的拉氏變換為 1/s2。因此，當求系統(tǒng) G(s)的單位斜坡響應時，可以先用 s 除G(s)，再利用階躍響應命令，就能求出系統(tǒng)的斜坡響應。例如，試求下列閉環(huán)系統(tǒng)的單位斜坡響應。 11)()(2sssRsC對于單位斜坡輸入量，R(s)=1/s2 ，因此 ssssssssC1) 1(1111)(222在 MATLAB 中輸入以下命令，得到如圖 2-5 所示的響應曲線： num=0 0 0 1; den=1 1 1 0;step(num,den)title(Unit-Ramp

9、Response Cuve for System G(s)=1/(s2+s+1)圖 2-4 單位脈沖響應的另一種表示法圖 2-5 單位斜坡響應文檔可能無法思考全面，請瀏覽后下載! 5 / 102. 特征參量和對二階系統(tǒng)性能的影響n標準二階系統(tǒng)的閉環(huán)傳遞函數(shù)為： 2222)()(nnnsssRsC二階系統(tǒng)的單位階躍響應在不同的特征參量下有不同的響應曲線。1）對二階系統(tǒng)性能的影響設定無阻尼自然振蕩頻率，考慮 5 種不同的值：)/( 1sradn=0,0.25,0.5,1.0 和 2.0，利用 MATLAB 對每一種求取單位階躍響應曲線，分析參數(shù)對系統(tǒng)的影響。為便于觀測和比較，在一幅圖上繪出 5 條

10、響應曲線（采用“hold”命令實現(xiàn)）。 num=0 0 1; den1=1 0 1; den2=1 0.5 1; den3=1 1 1; den4=1 2 1; den5=1 4 1;t=0:0.1:10; step(num,den1,t) grid text(4,1.7,Zeta=0); hold step(num,den2,t) text (3.3,1.5,0.25) step(num,den3,t) text (3.5,1.2,0.5) step(num,den4,t) text (3.3,0.9,1.0) step(num,den5,t) text (3.3,0.6,2.0) titl

11、e(Step-Response Curves for G(s)=1/s2+2(zeta)s+1)由此得到的響應曲線如圖 2-6 所示：文檔可能無法思考全面，請瀏覽后下載! 6 / 102）n對二階系統(tǒng)性能的影響同理，設定阻尼比時，25. 0當分別取 1,2,3 時，利用 MATLAB 求取單位階躍響應曲線，分析參數(shù)對系nn統(tǒng)的影響。num1=0 0 1; den1=1 0.5 1; t=0:0.1:10; step(num1,den1,t); grid; hold ontext(3.1,1.4,wn=1)num2=0 0 4; den2=1 1 4;step(num2,den2,t); hol

12、d ontext(1.7,1.4,wn=2)num3=0 0 9; den3=1 1.5 9;step(num3,den3,t); hold ontext(0.5,1.4,wn=3)由此得到的響應曲線如圖 2-7 所示：圖 2-6 不同時系統(tǒng)的響應曲線圖 2-7 不同時系統(tǒng)的響應曲線n文檔可能無法思考全面，請瀏覽后下載! 7 / 10三、實驗內(nèi)容三、實驗內(nèi)容1觀察函數(shù) step( )調(diào)用格式，假設系統(tǒng)的傳遞函數(shù)模型為 146473)(2342sssssssG繪制出系統(tǒng)的階躍響應曲線。num=0 0 1 3 7;den=1 4 6 4 1;t=0:0.1:10;step(num,den,t)gr

13、idxlabel(t/s),ylabel(c(t)title(Unit-impulse Response of G(s)=1/(s2+0.2s+1)2對典型二階系統(tǒng)2222)(nnnsssG文檔可能無法思考全面，請瀏覽后下載! 8 / 101）分別繪出，分別取 0,0.25,0.5,1.0 和 2.0 時的單位階)/(2sradn躍響應曲線，分析參數(shù)對系統(tǒng)的影響 num=0 0 4;den1=1 0 4;den2=1 1 4;den3=1 2 4;den4=1 4 4;den5=1 8 4;t=0:0.1:10; step(num,den1,t)gridtext(1,1.7,Zeta=0);

14、holdstep(num,den2,t)text(1.3,1.5,0.25); step(num,den3,t)text(1.34,1.2,0.5); step(num,den4,t)text(1.5,0.9,1.0); step(num,den5,t)text(1.8,0.6,2.0); title(Step-Response Curves for G(s)=1/(s2+2(zeta)s+1)文檔可能無法思考全面，請瀏覽后下載! 9 / 102）繪制出當=0.25, 分別取 1,2,4,6 時單位階躍響應曲線，分析參數(shù)n對系統(tǒng)的影響。nnum1=0 0 1; den1=1 0.5 1; t=

15、0:0.1:10; step(num1,den1,t); grid; hold on text(3.32,1.4,wn=1) num2=0 0 4; den2=1 1 4; step(num2,den2,t); hold ontext(3.26,0.8,wn=2) num3=0 0 16; den3=1 2 16; step(num3,den3,t); hold on text(3.23,0.962,wn=4) num4=0 0 36; den4=1 3 36; step(num4,den4,t); hold on text(1.4,1.09,wn=6) step(num3,den3,t); hold ontext(0.1,1.4,wn=6)文檔可能無法思考全面，請瀏覽后下載! 10 / 10四、