控制系統(tǒng)仿真 第五章控制系統(tǒng)CAD學(xué)習(xí)資料

自動(dòng)控制系統(tǒng)仿真Thesimulationofautomaticcontrolsystem第五章控制系統(tǒng)CAD第五章控制系統(tǒng)CAD控制系統(tǒng)根軌跡分析控制系統(tǒng)穩(wěn)定性分析控制系統(tǒng)時(shí)域分析控制系統(tǒng)頻域分析控制系統(tǒng)設(shè)計(jì)方法現(xiàn)代控制理論CAD第13次控制系統(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)定的。若連續(xù)時(shí)間系統(tǒng)的全部零極點(diǎn)都位于S左半平面；或若離散時(shí)間系統(tǒng)的全部零極點(diǎn)都位于Z平面單位圓內(nèi)，則系統(tǒng)是最小相位系統(tǒng)。一、系統(tǒng)穩(wěn)定及最小相位系統(tǒng)判據(jù)控制系統(tǒng)穩(wěn)定性分析2、直接判別MATLAB提供了直接求取系統(tǒng)所有零極點(diǎn)的函數(shù)，因此可以直接根據(jù)零極點(diǎn)的分布情況對(duì)系統(tǒng)的穩(wěn)定性及是否為最小相位系統(tǒng)進(jìn)行判斷。二、系統(tǒng)穩(wěn)定及最小相位系統(tǒng)的判別方法1、間接判別（工程方法）勞斯判據(jù)：勞斯表中第一列各值嚴(yán)格為正，則系統(tǒng)穩(wěn)定，如果勞斯表第一列中出現(xiàn)小于零的數(shù)值，系統(tǒng)不穩(wěn)定。赫爾維茨判據(jù)：當(dāng)且僅當(dāng)由系統(tǒng)分母多項(xiàng)式構(gòu)成的赫爾維茨矩陣為正定矩陣時(shí)，系統(tǒng)穩(wěn)定?？刂葡到y(tǒng)穩(wěn)定性分析例exp5_1.m已知某系統(tǒng)的模型如右所示：要求判斷系統(tǒng)的穩(wěn)定性及系統(tǒng)是否為最小相位系統(tǒng)。%n1>0,實(shí)部大于零的零點(diǎn)個(gè)數(shù)大于0，則為非最小相位系統(tǒng)%n2>0,實(shí)部大于零的極點(diǎn)個(gè)數(shù)大于0，則系統(tǒng)不穩(wěn)定，求出不穩(wěn)定極點(diǎn)ii=find(條件式)用來(lái)求取滿足條件的向量的下標(biāo)向量，以列向量表示。例如exp5_1.m中的條件式為jj=find(real(p)>0)，其含義就是找出極點(diǎn)向量p中滿足實(shí)部的值大于0的所有元素下標(biāo)，并將結(jié)果返回到j(luò)j向量中去。這樣如果找到了實(shí)部大于0的極點(diǎn)，則會(huì)將該極點(diǎn)的序號(hào)返回到j(luò)j下。如果最終的結(jié)果里jj的元素個(gè)數(shù)大于0，則認(rèn)為找到了不穩(wěn)定極點(diǎn)，因而給出系統(tǒng)不穩(wěn)定的提示，若產(chǎn)生的jj向量的元素個(gè)數(shù)為0，則認(rèn)為沒(méi)有找到不穩(wěn)定的極點(diǎn)，因而得出系統(tǒng)穩(wěn)定的結(jié)論?？刂葡到y(tǒng)穩(wěn)定性分析

例exp5_2.m系統(tǒng)模型如下所示，判斷系統(tǒng)的穩(wěn)定性，以及系統(tǒng)是否為最小相位系統(tǒng)。第五章控制系統(tǒng)CAD控制系統(tǒng)根軌跡分析控制系統(tǒng)穩(wěn)定性分析控制系統(tǒng)時(shí)域分析控制系統(tǒng)頻域分析控制系統(tǒng)設(shè)計(jì)方法現(xiàn)代控制理論CAD控制系統(tǒng)時(shí)域分析一個(gè)動(dòng)態(tài)系統(tǒng)的性能常用典型輸入作用下的響應(yīng)來(lái)描述。響應(yīng)是指零初始值條件下某種典型的輸入函數(shù)作用下對(duì)象的響應(yīng)，控制系統(tǒng)常用的輸入函數(shù)為單位階躍函數(shù)和脈沖激勵(lì)函數(shù)（即沖激函數(shù)）。在MATLAB的控制系統(tǒng)工具箱中提供了求取這兩種輸入下系統(tǒng)響應(yīng)的函數(shù)。一、時(shí)域分析的一般方法求取系統(tǒng)單位階躍響應(yīng)：step()求取系統(tǒng)的沖激響應(yīng)：impulse()控制系統(tǒng)時(shí)域分析求取系統(tǒng)單位階躍響應(yīng)：step()求取系統(tǒng)的沖激響應(yīng)：impulse()控制系統(tǒng)時(shí)域分析1、step()函數(shù)的用法控制系統(tǒng)時(shí)域分析2、impulse()函數(shù)的用法控制系統(tǒng)時(shí)域分析控制系統(tǒng)時(shí)域分析控制系統(tǒng)時(shí)域分析3、仿真時(shí)間t的選擇：對(duì)于典型二階系統(tǒng)根據(jù)其響應(yīng)時(shí)間的估算公式可以確定。對(duì)于高階系統(tǒng)往往其響應(yīng)時(shí)間很難估計(jì)，一般采用試探的方法，把t選大一些，看看響應(yīng)曲線的結(jié)果，最后再確定其合適的仿真時(shí)間。一般來(lái)說(shuō)，先不指定仿真時(shí)間，由MATLAB自己確定，然后根據(jù)結(jié)果，最后確定合適的仿真時(shí)間。在指定仿真時(shí)間時(shí)，步長(zhǎng)的不同會(huì)影響到輸出曲線的光滑程度，一般不易取太大?？刂葡到y(tǒng)時(shí)域分析例exp5_6_.m控制系統(tǒng)時(shí)域分析二、常用時(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)上反映出來(lái)。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)lsim()：連續(xù)系統(tǒng)對(duì)任意輸入的響應(yīng)

gensig()

