北郵計(jì)算機(jī)仿真設(shè)計(jì)報(bào)告（精）

1、計(jì)算機(jī)控制仿真課程設(shè)計(jì)姓名 學(xué)院 自動(dòng)化學(xué)院 專業(yè) 自動(dòng)化專業(yè) 班級(jí) 學(xué)號(hào) 10212006（03） 10212012（09） 10212014（12） 指導(dǎo)教師 2013年7月目錄一、引言11.1課程設(shè)計(jì)簡(jiǎn)介11.2選題及小組分工1二、數(shù)字PID閉環(huán)直流電機(jī)調(diào)速控制系統(tǒng)的設(shè)計(jì)和仿真實(shí)現(xiàn)22.1課程設(shè)計(jì)的基本理論知識(shí)2數(shù)字PID的控制算法2調(diào)節(jié)器參數(shù)對(duì)控制性能的影響4采樣時(shí)間T的選擇42.2 課程設(shè)計(jì)的設(shè)計(jì)要求52.3 課程設(shè)計(jì)的總體方案設(shè)計(jì)62.4 課程的詳細(xì)設(shè)計(jì)6模擬PID控制器的設(shè)計(jì)6數(shù)字PID控制器的設(shè)計(jì)9仿真結(jié)果分析122.5調(diào)試Simulink仿真12當(dāng)Kp變化時(shí)12當(dāng)Ti變化，即

2、Ki變化時(shí)14當(dāng)Td變化，即Kd變化時(shí)182.6設(shè)計(jì)總結(jié)21三、大林算法計(jì)算機(jī)控制系統(tǒng)設(shè)計(jì)及仿真實(shí)現(xiàn)213.1 設(shè)計(jì)要求213.2 總體方案設(shè)計(jì)223.3 基本理論知識(shí)223.4 詳細(xì)設(shè)計(jì)253.5 調(diào)試Simulink仿真26四、二階彈簧阻尼系統(tǒng)的PID控制器設(shè)計(jì)及其參數(shù)整定304.1 設(shè)計(jì)要求304.2 總體方案設(shè)計(jì)314.3 基本理論知識(shí)314.4 詳細(xì)設(shè)計(jì)344.5 分析與總結(jié)40五、二階系統(tǒng)串聯(lián)校正裝置的設(shè)計(jì)與分析405.1 設(shè)計(jì)要求405.2 基本理論知識(shí)405.3 詳細(xì)設(shè)計(jì)425.4 調(diào)試Simulink仿真455.5 設(shè)計(jì)總結(jié)46六、單級(jí)倒立擺的最優(yōu)控制器設(shè)計(jì)476.1 設(shè)計(jì)要

3、求476.2 總體方案設(shè)計(jì)476.3 基本理論知識(shí)486.4 詳細(xì)設(shè)計(jì)496.6 調(diào)試Simulink仿真566.6 設(shè)計(jì)總結(jié)57參考文獻(xiàn)58一、 引言1.1課程設(shè)計(jì)簡(jiǎn)介計(jì)算機(jī)仿真技術(shù)是一種借助高速、大存儲(chǔ)量數(shù)字計(jì)算機(jī)及相關(guān)技術(shù)，對(duì)復(fù)雜真實(shí)系統(tǒng)的運(yùn)行過程或狀態(tài)進(jìn)行數(shù)字化模擬的技術(shù)。隨著技術(shù)的發(fā)展，計(jì)算機(jī)仿真的應(yīng)用領(lǐng)域越來越廣，在通信，控制，航天等領(lǐng)域得到了廣泛的應(yīng)用。MATLAB是一套高性能的數(shù)值計(jì)算和可視化軟件，它集數(shù)值分析、矩陣運(yùn)算和圖形顯示于一體，構(gòu)成了一個(gè)方便的界面友好的用戶環(huán)境。由控制領(lǐng)域?qū)＜彝瞥龅腗ATLAB工具箱之一的控制系統(tǒng)(Control System)，在控制系統(tǒng)計(jì)算機(jī)輔助

4、分析與設(shè)計(jì)方面獲得了廣泛的應(yīng)用，并且MATLAB工具箱的內(nèi)容還在不斷增加，應(yīng)用范圍也越來越寬?？刂葡到y(tǒng)的分析與設(shè)計(jì)方法，不論是古典的還是現(xiàn)代的，都是以數(shù)學(xué)模型為基礎(chǔ)進(jìn)行的。MATLAB可以用于以傳遞函數(shù)形式描述的控制系統(tǒng)。在本文中，我們將以數(shù)字PID閉環(huán)直流電機(jī)調(diào)速控制系統(tǒng)的設(shè)計(jì)和仿真實(shí)現(xiàn)、大林算法計(jì)算機(jī)控制系統(tǒng)設(shè)計(jì)、二階彈簧阻尼系統(tǒng)的PID控制器設(shè)計(jì)及其參數(shù)整定、二階系統(tǒng)串聯(lián)校正裝置的設(shè)計(jì)與分析、單級(jí)倒立擺的最優(yōu)控制器設(shè)計(jì)為例，說明如何使用MATLAB進(jìn)行輔助分析。Simulink仿真、計(jì)算傳遞函數(shù)零極點(diǎn)、計(jì)算閉環(huán)系統(tǒng)對(duì)單位階躍輸入的響應(yīng)等等，都會(huì)有所涉及。1.2選題及小組分工1.2.1選題

5、選題：第一題：數(shù)字PID閉環(huán)直流電機(jī)調(diào)速控制系統(tǒng)的設(shè)計(jì)和仿真實(shí)現(xiàn) 第三題：大林算法計(jì)算機(jī)控制系統(tǒng) 第四題：二階彈簧阻尼系統(tǒng)的PID控制器設(shè)計(jì)及其參數(shù)整定 第五題：二階系統(tǒng)串聯(lián)校正裝置的設(shè)計(jì)與分析 第六題：?jiǎn)渭?jí)倒立擺的最優(yōu)控制器設(shè)計(jì)1.2.2分工分工：組長(zhǎng) 10212006 負(fù)責(zé)第一、六題； 組員 10212014 負(fù)責(zé)第三、四題； 10212012 負(fù)責(zé)第五題。二、數(shù)字PID閉環(huán)直流電機(jī)調(diào)速控制系統(tǒng)的設(shè)計(jì)和仿真實(shí)現(xiàn)2.1課程設(shè)計(jì)的基本理論知識(shí)2.1.1 數(shù)字PID的控制算法在模擬調(diào)節(jié)系統(tǒng)中，PID算法的表達(dá)式為：式中：調(diào)節(jié)器的輸出信號(hào)；調(diào)節(jié)器輸入的偏差信號(hào)，它等于測(cè)量值與給定值之差；調(diào)節(jié)器的比

6、例系數(shù)；調(diào)節(jié)器的時(shí)間積分常數(shù)；調(diào)節(jié)器的時(shí)間微分常數(shù)。對(duì)上式進(jìn)行離散化，用數(shù)字形式的差分方程代替連續(xù)系統(tǒng)的微分方程，此時(shí)積分項(xiàng)和微分項(xiàng)可用求和及增量式表示：式中：采樣周期；e(k)第k次采樣時(shí)的偏差值；e(k-1)第k-1次采樣時(shí)的偏差值；k采樣序號(hào)，k=0,1,2u(k)第k次采樣時(shí)調(diào)節(jié)器的輸出。則有：上面兩式相減：式中：積分系數(shù)；微分系數(shù)。位置式PID控制算法：圖2-1 位置式PID控制算法結(jié)構(gòu)框圖增量式PID控制算法圖2-2 增量式PID控制算法結(jié)構(gòu)框圖2.1.2 PID調(diào)節(jié)器參數(shù)對(duì)控制性能的影響1、不同Kp對(duì)控制性能的影響(1) 對(duì)動(dòng)態(tài)性能的影響比例控制參數(shù)Kp加大，使系統(tǒng)的動(dòng)作靈敏，速

7、度加快；Kp偏大，振蕩次數(shù)加多，調(diào)節(jié)時(shí)間加長(zhǎng)；當(dāng)Kp太大時(shí)，系統(tǒng)會(huì)趨于不穩(wěn)定；若Kp太小，又會(huì)使系統(tǒng)的動(dòng)作緩慢。(2) 對(duì)穩(wěn)態(tài)性能的影響加大比例控制系數(shù)Kp，在系統(tǒng)穩(wěn)定的情況下，可以減小穩(wěn)態(tài)誤差提高控制精度，卻不能完全消除穩(wěn)態(tài)誤差。2、積分控制參數(shù)TI對(duì)控制性能的影響(1) 對(duì)動(dòng)態(tài)性能的影響積分控制參數(shù)TI通常使系統(tǒng)的穩(wěn)定性下降。TI太小系統(tǒng)將不穩(wěn)定，TI偏小，振蕩次數(shù)較多，TI太大，對(duì)系統(tǒng)性能的影響減少。(2) 對(duì)穩(wěn)態(tài)性能的影響積分控制參數(shù)TI能消除系統(tǒng)的穩(wěn)態(tài)誤差，提高控制系統(tǒng)的控制精度。但TI太大時(shí)，積分作用呆太弱，以至不能減少穩(wěn)態(tài)誤差。3、微分控制參數(shù)TD對(duì)控制性能的影響微分控制可以改善

8、動(dòng)態(tài)特性，如超調(diào)量減少，調(diào)節(jié)時(shí)間縮短，允許加大比例控制，使穩(wěn)態(tài)誤差減小，提高控制精度。當(dāng)TD偏大時(shí)，超調(diào)量較大，當(dāng)TD偏小時(shí)，超調(diào)量也較大，調(diào)節(jié)時(shí)間也較長(zhǎng)。2.1.3 采樣時(shí)間T的選擇（1）T越小，隨動(dòng)性和抗干擾性能越好。（2）必須滿足采樣定理的要求，對(duì)于隨動(dòng)系統(tǒng)，為系統(tǒng)的開環(huán)截止頻率。若單路采樣時(shí)間為，則采樣周期，N為測(cè)量控制回路數(shù)。（3）選擇采樣周期T太小，將使微分積分作用不明顯。（4）快速系統(tǒng)的T應(yīng)取小，反之，T可取大些。（5）執(zhí)行機(jī)構(gòu)動(dòng)作慣性大時(shí)，T應(yīng)取大些。2.2 課程設(shè)計(jì)的設(shè)計(jì)要求已知某晶閘管直流單閉環(huán)調(diào)速系統(tǒng)的轉(zhuǎn)速控制器選用 PID控制器，結(jié)構(gòu)如圖1-3所示。圖2-3 某晶閘管直

9、流單閉環(huán)調(diào)速系統(tǒng)結(jié)構(gòu)框圖要求：1、運(yùn)用 MATLAB/Simulink 軟件對(duì)控制系統(tǒng)進(jìn)行建模并對(duì)模塊進(jìn)行參數(shù)設(shè)置；2、封裝PID模塊的控制圖；3、使用期望特性法來確定 Kp、 Ti、Td以及采樣周期 T，期望系統(tǒng)對(duì)應(yīng)的閉環(huán)特征根為：-300，-300，-30+j30和-30-j30，觀察其單位階躍響應(yīng)曲線，得出仿真結(jié)果并進(jìn)行仿真分析；4、記錄在改變PID控制某一控制參數(shù)（比例系數(shù)或積分系數(shù)或微分系數(shù)）時(shí)，該系統(tǒng)對(duì)應(yīng)的階躍響應(yīng)曲線的變化，并觀察闡述發(fā)生這種變化的規(guī)律；5、總結(jié)P、I、D控制參數(shù)的改變對(duì)系統(tǒng)控制效果的影響。2.3 課程設(shè)計(jì)的總體方案設(shè)計(jì)首先，運(yùn)用MATLAB/Simulink軟件

10、對(duì)控制系統(tǒng)進(jìn)行建模并對(duì)模塊進(jìn)行參數(shù)設(shè)置，然后使用期望特性法來確定Kp、Ti、Td以及采樣周期T，得到期望系統(tǒng)對(duì)應(yīng)的閉環(huán)特征根，最后通過改變PID的某一個(gè)參數(shù)，觀察此參數(shù)變化對(duì)系統(tǒng)的影響。2.4 課程的詳細(xì)設(shè)計(jì)2.4.1 模擬PID控制器的設(shè)計(jì)假設(shè)PID控制器的傳遞函數(shù)為：閉環(huán)系統(tǒng)的傳遞函數(shù)為：希望的閉環(huán)極點(diǎn)為-300，-300，-30+j30，-30-j30，得到的期望特征方程為：兩極點(diǎn)方程對(duì)應(yīng)項(xiàng)系數(shù)相等，可解得：1、設(shè)計(jì)框圖圖2-4 Simulink設(shè)計(jì)框圖封裝PID模塊的控制圖如圖2-5所示圖2-5 封裝PID模塊的控制圖2、系統(tǒng)仿真期望的系統(tǒng)的參數(shù)設(shè)置如圖2-6所示：圖2-6 期望的系統(tǒng)

