| 查看: 24770 | 回复: 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 ] |
回复此楼» 猜你喜欢
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有218人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
|
我自己尝试了下单方程最优化：Parameters X1[1,7], X2[1,7], X3[1,7]; Variable X1,X2,X3,Y; Minimum = Y; Function Y=355.45-38.35*X1-13.99*X2-85.78*X3+0.825*X1*X2+6.234*X1*X3+1.959*X2*X3+1.278*X1^2+0.272*X2^2+6.59*X3^2; Data; 1 2 4 96.15 1 4 7 89.07 2 6 3 95.44 2 1 7 96.08 3 3 3 95.42 3 5 6 96.37 4 7 2 92.66 4 1 6 96.49 5 3 2 95.73 5 5 5 93.93 6 7 1 92.53 6 2 5 95.68 7 4 1 95.37 7 6 4 93.97 得到X1 X2 X3最佳估算 没有Y值 请问这样准确吗？ |
36楼2015-01-21 22:00:59
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 计算强帖: 1
- 应助: 454 (硕士)
- 金币: 20659.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.7小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2013-11-05 11:23:25
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3891
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24217
- 在线: 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
- 金币: 3891
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24217
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
4楼2013-11-05 11:56:53
100