：產(chǎn)生信號(hào)對(duì)于離散系統(tǒng)只需在連續(xù)系統(tǒng)對(duì)應(yīng)函數(shù)前加d就可以，如dstep，dimpulse等。它們的調(diào)用格式與step、impulse類(lèi)似，可以通過(guò)help命令來(lái)察看自學(xué)?？刂葡到y(tǒng)時(shí)域分析例：繪制下列系統(tǒng)的零輸入響應(yīng)，初始狀態(tài)x0=[1;0]。控制系統(tǒng)時(shí)域分析例：繪制下列系統(tǒng)的方波響應(yīng)，其中方波的周期為4秒，持續(xù)時(shí)間為10秒，采樣周期為0.01秒?？刂葡到y(tǒng)時(shí)域分析三、時(shí)域分析應(yīng)用實(shí)例MATLAB的step()和impulse()函數(shù)本身可以處理多輸入多輸出的情況，因此編寫(xiě)MATLAB程序并不因?yàn)橄到y(tǒng)輸入輸出的增加而變得復(fù)雜?？刂葡到y(tǒng)時(shí)域分析控制系統(tǒng)時(shí)域分析pos=40,tp=0.8例exp5_8.m控制系統(tǒng)時(shí)域分析alph=0.5,wn=5pos=16.1301tp=0.3313tr=1.6123ts2=0.7288%退出循環(huán)%退出循環(huán)%退出循環(huán)第五章控制系統(tǒng)CAD控制系統(tǒng)根軌跡分析控制系統(tǒng)穩(wěn)定性分析控制系統(tǒng)時(shí)域分析控制系統(tǒng)頻域分析控制系統(tǒng)設(shè)計(jì)方法現(xiàn)代控制理論CAD控制系統(tǒng)頻域分析頻率響應(yīng)是指系統(tǒng)對(duì)正弦輸入信號(hào)的穩(wěn)態(tài)響應(yīng)，從頻率響應(yīng)中可以得出帶寬、增益、轉(zhuǎn)折頻率、閉環(huán)穩(wěn)定性等系統(tǒng)特征。頻率特性是指系統(tǒng)在正弦信號(hào)作用下，穩(wěn)態(tài)輸出與輸入之比對(duì)頻率的關(guān)系特性。頻率特性函數(shù)與傳遞函數(shù)有直接的關(guān)系，記為：一、頻域分析的一般方法求取系統(tǒng)對(duì)數(shù)頻率特性圖（波特圖）：bode()求取系統(tǒng)奈奎斯特圖（幅相曲線圖或極坐標(biāo)圖）：nyquist()頻域分析法是應(yīng)用頻率特性研究控制系統(tǒng)的一種典型方法。采用這種方法可直觀地表達(dá)出系統(tǒng)的頻率特性，分析方法比較簡(jiǎn)單，物理概念比較明確，對(duì)于諸如防止結(jié)構(gòu)諧振、抑制噪聲、改善系統(tǒng)穩(wěn)定性和暫態(tài)性能等問(wèn)題，都可以從系統(tǒng)的頻率特性上明確地看出其物理實(shí)質(zhì)和解決途經(jīng)。通常將頻率特性用曲線的形式進(jìn)行表示，包括對(duì)數(shù)頻率特性曲線和幅相頻率特性曲線簡(jiǎn)稱(chēng)幅相曲線，MATLAB提供了繪制這兩種曲線的函數(shù)?？刂葡到y(tǒng)頻域分析1、對(duì)數(shù)頻率特性圖（波特圖）bode()2、奈奎斯特圖（幅相頻率特性圖）nyquist()Figure(1)Figure(2)控制系統(tǒng)頻域分析二、常用頻域分析函數(shù)MATLAB除了提供前面介紹的基本頻域分析函數(shù)外，還提供了大量在工程實(shí)際中廣泛應(yīng)用的庫(kù)函數(shù)，由這些函數(shù)可以求得系統(tǒng)的各種頻率響應(yīng)曲線和特征值。如：margin：求幅值裕度和相角裕度及對(duì)應(yīng)的轉(zhuǎn)折頻率freqs：模擬濾波器特性nichols：求連續(xù)系統(tǒng)的尼科爾斯頻率響應(yīng)曲線（即對(duì)數(shù)幅相曲線）ngrid：尼科爾斯方格圖控制系統(tǒng)頻域分析margin()函數(shù)margin函數(shù)可以從頻率響應(yīng)數(shù)據(jù)中計(jì)算出幅值裕度、相角裕度以及對(duì)應(yīng)的頻率。幅值裕度和相角裕度是針對(duì)開(kāi)環(huán)SISO系統(tǒng)而言，它指示出系統(tǒng)閉環(huán)時(shí)的相對(duì)穩(wěn)定性。當(dāng)不帶輸出變量引用時(shí)，margin可在當(dāng)前圖形窗口中繪制出帶有裕量及相應(yīng)頻率顯示的Bode圖，其中幅值裕度以分貝為單位。幅值裕度是在相角為-180度處使開(kāi)環(huán)增益為1的增益量，如在-180度相頻處的開(kāi)環(huán)增益為g，則幅值裕度為1/g；若用分貝值表示幅值裕度，則等于：-20*log10(g)。類(lèi)似地，相角裕度是當(dāng)開(kāi)環(huán)增益為1.0時(shí)，相應(yīng)的相角與180度角的和。margin(mag,phase,w)：由bode指令得到的幅值mag（不是以dB為單位）、相角phase及角頻率w矢量繪制出帶有裕量及相應(yīng)頻率顯示的bode圖。margin(num,den)：可計(jì)算出連續(xù)系統(tǒng)傳遞函數(shù)表示的幅值裕度和相角裕度并繪制相應(yīng)波特圖。margin(a,b,c,d):可以計(jì)算出連續(xù)狀態(tài)空間系統(tǒng)表示的幅值裕度和相角裕度并繪制相應(yīng)波特圖。[gm,pm,wcg,wcp]=margin(mag,phase,w)：由幅值mag（不是以dB為單位）、相角phase及角頻率w矢量計(jì)算出系統(tǒng)幅值裕度和相角裕度及相應(yīng)的相角交界頻率wcg、截止頻率wcp，而不直接繪出Bode圖曲線?？刂葡到y(tǒng)頻域分析freqs()函數(shù)freqs用于計(jì)算由矢量a和b構(gòu)成的模擬濾波器H(s)=B(s)/A(s)的幅頻響應(yīng)。h=freqs(b,a,w)用于計(jì)算模擬濾波器的幅頻響應(yīng)，其中實(shí)矢量w用于指定頻率值，返回值h為一個(gè)復(fù)數(shù)行向量，要得到幅值必須對(duì)它取絕對(duì)值，即求模。[h,w]=freqs(b,a)自動(dòng)設(shè)定200個(gè)頻率點(diǎn)來(lái)計(jì)算頻率響應(yīng)，這200個(gè)頻率值記錄在w中。[h,w]=freqs(b,a,n)設(shè)定n個(gè)頻率點(diǎn)計(jì)算頻率響應(yīng)。不帶輸出變量的freqs函數(shù)，將在當(dāng)前圖形窗口中繪制出幅頻和相頻曲線，其中幅相曲線對(duì)縱坐標(biāo)與橫坐標(biāo)均為對(duì)數(shù)分度?？刂葡到y(tǒng)頻域分析三、頻域分析應(yīng)用實(shí)例Nyquist曲線是根據(jù)開(kāi)環(huán)頻率特性在復(fù)平面上繪出的幅相軌跡，根據(jù)開(kāi)環(huán)的Nyquist曲線，可以判斷閉環(huán)系統(tǒng)的穩(wěn)定性。系統(tǒng)穩(wěn)定的充要條件為：Nyquist曲線按逆時(shí)針包圍臨界點(diǎn)(-1,j0)的圈數(shù)R，等于開(kāi)環(huán)傳遞函數(shù)位于s右半平面的極點(diǎn)數(shù)P，否則閉環(huán)系統(tǒng)不穩(wěn)定，閉環(huán)正實(shí)部特征根個(gè)數(shù)Z=P-R。若剛好過(guò)臨界點(diǎn)，則系統(tǒng)臨界穩(wěn)定。例exp5_15.m線性時(shí)不變系統(tǒng)如下所示：要求繪制系統(tǒng)的波特圖和奈奎斯特圖，判斷系統(tǒng)穩(wěn)定性，如果系統(tǒng)穩(wěn)定，求出系統(tǒng)穩(wěn)定裕度，并繪制系統(tǒng)的單位沖激響應(yīng)以驗(yàn)證判斷結(jié)論。控制系統(tǒng)頻域分析Pade函數(shù)可以近似表示延時(shí)環(huán)節(jié)e^(-st)，它的調(diào)用格式為：(num,den)=pade(t,n)，產(chǎn)生最佳逼近時(shí)延t秒的n階傳遞函數(shù)形式。(a,b,c,d)=pade(t,n)，則產(chǎn)生的是n階SISO的狀態(tài)空間模型。10G(s)R(s)C(s)__例exp5_17.m系統(tǒng)結(jié)構(gòu)圖如下所示，試用nyquist頻率曲線判斷系統(tǒng)的穩(wěn)定性。第五章控制系統(tǒng)CAD控制系統(tǒng)根軌跡分析控制系統(tǒng)穩(wěn)定性分析控制系統(tǒng)時(shí)域分析控制系統(tǒng)頻域分析控制系統(tǒng)設(shè)計(jì)方法現(xiàn)代控制理論CAD控制系統(tǒng)根軌跡分析所謂根軌跡是指，當(dāng)開(kāi)環(huán)系統(tǒng)某一參數(shù)從零變到無(wú)窮大時(shí)，閉環(huán)系統(tǒng)特征方程的根在s平面上的軌跡。一般來(lái)說(shuō)，這一參數(shù)選作開(kāi)環(huán)系統(tǒng)的增益K，而在無(wú)零極點(diǎn)對(duì)消時(shí)，閉環(huán)系統(tǒng)特征方程的根就是閉環(huán)傳遞函數(shù)的極點(diǎn)。根軌跡分析方法是分析和設(shè)計(jì)線性定?？刂葡到y(tǒng)的圖解方法，使用十分簡(jiǎn)便。利用它可以對(duì)系統(tǒng)進(jìn)行各種性能分析。一、根軌跡分析方法的概念控制系統(tǒng)根軌跡分析（1）穩(wěn)定性當(dāng)開(kāi)環(huán)增益K從零到無(wú)窮大變化時(shí)，圖中的根軌跡不會(huì)越過(guò)虛軸進(jìn)入右半s平面，因此這個(gè)系統(tǒng)對(duì)所有的K值都是穩(wěn)定的。如果根軌跡越過(guò)虛軸進(jìn)入右半s平面，則其交點(diǎn)的K值就是臨界穩(wěn)定開(kāi)環(huán)增益。（2）穩(wěn)態(tài)性能開(kāi)環(huán)系統(tǒng)在坐標(biāo)原點(diǎn)有一個(gè)極點(diǎn)，因此根軌跡上的K值就是靜態(tài)速度誤差系數(shù)，如果給定系統(tǒng)的穩(wěn)態(tài)誤差要求，則可由根軌跡確定閉環(huán)極點(diǎn)容許的范圍。（3）動(dòng)態(tài)性能當(dāng)0<K<0.5時(shí)，所有閉環(huán)極點(diǎn)位于實(shí)軸上，系統(tǒng)為過(guò)阻尼系統(tǒng)，單位階躍響應(yīng)為非周期過(guò)程；當(dāng)K=0.5時(shí)，閉環(huán)兩個(gè)極點(diǎn)重合，系統(tǒng)為臨界阻尼系統(tǒng)，單位階躍響應(yīng)仍為非周期過(guò)程，但速度更快；當(dāng)K>0.5時(shí)，閉環(huán)極點(diǎn)為復(fù)數(shù)極點(diǎn)，系統(tǒng)為欠阻尼系統(tǒng)，單位階躍響應(yīng)為阻尼振蕩過(guò)程，且超調(diào)量與K成正比?？刂葡到y(tǒng)根軌跡分析二、根軌跡分析函數(shù) 通常來(lái)說(shuō)，繪制系統(tǒng)的根軌跡是很繁瑣的事情，因此在教科書(shū)中介紹的是一種按照一定規(guī)則進(jìn)行繪制的概略根軌跡。在MATLAB中，專(zhuān)門(mén)提供了繪制根軌跡的有關(guān)函數(shù)。pzmap：繪制線性系統(tǒng)的零極點(diǎn)圖rlocus：求系統(tǒng)根軌跡。rlocfind：計(jì)算給定一組根的根軌跡增益。sgrid：在連續(xù)系統(tǒng)根軌跡圖和零極點(diǎn)圖中繪制出阻尼系數(shù)和自然頻率柵格。控制系統(tǒng)根軌跡分析1、零極點(diǎn)圖繪制MATLAB提供了函數(shù)pzmap()來(lái)繪制系統(tǒng)的零極點(diǎn)圖，其用法如下：[p,z]=pzmap(a,b,c,d)：返回狀態(tài)空間描述系統(tǒng)的極點(diǎn)矢量和零點(diǎn)矢量，而不在屏幕上繪制出零極點(diǎn)圖。[p,z]=pzmap(num,den)：返回傳遞函數(shù)描述系統(tǒng)的極點(diǎn)矢量和零點(diǎn)矢量，而不在屏幕上繪制出零極點(diǎn)圖。pzmap(a,b,c,d)或pzmap(num,den)：不帶輸出參數(shù)項(xiàng)，則直接在s復(fù)平面上繪制出系統(tǒng)對(duì)應(yīng)的零極點(diǎn)位置，極點(diǎn)用×表示，零點(diǎn)用o表示。pzmap(p,z)：根據(jù)系統(tǒng)已知的零極點(diǎn)列向量或行向量直接在s復(fù)平面上繪制出對(duì)應(yīng)的零極點(diǎn)位置，極點(diǎn)用×表示，零點(diǎn)用o表示?？刂葡到y(tǒng)根軌跡分析2、根軌跡圖繪制MATLAB提供了函數(shù)rlocus()來(lái)繪制系統(tǒng)的根軌跡圖，其用法如下：rlocus(a,b,c,d)或者rlocus(num,den)：根據(jù)SISO開(kāi)環(huán)系統(tǒng)的狀態(tài)空間描述模型和傳遞函數(shù)模型，直接在屏幕上繪制出系統(tǒng)的根軌跡圖。開(kāi)環(huán)增益的值從零到無(wú)窮大變化。rlocus(a,b,c,d,k)或rlocus(num,den,k)：通過(guò)指定開(kāi)環(huán)增益k的變化范圍來(lái)繪制系統(tǒng)的根軌跡圖。r=rlocus(num,den,k)或者[r,k]=rlocus(num,den)：不在屏幕上直接繪出系統(tǒng)的根軌跡圖，而根據(jù)開(kāi)環(huán)增益變化矢量k，返回閉環(huán)系統(tǒng)特征方程1＋k*num(s)/den(s)=0的根r，它有l(wèi)ength(k)行，length(den)-1列，每行對(duì)應(yīng)某個(gè)k值時(shí)的所有閉環(huán)極點(diǎn)?；蛘咄瑫r(shí)返回k與r。若給出傳遞函數(shù)描述系統(tǒng)的分子項(xiàng)num為負(fù)，則利用rlocus函數(shù)繪制的是系統(tǒng)的零度根軌跡。（正反饋系統(tǒng)或非最小相位系統(tǒng)）控制系統(tǒng)根軌跡分析3、rlocfind()函數(shù)MATLAB提供了函數(shù)rlocfind()來(lái)找出給定的一組根（閉環(huán)極點(diǎn)）對(duì)應(yīng)的根軌跡增益。其用法如下：[k,p]=rlocfind(a,b,c,d)或者[k,p]=rlocfind(num,den)它要求在屏幕上先已經(jīng)繪制好有關(guān)的根軌跡圖。然后，此命令將產(chǎn)生一個(gè)光標(biāo)以用來(lái)選擇希望的閉環(huán)極點(diǎn)。命令執(zhí)行結(jié)果：k為對(duì)應(yīng)選擇點(diǎn)處根軌跡開(kāi)環(huán)增益；p為此點(diǎn)處的系統(tǒng)閉環(huán)特征根。

