| 查看: 2888 | 回复: 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的值吗？ 如果能拟合出来，请问该用什么方法啊？ |
回复此楼» 猜你喜欢
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
今天务委会开完了，明天出结果吗
已经有11人回复
-
有没有大神帮我看看基金代码
已经有9人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
放榜前的不淡定
已经有18人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 抛物型偏微分方程的数值解求解问题，谢谢！
已经有14人回复
- 请问这个4阶的非线性偏微分方程组在PDE中怎么解？等待高手解答。
已经有17人回复
- 请教mathematic问题？关于解微分方程组的和画图的一些问题。
已经有19人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- 关于微分方程组的求解
已经有4人回复
- 一阶微分方程求解，并优化参数值
已经有14人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- 这样一个偏微分方程有没有解析解？
已经有4人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab解微分方程
已经有10人回复
- 请问如何求解二元一阶微分方程组
已经有24人回复
- 复杂偏微分方程，望大侠指教……
已经有4人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】如何求解下列偏微分方程
已经有13人回复
- 【求助】急请微分方程高手看下这个方程组！
已经有14人回复
- 【求助】偏微分方程和常微分方程组成的方程组怎么求，大家帮帮忙啊
已经有3人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】一个二阶常微分方程的解法（重奖金币！）
已经有16人回复
- 【求助】有限元法求解二阶偏微分方程组
已经有6人回复
零下夏天
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1442.5
- 红花: 2
- 帖子: 391
- 在线: 134.1小时
- 虫号: 1979021
- 注册: 2012-09-06
- 专业: 其他无机非金属材料
2楼2013-03-26 23:41:03
零下夏天
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1442.5
- 红花: 2
- 帖子: 391
- 在线: 134.1小时
- 虫号: 1979021
- 注册: 2012-09-06
- 专业: 其他无机非金属材料
3楼2013-03-27 00:05:18
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2013-03-27 08:39:57
零下夏天
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1442.5
- 红花: 2
- 帖子: 391
- 在线: 134.1小时
- 虫号: 1979021
- 注册: 2012-09-06
- 专业: 其他无机非金属材料
5楼2013-03-27 09:42:02
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
6楼2013-03-27 10:02:25
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
零下夏天: 金币+7, ★★★很有帮助, 谢谢，提供了一个很好的思路 2013-03-28 11:43:19
感谢参与，应助指数 +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
零下夏天
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1442.5
- 红花: 2
- 帖子: 391
- 在线: 134.1小时
- 虫号: 1979021
- 注册: 2012-09-06
- 专业: 其他无机非金属材料
8楼2013-03-28 11:45:40
9楼2013-03-28 11:53:24
★
感谢参与，应助指数 +1
零下夏天: 金币+1, ★有帮助, 感谢回复，据说1stopt的拟合效果更好些 2013-03-29 00:13:51
感谢参与，应助指数 +1
零下夏天: 金币+1, ★有帮助, 感谢回复，据说1stopt的拟合效果更好些 2013-03-29 00:13:51
|
本帖内容被屏蔽 |
10楼2013-03-28 11:56:03
28