MATLAB 自动控制原理 传递函数的单位阶跃响应曲线绘制

    wc=1000;
    zuni1=0.5;zuni2=0.6;zuni3=0.7;zuni4=0.8;
    num1=[wc^2];
    den1=[1 2*zuni1*wc wc^2];
    den2=[1 2*zuni2*wc wc^2];
    den3=[1 2*zuni3*wc wc^2];
    den4=[1 2*zuni4*wc wc^2];
    sys1=tf(num1,den1);sys2=tf(num1,den2);sys3=tf(num1,den3);sys4=tf(num1,den4);
    t=0:0.0001:0.04;
    figure(1);
    c1=step(sys1,t);c2=step(sys2,t);c3=step(sys3,t);c4=step(sys4,t);
    plot(t,c1);hold on;plot(t,c2);plot(t,c3);plot(t,c4);
    legend('zuni1=0.5','zuni2=0.6','zuni3=0.7','zuni4=0.8');
    xlabel('时间/单位:秒');ylabel('单位阶跃响应幅值');
    title('wc=1000验证阻尼系数对系统