不帶輸出參數(shù)項(xiàng)[k,p]時(shí)，同樣可以執(zhí)行，只是此時(shí)只將k的值返回到缺省變量ans中。4、sgrid()函數(shù)sgrid：在現(xiàn)存的屏幕根軌跡或零極點(diǎn)圖上繪制出自然振蕩頻率wn、阻尼比矢量z對(duì)應(yīng)的格線。sgrid(‘new’)：是先清屏，再畫(huà)格線。sgrid(z,wn)：則繪制由用戶指定的阻尼比矢量z、自然振蕩頻率wn的格線。控制系統(tǒng)根軌跡分析三、根軌跡分析應(yīng)用實(shí)例控制系統(tǒng)根軌跡分析例exp5_22.m某開(kāi)環(huán)系統(tǒng)傳遞函數(shù)如下所示：要求繪制系統(tǒng)的閉環(huán)根軌跡，分析其穩(wěn)定性，并繪制出當(dāng)k=55和k=56時(shí)系統(tǒng)的閉環(huán)沖激響應(yīng)。k=55.6381控制系統(tǒng)根軌跡分析控制系統(tǒng)根軌跡分析系統(tǒng)開(kāi)環(huán)傳遞函數(shù)為G(s)=k/s(s+1)(s+2)試尋找一個(gè)合適的k值使得閉環(huán)系統(tǒng)具有較理想的階躍響應(yīng)。例exp5_23.m控制系統(tǒng)根軌跡分析例exp5_24.m某控制系統(tǒng)的開(kāi)環(huán)傳遞函數(shù)，G(s)=k(s+1)/s^2(s+2)(s+4)要求分別繪制正反饋系統(tǒng)和負(fù)反饋系統(tǒng)的根軌跡指出它們的穩(wěn)定性情況有什么不同控制系統(tǒng)根軌跡分析控制系統(tǒng)的分析是進(jìn)行控制系統(tǒng)設(shè)計(jì)的基礎(chǔ)，同時(shí)也是工程實(shí)際當(dāng)中解決問(wèn)題的主要方法，因而對(duì)控制系統(tǒng)的分析在控制系統(tǒng)仿真中具有舉足輕重的作用。通過(guò)求取系統(tǒng)的零極點(diǎn)增益模型直接獲得系統(tǒng)的零極點(diǎn)，從而可以直接對(duì)控制系統(tǒng)的穩(wěn)定性及是否為最小相位系統(tǒng)作出判斷。控制系統(tǒng)的經(jīng)典分析方法（時(shí)域、頻域分析）是目前控制系統(tǒng)界進(jìn)行科學(xué)研究的主要方法，是進(jìn)行控制系統(tǒng)設(shè)計(jì)的基礎(chǔ)，要求熟練掌握單位階躍響應(yīng)、波特圖等常用命令的使用。根軌跡分析是求解閉環(huán)特征方程根的簡(jiǎn)單的圖解方法，要求熟練掌握根軌跡的繪制。本章小結(jié)第五章控制系統(tǒng)CAD控制系統(tǒng)根軌跡分析控制系統(tǒng)穩(wěn)定性分析控制系統(tǒng)時(shí)域分析控制系統(tǒng)頻域分析控制系統(tǒng)設(shè)計(jì)方法現(xiàn)代控制理論CAD第14次控制系統(tǒng)設(shè)計(jì)方法（一）超前校正：超前校正是用超前相角對(duì)系統(tǒng)實(shí)現(xiàn)校正的，可改善系統(tǒng)的動(dòng)態(tài)性能，即加寬系統(tǒng)頻帶，提高系統(tǒng)的響應(yīng)速度。超前校正會(huì)使系統(tǒng)的相角裕量增加，從而提高系統(tǒng)的相對(duì)穩(wěn)定性，致使閉環(huán)系統(tǒng)的頻帶擴(kuò)寬。系統(tǒng)頻帶的加寬也會(huì)帶來(lái)一定的噪聲信號(hào)，這是系統(tǒng)所不希望的。

