自動(dòng)控制原理課程設(shè)計(jì)基于MATLAB軟件的自動(dòng)控制系統(tǒng)仿真

1、 自動(dòng)控制原理 課程設(shè)計(jì)說(shuō) 明 書 日期： 2012.5.21-2012.6.1 自動(dòng)控制原理 課程設(shè)計(jì)評(píng)閱書題目基于MATLAB軟件的自動(dòng)控制系統(tǒng)仿真學(xué)生姓名學(xué)號(hào)指導(dǎo)教師評(píng)語(yǔ)及成績(jī)指導(dǎo)教師簽名： 年 月 日答辯評(píng)語(yǔ)及成績(jī)答辯教師簽名： 年 月 日教研室意見總成績(jī)： 教研室主任簽名： 年 月 日摘 要 MATLAB 是一種直觀、高效的計(jì)算機(jī)語(yǔ)言，同時(shí)也是一個(gè)科學(xué)計(jì)算平臺(tái)。MATLAB軟件技術(shù)進(jìn)行控制系統(tǒng)建模，時(shí)域、負(fù)域和頻域分析，前饋校正設(shè)計(jì)，離散系統(tǒng)與非線性系統(tǒng)分析，以及現(xiàn)代控制理論中的可控性、可觀性和穩(wěn)定性的判別。在MATLAB中，常用的系統(tǒng)建模方法有傳遞函數(shù)模型，零極點(diǎn)模型以及狀態(tài)空間模

2、型等。在這次課程設(shè)計(jì)中，就使用了MATLAB軟件對(duì)我們過(guò)去學(xué)過(guò)的一些簡(jiǎn)單的函數(shù)進(jìn)行仿真。此次學(xué)會(huì)了使用簡(jiǎn)單的MATLAB語(yǔ)言編寫簡(jiǎn)單的傳遞函數(shù)，還有使用Simulink進(jìn)行仿真。 我們可將綜合性和設(shè)計(jì)性實(shí)驗(yàn)項(xiàng)目通過(guò)MATLAB 在計(jì)算機(jī)上仿真，使系統(tǒng)的觀察實(shí)驗(yàn)的動(dòng)態(tài)過(guò)程。目前，MATLAB 已經(jīng)成為我們當(dāng)代大學(xué)生必須掌握的基本技能，在設(shè)計(jì)研究單位和工業(yè)部門，MATLAB 已經(jīng)成為研究和解決各種具體工程問(wèn)題的一種標(biāo)準(zhǔn)軟件。在完成了驗(yàn)證性、綜合性和設(shè)計(jì)性實(shí)驗(yàn)后，課程設(shè)計(jì)必不可少。課程設(shè)計(jì)是工科實(shí)踐教學(xué)的一個(gè)重要的環(huán)節(jié)，目的是培養(yǎng)我們綜合運(yùn)用理論知識(shí)分析和解決實(shí)際問(wèn)題的方法和能力,實(shí)現(xiàn)由知識(shí)向技能的

3、初步化。所以課程設(shè)計(jì)是培養(yǎng)我們思維創(chuàng)造能力最有效的圖徑。 關(guān)鍵詞： 控制系統(tǒng)；MATLAB;模擬仿真；思維創(chuàng)造目 錄1 課題描述12 仿真過(guò)程2 2.1 控制系統(tǒng)建模22.2 線性系統(tǒng)時(shí)域分析22.3 線性系統(tǒng)根軌跡分析42.4 線性系統(tǒng)頻域分析52.5 線性系統(tǒng)校正63 Simulink仿真7總結(jié)11參考文獻(xiàn)121 課題描述MATLAB 是一種直觀、高效的計(jì)算機(jī)語(yǔ)言，同時(shí)也是一個(gè)科學(xué)計(jì)算平臺(tái)。它的伴隨工具Simulink 是用來(lái)對(duì)真實(shí)世界的動(dòng)力學(xué)系統(tǒng)建模、模擬仿真和分析的軟件。我們可將綜合性和設(shè)計(jì)性實(shí)驗(yàn)項(xiàng)目通過(guò)MATLAB 在計(jì)算機(jī)上仿真，使系統(tǒng)的觀察實(shí)驗(yàn)的動(dòng)態(tài)過(guò)程。目前，MATLAB 已經(jīng)

