| 查看: 2490 | 回复: 6 | ||
| 【悬赏金币】回答本帖问题，作者自制奶酪将赠送您 100 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
用matlab求动力学参数 已有3人参与
|
||
|
本人MATLAB小白，最近审稿人要求做个动力学模型，我不会用MATLAB拟合求速率常数，希望木虫的大神可以教我，谢谢。 (dC_A)/dt=-k1*CA-k7*CA (dC_B)/dt=k1*CA-k2*CB- k5*CB (dC_C)/dt=k2*CB-k3*CC- k6*CC (dC_D)/dt=k3*CC-k4*CD (dC_E)/dt= k4*CD (dC_F)/dt=k5*CB+k6*CC+ k7*CA 此外还采用Nelder-Mead 方法通过最小化目标函数优化速率常数 F_obj=∑_m▒{∑_(n=1)^(n_p)▒[C_(m,n,calc)-C_(m,n,exp) ] ^2 } C_(m,n,calc)是给定的时间n下组分的实验浓度；C_(m,n,exp)是给定的时间n下组分的计算浓度 第一组： T= (0.5, 1, 1.5, 2, 3, 4, 5, 6) A= (18.7, 6.5, 0.7, 2.1, 1.7, 1.5, 1.2, 0.2) B= (48.9, 31.3, 15.3, 6.9, 2.2, 1.8, 1.7, 0.5) C= (19.3, 34.7, 41.2, 40.9, 37.2, 31.7, 26.1, 19.4) D= (9.7, 19.6, 28.7, 34.5, 38.4, 38.2, 36.3, 33.4) E= (0.2, 2.6, 5.7, 9.3, 17.1, 24.8, 30.5, 36.7) 第二组： T= (0.5, 1, 1.5, 2, 3, 4, 5, 6) A= (23.7, 7.2, 1.7, 1.1, 0.7, 0.5, 0.2, 0.2) B= (55.4, 44.3, 27.3, 11.9, 1.2, 1.1, 0.7, 0.5) C= (14.3, 33.3, 43.2, 49.3, 51.5, 44.7,36.1, 30.4) D= (3.7, 11.3, 19.6, 26.4, 33.4, 35.2, 34.1, 32.4) E= (0.2, 1.2, 3.7, 7.5, 13.7, 18.3, 23.2, 27.3) 第三组： T= (0.5, 1, 1.5, 2, 3, 4, 5, 6) A= (33.7, 15.2, 3.7, 1.1, 0.7, 0.5, 0.2, 0.2) B= (49.7, 53.4, 44.3, 28.3, 8.7, 1.1, 0.7, 0.5) C= (9.3, 20.8, 32.2, 43.7, 53.5, 55.7, 51.4, 46.9) D= (2.7, 7.9, 13.7, 18.4, 25.4, 28.2, 29.4, 30.4) E= (0, 0.2, 2.1, 4.3, 8.7, 12.3, 15.7, 19.3) |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有123人回复
-
物理学I论文润色/翻译怎么收费?
已经有84人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有0人回复
5楼2021-01-05 21:42:43
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31030.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5603
- 在线: 856.8小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
2楼2021-01-05 21:04:32
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
3楼2021-01-05 21:21:18
4楼2021-01-05 21:41:39