在頻域法中，采用伯德圖進(jìn)行設(shè)計(jì)是最常見(jiàn)的。其基本思想是：改變?cè)邢到y(tǒng)開(kāi)環(huán)頻率特性的形狀，使其具有希望的低頻增益、希望的增益穿越頻率和充分的穩(wěn)定裕量。為方便起見(jiàn)，常常選用如下形式的傳遞函數(shù)的校正裝置：超前校正超前校正超前校正——頻率響應(yīng)法1.根據(jù)穩(wěn)態(tài)誤差要求確定開(kāi)環(huán)增益K2.繪制待校正系統(tǒng)對(duì)數(shù)頻率特性，計(jì)算相角裕度3.根據(jù)截止頻率要求，計(jì)算超前網(wǎng)絡(luò)參數(shù)和

選，由確定

或由期望相角裕度得出確定

再由確定得4.繪制校正后對(duì)數(shù)頻率特性5.驗(yàn)算已校正系統(tǒng)的相角裕度，若不滿足要求，重選控制系統(tǒng)設(shè)計(jì)方法解：（1）根據(jù)給定的系統(tǒng)穩(wěn)態(tài)指標(biāo)，確定系統(tǒng)的開(kāi)環(huán)增益kc，使速度誤差常數(shù)為10例4-4已知被控對(duì)象的傳函為速度誤差常數(shù)為10，相角裕量為45°

