matlab在控制系統(tǒng)中的應(yīng)用

1第五章Matlab在控制理論中的應(yīng)用

5-1控制系統(tǒng)數(shù)學(xué)模型的Matlab描述在線性系統(tǒng)理論中，一般常用的數(shù)學(xué)模型形式有：傳遞函數(shù)模型（系統(tǒng)的外部模型）;狀態(tài)方程模型（系統(tǒng)的內(nèi)部模型）;零極點(diǎn)增益模型和部分分式模型等。這些模型之間都有著內(nèi)在的聯(lián)系，可以相互進(jìn)行轉(zhuǎn)換。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第1頁!2一、連續(xù)系統(tǒng)的傳遞函數(shù)模型連續(xù)系統(tǒng)的傳遞函數(shù)如下：對(duì)線性定常系統(tǒng)，式中s的系數(shù)均為常數(shù)，且a1不等于零，這時(shí)系統(tǒng)在MATLAB中可以方便地由分子和分母系數(shù)構(gòu)成的兩個(gè)向量唯一地確定出來，這兩個(gè)向量分別用num和den表示。

num=[b1,b2,…,bm,bm+1] den=[a1,a2,…,an,an+1]

注意：它們都是按s的降冪進(jìn)行排列的。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第2頁!3零極點(diǎn)模型實(shí)際上是傳遞函數(shù)模型的另一種表現(xiàn)形式，其原理是分別對(duì)原系統(tǒng)傳遞函數(shù)的分子、分母進(jìn)行分解因式處理，以獲得系統(tǒng)的零點(diǎn)和極點(diǎn)的表示形式。在MATLAB中零極點(diǎn)增益模型用[z,p,K]矢量組表示。即：z=[z1,z2,…,zm]p=[p1,p2,...,pn]K=[k]函數(shù)tf2zp()可以用來求傳遞函數(shù)的零極點(diǎn)和增益。二、零極點(diǎn)增益模型K為系統(tǒng)增益，zi為零點(diǎn)，pj為極點(diǎn)matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第3頁!4控制系統(tǒng)常用到并聯(lián)系統(tǒng)，這時(shí)就要對(duì)系統(tǒng)函數(shù)進(jìn)行分解，使其表現(xiàn)為一些基本控制單元的和的形式。函數(shù)[r,p,k]=residue(b,a)對(duì)兩個(gè)多項(xiàng)式的比進(jìn)行部分展開，以及把傳函分解為微分單元的形式。向量b和a是按s的降冪排列的多項(xiàng)式系數(shù)。部分分式展開后，余數(shù)返回到向量r，極點(diǎn)返回到列向量p，常數(shù)項(xiàng)返回到k。[b,a]=residue(r,p,k)可以將部分分式轉(zhuǎn)化為多項(xiàng)式比p(s)/q(s)。四、部分分式展開matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第4頁!5例4-2零極點(diǎn)增益模型：num=[1,11,30,0];den=[1,9,45,87,50];[z,p,k]=tf2zp(num,den)》z=0-6-5p=-3.0000+4.0000i-3.0000-4.0000i-2.0000-1.0000k=1結(jié)果表達(dá)式：matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第5頁!6例4-4：系統(tǒng)為一個(gè)兩輸入兩輸出系統(tǒng)A=[16910;31268;47911;5121314];B=[46;24;22;10];C=[0021;8022];D=zeros(2,2);matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第6頁!7例4-5：1）已知系統(tǒng)狀態(tài)空間模型為：A=[01;-1-2];B=[0;1];C=[1,3];D=[1];[num,den]=ss2tf(A,B,C,D)num=152;den=121;[z,p,k]=ss2zp(A,B,C,D)z=-4.5616p=-1k=1-0.4384-1matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第7頁!83）系統(tǒng)的零極點(diǎn)增益模型：z=[-3];p=[-1,-2,-5];k=6;[num,den]=zp2tf(z,p,k)num=00618den=181710[a,b,c,d]=zp2ss(z,p,k)a=-1.000000b=12.0000-7.0000-3.1623103.162300c=001.8974d=0注意：零極點(diǎn)的輸入可以寫出行向量，也可以寫出列向量。

matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第8頁!91、并聯(lián)：parallel格式：[a,b,c,d]=parallel(a1,b1,c1,d1,a2,b2,c2,d2)％并聯(lián)連接兩個(gè)狀態(tài)空間系統(tǒng)。[num,den]=parallel(num1,den1,num2,den2)％將并聯(lián)連接的傳遞函數(shù)進(jìn)行相加。5.2.2模型的連接matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第9頁!103、反饋：feedback格式：[a,b,c,d]=feedback(a1,b1,c1,d1,a2,b2,c2,d2,sign)％系統(tǒng)1的所有輸出連接到系統(tǒng)2的輸入，系統(tǒng)2的所有輸出連接到系統(tǒng)1的輸入，sign用來指示系統(tǒng)2輸出到系統(tǒng)1輸入的連接符號(hào)，sign缺省時(shí)，默認(rèn)為負(fù)，即sign=-1?？傁到y(tǒng)的輸入/輸出數(shù)等同于系統(tǒng)1。[num,den]=feedback(num1,den1,num2,den2,sign)％可以得到類似的連接，只是子系統(tǒng)和閉環(huán)系統(tǒng)均以傳遞函數(shù)的形式表示。sign的含義與前述相同。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第10頁!11在進(jìn)行控制系統(tǒng)的仿真之前，建立系統(tǒng)的模型表達(dá)式是關(guān)鍵的一步。對(duì)于控制系統(tǒng)，有不同的分類，在本課程中主要討論的是線性定常連續(xù)系統(tǒng)系統(tǒng)的描述有不同的方法：微分方程；傳遞函數(shù)；零極點(diǎn)增益模式；部分分式展開；狀態(tài)空間模型等。系統(tǒng)的模型之間可以相互轉(zhuǎn)換，要求熟練掌握各種模型之間轉(zhuǎn)換的命令。模型之間可以進(jìn)行連接，常用的模型連接命令：串聯(lián)、并聯(lián)、反饋。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第11頁!12控制系統(tǒng)的分析方法控制系統(tǒng)的穩(wěn)定性分析控制系統(tǒng)的時(shí)域分析控制系統(tǒng)的頻域分析控制系統(tǒng)的根軌跡分析matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第12頁!13一、時(shí)域分析的一般方法控制系統(tǒng)仿真常用的幾個(gè)時(shí)間響應(yīng)函數(shù)如下：5.3.2控制系統(tǒng)的時(shí)域分析

一個(gè)動(dòng)態(tài)系統(tǒng)的性能常用典型輸入作用下的響應(yīng)來描述。響應(yīng)是指零初始值條件下某種典型的輸入函數(shù)作用下對(duì)象的響應(yīng)，控制系統(tǒng)常用的輸入函數(shù)為單位階躍函數(shù)和脈沖激勵(lì)函數(shù)（即沖激函數(shù)）。在MATLAB的控制系統(tǒng)工具箱中提供了求取這兩種輸入下系統(tǒng)響應(yīng)的函數(shù)。求取系統(tǒng)單位階躍響應(yīng)：step()求取系統(tǒng)的沖激響應(yīng)：impulse()matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第13頁!14例1、系統(tǒng)無零點(diǎn)與系統(tǒng)有零點(diǎn)表達(dá)式分別為試比較系統(tǒng)階躍響應(yīng)的差別。num1=[10];num2=[210];den=[1210];step(num1,den)holdonstep(num2,den)

解：matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第14頁!15impulse(num,den)

給定系統(tǒng)多項(xiàng)式模型num，den，求系統(tǒng)的單位脈沖響應(yīng)。繪制時(shí)間響應(yīng)曲線；也可以使用狀態(tài)空間模型impulse(A,B,C,D)。

[y,x,t]=impulse(num,den)返回變量格式，不繪圖。

y=y(t)為輸出響應(yīng)，x=x(t)為狀態(tài)變量，t為時(shí)間相量。線性系統(tǒng)的穩(wěn)態(tài)值可以通過函數(shù)dcgain()來求取，其調(diào)用格式為：dc=dcgain(num,den)或dc=dcgain(a,b,c,d)2、impulse()函數(shù)的用法matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第15頁!【例3】已知二階閉環(huán)控制系統(tǒng)如圖所示，試在4個(gè)子圖中繪出當(dāng)無阻尼自然振蕩頻率，阻尼比ζ

