24小时热门版块排行榜    

查看: 2909  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

零下夏天

木虫 (正式写手)

[求助] 微分方程组求参数问题,求高人指教,感谢

问题如下:
dX/dt=-a*X*Y+b*Z*Z/2;
dY/dt=-a*X*Y-c*Y*Z/2;
dZ/dt=dY/dt-2*dX/dt;
X0=3,Y0=3,Z0=0;
如果只知道一系列 X 和 t 的关系,可以拟合出参数a,b,c的值吗?
如果能拟合出来,请问该用什么方法啊?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

拉布拉多犬

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
零下夏天: 金币+7, ★★★很有帮助, 谢谢,提供了一个很好的思路 2013-03-28 11:43:19
1、将第一和第2个方程带到第三个方程,这样第三个方程也变成了左边导数,右边函数的形式,这样就是三个标准形式的方程。
2.这个方程组是非线性的,在一定的a、b、c取值范围内,系统的运动可能是混沌的,所以楼主应该先借助数值方法研究该方程解的特性,如果确实存在混沌区,那么在混沌区的x是没有办法通过拟合而确定参数abc的,因为混沌运动有初值敏感性。而对于非混沌区的x则可以通过拟合确定abc
3.至于用什么方法拟合,则方法比较多,如最小二乘等等。如果处在非混沌区,我认为楼主应该首先对abc的取值有个范围估计,然后在此范围内数值方法(如龙格库塔法)进行运算,找到最逼近你数据的abc。
7楼2013-03-27 20:18:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

零下夏天

木虫 (正式写手)

对了,X,Y,Z>=0的
2楼2013-03-26 23:41:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

零下夏天

木虫 (正式写手)

另外,一般的a>b,a>c,a,b,c>0
3楼2013-03-27 00:05:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
零下夏天: 金币+2, 有帮助, 我昨天搜了许多相关的帖子,知道能用这个软件拟合,但是我很想知道原理,您知道吗,方便解释下吗? 2013-03-27 09:39:30
用下1stOpt,很容易求解这种微分方程拟合问题。
4楼2013-03-27 08:39:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:23 by hZiFeudZyoGR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:13 by hZiFeudZyoGR
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 19:01 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:51 by hZiFeudZyoGR
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:50 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +8 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 17:11 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 14/700 2026-08-30 17:01 by hZiFeudZyoGR
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 4/200 2026-08-30 16:03 by sundiv
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 14:03 by GkCVyODlx4mR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 13:54 by GkCVyODlx4mR
[基金申请] 国社科又开始会评了,不知道这次命运如何 +3 雨打竹帘 2026-08-30 4/200 2026-08-30 12:35 by 白大多亮
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:14 by l0VvVHGBGRLv
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 06:00 by ZPa0EcMwuECS
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
信息提示
请填处理意见