自動(dòng)控制理論實(shí)驗(yàn)報(bào)告

1、精品實(shí)驗(yàn)一 典型環(huán)節(jié)的時(shí)域響應(yīng)一、實(shí)驗(yàn)?zāi)康亩?shí)驗(yàn)設(shè)備三、實(shí)驗(yàn)原理及內(nèi)容1典型環(huán)節(jié)的方框圖及傳遞函數(shù)2典型環(huán)節(jié)的模擬電路圖及輸出響應(yīng)四、實(shí)驗(yàn)結(jié)果比例環(huán)節(jié) 取 R0 = 200K；R1 = 100K積分環(huán)節(jié) 取 R0 = 200K；C = 1uF 取 R0 = 200K；C = 2uF比例積分環(huán)節(jié) 取 R0 = R1 = 200K；C = 1uF 取 R0=R1=200K；C=2uF慣性環(huán)節(jié) 取 R0=R1=200K；C=1uF 取 R0=R1=200K；C=2uF比例微分環(huán)節(jié) 取 R0 = R2 = 100K，R3 = 10K，C = 1uF；R1 = 100K 取 R0=R2=100K，R3

2、=10K，C=1uF；R1=200K五、心得體會(huì)實(shí)驗(yàn)二 典型系統(tǒng)的時(shí)域響應(yīng)和穩(wěn)定分析一、實(shí)驗(yàn)?zāi)康亩?、?shí)驗(yàn)設(shè)備三、實(shí)驗(yàn)原理及內(nèi)容 1、典型的二階系統(tǒng)穩(wěn)定性分析（1） 結(jié)構(gòu)框圖 圖1-2是典型二階系統(tǒng)的原理方框圖，其中T0=1s，T1=0.1s，K1。（2）模擬電路圖見圖1-3。（3）理論分析 開環(huán)傳函： 其中：開環(huán)增益。 （4）實(shí)驗(yàn)內(nèi)容 先算出臨界阻尼、欠阻尼、過阻尼時(shí)電阻R的理論值，再將理論值應(yīng)用于模擬電路中，觀察二階系統(tǒng)的動(dòng)態(tài)性能及穩(wěn)定性，應(yīng)與理論值分析基本吻合。在此實(shí)驗(yàn)中T0=1s,T1=0.2s,K1=200/R =K=200/R 閉環(huán)傳函： 其中：四、實(shí)驗(yàn)步驟 “ST”端插針與“S”端

3、插針用“短路板”短接。由于每個(gè)運(yùn)放單元均設(shè)置了鎖零場效應(yīng)管，所以運(yùn)放具有鎖零功能。將開關(guān)設(shè)在“方波”檔，分別調(diào)節(jié)調(diào)幅和調(diào)頻電位器，使得“OUT”端輸出的方波幅值為1V，周期為10s左右。 二階系統(tǒng)瞬態(tài)性能指標(biāo)的測試 按圖1-3接線，將1中的方波信號(hào)接至輸入端，取R=10K。 用示波器觀察系統(tǒng)階躍響應(yīng)C(t)，測量并記錄超調(diào)量Mp，峰值時(shí)間Tp和調(diào)節(jié)時(shí)間ts，并記錄在表1-3中。 分別按R=50K;160K;200K改變系統(tǒng)開環(huán)增益，觀察響應(yīng)的階躍響應(yīng)C(t)，測量并記錄性能指標(biāo)Mp，Tp和ts，及系統(tǒng)的穩(wěn)定性。并將測量值和計(jì)算值（實(shí)驗(yàn)前必須按公式計(jì)算出）進(jìn)行比較，參數(shù)取值及響應(yīng)曲線，詳見表1-

4、3。五、實(shí)驗(yàn)結(jié)果系統(tǒng)響應(yīng)曲線 R=10 R= R=160K R=200K 六、實(shí)驗(yàn)現(xiàn)象分析七、心得體會(huì).精品實(shí)驗(yàn)三 控制系統(tǒng)的穩(wěn)定性和穩(wěn)態(tài)誤差一、實(shí)驗(yàn)?zāi)康亩?、?shí)驗(yàn)設(shè)備三、實(shí)驗(yàn)內(nèi)容1利用MATLAB描述系統(tǒng)數(shù)學(xué)模型如果系統(tǒng)的的數(shù)學(xué)模型可用如下的傳遞函數(shù)表示則在MATLAB下，傳遞函數(shù)可以方便的由其分子和分母多項(xiàng)式系數(shù)所構(gòu)成的兩個(gè)向量惟一確定出來。即num=b0,b1 , bm; den=1,a1,a2 ,an例2-1 若系統(tǒng)的傳遞函數(shù)為試?yán)肕ATLAB表示。當(dāng)傳遞函數(shù)的分子或分母由若干個(gè)多項(xiàng)式乘積表示時(shí)，它可由MATLAB 提供的多項(xiàng)式乘法運(yùn)算函數(shù)conv( )來處理，以獲得分子和分母多項(xiàng)式向