系統(tǒng)的設(shè)計(jì)指標(biāo)要求如下：控制系統(tǒng)設(shè)計(jì)方法

(2)繪制未校正系統(tǒng)的伯德圖，并計(jì)算其相角裕量>>num=2000;>>den=[1302000];>>margin(num,den)控制系統(tǒng)設(shè)計(jì)方法（3)根據(jù)給定的相角裕量要求，計(jì)算所需的相角超前量（4）計(jì)算系數(shù)（5）計(jì)算校正后系統(tǒng)的剪切頻率，即未校正系統(tǒng)增益為處的頻率，由圖可知由圖可知控制系統(tǒng)設(shè)計(jì)方法（6）確定超前校正裝置（7）畫(huà)出校正后的系統(tǒng)伯德圖，校正裝置的傳遞函數(shù)為控制系統(tǒng)設(shè)計(jì)方法num=2000*[0.151];den=conv([0.081],[1302000]);margin(num,den)控制系統(tǒng)設(shè)計(jì)方法控制系統(tǒng)設(shè)計(jì)方法（二）滯后校正：

由于滯后校正裝置給系統(tǒng)加入了滯后的相角，因而將會(huì)使得系統(tǒng)的動(dòng)態(tài)穩(wěn)定性變差。滯后校正可降低系統(tǒng)穩(wěn)態(tài)誤差，并使得閉環(huán)系統(tǒng)的帶寬降低，從而使系統(tǒng)的動(dòng)態(tài)響應(yīng)速度變慢，這有利于減小外部噪聲信號(hào)對(duì)系統(tǒng)的影響。校正裝置中零點(diǎn)出現(xiàn)在極點(diǎn)之后，可降低系統(tǒng)的穩(wěn)態(tài)誤差，改善系統(tǒng)的穩(wěn)態(tài)性能，但同時(shí)會(huì)使系統(tǒng)的動(dòng)態(tài)穩(wěn)定性變差，因滯后相位使系統(tǒng)的頻帶有所減小，降低了系統(tǒng)的動(dòng)態(tài)響應(yīng)速度，然而卻增加了系統(tǒng)的抗干擾的能力。滯后校正β>1滯后校正校正后校正前-90°0°90°-180°滯后校正——頻率響應(yīng)法根據(jù)穩(wěn)態(tài)誤差要求確定開(kāi)環(huán)增益K繪制待校正系統(tǒng)對(duì)數(shù)頻率特性，確定剪切頻率和幅值裕度，相角裕度根據(jù)相角裕度要求，選擇由控制系統(tǒng)設(shè)計(jì)方法例4-5已知單位負(fù)反饋系統(tǒng)固有的傳遞函數(shù)為若要求系統(tǒng)滿足如下性能指標(biāo)：開(kāi)環(huán)放大倍數(shù)相角裕量，幅值裕量解：原有系統(tǒng)的伯德圖如圖，試設(shè)計(jì)校正裝置。系統(tǒng)不穩(wěn)定。滯后校正裝置傳遞函數(shù)查表r=50°時(shí)w=0.4877，mag=8.9526num=5;den=conv([110],[0.51]);[mag,phase,w]=bode(num,den);r=180+phase;margin(num,den)控制系統(tǒng)設(shè)計(jì)方法校正后的系統(tǒng)伯德圖控制系統(tǒng)設(shè)計(jì)方法控制系統(tǒng)設(shè)計(jì)方法滯后超前單純采用超前校正或滯后校正均只能改善系統(tǒng)暫態(tài)或穩(wěn)態(tài)一個(gè)方面的性能。若未校正系統(tǒng)不穩(wěn)定，并且對(duì)校正后系統(tǒng)的穩(wěn)態(tài)和暫態(tài)都有較高要求時(shí)，宜于采用串聯(lián)滯后-超前校正裝置。利用校正網(wǎng)絡(luò)中的超前部分改善系統(tǒng)的暫態(tài)性能，而校正網(wǎng)絡(luò)的滯后部分則可提高系統(tǒng)的穩(wěn)態(tài)精度。用頻率特性法設(shè)計(jì)串聯(lián)校正裝置除可按本節(jié)前述分析方法和步驟之外，還可以按期望特性法去設(shè)計(jì)。滯后超前[-20][+20]滯后-超前校正例4-6要求開(kāi)環(huán)放大倍數(shù)Kv=10，相角裕量≥40°，幅值裕量≥10dB解：num=10;den=conv([110],[0.51]);