分別為0.2、0.5、1.0和2.5等不同值時(shí)，系統(tǒng)的單位階躍響應(yīng)曲線。解建立一個(gè)m文件，不妨命名為step4re.m如下將該m文件保存在work文件夾中，然后在CommandWindow中鍵入step4re，回車。計(jì)算機(jī)就分別在4個(gè)子圖中繪出4個(gè)單位階躍響應(yīng)曲線。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第16頁!matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第17頁!183.求特定輸入u下的時(shí)間響應(yīng)：lsim(num,den,u,T)

給定系統(tǒng)多項(xiàng)式模型num，den，任意控制輸入向量u，時(shí)間向量T，繪制時(shí)間響應(yīng)曲線；也可以使用狀態(tài)空間模型lsim(A,B,C,D,u,T)。

[y,x,t]=lsim(num,den)

返回變量格式，不繪圖。

y=y(t)為輸出響應(yīng)，x=x(t)為狀態(tài)變量，t為時(shí)間向量。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第18頁!19（a）I型系統(tǒng)（b）II型系統(tǒng)

兩個(gè)系統(tǒng)對(duì)斜坡信號(hào)的跟蹤曲線如上圖所示。I型系統(tǒng)在跟蹤速率信號(hào)時(shí)，有明顯的跟蹤誤差，即ess=1，如圖(a)所示；但是II型系統(tǒng)由圖(b)可見，穩(wěn)態(tài)跟蹤誤差為零，即ess=0。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第19頁!20常用時(shí)域分析函數(shù)

時(shí)間響應(yīng)探究系統(tǒng)對(duì)輸入和擾動(dòng)在時(shí)域內(nèi)的瞬態(tài)行為，系統(tǒng)特征如：上升時(shí)間、調(diào)節(jié)時(shí)間、超調(diào)量和穩(wěn)態(tài)誤差都能從時(shí)間響應(yīng)上反映出來。MATLAB除了提供前面介紹的對(duì)系統(tǒng)階躍響應(yīng)、沖激響應(yīng)等進(jìn)行仿真的函數(shù)外，還提供了大量對(duì)控制系統(tǒng)進(jìn)行時(shí)域分析的函數(shù)，如：covar：連續(xù)系統(tǒng)對(duì)白噪聲的方差響應(yīng)initial：連續(xù)系統(tǒng)的零輸入響應(yīng)對(duì)于離散系統(tǒng)只需在連續(xù)系統(tǒng)對(duì)應(yīng)函數(shù)前加d就可以，如dstep，dimpulse等。它們的調(diào)用格式與step、impulse類似，可以通過help命令來察看自學(xué)。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第20頁!21nichols(num,den)

給定控制系統(tǒng)的開環(huán)多項(xiàng)式模型num，den，繪制尼柯爾斯圖；也可以使用狀態(tài)空間模型nichols(A,B,C,D)。

[m,p,w]=nichols(num,den)返回變量格式，不繪圖。

m為模矩陣，p為相位角矩陣，w為頻率向量。nyquist(num,den)

給定控制系統(tǒng)的開環(huán)多項(xiàng)式模型num，den，繪制奈奎斯特圖；也可以使用狀態(tài)空間模型nyquist(A,B,C,D)。

[re,im,w]=nyquist(num,den)返回變量格式，不繪圖。

re為返回的實(shí)部向量，im為返回的虛部向量，w為頻率向量。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第21頁!22采用輸出裕量的方法：num=[15];den=[12100];margin(num,den)

得到圖線如右圖所示，從圖上可以讀到開環(huán)截止頻率為wc=2.24，相位裕度為Pm=48.2度，相位裕度Pm較大，但是由于二階環(huán)節(jié)諧振峰值的影響，幅值裕度Gm較小，且Gm=2.5dB。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第22頁!23num=[45];den=[16.0234.12104.682.08];nyquist(num,den)例6、控制系統(tǒng)的開環(huán)傳遞函數(shù)為繪制其nyquist圖。

由于坐標(biāo)軸范圍的設(shè)置是自動(dòng)的，為了表現(xiàn)全圖的概貌，原點(diǎn)處的圖線就變的不夠清晰。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第23頁!243、繪制根軌跡圖rlocus(num,den)

給定控制系統(tǒng)的開環(huán)多項(xiàng)式模型num，den，繪制閉環(huán)根的根軌跡圖；也可以使用狀態(tài)空間模型rlocus(A,B,C,D)。根軌跡繪圖的幾個(gè)常用的函數(shù)：[k,r]=rlocfind(num,den)

