动力工程计算机控制matlab上机作业

《动力工程计算机控制matlab上机作业》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、2——12题已知系统如图所示►►HG(z)r⑺r(.v)n+、4、0.543(1-0.5z_,)(1-0.368z_,)*仏奸此丄〜、10(l-z，(l+0.718z，5(5+1)采用零阶保持器，单位阶跃输入，采样周期1秒。(1)求出系统的闭环Z传递函数⑵、误差Z传递函数£(2)、控制命令Z传递函数/7(z)和输出的Z变换K(z)；(2)编程上机，算出y(Z:7)(々=0,1，…，20),画出响应曲线，求出超调量％、过渡时间Matlab程序如下:»Ts=l;z=0;p=l;k=l;〉〉Rz=zpk(z，p，k，Ts)%按零点、极点、比例系数生成

2、R⑵表达式(注意含采样周期Ts)»z=[];p=[0-1];k=10;〉〉Gs=zpk(z，p，k);%按零点、极点、比例系数生成G(s)表达式»Gd=c2d(Gs,Ts；zoh')%将G⑻按H0⑻离散化»z=[0.50.368];p=[l-0.718];k=0.543;〉〉Dd=zpk(z，p，k，Ts)%按零点、极点、比例系数生成D(z)表达式(注意含采样周期Ts)»G0d=Dd*Gd»Gcd=feedback(GOd，l)%按单位负反馈生成W环传递函数Gc(z)表达式»Yz=Rz*Gcd»[ykt]=step(Gcd,0:1:20)%求G

3、c⑵对阶跃输入的响应，即y(kT)yk==01.99761.00390.99781.00110.99921.00060.99961.00030.99981.00020.99991.00010.99991.00001.00001.00001.00001.00001.00001.0000t=01234567891011121314151617181920»pp=spline(t,yk);»yy=ppval(pp,linspace(0,20));»plot(linspace(0,20),yy,*b');»holdon»xlabel(’t’);ylab

4、el(’y’)；»Ez=l/(l+GOd)»[ekt]=step(Ez,0:1:20)%求E(z)对阶跃输入的响应ek=1.0000-0.9976-0.00390.0022-0.00110.0008-0.00060.0004-0.00030.0002-0.00020.0001-0.00010.0001-0.00000.0000-0.00000.0000-0.00000.0000-0.0000t=01234567891011121314151617181920»Uz=Dd/(l+G0d)»[ukt]=step(Uz,0:1:20)%求U(z)对阶

5、跃输入的响应uk=0.5430-0.85990.7153-0.51230.3671-0.26310.1886-0.13510.0968-0.06940.0497-0.03560.0255-0.01830.0131-0.00940.0067-0.00480.0035-0.00250.0018t=01234567891011121314151617181920»pp=spline(t，uk);»yy=ppval(pp，linspace(0,20));»plot(linspace(0,20)，yy，’b’);»holdon»Ez=l/(l+G0d»[

6、ekt]=step(Ez,0:1:20)%求E(z)对阶跃输入的响应ek=1.0000-0.9976-0.00390.0022-0.00110.0008-0.00060.0004-0.00030.0002-0.00020.0001-0.00010.0001-0.00000.0000-0.00000.0000-0.00000.0000-0.0000t=01234567891011121314151617181920»pp=spline(t，ek);»yy=ppval(pp，linspace(0,20));»plot(linspace(0,20)，

7、yy，'b');»[ykmaxtp]=max(yk);»tp=tp-1»ykp=ykmax-l»r=21;while(yk(r)>0.98&yk(r)<1.02)，r=r-1;end;»ts=r+1»gtext(’y(kT)’);gtextCe(kT)’);gtextCu(kT)’)»grid;»title(’StepResponse’)；StepResponse波形图: