控制基礎(chǔ)第三次實(shí)驗(yàn)——MatlabSimulink仿真實(shí)驗(yàn)

1、東南大學(xué)自動(dòng)化學(xué)院自動(dòng)控制原理實(shí)驗(yàn)實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)五 Matlab/Simulink仿真實(shí)驗(yàn)姓 名： 學(xué) 號(hào)： 專 業(yè)： 自動(dòng)化 實(shí) 驗(yàn) 室： 組 別： 同組人員： 設(shè)計(jì)時(shí)間： 2014年12 月 12日評(píng)定成績(jī)： 審閱教師： 1、 實(shí)驗(yàn)?zāi)康?（1） 學(xué)習(xí)系統(tǒng)數(shù)學(xué)模型的多種表達(dá)方法，并會(huì)用函數(shù)相互轉(zhuǎn)換。（2） 學(xué)習(xí)模型串并聯(lián)及反饋連接后的系統(tǒng)傳遞函數(shù)。（3） 掌握系統(tǒng)BODE圖，根軌跡圖及奈奎斯特曲線的繪制方法。并利用其對(duì)系統(tǒng)進(jìn)行分析。（4） 掌握系統(tǒng)時(shí)域仿真的方法，并利用其對(duì)系統(tǒng)進(jìn)行分析。2、 預(yù)習(xí)與回答 借閱相關(guān)Matlab/Simulink參考書，熟悉能解決題目問題的相關(guān)Matlab函數(shù)。構(gòu)

2、建傳遞函數(shù)：H = tf(num,den);狀態(tài)空間：A,B,C,D = tf2ss(num,den);部分分式：r,p,k = residue(num,den);聯(lián)接：串聯(lián)H1*H2;并聯(lián)H1+H2;反饋feedback(H1,H2,-1);伯德圖：bode(H);根軌跡：rlocus(H);rlocus(-H);奈奎斯特：Nyquist(H);時(shí)域階躍響應(yīng)：Step(H);Grid on3、 實(shí)驗(yàn)內(nèi)容1已知H（s）=，求H（s）的零極點(diǎn)表達(dá)式和狀態(tài)空間表達(dá)式。答：先創(chuàng)建傳遞函數(shù)H(s)>> num = 0.05 1;>> den = 0.02 0.3 1;>&

3、gt; H = tf(num,den)H = 0.05 s + 1 - 0.02 s2 + 0.3 s + 1 Continuous-time transfer function.H（s）的零極點(diǎn)表達(dá)式如下：>> Z,P,K = tf2zp(num,den)Z = -20P = -10 -5K = 2.5000H（s）的狀態(tài)空間表達(dá)式如下：>> A,B,C,D = tf2ss(num,den)A = -15 -50 1 0B = 1 0C = 2.5000 50.0000D = 02已知，。（1） 求兩模型串聯(lián)后的系統(tǒng)傳遞函數(shù)。答：先創(chuàng)建傳遞函數(shù)H(s)>>

4、 num1 = 1 5;>> den1 = 1 3 2 0;>> H1 = tf(num1,den1)H1 = s + 5 - s3 + 3 s2 + 2 s Continuous-time transfer function.>> num2 = 1;>> den2 = 1 1;>> H2 = tf(num2,den2)H2 = 1 - s + 1 Continuous-time transfer function.兩模型串聯(lián)后的系統(tǒng)傳遞函數(shù)如下：>> H3 = H1*H2H3 = s + 5 - s4 + 4 s3 +

5、5 s2 + 2 s Continuous-time transfer function.（2） 求兩模型并聯(lián)后的系統(tǒng)傳遞函數(shù)。答：兩模型并聯(lián)后的系統(tǒng)傳遞函數(shù)如下：>> H4 = H1+H2H4 = s3 + 4 s2 + 8 s + 5 - s4 + 4 s3 + 5 s2 + 2 s Continuous-time transfer function.（3） 求兩模型在負(fù)反饋連接下的系統(tǒng)傳遞函數(shù)。答：兩模型在負(fù)反饋連接下的系統(tǒng)傳遞函數(shù)如下：>> H5 = feedback(H1,H2,-1)H5 = s2 + 6 s + 5 - s4 + 4 s3 + 5 s2 +

6、 3 s + 5 Continuous-time transfer function.3 作出上題中（1）的BODE圖，并求出幅值裕度與相位裕度。答：作出bode圖命令如下：>> bode(H3)Bode圖如下：標(biāo)定出如下：從上可以看出：幅值裕度為：4.95dB相位裕度為：-19°使用margin函數(shù)直接可求取得（l為幅值裕度，f為相位裕度）：>> l,f = margin(H3)Warning: The closed-loop system is unstable. > In warning at 26 In DynamicSystem.margin

7、at 63 l = 0.5576f = -19.36624 給定系統(tǒng)開環(huán)傳遞函數(shù)為，繪制系統(tǒng)的根軌跡圖與奈奎斯特曲線，并求出系統(tǒng)穩(wěn)定時(shí)的增益K的范圍。答：構(gòu)建傳遞函數(shù)>> num = 1;>> den = conv(1 2,1 2 5);>> G = tf(num,den)G = 1 - s3 + 4 s2 + 9 s + 10 Continuous-time transfer function.繪制系統(tǒng)的根軌跡圖（0至無(wú)窮大）如下：>> rlocus(G)由第二張圖可知K的范圍是：繪制系統(tǒng)的根軌跡圖（負(fù)無(wú)窮大至0）如下：>> rlo

8、cus(-G)由第二張圖可知K的范圍是：將二張圖放在一起，即得：>> rlocus(G);>> hold on>> rlocus(-G)>> hold on綜上，可知K的范圍是：繪制k=1時(shí)奈奎斯特曲線如下：>> num = 1;>> den = conv(1 2,1 2 5);>> G = tf(num,den);由穩(wěn)定條件判別方法可知，故不穩(wěn)定。繪制k=30時(shí)奈奎斯特曲線如下：由穩(wěn)定條件判別方法可知，故不穩(wěn)定。繪制k=-3時(shí)奈奎斯特曲線如下：由穩(wěn)定條件判別方法可知，故穩(wěn)定。繪制k=-11時(shí)奈奎斯特曲線如下：

9、由穩(wěn)定條件判別方法可知，故不穩(wěn)定。5 對(duì)內(nèi)容4中的系統(tǒng)，當(dāng)K=10和40時(shí)，分別作出閉環(huán)系統(tǒng)的階躍響應(yīng)曲線，要求用Simulink實(shí)現(xiàn)。答：Simulink模擬原理圖如下：K = 10 時(shí)的階躍響應(yīng)曲線如下：可見時(shí)最終會(huì)達(dá)到穩(wěn)定。K = 40時(shí)的階躍響應(yīng)曲線如下：可見時(shí)最終振蕩，不會(huì)達(dá)到穩(wěn)定。設(shè)置正反饋的Simulink模擬原理圖如下：K = 5時(shí)的階躍響應(yīng)曲線如下：可見時(shí)最終會(huì)達(dá)到穩(wěn)定。4、 實(shí)驗(yàn)總結(jié)1，這次實(shí)驗(yàn)讓我更進(jìn)一步的理解了bode圖和奈奎斯特圖，尤其深化理解了bode圖切線位置為我們手工畫bode圖的拐點(diǎn)處，即可以近似為特征根的解。2，指標(biāo) - 模型 = 調(diào)節(jié)規(guī)律。3，本次實(shí)驗(yàn)是系統(tǒng)頻率特性的測(cè)試，這章