自动控制原理中传递函数的Matlab编程求解性能指标

    %%  系统建模
    num=[67 4596];
    den=[1 104 903 4144];
    G=tf(num,den);
    %% 单位阶跃响应表达式计算
    num1=[1];
    den=[1 0];
    shuru1=tf(num,den);
    G_jieyue=shuru1*G
    syms s;
    G_jieyues=(4489*s^2+615864*s+2.112e07)/(s^4+104*s^3+903*s^2+4144*s)
    G_jieyuet=ilaplace(G_jieyues)
    %% 单位脉冲响应表达式计算
    G_maichong=G
    G_maichongs=(67*s+4596)/(s^4+104*s^3+903*s^2+4144*s);
    G_maichongt=ilaplace(G_maichongs)
    %% 单位阶跃响应绘图
    figure(1);
    t=0:0.01:5;
    c=step(G,t);
    plot(t,c);
    Gss=dcgain(G)
    grid
    title('Unit-Step Response of G(s)=(67s+4596)/(s^3+104s^2+903s+414