【正文】 10)(1 )()(23 ??????? ssssGsGs () 由此可得系統(tǒng)的閉環(huán)特征方程為： )( 23 ???? ssssD () 列出勞斯表，如下表所示： 武漢理工大學《自動控制原理》課程設計說明書 5 ns 第一列 第二列 3s 1 2s 10 1s 1 0 0s 10 0 由于勞斯表第一列中有一個系數(shù)為 10，因此第一列系數(shù)變化兩次，說明系統(tǒng)閉環(huán)特征方程有兩個正實部的根，系統(tǒng)不穩(wěn)定。 ) +20l gT 39。39。) 結果如下圖所示： 武漢理工大學《自動控制原理》課程設計說明書 18 已校正系統(tǒng)單位階躍響應T i m e ( s e c o n d s )Amplitude0 5 10 1500 . 20 . 40 . 60 . 811 . 21 . 4S y s t e m : s y s 3 _ s t e pR i s e t i m e ( s e c o n d s ) : 0 . 9 8 1S y s t e m : s y s 3 _ s t e pP e a k a m p l i t u d e : 1 . 2 2O v e r s h o o t ( % ) : 2 2A t t i m e ( s e c o n d s ) : 2 . 3 5S y s t e m : s y s 3 _ s t e pS e t t l i n g t i m e ( s e c o n d s ) : 8 . 1 8S y s t e m : s y s 3 _ s t e pF i n a l v a l u e : 1 圖 14 已校正系統(tǒng)的單位階躍響應曲線 從圖中可以看出上升時間 ，峰值時間 ，超調量 22%，調節(jié)時間（ ? =? 2%）。 總之，通過這次做課程設計，我學到了很多，無論是課內的還是課外的，這些都為今后繼續(xù)學習奠定了基礎。 繪制已校正系統(tǒng)的根軌跡圖 用 MATLAB 繪制出已校正系統(tǒng)的根軌跡圖，命令如下： rlocus(sys3) %已校正系統(tǒng)根軌跡圖 title(39。39。39。c b cL ??可由帶校正系統(tǒng)幅頻特性曲線斜率為 40dB/dec 的部分在39。39。每種方法都有不同的適用范圍，應當根據(jù)實際要求恰當?shù)倪x擇，由于本題要求采用串聯(lián)滯后 超前校正，下面將著重介紹這種方法。 用 MATLAB分析未校正和已校正系統(tǒng)的根軌跡，單位階躍響應和單位斜坡響應。學會利用MATLAB 進行建模仿真應當是大學生的一項基本技能。c? ，相角裕度 39。)cL? 為待校正系統(tǒng)的幅頻特性曲線在 39。由于相角裕度小于零，幅值裕度大于零 ,亦證明系統(tǒng)不穩(wěn)定。)= 39。) grid on 結果如下圖所示： 2 0 1 5 1 0505Magnitude (dB)104103102101100101102103 9 0 4 504590Phase (deg)校正裝置 bode 圖F r e q u e n c y ( r a d / s ) 圖 9 校正網(wǎng)絡的伯德圖 對已校正系統(tǒng)的驗證及穩(wěn)定性分析 繪制已校正系統(tǒng)的伯德圖 將校正網(wǎng)絡與待校正系統(tǒng)的串聯(lián)后，利用 MATLAB 繪出已校正系統(tǒng)的伯德圖，鍵入如下命令： 武漢理工大學《自動控制原理》課程設計說明書 13 %校正裝置 num1=conv([,1],[1,1]) den1=conv([,1],[,1]) sys1=tf(num1,den1) num2=10 den2=conv(conv([1,0],[1,1]),[,1]) sys2=tf(num2,den2) sys3=series(sys1,sys2) margin(sys3)%已校正系統(tǒng) bode 圖 [hdb3,r3,wx3,wc3]=margin(sys3) sys3_step=feedback(sys3,1) %求校正后系統(tǒng)閉環(huán)傳遞函數(shù) sys3_bandwidth=bandwidth(sys3_step)%求校正后閉環(huán)系統(tǒng)帶寬頻率 結果如下圖所示 ： 圖 10 已校正系統(tǒng)的伯德圖 命令行中的結果為： hdb3 =； r3 =； wx3 =； wc3 = 因此校正后系統(tǒng)穿越頻率為 ，截止武漢理工大學《自動控制原理》課程設計說明書 14 頻率為 。其中相角裕度的 增加意味著阻尼比增大，超調量減小，系統(tǒng)動態(tài)性能變好； 3) 校正后系統(tǒng)截止頻率從 。 從前后對比中可以看出，所得出的結論與串聯(lián)滯后 超前校正的特點一致，同時可以看出系統(tǒng)的各項性能指標之間存在著矛盾的，比如超調量和調節(jié)時間，因此在實際的操作過程中應綜合考慮校正網(wǎng)絡的特點 和所要達到的目標，合理的選取校正網(wǎng)絡及參數(shù)，通常無法一次性達到要求，需反復的計算嘗試，這也是分析法 (又稱試探法 )的缺點之一。 2345 2 ????? ???? sssss sss () 所以已校正系統(tǒng)的閉環(huán)特征方程為： 102 3 53 8 2 9 )(39。 ( 1. 41 ) +20l g( 1 1. 41 ) =0L? ? () 可以解得 ?? ，帶入到滯后 超前校正網(wǎng)絡的傳遞函數(shù)中得： ( 1+ T s) ( 1+ s)( s) = 1( 1+ T s) ( 1+ s)7. 08acaG () 式中只有 aT 一個未知參數(shù)，因此校正后系統(tǒng)的開環(huán)傳遞函數(shù)為： ))()(( )1(10)()()(39。在 MATLAB 命令窗口中輸入如下命令： sys2_step=feedback(sys2,1) %求未校正系統(tǒng)閉環(huán)傳遞函數(shù) 武漢理工大學《自動控制原理》課程設計說明書 9 step(sys2_step) %未校正系統(tǒng)單位階躍響應 H(S)=1 hold on title(‘未校正系統(tǒng)單位階躍響應 ’) 結果如下圖所示： 圖 6 待校正系統(tǒng)的單位階躍 響應曲線 從圖中可以看出系統(tǒng)的單位階躍響應呈發(fā)散震蕩形式，系統(tǒng)嚴重不穩(wěn)定。c? 處的值，39。(dB)h ； 3) 在待校正系統(tǒng)對數(shù)幅頻特性曲線上，選擇斜率從 20dB/dec 變?yōu)?40dB/dec 的交接頻率作為校正網(wǎng)絡超前部分的交接頻率 b? 。如果性能指標以單位階躍響應的峰值時間、調節(jié)時間、超調量、阻尼比、穩(wěn)態(tài)誤差等時域特征量給出時，一般采用時域法校正；如果性能指標以系統(tǒng)的相角裕度、諧振峰值、閉環(huán)帶寬、靜態(tài)誤差系數(shù)等頻域特征量給出，如本題，一般采用頻率