ex742a1.m

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

syms s T K z t
Gs=(1/s)*K/s/(0.1*s+1);
ft=ilaplace(Gs);ftt=subs(ft,t,t*T);
GoGpZ=(1-z^-1)*ztrans(ftt);GoGpZ=simplify(GoGpZ);
GoGpZ=subs(GoGpZ,T,0.1);GoGpZ=subs(GoGpZ,K,10);
[dnum,dden]=numden(GoGpZ);
dnum=sym2poly(dnum);dden=sym2poly(dden);
disp('广义对象的Z传递函数GoGp(z)=');
[z p k]=tf2zp(dnum,dden);GoGpZ=zpk(z,p,k,0.1)