11、的參數(shù)設(shè)置仿真結(jié)果如圖2-7所示：圖2-7 仿真結(jié)果超調(diào)量為：6.82% 上升時(shí)間：0.0121s 調(diào)節(jié)時(shí)間：0.0784s位置式PID控制算法：兩式相減，得：式中積分系數(shù)；微分系數(shù)。由前面算得：由經(jīng)驗(yàn)法得到T為0.002s，則：1、設(shè)計(jì)框圖圖2-8 Simulink設(shè)計(jì)框圖封裝PID模塊的控制圖如圖2-9所示：圖2-9封裝PID模塊的控制圖2、系統(tǒng)仿真期望的系統(tǒng)的參數(shù)設(shè)置如圖2-10所示：圖2-10 期望的系統(tǒng)的參數(shù)設(shè)置仿真結(jié)果如圖2-11所示：圖2-11 仿真結(jié)果超調(diào)量為，調(diào)節(jié)時(shí)間，上升時(shí)間。2.4.3 仿真結(jié)果分析根據(jù)仿真結(jié)果，我們可以得到數(shù)字PID控制器控制的系統(tǒng)輸出，超調(diào)量為，調(diào)節(jié)時(shí)

12、間，上升時(shí)間。系統(tǒng)可快速達(dá)到穩(wěn)態(tài)，但是相對(duì)于模擬PID控制器的輸出，系統(tǒng)的快速性、穩(wěn)定性均有一定程度的降低。2.5調(diào)試Simulink仿真2.5.1 當(dāng)Kp變化時(shí)當(dāng)Kp=1時(shí)：圖2-12 當(dāng)Kp=1時(shí)的仿真結(jié)果超調(diào)量:27.5% 上升時(shí)間：1.02s 調(diào)節(jié)時(shí)間:0.85s，。當(dāng)Kp=8時(shí)：圖2-13 當(dāng)Kp=8時(shí)的仿真結(jié)果超調(diào)量：88.7% 上升時(shí)間：1.017s 調(diào)節(jié)時(shí)間：4.5s當(dāng)Kp=15時(shí)：圖2-14 當(dāng)Kp=15時(shí)的仿真結(jié)果此時(shí)系統(tǒng)已不穩(wěn)定結(jié)論：比例系數(shù)只改變系統(tǒng)的增益，對(duì)系統(tǒng)的影響主要反映在系統(tǒng)的穩(wěn)態(tài)誤差和穩(wěn)定性上。比例系數(shù)KP太小，系統(tǒng)動(dòng)作緩慢，增大比例系數(shù)可提高系統(tǒng)的開環(huán)增益，

13、減小系統(tǒng)的穩(wěn)態(tài)誤差，提高響應(yīng)速度但會(huì)降低系統(tǒng)的相對(duì)穩(wěn)定性，過大的比例系數(shù)會(huì)使系統(tǒng)有比較大的超調(diào)，并產(chǎn)生震蕩，使穩(wěn)定性變差。2.5.2 當(dāng)Ti變化，即Ki變化時(shí)當(dāng)Ki=1時(shí)：圖2-15 當(dāng)Ki=1時(shí)的仿真結(jié)果超調(diào)量：45.75% 上升時(shí)間：0.017s 調(diào)節(jié)時(shí)間：9.997s當(dāng)Ki=10時(shí)：圖2-16 當(dāng)Ki=10時(shí)的仿真結(jié)果超調(diào)量：45.29% 上升時(shí)間：0.0165s 調(diào)節(jié)時(shí)間：2.9s當(dāng)Ki=50時(shí)：圖2-17 當(dāng)Ki=60時(shí)的仿真結(jié)果超調(diào)量：48.86% 上升時(shí)間：0.0171s 調(diào)節(jié)時(shí)間：0.5s當(dāng)Ki=100時(shí)：圖2-18 當(dāng)Ki=100時(shí)的仿真結(jié)果超調(diào)量：56.85% 上升時(shí)間：0

14、.017s 調(diào)節(jié)時(shí)間：0.35s當(dāng)Ki=200時(shí)：圖2-19 當(dāng)Ki=200時(shí)的仿真結(jié)果超調(diào)量：69.93% 上升時(shí)間：0.0175s 調(diào)節(jié)時(shí)間：0.4s當(dāng)Ki=400時(shí)：圖2-20 當(dāng)Ki=400時(shí)的仿真結(jié)果超調(diào)量：96.81% 上升時(shí)間：0.0175s 調(diào)節(jié)時(shí)間：0.4s當(dāng)Ki=500時(shí)：圖2-21 當(dāng)Ki=500時(shí)的仿真結(jié)果超調(diào)量：110% 上升時(shí)間：0.0175s 調(diào)節(jié)時(shí)間：0.4s結(jié)論：積分控制主要目的使系統(tǒng)無穩(wěn)態(tài)誤差。主要改善系統(tǒng)的穩(wěn)態(tài)性能。減小積分時(shí)間有利于減小超調(diào)，減小震蕩，使系統(tǒng)的穩(wěn)定性增加，但是系統(tǒng)靜差消除時(shí)間變長(zhǎng)。2.5.3 當(dāng)Td變化，即Kd變化時(shí)當(dāng)Kd=0.001時(shí)：

15、圖2-22 當(dāng)Kd=0.001時(shí)的仿真結(jié)果超調(diào)量：88.43% 上升時(shí)間：0.0278s 調(diào)節(jié)時(shí)間：2.21s當(dāng)Kd=0.01時(shí)：圖2-23 當(dāng)Kd=0.01時(shí)的仿真結(jié)果超調(diào)量：78.94% 上升時(shí)間：0.026s 調(diào)節(jié)時(shí)間：0.92s當(dāng)Kd=0.05時(shí)：圖2-24 當(dāng)Kd=0.05時(shí)的仿真結(jié)果超調(diào)量：60.43% 上升時(shí)間：0.0193s 調(diào)節(jié)時(shí)間：0.43s當(dāng)Kd=0.10時(shí)：圖2-25 當(dāng)Kd=0.10時(shí)的仿真結(jié)果超調(diào)量：62.02% 上升時(shí)間：0.016s 調(diào)節(jié)時(shí)間：0.41s當(dāng)Kd=0.2時(shí)圖2-26 當(dāng)Kd=0. 5時(shí)的仿真結(jié)果超調(diào)量：79.72% 上升時(shí)間：0.0125s 調(diào)節(jié)時(shí)

16、間：1.82s當(dāng)Kd=0.3時(shí)圖2-27 當(dāng)Kd=0. 3時(shí)的仿真結(jié)果結(jié)論：微分系數(shù)能預(yù)測(cè)誤差變化得趨勢(shì)，能抑制誤差的控制作用等于0，避免被控量的嚴(yán)重超調(diào)。主要改善系統(tǒng)的動(dòng)態(tài)性能。增大微分時(shí)間常數(shù)Td（即減小微分系數(shù)Kd=Kp/Kd）有利于加快系統(tǒng)的響應(yīng)速度，使系統(tǒng)超調(diào)量減小，穩(wěn)定性增加，但系統(tǒng)對(duì)擾動(dòng)的抑制能力減弱。微分系數(shù)Kd對(duì)系統(tǒng)超調(diào)的抑制作用非常明顯。2.6設(shè)計(jì)總結(jié)P、I、D控制是經(jīng)典常用的控制方法，在現(xiàn)實(shí)中有著廣泛的應(yīng)用。了解P、I、D參數(shù)的改變對(duì)系統(tǒng)控制效果的影響十分重要。PID控制的結(jié)構(gòu)比較簡(jiǎn)單，但三個(gè)系數(shù)有著比較明顯的意義：比例控制其直接響應(yīng)與當(dāng)前的誤差信號(hào)，一旦發(fā)生誤差信號(hào)，則

