ex3_4.m

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

% ex3_4
% 二阶环节的斜坡响应
figure('pos',[120,130,300,250],'color','w')
axes('pos',[0.1 0.15 0.8 0.72]);
for zeta=[0:0.25:0.5,1:0.5:2];
wn=0.15; num=wn^2;
den1=[1 2*zeta*wn wn^2];
den=[den1 0 0];
[r,p,k]=residue(num,den) 
syms  s t;
if zeta==1,
    tt=[1 t 1 t]';
else 
    tt=[1 1 1 t]';
end;
y=sum(tt.*ilaplace(r./(s-p)))
ezplot(y,[0 100]);
hold on;
end;
grid; 
axis([0 100 0 100]);
title(['(',char(poly2sym(num,'s')),') / (',char(poly2sym(den1,'s')),')']);
line([0 100],[0 100],'linestyle','--');