margin(num,den)求得相角裕量為-28°，幅值裕量為-10.5dB，Wg=1.41rad/s，系統(tǒng)不穩(wěn)定。取期望剪切頻率wc=1.4rad/s此時(shí)相角=-180°滯后-超前校正r=40°，滯后校正裝置傳函：超前校正裝置傳函：控制系統(tǒng)設(shè)計(jì)方法例5-6（三）滯后-超前校正：控制系統(tǒng)設(shè)計(jì)方法例5-7（四）反饋校正：基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)如圖所示的“一階倒立擺控制系統(tǒng)”中，通過(guò)檢測(cè)小車(chē)位置與擺桿的擺動(dòng)角，來(lái)適當(dāng)控制驅(qū)動(dòng)電動(dòng)機(jī)拖動(dòng)力的大小，控制器由一臺(tái)工業(yè)控制計(jì)算機(jī)（IPC）完成。一階倒立擺控制系統(tǒng)控制系統(tǒng)設(shè)計(jì)方法1系統(tǒng)模型2模型驗(yàn)證3雙閉環(huán)PID控制器設(shè)計(jì)4仿真實(shí)驗(yàn)5結(jié)論基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)

一階倒立擺系統(tǒng)模型問(wèn)題的提出：獨(dú)輪自行車(chē)，機(jī)器人行走過(guò)程中的平衡控制，火箭發(fā)射中的垂直度控制，衛(wèi)星飛行中的姿態(tài)控制，海上鉆井平臺(tái)的穩(wěn)定控制，飛機(jī)的安全著陸控制。單一剛性鉸鏈，兩自由度動(dòng)力學(xué)問(wèn)題一階倒立擺系統(tǒng)模型1）擺桿繞其重心的轉(zhuǎn)動(dòng)方程2）擺桿重心的水平運(yùn)動(dòng)方程3）擺桿重心的垂直運(yùn)動(dòng)方程4）小車(chē)的水平運(yùn)動(dòng)方程一階倒立擺系統(tǒng)模型在平衡點(diǎn)附近線性化：對(duì)上述微分方程模型通過(guò)上一節(jié)的變換可以得到一階直線倒立擺系統(tǒng)的微分方程，傳遞函數(shù)和狀態(tài)方程三種形式。由得：基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)（一）對(duì)象模型一階倒立擺的精確模型為：代入具體參數(shù)后，得到模型為：1系統(tǒng)模型基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)1系統(tǒng)模型若只考慮在其工作點(diǎn)附近的細(xì)微變化，這時(shí)可以將模型線性化，得到近似模型為其等效動(dòng)態(tài)結(jié)構(gòu)圖如下圖所示：基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)1系統(tǒng)模型（二）電動(dòng)機(jī)、驅(qū)動(dòng)器及機(jī)械傳動(dòng)裝置的模型若忽略電動(dòng)機(jī)的空載轉(zhuǎn)矩和系統(tǒng)摩擦，就可認(rèn)為驅(qū)動(dòng)器和機(jī)械傳動(dòng)裝置均為純比例環(huán)節(jié)，并假設(shè)這兩個(gè)環(huán)節(jié)的增益為和.對(duì)于交流伺服電動(dòng)機(jī)，其傳遞函數(shù)可近似為由于選擇小慣性電動(dòng)機(jī)，其時(shí)間常數(shù)相對(duì)都很小，這樣就可以將電動(dòng)機(jī)模型近似等效為比例環(huán)節(jié)。綜上所述，電動(dòng)機(jī)、驅(qū)動(dòng)器及機(jī)械傳動(dòng)裝置三個(gè)環(huán)節(jié)就可以合成一個(gè)比例環(huán)節(jié)，?；陔p閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)2模型驗(yàn)證（一）SIMULINK子系統(tǒng)子系統(tǒng)通過(guò)將大的復(fù)雜的模型分割成幾個(gè)小的模型系統(tǒng)，使得整個(gè)系統(tǒng)模型更加簡(jiǎn)捷，可讀性更高。把已存在的Simulink模型中的某個(gè)部分或全部“封裝”成子系統(tǒng)的操作程序如下：1、首先使用范圍框?qū)⒁胺庋b”成子系統(tǒng)的部分選中，包括模塊和信號(hào)線。為了使范圍框圈住所需要的模塊，常常需要事先重新安排各模塊的位置（注意：這里只能用范圍框，而不能用Shift逐個(gè)選定）。2、在模塊窗口菜單選項(xiàng)中選擇[Edit>>CreateSubsystem],

Simulink將會(huì)用一個(gè)子系統(tǒng)模塊代替選中的模塊組?；陔p閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)2模型驗(yàn)證（一）SIMULINK子系統(tǒng)3、所得子系統(tǒng)模塊將有默認(rèn)的輸入和輸出端口。輸入端口和輸出端口的默認(rèn)名稱(chēng)分別為In1和Out1。調(diào)整子系統(tǒng)和模型窗口的大小使之更加美觀。若想查看子系統(tǒng)的內(nèi)容或?qū)ψ酉到y(tǒng)進(jìn)行再編輯，可以雙擊子系統(tǒng)模塊，則會(huì)出現(xiàn)一個(gè)顯示子系統(tǒng)內(nèi)容的新窗口。在窗口內(nèi)除了原始模塊外，Simulink自動(dòng)添加輸入模塊和輸出模塊，分別代表子系統(tǒng)的輸入端口和輸出端口。改變其標(biāo)簽會(huì)使子系統(tǒng)的輸入輸出端口的標(biāo)簽也隨之變化?；陔p閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)（二）仿真驗(yàn)證

1模型封裝我們采用仿真實(shí)驗(yàn)的方法在Matlab的Simulink圖形仿真環(huán)境下進(jìn)行模型驗(yàn)證實(shí)驗(yàn)。其原理如下圖所示。其中，上半部分為精確模型仿真圖，下半部分為簡(jiǎn)化模型仿真圖。2模型驗(yàn)證基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)（二）仿真驗(yàn)證