17、控制其立即發(fā)生作用以減少偏差，Kp的值增大則偏差將減小，然而這不是絕對(duì)的，考慮根軌跡分析，Kp無限增大會(huì)使閉環(huán)系統(tǒng)不穩(wěn)定。積分控制器對(duì)以往的誤差信號(hào)發(fā)生作用，引入積分環(huán)節(jié)能消除控制中的靜態(tài)誤差，但Ki的值增大可能增加系統(tǒng)的超調(diào)量。積分作用太強(qiáng)也會(huì)引起振蕩，太弱會(huì)使系統(tǒng)存在余差。微分控制對(duì)誤差的導(dǎo)數(shù)，即變化率發(fā)生作用，有定的預(yù)報(bào)功能，有超前調(diào)節(jié)的作用，能在誤差有大的變化趨勢(shì)時(shí)施加合適的控制，Kd的值增大能加快系統(tǒng)的響應(yīng)速度，減少調(diào)節(jié)時(shí)間。對(duì)滯后大的對(duì)象有很好的效果。但不能克服純滯后，使用微分調(diào)節(jié)可使系統(tǒng)收斂周期的時(shí)間縮短。微分時(shí)間太長(zhǎng)也會(huì)引起振蕩。三、大林算法計(jì)算機(jī)控制系統(tǒng)設(shè)計(jì)及仿真實(shí)現(xiàn)3.1

18、設(shè)計(jì)要求已知被控對(duì)象的傳遞函數(shù)為：采樣周期為T=0.5s，用大林算法設(shè)計(jì)數(shù)字控制器D(z)，并分析是否會(huì)產(chǎn)生振鈴現(xiàn)象。要求：1、用大林算法設(shè)計(jì)數(shù)字控制器D(z)；2、在 Simulink 仿真環(huán)境畫出仿真框圖及得出仿真結(jié)果，畫出數(shù)字控制；3、繪制并分析數(shù)字控制器的振鈴現(xiàn)象；4、對(duì)振鈴現(xiàn)象進(jìn)行消除；5、得出仿真結(jié)果并進(jìn)行仿真分析；6、程序清單及簡(jiǎn)要說明；7、撰寫設(shè)計(jì)報(bào)告（列出參考文獻(xiàn)，以及仿真結(jié)果及分析）。3.2 總體方案設(shè)計(jì)圖3-1 系統(tǒng)控制框圖3.3 基本理論知識(shí)3.3.1大林控制算法的設(shè)計(jì)目標(biāo)大林控制算法的設(shè)計(jì)目標(biāo)是使整個(gè)閉環(huán)系統(tǒng)所期望的傳遞函數(shù)相當(dāng)于一個(gè)延遲環(huán)節(jié)和一個(gè)慣性環(huán)節(jié)相串聯(lián)，即：

19、整個(gè)閉環(huán)系統(tǒng)的純滯后時(shí)間和被控對(duì)象的純滯后時(shí)間相同。閉環(huán)系統(tǒng)的時(shí)間常數(shù)為，純滯后時(shí)間y與采樣周期T有整數(shù)倍的關(guān)系，。3.3.2大林算法的D(z)基本形式帶有純滯后的一階慣性環(huán)節(jié)：其與零階保持器相串聯(lián)的脈沖傳遞函數(shù)為：于是相應(yīng)的控制器形式為：3.3.3振鈴現(xiàn)象及其抑制（1）定義：控制量以1/2的采樣頻率（即二倍采樣周期）振蕩的現(xiàn)象稱為“振鈴”。這種振蕩一般是衰減的。（2）產(chǎn)生原因：如果在U(z)的脈沖傳遞函數(shù)表達(dá)式中，包含有在z 平面單位圓內(nèi)接近-1的實(shí)數(shù)極點(diǎn)，則會(huì)產(chǎn)生振鈴現(xiàn)象。 （3）解決辦法：令數(shù)字控制器中產(chǎn)生振鈴現(xiàn)象的極點(diǎn)（左半平面上接近-1的極點(diǎn)）的因子中的z=1，就可以消除振鈴現(xiàn)象。3

20、.3.4振鈴現(xiàn)象的消除被控對(duì)象為一階慣性環(huán)節(jié)。當(dāng)被控對(duì)象為純滯后一階慣性環(huán)節(jié)是，數(shù)字控制器D(z)為：由此可以得到振鈴幅度為：于是，如果選擇 ，則，無振鈴現(xiàn)象；如果選擇，則有振鈴現(xiàn)象。由此可見，當(dāng)系統(tǒng)的時(shí)間常數(shù)大于或是等于被控對(duì)象的時(shí)間常數(shù)時(shí)，即可消除振鈴現(xiàn)象。將D(z)的分母進(jìn)行分解可得：由上式，z=1處的極點(diǎn)不會(huì)引起振鈴現(xiàn)象?？赡芤鹫疋彫F(xiàn)象的因子為：當(dāng)N=0時(shí)，此因子消失，無振鈴可能。當(dāng)N=1時(shí)，有一個(gè)極點(diǎn)在，當(dāng)遠(yuǎn)小于T時(shí)，即當(dāng)遠(yuǎn)小于T時(shí)，將產(chǎn)生嚴(yán)重的振鈴現(xiàn)象。當(dāng)N=時(shí)，極點(diǎn)為當(dāng)遠(yuǎn)小于T時(shí)，于是，有嚴(yán)重的振鈴現(xiàn)象。以N=2時(shí)為例，數(shù)字控制器的形式為：當(dāng)遠(yuǎn)小于T時(shí)，有嚴(yán)重的振鈴現(xiàn)象，產(chǎn)生振

21、鈴現(xiàn)象的極點(diǎn)為式中的，于是令式中z=1，于是該因子變?yōu)椋汗氏疋彫F(xiàn)象后，D(z)的形式為：3.4 詳細(xì)設(shè)計(jì)3.4.1 設(shè)計(jì)大林控制器已知某控制系統(tǒng)被控對(duì)象的傳遞函數(shù)為則可得，連同零階保持器在內(nèi)的系統(tǒng)廣義被控對(duì)象的傳遞函數(shù)：其廣義對(duì)象的z傳遞函數(shù)：大林算法的設(shè)計(jì)目標(biāo)就是設(shè)計(jì)一個(gè)數(shù)字控制器，使整個(gè)閉環(huán)系統(tǒng)的z傳遞函數(shù)相當(dāng)于一個(gè)帶有純滯后的一階慣性環(huán)節(jié)，若，則可得：則大林控制器為：3.4.2 振鈴現(xiàn)象分析及消除由D(z)得到振鈴幅度為：由于RA0，則有振鈴現(xiàn)象。從D(z)得其三個(gè)極點(diǎn)為：根據(jù)判定結(jié)論，處的極點(diǎn)不會(huì)引起振鈴現(xiàn)象。所以，引起振鈴現(xiàn)象的極點(diǎn)為：依據(jù)大林消除振鈴現(xiàn)象的方法，應(yīng)去掉分母中的因

