| 查看: 1339 | 回复: 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
» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
新西兰Robinson研究所招收全奖PhD
已经有0人回复
物理学I论文润色/翻译怎么收费?
已经有298人回复
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
笼目材料中量子自旋液体基态的证据
已经有0人回复
数学教学论硕士可以读数学物理博士吗?
已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
求助NH4V4O10晶体的CIF文件
已经有0人回复
» 本主题相关价值贴推荐,对您同样有帮助:
matlab拟合方程参数时初值的选择
已经有15人回复
matlab拟合求值
已经有17人回复
在matlab中拟合曲面的表达式
已经有29人回复
matlab拟合拟均相动力学参数
已经有6人回复
需要用matlab拟合数据的外行人急需帮助
已经有10人回复
matlab非线性参数拟合问题
已经有7人回复
【求助】用matlab最优化方法进行参数拟合
已经有17人回复
【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
如何用MATLAB拟合指数函数
已经有3人回复
【求助】matlab曲面拟合表达式
已经有10人回复
【求助】如何使得matlab拟合的曲线强制经过一个点
已经有8人回复
【求助】使用MATLAB怎么实现拟合力场参数的程序?【已解决】
已经有15人回复
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
5楼2012-04-14 10:15:59
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
- 专业: 化学环境污染与健康
7楼2012-04-14 23:33:44












回复此楼
v1iv1ia1n