4、成為我們當(dāng)代大學(xué)生必須掌握的基本技能，在設(shè)計(jì)研究單位和工業(yè)部門，MATLAB 已經(jīng)成為研究和解決各種具體工程問(wèn)題的一種標(biāo)準(zhǔn)軟件。在完成了驗(yàn)證性、綜合性和設(shè)計(jì)性實(shí)驗(yàn)后，課程設(shè)計(jì)必不可少。課程設(shè)計(jì)是工科實(shí)踐教學(xué)的一個(gè)重要的環(huán)節(jié)，目的是培養(yǎng)我們綜合運(yùn)用理論知識(shí)分析和解決實(shí)際問(wèn)題的方法和能力,實(shí)現(xiàn)由知識(shí)向技能的初步化。所以課程設(shè)計(jì)是培養(yǎng)我們思維創(chuàng)造能力最有效的途徑。 在學(xué)習(xí)自動(dòng)控制原理這個(gè)學(xué)科之后，這次課程設(shè)計(jì)又針對(duì)這個(gè)學(xué)科做了MATLAB軟件的仿真學(xué)習(xí)。這次課程設(shè)計(jì)的基本要求是學(xué)會(huì)使用MATLAB軟件并進(jìn)行仿真，學(xué)會(huì)建立簡(jiǎn)單的數(shù)學(xué)模型，例如，傳遞函數(shù)模型，零極點(diǎn)模型以及狀態(tài)空間模型等。然后對(duì)這個(gè)軟件

5、進(jìn)行一定的了解與認(rèn)識(shí)。而且還要學(xué)會(huì)使用簡(jiǎn)單的MARLAB語(yǔ)言和Simulink構(gòu)建簡(jiǎn)單的數(shù)學(xué)模型。通過(guò)學(xué)習(xí)培養(yǎng)理論聯(lián)系實(shí)際的科學(xué)態(tài)度，訓(xùn)練綜合運(yùn)用經(jīng)典控制理論和相關(guān)課程知識(shí)的能力。掌握自動(dòng)控制原理的時(shí)域分析法，根軌跡法，頻域分析法，以及各種（矯正）裝置的作用及用法，能夠利用不同的分析方法對(duì)給定系統(tǒng)進(jìn)行性能分析，能根據(jù)不同的系統(tǒng)性能指標(biāo)要求進(jìn)行合理的系統(tǒng)設(shè)計(jì)，并調(diào)試滿足系統(tǒng)的指標(biāo)。學(xué)會(huì)使用MATLAB語(yǔ)言及Simulink動(dòng)態(tài)仿真工具進(jìn)行系統(tǒng)的仿真與調(diào)試。鍛煉獨(dú)立思考和動(dòng)手解決控制系統(tǒng)實(shí)際問(wèn)題的能力。2 設(shè)計(jì)過(guò)程 2.1 控制系統(tǒng)建模已知多回路反饋系統(tǒng)的結(jié)構(gòu)如圖，求閉環(huán)系統(tǒng)的傳遞函數(shù)C(S)/R

6、（S),。解 MATLAB文本如下G1=tf(1,1 10);G2=tf(1,1 1);G3=tf(1 0 1,1 4 4);numg4=1 1;deng4=1 6;G4=tf(numg4,deng4);H1=zpk(-1,-2,1);numh2=2;denh2=1;H3=1;nh2=conv(numh2,deng4);dh2=conv(denh2,numg4);H2=tf(nh2,dh2);sys1=series(G3,G4);sys2=feedback(sys1,H1,+1);sys3=series(G2,sys2);sys4=feedback(sys3,H2);sys5=series(G

7、1,sys4);sys=feedback(sys5,H3)程序執(zhí)行結(jié)果Zero/pole/gain: 0.083333 (s+1)2 (s+2) (s2 + 1)-(s+10.12) (s+2.44) (s+2.349) (s+1) (s2 + 1.176s + 1.023)2.2 線性系統(tǒng)時(shí)域分析 已知系統(tǒng)的閉環(huán)傳遞函數(shù)為，其中&=0.707，求二階系統(tǒng)的單位脈沖響應(yīng)，單位脈沖響應(yīng)和單位斜坡響應(yīng)。解 MATLAB 文本如下zeta=0.707;num=16;den=1 8 * zeta 16;sys=tf(num,den);p=roots(den)t=0:0.01:3;figure(

8、1)impulse(sys,t);gridxlabel('t');ylabel('c(t)');title('impulse response');figure(2)step(sys,t);gridxlabel('t');ylabel('c(t)');title('step response');figure(3)u=t;lsim(sys,u,t,0);gridxlabel('t');ylabel('c(t)');title('ramp response

9、9;);程序執(zhí)行結(jié)果單位脈沖響應(yīng)單位階躍響應(yīng)單位斜坡響應(yīng)分析系統(tǒng)對(duì)輸入信號(hào)導(dǎo)數(shù)的響應(yīng)，就等于系統(tǒng)對(duì)該信號(hào)響應(yīng)的導(dǎo)數(shù)；或者，系統(tǒng)對(duì)輸入信號(hào)的積分的響應(yīng)，就等于系統(tǒng)對(duì)該輸入信號(hào)響應(yīng)的積分，而積分常數(shù)由零輸入初始條件確定。2.3 線性系統(tǒng)根軌跡分析已知單位負(fù)反饋系統(tǒng)的開環(huán)傳遞函數(shù)為試畫出K從零點(diǎn)變化到無(wú)窮時(shí)的根軌跡圖，并求出系統(tǒng)臨界阻尼時(shí)對(duì)應(yīng)的K值及其閉環(huán)極點(diǎn)。解 系統(tǒng)閉環(huán)特征多項(xiàng)式為D（s)=+4s+Ks+20=+4s+20+K(s+4)=0等效開環(huán)傳遞函數(shù)MATLAB文本如下G=TF(1 4,1 4 20);figure(1)pzmap(G);figure(2)rLocus(G);rlocfin