22、子，即令z=1，代入上式即可。此時(shí)，3.5 調(diào)試Simulink仿真3.5.1大林控制器設(shè)計(jì)通過前面的分析，我們得到大林控制器為：則Simulink仿真程序如圖3-2所示：圖3-2 仿真程序?yàn)榱蓑?yàn)證是否有振鈴現(xiàn)象，我們將控制器的輸入端加單位階躍信號(hào)，觀察其輸出控制量波形，如圖3-3所示：圖3-3 控制器在單位階躍輸入作用下的輸出控制量波形從上圖，我們觀察到，數(shù)字控制器的輸出以接近2T的周期大幅度上下擺動(dòng)。因此該控制器有振鈴現(xiàn)象，需要消除振鈴。系統(tǒng)的輸出曲線如圖3-4所示：圖3-4 輸出曲線3.5.2 消除振鈴現(xiàn)象根據(jù)前面分析的結(jié)果，為了消除振鈴，我們將控制器設(shè)置為：則Simulink仿真程序如

23、圖3-5所示：圖3-5 修正后的仿真程序此時(shí)我們觀察修正后的控制器在單位階躍輸入作用下的輸出控制量波形，如圖3-6所示，我們發(fā)現(xiàn)，控制器的輸出不再以接近2T的周期大幅度上下擺動(dòng)，振鈴現(xiàn)象被消除。圖3-6 修正后的控制器在單位階躍輸入作用下的輸出控制量波形修正后系統(tǒng)的輸出曲線如圖3-7所示：圖3-7 修正后的輸出曲線系統(tǒng)的超調(diào)量%=1.5%，調(diào)節(jié)時(shí)間Ts=7.14s。從本題我們得到以下結(jié)論，雖然振鈴現(xiàn)象已消除，但是系統(tǒng)的快速性有所降低，系統(tǒng)需要更長(zhǎng)的時(shí)間達(dá)到穩(wěn)態(tài)，而未消除振鈴之前，系統(tǒng)能夠迅速達(dá)到穩(wěn)態(tài)，但是由于存在振鈴，會(huì)降低設(shè)備的壽命，因此實(shí)際生產(chǎn)過程中，我們往往需要多方面考慮問題，避免顧此失

24、彼。四、二階彈簧阻尼系統(tǒng)的PID控制器設(shè)計(jì)及其參數(shù)整定4.1 設(shè)計(jì)要求考慮彈簧阻尼系統(tǒng)如圖4-1所示，其被控對(duì)象為二階環(huán)節(jié)，傳遞函數(shù)G(s)如下，參數(shù)為M=1kg，b=2N.s/m，k=25N/m，F(xiàn)(s)=1。圖4-1 彈簧阻尼系統(tǒng)示意圖彈簧阻尼系統(tǒng)的微分方程和傳遞函數(shù)為：圖4-2 閉環(huán)控制系統(tǒng)結(jié)構(gòu)圖要求完成：1、控制器為P控制器時(shí)，改變比例帶或比例系數(shù)大小，分析對(duì)系統(tǒng)性能的影響并繪制相應(yīng)曲線；2、控制器為PI控制器時(shí)，改變積分時(shí)間常數(shù)大小，分析對(duì)系統(tǒng)性能的影響并繪制相應(yīng)曲線；(當(dāng)kp=50時(shí)，改變積分時(shí)間常數(shù))3、設(shè)計(jì)PID控制器，選定合適的控制器參數(shù)，使階躍響應(yīng)曲線的超調(diào)量%20%,過渡

25、過程時(shí)間ts2s, 并繪制相應(yīng)曲線。4.2 總體方案設(shè)計(jì)圖4-3閉環(huán)控制系統(tǒng)結(jié)構(gòu)圖將控制器分別設(shè)計(jì)成P、PI、PID三種控制器，并比較觀察其各自特點(diǎn)。4.3 基本理論知識(shí)4.3.1比例（P）控制比例(P)控制是一種最簡(jiǎn)單的控制方式，其控制器的輸出與輸入誤差信號(hào)成比例關(guān)系。當(dāng)僅有比例控制時(shí)系統(tǒng)輸出存在穩(wěn)定誤差。比例控制器的傳遞函數(shù)為：式中，Kp稱為比例系數(shù)或增益（視情況可設(shè)置為正或負(fù)），一些傳統(tǒng)的控制器又常用比例帶（Proportional Band，PB），來取代比例系數(shù)Kp，比例帶是比例系數(shù)的倒數(shù)，比例帶也稱為比例度。對(duì)于單位反饋系統(tǒng)，0型系統(tǒng)響應(yīng)實(shí)際階躍信號(hào)的穩(wěn)態(tài)誤差與其開環(huán)增益K近視成反

26、比，即：對(duì)于單位反饋系統(tǒng)，I型系統(tǒng)響應(yīng)勻速信號(hào)的穩(wěn)態(tài)誤差與其開環(huán)增益Kv近視成反比，即：P控制只改變系統(tǒng)的增益而不影響相位，它對(duì)系統(tǒng)的影響主要反映在系統(tǒng)的穩(wěn)態(tài)誤差和穩(wěn)定性上，增大比例系數(shù)可提高系統(tǒng)的開環(huán)增益，減小系統(tǒng)的穩(wěn)態(tài)誤差，從而提高系統(tǒng)的控制精度，但這會(huì)降低系統(tǒng)的相對(duì)穩(wěn)定性，甚至可能造成閉環(huán)系統(tǒng)的不穩(wěn)定，因此，在系統(tǒng)校正和設(shè)計(jì)中P控制一般不單獨(dú)使用。4.3.2 比例積分(PI)控制比例積分(PI)控制具有比例加積分控制規(guī)律的控制稱為比例積分控制器，即PI控制，PI控制的傳遞函數(shù)為：其中，Kp 為比例系數(shù)，Ti稱為積分時(shí)間常數(shù)，兩者都是可調(diào)的參數(shù).控制器的輸出信號(hào)為：PI控制器可以使系統(tǒng)在進(jìn)

