CAD实验五非线性系统仿真.doc

《CAD实验五非线性系统仿真.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验五非线性系统仿真【实验目的】1，了解三种典型非线性环节对系统性能的影响2，了解三种典型非线性环节的特性【实验内容】如图所示系统，分别针对以下三种情况加入非线性环节，采用离散相似法进行仿真在A处加入饱和非线性环节在A处加入死区非线性环节在A处加入滞环非线性环节要求：适当改变S同时与非线性环节的情况比较,分析各非线性环节对系统动态性能的影响。程序如下：%------filename:satu.m------饱和非线性子函数functionuc=satu(ur,s1)if(abs(ur)>=s1)if(ur>0)uc=s

2、1;elseuc=-s1;end;elseuc=ur;end;%------filename:dead.m------死区非线性子函数functionuc=dead(ur,s1)if(abs(ur)>=s1)if(ur>0)uc=ur-s1;elseuc=ur+s1;end;elseuc=0;end;%------filename:backlash.m------滞环非线性子函数function[Uc,Ubb]=backlash(Urb,Ur,Ucb,S1)if(Ur>Urb)if((Ur-S1)>=Ucb)Uc=Ur

3、-S1;elseUc=Ucb;end;elseif(Ur

4、:,3);D=P(:,4);m=length(WIJ(:,1));W0=zeros(n,1);W=zeros(n,n);fork=1:mif(WIJ(k,2)==0);W0(WIJ(k,1))=WIJ(k,3);elseW(WIJ(k,1),WIJ(k,2))=WIJ(k,3);end;end;fori=1:n%按环节离散化if(A(i)==0);%积分或积分环节各系数FI(i)=1;FIM(i)=h*C(i)/B(i);FIJ(i)=h*h*C(i)/B(i)/2;FIC(i)=1;FID(i)=0;if(D(i)~

5、=0);FID(i)=D(i)/B(i);endelseFI(i)=exp(-h*A(i)/B(i));%惯性环节个比例系数FIM(i)=(1-FI(i))*(C(i)/A(i));FIJ(i)=h*C(i)/A(i)-FIM(i)*B(i)/A(i);FIC(i)=1;FID(i)=0;if(D(i)~=0);FIC(i)=C(i)/D(i)-A(i)/B(i);FID(i)=D(i)/B(i);endendendZ=[0000];S=[0000];Y=zeros(n,1);X=Y;Uk=zeros(n,1);Ubb

6、=Uk;y=0;t=T0:h:Tf;N=length(t);fork=1:N-1Ub=Uk;Uk=W*Y+W0*Y0;fori=1:nif(Z(i)~=0)if(Z(i)==1)Uk(i)=satu(Uk(i),S(i));endif(Z(i)==2)Uk(i)=dead(Uk(i),S(i));endif(Z(i)==3)[Uk(i),Ubb(i)]=backlash(Ubb(i),Uk(i),Ub(i),S(i));endendendUdot=(Uk-Ub)/h;Uf=2*Uk-Ub;X=FI'.*X+FIM'.*

7、Uk+FIJ'.*Udot;%各环节状态Yb=Y;Y=FIC'.*X+FID'.*Uf;fori=1:nif(Z(i)~=0)if(Z(i)==4)Y(i)=satu(Y(i),S(i));endif(Z(i)==5)Y(i)=dead(Y(i),S(i));endif(Z(i)==6)[Y(i),Ubb(i)]=backlash(Ubb(i),Y(i),Yb(i),S(i));endendendy=[y,Y(nout)];endplot(t,y,'y')holdon仿真结果如下：S=1饱和非线性时，不同的s值S=0与

8、s=5比较接近；s值越大越不容易进入非线性S=5S=3死区非线性时，不同的s值：S=0S=5S=3S=1滞环非线性时，不同的s值S=0S=3S=5S=1结果分析：饱和非线性对系统动态性能的影响：(1)饱和特性使系统开环增益下降，对动态响应的平稳性有利。(2)如果饱和点过低，则在提高系统平稳性的同时，将使系统的快速性和稳态跟踪精度有