已知系统结构图如下，试采用四阶龙格－库塔法求系统输出响应。 
正确答案:
程序清单如下：
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）；
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）；
答案解析:有
微信扫一扫手机做题