《线性系统串联校正》由会员分享，可在线阅读，更多相关《线性系统串联校正（14页珍藏版）》请在金锄头文库上搜索。
1、实验八 线性系统串联校正一、实验目的1熟练掌握用 MATLAB 语句绘制频域曲线。2掌握控制系统频域范围内的分析校正方法。3掌握用频率特性法进行串联校正设计的思路和步骤。二、基础知识控制系统设计的思路之一就是在原系统特性的基础上，对原特性加以校正，使之达到要求的性能指标。最常用的经典校正方法有根轨迹法和频域法。而常用的串联校正装置有超前校正、滞后校正和超前滞后校正装置。本实验主要讨论在 MATLAB 环境下进行串联校正设计。1基于频率法的串联超前校正超前校正装置的主要作用是通过其相位超前效应来改变频率响应曲线的形状，产生足够大的相位超前角，以补偿原来系统中元件造成的过大的相位滞后。因此校正时应
2、使校正装置的最大超前相位角出现在校正后系统的开环截止频率处。c例 5-1：单位反馈系统的开环传递函数为 ，试确定串联校正)1()sKG装置的特性，使系统满足在斜坡函数作用下系统的稳态误差小于 0.1，相角裕度。045r解：根据系统静态精度的要求，选择开环增益 10.)1()(200 KskLimsEies取 ，求原系统的相角裕度。12Knum0=12; den0=2,1,0; w=0.1:1000;gm1,pm1,wcg1,wcp1=margin(num0,den0);mag1,phase1=bode(num0,den0,w);gm1,pm1,wcg1,wcp1margin(num0,den0
3、) %计算系统的相角裕度和幅值裕度，并绘制出Bode 图grid;ans =Inf 11.6548 Inf 2.4240由结果可知，原系统相角裕度 ， ，不满足指标要求，系06.1rsradc/4.2统的 Bode 图如图 5-1 所示。考虑采用串联超前校正装置，以增加系统的相角裕度。确定串联装置所需要增加的超前相位角及求得的校正装置参数。 ),5,45( 000 cmc 令取为 原 系 统 的 相 角 裕 度 msin1e=5; r=45; r0=pm1;phic=(r-r0+e)*pi/180;alpha=(1+sin(phic)/(1-sin(phic);将校正装置的最大超前角处的频率
4、作为m校正后系统的剪切频率 。则有：c 1)(0)(lg200 ccc jGjjG即原系统幅频特性幅值等于 时的频率，选为 。lg2c根据 = ，求出校正装置的参数 。即 。mc Tc1il,ii=min(abs(mag1-1/sqrt(alpha);wc=w( ii); T=1/(wc*sqrt(alpha);numc=alpha*T,1; denc=T,1;num,den=series(num0,den0,numc,denc); %原系统与校正装置串联图 5-1 原系统的 Bode 图 gm,pm,wcg,wcp=margin(num,den); %返回系统新的相角裕度和幅值裕度print
5、sys(numc,denc) %显示校正装置的传递函数disp(校正之后的系统开环传递函数为: );printsys(num,den) %显示系统新的传递函数mag2,phase2=bode(numc,denc,w); %计算指定频率内校正装置的相角范围和幅值范围mag,phase=bode(num,den,w); %计算指定频率内系统新的相角范围和幅值范围 subplot(2,1,1);semilogx(w,20*log10(mag),w,20*log10(mag1),-,w,20*log10(mag2),-.);grid; ylabel(幅值(db); title(-Go,-Gc,GoGc
6、);subplot(2,1,2); semilogx(w,phase,w,phase1,-,w,phase2,-,w,(w-180-w),:);grid; ylabel(相位( 0); xlabel(频率(rad/sec);title(校正前：幅值裕量=,num2str(20*log10(gm1),db,相位裕量=,num2str(pm1),0;校正后：幅值裕量=,num2str(20*log10(gm),db,相位裕量=,num2str(pm),0);图 5-2 系统校正前后的传递函数及 Bode 图2基于频率法的串联滞后校正滞后校正装置将给系统带来滞后相角。引入滞后装置的真正目的不是为了提
7、供一个滞后相角，而是要使系统增益适当衰减，以便提高系统的稳态精度。滞后校正的设计主要是利用它的高频衰减作用，降低系统的截止频率，以便能使得系统获得充分的相位裕量。例 5-2：单位反馈系统的开环传递函数为 ，试确)12.0)(1.()ssKG定串联校正装置的特性，使校正后系统的静态速度误差系数等于 30/s，相角裕度 ，幅值裕量不小于 10dB，截止频率不小于 2.3rad/s。04r解：根据系统静态精度的要求，选择开环增益 30)12.0)(1.()(00 KssKLimsGiKsv利用 MATLAB 绘制原系统的 bode 图和相应的稳定裕度。num0=30; den0=conv(1,0,c
8、onv(0.1,1,0.2,1); w=logspace(-1,1.2);gm1,pm1,wcg1,wcp1=margin(num0,den0);mag1,phase1=bode(num0,den0,w);gm1,pm1,wcg1,wcp1margin(num0,den0) grid;ans =0.5000 -17.2390 7.0711 9.7714由结果可知，原系统不稳定，且截止频率远大于要求值。系统的 Bode 图如图 5-3 所示，考虑采用串联超前校正无法满足要求，故选用滞后校正装置。根据对相位裕量的要求，选择相角为处的频率作为校正后系统的截止频率 。)40,15(1800 c确定原系
9、统在新 处的幅值衰减到 0dB 时所需的衰减量为 。一般取校clg20正装置的转折频率分别为 和 。cT)105(T图 5-3 原系统的 Bode 图e=10; r=40; r0=pm1;phi=(-180+r+e);il,ii=min(abs(phase1-phi);wc=w( ii); beit=mag1(ii); T=10/wc;numc= T,1; denc= beit*T,1;num,den=series(num0,den0,numc,denc); %原系统与校正装置串联gm,pm,wcg,wcp=margin(num,den); %返回系统新的相角裕度和幅值裕度printsys(n
10、umc,denc) %显示校正装置的传递函数disp(校正之后的系统开环传递函数为: );printsys(num,den) %显示系统新的传递函数mag2,phase2=bode(numc,denc,w); %计算指定频率内校正装置的相角范围和幅值范围mag,phase=bode(num,den,w); %计算指定频率内系统新的 相角范围和幅值范围 subplot(2,1,1);semilogx(w,20*log10(mag),w,20*log10(mag1),-,w,20*log10(mag2),-.);grid; ylabel(幅值(db); title(-Go,-Gc,GoGc);su
11、bplot(2,1,2); semilogx(w,phase,w,phase1,-,w,phase2,-,w,(w-180-w),:);grid; ylabel(相位( 0); xlabel(频率(rad/sec);title(校正前：幅值裕量=,num2str(20*log10(gm1),db,相位裕量=,num2str(pm1),0;校正后：幅值裕量=,num2str(20*log10(gm),db,相位裕量=,num2str(pm),图 5-4 系统校正前后的传递函数及 Bode 图 0);3基于频率法的串联滞后-超前校正滞后-超前校正装置综合了超前校正和滞后校正的优点，从而改善了系统的
12、性能。例 5-3：单位反馈系统的开环传递函数为 ，若要求)14.0)()(ssKG相角裕度 ，幅值裕量大于 10dB， ，试确定串联校正装置的045r /10Kv特性。解：根据系统静态精度的要求，选择开环增益 10)(0sGLimsv利用 MATLAB 绘制原系统的 bode 图和相应的稳定裕度，如图 5-5 所示。num0=10; den0=conv(1,0,conv(1,1,0.4,1); w=logspace(-1,1.2);gm1,pm1,wcg1,wcp1=margin(num0,den0);mag1,phase1=bode(num0,den0,w);gm1,pm1,wcg1,wcp
13、1margin(num0,den0) grid;ans =0.3500 -24.1918 1.5811 2.5520由结果可以看出，单级超前装置难以满足要求，故设计一个串联滞后-超前装置。选择原系统 的频率为新的截止频率 ，则可以确定滞后部分的 和018c2T。其中， 。由原系统， ，此时的幅值为ccT1.01220sradc/58.19.12dB。根据校正后系统在新的幅值交接频率处的幅值必须为 0dB，确定超前校正部分的 。在原系统 ，即（1.58,-9.12）处画一条斜率为1T)(lg20,(ccjG图 5-5 原系统的 Bode 图的直线，此直线与 0dB 线及-20dB 线的交点分别为
14、超前校正部分的两decB/20个转折频率。wc=1.58; beit=10; T2=10/wc;lw=20*log10(w/1.58)-9.12; il,ii=min(abs(lw+20); w1=w(ii);numc1=1/w1,1;denc1=1/ (beit*w1),1;numc2= T2,1;denc2= beit*T2,1;numc,denc=series(numc1,denc1,numc2,denc2);num,den=series(num0,den0,numc,denc);printsys(numc,denc) disp(校正之后的系统开环传递函数为: );printsys(num,den) mag2,phase2=bode(numc,denc,w); mag,phase=bode(num,den,w); gm,pm,wcg,wcp=margin(num,den); subplot(2,1,1);semilogx(w,20*log10(mag),w,20*log10(mag1),-,w,20*log10(mag2),-.);grid; ylabel(幅值(db); title(-Go,-Gc,GoGc);subplot(2,1,2); sem
