ex7_6.m

来自「反馈控制系统分析与设计,虽然有点抽象」· M 代码 · 共 8 行

G=tf(10,[1,10,35,50,24]); N=10;
[Kc,Pm,wc,wp]=margin(G); Tc=2*pi/wc;
[Gc1,Kp1,Ti1,Td1]=ziegler(3,[Kc,Tc,N]); 
[Gc2,Kp2,Ti2,Td2,H]=ziegler(4,[Kc,Tc,N]); [Kp2,Ti2,Td2]

G_c1=feedback(G*Gc1,1); step(G_c1); hold on
G_c2=feedback(G*Gc2,H); step(G_c2)