关于微分方程分段求解
解微分方程时,如果设t0=0,tf=1000,
为什么用方法一:直接求t=t0~tf的值
与用方法二分段求解t=0~500,t=500~1000算出来的值不一样呢?
方法二种将第一时间段算出的t=500时的值作为第二时间段的初始值。
哪位大侠解释一下,谢谢!
function dy=rigid(t,y)
dy=zeros(3,1);
dy(1)=y(2)*y(3);
dy(2)=-y(1)*y(3);
dy(3)=-0.51*y(1)*y(2);
方法一:
[T,Y]=ode45(@rigid,[0:10:1000],[0 1 1]);
方法二:
[T1,Y1]=ode45(@rigid,[0:10:500],[0 1 1]);
[T2,Y2]=ode45(@rigid,[500:10:1000],[Y1(51,1) Y1(51,2) Y1(51,3)]);