基于matlab的系统校正.docVIP

基于matlab的系统校正 1 num=1000; den=conv([1,0],conv([0.1,1],[0.001,1])); G0=tf(num,den); [Gm,Pm,Wcg,Wcp]=margin(G0); w=0.1:0.1:10000; [mag,phase]=bode(G0,w); magdb=20*log10(mag); phiml=45;deta=8;phim=phiml-Pm+deta; bita=(1-sin(phim*pi/180))/(1+sin(phim*pi/180)); n=find(magdb+10*log10(1/bita)=0.0001); wc=n(1); w1=(wc/10)*sqrt(bita); w2=(wc/10)/sqrt(bita); numc=[1/w1,1]; denc=[1/w2,1]; Gc=tf(numc,denc); G=Gc*G0; [Gmc,Pmc,Wcgc,Wcpc]=margin(G); CmcdB=20*log10(Gmc); disp(校正装置传递函数和校正后系统开环传递函数),Gc,G; disp(校正后系统的频域性能指标Kg,γ,wc),[Gmc,Pmc,Wcpc]; disp(校正装置的参数T和β值：),T=1/w1;[T,bita], bode(G0,G); hold on,margin(G); 校正装置传递函数和校正后系统开环传递函数 Transfer function: 0.01741 s + 1 -------------- 0.001955 s + 1 校正后系统的频域性能指标Kg,γ,wc 校正装置的参数T和β值： ans = 0.1123 2、 num=10; den=[0.2,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); phiml=35; deta=18; phim=(phiml-Pm+deta)/2; alpha=(1+sin(phim*pi/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]; Gcl=tf(numc,denc); Gc=Gcl*Gcl; G=(alpha)^2*Gc*G0; disp(显示单级校正网络传递函数，2级串联校正网络传递函数及α，T值), T=1/w2;Gcl,Gc,[alpha,T] t=[0:0.001:4]; bode(G0,G); hold on,margin(G),figure(2); sys0=feedback(G0,1);step(sys0,t);hold on, sys=feedback(G,1);step(sys,t); 显示单级校正网络传递函数，2级串联校正网络传递函数及α，T值 Transfer function: 0.07348 s + 0.2009 ------------------ 0.07348 s + 1 Transfer function: 0.0054 s^2 + 0.02953 s + 0.04037 -------------------------------- 0.0054 s^2 + 0.147 s + 1