24小时热门版块排行榜    

查看: 1075  |  回复: 3

小东东喵

新虫 (初入文坛)

[求助] 求助电脑大神 软件拟合1stopt二元体系的Wilson方程 为什么总是得不到拟合的结果 已有1人参与

求助电脑大神。我是用了1.5版,发现做不了,哪位大神可以指点一下   用5.5版共享参数拟合模式要怎么做?或者帮我运行一下。因为要毕业了,比较着急,谢谢啦、
itle "Wilson Model";                                       
Parameters a,b,c,d;                                       
Variable T,x1;                                       
Constant R=8.314;                                       
ConstStr x2=1-x1;                                       
ConstStr M=exp*(-c)/R*T;                                       
ConstStr N=exp*(-d)/R*T;                                       
Function a=exp(-ln(x1+M*x2)+x2*(M/(x1+M*x2)-N/(x2+N*x1)));                                       
Function b=exp(-ln(x2+N*x1)-x1*(M/(x1+M*x2)-N/(x2+N*x1)));                                       
Data;                                       
T,x1=                                       
305.15        0                               
301.15        0.137                               
299.15        0.186                               
296.8        0.266                               
295        0.302                               
286.65        0.481                               
286.65        0.481                               
295.15        0.571                               
299.15        0.608                               
303.15        0.676                               
307.9        0.781                               
313.15        0.9                               
316.15        1@月只蓝@beefly
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
a=exp(-ln(x1+M*x2)+x2*(M/(x1+M*x2)-N/(x2+N*x1)));                                       
b=exp(-ln(x2+N*x1)-x1*(M/(x1+M*x2)-N/(x2+N*x1)));         

将这两行转化一下,都写x1=f(x1,..)形式试试
数值计算
2楼2018-12-06 18:37:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wonderful899

木虫 (小有名气)

引用回帖:
2楼: Originally posted by 独孤神宇 at 2018-12-06 18:37:09
a=exp(-ln(x1+M*x2)+x2*(M/(x1+M*x2)-N/(x2+N*x1)));                                       
b=exp(-ln(x2+N*x1)-x1*(M/(x1+M*x2)-N/(x2+N*x1)));         

将这两行转化一下,都写x1=f(x1,..)形式试试

我也碰到楼主一样的问题,关键转化成x1=f(x1,..)形式,写不出来啊
3楼2022-08-09 22:13:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by wonderful899 at 2022-08-09 22:13:25
我也碰到楼主一样的问题,关键转化成x1=f(x1,..)形式,写不出来啊...

怎么会写不出来,不需要显示表达,隐式也可以。

发自小木虫Android客户端
数值计算
4楼2022-08-10 14:08:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小东东喵 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-02 02:58 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 20/1000 2026-11-02 02:52 by xzcv0Tehd8lQ
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-02 02:45 by xzcv0Tehd8lQ
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 X5bHxOULxFqb 2026-10-30 24/1200 2026-11-02 02:39 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 3/150 2026-11-02 01:00 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 4/200 2026-11-02 00:50 by xzcv0Tehd8lQ
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:24 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 21:45 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 16:58 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
信息提示
请填处理意见