m6_9_1.m

来自「MATLAB控制系统仿真与实例详解（附书光盘）」· M 代码 · 共 13 行

wn=1;
zetas=[0:0.1:1,2,3,5];
t=0:0.1:12; 
hold on 
for i=1:length(zetas) 
Gc=tf(wn^2,[1,2*zetas(i)*wn,wn^2]); 
step(Gc,t) 
end 
hold off
title('阶跃响应曲线');
grid