| 查看: 23956 | 回复: 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人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有186人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有11人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68129.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
|
运行代码： 结果： 相关系数之平方(R^2): 0.910015554060544 决定系数(DC): 0.910015554060543 卡方系数(Chi-Square): 6.30625010587871 F统计(F-Statistic): 121.356380369046 参数 最佳估算 ---------- ------------- b0 274.297442101538 b1 -44.818446582328 b2 21.7155760160432 b3 -78.7219393455877 b4 0.820652156345272 b5 6.79527940302928 b6 0.732950301205965 b7 1.89990679569777 b8 -4.33596272234971 b9 6.34139744840854 ====== 结果输出 ===== No 实测值y 计算值y 1 76.84 78.7268296 2 77.85 79.9128753 3 62.19 51.1783409 4 70.77 71.2205570 5 73.3 79.1110228 6 66.23 69.6920682 7 10.66 20.4925341 8 79.36 69.5274660 9 79.87 76.4079319 10 77.34 71.5289771 11 21.77 21.3194432 12 74.82 85.8316589 13 74.31 72.2471250 14 58.65 56.7631700 可以看到上面的三列数据，实测值y和计算值y 都在[0,100]的范围。 |
48楼2015-01-23 09:16:14
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
- 金币: 3880
- 散金: 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
- 金币: 3880
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
4楼2013-11-05 11:56:53
