| 查看: 1730 | 回复: 6 | ||
| 本帖产生 1 个 计算强帖 ，点击这里进行查看 | ||
v1iv1ia1n至尊木虫 (知名作家)
|
[求助]
matlab参数拟合
|
|
|
求助： 不太会用matlab，请路过的大侠拔刀相助！！ 想求ku、kep、kem和km这几个参数，微分方程如下： dCtot/dt = (ku*1.5)-(kep*Cp)-(kem*Cm) dCp/dt = (ku*1.5) -(km*Cp)-(kep*Cp) dCm/dt = (km*Cp)-(kem*Cm) 已知数据为Ctot和t： t Ctot 0.5h 1 0.326 2 0.289 3 0.443 AVE 0.352666667 1h 1 0.362 2 0.38 3 0.417 AVE 0.386333333 2h 1 0.43 2 0.401 3 AVE 0.4155 4h 1 0.292 2 0.368 3 0.538 AVE 0.399333333 8h 1 0.507 2 0.492 3 0.513 AVE 0.504 12h 1 0.275 2 0.183 3 0.202 AVE 0.22 24h 1 0.111 2 0.122 3 0.094 AVE 0.109 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : data.xls
2012-04-13 16:06:47, 39.5 K
» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有167人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有11人回复
-
散金币祈福
已经有94人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab拟合方程参数时初值的选择
已经有15人回复
- matlab拟合求值
已经有17人回复
- 在matlab中拟合曲面的表达式
已经有29人回复
- matlab拟合拟均相动力学参数
已经有6人回复
- 需要用matlab拟合数据的外行人急需帮助
已经有10人回复
- matlab非线性参数拟合问题
已经有7人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
- 如何用MATLAB拟合指数函数
已经有3人回复
- 【求助】matlab曲面拟合表达式
已经有10人回复
- 【求助】如何使得matlab拟合的曲线强制经过一个点
已经有8人回复
- 【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】
已经有15人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
cenwanglai: 金币+5, 计算强帖+1, 非常感谢热心应助，欢迎常来交流探讨~ 2012-04-14 21:17:21
v1iv1ia1n: 金币+30, ★★★★★最佳答案, 谢谢大侠，你太牛了，我好好研究研究！！！ 2012-04-14 23:32:28
感谢参与，应助指数 +1
cenwanglai: 金币+5, 计算强帖+1, 非常感谢热心应助，欢迎常来交流探讨~ 2012-04-14 21:17:21
v1iv1ia1n: 金币+30, ★★★★★最佳答案, 谢谢大侠，你太牛了，我好好研究研究！！！ 2012-04-14 23:32:28
|
应该属于微分方程拟合了，Cp、Cm的初值都知道吗？下面是1stOpt求解代码： 结果（有多组）： 均方差(RMSE): 7.21477021187873E-12 残差平方和(SSE): 7.28740728942977E-22 相关系数(R): 1 相关系数之平方(R^2): 1 决定系数(DC): 1 参数 最佳估算 -------------------- ------------- ku 0.0276775305197549 kep -0.0456535422528521 kem 0.00133632074107808 km 0.514129491172468 cp初值(数据文件 - 1) -13.5280340784502 cm初值(数据文件 - 1) -306.857892086343 cp初值(数据文件 - 2) -1.26798177439112 cm初值(数据文件 - 2) -16.1124930625247 cp初值(数据文件 - 3) -3.11756100776893 cm初值(数据文件 - 3) -30.9936055993614 cp初值(数据文件 - 4) -6.62174164284473 cm初值(数据文件 - 4) -204.492835490761 cp初值(数据文件 - 5) -2.4771057659164 cm初值(数据文件 - 5) -24.1373425222873 cp初值(数据文件 - 6) -7.81801022205547 cm初值(数据文件 - 6) -111.07846956739 cp初值(数据文件 - 7) 2.86865604413812 cm初值(数据文件 - 7) 101.125291957442 ====== 输出结果 ===== 文件：数据文件 - 1 No t 目标ctot 计算ctot 1 2 0.289 0.288999999981285 2 3 0.443 0.443000000008992 文件：数据文件 - 2 No t 目标ctot 计算ctot 1 2 0.38 0.37999999999826 2 3 0.417 0.417000000000818 文件：数据文件 - 3 No t 目标ctot 计算ctot 1 2 0.401 0.400999999995707 2 3 0.417 0.417000000002048 文件：数据文件 - 4 No t 目标ctot 计算ctot 1 2 0.368 0.367999999990831 2 3 0.538 0.538000000004393 文件：数据文件 - 5 No t 目标ctot 计算ctot 1 2 0.492 0.491999999996592 2 3 0.513 0.513000000001621 文件：数据文件 - 6 No t 目标ctot 计算ctot 1 2 0.183 0.182999999989207 2 3 0.202 0.202000000005179 文件：数据文件 - 7 No t 目标ctot 计算ctot 1 2 0.122 0.12200000000272 2 3 0.094 0.0939999999982719 |
6楼2012-04-14 20:29:52
v1iv1ia1n
至尊木虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 30721.4
- 散金: 18
- 红花: 14
- 帖子: 5617
- 在线: 557.7小时
- 虫号: 1193844
- 注册: 2011-01-21
- 专业: 化学环境污染与健康
2楼2012-04-13 21:37:20
3楼2012-04-14 09:59:06
4楼2012-04-14 10:05:00
v1iv1ia1n5楼2012-04-14 10:15:59
v1iv1ia1n
至尊木虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 30721.4
- 散金: 18
- 红花: 14
- 帖子: 5617
- 在线: 557.7小时
- 虫号: 1193844
- 注册: 2011-01-21
- 专业: 化学环境污染与健康
7楼2012-04-14 23:33:44
40