24小时热门版块排行榜    

查看: 1637  |  回复: 5

嬛嬛爱豆豆

木虫 (正式写手)

[求助] 求助高手:1stopt软件拟合二元体系的Wilson 方程,为什么拟合结果数值好小 已有1人参与

拟合结果好小,文献中都几百几千的,我想知道是为什么?
Title "Wilson Model";
Parameter a,b,c;
Variable T,x1;
Constant R=8.314;
Constant V1=58.29;
Constant V2=129.57;
ConstStr x2=1-x1;
ConstStr M=(V1/V2)exp(-b/(R*T));
ConstStr N=(V2/V1)exp(-c/(R*T));
Function a=exp(-ln(x1+M*x2)+(M/(x1+M*x2)-N/(x2+N*x1)));
Data;
T,x1=
293.15         0.1972
298.15         0.2008
303.15         0.2086
308.15         0.2067
313.15         0.2209
318.15         0.2210
323.15         0.2420
328.15         0.2496
333.15         0.2547
338.15         0.2835
拟合结果:
----------        -------------
a                  0.999999999999938
b                  1.29359453327508E-10
c                  -3.00010275167105E-10


求助高手:1stopt软件拟合二元体系的Wilson 方程,为什么拟合结果数值好小
20170323103623.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
代码不对。
Parameter关键词后面定义的是待求的参数;Variable关键词后面定义的是自变量和因变量;混淆致错。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2017-03-23 16:04:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

嬛嬛爱豆豆

木虫 (正式写手)

引用回帖:
2楼: Originally posted by 月只蓝 at 2017-03-23 16:04:34
代码不对。
Parameter关键词后面定义的是待求的参数;Variable关键词后面定义的是自变量和因变量;混淆致错。

区长大人:你好,感谢回复。
Parameter后面的a,b,c的确的要求的参数,Variable关键词后面也是自变量温度T,因变量x1。
现在只有T、x1的数据,希望得到b,c参数值。

附图中两个公式中,是不是我只编写了含有γ1的第一个公式来进行拟合而产生错误。
γ2也是未知的,若我写两个Function软件就会提示有错误无法继续计算运行。
不知道我描述的区长是否能明白,Help me
3楼2017-03-24 10:02:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 嬛嬛爱豆豆 at 2017-03-24 10:02:50
区长大人:你好,感谢回复。
Parameter后面的a,b,c的确的要求的参数,Variable关键词后面也是自变量温度T,因变量x1。
现在只有T、x1的数据,希望得到b,c参数值。

附图中两个公式中,是不是我只编写了含有γ1的 ...

两个方程共用参数,需要使用共享参数拟合模式。
普通1.5版本的1stOpt做不了,由于你的参数只有3个,可以用5.5试用版的来做。
相关的帖子,论坛里已经有很多,搜索了解之。

» 本帖已获得的红花(最新10朵)

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2017-03-24 10:06:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

嬛嬛爱豆豆

木虫 (正式写手)

送红花一朵
引用回帖:
4楼: Originally posted by 月只蓝 at 2017-03-24 10:06:57
两个方程共用参数,需要使用共享参数拟合模式。
普通1.5版本的1stOpt做不了,由于你的参数只有3个,可以用5.5试用版的来做。
相关的帖子,论坛里已经有很多,搜索了解之。...

谢谢,我去找找
5楼2017-03-24 10:35:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

付出+分享=?

新虫 (小有名气)

您好,请问您做出来了吗?

发自小木虫IOS客户端
6楼2023-03-10 23:02:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 嬛嬛爱豆豆 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 3/150 2026-11-10 01:40 by FChKX8meu1ds
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-10 01:29 by FChKX8meu1ds
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见