program2_1.m

此为matlab的代码。包括数值o45算法及o23算法

%MATLAB Program2_1a

%系统的状态转移矩阵
syms s t
A=[-23/3 -22/3 -2/3;1 0 0;0 1 0];
FS=inv(s*eye(3)-A);
eAt=ilaplace(FS,s,t);
eAt=simplify(eAt)

%系统阶跃响应解析解
S=dsolve('Dv=-23/3*v-22/3*w-2/3*z+1,Dw=v,Dz=w','v(0)=0,w(0)=0,z(0)=0');
x1=S.v
x2=S.w
x3=S.z