24小时热门版块排行榜    

查看: 2437  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 零下夏天 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 284求调剂 +6 天下熯 2026-02-28 6/300 2026-03-01 14:19 by Ducount.Y
[考研] 0805总分292,求调剂 +3 幻想之殇 2026-03-01 3/150 2026-03-01 14:11 by yc258
[考研] 0856材料求调剂 +3 麻辣鱿鱼 2026-02-28 3/150 2026-03-01 14:06 by yc258
[考研] 303求调剂 +3 今夏不夏 2026-03-01 3/150 2026-03-01 14:01 by yc258
[考研] 材料化工调剂 +8 今夏不夏 2026-03-01 9/450 2026-03-01 13:25 by yc258
[考研] 化工299分求调剂 一志愿985落榜 +4 嘻嘻(*^ω^*) 2026-03-01 4/200 2026-03-01 13:15 by wang_dand
[考研] 290求调剂 +8 材料专硕调剂; 2026-02-28 9/450 2026-03-01 12:46 by 闭眼看蓝天
[考研] 0856化工专硕求调剂 +6 董boxing 2026-03-01 6/300 2026-03-01 12:45 by houyaoxu
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考研] 求调剂 +5 repeatt?t 2026-02-28 5/250 2026-03-01 11:43 by 王伟要上岸啊
[考研] 高分子化学与物理调剂 +5 好好好1233 2026-02-28 9/450 2026-03-01 10:59 by fengyu211
[考研] 寻找调剂 +4 LYidhsjabdj 2026-02-28 4/200 2026-03-01 10:56 by sunny81
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[考研] 0856求调剂285 +6 吕仔龙 2026-02-28 6/300 2026-03-01 10:03 by wang_dand
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
[考研] 298求调剂 +5 axyz3 2026-02-28 5/250 2026-03-01 06:45 by 刘兵
[考研] 292求调剂 +3 yhk_819 2026-02-28 3/150 2026-02-28 21:57 by gaoxiaoniuma
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
信息提示
请填处理意见