文档介绍：仿真实验报告
3、采用四阶龙格库塔法求高阶系统阶单位跃响应曲线的数值解。
，，
4、 自学OED45指令用法，并求解题2中二阶系统的单位阶跃响应。
程序：y=RKT2(,5,[0 1;-100 -10],[0;100],[0;0])；
y=RKT3(,20,[0,1,0;0,0,1;-20,-102,-],[0;0;20],[0;0;0])；
[t,y]=ode45('func',[0,1],[0,0]);plot(t,y)；
M文件
function y=RKT2(h,T,A,B,Y)
YY=Y;
for t=0:h:T;
if(t<=1)
u=1;
else
u=0;
end
K1=A*YY+B*u;
K2=A*[YY+h/2*K1]+B*u;
K3=A*[YY+h/2*K2]+B*u;
K4=A*[YY+h*K3]+B*u;
YY=YY+h/6*(K1+2*K2+2*K3+K4);
Y=[Y YY];
end
y=Y(1,:);
t=0:h:(T+h);
function y=RKT3(h,T,A,B,Y)
YY=Y;
for t=0:h:T;
K1=A*YY+B;
K2=A*[YY+h/2*K1]+B;
K3=A*[YY+h/2*K2]+B;
K4=A*[YY+h*K3]+B;
YY=YY+h/6*(K1+2*K2+2*K3+K4);
Y=[Y YY];
end
y=Y(1,:);
t=0:h:(T+h);
plot(t,y)
function yp=func(t,y)
yp=[0;0];
yp(1)=y(2);
yp(2)=100-10*y(2)-100*y(1);
end
实验结果图像：
第二次上机任务
1、试用simulink方法解微分方程，并封装模块，输出为。得到各状态变量的时间序列，以及相平面上的吸引子。
参数入口为的值以及的初值。（其中，以及初值分别为） 提示：模块输入是输出量的微分。
模块的封装图：
封装内部结构图：
各状态的时间曲线图：
吸引子：