24小时热门版块排行榜    

查看: 2395  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

清晨清晨

银虫 (初入文坛)

[求助] 好心虫友,帮忙看看这个matlab程序那里出错了 已有1人参与

用matlab拟合数据,进行NRTL模型拟合:我编写了如下的程序,一直出错,麻烦有懂的虫友帮忙看看,怎么改错??
(一)定义.m文件。
function g4=f(beta0,x1,x2,T,R)
g1=beta0(1);
g2=beta0(2);
g3=beta0(3);
g4=beta0(4);
g4=x2.*x2.*((((g(1)./(R.*T)).*(exp(-g(3).*g(1)./(R.*T))).^2)./((x1+(exp(-g(3).*g(1)./(R.*T))).*x2).^2))+(((g(2)./(R.*T)).*(exp(-g(3).*g(2)./(R.*T))).^2)./((x2+(exp(-g(3).*g(2)./(R.*T))).*x1).^2)))
(二)编写程序。
x1=[0.01525        0.02059        0.02685        0.03417        0.04221        0.04958];
x2=[0.98475        0.97941        0.97315        0.96583        0.95779        0.95042];
T=[298.15        303.15        308.15        313.15        318.15        323.15];
R=[8.314 8.314 8.314 8.314 8.314 8.314];
beta0=[0,0,0,0];
[g,resnorm]=lsqcurvefit(@f,beta0,x1,x2,T,R)
运行后出错:
Warning: Length of lower bounds is > length(x); ignoring extra bounds.
> In checkbounds at 27
  In optim\private\lsqncommon at 39
  In lsqcurvefit at 149
  In Untitled2 at 6
Warning: Length of upper bounds is > length(x); ignoring extra bounds.
> In checkbounds at 37
  In optim\private\lsqncommon at 39
  In lsqcurvefit at 149
  In Untitled2 at 6
Exiting due to infeasibility:  4 lower bounds exceed the corresponding upper bounds.

g =

     0     0     0     0


resnorm =

     []
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清晨清晨

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by somomo91 at 2017-07-09 17:29:24
这里的 beta0 有4个元素,但是 上下边界 有 6个元素,所以会出错
beta0 元素的个数应该和 R,T 一致。
lsqcurvefit(fun,x0,xdata,ydata,lb,ub)

但即使beta0 有6个0,你会看到结果应该还不对。
没看方程,没法 ...

这个是需要拟合的公式,只需要求4个变量,所以我beta0只给了4个值
好心虫友,帮忙看看这个matlab程序那里出错了
IMG_20170709_174704.jpg

3楼2017-07-09 17:49:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

somomo91

专家顾问 (职业作家)

【答案】应助回帖

★
感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2017-09-28 23:50:04
这里的 beta0 有4个元素,但是 上下边界 有 6个元素,所以会出错
beta0 元素的个数应该和 R,T 一致。
CODE:
lsqcurvefit(fun,x0,xdata,ydata,lb,ub)

但即使beta0 有6个0,你会看到结果应该还不对。
没看方程,没法确定到底你的程序哪里出错
2楼2017-07-09 17:29:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清晨清晨

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by somomo91 at 2017-07-09 17:29:24
这里的 beta0 有4个元素,但是 上下边界 有 6个元素,所以会出错
beta0 元素的个数应该和 R,T 一致。
lsqcurvefit(fun,x0,xdata,ydata,lb,ub)

但即使beta0 有6个0,你会看到结果应该还不对。
没看方程,没法 ...

这个是需要拟合的公式,只需要求4个变量,所以我beta0只给了4个值
好心虫友,帮忙看看这个matlab程序那里出错了-1
IMG_20170709_174704.jpg

4楼2017-07-09 17:49:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

somomo91

专家顾问 (职业作家)

【答案】应助回帖

★
jjdg: 金币+1, 感谢参与 2017-09-28 23:50:16
方程里面有 ln(gamma1),所以你的程序,正确的应该是 g4 = exp(x2.*x2.* ......)
其次,tau12 和 tau21 的值应该不是零吧,但你的程序里面 设定 tau12 = tau21 = 0
这也是为什么 即使没有边界错误,输出为零的原因
5楼2017-07-09 18:04:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 申请科研项目纯靠实力能成功吗? 200+7 David1982 2026-10-13 11/550 2026-10-15 10:30 by zyqchem
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见