ex7_12.m

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

G=tf(10,[1,10,35,50,24]);
[k,L,T]=getfod(G); 
[Kc,p,wc,m]=margin(G);  Tc=2*pi/wc;
[Gc,Kp,Ti,Td,beta,H]=rziegler([k,L,T,10,Kc,Tc]);
[Kp,Ti,Td,beta]

G_c=feedback(G*Gc,H); step(G_c); hold on
[Gc,Kp1,Ti1,Td1]=ziegler(3,[k,L,T,10]); [Kp1,Ti1,Td1]

G_c=feedback(G*Gc,1); step(G_c);