资源描述:
《基于matlab的系统校正》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、基于matlab的系统校正1num=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=
2、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);holdon,marg
3、in(G);校正装置传递函数和校正后系统开环传递函数Transferfunction:0.01741s+1--------------0.001955s+1校正后系统的频域性能指标Kg,γ,wc校正装置的参数T和β值:ans=0.1740.11232、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;
4、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.
5、001:4];bode(G0,G);holdon,margin(G),figure(2);sys0=feedback(G0,1);step(sys0,t);holdon,sys=feedback(G,1);step(sys,t);显示单级校正网络传递函数,2级串联校正网络传递函数及α,T值Transferfunction:0.07348s+0.2009------------------0.07348s+1Transferfunction:0.0054s^2+0.02953s+0.04037--------------------------------0.0054s^
6、2+0.147s+1