控制工程基礎(chǔ)(經(jīng)典控制部分)MATLAB分析.doc

控制工程基礎(chǔ)（經(jīng)典控制部分）的MATLAB分析第一章MATLAB的基本使用1-1MATLAB語言簡(jiǎn)介MATLAB是一種高級(jí)矩陣語言，它由Math Works公司于1984年正式推出，它的基本處理對(duì)象是矩陣，即使是一個(gè)標(biāo)量純數(shù)，MATLAB也認(rèn)為它是只有一個(gè)元素的矩陣。隨著MATLAB的發(fā)展，特別是它所包含的大量工具箱（應(yīng)用程序集）的集結(jié)，使MATLAB已經(jīng)成為帶有獨(dú)特?cái)?shù)據(jù)結(jié)構(gòu)、輸入輸出、流程控制語句和函數(shù)、并且面向?qū)ο蟮母呒?jí)語言。MATLAB語言被稱為一種“演算紙式的科學(xué)計(jì)算語言”，它在數(shù)值計(jì)算、符號(hào)運(yùn)算、數(shù)據(jù)處理、自動(dòng)控制、信號(hào)處理、神經(jīng)網(wǎng)絡(luò)、優(yōu)化計(jì)算、模糊邏輯、系統(tǒng)辨識(shí)、小波分析、圖象處理、統(tǒng)計(jì)分析、甚至于金融財(cái)會(huì)等廣大領(lǐng)域有著十分廣泛的用途。MATLAB語言在工程計(jì)算與分析方面具有無可比擬的優(yōu)異性能。它集計(jì)算、數(shù)據(jù)可視化和程序設(shè)計(jì)于一體，并能將問題和解決方案以使用者所熟悉的數(shù)學(xué)符號(hào)或圖形表示出來。MATLAB語言和C語言的關(guān)系與C語言和匯編語言的關(guān)系類似。例如當(dāng)我們需要求一個(gè)矩陣的特征值時(shí),在MATLAB下只需由幾個(gè)字符組成的一條指令即可得出結(jié)果，而不必去考慮用什么算法以及如何實(shí)現(xiàn)這些算法等低級(jí)問題，也不必深入了解相應(yīng)算法的具體內(nèi)容。就象在C語言下不必象匯編語言中去探究乘法是怎樣實(shí)現(xiàn)的，而只需要采用乘積的結(jié)果就可以了。MATLAB語言還有一個(gè)巨大的優(yōu)點(diǎn)是其高度的可靠性。例如對(duì)于一個(gè)病態(tài)矩陣的處理，MATLAB不會(huì)得出錯(cuò)誤的結(jié)果，而用C或其它高級(jí)語言編寫出來的程序可能會(huì)得出錯(cuò)誤的結(jié)果。這是因?yàn)镸ATLAB函數(shù)集及其工具箱都是由一些在該領(lǐng)域卓有研究成果，造詣很深的權(quán)威學(xué)者經(jīng)過反復(fù)比較所得出來的最優(yōu)方法，而且經(jīng)過多年的實(shí)踐檢驗(yàn)被證明是正確可靠的。1-2MATLAB的工作窗口下面以MATLAB6.1為例介紹。從實(shí)用的角度MATLAB的工作窗口包括命令窗口、M文件編輯器窗口、圖形編輯窗口、數(shù)學(xué)函數(shù)庫、應(yīng)用程序接口及在線窗口。下面首先介紹MATLAB的命令窗口及M文件編輯器。一、命令窗口啟動(dòng)MATLAB之后，屏幕上自動(dòng)出現(xiàn)命令窗口MATLAB，它是MATLAB提供給用戶的操作界面，用戶可以在命令窗口內(nèi)提示符“”之后（有的MATLAB版本命令窗口沒有提示符）鍵入MATLAB命令，回車即獲得該命令的答案。命令窗口內(nèi)有File、Edit、View、Web、Window、Help等菜單條。二、M文件編輯窗口M文件是MATLAB語言所特有的文件。用戶可以在M文件編輯窗口內(nèi)，編寫一段程序，調(diào)試，運(yùn)行并存盤，所保存的用戶程序即是用戶自己的M文件。MATLAB工具箱中大量的應(yīng)用程序也是以M文件的形式出現(xiàn)的，這些M文件可以打開來閱讀，甚至修改，但應(yīng)注意，不可改動(dòng)工具箱中的M文件！1進(jìn)入M文件窗口有兩種方法1) 命令窗口 File New M-File；2) 命令窗口 點(diǎn)擊“File”字樣下面的圖標(biāo)。M文件編輯窗口的標(biāo)記是“Untitled”（無標(biāo)題的）。當(dāng)用戶編寫的程序要存盤時(shí)，Untitled作為默認(rèn)文件名提供給用戶，自然，用戶可以，也應(yīng)當(dāng)自己命名。若用戶不自己命名，則MATLAB會(huì)對(duì)Untitled進(jìn)行編號(hào)。2M文件的執(zhí)行：返回命令窗口，在當(dāng)前目錄（Current Directory）內(nèi)選擇所要運(yùn)行的M文件的目錄，在命令窗口提示符“”后，直接鍵入文件名（不加后綴）即可運(yùn)行。注意：（1）機(jī)器默認(rèn)路徑為一級(jí)子目錄MATLAB6 p1work；（2）MATLAB 6.1以前的版本，運(yùn)行M文件的方法稍有不同，它必須在File菜單下，打開“Run Script”子菜單，鍵入需要運(yùn)行的文件路徑及名稱再回車，在這種情況下，work作為根目錄對(duì)待，不出現(xiàn)在M文件的路徑之中。本講義的參考程序都是在M文件窗口下編制的。三、在線幫助窗口在命令窗口中鍵入Help(空格) 函數(shù)名，可以立即獲得該函數(shù)的使用方法。1-3 MATLAB最基本的矩陣操作作為命令窗口及M文件編輯器的應(yīng)用實(shí)例，介紹幾個(gè)最基本的矩陣運(yùn)算命令。一、矩陣的輸入在方括號(hào)內(nèi)依次按行鍵入矩陣元素，在一行內(nèi)的各元素之間用空格或逗號(hào)分開，每行之間用分號(hào)分開。例如，在命令窗內(nèi)輸入A=2 2 3;4 5 4;7 8 9 (注意：方括號(hào)，分號(hào)為矩陣行標(biāo)記)B=1,3,5;6,-4,2;3,5,1 (逗號(hào)與空格功能相同)A223B135 4546-42 789 351同理：輸入A1246得到行矢量，輸入A22；4；6得到列矢量，于是，當(dāng)輸入C=A；A1有C=123456789246A1作為矩陣C的最后一行，C和A相比，增加了一行。二、矩陣的轉(zhuǎn)置矩陣A的轉(zhuǎn)置用A表示，顯然，A1與A2互為轉(zhuǎn)置，即A1會(huì)得到以2,4，6為元素的列矢量。思考一下輸入C1 AA2 C2 AA1有什么結(jié)果？而輸入A；A1有無意義？三、矩陣的四則運(yùn)算1矩陣的加減法：當(dāng)兩個(gè)矩陣維數(shù)相同時(shí)可以直接進(jìn)行“”或“”運(yùn)算。如 D1AB，D2AB2 矩陣的乘法：當(dāng)矩陣A，B維數(shù)相容時(shí)C3AB：普通意義下的矩陣相乘C4A .B：矩陣A與B的對(duì)應(yīng)元素相乘顯然，ABBA（一般情況），而A .BB .A。A .B稱為數(shù)列型乘法，它要求參加運(yùn)算的矩陣或數(shù)列具有相同的行列數(shù)，這是MATLAB語言中的一種特殊運(yùn)算，它在今后求取函數(shù)值等運(yùn)算時(shí)是很重要的。實(shí)際上，前面所述的矩陣加、減法就是一種數(shù)列型運(yùn)算。3 矩陣的除法D4AB：表示A-1B或inv(A)*B，即A的逆矩陣左乘矩陣B。D5B/A：表示BA-1或Binv(A)，即A的逆矩陣右乘B。D6A .B：表示B的每一個(gè)元素被A的對(duì)應(yīng)元素除。D7A ./B：表示A的每一個(gè)元素被B的對(duì)應(yīng)元素除。顯然，A .B與A ./B的各對(duì)應(yīng)元素互為倒數(shù)。讀者可以思考一下，D6 .D7等于什么？D8inv(A)：A的逆矩陣。打開M文件編輯窗口，將上述命令依次鍵入，得到fanli001如下：參考程序fanli001:矩陣的四則運(yùn)算A=2 2 3;4 5 4;7 8 9 % 三階矩陣輸入B=1,3,5;6,-4,2;3,5,1 % 三階矩陣輸入A1=2 4 6 %行向量A2=2;4;6 %列向量C=A;A1 %矩陣A增加一行C1=A A2%矩陣A增加一列C2=A A1 %矩陣A增加一列D1=A+B %矩陣相加D2=A-B%矩陣相減C3=A*B%矩陣與矩陣相乘C4=A.*B%矩陣的對(duì)應(yīng)元素相乘D3=AB%A的逆左乘BD4=B/A%A的逆右乘BD6=A.B%B的各元素被A的對(duì)應(yīng)元素除D7=A./B%A的各元素被B的對(duì)應(yīng)元素除D8=inv(A) %A的逆矩陣語句后面的%為語句說明符。MATLAB中矩陣運(yùn)算的其它主要命令可通過在線幫助獲得。1-4MATLAB的符號(hào)運(yùn)算操作一、進(jìn)入符號(hào)運(yùn)算功能在命令窗口鍵入syms x y z t 此后，即可以使用x，y，z，t等作自變量定義函數(shù)。syms x y z t real 規(guī)定所定義的變量為實(shí)型。二、代數(shù)方程求解使用命令solve 可以求解代數(shù)方程，如求下例方程組的解，命令為x,y=solve(x2+x*y+y-3=0,x2-4*x-2*y+3=0)程序見范例程序fanli002。參考程序fanli002:代數(shù)方程求解syms x y%進(jìn)入符號(hào)運(yùn)算功能f1=x2+x*y+y-3%函數(shù)f1f2=x2-4*x-2*y+3 %函數(shù)f2x,y=solve(x2+x*y+y-3=0,x2-4*x-2*y+3=0)%求解方程組f1a=simplify(subs(f1) %用求解出的x,y檢驗(yàn)方程1f2a=simplify(subs(f2) %用求解出的x,y檢驗(yàn)方程2x=double(x) %將符號(hào)變量轉(zhuǎn)換成浮點(diǎn)數(shù)y=double(y) %將符號(hào)變量轉(zhuǎn)換成浮點(diǎn)數(shù)f1=subs(f1) %用浮點(diǎn)數(shù)x,y檢驗(yàn)方程1f2=subs(f2) %用浮點(diǎn)數(shù)x,y檢驗(yàn)方程2由solve求出的根是根的符號(hào)表達(dá)形式，是準(zhǔn)確解。命令simplify(f)表示化簡(jiǎn)，subs(f1)表示將求解出x，y代回f1中；double是將符號(hào)變量轉(zhuǎn)換成浮點(diǎn)數(shù)，是準(zhǔn)確解x，y的近似值。這從程序運(yùn)行f1a=0，f2a=0，f10，f20可以確認(rèn)。三、符號(hào)矩陣運(yùn)算符號(hào)矩陣可以和數(shù)值矩陣一樣進(jìn)行運(yùn)算，例如：求矩陣特征值eig，求矩陣的逆inv等命令都支持符號(hào)運(yùn)算。設(shè)計(jì)算其特征值eig(A)，逆矩陣inv(A)，程序見fanli003。參考程序fanli003:符號(hào)矩陣的特征值syms t real %定義為實(shí)型變量A=sin(t) -cos(t);cos(t) sin(t) %定義矩陣AB1=eig(A) %求矩陣A的特征值B1=simple(B1) %化簡(jiǎn)A的特征值表達(dá)式B2=inv(A) %求矩陣A的逆矩陣B2=simple(B2) %化簡(jiǎn)A的逆矩陣表達(dá)式C1=A*B2 %檢驗(yàn)A的逆矩陣C1=simple(C1) %C1為單位矩陣注意函數(shù)的輸入方法，自變量用圓括號(hào)括起來。四、微積分運(yùn)算設(shè)函數(shù)，則MATLAB中微積分運(yùn)算命令為fiff(f)：求函數(shù)f對(duì)自變量x的一階導(dǎo)數(shù)；diff(f，2)：求函數(shù)f對(duì)自變量x的二階導(dǎo)數(shù)；int(f)：求函數(shù)f的不定積分例1.1，設(shè)試計(jì)算其一階，二階導(dǎo)數(shù)，積分運(yùn)算，并作出函數(shù)圖象，見范例fanli004。參考程序fanli004:函數(shù)的微分與積分syms xf=1/(5+4*cos(x)ezplot(f) %函數(shù)f的曲線f1=diff(f) %函數(shù)f的一階導(dǎo)數(shù)figure,ezplot(f1) %函數(shù)f一階導(dǎo)數(shù)的曲線f2=diff(f,2) %函數(shù)f的二階導(dǎo)數(shù)figure,ezplot(f2) %函數(shù)f二階導(dǎo)數(shù)的曲線g=int(int(f2) %函數(shù)f的二階導(dǎo)數(shù)f2的二重積分figure,ezplot(g) %函數(shù)f2二重積分的曲線e=f-g %二階導(dǎo)數(shù)的二重積分與原函數(shù)的差e=simple(e)figure,ezplot(e)程序中ezplot(f)：作函數(shù)的圖形，x的取值范圍默認(rèn)值為-2 x 2fanli004的函數(shù)曲線fanli004的一階導(dǎo)函數(shù)曲線fanli004的二階導(dǎo)函數(shù)曲線fanli004的二階導(dǎo)函數(shù)的二重積分曲線ezplot(f)是一個(gè)很有用的作圖命令，它的其它應(yīng)用形式，請(qǐng)查在線幫助。細(xì)心的讀者會(huì)發(fā)現(xiàn)，一個(gè)函數(shù)求二階導(dǎo)數(shù)后再對(duì)二階導(dǎo)數(shù)進(jìn)行二重積分，其結(jié)果與原函數(shù)相差一個(gè)常數(shù)。相當(dāng)于縱坐標(biāo)發(fā)生平移。第二章系統(tǒng)的時(shí)域特性2-1傳遞函數(shù)一、傳遞函數(shù)的兩種形式傳遞函數(shù)通常表達(dá)成s的有理分式形式及零極點(diǎn)增益形式。設(shè)傳遞函數(shù)1有理分式形式分別將分子、分母中s多項(xiàng)式的系數(shù)按降冪排列成行矢量，缺項(xiàng)的系數(shù)用0補(bǔ)齊。上述函數(shù)可表示為num1=2 1 %（注意：方括號(hào)，同一行的各元素間留空格或逗號(hào)）。den1=1 2 2 1syss1=tf(num1,den1)運(yùn)行后，返回傳遞函數(shù)的形式。這種形式不能直接進(jìn)行符號(hào)運(yùn)算！2零極點(diǎn)增益形式Z,P,K = tf2zp(num1,den1)sys2 = zpk(Z,P,K)返回零、極點(diǎn)、增益表達(dá)式，其Z，P分別將零點(diǎn)和極點(diǎn)表示成列向量，若無零點(diǎn)或極點(diǎn)用 （空矩陣）代替。運(yùn)行得到的點(diǎn) Z = -0.5極點(diǎn)P= -1，-0.5j0.866增益K = 2指令zp2tf(Z，P，K)將零極點(diǎn)增益變換成有理分式形式，見程序fanli005。參考程序fanli005:傳遞函數(shù)的有理分式及零極點(diǎn)增益模型num1=2 1 % 傳遞函數(shù)的分子系數(shù)向量den1=1 2 2 1 % 傳遞函數(shù)的分母系數(shù)向量sys1=tf(num1,den1) % 傳遞函數(shù)的有理分式模型Z,P,K=tf2zp(num1,den1) % 有理分式模型轉(zhuǎn)換成零極點(diǎn)增益模型num2,den2=zp2tf(Z,P,K) % 零極點(diǎn)增益模型轉(zhuǎn)換成有理分式模型sys2=zpk(Z,P,K) % 傳遞函數(shù)的零極點(diǎn)增益模型A1,B1,C1,D1=tf2ss(num1,den1) % 有理分式模型轉(zhuǎn)換成狀態(tài)空間模型A2,B2,C2,D2=zp2ss(Z,P,K)% 零極點(diǎn)及增益模型轉(zhuǎn)換成狀態(tài)空間模型num1,den1=ss2tf(A1,B1,C1,D1)% 狀態(tài)空間模型轉(zhuǎn)換成有理分式模型Z,P,K=ss2zp(A2,B2,C2,D2) % 狀態(tài)空間模型轉(zhuǎn)換成零極點(diǎn)增益模型程序中，命令tf2ss，zp2ss及ss2tf，ss2zp是狀態(tài)空間模型與有理分式及零、極點(diǎn)、增益模型之間的相互轉(zhuǎn)換。二、傳遞函數(shù)框圖的處理用框圖可以方便地表示傳遞函數(shù)的并聯(lián)，串聯(lián)及反饋。為簡(jiǎn)潔，僅以有理分式模型為例。G1G1G1+G21 并聯(lián)sysp = parallel(sys1,sys2)num,den=parallel(num1,den1,num2,den2)2 串聯(lián)G1(s)G2(s)G1(s)G2(s)syss = series(sys1,sys2)nums,dens = series(num1,den1,num2,den2)3 反饋G1(s)G2(s)G1(s)G2(s)1+G1(s)G2(s) G3(s)G3(s)sysc=feedback(syss,sys3,1) %默認(rèn)值（-1）numc, denc = feedback(nums, dens, num3, den3)4 單位反饋G1(s)G2(s)G1(s)G2(s)1+G1(s)G2(s) sysd = feedback(syss, 1)numd, dend = feedback(nums, dens, 1, 1) %(單位反饋)上面給出了同一指令的兩種形式，相當(dāng)于兩套平行指令。對(duì)于零極點(diǎn)增益形式，書寫稍復(fù)雜一些，可先用zpk轉(zhuǎn)換成系統(tǒng)形式，或用zp2tf轉(zhuǎn)折換成有理分式形式后再進(jìn)行框圖化簡(jiǎn)操作。三、簡(jiǎn)單函數(shù)的拉普拉斯變換在MATLAB的符號(hào)功能中，可以對(duì)簡(jiǎn)單函數(shù)進(jìn)行拉普拉斯正、逆變換。拉氏正變換：laplace(f(t)拉氏逆變換：ilaplace(L(s)其中為原函數(shù)，為象函數(shù)。命令格式參見fanli007。參考程序fanli007:拉普拉斯變換syms s t w a b cf1=sqrt(b-a)2+w2)/w*exp(-a*t)*sin(w*t+atan(w/(b-a)%原函數(shù)f1L1=laplace(f1) %f1的拉氏變換（象函數(shù)）L1=simple(L1) %化簡(jiǎn)f2=ilaplace(L1)%L1的拉氏逆變換f2=simple(f2) %化簡(jiǎn)在MATLAB中使用laplace及ilaplace命令時(shí)，要注意象、原函數(shù)的符號(hào)，特別是對(duì)初相不等于零的振蕩系統(tǒng)，運(yùn)行結(jié)果常常同手冊(cè)上的結(jié)果相差一個(gè)符號(hào)，這要注意函數(shù)表達(dá)式成立的條件。保險(xiǎn)的辦法是再使用拉氏變換的初值定理確定象、原函數(shù)的符號(hào)。2-2系統(tǒng)時(shí)域特性曲線在MATLAB中，當(dāng)傳遞函數(shù)已知時(shí)，可以方便地求出系統(tǒng)的單位脈沖響應(yīng)、單位階躍響應(yīng)等曲線。一、系統(tǒng)的單位階躍響應(yīng)stepstep有以下幾種格式step(sys)：直接作出sys的單位階躍響應(yīng)曲線。其中sys = tf(num, den) 或sys = zpk(z, p, k)，MATLAB自動(dòng)決定響應(yīng)時(shí)間。step(sys, t)設(shè)定響應(yīng)時(shí)間的單位階躍響應(yīng)。t可以設(shè)定為最大響應(yīng)時(shí)間 t = t終值(秒)，也可以設(shè)置為一個(gè)向量t = 0 : t : t終值注意冒號(hào)的使用。它產(chǎn)生一個(gè)從0到t終值的行矢量，元素之間的間隔為t。step(sys1, sys2, , sysn)在同一幅圖上畫出幾個(gè)系統(tǒng)的單位階躍響應(yīng)。y, t = step(sys)；命令輸出對(duì)應(yīng)時(shí)刻t的各個(gè)單位階躍響應(yīng)值，不畫圖。語句后的分號(hào)控制數(shù)據(jù)的屏幕顯示。如果要查看機(jī)器計(jì)算了多少個(gè)數(shù)據(jù)，可以使用命令size(y)得出的結(jié)果也表明數(shù)據(jù)作為列矢量的行數(shù)。要將計(jì)算出的y, t作成曲線，使用一般的作圖命令plot(t, y)plot后面跟的兩個(gè)參數(shù)橫坐標(biāo)在前，縱坐標(biāo)在后。參考程序見fanli008：參考程序fanli008:系統(tǒng)的單位階躍響應(yīng)num1=4 2den1=2 8 14 11 4sys1=tf(num1,den1) %系統(tǒng)G1(s)num2=2 1den2=1 4 6 7 3sys2=tf(num2,den2) %系統(tǒng)G2(s)y1,t1=step(sys1); %系統(tǒng)G1(s)的單位階躍響應(yīng)數(shù)據(jù)y2,t2=step(sys2); %系統(tǒng)G2(s)的單位階躍響應(yīng)數(shù)據(jù)step(sys1,sys2) %系統(tǒng)G1(s)、G2(s)的單位階躍響曲線figure,step(sys1,sys2,20) %系統(tǒng)G1(s)、G2(s)在自選時(shí)間（20秒）內(nèi)的單位階躍響曲線figure,plot(t1,y1) %系統(tǒng)G1(s)的單位階躍響應(yīng)曲線figure,plot(t2,y2) %系統(tǒng)G2(s)的單位階躍響應(yīng)曲線fanli008:step(sys1,sys2,t)單位階躍曲線二、系統(tǒng)的單位脈沖響應(yīng)impulse命令格式與單位階躍響應(yīng)step的命令格式完全相同，只需將語句中的step用impulse代替即可。針對(duì)同樣的系統(tǒng)，其單位脈沖響應(yīng)的參考程序見fanli009。參考程序fanli009:系統(tǒng)的單位脈沖響應(yīng)num1=4 2den1=2 8 14 11 4sys1=tf(num1,den1)num2=2 1den2=1 4 6 7 3sys2=tf(num2,den2)y1,t1=impulse(sys1); %系統(tǒng)G1(s)的單位脈沖響應(yīng)數(shù)據(jù)y2,t2=impulse(sys2); %系統(tǒng)G2(s)的單位脈沖響應(yīng)數(shù)據(jù)impulse(sys1,sys2) %系統(tǒng)G1(s)、G2(s)的單位脈沖響應(yīng)曲線figure,impulse(sys1,sys2,20)%系統(tǒng)G1(s)、G2(s)在自選時(shí)間（20秒）內(nèi)的單位脈沖響應(yīng)曲線figure,plot(t1,y1) %系統(tǒng)G1(s)的單位脈沖響應(yīng)曲線figure,plot(t2,y2) %系統(tǒng)G2(s)的單位脈沖響應(yīng)曲線hold on,step(sys2) %系統(tǒng)G2(s)的單位階躍和單位脈沖響應(yīng)曲線fanli009:impulse(sys1,sys2)單位脈沖響應(yīng)曲線程序的最后一句hold on是當(dāng)前圖形保護(hù)模式。當(dāng)要將新圖形作在當(dāng)前圖形上時(shí)，必須使用hold on。而figure的含意是另開一個(gè)新的圖形窗口，如果不用figure或hold on，則新的圖形會(huì)占用原圖形窗口，始終只保留一個(gè)最新的圖形窗口。三、一階系統(tǒng)及二階系統(tǒng)的時(shí)域特性一階系統(tǒng)及二階系統(tǒng)是最基本也是最重要的系統(tǒng)，高階系統(tǒng)總可以視為由若干個(gè)一階和（或）二階系統(tǒng)組合構(gòu)成。1 一階系統(tǒng)（設(shè)增益為1）影響系統(tǒng)特性的參數(shù)是其時(shí)間常數(shù)T，T越大，系統(tǒng)慣性越大，響應(yīng)越慢。參考程序fanli010給出了T0.4, 1.2, 2.0, 2.8, 3.6, 4,4六條單位階躍響應(yīng)曲線。參考程序fanli010:一階系統(tǒng)的單位階躍響應(yīng)曲線num=1;i=1;for del=0.1:0.2:1.1 %一階系統(tǒng)時(shí)間常數(shù)遞增間隔den=4*del 1; %一階系統(tǒng)分母向量step(tf(num,den) %一階系統(tǒng)單位階躍響應(yīng)曲線hold on, %不同時(shí)間常數(shù)的一階系統(tǒng)單位階躍響應(yīng)曲線簇i=i+1;end同理，可以作出對(duì)應(yīng)的單位脈沖響應(yīng)曲線，參考程序fanli011。參考程序fanli011:一階系統(tǒng)的單位脈沖響應(yīng)曲線num=1;i=1;fanli010一階系統(tǒng)時(shí)間常數(shù)對(duì)單位階躍響應(yīng)的影響fanli011一階系統(tǒng)時(shí)間常數(shù)對(duì)單位脈沖響應(yīng)的影響for del=0.1:0.2:1.1den=4*del 1;impulse(tf(num,den),10) %一階系統(tǒng)單位階脈沖應(yīng)曲線hold on, %不同時(shí)間常數(shù)的一階系統(tǒng)單位脈沖響應(yīng)曲線簇 i=i+1;end注意MATLAB中for語句的結(jié)構(gòu)。讀者可以改變不同的增益，看看圖形有何變化。2 二階系統(tǒng)（設(shè)01）設(shè)二階系統(tǒng)為二階系統(tǒng)的特征參數(shù)為固有頻率及阻尼比。當(dāng)增大，系統(tǒng)振動(dòng)頻率加快，振蕩加?。欢S著減小，系統(tǒng)振蕩加劇，振蕩峰尖銳。參考程序fanli012示出了當(dāng)0.5, =1, 2, 3, 4, 5 rad/s時(shí)的間接階躍曲線簇。參考程序fanli012: 不同固有頻率的二階系統(tǒng)的單位階躍響應(yīng)曲線（=0.5）i=1;for del=1:1:5; % 二階系統(tǒng)固有頻率遞增間隔num=del2; % 二階系統(tǒng)傳遞函數(shù)分子系數(shù)向量den=1 del del2; % 不同固有頻率的二階系統(tǒng)分母系數(shù)向量step(tf(num,den),6) %二階系統(tǒng)單位階躍響應(yīng)曲線hold on, %不同固有頻率的二階系統(tǒng)單位階躍響應(yīng)曲線簇i=i+1;endfanli012二階系統(tǒng)固有頻率對(duì)單位階躍響應(yīng)的影響參考程序fanli013示出了同一二階系統(tǒng)的單位脈沖響應(yīng)曲線簇。參考程序fanli013: 不同固有頻率的二階系統(tǒng)的單位脈沖響應(yīng)曲線（=0.5）i=1;for del=1:1:5;num=del2;den=1 del del2; % 不同固有頻率的二階系統(tǒng)分母系數(shù)向量impulse(tf(num,den),6) %二階系統(tǒng)單位脈沖響應(yīng)曲線hold on, %不同固有頻率的二階系統(tǒng)單位脈沖響應(yīng)曲線簇i=i+1;endfanli013二階系統(tǒng)固有頻率對(duì)單位脈沖響應(yīng)的影響參考程序fanlio14示出了當(dāng)1，0.1，0.3，0.5，0.7，0.9的二階系統(tǒng)的單位階躍響應(yīng)曲線簇。參考程序fanli014: 不同阻尼比的二階系統(tǒng)的單位階躍響應(yīng)曲線（=1）i=1;for del=0.1:0.2:0.9; % 二階系統(tǒng)阻尼比遞增間隔 num=1; den=1 2*del 1; % 不同阻尼比的二階系統(tǒng)分母系數(shù)向量 step(tf(num,den),30) hold on, %不同阻尼比的二階系統(tǒng)單位階躍響應(yīng)曲線簇 i=i+1;endfamli014二階系統(tǒng)阻尼比對(duì)單位階躍響應(yīng)的影響參考程序fanli015示出了同一二階系統(tǒng)當(dāng)取0.1，0.3，0.5，0.7，0.9的單位脈沖響應(yīng)曲線簇。參考程序fanli015: 不同阻尼比的二階系統(tǒng)的單位脈沖響應(yīng)曲線（=1）i=1;for del=0.1:0.2:0.9; num=1; den=1 2*del 1; impulse(tf(num,den),30) hold on, %不同阻尼比的二階系統(tǒng)單位脈沖響應(yīng)曲線簇 i=i+1;endfanli015二階系統(tǒng)阻尼比對(duì)單位脈沖響應(yīng)的影響四、帶延時(shí)環(huán)節(jié)系統(tǒng)的典型響應(yīng)設(shè)具有純延時(shí)環(huán)節(jié)的傳遞函數(shù)為計(jì)算這種系統(tǒng)的單位階躍響應(yīng)不能使用一般方法。首先應(yīng)使用Pade法對(duì)延時(shí)環(huán)節(jié)進(jìn)行近似展開，numT, dent = pade(T, 5) （表示使用5階Pade級(jí)數(shù)）得出的分子、分母系數(shù)的向量，然后將系統(tǒng)視為慣性環(huán)節(jié)與延遲環(huán)節(jié)的串聯(lián)，再求其階躍脈沖響應(yīng)。設(shè)K = 2.5，T1，程序見fanli016。參考程序fanli016: 帶延時(shí)環(huán)節(jié)系統(tǒng)的單位階躍響應(yīng)K=2.5;T=1; T1=0.5;numt,dent=pade(T,5);%5階Pade法近似延時(shí)環(huán)節(jié)syst=tf(numt,dent); %延時(shí)環(huán)節(jié)的近似傳遞函數(shù)num1=K;den1=T1 1;sys1=tf(num1,den1);%慣性環(huán)節(jié)sys=series(sys1,syst);% 帶延時(shí)環(huán)節(jié)的慣性系統(tǒng)step(sys,6),grid % 帶延時(shí)環(huán)節(jié)慣性系統(tǒng)的單位階躍響應(yīng)figure,impulse(sys,6),grid% 帶延時(shí)環(huán)節(jié)慣性系統(tǒng)的單位脈沖響應(yīng)fanli016帶延遲的一階系統(tǒng)的單位階躍響應(yīng)fanli016帶延遲的一階系統(tǒng)的單位脈沖響應(yīng)從響應(yīng)曲線上可見，具有延時(shí)環(huán)節(jié)的慣性系統(tǒng)當(dāng)加上單位階躍或單位脈沖輸入前后，其響應(yīng)初始有振蕩特性，脈沖響應(yīng)更為明顯。這種振蕩現(xiàn)象使得對(duì)這種系統(tǒng)進(jìn)行校正需要使用特殊的方法。在作圖命令后，加上格線命令grid，機(jī)器自動(dòng)給曲線圖加上格線。2-3響應(yīng)曲線的動(dòng)態(tài)分析一、MATLAB圖形的編輯在MATLAB中，用于編輯圖形的命令很多，因?yàn)榭梢暬荕ATLAB的一種強(qiáng)大而優(yōu)越的性能。從工程實(shí)用角度，僅在圖形窗口進(jìn)行編輯已經(jīng)夠用了。在圖形窗口頂部的第二排依次排有下列命令按鈕：Edit Plot（）：圖形編輯。點(diǎn)擊進(jìn)入圖形編輯功能。點(diǎn)擊選擇編輯對(duì)象。例如，曲線，標(biāo)題，縱、橫坐標(biāo)的說明等，可對(duì)選中項(xiàng)進(jìn)行編輯。Insert Text（）：插入文本。點(diǎn)擊可在圖區(qū)內(nèi)插入文字，并可選擇字體及大小等。Insert Arrow（）：插入箭頭。點(diǎn)擊可在圖區(qū)內(nèi)畫箭頭，并可對(duì)該箭頭進(jìn)行編輯。Insert Line()：插入直線。點(diǎn)擊可在圖區(qū)內(nèi)插入直線，并可對(duì)該直線進(jìn)行編輯，改變線寬及顏色。Zoom In（）：曲線放大。激活該按鈕，可以逐次放大曲線，以觀察曲線某些部分的細(xì)節(jié)，這對(duì)觀察變化劇烈的曲線部分很有幫助。Zoom Out（）：曲線縮小。激活該按鈕，可以逐次從放大狀態(tài)返回到原狀態(tài)。Rotate 3D（）：三維旋轉(zhuǎn)。激活該按鈕，在圖區(qū)內(nèi)按住鼠標(biāo)左鍵，拖曳鼠標(biāo)，原二維曲線可三維空間內(nèi)任意旋轉(zhuǎn)。此外，還有平面旋轉(zhuǎn)、三維動(dòng)畫、前后放大、平移、飛逸等。例如，將曲線加粗，以利作圖（曲線默認(rèn)寬度為0.5）。按Edit Plot 選擇曲線雙擊左鍵 出現(xiàn)線寬選擇框（line width） 選擇線寬 OK。二、對(duì)曲線進(jìn)行數(shù)據(jù)分析不激活Edit Plot按鈕（圖形編輯不使能）， 讀取曲線上任一點(diǎn)的坐標(biāo)值：鼠標(biāo)指向選擇點(diǎn)，單擊左鍵，顯示曲線編號(hào)，橫坐標(biāo)及縱坐標(biāo)值； 研究曲線擬合情況：按Tools Basic Fitting，出現(xiàn)Plot fitt（圖形擬合）窗口，選擇需要擬合的曲線，點(diǎn)擊擬合方式，即可對(duì)所選圖形進(jìn)行曲線擬合。MATLAB提供了樣條插值、保形插值（shape-preserving interpolant）及直到10階的多項(xiàng)式擬合。當(dāng)選擇多項(xiàng)式擬合時(shí)，機(jī)器會(huì)給出曲線及擬合方程。在所給出的擬合方式中，三次樣條插值（Cubic spline interpolant）及保形插值（shape-preserving interpolant）效果最好。第三章 系統(tǒng)的頻率特性系統(tǒng)的頻率特性就是，這是一個(gè)復(fù)變函數(shù)。3-1MATLAB中的復(fù)數(shù)及復(fù)變量在MATLAB中復(fù)數(shù)x的輸入方法：x = 2 + 3*i或2 + 3*j也可以s = 2 + 3i或s = 3 + 2j在符號(hào)運(yùn)算功能下，復(fù)變函數(shù)表示為（設(shè)a，b都是實(shí)數(shù)）：syms a b realx = ab*I或ab*j（*不能省略）在MATLAB中與復(fù)數(shù)有關(guān)的命令有abs(x)：求實(shí)數(shù)的絕對(duì)值及復(fù)數(shù)的模，支持符號(hào)運(yùn)算；angle(x)：求復(fù)數(shù)的相角，單位為弧度，不支持符號(hào)運(yùn)算。conj(x)：求復(fù)數(shù)的共軛復(fù)數(shù)，支持符號(hào)運(yùn)算；imag(x)：求復(fù)數(shù)的虛部，支持符號(hào)運(yùn)算；real(x)：求復(fù)數(shù)的實(shí)部，支持符號(hào)運(yùn)算。復(fù)數(shù)的運(yùn)算見參考程序fanli017。參考程序fanli017:復(fù)數(shù)運(yùn)算（含符號(hào)功能）syms a b c d real%符號(hào)均為實(shí)變量X1=2-3i X2=a+c/b*jX3=c/a-2*d*jS1=X2*X3/X1 %復(fù)數(shù)的乘除S2=simple(abs(S1) %復(fù)數(shù)的模S3=simple(real(S1) %復(fù)數(shù)的實(shí)部S4=simple(imag(S1) %復(fù)數(shù)的虛部S5=simple(conj(S1) %復(fù)數(shù)的共軛S6=simple(expand(S1*S5) %復(fù)數(shù)與其共軛的積S7=simple(S32+S42) %復(fù)數(shù)實(shí)部、虛部的平方和S8=simple(S6-S7) %檢驗(yàn)復(fù)數(shù)運(yùn)算的正確性fanli018給出了系統(tǒng)的頻率特性。參考程序fanli018:系統(tǒng)頻率特性的解析表達(dá)syms w real %定義頻率w 為實(shí)型變量g=50*(0.6*w*i+1)/(i*w)2/(4*w*j+1) %系統(tǒng)頻率特性g1=simple(conj(g) %系統(tǒng)頻率特性的復(fù)共軛gr=simple(real(g) %系統(tǒng)實(shí)頻特性gi=simple(imag(g) %系統(tǒng)虛頻特性ga=simple(abs(g) %系統(tǒng)幅頻特性（模）gb=simple(sqrt(gr2+gi2) %系統(tǒng)幅頻特性的另一計(jì)算方法gc=simple(ga-gb) %檢驗(yàn)系統(tǒng)幅頻特性的兩種計(jì)算方法3-2頻率特性的Nyquist圖將系統(tǒng)傳遞函數(shù)中的復(fù)變數(shù)s用純虛數(shù)j（為角頻率，rad/s）代替，即得到系統(tǒng)的頻率特性，又稱為諧波傳遞函數(shù)，有三種表示方法：式中：實(shí)頻特性：虛頻特性：幅頻特性：相頻特性在研究控制系統(tǒng)的特性時(shí)，采用圖示方法比采用解析方法直觀、簡(jiǎn)單，便于對(duì)動(dòng)態(tài)特性進(jìn)行更深入的研究。以下討論設(shè)系統(tǒng)傳遞函數(shù)為sys。一、Nyquist圖1 Nyquist 圖的有關(guān)命令nyquist(sys)直接返回系統(tǒng)sys的Nyquist圖，機(jī)器自己確定頻率范圍（通常為- +）nyquist(sys, w)返回在規(guī)定頻率（rad/s）范圍內(nèi)的Nyquist圖，頻率范圍的格式為w = wmin, wmax (0wminwmin)nyquist(sys1,sys2,sysn)nyquist(sys1,sys2,sysn,w)在同一幅度圖上畫出幾個(gè)系統(tǒng)的Nyquist圖，可以由機(jī)器確定頻率范圍，也可以自選頻率范圍。re,im,w=nyquist(sys)分別返回系統(tǒng)sys對(duì)應(yīng)頻率w的實(shí)部（re）及虛部（im）值。讀取實(shí)部，虛部及對(duì)應(yīng)頻率點(diǎn)的方法：re(1,1,k)im(1,1,k)w(k) (rad/s)其中，k為機(jī)器計(jì)算Nyquist響應(yīng)的序列編號(hào)，必須是正整數(shù)而且應(yīng)當(dāng)在計(jì)算范圍內(nèi)?？梢允褂妹顂ize(w)查看。k值隨問題的簡(jiǎn)、繁而變化。這是一個(gè)非常有用的命令，當(dāng)然也可以采用直接在Nyquist圖上用鼠標(biāo)讀取的方法。re,im = nyquist(sys,w)返回在某一頻率w（rad/s）或某一頻率范圍w = wmin, wmax 系統(tǒng)的實(shí)部及虛部值。當(dāng)選定為頻率范圍時(shí)，調(diào)出實(shí)部，虛部及對(duì)應(yīng)頻率值的方法同前一指令。例3.1 分析系統(tǒng)的頻率特性。該傳遞函數(shù)是一個(gè)零極點(diǎn)及有理分式的混合形式。首先將子系統(tǒng)變換成零極點(diǎn)增益形式z1，p1，k1，再將系統(tǒng)變換成零極點(diǎn)增益模型（z，p，k）再作Nyquist圖。使用re，im，wnyquist(sys)語句，可計(jì)算出若干頻率點(diǎn)所對(duì)應(yīng)的實(shí)部、虛部值，而且在w(25)=0.8802rad/sw(26)=0.8951rad/s之間，曲線將從第三象限穿過負(fù)實(shí)軸進(jìn)入第二象限，在w(0.8802，0.8951)之間使用語句re，im=nyquist(sys，w)0.8802，w，0.8951進(jìn)行更細(xì)致的計(jì)算，在需要的精度內(nèi)求出穿越負(fù)實(shí)軸的參數(shù)。參考程序見fanli019。參考程序fanli019:用Nyquist圖分析系統(tǒng)頻率特性num=5;den=1 0.2 1;%子系統(tǒng)的傳遞函數(shù)Z1,P1,K1=tf2zp(num,den) %子系統(tǒng)的極點(diǎn)Z=-2.5;-2 %系統(tǒng)的零點(diǎn)P=P1;1/1.2;-1/1.4 %系統(tǒng)的極點(diǎn)K=num %系統(tǒng)的增益sys=zpk(Z,P,K) %系統(tǒng)的零、極點(diǎn)、增益模型nyquist(sys) %系統(tǒng)的Nyquist圖re,im,w=nyquist(sys) %計(jì)算系統(tǒng)若干頻率點(diǎn)所對(duì)應(yīng)的實(shí)部與虛部值w(25) %w(25)=0.8802,曲線處于第三象限靠近負(fù)實(shí)軸的頻率點(diǎn)w(26) %w(26)=0.8951,曲線處于第二象限靠近負(fù)實(shí)軸的頻率點(diǎn)re1,im1=nyquist(sys,0.8802:0.00005:0.8951) %在（0.8802, 0.8951）(rad/s)范圍內(nèi)，以0.00005rad/s為間隔計(jì)算實(shí)部及虛部值fanli019的Nyquist圖2 帶延時(shí)環(huán)節(jié)系統(tǒng)的Nyquist圖對(duì)帶有延時(shí)環(huán)節(jié)的系統(tǒng)，由于不能表示成傳遞函數(shù)的有理分式，所以MATLAB中沒有直接命令可用。但由于Nyquist圖實(shí)際上是以頻率為參變量，橫坐標(biāo)為實(shí)頻特性，縱坐標(biāo)為虛頻特性的直角坐標(biāo)系圖。所以，可以通過計(jì)算實(shí)頻特性與虛頻特性，再使用plot(x，y)作出Nyquist圖。例3.2 求作具有延時(shí)環(huán)節(jié)的系統(tǒng)的Nyquist圖。程序見fanli020。參考程序fanli020:帶延時(shí)環(huán)節(jié)系統(tǒng)的Nyquist圖w=0:0.01:30;%設(shè)定頻率變化范圍gjw=5*exp(-1.5*j*w)./(1+j*2*w);%計(jì)算頻率特性x=real(gjw); %實(shí)頻特性y=imag(gjw); %虛頻特性plot(x,y),grid %帶延時(shí)環(huán)節(jié)系統(tǒng)的Nyquist圖fanli020帶延時(shí)環(huán)節(jié)的慣性系統(tǒng)的Nyquist圖注意，程序中兩項(xiàng)相除時(shí)采用了對(duì)應(yīng)元素相除的方法，格式為：分子./分母。3-3系統(tǒng)的Bode圖Bode圖是系統(tǒng)的對(duì)數(shù)頻率特性曲線，由兩幅圖構(gòu)成：對(duì)數(shù)幅頻特性圖和對(duì)數(shù)相頻特性圖?？v坐標(biāo)以線性分度，幅值，以分貝表示；以度表示。兩張圖共用橫軸。橫坐標(biāo)為頻率，采用對(duì)數(shù)分度（），但習(xí)慣上仍標(biāo)真數(shù)。因而，橫坐標(biāo)是按十倍頻程均勻分度的。對(duì)于寫成典型環(huán)節(jié)形式的傳遞函數(shù)，其對(duì)應(yīng)頻率特性為幅頻特性相頻特性一、Bode圖的有關(guān)命令Bode圖的命令格式與Nyquist圖的命令格式十分相似，其作用對(duì)象仍然是由有理分式或零極點(diǎn)模型構(gòu)成的系統(tǒng)傳遞函數(shù)sys。Bode圖的有關(guān)命令如下：Bode(sys) 直接返回bode圖bode(sys, w) 在設(shè)定頻率范圍w=wmin, wmax或w=wmin: w:wmax內(nèi)的Bode圖。Bode(sys1，sys2，sysN)Bode(sys1，sys2，sysN，w)在同一幅圖上作幾個(gè)系統(tǒng)的Bode圖。mag，phase，wbode(sys)其中，w表示一系列計(jì)算點(diǎn)的頻率（rad/s）phase：對(duì)應(yīng)計(jì)算頻率點(diǎn)的相角（度）mag：對(duì)應(yīng)計(jì)算頻率點(diǎn)的幅值（非分貝）若要用分貝表示mag(db)=20*log10(mag)讀取mag，phase，w的方法與Nyquistg命令中讀取實(shí)部、虛部及頻率的方法相同。mag，phasebode(sys，w)返回某一頻率點(diǎn)w，或某一頻率范圍w =wmin, wmax或w = wmin: w:wmax的幅值及相角。當(dāng)選它為頻率范圍時(shí)，讀取方法相同。例3.3 作系統(tǒng)的Bode圖。注意該系統(tǒng)是一個(gè)II型系統(tǒng)。程序見fanli021。參考程序fanli021:型系統(tǒng)的Bode圖num=1 0.1 7.5;den=1 0.12 9 0 0;sys=tf(num,den)bode(sys,1,5),gridmag,phase,w=bode(sys); %計(jì)算系統(tǒng)若干頻率點(diǎn)所對(duì)應(yīng)的幅值與相角mag,phase=bode(sys,2.8479:0.001:2.8719) %計(jì)算小頻率范圍的幅值與相角fanli021的Bode圖程序的最后一句，用以計(jì)算相頻特性曲線峰點(diǎn)所對(duì)應(yīng)的相角，頻率及幅值。例3.4對(duì)比Nyquist圖及Bode圖。fanli022系統(tǒng)的Nyquist圖fanli022系統(tǒng)的Bode圖參考程序fanli022:系統(tǒng)Nyquist 圖與Bode圖的對(duì)比num=50;den=0.01 0.14 1;Z1,P1,K1=tf2zp(num,den) %子系統(tǒng)的極點(diǎn)Z=-2.5;-1/0.2P=P1;1/0.9;-1/0.3K=numsys=zpk(Z,P,K) %系統(tǒng)的零極點(diǎn)增益模型bode(sys),grid %系統(tǒng)的Bode圖figure,nyquist(sys) %系統(tǒng)的Nyquist圖mag,phase,w=bode(sys);%系統(tǒng)幅值及相角的計(jì)算二、典型環(huán)節(jié)的Nyquist圖和Bode圖根據(jù)典型環(huán)節(jié)的傳遞函數(shù)，使用前述的Nyquist圖及Bode圖的相關(guān)命令，很容易作出它們的這兩種頻率特性圖。讀者應(yīng)重視典型環(huán)節(jié)Nyquist圖及Bode圖的作法及走向，以及它們所反映的系統(tǒng)動(dòng)態(tài)特性，研究在原有系統(tǒng)增加或減少一個(gè)典型環(huán)節(jié)時(shí)對(duì)系統(tǒng)的特性曲線所引起的變化以及這種變化的物理及工程含義。例3.5：對(duì)比及的Nyquist圖及Bode圖，特別注意Bode圖中的相角變化，程序見fanli023。參考程序fanli023:系統(tǒng)Nyquist 圖與Bode圖的對(duì)比num1=1 -2;num2=-1 2;den=1 0.1 4;sys1=tf(num1,den) %系統(tǒng)1sys2=tf(num2,den) %系統(tǒng)2bode(sys1),grid %系統(tǒng)1的Bode圖figure,bode(sys2),grid %系統(tǒng)2的Bode圖figure,nyquist(sys1),grid %系統(tǒng)1的Nyquist圖figure,nyquist(sys2),grid %系統(tǒng)2的Nyquist圖fanli023:系統(tǒng)1的Bode圖fanli023:系統(tǒng)2的Bode圖fanli023:系統(tǒng)1的Nyquist圖fanli023:系統(tǒng)2的Nyquist圖Nyquist圖中的虛線是等幅值圓。第四章系統(tǒng)的校正4-1系統(tǒng)的穩(wěn)定性分析系統(tǒng)穩(wěn)定的基本條件是其極點(diǎn)具有負(fù)實(shí)部。判定系統(tǒng)穩(wěn)定性的主要方法有1. 若已知系統(tǒng)的特征多項(xiàng)式 可以用求解方程的辦法，求出特征根，常用如下指令：solve(f(s)返回特征根的值。Solve命令必須在符號(hào)運(yùn)算功能下運(yùn)行。2. 設(shè)特征多項(xiàng)式為其降排列的系數(shù)向量為CC1C2Cn-1 則指令C1roots(C)返回特征根的數(shù)值。若使用命令poly(C1)則返回原多項(xiàng)式系數(shù)向量，即roots與poly互為逆操作。Roots命令不需要符號(hào)功能支持。若給出的是開環(huán)傳遞函數(shù)，需要判斷對(duì)應(yīng)的單位負(fù)反饋系統(tǒng)的穩(wěn)定性，必須用feedback(sys，1)命令求出對(duì)應(yīng)的閉環(huán)系統(tǒng)，再求閉環(huán)特征根。2 若已知傳遞函數(shù)。可以用tf2zp轉(zhuǎn)換成零、極點(diǎn)模型判斷其穩(wěn)定性。3 根據(jù)Nyquist穩(wěn)定性判據(jù)，作出開環(huán)系統(tǒng)的Nyquist圖，通過該系統(tǒng)圍繞(-1，j0)點(diǎn)的情況判斷相應(yīng)的單位負(fù)反饋閉環(huán)系統(tǒng)的穩(wěn)定性。4 Bode判據(jù)，通過開環(huán)系統(tǒng)的Bode圖判斷相應(yīng)的單位負(fù)反饋閉環(huán)系統(tǒng)穩(wěn)定性。對(duì)于復(fù)雜系統(tǒng)，可以綜合使用上述幾種方法，還可用step(sys)命令檢驗(yàn)系統(tǒng)的單位階躍響應(yīng)是否收斂來映證判定結(jié)果的正確性。例4.1在系統(tǒng)的Nyquist圖上，疊加上單位圓，以便使用Nyquist判據(jù)。設(shè)開環(huán)系統(tǒng)作出其Nyquist圖并疊加上單位圓圖。由圖可見，該開環(huán)系統(tǒng)的Nyquist圖圍繞（1，j0）點(diǎn)逆時(shí)針及順時(shí)針各一次，因而相應(yīng)的閉環(huán)系統(tǒng)穩(wěn)定。程序見fanli024。參考程序fanli024: Nyquist圖與單位圓圖w2=linspace(0