ex741.m

来自「详细介绍了Mathlab的使用方法与技巧」· M 代码 · 共 14 行

syms s T K z t k
Gs=(1/s)*K/s/(0.1*s+1);                                  % 
ft=ilaplace(Gs);ftt=subs(ft,t,k*T);                           %求Gs的反变换， 
GoGpZ=(1-z^-1)*ztrans(ftt);GoGpZ=simplify(GoGpZ);         %求 ，并化简
Ez=(1-z^-1)^2;                                         % 
Dcz=(1-Ez)/(Ez*GoGpZ);Dcz=simplify(Dcz);                %(7.4.4)式，并化简
Dcz=subs(Dcz,T,0.1);Dcz=subs(Dcz,K,10);                  % ， 
[dnum,dden]=numden(Dcz);                           %取得 的分子和分母的多项式
dnum=sym2poly(dnum);dden=sym2poly(dden);           %转化为多项式系数向量
disp('单位速度输入时的有限拍控制器Dc(z)=');
[z p k]=tf2zp(dnum,dden);Dcz=zpk(z,p,k)
disp('系统脉冲传递函数D(z)=C(z)/R(z)');
Dz=1-Ez;Dz=simplify(Dz)