| 查看: 2457 | 回复: 8 | ||||
winterhao铁杆木虫 (正式写手)
|
[求助]
matlab常微分方程多参数拟合的参数之间cross correlations 求助@月只蓝以及matlab高手 已有1人参与
|
|
用lsqnonlin进行参数拟合后 怎样求参数 k1 k2 k3 之间的相关性 cross correlations lsqnonlin有这个功能吗 还是需要我用别的函数去评价 我看了您以前列出的例子 如果能结合任何一个简单的例子的code给我讲讲 那真是太感谢了 |
回复此楼» 猜你喜欢
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有247人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
-
3月 中国石油和化学工业联合会科学技术奖 马上开始，申获奖要点
已经有15人回复
-
中国石油和化学工业联合会科学技术奖的评审标准是什么？
已经有5人回复
-
经费如何使用，花不完怎么办？
已经有5人回复
-
兼职
已经有12人回复
-
海外Gustau Catalan课题组情况了解
已经有1人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68184.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.5小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2015-03-10 08:40:23
winterhao
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4492.7
- 散金: 20
- 红花: 2
- 帖子: 504
- 在线: 187.2小时
- 虫号: 666604
- 注册: 2008-12-02
- 性别: GG
- 专业: 流体力学
3楼2015-03-10 18:33:20
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68184.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.5小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2015-03-10 19:01:13
winterhao
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4492.7
- 散金: 20
- 红花: 2
- 帖子: 504
- 在线: 187.2小时
- 虫号: 666604
- 注册: 2008-12-02
- 性别: GG
- 专业: 流体力学
5楼2015-03-10 19:27:31
winterhao
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4492.7
- 散金: 20
- 红花: 2
- 帖子: 504
- 在线: 187.2小时
- 虫号: 666604
- 注册: 2008-12-02
- 性别: GG
- 专业: 流体力学
|
这是一个关于Mg(OH)2颗粒与水和CO2反应生成 MgCO3晶体的 多相流反应 包含了Mg(OH)2颗粒溶解模型 CO2吸收模型 以及MgCO3结晶模型 前两个模型参数已经确定 现在只需要用实验数据 拟合 MgCO3结晶模型的参数 我已经写好了 lsqnonlin的code 请见matlab文件 可以直接运行 运行发现： B = k1*(S^k3); % 成核速率 G = k2*(S^k4)；% 生长速率 第一个参数 k1 由于数量级很大（物理模型原因） 在拟合时变化很小 单独拟合这个参数变动也不大 拟合结果R2=0.90 我现在不确定是模型本身原因 还是我的初值选的不好的原因 我正在看您写的遗传算法 请您运行一下我的文件 帮我找找原因 非常感谢 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : fitting_lsqnonlin_Crystallization_kinetics.m
2015-03-10 20:11:18, 13.29 K
6楼2015-03-10 20:12:09
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68184.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.5小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
winterhao: 金币+50, ★★★很有帮助, 遗传算法我也试过了 R2能到0.93 我再试试1stopt 非常感谢 2015-03-10 22:30:32
winterhao: 金币+50, ★★★很有帮助, 遗传算法我也试过了 R2能到0.93 我再试试1stopt 非常感谢 2015-03-10 22:30:32
|
您好，m文件我看了，程序是可以运行的，计算结果不太好。 我觉得还是初值的问题，4个k初值之间的数量级相差太大，MATLAB在比较复杂的常微分方程拟合方面，初值的选取一直是一个问题。 如果遗传算法效果还是不理想，试试用1stopt软件，由于只有4个参数，你试试高版本的1stopt软件吧：http://muchong.com/bbs/viewthread.php?tid=8050245 |
7楼2015-03-10 22:01:33
winterhao
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4492.7
- 散金: 20
- 红花: 2
- 帖子: 504
- 在线: 187.2小时
- 虫号: 666604
- 注册: 2008-12-02
- 性别: GG
- 专业: 流体力学
8楼2015-03-11 04:15:46
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68184.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.5小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
9楼2015-03-11 08:32:18
