| 查看: 730 | 回复: 4 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
求助matlab高手看一下菜鸟写的非线性拟合程序，有关苯胺络合萃取的 已有2人参与
|
|||||
待拟合方程 ，需得到b1,b2的值原始数据： x1 x2 y 0.10 3.27 6.49 0.20 2.99 9.10 0.30 2.86 12.65 0.40 2.77 16.96 0.50 2.73 22.88 0.60 2.70 28.12 已写的程序，一直出错，求各位大神帮我修改一下： m文件： function y=aniline_1fun(b,x) x1=x(0:0.01:1.00); x2=x(4.00:-0.02:2.00); b1=m; b2=KE; y=b1*b2*0.0104^(b1-1)*x1/(1+b2*0.0104^b1)/(1+10^(4.67-x2)); >> x=[0.10,0.20,0.30,0.40,0.50,0.60; 3.27,2.99,2.86,2.77,2.73,2.70]'; >> y=[6.49,9.10,12.65,16.96,22.88,28.12]; >> b0=[0,5]'; >> [b,Q]=lsqcurvefit(aniline_1fun,b0,x,y) |
回复此楼» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有20人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有76人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有9人回复
-
桂林理工大学物理学专业招收调剂
已经有6人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68129.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Charles201: 金币+20, ★★★很有帮助 2015-08-05 15:14:38
感谢参与，应助指数 +1
Charles201: 金币+20, ★★★很有帮助 2015-08-05 15:14:38
|
dingd大神这么耐心地讲解，我直接给你结果吧。 我这些结果是用1stOpt软件做的，MATLAB肯定也可以做，在这里不我讨论MATLAB代码。 dingd大神说的过拟合很有道理的。 如果按照双参数公式： y=b1*b2*0.0104^(b1-1)*x1/(1+b2*0.0104^b1)/(1+10^(4.67-x2)) 将会得到无限组结果： b1 -0.480056251311531 b2 -0.110217887023378 ---------- ------------- b1 -1.43102339222656 b2 -0.0013970972484279 ---------- ------------- b1 -0.594292430579858 b2 -0.0652159450286256 ---------- ------------- b1 -0.438212028430683 b2 -0.13357689181283 ---------- ------------- b1 1.03961552899074 b2 -118.681614122118 …… …… 尽管b1和b2变化不定，但是两个参数按照下式计算： b3=b1*b2*0.0104^(b1-1)/(1+b2*0.0104^b1) b3一定都是等于3425.94298397334。 也就是说，楼主最原始的公式是过拟合的。 |
5楼2015-08-04 11:17:06
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Charles201: 金币+30, ★★★★★最佳答案 2015-08-05 15:14:26
感谢参与，应助指数 +1
Charles201: 金币+30, ★★★★★最佳答案 2015-08-05 15:14:26
|
公式有问题吧，过拟合： y=b1*b2*0.0104^(b1-1)*x1/(1+b2*0.0104^b1)/(1+10^(4.67-x2)); 可写成： y=b1*b2*0.0104^(b1-1)/(1+b2*0.0104^b1)/(1+10^(4.67-x2))*x1; 上式中：b1*b2*0.0104^(b1-1)/(1+b2*0.0104^b1)可以用一个参数代替，如b3 等价公式：y=b3/(1+10^(4.67-x2))*x1; 仅有一个参数 |
» 本帖已获得的红花（最新10朵）
Charles2012楼2015-08-03 18:22:37
3楼2015-08-03 18:26:06
|
我将m文件改成了如下形式 function y=aniline_1fun(b,x) x1=x(0:0.01:1.00); x2=x(4.00:-0.02:2.00); b3=b1*b2*0.0104^(b1-1)/(1+b2*0.0104^b1); y=b3*x1/(1+10^(4.67-x2)); 主程序： x=[0.10,0.20,0.30,0.40,0.50,0.60; 3.27,2.99,2.86,2.77,2.73,2.70]'; y=[6.49,9.10,12.65,16.96,22.88,28.12]; b3=[0]'; [b,Q]=lsqcurvefit(aniline_1fun,b0,x,y) 但是运行后提示输入参数的数目不足。请问这是什么原因呢？ |
4楼2015-08-03 19:39:23