給定控制系統(tǒng)的開環(huán)多項(xiàng)式模型num，den，作根軌跡圖，并用鼠標(biāo)在圖上選定閉環(huán)極點(diǎn)的位置，返回對(duì)應(yīng)的增益k和閉環(huán)極點(diǎn)的值。pzmap(num,den)

在復(fù)數(shù)平面上畫出零極點(diǎn)的位置。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第24頁!25>>[k,r]=rlocfind(num,den)Selectapointinthegraphicswindowselected_point=0+1.4247ik=4.0595r=-2.00990.0049+1.4212i0.0049-1.4212i>>輸入指令以后，用鼠標(biāo)在圖中選取特征點(diǎn)（本例中為根軌跡與虛軸的交點(diǎn)），即解得臨界增益約為k=4.06。選擇點(diǎn)以后的圖如右圖所示。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第25頁!為了避免純微分運(yùn)算，通常采用近似的PID控制器，其傳遞函數(shù)為2.離散PID控制器離散PID控制器的表達(dá)式為簡化為離散PID控制器的脈沖傳遞函數(shù)為matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第26頁!解使用期望特性法來設(shè)計(jì)PID控制器。假設(shè)PID控制器的傳遞函數(shù)為系統(tǒng)閉環(huán)的傳遞函數(shù)為不妨假設(shè)希望閉環(huán)極點(diǎn)為：，，和則期望特征多項(xiàng)式為對(duì)應(yīng)系數(shù)相等，可求得：在CommandWindow中輸入這3個(gè)參數(shù)值matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第27頁!3.3使用Ziegler-Nichols經(jīng)驗(yàn)整定公式進(jìn)行PID控制器設(shè)計(jì)Ziegler-Nichols經(jīng)驗(yàn)整定公式是針對(duì)被控對(duì)象模型為帶有延遲的一階慣性傳遞函數(shù)提出的matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第28頁!MATLAB/Simulink模型如圖所示，仿真結(jié)果如圖所示，可見系統(tǒng)可以穩(wěn)定工作。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第29頁!二、實(shí)驗(yàn)內(nèi)容1采用求取特征值的方法判別下面系統(tǒng)的穩(wěn)定性。2某單位負(fù)反饋系統(tǒng)的開環(huán)控制系統(tǒng)的傳遞函數(shù)為（1）繪制系統(tǒng)的根軌跡；（2）當(dāng)時(shí)，繪制系統(tǒng)的Bode圖，判斷系統(tǒng)的穩(wěn)定性，并且求出幅值裕度和相角裕度。。某的開環(huán)控制系統(tǒng)的傳遞函數(shù)為

試?yán)L制奈奎斯特曲線，并利用奈氏判據(jù)判別閉環(huán)系統(tǒng)的穩(wěn)定性。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第30頁!31狀態(tài)方程與輸出方程的組合稱為狀態(tài)空間表達(dá)式，又稱為動(dòng)態(tài)方程，經(jīng)典控制理論用傳遞函數(shù)將輸入—輸出關(guān)系表達(dá)出來，而現(xiàn)代控制理論則用狀態(tài)方程和輸出方程來表達(dá)輸入—輸出關(guān)系，揭示了系統(tǒng)內(nèi)部狀態(tài)對(duì)系統(tǒng)性能的影響。三狀態(tài)方程描述在MATLAB中，系統(tǒng)狀態(tài)空間用（A,B,C,D)矩陣組表示。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第31頁!32例4-1：傳遞函數(shù)描述1）》num=[12,24,0,20];den=[24622];2）借助多項(xiàng)式乘法函數(shù)conv來處理：》num=4*conv([1,2],conv([1,6,6],[1,6,6]));》den=conv([1,0],conv([1,1],conv([1,1],conv([1,1],[1,3,2,5]))));matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第32頁!33例4-3部分分式展開：num=[2,0,9,1];den=[1,1,4,4];[r,p,k]=residue(num,den)p=0.0000+2.0000i0.0000-2.0000i-1.0000k=2r=0.0000-0.2500i0.0000+0.2500i-2.0000結(jié)果表達(dá)式：matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第33頁!34在一些場合下需要用到某種模型，而在另外一些場合下可能需要另外的模型，這就需要進(jìn)行模型的轉(zhuǎn)換。模型轉(zhuǎn)換的函數(shù)包括：residue：傳遞函數(shù)模型與部分分式模型互換ss2tf：狀態(tài)空間模型轉(zhuǎn)換為傳遞函數(shù)模型ss2zp：狀態(tài)空間模型轉(zhuǎn)換為零極點(diǎn)增益模型tf2ss：傳遞函數(shù)模型轉(zhuǎn)換為狀態(tài)空間模型tf2zp：傳遞函數(shù)模型轉(zhuǎn)換為零極點(diǎn)增益模型zp2ss：零極點(diǎn)增益模型轉(zhuǎn)換為狀態(tài)空間模型zp2tf：零極點(diǎn)增益模型轉(zhuǎn)換為傳遞函數(shù)模型

