实验8线性系统串联校正

实验八线性系统串联校正一、实验目的熟练掌握用MATLAB语句绘制频域曲线。掌握控制系统频域范围内的分析校正方法。掌握用频率特性法进行串联校正设计的思路和步骤。二、基础知识控制系统设计的思路之一就是在原系统特性的基础上，对原特性加以校正，使之达到要求的性能指标。最常用的经典校正方法有根轨迹法和频域法。而常用的串联校正装置有超前校正、滞后校正和超前滞后校正装置。本实验主要讨论在MATLAB^境下进行串联校正设计。1.基于频率法的串联超前校正超前校正装置的主要作用是通过其相位超前效应来改变频率响应曲线的形状，产生足够大的相位超前角，以补偿原来系统中元件造成的过大的相位滞后。因此校正时应使校正装置的最大超前相位角出现在校正后系统的开环截止频率c处。例5-1:单位反馈系统的开环传递函数为G(s) ，试确定串联校正装S(s1)置的特性，使系统满足在斜坡函数作用下系统的稳态误差小于0.1，相角裕度r450。essessLimsE(s)Lims——sss0.1K10ks(s1)取K12，求原系统的相角裕度。>>num0=12;den0=[2,1,0];w=0.1:1000;[gm1,pm1,wcg1,wcp1]=margin(num0,den0);[mag1,phase1]=bode(numO,den0,w);Bode图grid;%Bode图grid;%十算系统的相角裕度和幅值裕度，并绘制出[gm1,pm1,wcg1,wcp1]margin(numO,denO)图图5-1原系统的Bode图ans=Inf11.6548Inf2.4240由结果可知，原系统相角裕度 r11.6o,c2.4rad/s,不满足指标要求，系统的Bode图如图5-1所示。考虑采用串联超前校正装置，以增加系统的相角裕度。确定串联装置所需要增加的超前相位角及求得的校正装置参数。(45o,0为原系统的相角裕度，取5o,令mJ>1rfi'IrtnmfrRr-11'riAglit:sinsine=5;r=45; r0=pm1;phic=(r-r0+e)*pi/180;alpha=(1+sin(phic))/(1-sin(phic));将校正装置的最大超前角处的频率 m作为校正后系统的剪切频率 c。则有:120lgq(jJG°(jc)0|G°(jJ 即原系统幅频特性幅值等于20lg.时的频率，选为根据m=c，求出校正装置的参数T即T=min(abs(mag1-1/sqrt(alpha)));wc=w(II);T=1/(wc*sqrt(alpha));numc=[alpha*T,1];denc=[T,1];[num,den]=series(numO,denO,numc,denc);%原系统与校正装置串联[gm,pm,wcg,wcp]=margin(num,den);%返回系统新的相角裕度和幅值裕度printsys(numc,denc) %显示校正装置的传递函数disp('校正之后的系统开环传递函数为 ：'printsys(num,den) %显示系统新的传递函数[mag2,phase2]=bode(numCjdenc,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');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),；';

'校正后：幅值裕量=；num2str(20*log10(gm)),'db；'相位裕量=',num2str(pm),'0']);nunydsn=0.503SI2E+1—G叫-(5—G叫-(5匚Gc2基于频率法的串联滞后校正滞后校正装置将给系统带来滞后相角。引入滞后装置的真正目的不是为了提供一个滞后相角，而是要使系统增益适当衰减，以便提高系统的稳态精度。K,K,试确定30/s,相角裕度30K30例5-2:单位反馈系统的开环传递函数为G(s)s(0.1s串联校正装置的特性，使校正后系统的静态速度误差系数等于r40°，幅值裕量不小于10dB，截止频率不小于2.3rad/s解：根据系统静态精度的要求，选择开环增益KLimsG(s)Limsvs0s0s(0.1s1)(0.2s1)利用MATLAB绘制原系统的bode图和相应的稳定裕度>>num0=30;denO=conv([1,0],conv([O・1,1],[0.2,1]));w=logspace(-1,1.2);[gm1,pm1,wcg1,wcp1]=margin(num0,den图5-3原系统的Bode图

0);0);[mag1,phase1]=bode(num0,den0,w);[gm1,pm1,wcg1,wcp1]margin(num0,den0)grid;ans=0.5000-17.23907.07119.7714由结果可知，原系统不稳定，且截止频率远大于要求值。系统的Bode图如图5-3所示，考虑采用串联超前校正无法满足要求，故选用滞后校正装置。根据对相位裕量的要求，选择相角为 1800 (50~100, 400)处的频率作为校正后系统的截止频率0dB处的频率作为校正后系统的截止频率0dB时所需的衰减量为20lgC。确定原系统在新 c处的幅值衰减到一般取校正装置的转折频率分别为e=10;r=40;r0=pm1;e=10;r=40;r0=pm1;phi=(-180+r+e);=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(numc,denc) %显示校正装置的传递函数disp('校正之后的系统开环传递函数为 ：'；printsys(num,den) %显示系统新的传递函数

[mag2,phase2]=bode(numcQenc,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');subplot(2,1,2);semilogx(w,phase,w,phase'-',w,phase2'；w,(w-180-w),)grid; ylabel()相位(0)'xlabel()频率(rad/sec)');title(['校正前：幅值裕量=',num2str(20*log10(gm1)),'db)‘相位裕量=',num2str(pm1),')nun/den=4.0566s+1■L<-K'校正后：幅值裕量=nun/den=4.0566s+1■L<-K42.3922s+1校正Z后的系统开环传递函数为:FiTii-i-i'v-.Inun/den.=FiTii-i-i'v-.I121.6983s+30图5-4系统校正前后的传递函数及 Bode图0A5934s4+12*9177s3+4展・2&22图5-4系统校正前后的传递函数及 Bode图=；num2str(pm),'));3.基于频率法的串联滞后-超前校正滞后-超前校正装置综合了超前校正和滞后校正的优点，从而改善了系统的性能例5-3:单位反馈系统的开环传递函数为例5-3:单位反馈系统的开环传递函数为G(s)s(s1)(0.4s1)，若要求相角裕度r45o，幅值裕量大于10dB,心10(1/s),试确定串联校正装置的特性。解：根据系统静态精度的要求，选择开环增益KvLimsG(s)K10s0利用MATLAB绘制原系统的bode图和相应的稳定裕度，如图5-5所示>>num0=10;den0=conv([1,0],conv([1,1],[0.4,1]));w=logspace(-1,1.2);I■IL：八，.“■ I■IL：八，.“■ i A1 .a图5-5原系统的Bode图grid;[mag1,phase1]=bode(num0,den0,w);[gm1,pm1,wcg1,wcp1]margin(num0,den0)ans=0.3500-24.19181.58112.5520由结果可以看出，单级超前装置难以满足要求，故设计一个串联滞后 -超前装置。选择原系统180。的频率为新的截止频率 c，则可以确定滞后部分的T2和。其中1 T2 10。由原系统，c1.58rad/s，此时的幅值为T10 2 0：T c29.12dB。根据校正后系统在新的幅值交接频率处的幅值必须为 0dB，确定超前校正部分的Tj。在原系统(c,20lgG°(jJ)，即(1・58,-9.12)处画一条斜率为20dB/dec的直线，此直线与OdB线及-20dB线的交点分别为超前校正部分的两个转折频率。wc=1.58;beit=10;T2=10/wc;lw=20*log10(w/1.58)-9.12;=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);semilogx(w,phase,w,phase)-',w,phase2)；w,(w-180-w),':);

grid; ylabel('相位(0));xlabel('频率(rad/sec)');title(['校正后：幅值裕量=；num2str(20*log10(gm)),'db','相位裕量=',num2str(pm),'0']);num/den=134237八2+84501s+1B.4237sA2+615032校正N后的慕统开环传谨函数为134.2374八2+84.5006s+104U-■GOjSCjD^Oc S~~!!!53和4U-■GOjSCjD^Oc S~~!!!53和5尹,*441345严4十W2.72S3sA3+64.P032洱十?粧率:radome)Bode图图5-6系统校正前后的传递函数及、实验内容1某单位负反馈控制系统的开环传递函数为 G(s)s(s1)，试设计超前校正装置，使校正后系统的静态速度误差系数Kv20s1，相位裕量 500，增益裕量20lgKg10dB。Matlab程序：num0=4;den0=[1,1,0];w=0.1:1000;[gm1,pm1,wcg1,wcp1]=margin(num0,den0);[mag1,phase1]=bode(num0,den0,w);[gm1,pm1,wcg1,wcp1]margin(num0,den0)ans=Inf28.0202 Inf 1.8791由结果可知，原系统相角裕度 为28・02°wc=1・879rad/s,不满足指标要求，系统的Bode图如下所示。考虑采用串联超前校正装置，以增加系统的相角裕度。BodeDiagramGm=InfdB(atInfrad/sec),Pm=28deg(at1.88rad/sec)100校正后的传递函数及Bode图：Matlab：e=5.5275;r=50; r0=pm1;phic=(r-r0+e)*pi/180;alpha=(1+sin(phic))/(1-sin(phic));=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);[gm,pm,wcg,wcp]=margin(num,den);%返回系统新的相角裕度和幅值裕度printsys(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');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),'O']);结果：num/den=2.7973s+10.0060511s+1校正之后的系统开环传递函数为num/den=67.1354s+240.0060511sA4+1.0182sA3+3.0182sA2+3.0061s+1校正后的传递函数及Bode图：--Go,-Gc,GoGc校正前：幅值裕量=lnfdb相位裕量=28.02020校正后:幅值裕量=lnfdb相位裕量=50.00010

2某单位负反馈控制系统的开环传递函数为 G(s)E,试设计 个合适的滞后校正网络，使系统阶跃响应的稳态误差约为 0.04,相角裕量约为45°由稳态误差可求出，1/(k+1)=0.04得出k=24Matlab程序：num0=24;den0=conv([1,1],conv([1,1],[1,1]));w=logspace(-1,1.2);[gm1,pm1,wcg1,wcp1]=margin(num0,den0);[mag1,phase1]=bode(num0,den0,w);[gm1,pm1,wcg1,wcp1]margin(num0,den0)grid;结果：ans=0.3334Bode-29.14671.73222.7056图--Go,-Gc,GoGcDiagramGm=-954--Go,-Gc,GoGcDiagramGm=-954dB(at173血msec).Pm=-291血g間271Frequency{r#c^WG|-校正Matlab程序:e=10;r=45;r0=pm1;phi=(-180+r+e);=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(numc,denc) %显示校正装置的传递函数disp(校正之后的系统开环传递函数为:')；printsys(num,den) %显示系统新的传递函数

nag2細a=22j=bod0=num^nconv%计算指指鞭率内校正装置的相角范围和幅值范围[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');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)),'dbT相位裕量=',num2str(pm1),'0';校正后：幅值裕量=',num2str(20*log10(gm)),'dbT相位裕量=',num2str(pm),'0']);结果：num/den=11.4062s+116.386s+1校正之后的系统开环传递函数为：num/den=273.75s+24116.386sA4+350.1579sA3+352.1579sA2+119.386s+1图如下:•Go,-Gc,GoGc道蝎$-10(_2p-10(相-20C-30CIU■_2IU频率rad/sec)3某单位负反馈控制系统的开环传递函数为3某单位负反馈控制系统的开环传递函数为G(s)KF'试设计一滞后-超前校正装置，使校正后系统的静态速度误差系数Kv10s1,相位裕量500，增益裕量20lgKg10dB由Kv10s1可知k=20Matlab程序：

[gm1,pm1,wcg1,wcp1]=margin(num0,den0);[mag1,phase1]=bode(num0,den0,w);[gm1,pm1,wcg1,wcp1]margin(num0,den0)grid;结果：ans=0.3000-28.08141.41422.4253Bode图：BodeDiagram-10.5dB(at1.41rad/sec),Pm=-28.1deg(at2.43rad/sec)-10101010Frequency(rad/sec)210选择原系统BodeDiagram-10.5dB(at1.41rad/sec),Pm=-28.1deg(at2.43rad/sec)-10101010Frequency(rad/sec)210选择原系统1800的频率为新的截止频率c,则可以确定滞后部分的T2和。其中10。由原系统，c1.58rad/s,10。由原系统，c1.58rad/s,此时的幅值为T2 10c20.19.12dB。根据校正后系统在新的幅值交接频率处的幅值必须为 0dB，确定超前校正部分的「。在原系统(c,20lgG°(jJ)，即(1・58,-9.12)处画一条斜率为20dB/dec的直线，此直线与OdB线及-20dB线的交点分别为超前校正部分的两个转折频率。num0=20;den0=conv([1,0],conv([1,1],[1,2]));w=logspace(-1,1.2);[gm1,pm1,wcg1,wcp1]=margin(num0,den0);[mag1,phase1]=bode(num0,den0,w);[gm1,pm1,wcg1,wcp1]margin(num0,den0)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/