1模型封裝利用前面介紹的Simulink壓縮子系統(tǒng)功能可將原理圖更加簡(jiǎn)捷的表示為如下形式:2模型驗(yàn)證基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)（二）仿真驗(yàn)證

2實(shí)驗(yàn)設(shè)計(jì)假定使倒立擺在（）初始狀態(tài)下突加微小沖擊力作用，則依據(jù)經(jīng)驗(yàn)知：小車(chē)將向前移動(dòng)，擺桿將倒下。下面利用仿真實(shí)驗(yàn)來(lái)驗(yàn)證正確數(shù)學(xué)模型的這一“必要性質(zhì)”。

3繪制繪圖子程序

具體程序請(qǐng)參見(jiàn)課本。

4仿真實(shí)驗(yàn)2模型驗(yàn)證基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)4仿真實(shí)驗(yàn)從下圖中可見(jiàn)：在0.1N的沖擊力作用下，擺桿倒下（由零逐步增大），小車(chē)位置逐漸增加；這一結(jié)果符合前述的實(shí)驗(yàn)設(shè)計(jì)，故可以在一定程度上確認(rèn)該“一階倒立擺系統(tǒng)”的數(shù)學(xué)模型是有效的。（二）仿真驗(yàn)證

2模型驗(yàn)證基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)從一階倒立擺系統(tǒng)動(dòng)態(tài)結(jié)構(gòu)圖中不難看出，該系統(tǒng)為“自不穩(wěn)定的非最小相位系統(tǒng)”。由于“一階倒立擺系統(tǒng)位置伺服控制”的核心是“在保證擺桿不倒的條件下，使小車(chē)位置可控,因此依據(jù)“負(fù)反饋閉環(huán)控制原理”，將系統(tǒng)小車(chē)位置作為“外環(huán)”，而將擺桿擺角作為“內(nèi)環(huán)”，則擺角作為外環(huán)內(nèi)的一個(gè)擾動(dòng)，能夠得到閉環(huán)系統(tǒng)的有效抑制。綜上，設(shè)計(jì)“一階倒立擺位置伺服控制系統(tǒng)”如下圖所示,剩下的問(wèn)題就是如何確定控制器的結(jié)構(gòu)與參數(shù)。3雙閉環(huán)PID控制器設(shè)計(jì)基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)（一）內(nèi)環(huán)控制器的設(shè)計(jì)

1、控制器結(jié)構(gòu)的選擇下圖為采用反饋校正控制的系統(tǒng)內(nèi)環(huán)框圖,反饋控制器選用PD形式。

3雙閉環(huán)PID控制器設(shè)計(jì)基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)（一）內(nèi)環(huán)控制器的設(shè)計(jì)

2、控制器參數(shù)的整定首先暫定比例環(huán)節(jié)的增益,又已知。這樣我們可以求出內(nèi)環(huán)的傳遞函數(shù)為：

3雙閉環(huán)PID控制器設(shè)計(jì)基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)（一）內(nèi)環(huán)控制器的設(shè)計(jì)

2、控制器參數(shù)的整定對(duì)于這一典型的二階系統(tǒng)我們采取典型參數(shù)整定辦法，即以保證內(nèi)環(huán)系統(tǒng)具有“快速跟隨性能特性”為條件來(lái)確定反饋控制器的參數(shù)，這樣就有：系統(tǒng)內(nèi)環(huán)的閉環(huán)傳遞函數(shù)為：3雙閉環(huán)PID控制器設(shè)計(jì)基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)（一）內(nèi)環(huán)控制器的設(shè)計(jì)

2、系統(tǒng)內(nèi)環(huán)的動(dòng)態(tài)跟隨性能指標(biāo)

(1)理論分析

(2)仿真實(shí)驗(yàn)仿真曲線如下圖所示,從曲線中可以很清楚地得知，其響應(yīng)時(shí)間和超調(diào)量與理論分析值相符合。

3雙閉環(huán)PID控制器設(shè)計(jì)基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)（二）外環(huán)控制器的設(shè)計(jì)

1、系統(tǒng)外環(huán)模型的降階

(1)對(duì)內(nèi)環(huán)等效閉環(huán)傳遞函數(shù)的近似處理

(2)對(duì)象模型的近似處理

近似條件為3雙閉環(huán)PID控制器設(shè)計(jì)基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)（二）外環(huán)控制器的設(shè)計(jì)

2控制器設(shè)計(jì)下圖為系統(tǒng)外環(huán)前向通道上傳遞函數(shù)的等效過(guò)程,我們可以將外環(huán)系統(tǒng)設(shè)計(jì)成典型Ⅱ型的結(jié)構(gòu)形式.3雙閉環(huán)PID控制器設(shè)計(jì)基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)（二）外環(huán)控制器的設(shè)計(jì)

2控制器設(shè)計(jì)系統(tǒng)的閉環(huán)結(jié)構(gòu)圖如下所示,調(diào)節(jié)器仍選擇PD形式,并采用單位反饋來(lái)構(gòu)成外環(huán)反饋通道.3雙閉環(huán)PID控制器設(shè)計(jì)（二）外環(huán)控制器的設(shè)計(jì)

2控制器設(shè)計(jì)根據(jù)典型Ⅱ型系統(tǒng)設(shè)計(jì)方法,確定外環(huán)調(diào)節(jié)器的兩個(gè)參數(shù)為,這樣可得到完整的系統(tǒng)仿真結(jié)構(gòu)如下圖所示:基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)3雙閉環(huán)PID控制器設(shè)計(jì)4仿真實(shí)驗(yàn)1畫(huà)圖子程序綜合上述內(nèi)容，可得到下圖所示的Simulink仿真系統(tǒng)結(jié)構(gòu)圖。需要強(qiáng)調(diào)的是：其中的對(duì)象模型為精確模型的封裝子系統(tǒng)形式。畫(huà)圖子程序參見(jiàn)課本?；陔p閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)2仿真結(jié)果仿真實(shí)驗(yàn)結(jié)果如下圖所示,從中可見(jiàn),雙閉環(huán)PID控制方案是有效的.基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)4仿真實(shí)驗(yàn)2仿真結(jié)果我們還可以改變倒立擺系統(tǒng)的部分參數(shù)來(lái)檢驗(yàn)系統(tǒng)是否具有一定的魯棒性。例如，我們將倒立擺的擺桿質(zhì)量改為1.1kg，此時(shí)的仿真結(jié)果如下圖所示。從仿真結(jié)果可見(jiàn)：控制系統(tǒng)仍能有效的控制其保持倒擺直立并使小車(chē)移動(dòng)到指定位置?；陔p閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)4仿真實(shí)驗(yàn)2仿真結(jié)果