10、d(G)結(jié)果得到如下運(yùn)行結(jié)果零極點(diǎn)分布圖根軌跡圖分析 根軌跡的起點(diǎn)和終點(diǎn)，根軌跡起于開環(huán)極點(diǎn)，終于開環(huán)零點(diǎn)；根軌跡的分支數(shù)，對(duì)稱性和連續(xù)性，根軌跡的分支數(shù)與開環(huán)有限零點(diǎn)數(shù)m和有限極點(diǎn)數(shù)n中大者相等，它們是連續(xù)并且對(duì)稱于實(shí)軸；當(dāng)開環(huán)有限極點(diǎn)數(shù)n大于有限零點(diǎn)數(shù)m時(shí)，有n-m條根軌跡分支沿著與實(shí)軸交角為; k=0,1,2,.n-m-12.4 線性系統(tǒng)頻域分析已知單位負(fù)反饋系統(tǒng)的開環(huán)傳遞函數(shù)試?yán)L制其伯德圖和奈奎斯圖，并判斷閉環(huán)系統(tǒng)的穩(wěn)定性解 MATLAB 文本如下G=tf(1280 640,1 2402 1604.81 320.24 16);figure(1)margin(G);figure(2)ny

11、quist(G);axis equal運(yùn)行結(jié)果伯德圖奈奎斯特圖分析由于系統(tǒng)無(wú)右半邊開環(huán)極點(diǎn)，從圖可以看出，開環(huán)幅相曲線不包括（-1，j0)點(diǎn)，閉環(huán)系統(tǒng)穩(wěn)定，另外在伯德圖可以得出系統(tǒng)的幅值裕度h=29.5dB,相角裕度，相應(yīng)的截止頻率，穿越頻率。有奈式判知，系統(tǒng)閉環(huán)穩(wěn)定2.5 線性系統(tǒng)校正設(shè)單位負(fù)反饋系統(tǒng)的開環(huán)傳遞函數(shù)為若要求系統(tǒng)在單位斜坡輸入信號(hào)作用時(shí)，位置輸出穩(wěn)態(tài)誤差，開環(huán)系統(tǒng)截止頻率，相角裕度，幅值裕度，式設(shè)計(jì)串聯(lián)無(wú)緣超前網(wǎng)絡(luò)。解 1） 根據(jù)穩(wěn)態(tài)誤差要求，確定開環(huán)增益K。2） 利用以確定的開環(huán)增益，計(jì)算待校正系統(tǒng)的幅值裕度、相角裕度及其對(duì)應(yīng)的截止頻率、穿越頻率。3） 根據(jù)截止頻率的要求，計(jì)