5、量，此函數(shù)的調(diào)用格式為 p=conv(p1,p2)其中，p1和p2分別為由兩個(gè)多項(xiàng)式系數(shù)構(gòu)成的向量，而p為p1和p2多項(xiàng)式的乘積多項(xiàng)式系數(shù)向量。conv( )函數(shù)的調(diào)用是允許多級(jí)嵌套的。例2-2 若系統(tǒng)的傳遞函數(shù)為試?yán)肕ATLAB求出其用分子和分母多項(xiàng)式表示的傳遞函數(shù)。2利用MATLAB分析系統(tǒng)的穩(wěn)定性在分析控制系統(tǒng)時(shí)，首先遇到的問題就是系統(tǒng)的穩(wěn)定性。判斷一個(gè)線性系統(tǒng)穩(wěn)定性的一種最有效的方法是直接求出系統(tǒng)所有的極點(diǎn)，然后根據(jù)極點(diǎn)的分布情況來確定系統(tǒng)的穩(wěn)定性。對(duì)線性系統(tǒng)來說，如果一個(gè)連續(xù)系統(tǒng)的所有極點(diǎn)都位于左半s平面，則該系統(tǒng)是穩(wěn)定的。MATLAB中根據(jù)特征多項(xiàng)式求特征根的函數(shù)為roots(

6、)，其調(diào)用格式為r=roots(p)其中，p為特征多項(xiàng)式的系數(shù)向量；r為特征多項(xiàng)式的根。另外，MATLAB中的pzmap( )函數(shù)可繪制系統(tǒng)的零極點(diǎn)圖，其調(diào)用格式為p,z=pzmap(num,den)其中，num和den分別為系統(tǒng)傳遞函數(shù)的分子和分母多項(xiàng)式的系數(shù)按降冪排列構(gòu)成的系數(shù)行向量。當(dāng)pzmap( )函數(shù)不帶輸出變量時(shí)，可在當(dāng)前圖形窗口中繪制出系統(tǒng)的零極點(diǎn)圖；當(dāng)帶有輸出變量時(shí)，也可得到零極點(diǎn)位置，如需要可通過pzmap(p,z)繪制出零極點(diǎn)圖，圖中的極點(diǎn)用“”表示，零點(diǎn)用“o”表示。例2-3 已知系統(tǒng)的傳遞函數(shù)為給出系統(tǒng)的零極點(diǎn)圖，并判定系統(tǒng)的穩(wěn)定性。 3利用MATLAB計(jì)算系統(tǒng)的穩(wěn)態(tài)誤

7、差對(duì)于圖2-2所示的反饋控制系統(tǒng)，根據(jù)誤差的輸入端定義，利用拉氏變換終值定理可得穩(wěn)態(tài)誤差ess圖2-2 反饋控制系統(tǒng)在MATLAB中，利用函數(shù)dcgain( )可求取系統(tǒng)在給定輸入下的穩(wěn)態(tài)誤差，其調(diào)用格式為ess=dcgain (nume,dene)其中，ess為系統(tǒng)的給定穩(wěn)態(tài)誤差；nume和dene分別為系統(tǒng)在給定輸入下的穩(wěn)態(tài)傳遞函數(shù)的分子和分母多項(xiàng)式的系數(shù)按降冪排列構(gòu)成的系數(shù)行向量例2-4 已知單位反饋系統(tǒng)的開環(huán)傳遞函數(shù)為試求該系統(tǒng)在單位階躍和單位速度信號(hào)作用下的穩(wěn)態(tài)誤差。解 (1) 系統(tǒng)在單位階躍和單位速度信號(hào)作用下的穩(wěn)態(tài)傳遞函數(shù)分別為四、實(shí)驗(yàn)結(jié)果例2-1例2-2例2-3例2-4課本題目

8、驗(yàn)證:例3-13單位階躍響應(yīng)曲線例3-14例3-15例3-16subplot(1,2,1);num0=44;den0=1 10 24 44;step(num0,den0,-.);hold on;num=5.78;den=1 2.4 5.78;step(num,den);legend(原系統(tǒng)1,降價(jià)后2階系統(tǒng))subplot(1,2,2);num0=44*1 7.8;den0=1 10 24 44;step(num0,den0,-.);hold on;num=45.08;den=1 2.4 5.78;step(num,den);legend(原系統(tǒng)2,降價(jià)后2階系統(tǒng))單位階躍響應(yīng)曲線五、實(shí)驗(yàn)心得

