自動(dòng)控制原理實(shí)驗(yàn)七基于MATLAB控制系統(tǒng)頻域法串聯(lián)校正設(shè)計(jì)

1、實(shí)驗(yàn)七 基于MATLAB控制系統(tǒng)頻域法串聯(lián)校正設(shè)計(jì)一、實(shí)驗(yàn)?zāi)康模?）對(duì)給定系統(tǒng)設(shè)計(jì)滿足頻域或時(shí)域指標(biāo)的串聯(lián)校正裝置；（2）掌握頻域法設(shè)計(jì)串聯(lián)校正的方法；（3）掌握串聯(lián)校正環(huán)節(jié)對(duì)系統(tǒng)穩(wěn)定性及過(guò)渡過(guò)程的影響。二、實(shí)驗(yàn)原理及內(nèi)容利用MATLAB可以方便的畫(huà)出Bode圖并求出幅值裕量和相角裕量。將MATLAB應(yīng)用到經(jīng)典理論的校正方法中，可以方便的校驗(yàn)系統(tǒng)校正前后的性能指標(biāo)。通過(guò)反復(fù)試探不同校正參數(shù)對(duì)應(yīng)的不同性能指標(biāo)，能夠設(shè)計(jì)出最佳的校正裝置。1、串聯(lián)超前校正用頻域法對(duì)系統(tǒng)進(jìn)行超前校正的基本原理，是利用超前校正網(wǎng)絡(luò)的相位超前特性來(lái)增大系統(tǒng)的相位裕量，以達(dá)到改善系統(tǒng)瞬態(tài)響應(yīng)的目標(biāo)。為此，要求校正網(wǎng)絡(luò)最大的

2、相位超前角出現(xiàn)在系統(tǒng)的截止頻率（剪切頻率）處。串聯(lián)超前校正的特點(diǎn)：主要對(duì)未校正系統(tǒng)中頻段進(jìn)行校正，使校正后中頻段幅值的斜率為-20dB/dec，且有足夠大的相位裕度；超前校正會(huì)使系統(tǒng)瞬態(tài)響應(yīng)的速度變快，校正后系統(tǒng)的截止頻率增大。這表明校正后，系統(tǒng)的頻帶變寬，瞬態(tài)響應(yīng)速度變快，相當(dāng)于微分效應(yīng)；但系統(tǒng)抗高頻噪聲的能力變差。用頻率法對(duì)系統(tǒng)進(jìn)行串聯(lián)超前校正的一般步驟為：1）根據(jù)穩(wěn)態(tài)誤差的要求，確定開(kāi)環(huán)增益K。2）根據(jù)所確定的開(kāi)環(huán)增益K，畫(huà)出未校正系統(tǒng)的波特圖，計(jì)算未校正系統(tǒng)的相位裕度。3）計(jì)算超前網(wǎng)絡(luò)參數(shù)a和T。 4）確定校正網(wǎng)絡(luò)的轉(zhuǎn)折頻率。5）畫(huà)出校正后系統(tǒng)的波特圖，驗(yàn)證已校正系統(tǒng)的相位裕度。【7-

