自動(dòng)控制原理試驗(yàn)時(shí)域分析

1、實(shí)驗(yàn)報(bào)告課程名稱(chēng)：自動(dòng)控制原理課題名稱(chēng)：實(shí)驗(yàn)三線(xiàn)性系統(tǒng)時(shí)域分析-1專(zhuān)業(yè)班級(jí)：14通信工程 一班實(shí)險(xiǎn)成員：呂傳磊（）一、實(shí)驗(yàn)工具1、PC2、MATLAB 套件二、實(shí)驗(yàn)?zāi)康?、熟練掌握式印0函數(shù)，impulse。函數(shù)，Istm。函數(shù)和initial。函數(shù)的使用方法，研究線(xiàn)性系統(tǒng)在單位階躍、單位脈沖及任意激勵(lì)下的響應(yīng)，以及波形圖的繪制;2、通過(guò)響應(yīng)曲線(xiàn)觀測(cè)特征參量網(wǎng)二階系統(tǒng)動(dòng)態(tài)性能的影響;三、實(shí)驗(yàn)任務(wù)1.編程完成例2-12的單位脈沖響應(yīng)、階躍響應(yīng)、零輸入響應(yīng)與零狀態(tài)響應(yīng)的求解, 并畫(huà)出各響應(yīng)的波形圖。根據(jù)電路關(guān)系，列寫(xiě)方程，并做拉氏變換，得4 (s)0.1s + 0.2令Ui(s)=l河求得單位脈沖

2、響應(yīng)；令Ui(s)=l/s,可求得階躍響應(yīng)；零狀態(tài)響應(yīng)時(shí)，寫(xiě)出傳遞函數(shù)，與輸入相乘可得零輸入響應(yīng)；零輸入響應(yīng)可以利用初始條件，通過(guò)等效電壓電流源的方法，寫(xiě)出系統(tǒng)函數(shù), 從而求得響應(yīng)，也可利用階躍響應(yīng)零狀態(tài)響應(yīng)獲得。波形繪制可以利用impulse。函數(shù)實(shí)現(xiàn)，結(jié)果求解可以利用ilaplaceO函數(shù)，即拉 氏反變換完成。程序如下:%實(shí)驗(yàn)三例2-12% 201400121100 呂傳磊%14通信一班clc;clear;%單位脈沖響應(yīng),Ui=l;sysA=tf(0.1 1.2 JI 1 1) t=0:0.01:12;figure,inipulse(sysA,t/rr)Jiold on, grid on;

3、%階躍響應(yīng),Ui=l/s;sysl=tf(l,l 1 1 0);sys2=tf(0.1 0.2,l 1 1);sysB=parallel(sysl ,sys2) inipulse(sysBJiol d on;%零狀態(tài)響應(yīng)JI=1/sFZ+s+D;sysC=tf(l,l 1 1 0)inipulse(sysC,t/gf)Jiold on;%也可使用step求解，此時(shí)H=l/(8A2+s+l) %step(sysC,t/g%hold on;%零輸入響應(yīng)，利用全響應(yīng)零狀態(tài)響應(yīng) sysD=parallel(sysBrsysC) impulselsysD.tM)9bold on;%零輸入響應(yīng)，也可利用初

4、始條件與拉氏變換%H=(0.1 -0.1 /s) * (1)+0.1 /ssysl=tf(0,1,l 1 1);sys2=tf(0.1,l 1 1 0);sys3=tf(0.141 0);sysE=par all el (sysl ,sys2);sysE=parallel(sysE,sys3)impul seCsysE,1) ,hol d on;%驗(yàn)證零輸入響應(yīng)%與時(shí)域結(jié)果對(duì)比，驗(yàn)證結(jié)果%曲線(xiàn)重合即正確y=0.2.*exp(-0.5.*t).*sin(0.866.*t+pi./6);plot(t,y)Jiold on;%驗(yàn)證階躍響應(yīng)，即全響應(yīng)%與時(shí)域結(jié)果對(duì)比，驗(yàn)證結(jié)果%曲線(xiàn)重合即正確y=l+L

