ex3_2a.m

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

% ex3_2a
% 二阶环节的脉冲响应
figure('pos',[120,130,300,250],'color','w')
axes('pos',[0.1 0.15 0.8 0.72]);
syms  s 
for zeta=[0:0.2:0.8,1:.5:2.0];
wn=0.15; 
wn=sym(num2str(wn));
zet=sym(num2str(zeta));
   if zeta==1
   ezplot(ilaplace(wn^2/(s+wn)^2),[0 100]); 
   hold on;
   else
   ezplot(ilaplace(wn^2/(s^2+2*zet*wn*s+wn^2)),[0 100]); 
   hold on;
end;
end;
grid; 
axis([0 100 -.15 .15]);
title('ζ:0,0.2,0.4,0.6,0.8,1,1.5,2');