| 查看: 2281 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
动力学方程组参数求解问题 已有3人参与
|
|||
|
本人动力学方程如下： dx1/dt=k21x2+k31x3+k41x4+k51x5+k61x6+k71x7+k81x8 dx2/dt=-k21x2+k32x3+k42x4+k52x5+k62x6+k72x7+k82x8 dx3/dt=-k31x3-k32x3+k43x4+ k53x5+k63x6+k73x7+k83x8 dx4/dt=-k41x4-k42x4-k43x4+ k54x5+k64x6+k74x7+k84x8 dx5/dt=-k51x5-k52x5-k53x5-k54x5-k56x52-k57x52+k65x6+k75x7+k85x8 dx6/dt=-k61x6-k62x6-k63x6-k64x6-k65x6-k67x62+k56x52+k76x7+k86x8 dx7/dt=-k71x7-k72x7-k73x7-k74x7-k75x7-k76x7-k78x72+k57x52+k67x62+k87x8 dx8/dt=-k81x8-k82x8-k83x8-k84x8-k85x8-k86x8-k87x8-k89x82+k78x72 dx9/dt=k89x82 方程中有多个缩合反应为表观二级（即平方），已上传图片供参考 k为动力学参数，k大于0，共33个 x为质量分数，x＞0，九个x之和为1 已知数据如下： 温度 时间 x1 x2 x3 x4 x5 x6 x7 x8 x9 430 0.00 0.0192 0.0236 0.1006 0.3284 0.1647 0.1327 0.1981 0.0319 0.0008 430 20.00 0.0287 0.0503 0.1313 0.3156 0.1391 0.1206 0.1755 0.0366 0.0022 430 40.00 0.0349 0.0858 0.1827 0.3058 0.1203 0.0805 0.1476 0.0384 0.0039 430 60.00 0.0423 0.1250 0.2086 0.2950 0.1005 0.0624 0.1180 0.0428 0.0053 本人想编写一个MATLAB程序实现以下功能： 设k初值均为0，求解方程中的33个参数，将求得的参数代入方程中，根据t=0min的数据求得t=20、40、60min的x计算值，将x计算值与x实验值比较，如果误差＜5%，得到的参数就是最终参数，如果误差大于5%，就将得到的参数最为初值 循环上述计算，直到 误差＜5%。 哪位大侠能帮忙编写一下，小女子不胜感激呢~  动力学方程组.jpg |
回复此楼» 猜你喜欢
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有267人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
广州大学光电信息工程专业调剂，招收物理学专业学生
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用matlab最优化方法拟合获得多个动力学参数中的问题
已经有4人回复
- 有关伪一级动力学和伪二级动力学方程
已经有9人回复
- 用matlab求解方程
已经有6人回复
- 求助动力学方程？？？
已经有3人回复
- 求助用matlab拟合动力学方程
已经有9人回复
- 如何求解带有固定参数的超定方程组？
已经有21人回复
- 求助，matlab自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 【求助】非线性方程组的求解问题
已经有6人回复
- 求助matlab拟合方程组
已经有12人回复
- 动力学方程参数估计方法
已经有14人回复
- SCR反应动力学方程求解的问题
已经有4人回复
- 动力学参数拟合
已经有26人回复
- 请教matlab反应动力学参数估计遇到的问题，谢谢
已经有15人回复
- 车辆动力学模型的解法
已经有14人回复
- matlab m文件运行出错
已经有4人回复
- matlab-常微分方程参数估计
已经有12人回复
- lsqnonlin函数拟合微分方程组参数拟合问题
已经有10人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- matlab非线性参数拟合问题
已经有7人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】多元方程组求解及拟合
已经有9人回复
- 【求助】催化反应动力学
已经有5人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
6楼2014-09-03 10:17:13
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68129.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2014-09-02 14:26:53
Cle_key
铁虫 (小有名气)
- 应助: 22 (小学生)
- 金币: 1560.8
- 红花: 2
- 帖子: 178
- 在线: 28.9小时
- 虫号: 2625972
- 注册: 2013-08-31
- 性别: GG
- 专业: 内流流体力学
3楼2014-09-02 17:38:33
polypro
木虫 (正式写手)
- 应助: 197 (高中生)
- 金币: 3921.3
- 散金: 160
- 红花: 22
- 帖子: 512
- 在线: 62.6小时
- 虫号: 1809438
- 注册: 2012-05-11
- 性别: GG
- 专业: 理论和计算化学
4楼2014-09-03 08:20:26
