试题详情
简答题 已知系统结构图如下，试采用四阶龙格－库塔法求系统输出响应。
  • 程序清单如下：
    r=2；
    [numo，deno]=series（[1，1]，[1，4]，4，[1，3，2]）；
    [numc，denc]=cloop（numo，deno）；
    [A，B，C，D]=tf2ss（numc，denc）；
    Tf=input（‘仿真时间Tf=’）；
    h=input（‘计算步长h=’）；
    X=zeros（size（A，1），1）；
    y=0；t=0；
    fori=1：Tf/h
    K1=A*X+B*r；
    K2=A*（X+h*K1/2）+B*r；
    K3=A*（X+h*K2/2）+B*r；
    K4=A*（X+h*K3）+B*r；
    X=X+h*（K1+2*K2+2*K3+K4）/6；
    y=[y，C*X]；t=[t，t（i）+h]；
    end
    plot（t，y）；

  • 关注下方微信公众号，在线模考后查看

热门试题