| 查看: 2922 | 回复: 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的值吗？ 如果能拟合出来，请问该用什么方法啊？ |
回复此楼» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有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
投票: