| 查看: 23796 | 回复: 73 | |||||||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||||||||||
dshuanger木虫 (小有名气)
|
[求助]
matlab非线性拟合模型参数，大神你在哪！
|
||||||||||||||
|
本人非计算模拟专业，小白一枚，但最近需要模拟参数，不知如何运用matlab求参数，顾求助各位大神。 模型表达式： y=a/E1 +a/E2[1-exp(-x/t)]+a*x/E3 其中：x=[2，5，10，15，30，60，120，200，300，550，700]； y=[8.66667，10，10.66667，12.66667，13.33333，15，15.33333， 16.66667，17.33333，18.33333，20]； a=2 未知参数：E1，E2，t，E3 不知道能否通过matlab拟合参数 请各位大神能鼎力相助，最好能给出代码或者计算方法，谢谢大家！ [ Last edited by dshuanger on 2013-11-5 at 11:06 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 matlab典型案例及小技巧 | 宝贵经验 | Polyflow反应模拟 | 分子动力学 |
| 程序 | 第一性原理杂 | 研究僧 | 生物试验统计分析 |
| 动力学拟合 | 技术贴 | 考博 | 科学计算tool疑难问题 |
| 软件 |
» 猜你喜欢
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有236人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 非线性拟合得到的模型可靠性怎么验证？
已经有8人回复
- 用Matlab做多变量非线性拟合求参数值，求各位大侠鼎力相助！
已经有3人回复
- 非线性神经网络模型的模型预测控制 程序问题
已经有4人回复
- origin非线性拟合
已经有5人回复
- 关于matlab微分、及非线性拟合问题
已经有5人回复
- 用spss拟合生长曲线（非线性回归）其中参数初始值怎么确定？
已经有10人回复
- Herschel-Bulkley模型怎么拟合
已经有3人回复
- 如何用origin进行非线性拟合？
已经有9人回复
- 请问matlab拟合二元方程的问题
已经有3人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- 求助如何用matlab取出直线上的点
已经有10人回复
- 是否可以用求模型参数的数据来验证求模型参数方法的可靠性？（van Genuchten）
已经有24人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- 用matlab拟合双因变量实验结果求助
已经有7人回复
- matlab非线性参数拟合问题
已经有7人回复
- 拟合模型的好坏
已经有4人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
- 【求助】急求关于matlab做非线性拟合的函数
已经有8人回复
- 【求助】求助origin数学模型
已经有6人回复
- 【求助】多变量非线性拟合【已完结】
已经有8人回复
- 求助： 1stopt 进行非线性拟合
已经有5人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
dshuanger: 金币+100, ★★★很有帮助, 非常感谢，呵呵，虽然有些看不懂，但感觉比较相符 2013-11-05 20:21:50
感谢参与，应助指数 +1
dshuanger: 金币+100, ★★★很有帮助, 非常感谢，呵呵，虽然有些看不懂，但感觉比较相符 2013-11-05 20:21:50
|
结果不唯一。另外一组结果： e1 = 54.5440439848079 e2 = 241.457452201238 t = 0.104301409896579 e3 = 2.35332457537677E17 ====== 输出结果 ===== No. Observed y Calculated y 1 8.6667 8.66670000000005 2 10 9.9999999999611 3 10.6667 10.6666999999611 4 12.6667 12.6666999999611 5 13.3333 13.3332999999611 6 15 14.9999999999611 7 15.3333 15.3332999999611 8 16.6667 16.6666999999611 9 17.3333 17.3332999999611 10 18.3333 18.3332999999611 11 20 19.9999999999611 |
9楼2013-11-05 14:42:42
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 计算强帖: 1
- 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2013-11-05 11:23:25
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
★
dingd: 金币+1, 谢谢分享。 2013-11-05 16:37:58
dingd: 金币+1, 谢谢分享。 2013-11-05 16:37:58
|
y=2/e1 +2/e2/(1-exp(-x/t))+a*x/3; // a=2 x=[2，5，10，15，30，60，120，200，300，550，700]； y=[8.66667，10，10.66667，12.66667，13.33333，15，15.33333， 16.66667，17.33333，18.33333，20]； 相关系数(R): 1 相关系数之平方(R^2): 1 决定系数(DC): 0.999999999996653 卡方系数(Chi-Square): 1.64243431014082E-11 参数 最佳估算 ---------- -------------- e1 =1.2090897803241 e2 =3.11591327892529 t =0.00209072546506921 ====== 输出结果 ===== No. Observed y Calculated y 1 8.66667 8.66666365993411 2 10 9.99999365993411 3 10.66667 10.6666636599341 4 12.66667 12.6666636599341 5 13.33333 13.3333236599341 6 15 14.9999936599341 7 15.33333 15.3333236599341 8 16.66667 16.6666636599341 9 17.33333 17.3333236599341 10 18.33333 18.3333236599341 11 20 19.9999936599341 |
3楼2013-11-05 11:54:42
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
4楼2013-11-05 11:56:53
20