program2_2c_2.asv

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

%MATLAB Program2_2c_2
A=[0 1 0;0 0 1;-2 -4 -3];
B=[1 0;0 1;-1 1];
C=[0 1 -1;1 2 1];
D=[0 0; 0 0];
T=0:0.01:4;
[YY1,XX1]=step(A,B,C,D,1,t);
[YY2,XX2]=step(A,B,C,D.2,t);
y1=YY1(:,1)+YY2(:,1);
y2=YY1(:,1)+YY2(:,1);
subplot(2,1,1),plot(t,y1);ylabel('To:Y(1)'),grid
subplot(2,1,2),plot(t,y2);ylabel('To:Y(2)'),xlabel('time(secgrid