9、實(shí)驗(yàn)四 控制系統(tǒng)的根軌跡和頻域特性分析一、實(shí)驗(yàn)?zāi)康亩?、?shí)驗(yàn)設(shè)備三、實(shí)驗(yàn)內(nèi)容1基于MATLAB的控制系統(tǒng)根軌跡分析1）利用MATLAB繪制系統(tǒng)的根軌跡利用rlocus( )函數(shù)可繪制出當(dāng)根軌跡增益k由0至+變化時(shí)，閉環(huán)系統(tǒng)的特征根在s平面變化的軌跡，該函數(shù)的調(diào)用格式為r,k=rlocus(num,den) 或 r,k=rlocus(num,den,k)其中，返回值r為系統(tǒng)的閉環(huán)極點(diǎn)，k為相應(yīng)的增益。rlocus( )函數(shù)既適用于連續(xù)系統(tǒng)，也適用于離散系統(tǒng)。rlocus(num,den)繪制系統(tǒng)根軌跡時(shí)，增益k是自動(dòng)選取的，rlocus(num,den, k)可利用指定的增益k來繪制系統(tǒng)的根軌跡。

10、在不帶輸出變量引用函數(shù)時(shí)，rolcus( )可在當(dāng)前圖形窗口中繪制出系統(tǒng)的根軌跡圖。當(dāng)帶有輸出變量引用函數(shù)時(shí)，可得到根軌跡的位置列向量r及相應(yīng)的增益k列向量，再利用plot(r,x)可繪制出根軌跡。2）利用MATLAB獲得系統(tǒng)的根軌跡增益在系統(tǒng)分析中，常常希望確定根軌跡上某一點(diǎn)處的增益值k，這時(shí)可利用MATLAB中的rlocfind( )函數(shù)，在使用此函數(shù)前要首先得到系統(tǒng)的根軌跡，然后再執(zhí)行如下命令k,poles=rlocfind(num,den) 或 k,poles=rlocfind(num,den,p)其中，num和den分別為系統(tǒng)開環(huán)傳遞函數(shù)的分子和分母多項(xiàng)式的系數(shù)按降冪排列構(gòu)成的系數(shù)向

11、量；poles為所求系統(tǒng)的閉環(huán)極點(diǎn)；k為相應(yīng)的根軌跡增益；p為系統(tǒng)給定的閉環(huán)極點(diǎn)。例3-1 已知某反饋系統(tǒng)的開環(huán)傳遞函數(shù)為試?yán)L制該系統(tǒng)根軌跡，并利用根軌跡分析系統(tǒng)穩(wěn)定的k值范圍。由此可見根軌跡與虛軸交點(diǎn)處的增益k=6，這說明當(dāng)k6時(shí)，系統(tǒng)不穩(wěn)定；利用rlocfind( )函數(shù)也可找出根軌跡從實(shí)軸上的分離點(diǎn)處的增益k =0.38, 這說明當(dāng)0k0.38時(shí)，系統(tǒng)為單調(diào)衰減穩(wěn)定，當(dāng)0.38ksubplot(2,1,1);semilogx(w,20*log10(mag);subplot(2,1,2);semilogx(w,phase)如果只想繪制出系統(tǒng)的Bode圖，而對(duì)獲得幅值和相位的具體數(shù)值并不感興

12、趣，則可以采用如下簡單的調(diào)用格式 bode(num,den)例3-3 已知二階系統(tǒng)的開環(huán)傳遞函數(shù)為繪制出當(dāng)wn=3和=0.3時(shí)系統(tǒng)的Bode圖。在曲線窗口中，通過利用鼠標(biāo)單擊曲線上任意一點(diǎn)，可以獲得此點(diǎn)所對(duì)應(yīng)的系統(tǒng)在該點(diǎn)的頻率與幅值或頻率與相位等有關(guān)信息。2）利用MATLAB繪制系統(tǒng)的Nyquist圖MATLAB提供的函數(shù)nyquist( )可以繪制系統(tǒng)Nyquist圖，該函數(shù)的調(diào)用格式為Re,Im,w=nyquist(num,den)其中，num和den分別為系統(tǒng)開環(huán)傳遞函數(shù)的分子和分母多項(xiàng)式的系數(shù)按降冪排列構(gòu)成的系數(shù)行向量；Re,Im和w分別為頻率特性的實(shí)部向量、虛部向量和對(duì)應(yīng)的頻率向量。有了這些值就可利用命令plot(Re,Im)來直接繪出系統(tǒng)的奈奎斯特圖。當(dāng)然，Nyquist圖也可采用與Bode圖類似的簡單命令來直接繪制。例3-4 已知系統(tǒng)的開環(huán)傳遞函數(shù)為繪制Nyquist圖，并判斷系統(tǒng)的穩(wěn)定性。在Nyquist曲線窗口中，也可利用鼠標(biāo)通過單擊曲線上任意一點(diǎn)，獲得此點(diǎn)所對(duì)應(yīng)的系統(tǒng)的開環(huán)頻率特性，在該點(diǎn)的實(shí)部和虛部及其頻率的值。四、實(shí)驗(yàn)結(jié)果例3-1根軌跡圖例3-2根軌跡圖例3-3wn=3;zeta=0.3; w=lo