24小时热门版块排行榜    

查看: 2958  |  回复: 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的回帖

零下夏天

木虫 (正式写手)

对了,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的回帖

零下夏天

木虫 (正式写手)

引用回帖:
4楼: Originally posted by dingd at 2013-03-27 08:39:57
用下1stOpt,很容易求解这种微分方程拟合问题。

就是他用的数值分析或者最优化里的相关的定理什么的?现在我觉得自己知识有点不够用,想再充充电
5楼2013-03-27 09:42:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


零下夏天: 金币+1, 有帮助, 谢谢,看来还是知识太欠缺了 2013-03-28 11:22:59
全局最优化不是一两句能说的清,算法实现、效果等千差万别,先找相关书籍看看。
6楼2013-03-27 10:02:25
已阅   回复此楼   关注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的回帖

零下夏天

木虫 (正式写手)

引用回帖:
7楼: Originally posted by 拉布拉多犬 at 2013-03-27 20:18:41
1、将第一和第2个方程带到第三个方程,这样第三个方程也变成了左边导数,右边函数的形式,这样就是三个标准形式的方程。
2.这个方程组是非线性的,在一定的a、b、c取值范围内,系统的运动可能是混沌的,所以楼主应 ...

就是还有一点不明白,怎么用数值方法研究方程解的特性,是不是存在混沌区,什么书上有这样的内容,可以推荐一下吗?
8楼2013-03-28 11:45:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

拉布拉多犬

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
零下夏天: 金币+5, 有帮助, 谢谢了,先学学看 2013-03-29 00:12:52
参考《非线性动力学》方面的书,有混沌的介绍,再参考《数值计算方法》方面的书,看一阶常微分方程组的数值解法,龙格库塔是个不错的方法。
9楼2013-03-28 11:53:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

456895

禁虫 (正式写手)


感谢参与,应助指数 +1
零下夏天: 金币+1, 有帮助, 感谢回复,据说1stopt的拟合效果更好些 2013-03-29 00:13:51
本帖内容被屏蔽

10楼2013-03-28 11:56:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 零下夏天 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +6 两三岁ss 2026-09-14 6/300 2026-09-15 20:17 by 只吃肉葡萄
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:23 by CgyNCDVNhGVg
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:11 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:35 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 16:15 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-15 15:50 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-15 14:49 by VMDqgnlKimtZ
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:38 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 05:22 by 5BDX0d0WFp7t
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
信息提示
请填处理意见