ex4_26.m

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

% ex4_26
% 时滞系统的观测补偿控制系统稳定性分析
figure('pos',[50,50,200,200],'color','w');
axes('pos',[0.12,0.14,0.75,0.75]);
pn=[1];
pd=[10 1];
Gp=tf(pn,pd);
Gp.inputdelay=2;
cn=[10 1];
cd=[10 0];
Gc=tf(cn*9,cd);
mn=[1.1];
md=[11 1];
Gm=tf(mn,md);
w=logspace(-1,0.8,100);
[r1,i1]=nyquist(Gc*Gp,w);
r11(1,:)=r1(1,:);
i11(1,:)=i1(1,:);
plot(r11,i11);
hold on;
for kk=[0.1 0.5 5];
    kn=[50 1];
    kd=[50 0];
    Gk=tf(kn*kk,kd);
    nyquist(-(1+inv(Gk*Gm)));
end;
axis([-2 0.5 -5 1])