5-2控制模型的轉(zhuǎn)換與連接5.2.1、模型的轉(zhuǎn)換matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第34頁!352）已知一個(gè)單輸入三輸出系統(tǒng)的傳遞函數(shù)模型為：num=[00-2;0-1-5;120];den=[16116];[A,B,C,D]=tf2ss(num,den)A=-6-11-6B=1C=00-2D=010000-1-5001001200matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第35頁!364）已知部分分式：r=[-0.25i,0.25i,-2];p=[2i,-2i,-1];k=2;[num,den]=residue(r,p,k)num=2091den=1144注意余式一定要與極點(diǎn)相對(duì)應(yīng)。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第36頁!372、串聯(lián)：series格式：[a,b,c,d]=series(a1,b1,c1,d1,a2,b2,c2,d2)％串聯(lián)連接兩個(gè)狀態(tài)空間系統(tǒng)。[num,den]=series(num1,den1,num2,den2)％將串聯(lián)連接的傳遞函數(shù)進(jìn)行相乘。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第37頁!38ctrb和obsv函數(shù)可以求出狀態(tài)空間系統(tǒng)的可控性和可觀性矩陣。 格式：co=ctrb(a,b)ob=obsv(a,c)對(duì)于n×n矩陣a，n×m矩陣b和p×n矩陣cctrb(a,b)可以得到n×nm的可控性矩陣co=[baba2b…an-1b]obsv(a,c)可以得到nm×n的可觀性矩陣ob=[ccaca2…can-1]’當(dāng)co的秩為n時(shí)，系統(tǒng)可控；當(dāng)ob的秩為n時(shí)，系統(tǒng)可觀。5.2.3模型的屬性matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第38頁!395-3控制系統(tǒng)的分析方法早期的控制系統(tǒng)分析過程復(fù)雜而耗時(shí)，如想得到一個(gè)系統(tǒng)的沖激響應(yīng)曲線，首先需要編寫一個(gè)求解微分方程的子程序，然后將已經(jīng)獲得的系統(tǒng)模型輸入計(jì)算機(jī)，通過計(jì)算機(jī)的運(yùn)算獲得沖激響應(yīng)的響應(yīng)數(shù)據(jù)，然后再編寫一個(gè)繪圖程序，將數(shù)據(jù)繪制成可供工程分析的響應(yīng)曲線。MATLAB控制系統(tǒng)工具箱和SIMULINK輔助環(huán)境的出現(xiàn)，給控制系統(tǒng)分析帶來了福音?？刂葡到y(tǒng)的分析包括系統(tǒng)的穩(wěn)定性分析、時(shí)域分析、頻域分析及根軌跡分析。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第39頁!405.3.1控制系統(tǒng)的穩(wěn)定性分析對(duì)于連續(xù)時(shí)間系統(tǒng)，如果閉環(huán)極點(diǎn)全部在S平面左半平面，則系統(tǒng)是穩(wěn)定的。對(duì)于離散時(shí)間系統(tǒng)，如果系統(tǒng)全部極點(diǎn)都位于Z平面的單位圓內(nèi)，則系統(tǒng)是穩(wěn)定的。MATLAB提供了直接求取系統(tǒng)所有零極點(diǎn)的函數(shù)，因此可以直接根據(jù)零極點(diǎn)的分布情況對(duì)系統(tǒng)的穩(wěn)定性matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第40頁!411、step()函數(shù)的用法y=step(num,den,t)：其中num和den分別為系統(tǒng)傳遞函數(shù)描述中的分子和分母多項(xiàng)式系數(shù)，t為選定的仿真時(shí)間向量，一般可以由t=0:step:end等步長地產(chǎn)生出來。該函數(shù)返回值y為系統(tǒng)在仿真時(shí)刻各個(gè)輸出所組成的矩陣。如果對(duì)具體的響應(yīng)值不感興趣，而只想繪制系統(tǒng)的階躍響應(yīng)曲線，可調(diào)用以下的格式：step(num,den)；step(num,den,t)；[y,x,t]=step(num,den)：此時(shí)時(shí)間向量t由系統(tǒng)模型的特性自動(dòng)生成。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第41頁!42得到的圖線：計(jì)算超調(diào)量的命令：[y1,t1]=step(num1,den);[y2,t2]=step(num2,den);[max(y1),max(y2)]ans=1.35091.4348>>可見，超調(diào)量分別是35.09%和43.48%。

