計算機仿真第四章控制系統(tǒng)的分析方法-wenkub
2023-05-21 14:33:58 本頁面
【正文】 ,?? 5?nw ，求系統(tǒng)的階躍響應(yīng)曲線。 x為系統(tǒng)返回的狀態(tài)軌跡。在 MATLAB的控制系統(tǒng)工具箱中提供了求取這兩種輸入下系統(tǒng)響應(yīng)的函數(shù)。這樣如果找到了實部大于 0的極點，則會將該極點的序號返回到 ii下。 例 已知某系統(tǒng)的模型如右所示： ? ? uxyuxx7165210016127587403622121??????????????? ???????????????????要求判斷系統(tǒng)的穩(wěn)定性及系統(tǒng)是否為最小相位系統(tǒng)。 ?若連續(xù)時間系統(tǒng)的 全部零極點 都位于 S左半平面；或若離散時間系統(tǒng)的全部零極點都位于 Z平面單位圓內(nèi)，則系統(tǒng)是最小相位系統(tǒng)。 ? MATLAB控制系統(tǒng)工具箱和 SIMULINK輔助環(huán)境的出現(xiàn)，給控制系統(tǒng)分析帶來了福音。 ?控制系統(tǒng)的分析包括系統(tǒng)的穩(wěn)定性分析、時域分析、頻域分析及根軌跡分析。 一、系統(tǒng)穩(wěn)定及最小相位系統(tǒng)判據(jù) 直接判別 MATLAB提供了直接求取系統(tǒng)所有零極點的函數(shù)，因此可以直接根據(jù)零極點的分布情況對系統(tǒng)的穩(wěn)定性及是否為最小相位系統(tǒng)進行判斷。 例 系統(tǒng)模型如下所示，判斷系統(tǒng)的穩(wěn)定性，以及系統(tǒng)是否為最小相位系統(tǒng)。如果最終的結(jié)果里 ii的元素個數(shù)大于 0，則認為找到了不穩(wěn)定極點，因而給出系統(tǒng)不穩(wěn)定的提示，若產(chǎn)生的 ii向量的元素個數(shù)為 0，則認為沒有找到不穩(wěn)定的極點，因而得出系統(tǒng)穩(wěn)定的結(jié)論。 一、時域分析的一般方法 ?求取系統(tǒng)單位階躍響應(yīng)： step() ?求取系統(tǒng)的沖激響應(yīng)： impulse() step()函數(shù)的用法 ?y=step(num,den,t)：其中 num和 den分別為系統(tǒng)傳遞函數(shù)描述中的分子和分母多項式系數(shù)， t為選定的仿真時間向量，一般可以由 t=0:step:end等步長地產(chǎn)生出來。 ?如果對具體的響應(yīng)值不感興趣，而只想繪制系統(tǒng)的階躍響應(yīng)曲線，可調(diào)用以下的格式： step(num,den)； step(num,den,t)； step(A,B,C,D,iu,t)； step(A,B,C,D,iu)； 例 e xp 4_ 3. m 已知系統(tǒng)的開環(huán)傳遞函數(shù)為：sssssG o4036820)(234 ???? 求系統(tǒng)在單位負反饋下的階躍響應(yīng)曲線。 例 e xp 4_ 6 . m 已知某閉環(huán)系統(tǒng)的傳遞函數(shù)為：2510)(23 ?????sssssG 求其階躍響應(yīng)曲線。 ?一般來說，先不指定仿真時間，由 MATLAB自己確定，然后根據(jù)結(jié)果，最后確定合適的仿真時間。 它們的調(diào)用格式與 step、 impulse類似，可以通過 help命令來察看自學。 )1( ?ssd 1 ＋ e s + _ R ( s) C ( s) 由圖可得閉環(huán)傳遞函數(shù)為：dsedsdsGc?????)1()(2，其為典型二階系統(tǒng)。采用這種方法可直觀地表達出系統(tǒng)的頻率特性，分析方法比較簡單，物理概念比較明確，對于諸如防止結(jié)構(gòu)諧振、抑制噪聲、改善系統(tǒng)穩(wěn)定性和暫態(tài)性能等問題，都可以從系統(tǒng)的頻率特性上明確地看出其物理實質(zhì)和解決途經(jīng)。 MATLAB提供了函數(shù) bode()來繪制系統(tǒng)的波特圖，其用法如下： ?bode(a,b,c,d)：自動繪制出系統(tǒng)的一組 Bode圖，它們是針對連續(xù)狀態(tài)空間系統(tǒng) [a,b,c,d]的每個輸入的 Bode圖。 ?bode(a,b,c,d,iu,w)或 bode(num,den,w)：可利用指定的角頻率矢量繪制出系統(tǒng)的波特圖。 MATLAB提供了函數(shù) nyquist()來繪制系統(tǒng)的極坐標圖，其用法如下： ?nyquist(a,b,c,d)：繪制出系統(tǒng)的一組 Nyquist曲線，每條曲線相應(yīng)于連續(xù)狀態(tài)空間系統(tǒng) [a,b,c,d]的輸入 /輸出組合對。 ?nyquist(a,b,c,d,iu,w)或 nyquist(num,den,w)：可利用指定的角頻率矢量繪制出系統(tǒng)的極坐標圖。 二、常用頻域分析函數(shù) MATLAB除了提供前面介紹的基本頻域分析函數(shù)外，還提供了大量在工程實際中廣泛應(yīng)用的庫函數(shù)，由這些函數(shù)可以求得系統(tǒng)的各種頻率響應(yīng)曲線和 特征值。 ?幅值裕度是在相角為 180度處使開環(huán)增益為 1的增益量，如在 180度相頻處的開環(huán)增益為 g，則幅值裕度為 1/g；若用分貝值表示幅值裕度，則等于：20*log10(g)。類似， margin(a,b,c,d)可以計算出連續(xù)狀態(tài)空間系統(tǒng)表示的幅值裕度和相角裕度并繪制相應(yīng)波特圖。 ?[h,w]=freqs(b,a)自動設(shè)定 200個頻率點來計算頻率響
公司管理相關(guān)推薦
鄂ICP備17016276號-1