声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 802|回复: 1

[综合讨论] 关于微分方程分段求解问题

[复制链接]
发表于 2007-4-10 12:02 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
解微分方程时,如果设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)]);
回复
分享到:

使用道具 举报

发表于 2007-4-10 17:12 | 显示全部楼层

回复 #1 thuwb 的帖子

这就是非线性系统的响应中对初值的敏感性啊!这个问题如果你能早于Lorenz发现的话,你就是混沌之父了!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 04:22 , Processed in 0.067364 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表