27、入穩(wěn)態(tài)后無穩(wěn)態(tài)誤差。PI控制器在與被控對(duì)象串聯(lián)時(shí)，相當(dāng)于在系統(tǒng)中增加了一個(gè)位于原點(diǎn)的開環(huán)極點(diǎn)，同時(shí)也增加了一個(gè)位于s左半平面的開環(huán)零點(diǎn)。位于原點(diǎn)的極點(diǎn)可以提高系統(tǒng)的型別，以消除或減小系統(tǒng)的穩(wěn)態(tài)誤差,改善系統(tǒng)的穩(wěn)態(tài)性能，而增加的負(fù)實(shí)部零點(diǎn)則可減小系統(tǒng)的阻尼程度，緩和PI控制器極點(diǎn)對(duì)系統(tǒng)穩(wěn)定性及動(dòng)態(tài)過程產(chǎn)生的不利影響.在實(shí)際工程中，PI控制器通常用來改善系統(tǒng)的穩(wěn)定性能。4.3.3 比例積分微分(PID)控制具有比例+積分+微分控制規(guī)律的控制稱為比例積分微分控制，即PID控制，PID控制的傳遞函數(shù)為：其中， Kp 為比例系數(shù)，Ti為微分時(shí)間常數(shù)，為微分時(shí)間常數(shù)，三者都是可調(diào)的參數(shù)。PID控制器的輸出

28、信號(hào)為：PID控制器的傳遞函數(shù)可寫成：PI控制器與被控對(duì)象串聯(lián)連接時(shí)，可以使系統(tǒng)的型別提高一級(jí)，而且還提供了兩個(gè)負(fù)實(shí)部的零點(diǎn)。與PI控制器相比，PID控制器除了同樣具有提高系統(tǒng)穩(wěn)定性能的優(yōu)點(diǎn)外，還多提供了一個(gè)負(fù)實(shí)部零點(diǎn)。因此在提高系統(tǒng)動(dòng)態(tài)系統(tǒng)方面提供了很大的優(yōu)越性。在實(shí)際過程中，PID控制器被廣泛應(yīng)用。PID控制通過積分作用消除誤差，而微分控制可縮小超調(diào)量，加快反應(yīng)，是綜合了PI控制與PD控制長(zhǎng)處并去除其短處的控制。從頻域角度看，PID控制通過積分作用于系統(tǒng)的低頻段，以提高系統(tǒng)的穩(wěn)定性，而微分作用于系統(tǒng)的中頻段，以改善系統(tǒng)的動(dòng)態(tài)性能。PID控制器的參數(shù)整定是控制系統(tǒng)設(shè)計(jì)的核心內(nèi)容。它是根據(jù)被控

29、過程的特性確定PID控制器的比例系數(shù)、積分時(shí)間和微分時(shí)間的大小。PID控制器參數(shù)整定的方法很多，概括起來有兩大類：一是理論計(jì)算整定法。它主要是依據(jù)系統(tǒng)的數(shù)學(xué)模型，經(jīng)過理論計(jì)算確定控制器參數(shù)。這種方法所得到的計(jì)算數(shù)據(jù)未必可以直接用，還必須通過工程實(shí)際進(jìn)行調(diào)整和修改。二是工程整定方法，它主要依賴工程經(jīng)驗(yàn)，直接在控制系統(tǒng)的試驗(yàn)中進(jìn)行，且方法簡(jiǎn)單、易于掌握，在工程實(shí)際中被廣泛采用。PID控制器參數(shù)的工程整定方法，主要有臨界比例法、反應(yīng)曲線法和衰減法。三種方法各有其特點(diǎn)，其共同點(diǎn)都是通過試驗(yàn)，然后按照工程經(jīng)驗(yàn)公式對(duì)控制器參數(shù)進(jìn)行整定。但無論采用哪一種方法所得到的控制器參數(shù)，都需要在實(shí)際運(yùn)行中進(jìn)行最后調(diào)整

30、與完善。現(xiàn)在一般采用的是臨界比例法。利用該方法進(jìn)行 PID控制器參數(shù)的整定步驟如下：(1)首先預(yù)選擇一個(gè)足夠短的采樣周期讓系統(tǒng)工作；(2)僅加入比例控制環(huán)節(jié)，直到系統(tǒng)對(duì)輸入的階躍響應(yīng)出現(xiàn)臨界振蕩，記下這時(shí)的比例放大系數(shù)和臨界振蕩周期；(3)在一定的控制度下通過公式計(jì)算得到PID控制器的參數(shù)。 4.4 詳細(xì)設(shè)計(jì)4.4.1 P控制器P控制器的傳遞函數(shù)為：系統(tǒng)傳遞函數(shù)：用matlab作出Kp分別為20、40、60、80的階躍響應(yīng)曲線：z=20:20:80;i=1;hold onfor kp=zsubplot(2,2,i);num=kp;den=1,2,25+kp;step(num,den); i=i

31、+1; %依次為Kp等于20,40,60,80的階躍響應(yīng)圖endtitle(Step Response);hold off仿真圖如下：圖4-4 依次為Kp=20、40、60、80的情況由圖可看出，隨著Kp加大，曲線圖波動(dòng)變大，穩(wěn)定性變差，準(zhǔn)確性變好，系統(tǒng)響應(yīng)時(shí)間加快。系統(tǒng)的超調(diào)量%和調(diào)節(jié)時(shí)間Ts：Kp=20，%=62.2%，Ts=3.86s；Kp=40，%=67.5%，Ts=3.83s；Kp=60，%=71.0%，Ts=3.76s；Kp=80，%=73.5%，Ts=3.64s；4.4.2 PI控制器PI控制器的傳遞函數(shù)為： Simulink建立的仿真模型為：圖4-5 PI控制的仿真建模圖4-6

32、 在Kp=50的情況下Ki為20的情況圖4-7 在Kp=50的情況下Ki為40的情況圖4-8 在Kp=50的情況下Ki為60的情況圖4-9 在Kp=50的情況下Ki為80的情況系統(tǒng)傳遞函數(shù)： 用matlab作出Kp=50時(shí)，Ki分別為20、40、60、80的階躍響應(yīng)曲線：kp=50;z=20:20:80;i=1;hold onfor ki=zsubplot(2,2,i); num=kp, ki;den=1,2,25+kp,ki; step(num,den);i=i+1; % 依次為Ki等于20,40,60,80的階躍響應(yīng)圖endtitle(Step Response);hold off圖4-1

33、0 在Kp=50的情況下Ki分別為20、40、60、80的情況由圖可看出，在積分控制中，控制器的輸出與輸入誤差信號(hào)的積分成正比關(guān)系。積分項(xiàng)對(duì)誤差取決于時(shí)間的積分，隨著時(shí)間的增加，積分項(xiàng)會(huì)增大。這樣，即便誤差很小，積分項(xiàng)也會(huì)隨著時(shí)間的增加而加大，它推動(dòng)控制器的輸出增大使穩(wěn)態(tài)誤差進(jìn)一步減小，直到等于零。因此，PI控制器可以使系統(tǒng)在進(jìn)入穩(wěn)態(tài)后無穩(wěn)態(tài)誤差。幾個(gè)曲線圖比較發(fā)現(xiàn)后，Ki變大,波形圖略微變化，穩(wěn)定性變差，準(zhǔn)確性變好，響應(yīng)時(shí)間稍微變快。系統(tǒng)的超調(diào)量%和調(diào)節(jié)時(shí)間Ts：Ki=20，%=18.7%，Ts=10.5s；Ki=40，%=24.4%，Ts=7.39s；Ki=60，%=29.9%，Ts=5.