顯然，系統(tǒng)有s=-5的零點(diǎn)時(shí)，其階躍響應(yīng)疊加了微分項(xiàng)，超調(diào)量稍有增加。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第42頁!43例2、控制系統(tǒng)的傳遞函數(shù)分別為作系統(tǒng)的脈沖響應(yīng)曲線。num1=[001];num2=[010];num3=[100];den=[1210];impulse(num1,den);holdon;impulse(num2,den);impulse(num3,den);也可以簡單的將參數(shù)填到命令中：>>impulse([001],[1210]);holdon;>>impulse([010],[1210]);impulse([100],[1210]);解：matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第43頁!matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第44頁!45得到的圖線結(jié)果：（在matlab彈出的圖形窗口中，選擇Edit/CopyFigure，可以實(shí)現(xiàn)拷貝圖像功能）matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第45頁!46例4、系統(tǒng)結(jié)構(gòu)如圖所示，考察I型系統(tǒng)與II型系統(tǒng)對(duì)速率信號(hào)的跟蹤能力。+-+-解：繪制圖a：num1=[2];den1=[122];t=0:0.1:10;u=t;lsim(num1,den1,u,t);holdonplot(t,u,‘r:’)%將t-u直線描紅色繪制圖b：num2=[22];den2=[1222];t=0:0.1:10;u=t;lsim(num2,den2,u,t);holdonplot(t,u,'r:')%將t-u直線描紅色圖aI型系統(tǒng)圖bII型系統(tǒng)matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第46頁!47>>num1=[2];>>den1=[122];>>t=0:0.1:10;u=t;>>[y]=lsim(num1,den1,u,t);>>y=y‘;>>plot(t,u-y);>>num2=[22];>>den2=[1222];>>t=0:0.1:50;u=t;>>[y]=lsim(num2,den2,u,t);>>y=y‘;>>plot(t,u-y);誤差響應(yīng)曲線圖：matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第47頁!485.3.3頻率特性繪圖控制系統(tǒng)頻率特性分析常用的幾個(gè)繪圖函數(shù)為：bode(num,den)

給定控制系統(tǒng)的開環(huán)多項(xiàng)式模型num，den，繪制伯德圖；也可以使用狀態(tài)空間模型bode(A,B,C,D)。

[m,p,w]=bode(num,den)返回變量格式，不繪圖。

m=m(w)為幅值向量，p=p(w)為相位向量，w為頻率向量。margin(num,den)

給定控制系統(tǒng)的開環(huán)多項(xiàng)式模型num，den，繪制帶有穩(wěn)定裕度的伯德圖；也可以使用狀態(tài)空間模型margin(A,B,C,D)。

[gm,pm,wg,wp]=margin(num,den)返回變量格式，不繪圖。

gm為幅值裕度，pm為相位裕度，wg是相角等于-π時(shí)的頻率，wp是模等于1時(shí)的頻率。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第48頁!49例4、控制系統(tǒng)的開環(huán)傳遞函數(shù)為作伯德圖并確定系統(tǒng)的穩(wěn)定裕度。解：如果單純做bode圖：num=[15];den=[12100];bode(num,den)

得到圖線如右圖所示，圖中沒有各種裕度的顯示，不能滿足題目的要求。matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第49頁!50例5、控制系統(tǒng)的開環(huán)傳遞函數(shù)為繪制其nyquist圖。解：num=[26];den=[1-10];nyquist(num,den)matlab在控制系統(tǒng)中的應(yīng)用共58頁，您現(xiàn)在瀏覽的是第50頁!51下面我們加上一句語句，將坐標(biāo)軸