| 查看: 23958 | 回复: 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 ] |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有208人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有11人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有1人回复
|
Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9,X1[1,7],X2[1,7],X3[1,7]; Function Y=b0+b1*X1+b2*X2+b3*X3+b4*X1*X2+b5*X1*X3+b6*X2*X3+b7*X1^2+b8*X2^2+b9*X3^2; Minimum Z Z=(b0+b1*1+b2*2+b3*4+b4*1*2+b5*1*4+b6*2*4+b7*1^2+b8*2^2+b9*4^2-76.84)^2 +(b0+b1*1+b2*4+b3*7+b4*1*4+b5*1*7+b6*4*7+b7*1^4+b8*4^2+b9*7^2-77.85)^2 +(b0+b1*2+b2*6+b3*3+b4*2*6+b5*2*3+b6*6*3+b7*2^4+b8*6^2+b9*3^2-62.19)^2 +(b0+b1*2+b2*1+b3*7+b4*2*1+b5*2*7+b6*1*7+b7*2^4+b8*1^2+b9*7^2-70.77)^2+(b0+b1*3+b2*3+b3*3+b4*3*3+b5*3*3+b6*3*3+b7*3^3+b8*3^2+b9*3^2-73.30)^2 +(b0+b1*3+b2*5+b3*6+b4*3*5+b5*3*6+b6*5*6+b7*3^4+b8*5^2+b9*6^2-66.23)^2 +(b0+b1*4+b2*7+b3*2+b4*4*7+b5*4*2+b6*7*2+b7*4^4+b8*7^2+b9*2^2-10.66)^2 +(b0+b1*4+b2*1+b3*6+b4*4*1+b5*4*6+b6*1*6+b7*4^4+b8*1^2+b9*6^2-79.36)^2 +(b0+b1*5+b2*3+b3*2+b4*5*3+b5*5*2+b6*3*2+b7*5^4+b8*3^2+b9*2^2-79.87)^2 +(b0+b1*5+b2*5+b3*5+b4*5*5+b5*5*5+b6*5*5+b7*5^4+b8*5^2+b9*5^2-77.34)^2 +(b0+b1*6+b2*7+b3*1+b4*6*7+b5*6*1+b6*7*1+b7*6^4+b8*7^2+b9*1^2-21.77)^2 +(b0+b1*6+b2*2+b3*5+b4*6*2+b5*6*5+b6*2*5+b7*6^4+b8*2^2+b9*5^2-74.82)^2 +(b0+b1*7+b2*4+b3*2+b4*7*4+b5*7*1+b6*4*1+b7*7^4+b8*4^2+b9*1^2-74.31)^2 +(b0+b1*7+b2*6+b3*4+b4*7*6+b5*7*4+b6*6*4+b7*7^4+b8*6^2+b9*4^2-58.65)^2 Data; ... 然后就出现定义Z有误了... 能不能麻烦你给出可执行的命令啊？谢谢  16.png |
55楼2015-01-23 15:31:21
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
