应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助matlab高手看一下菜鸟写的非线性拟合程序，有关苯胺络合萃取的
51/1返回列表
查看: 690  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Charles201

新虫 (小有名气)

[求助] 求助matlab高手看一下菜鸟写的非线性拟合程序，有关苯胺络合萃取的 已有2人参与

待拟合方程求助matlab高手看一下菜鸟写的非线性拟合程序，有关苯胺络合萃取的，需得到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)
回复此楼

» 猜你喜欢
1楼 2015-08-03 18:13:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Charles201

新虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2015-08-03 18:22:37
公式有问题吧，过拟合：

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.010 ...

我将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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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朵）

Charles201
2楼2015-08-03 18:22:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Charles201

新虫 (小有名气)
送红花一朵
引用回帖:
2楼: Originally posted by dingd at 2015-08-03 18:22:37
公式有问题吧，过拟合：

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.010 ...

谢谢，我一会改一下，试试能不能搞定

[ 发自小木虫客户端 ]
赞一下 回复此楼
3楼2015-08-03 18:26:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Charles201: 金币+20, ★★★很有帮助 2015-08-05 15:14:38
引用回帖:
4楼: Originally posted by Charles201 at 2015-08-03 19:39:23
我将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=';
y=;
b3=';
=l ...

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。
也就是说，楼主最原始的公式是过拟合的。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2015-08-04 11:17:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