3、1】給定系統(tǒng)如圖7-1所示，試設(shè)計(jì)一個(gè)串聯(lián)校正裝置，使系統(tǒng)滿足幅值裕量大于10分貝，相位裕量45o為了滿足上述要求，試探地采用超前校正裝置Gc(s)，使系統(tǒng)變?yōu)閳D7-2的結(jié)構(gòu)。 圖7-1 校正前系統(tǒng) 圖7-2 校正后系統(tǒng) 用下面地MATLAB語(yǔ)句得出原系統(tǒng)的幅值裕量與相位裕量。 G=tf(100, 0.04, 1, 0);Gw, Pw, Wcg, Wcp=margin(G); Gw =InfPw =28.0243Wcg=InfWcp=46.9701可以看出，這個(gè)系統(tǒng)有無(wú)窮大的幅值裕量，并且其相位裕量28o，幅值穿越頻率Wcp=47rad/sec。引入一個(gè)串聯(lián)超前校正裝置：0.05s+1(0.0

4、2s+1) (自主修正) 通過(guò)下面的MATLAB語(yǔ)句得出校正前后系統(tǒng)的Bode圖如圖7-3,校正前后系統(tǒng)的階躍響應(yīng)圖如圖7-4。其中、ts1分別為校正前系統(tǒng)的幅值穿越頻率、相角裕量、調(diào)節(jié)時(shí)間，2、ts2分別為校正后系統(tǒng)的幅值穿越頻率、相角裕量、調(diào)節(jié)時(shí)間。代碼：G1=tf(100,0.04,1,0); G2=tf(100*0.05,1,conv(0.04,1,0,0.02,1) bode(G1)hold bode(G2, r-)figureG1_c=feedback(G1,1)G2_c=feedback(G2,1)Gw, Pw, Wcg, Wcp=margin(G2)step(G1_c)hold

5、 step(G2_c, r-)運(yùn)行結(jié)果： Gw =Inf，Pw = 38.9425，Wcg =Inf，Wcp = 70.7105圖7-3 校正前后系統(tǒng)的Bode圖圖7-4 校正前后系統(tǒng)的階躍響應(yīng)圖可以看出，在這樣的控制器下，校正后系統(tǒng)的相位裕量由增加到48o，調(diào)節(jié)時(shí)間由0.28s減少到0.08s。系統(tǒng)的性能有了明顯的提高，滿足了設(shè)計(jì)要求?！?-2】已知單位負(fù)反饋傳遞函數(shù)Ks2(0.5s+1) (自主修正)，試設(shè)計(jì)無(wú)源串聯(lián)超前校正網(wǎng)絡(luò)的傳遞函數(shù)，使系統(tǒng)的靜態(tài)加速度誤差系數(shù)ka=20，相角裕度 g35 。(1) k=10(2) 繪制未校正系統(tǒng)的bode圖, (g=-30.5 ，不穩(wěn)定）(3) 設(shè)計(jì)

6、串聯(lián)超前校正裝置，確定校正裝置提供的相位超前量jm（84，采用兩級(jí)超前校正網(wǎng)絡(luò)實(shí)現(xiàn)，每一級(jí)為42 ）(4) 確定校正網(wǎng)絡(luò)的轉(zhuǎn)折頻率w1和w2 ，然后確定校正器的傳遞函數(shù)代碼：num =10 ;den= 0.5,1, 0, 0 ; G0 =tf (num,den) ; Gm,Pm,Wcg,Wcp=margin (G0) w=0.1:1 : 10000 ; mag ,phase = bode (G0,w) ; magdb = 20*log10(mag) phim1 = 35 , deta= 18 ; phim = (phim1-Pm+deta)/2 ; alpha =(1+sin (phim*pi

7、/ 180) / (1-sin (phim*pi/ 180) , n = find (magdb+10*log10(alpha) = 0.0001) ; wc= n(1)+0.1, w1=wc/sqrt(alpha) , w2 = sqrt (alpha)*wc ; numc=(1/alpha)*1/w1,1,denc=1/w2,1, Gc1=tf(numc,denc),Gc = Gc1*Gc1, G=(alpha)2*Gc*G0, disp (顯示單級(jí)矯正網(wǎng)絡(luò)傳遞函數(shù)，2級(jí)串聯(lián)校正網(wǎng)絡(luò)傳遞函數(shù)及T值), T=1/w2;Gc1, Gc, alpha,T, bode(G0,G) ;hold on

8、,margin(G),figure (2), sys0 = feedback (G0 , 1); step(sys0), hold on, sys= feedback (G,1) ,step(sys)Gc1 = 0.07348 s + 0.2009 - 0.07348 s + 1Gc = 0.0054 s2 + 0.02953 s + 0.04037 - 0.0054 s2 + 0.147 s + 1ans = 4.9768 0.0735圖7-5 校正前（藍(lán)色）和校正后（橙色）系統(tǒng)的Bode圖圖7-6 校正前（藍(lán)色）和校正后（橙色）系統(tǒng)的單位階躍響應(yīng)2、串聯(lián)滯后校正滯后校正網(wǎng)絡(luò)具有低通濾波器的

9、特性，會(huì)使系統(tǒng)開(kāi)環(huán)頻率特性的中頻和高頻段增益降低和截止頻率減小，從而有可能使系統(tǒng)獲得足夠大的相位裕度，它不影響頻率特性的低頻段。但是校正后系統(tǒng)的截止頻率會(huì)減小，瞬態(tài)響應(yīng)的速度要變慢 。 單位反饋?zhàn)钚∠辔幌到y(tǒng)頻率法設(shè)計(jì)串聯(lián)滯后校正網(wǎng)絡(luò)的步驟如下：1）根據(jù)穩(wěn)態(tài)性能要求，確定開(kāi)環(huán)增益K；2）利用已確定的開(kāi)環(huán)增益，畫(huà)出未校正系統(tǒng)對(duì)數(shù)頻率特性曲線，確定未校正系統(tǒng)的截止頻率、相位裕度和幅值裕度；3）選擇不同的，計(jì)算或查出不同的值，在伯特圖上繪制曲線；4）根據(jù)相位裕度要求，選擇已校正系統(tǒng)的截止頻率； 5）確定滯后網(wǎng)絡(luò)參數(shù)b和T ；6）驗(yàn)算已校正系統(tǒng)的相位裕度和幅值裕度。 【范例7-3】已知單位負(fù)反饋系統(tǒng)Ks

10、(0. 05s+1)(0.25s+1) (自主修正)，設(shè)計(jì)要求該系統(tǒng)的相角裕度滿足，幅值裕度，靜態(tài)速度誤差。求串聯(lián)滯后校正裝置的傳遞函數(shù)。代碼：num=40,den=conv (1,0,0.05,1),den =conv(den,0.25,1);G0= tf(num,den) , margin (G0),gamma0=50,delta=6,gamma = gamma0 + delta ;w=0.01:0.01:1000,mag,phase= bode (G0,w),n =find (180 + phase-gamma 45，靜態(tài)誤差系數(shù)Kv=1000s-1.代碼如下：delta=2;s=tf(

11、s),G=1000/ (s*(0.1*s+1)*(0.001*s+1);margin(G)gm,pm=margin(G),phim1=50;phim=phim1 -pm+delta;phim=phim*pi/180;alfa=(1+sin(phim)/(1-sin(phim);a=10*log10(alfa);mag,phase,w=bode(G) ;adB=20*log10 (mag) ;Wm=spline(adB,w,-a) ;t=1/ (Wm*sqrt(alfa);Gc=(1+alfa*t*s)/ (1+t*s) ; gmc, pmc=margin(G*Gc);figure;margin (G*Gc)figure;step(feedback(G,1),figure;step(feedback(G*Gc,1)命令窗口結(jié)果：G=1000 -0.01 s3 + 0.101 s2 + sGc = 0.01716 s + 1 - 0.002042 s + 1G*Gc= 17.16 s + 1000 - 2.042e-07 s4 + 0.0003062 s3 + 0.103 s2 + s運(yùn)行結(jié)果：校正前（藍(lán)色）和校正后（紅色）系統(tǒng)的Bo