ex7_11.m

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

G=tf(10,[1,10,35,50,24]);
[k,L,T]=getfod(G); 
[Gc1,Kp1]=cohenpid(1,[k,L,T,10]); [Kp1]

[Gc2,Kp2,Ti2]=cohenpid(2,[k,L,T,10]); [Kp2,Ti2]

[Gc3,Kp3,Ti3,Td3]=cohenpid(3,[k,L,T,10]); [Kp3,Ti3,Td3]

[Gc4,Kp4,Ti4,Td4]=cohenPID(5,[k,L,T,10]); [Kp4,Td4]

t=0:.1:10;
G_c1=feedback(G*Gc1,1); y=step(G_c1,t); 
G_c2=feedback(G*Gc2,1); y=[y step(G_c2,t)]; 
G_c3=feedback(G*Gc3,1); y=[y step(G_c3,t)]; 
G_c4=feedback(G*Gc4,1); y=[y step(G_c4,t)]; 
plot(t,y)