5、15.*exp(-0.5.*t).*sin(0.866.*t2.*pi./3)+0.2.*exp(-0.5.*t).*sin(0.866.*Hpi./6); plot(t,y)Jiold on;程序運(yùn)行結(jié)果如下：其中階躍響應(yīng)為兩曲線(xiàn)重合曲線(xiàn)，即與時(shí)域表達(dá)式波形進(jìn)行了重合驗(yàn)證；零輸入響應(yīng)為三曲線(xiàn)重合曲線(xiàn)，分別為（階躍零狀態(tài)響應(yīng)）、等效法求系統(tǒng)函數(shù)、 時(shí)域表達(dá)式波形。Dpn 七-dmv部分傳遞函數(shù)表達(dá)式如下：sysA =0. 1 s + 1. 2sA2 + s + 1Continuous-time transfer function.sysB =0. 1 sa4 +0.3 sA3 + 1. 3 s

6、*2 + 1.2 s + 1s5 + 2s4+3s3 + 2s2+s可以利用ilaplace。函數(shù)求解其時(shí)域表達(dá)式，部分求解如下: test030102 TOC o 1-5 h z / sqrt(3) t Isqrt(3) sin| | 23 |/ t |/ sqrt(3) t 2/|exp|I I cos| I + I2 / 2/3/10/ sqrt(3) t Isqrt (3) sin| | 7 |/ t |/ sqrt (3) t 2/|explI I cos|I +I 92 / 2/27/10可見(jiàn)，運(yùn)行正確。.改變系統(tǒng)阻尼比，畫(huà)出不同階躍響應(yīng)的波形圖，并比較阻尼比對(duì)峰值時(shí)間等系 統(tǒng)性能

7、的影響。傳遞函數(shù)表達(dá)式：(S)= 2 :7二 s2+8s+16利用tfO函數(shù)構(gòu)建系統(tǒng)模型，利用Step。求階躍響應(yīng)，并畫(huà)出圖形程序如下：%實(shí)驗(yàn)三例3-16%201400121100 呂傳磊%14通信一班clc;dear;zeta=O 0.1 0.3 0.5 0.707 0.9 1 2;t=0:0.01:12;for i=l:length(zeta)sys=tf(16,l 8*zeta 16)；%subpl ot(334)ane(ce112mat(sys.num),.%cell2mat(sys.den),title(阻尼比=，nun12str(zeta(i)；step(sys,t),hold o

8、n;grid on;end運(yùn)行結(jié)果如下:Step Response8Pm-dujaTime (seconds)系統(tǒng)零極點(diǎn)圖如下:阻尼比U)-202Real Pat 明尼比=0 707阿尼比二0 1 2 0.25.0G6OE-ia后比032024Real PM里尼比R54-20Real Part尼比11 O-14-3 ted XJ2meE-6440246 Real Part 明房比=0 9Real Part阻尼比二2i -101RW Par 1 。 5 J 510 0 4XJeBE-4.155 0-5pfld TeQeE-Real Part列下面表格，研究欠阻尼系統(tǒng)中不同阻尼比對(duì)系統(tǒng)性能的影響:阻尼比0.10.30.50.7070.9上升時(shí)間 tr(s)0.2760.330.4090.5370.721峰值時(shí)間 tp(s)0.790.820.911.111.8超調(diào)量72.9%37.2%16.3%4.33%0.152%調(diào)節(jié)時(shí)間 ts(s)9.62.812.021.491.17可由上表數(shù)據(jù)得出結(jié)論：在振蕩頻率叫一定的情況下，阻尼比越小，超調(diào)量越大，上升 時(shí)間越短，調(diào)節(jié)時(shí)間越長(zhǎng)，通常取0.40.8為宜