24小时热门版块排行榜    

查看: 2440  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 零下夏天 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 江苏省农科院招调剂1名 +3 Qwertyuop 2026-03-01 3/150 2026-03-01 23:18 by aaadim
[考研] 0854复试调剂 276 +3 wmm9 2026-03-01 3/150 2026-03-01 23:13 by 热情沙漠
[考研] 0856调剂 +5 刘梦微 2026-02-28 5/250 2026-03-01 22:30 by wang_dand
[考研] 275求调剂 +3 明远求学 2026-03-01 3/150 2026-03-01 22:29 by 刘兵
[考研] 26考研报考西工大材料308分求调剂 +3 weizhong123 2026-03-01 3/150 2026-03-01 21:42 by 公瑾逍遥
[考研] 一志愿郑大材料学硕298分,求调剂 +3 wsl111 2026-03-01 3/150 2026-03-01 21:38 by ms629
[考研] 274求调剂 +3 cgyzqwn 2026-03-01 6/300 2026-03-01 21:24 by cgyzqwn
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
[考研] 材料学硕318求调剂 +9 February_Feb 2026-03-01 11/550 2026-03-01 19:47 by 无懈可击111
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考研] 0856材料求调剂 +11 hyf hyf hyf 2026-02-28 12/600 2026-03-01 18:57 by 18137688336
[考研] 材料学调剂 +9 提神豆沙包 2026-02-28 11/550 2026-03-01 18:15 by ms629
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 285求调剂 +8 满头大汗的学生 2026-02-28 8/400 2026-03-01 16:47 by caszguilin
[考研] 311求调剂 +6 亭亭亭01 2026-03-01 6/300 2026-03-01 15:41 by 324616
[考研] 295复试调剂 +3 简木ChuFront 2026-03-01 3/150 2026-03-01 14:27 by zzxw520th
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 材料284求调剂,一志愿郑州大学英一数二专硕 +10 想上岸的土拨鼠 2026-02-28 10/500 2026-03-01 14:12 by yc258
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
信息提示
请填处理意见