為了進(jìn)一步驗(yàn)證控制系統(tǒng)的魯棒性能，并便于進(jìn)行比較，我們不妨改變倒立擺的擺桿質(zhì)量和長(zhǎng)度多作幾組試驗(yàn),部分實(shí)驗(yàn)結(jié)果如下所示?？梢?jiàn)，所設(shè)計(jì)的雙閉環(huán)PID控制器在系統(tǒng)參數(shù)的一定變化范圍內(nèi)能有效的工作，保持?jǐn)[桿直立并使小車(chē)有效定位，控制系統(tǒng)具有一定的魯棒性?；陔p閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)4仿真實(shí)驗(yàn)擺桿長(zhǎng)度不變而擺桿質(zhì)量變化時(shí)系統(tǒng)仿真結(jié)果

基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)4仿真實(shí)驗(yàn)擺桿質(zhì)量不變而擺桿長(zhǎng)度變化時(shí)系統(tǒng)的仿真結(jié)果基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)4仿真實(shí)驗(yàn)5結(jié)論本節(jié)從理論上證明了所設(shè)計(jì)的“一階直線倒立擺”雙閉環(huán)PID控制方案是可行的。本節(jié)的結(jié)果在實(shí)際應(yīng)用時(shí)（實(shí)物仿真）還有如下問(wèn)題：

(1)微分控制規(guī)律易受“噪聲”干擾，具體實(shí)現(xiàn)時(shí)應(yīng)充分考慮信號(hào)的數(shù)據(jù)處理問(wèn)題。

(2)如采用“模擬式旋轉(zhuǎn)電位器”進(jìn)行擺角檢測(cè)，在實(shí)際應(yīng)用中檢測(cè)精度不佳。

(3)實(shí)際應(yīng)用中還需考慮初始狀態(tài)下的“起擺過(guò)程控制問(wèn)題”?！耙浑A直線倒立擺的控制問(wèn)題”是一個(gè)非常典型而具有明確物理意義的“運(yùn)動(dòng)控制系統(tǒng)問(wèn)題”，對(duì)其深入的分析與應(yīng)用研究，有助于提高我們的分析問(wèn)題與解決問(wèn)題的能力。基于雙閉環(huán)PID控制的一階倒立擺控制系統(tǒng)設(shè)計(jì)第五章控制系統(tǒng)CAD控制系統(tǒng)根軌跡分析控制系統(tǒng)穩(wěn)定性分析控制系統(tǒng)時(shí)域分析控制系統(tǒng)頻域分析控制系統(tǒng)設(shè)計(jì)方法現(xiàn)代控制理論CAD第15次控制工具箱——模型轉(zhuǎn)換5.ss2tf：系統(tǒng)狀態(tài)空間形式變傳遞函數(shù)形式6.ss2zp:狀態(tài)空間變零極點(diǎn)增益形式7.tf2ss:傳函形式變狀態(tài)空間形式8.tf2zp：傳函形式變零極點(diǎn)增益形式9.zp2ss:零極點(diǎn)增益形式為狀態(tài)空間形式10.zp2tf:零極點(diǎn)增益形式變傳函形式iu為指定變換所用的輸入量控制工具箱——系統(tǒng)特性1.ctrb,obsv：可控性和可觀性矩陣可求出狀態(tài)空間系統(tǒng)的可控性和可觀性矩陣可控性矩陣：可觀性矩陣：控制工具箱——模型實(shí)現(xiàn)2.ctrbf,obsvf:可控性和可觀性階梯形式3.ss2ss：相似變換ctrbf可將系統(tǒng)分解為可控/不可控兩部分,obsvf可將系統(tǒng)分為可觀/不可觀兩部分tol可指定誤差容限完成相似變換控制工具箱——方程求解1.care:代數(shù)Riccati（李卡提）方程求解2.Lyapunov（李亞普諾夫）方程求解Riccati（李卡提）方程Lyapunov（李亞普諾夫）方程特殊形式Lyapunov（李亞普諾夫）方程lyap2函數(shù)類(lèi)似于lyap，但是運(yùn)算速度大大加快?，F(xiàn)代控制理論CAD極點(diǎn)配置例：系統(tǒng)狀態(tài)方程通過(guò)狀態(tài)反饋控制u=-Kx，配置閉環(huán)極點(diǎn)在p1=-2+j4,p2=-2-j4,p3=-10現(xiàn)代控制理論CAD%反饋控制配置閉環(huán)極點(diǎn)%開(kāi)環(huán)系統(tǒng)狀態(tài)方程描述：A=[010;001;-1-5-6];B=[0;0;1];%期望閉環(huán)極點(diǎn)J=[-2+j*4-2-j*4-10];%求反饋增益矩陣KK=acker(A,B,J)現(xiàn)代控制理論CADACKERPoleplacementgainselectionusingAckermann'sformula.

K=ACKER(A,B,P)calculatesthefeedbackgainmatrixKsuchthatthesingleinputsystem.x=Ax+Bu

withafeedbacklawofu=-KxhasclosedlooppolesatthevaluesspecifiedinvectorP,i.e.,P=eig(A-B*K).

Note:ThisalgorithmusesAckermann'sformula.ThismethodisNOTnumericallyreliableandstartstobreakdownrapidlyforproblemsofordergreaterthan10,orforweaklycontrollablesystems.Awarningmessageisprintedifthenonzeroclosed-looppolesaregreaterthan10%fromthedesiredlocationsspecifiedinP.現(xiàn)代控制理論CAD%反饋控制配置閉環(huán)極點(diǎn)%開(kāi)環(huán)系統(tǒng)狀態(tài)方程描述：A=[010;001;-1-5-6];B=[0;0;1];%期望閉環(huán)極點(diǎn)J=[-2+j*4-2-j*4-10];%求反饋增益矩陣KK=place(A,B,J)現(xiàn)代控制理論CADPLACEPoleplacementtechnique

K=PLACE(A,B,P)computesastate-feedbackmatrixKsuchthattheeigenvaluesofA-B*KarethosespecifiedinvectorP.Noeigenvalueshouldhaveamultiplicitygreaterthanthenumberofinputs.

[K,PREC,MESSAGE]=PLACE(A,B,P)returnsPREC,anestimateofhowcloselytheeigenvaluesofA-B*Kmatchthespecifiedloca