12、算超前網(wǎng)絡(luò)參數(shù)a和T。為保證正系統(tǒng)的響應(yīng)速度，并充分利用網(wǎng)絡(luò)的橡膠超前特性，可選擇最大超前角頻率等于截止頻率，即，其中a由確定，然后再由，確定T值。4） 確定無(wú)源超前網(wǎng)絡(luò)和最大超前角；5） 驗(yàn)算已校正系統(tǒng)的幅值裕度，相角裕度以及對(duì)應(yīng)的截止頻率，穿越頻率。若驗(yàn)算結(jié)果不滿意指標(biāo)要求，需要新選擇，然后重復(fù)以上的設(shè)計(jì)步驟。 MTALAB 文本如下k=1/0.1;G0=zpk(,0 -1,k);h0,r,wx,ec=margin(G0)wm=4.4;L=bode(G0,wm);Lwc=20*log10(L)a=10(-0.1*Lwc)T=1/(wm*sqrt(a);phi=asin(a-1)/(a+1)

13、Gc=(1/a)*tf(a*T 1,T 1);Gc=a*Gc;G=Gc*G0;bode(G,'r',G0,'b-');grid;h,f,wx,wc=margin(G)運(yùn)行結(jié)果對(duì)數(shù)幅頻特性曲線分析圖中虛線部分為系統(tǒng)校正前的對(duì)數(shù)幅頻特性曲線，實(shí)線部分為系統(tǒng)校正后的對(duì)數(shù)幅頻特性曲線，如果存在誤差時(shí)，重新選取截止頻率，從而得到準(zhǔn)確的實(shí)驗(yàn)數(shù)據(jù)3 Simulink仿真已知系統(tǒng)的閉環(huán)傳遞函數(shù)為，其中&=0.707，求二階系統(tǒng)的單位脈沖響應(yīng)，單位脈沖響應(yīng)和單位斜坡響應(yīng)。解 MATLAB 文本如下zeta=0.707;num=16;den=1 8 * zeta 16;sy

14、s=tf(num,den);p=roots(den)t=0:0.01:3;figure(1)impulse(sys,t);gridxlabel('t');ylabel('c(t)');title('impulse response');figure(2)step(sys,t);gridxlabel('t');ylabel('c(t)');title('step response');figure(3)u=t;lsim(sys,u,t,0);gridxlabel('t');ylabel

15、('c(t)');title('ramp response');用Simulink仿真的環(huán)節(jié)模塊圖如圖。由于積分環(huán)節(jié)附帶的增益比較大（積分時(shí)間常數(shù)T0.1），Scope(示波器)繪出圖形的輻值顯示范圍并不是很理想。我們可以在Scope的顯示圖中點(diǎn)擊鼠標(biāo)右鍵，選Axes properties菜單，在彈出的對(duì)話框中設(shè)置Y-max屬性為100，則輸出結(jié)果如圖所示 總 結(jié) 本次課程設(shè)計(jì)，不僅是對(duì)前面所學(xué)知識(shí)的一種檢驗(yàn)，而且是對(duì)自己能力得一種提升。在設(shè)計(jì)過(guò)程中，我們通過(guò)到圖書管借閱了大量的自動(dòng)控制原理MATLAB實(shí)現(xiàn)的相關(guān)資料，自學(xué)的時(shí)候與同學(xué)、老師交流，一步一步的分析和

16、研究，最終完成了課程設(shè)計(jì)。在設(shè)計(jì)和分析過(guò)程中遇到不少問(wèn)題，首先是因?yàn)閷?duì)原先學(xué)過(guò)的知識(shí)有些遺忘，再次，MATLAB的初次見面也很陌生.不僅要復(fù)習(xí)自動(dòng)控制原理的相關(guān)知識(shí)，還要學(xué)習(xí)MATLAB軟件的使用。通過(guò)這次課程設(shè)計(jì)，我鞏固了自動(dòng)控制原理所學(xué)的基本知識(shí)。同時(shí)，最主要的是我對(duì)MATLAB軟件有了初步的了解，再以后的學(xué)習(xí)中可以輕松的解決MATLAB相關(guān)方面的問(wèn)題，也有助于自己的學(xué)習(xí)研究?？傊?，不論做好設(shè)計(jì)，還是學(xué)習(xí)，最主要的就是我們的態(tài)度，“努力造就實(shí)力，態(tài)度決定高度”。態(tài)度決定一切。初步使用MATLAB軟件時(shí)，不是很熟悉軟件的操作，出現(xiàn)了許多的問(wèn)題，雖然問(wèn)題不是很大，但影響了實(shí)驗(yàn)的準(zhǔn)確性，在輸入MATLAB文本時(shí)，疏漏了其中的標(biāo)點(diǎn)符號(hào)，甚至把函數(shù)的名字輸入錯(cuò)誤，倒是實(shí)驗(yàn)數(shù)據(jù)不準(zhǔn)確。當(dāng)?shù)貌坏綔?zhǔn)確的實(shí)驗(yàn)數(shù)據(jù)時(shí)，首先根絕軟件特有的功能，自己先進(jìn)行查詢，修改MATLAB文本中的錯(cuò)誤，自己不能查找解決的，通過(guò)老師和同學(xué)的幫助，與老師交流，一步一步的分析和研究，從而解決了其中所遇到的問(wèn)題，最終完成了其課程設(shè)計(jì)。 多鍛煉我們自己獨(dú)立思考和手解決控制系統(tǒng)實(shí)際問(wèn)題的能力。通過(guò)課程設(shè)計(jì)會(huì)讓我們學(xué)會(huì)