34、95s；Ki=80，%=35.4%，Ts=5.93s；4.4.3 PID控制器PID控制器的傳遞函數(shù)為： Simulink建立的仿真模型為：圖4-11 PID控制的仿真建模經(jīng)過多輪調(diào)試當(dāng)Kp=20，Ki=70，Kd=10滿足使階躍響應(yīng)曲線的超調(diào)量%20%，過渡過程時(shí)間Ts2s。圖4-12 PID控制仿真結(jié)果系統(tǒng)的超調(diào)量%=3.4999%，調(diào)節(jié)時(shí)間Ts=1.7421s。4.5 分析與總結(jié)（1）P控制器只改變系統(tǒng)的增益而不影響相位，它對(duì)系統(tǒng)的影響主要反映在系統(tǒng)的穩(wěn)態(tài)誤差和穩(wěn)定性上，增大比例系數(shù)可提高系統(tǒng)的開環(huán)增益，減小系統(tǒng)的穩(wěn)態(tài)誤差，從而提高系統(tǒng)的控制精度，但這會(huì)降低系統(tǒng)的相對(duì)穩(wěn)定性，甚至可能造成

35、閉環(huán)系統(tǒng)的不穩(wěn)定。（2）PI控制器消除或減小系統(tǒng)的穩(wěn)態(tài)誤差，改善系統(tǒng)的穩(wěn)態(tài)性能。（3）PID控制通過積分作用消除誤差，而微分控制可縮小超調(diào)量，加快反應(yīng),是綜合了PI控制與PD控制長(zhǎng)處并去除其短處的控制。從頻域角度看，PID控制通過積分作用于系統(tǒng)的低頻段，以提高系統(tǒng)的穩(wěn)定性，而微分作用于系統(tǒng)的中頻段，以改善系統(tǒng)的動(dòng)態(tài)性能。五、二階系統(tǒng)串聯(lián)校正裝置的設(shè)計(jì)與分析5.1 設(shè)計(jì)要求設(shè)某被控系統(tǒng)的傳遞函數(shù)G(s)如下：選用合適的方法設(shè)計(jì)一個(gè)串聯(lián)校正裝置K(s)，使閉環(huán)系統(tǒng)的階躍響應(yīng)曲線超調(diào)量，過渡過程時(shí)間，開環(huán)比例系數(shù)，并分析串聯(lián)校正裝置中增益、極點(diǎn)和零點(diǎn)對(duì)系統(tǒng)性能的影響。5.2 基本理論知識(shí)基于系統(tǒng)頻率

36、特性的串聯(lián)校正設(shè)計(jì)方法在工程中被大量采用，主要分為相位超前、相位滯后以及相位滯后超前校正設(shè)計(jì)等三種方法。相位超前校正主要用于改善閉環(huán)系統(tǒng)的動(dòng)態(tài)特性，對(duì)于系統(tǒng)的穩(wěn)態(tài)精度影響較小。相位滯后校正可以明顯地改善系統(tǒng)的穩(wěn)態(tài)性能，但會(huì)使動(dòng)態(tài)響應(yīng)過程變緩；而相位滯后超前校正則把兩者的特性結(jié)合起來，用于動(dòng)態(tài)、靜態(tài)特性要求較高的系統(tǒng)。如果性能指標(biāo)以單位階躍響應(yīng)的峰值時(shí)間、調(diào)整時(shí)間、超調(diào)量、阻尼系統(tǒng)、穩(wěn)態(tài)誤差等時(shí)域特征量給出時(shí)，一般采用根軌跡法校正。根軌跡法校正的基本思路為借助根軌跡曲線進(jìn)行校正。如果系統(tǒng)的期望主導(dǎo)極點(diǎn)往往不在系統(tǒng)的根軌跡上，由根軌跡的理論，添加上開環(huán)零點(diǎn)或極點(diǎn)可以使根軌跡曲線形狀改變。若期望主導(dǎo)

37、極點(diǎn)在原根軌跡的左側(cè)，則只要加上一對(duì)零、極點(diǎn)，使零點(diǎn)位置位于極點(diǎn)右側(cè)。如果適當(dāng)選擇零、極點(diǎn)的位置，就能夠使系統(tǒng)根軌跡通過期望主導(dǎo)極點(diǎn)s1，并且使主導(dǎo)極點(diǎn)在s1點(diǎn)位置時(shí)的穩(wěn)態(tài)增益滿足要求。此即相當(dāng)于相位超前校正。如果系統(tǒng)的期望主導(dǎo)極點(diǎn)若在系統(tǒng)的根軌跡上，但是在該點(diǎn)的靜態(tài)特性不滿足要求，即對(duì)應(yīng)的系統(tǒng)開環(huán)增益K太小。單純?cè)龃驥值將會(huì)使系統(tǒng)阻尼比變小，甚至于使閉環(huán)特征根跑到復(fù)平面S的右半平面去。為了使閉環(huán)主導(dǎo)極點(diǎn)在原位置不動(dòng)，并滿足靜態(tài)指標(biāo)要求，則可以添加上一對(duì)偶極子，其極點(diǎn)在其零點(diǎn)的右側(cè)。從而使系統(tǒng)原根軌跡形狀基本不變，而在期望主導(dǎo)極點(diǎn)處的穩(wěn)態(tài)增益得到加大。此即相當(dāng)于相位滯后校正。根軌跡超前校正的主

38、要步驟為：（1）依據(jù)要求的系統(tǒng)性能指標(biāo)，求出主導(dǎo)極點(diǎn)的期望位置。（2）觀察期望的主導(dǎo)極點(diǎn)是否位于校正前的系統(tǒng)根軌跡上。（3）如果需要設(shè)計(jì)校正網(wǎng)絡(luò)，設(shè)計(jì)校正網(wǎng)絡(luò)。（4）校正網(wǎng)絡(luò)零點(diǎn)的確定?？芍苯釉谄谕拈]環(huán)極點(diǎn)位置下方（或在頭兩個(gè)實(shí)極點(diǎn)的左側(cè)）增加一個(gè)相位超前網(wǎng)絡(luò)的實(shí)零點(diǎn)。（5）校正網(wǎng)絡(luò)極點(diǎn)的確定。確定校正網(wǎng)絡(luò)極點(diǎn)的位置，使期望的主導(dǎo)極點(diǎn)位于校正后的根軌跡上。利用校正網(wǎng)絡(luò)極點(diǎn)的相角，使得系統(tǒng)在期望主導(dǎo)極點(diǎn)上滿足根軌跡的相角條件。（6）估計(jì)在期望的閉環(huán)主導(dǎo)極點(diǎn)處的總的系統(tǒng)開環(huán)增益。計(jì)算穩(wěn)態(tài)誤差系數(shù)。如果穩(wěn)態(tài)誤差系數(shù)不滿足要求，重復(fù)不述步驟。（7）利用根軌跡設(shè)計(jì)相位超前網(wǎng)絡(luò)時(shí)，超前網(wǎng)絡(luò)的傳遞函數(shù)可表

