数字控制系统.docx

数字控制系统 分析、设计与实现 学院：电信学院 班级：电信（硕）81班 姓名：尹鑫 学号：08058029 日期：2011年3月22日数字控制系统一、题目：1.1 Matlab仿真Z传递函数；1.2 Matlab仿真系统性能分析。二、仿真过程：2.1 Matlab仿真Z传递函数：2.1.1 课本第五十三页例2.12：Matlab仿真程序如下： clearclcT0=1;G1=tf(0.02*2 1.4,1 -1.7 0.72,T0);zG1,pG1,kG1=zpkdata(G1,v);pp=0.9*exp(j*pi/5);pG2=pp;conj(pp);-0.8;G2=zpk(-0.2,pG2,0.5,T0);Tzpk=G2*G1;zT,pT,kT=zpkdata(Tzpk,v);%xy2p(pT);Ttf=tf(Tzpk);dt=0:T0:50;yzpk=step(Tzpk,dt);ytf=step(Ttf,dt);figurestem(dt,yzpk,:,o);gridhold onplot(dt,ytf,+);hold off仿真结果：2.1.2 课本第五十五页例2.13：Matlab仿真程序如下： clearclcT0=1;G1=tf(1 0.5,1 0.75,T0);denG2=conv(1 -0.9,conv(1 -0.8,1 0.8);G2=tf(0.15 0,denG2,T0);T=G1+G2;zT,pT,kT=zpkdata(T,v);%xy2p(zT);pG1=pole(G1);pG2=pole(G2);figure;y,k=step(T);stem(k,y,:,filled);grid；仿真结果：2.1.3 课本第五十六页例2.14：Matlab仿真程序如下： clearclcT0=1;G1=tf(0.02*2 1.4,1 -1.7 0.72,T0);pp=0.9*exp(j*pi/5);pG2=pp;conj(pp);-0.8;G2=zpk(-0.2,pG2,0.5,T0);G3=tf(1 0.5,1 0.75,T0);denG4=conv(1 -0.9,conv(1 -0.8,1 0.8);G4=tf(0.15 0,denG4,T0);T=G4*(G3+G2*G1);zG1,pG1,kG1=zpkdata(G1,v);zG2,pG2,kG2=zpkdata(G2,v);%xy2p(pG2);zG3,pG3,kG3=zpkdata(G3,v);zG4,pG4,kG4=zpkdata(G4,v);zT,pT,kT=zpkdata(T,v);disp(zeros of Tz);%xy2p(zT);disp(poles of Tz);%xy2p(pT);resp,k=step(T);figure;plot(k,resp,+);grid仿真结果：2.1.4 课本第五十九页例2.15：Matlab仿真程序如下： clearclcnumG=0.2*1 -0.85;a=0.9*exp(j*pi/4);denG=conv(1 -a,1 -conj(a);G=tf(numG,denG,1);zG,pG,kG=zpkdata(G,v);%xy2p(pG);numF=1 -0.3;denF=conv(1 -0.8,1 -0.1);F=tf(numF,denF,1);zF,pK,kF=zpkdata(F,v);T=feedback(G,F);zT,pT,kT=zpkdata(T,v);%xy2p(pT);resp,dt=step(T);figure;plot(dt(1:150),resp(1:150),-o),grid仿真结果：2.2 Matlab仿真系统性能分析：2.2.1课本第九十一页例2.26：Matlab仿真程序如下： Gp=tf(1,1 1 0);f0=10;T0=0.1;Gz=c2d(Gp,T0,zoh);G_cz=tf(5*1-0.9,1-0.7,T0);Ta=feedback(Gz*G_cz,1,-1);k=0:1:60;ya=step(Ta,k*T0);H=tf(10,1 10);HGz=c2d(H*Gp,T0,zoh);Tb=Gz*G_cz/(1+HGz*G_cz);Tb=minreal(Tb);figure;yb=step(Tb,k*T0);plot(k*T0,ya,*,k*T0,yb,o);grid on仿真结果：2.2.2课本第九十二页例2.27：Matlab仿真程序如下： T0=0.1;Gp=tf(1,1 1 0);Gz=c2d(Gp,T0,zoh);Gc=tf(5*1-0.9,1-0.7,T0);Ta=feedback(Gz*Gc,1,-1);k=0:1:60;ya=step(Ta,k*T0);%Moa,tpa,tra,tsa,essa=kstats(k*T0,ya,1);H=tf(10,1 10);HGz=c2d(H*Gp,T0,zoh);Tb=Gz*Gc/(1+HGz*Gc);Tb=minreal(Tb);yb=step(Tb,k*T0);%Mob,tpb,trb,tsb,essb=kstats(k*T0,yb,1);figure;plot(k*T0,ya,*);grid on;hold on;plot(k*T0,yb,o);仿真结果：2.2.3课本第九十三页例2.28：Matlab仿真程序如下： T0=0.1;Gp=tf(1,1 1 0);Gz=c2d(Gp,T0,zoh);Kz=tf(5*1-0.9,1-0.7,T0);Ta=feedback(