ziegler.m

PID控制器设计

%Ziegler-Nichols整定函数
function[Gc,Kp,Ti,Td]=ziegler(key,vars)
Kp=[];Ti=[];Td=[];
if length(vars)==3
    K=vars(1);L=vars(2);T=vars(3);
    a=K*L/T;
    if key==1   %P控制器
        Kp=1/a;
    elseif key==2   %PI控制器
        Kp=0.9/a;Ti=3.3*L;
    elseif key==3     %PID控制器
        Kp=1.2/a;Ti=2*L;;Td=L/2;
    end
elseif length(vars)==2
    Kc=vars(1);Tc=vars(2);
    if key==1
        Kp=0.5*Kc;
    elseif key==2
        Kp=0.4*Kc;Ti=0.8*Tc;
    elseif key==3
        Kp=0.6*Kc;Ti=0.5*Tc;;Td=0.12*Tc;
    end
end
switch key
    case 1,Gc=Kp;
    case 2,Gc=tf(Kp*[Ti,1],[Ti,0]);
    case 3
        nn=[1.1*Kp*Ti*Td,Kp*(Ti+0.1*Td),Kp];
        dd=Ti*[0.1*Td,1,0];
        Gc=tf(nn,dd);
end