paper1.m

来自「华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 」· M 代码 · 共 7 行

syms s kc ti tp kp torp Gc Gp G y yy
Gc=kc*(1+1/ti/s);Gp=kp*exp(-s*torp)/(1+tp*s);
G=Gc*Gp/(1+Gc*Gp);
y=ilaplace(G/s);
yy=subs(y,{kc,ti,kp,tp,torp},{1 1 1 10 2});
%yt=subs(yy,[0:100])
[num,den]=numden(yy)