ex3_2.m

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

% ex3_2
% 二阶环节的阶跃响应
figure('pos',[120,130,300,250],'color','w')
axes('pos',[0.1 0.15 0.8 0.72]);
syms  s 
for zeta=[0.1:0.2:0.9,1:.5:2.0];
    wn=0.15; 
    wn=sym(num2str(wn));
    zet=sym(num2str(zeta));
    if zeta==1
    ezplot(ilaplace(wn^2/s/(s+wn)^2),[0 100]); 
    hold on;
    text(26,0.85,'1.0');
    else
    ezplot(ilaplace(wn^2/s/(s^2+2*zet*wn*s+wn^2)),[0 100]); 
    hold on;
    if zeta==0.1,text(18,1.85,'0.1');end;
    if zeta==0.5,text(22,1.25,'0.5');end;
    if zeta==2,text(29.5,0.6,'2.0');end;
    end;
end;
grid; 
axis([0 100 0 2]);
title('ζ:0.1,0.3,0.5,0.7,0.9,1,1.5,2');