39、示為： 其中|z| num=20;den=1,2,0; Gs=tf(num,den); G=feedback(Gs,1);在控制器C取值為常數(shù)1的情況下，繪制此單位負(fù)反饋線性系統(tǒng)的根軌跡圖、系統(tǒng)的伯德圖以及閉環(huán)階躍響應(yīng)曲線： sisotool(Gs) step(G)圖5-1 校正前根軌跡圖、系統(tǒng)的伯德圖以及閉環(huán)階躍響應(yīng)曲線圖5-2 校正前的階躍響應(yīng)曲線由階躍響應(yīng)曲線可以看到，此時(shí)在沒有串聯(lián)校正裝置情況下，超調(diào)量為48.5% 20%，過渡過程時(shí)間Ts=3.78s1.5s，達(dá)不到指標(biāo)要求。5.3.2 設(shè)計(jì)根軌跡校正器在SISO Design Tool 窗口中，對(duì)控制器C進(jìn)行設(shè)置。利用增加和刪除零極

40、點(diǎn)的設(shè)置菜單，對(duì)控制器C的零極點(diǎn)任意設(shè)置。同時(shí)對(duì)控制器C的增益進(jìn)行設(shè)置。在系統(tǒng)設(shè)計(jì)完成后，需要對(duì)其做進(jìn)一步分析。在matlab中分析校正后閉環(huán)系統(tǒng)的階躍響應(yīng)，計(jì)算超調(diào)量、調(diào)整時(shí)間等指標(biāo)，以確保系統(tǒng)滿足設(shè)計(jì)要求。若不滿足設(shè)計(jì)要求，則重復(fù)上述步驟，不斷調(diào)整Gc的零極點(diǎn)，最終達(dá)到設(shè)計(jì)要求。重復(fù)執(zhí)行多次上述步驟后，可得出串聯(lián)校正裝置中增益、極點(diǎn)和零點(diǎn)對(duì)系統(tǒng)性能產(chǎn)生一定的影響。（1）校正裝置中增益對(duì)系統(tǒng)性能的影響：可以改變開環(huán)增益的大小，從而改善穩(wěn)態(tài)誤差。（2）校正裝置中極點(diǎn)對(duì)系統(tǒng)性能的影響：增加開環(huán)極點(diǎn)，使得原系統(tǒng)根軌跡的整體走向在S平面向右移，使系統(tǒng)穩(wěn)定性變壞。（3）校正裝置中零點(diǎn)對(duì)系統(tǒng)性能的影響：

41、增加開環(huán)零點(diǎn)，使得原系統(tǒng)根軌跡的整體走向在S平面向右移，使系統(tǒng)穩(wěn)定性得到改善。經(jīng)過多次調(diào)節(jié)后，最終得到理想結(jié)果。圖5-3 得到理想效果時(shí)的零極點(diǎn)配置5.3.3 分析校正后系統(tǒng)的穩(wěn)態(tài)性能Matlab如下： numc=0.25,1;denc=0.067,1; Gc=tf(numc,denc) G=Gc*Gs; G=feedback(G,1)Transfer function: 5 s + 20-0.04 s3 + 1.08 s2 + 7 s + 20 step(G)圖5-4 校正后的階躍響應(yīng)曲線經(jīng)計(jì)算，超調(diào)量12.1% 20%，過渡過程時(shí)間Ts = 1.04s %系統(tǒng)狀態(tài)方程A=0 1 0 0;0

42、 -0.1818 2.6727 0;0 0 0 1;0 -0.4545 31.1818 0;B=0; 1.8182;0; 4.5455;C=1 0 0 0;0 0 1 0;D=0;%求解系統(tǒng)的特征值p=eig(A)t=0:0.01:1;%階躍響應(yīng)step(A,B,C,D,1,t)系統(tǒng)的特征值p = 0-5.6041-0.14285.5651系統(tǒng)有一個(gè)位于右半平面的極點(diǎn)，故不穩(wěn)定。系統(tǒng)的階躍響應(yīng)曲線如圖6-3所示，上圖是小車坐標(biāo)x的階躍響應(yīng)曲線，下圖是倒立擺的垂直角度的階躍響應(yīng)曲線，所以必須加入校正裝置。圖6-3 開環(huán)系統(tǒng)階躍響應(yīng)曲線（2）線性二次型最優(yōu)控制器的設(shè)計(jì)設(shè)計(jì)線性二次型最優(yōu)控制器的設(shè)計(jì)

43、的關(guān)鍵是選擇加權(quán)矩陣Q。一般來說，Q選擇的越大，系統(tǒng)達(dá)到穩(wěn)定所需要的時(shí)間越短。首先選擇，R=1，然后根據(jù)實(shí)際情況調(diào)節(jié)。編程如下：%Q和R矩陣的選擇m=1;n=1Q=m 0 0 0;0 0 0 0;0 0 n 0;0 0 0 0;R=1;%求解線性二次型最優(yōu)狀態(tài)k,p,e=lqr(A,B,Q,R)%求解系統(tǒng)閉環(huán)狀態(tài)方程Ac=(A-B*k);Bc=B;Cc=C;Dc=D;%輸出系統(tǒng)階躍響應(yīng)T=0:0.02:10;U=ones(size(T);Y,X=lsim(Ac,Bc,Cc,Dc,U,T);y1=Y(:,1);y2=Y(:,2);plot(T,y1,r:,T,y2);title(Inverted

44、 Pendulum LQ Step Response);xlabel(時(shí)間（sec）);ylabel(響應(yīng));grid;legend(小車位移,擺桿傾角)結(jié)果為：k = -1.0000 -1.6567 18.6852 3.4594p = 1.5567 1.2067 -3.4594 -0.7027 1.2067 1.4554 -4.6826 -0.9466 -3.4594 -4.6826 31.6315 5.9838 -0.7027 -0.9466 5.9838 1.1397e = -0.8494 + 0.8323i -0.8494 - 0.8323i -5.5978 + 0.4070i -5.5978 - 0.4070i此時(shí)求得的線性二次型最優(yōu)狀態(tài)反饋矩陣為：其響應(yīng)曲線如圖6-4所示：圖6-4 m=1、n=1時(shí)階躍響應(yīng)曲線從圖6-4可以看出，超調(diào)量基本滿足要求，但穩(wěn)定值與系統(tǒng)期望值相差太大（小車坐標(biāo)的響應(yīng)曲線穩(wěn)態(tài)值為負(fù)值）；另一方