关于微分方程分段求解问题
解微分方程时，如果设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);
方法一：
=ode45(@rigid,,);
方法二：
=ode45(@rigid,,);
=ode45(@rigid,,);
回复 #1 thuwb 的帖子
这就是非线性系统的响应中对初值的敏感性啊！这个问题如果你能早于Lorenz发现的话，你就是混沌之父了！
页:
[1]