草稿本.m

来自「PID控制器设计」· M 代码 · 共 12 行

num=10;
den=[1,10,25,50,24];
G=tf(num,den);
vv=getfod(G);    %用模型转换
[Gc1,Kp1]=cohen(1,vv);   %可以将cohen换成ziegler，下同
[Gc2,Kp2,Ti2]=cohen(2,vv);
[Gc3,Kp3,Ti3,Td3]=cohen(3,vv);
[Gc4,Kp4,Ti4,Td4]=cohen(4,vv);
G_c1=feedback(G*Gc1,1);step(Gc1);
G_c2=feedback(G*Gc2,1);step(Gc2);
G_c3=feedback(G*Gc3,1);step(Gc3);
G_c4=feedback(G*Gc4,1);step(Gc4);