ex6_24.m

来自「华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 」· M 代码 · 共 23 行

% ex6_24
% 确定控制器放大倍数Kc和积分时间Ti，使系统有4：1衰减比输出响应曲线
figure('pos',[270,50,480,280],'color','w');
axes('pos',[0.1,0.14,.8,.8]);
G=zpk([],[-1,-2,-3],6); 
[Gmm1,Pmm,Wcgm,Wcpm]=freqm(G,0.221,'margin');
G=tf(G);
for i=1:2;
    k(1)=2/1.2;
    k(2)=Gmm1/1.2;
    num=[6*k(i)];
    den=[G.den{1}+[zeros(1,length(G.den{1})-1) 6*k(i)] 0];
    G1=tf(num,den);
    [Gmm(i),Pmm,Wcgm,Wcpm]=freqm(G1,0.221,'margin');
end;
Gmm(3)=1/pi;k(1,3)=k(1,2)
for i=1:3,
    tm=Gmm(i);
    Gc=tf(k(1,i)*[1 tm],[1 0]);
    Gloop=feedback(G*Gc,1);
    step(Gloop);
    hold on;
end;