版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

zxw1982

木虫 (正式写手)

应助: 2 (幼儿园)
金币: 3396.8
帖子: 342
在线: 127.9小时
虫号: 260464
注册: 2006-06-17
性别: GG
专业: 能源化工

[求助] 求助高手：1stopt软件拟合三元体系的NRTL 方程，为什么一运行软件就死？已有1人参与

我用1stopt软件计算三元体系的NRTL 方程参数，不知道为什么每次一运行计算，软件就无响应。
求助高手给看看，是哪里出了问题？
非常感谢！

Title "NRTL Model";
Parameters g12,g13,g21,g23,g31,g32;
Variable T,x1,y;
Constant R=8.314;
Constant a=0.3;
ConstStr x2=0.1*(1-x1);
ConstStr x3=0.9*(1-x1);
ConstStr G12=exp(-a*t12);
ConstStr G13=exp(-a*t13);
ConstStr G21=exp(-a*t21);
ConstStr G23=exp(-a*t23);
ConstStr G31=exp(-a*t31);
ConstStr G32=exp(-a*t32);
ConstStr t12=g12/(R*T);
ConstStr t13=g13/(R*T);
ConstStr t21=g21/(R*T);
ConstStr t23=g23/(R*T);
ConstStr t31=g31/(R*T);
ConstStr t32=g32/(R*T);
Function y=exp((G21*x2+G31*x3)*(t21*G21*x2+t31*G31*x3)/(x1+x2*G21+x3*G31)^2+(t13*G13*x3*x3+G13*G23*x2*x3*(t13-t23))/(x3+x1*G13+x2*G23)^2+(t12*G12*x2*x2+G12*G32*x2*x3*(t12-t32))/(x2+x1*G12+x3*G32)^2);
Data;
//T,x1,y=
293.15 0.0418 3.8504
298.15 0.052 3.5225
303.15 0.0644 3.2231
308.15 0.0792 2.9579
313.15 0.0986 2.6711
318.15 0.1196 2.4667
323.15 0.1468 2.2432
328.15 0.1792 2.0443
333.15 0.2205 1.8423
338.15 0.2708 1.6584
343.15 0.3334 1.4849@月只蓝@beefly

回复此楼

» 猜你喜欢

山西大同大学物理学硕士研究生招收调剂生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有135人回复
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有2人回复
070200求调剂，一志愿某211，288分已经有18人回复
法国博士后职位已经有0人回复
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助厄米多项式已经有0人回复
NIS 配置方面问题已经有1人回复
FORTRAN程序DO WHILE出现语句错误，请哪位高手修改一下（紧急求助）已经有1人回复
【求助】关于vasp结构优化脚本运行输出的疑问已经有1人回复
数据挖掘决策树C4.5算法matlab程序参数意义求助已经有1人回复
晶体学报E审稿意见求助已经有1人回复
transiesta出错，求助。已经有0人回复
求助ubuntu下maui的安装配置已经有5人回复
gromacs做mdrun中出现Segmentation fault 已经有3人回复
【求助】读文本数据出现段错误[fortran] 已经有14人回复
【求助】fortran一程序解读已经有2人回复
【求助】有关packmol的问题已经有5人回复
【求助】程序出错已经有4人回复
【求助】请液相高手指教已经有5人回复

1楼 2016-04-18 10:40:44

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

回帖支持 ( 显示支持度最高的前 50 名 )

lipenggg

专家顾问 (正式写手)

专家经验: +65
应助: 122 (高中生)
金币: 4841.3
红花: 43
帖子: 722
在线: 435.7小时
虫号: 3929201
注册: 2015-06-17
性别: GG
专业: 流体力学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zxw1982: 金币+50, ★★★★★最佳答案, 高手，太牛了，非常感谢啊，小小金币，请笑纳 2016-04-18 13:24:32

Gxx和txx有点像循环，将Gxx带入公式就可以运行。
"NRTL Model"

迭代数: 57
计算用时(时:分:秒:微秒): 00:00:34:778
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.00562784682002697
残差平方和(SSE): 0.000348399258126565
相关系数(R): 0.999971252921243
相关系数之平方(R^2): 0.99994250666888
决定系数(DC): 0.999942482492335
卡方系数(Chi-Square): 6.69010417634671E-5
F统计(F-Statistic): 156530.894709368

参数最佳估算
---------- -------------
g12         3006.09966005237
g13         129.553350585897
g21         -3684.92515292263
g23         2.28772252829897
g31         5836.85272493941
g32         40922.3158718221

====== 结果输出 =====

No 实测值y 计算值y
1 3.8504 3.8504756
2 3.5225 3.5228717
3 3.2231 3.2225105
4 2.9579 2.9501120
5 2.6711 2.6865438
6 2.4667 2.4625225
7 2.2432 2.2424896
8 2.0443 2.0396803
9 1.8423 1.8418459
10 1.6584 1.6583433
11 1.4849 1.4879393

赞一下(2人)

回复此楼

» 本帖已获得的红花（最新10朵）

嬛嬛爱豆豆

一团和气

5楼2016-04-18 13:01:07

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

lipenggg

专家顾问 (正式写手)

专家经验: +65
应助: 122 (高中生)
金币: 4841.3
红花: 43
帖子: 722
在线: 435.7小时
虫号: 3929201
注册: 2015-06-17
性别: GG
专业: 流体力学
管辖: 计算模拟

【答案】应助回帖

感谢参与，应助指数 +1

Title "NRTL Model";
Parameters g12,g13,g21,g23,g31,g32;
Variable T,x1,y;
Constant R=8.314;
Constant a=0.3;
ConstStr x2=0.1*(1-x1);
ConstStr x3=0.9*(1-x1);
ConstStr t12=g12/(R*T);
ConstStr t13=g13/(R*T);
ConstStr t21=g21/(R*T);
ConstStr t23=g23/(R*T);
ConstStr t31=g31/(R*T);
ConstStr t32=g32/(R*T);
Function y=exp((exp(-a*t21)*x2+exp(-a*t31)*x3)*(t21*exp(-a*t21)*x2+t31*exp(-a*t31)*x3)/(x1+x2*exp(-a*t21)+x3*exp(-a*t31))^2
         +(t13*exp(-a*t13)*x3*x3+exp(-a*t13)*exp(-a*t23)*x2*x3*(t13-t23))
         /(x3+x1*exp(-a*t13)+x2*exp(-a*t23))^2+(t12*exp(-a*t12)*x2*x2+exp(-a*t12)*exp(-a*t32)*x2*x3*(t12-t32))/(x2+x1*exp(-a*t12)+x3*exp(-a*t32))^2);
Data;
//T,x1,y=
293.15       0.0418       3.8504
298.15       0.0520       3.5225
303.15       0.0644       3.2231
308.15       0.0792       2.9579
313.15       0.0986       2.6711
318.15       0.1196       2.4667
323.15       0.1468       2.2432
328.15       0.1792       2.0443
333.15       0.2205       1.8423
338.15       0.2708       1.6584
343.15       0.3334       1.4849

赞一下(1人)

回复此楼

一团和气

4楼2016-04-18 12:54:43

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68128.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

小写g12和大写G12是一个参数吗，还是不同的参数？
须注意1stOpt是不区分字母的大小写的。
如果是同一个参数，改成如下，即可运行：
Parameters g12,g13,g21,g23,g31,g32;
Variable T,x1,y;
Constant R=8.314;
Constant a=0.3;
ConstStr x2=0.1*(1-x1);
ConstStr x3=0.9*(1-x1);

ConstStr t12=g12/(R*T);
ConstStr t13=g13/(R*T);
ConstStr t21=g21/(R*T);
ConstStr t23=g23/(R*T);
ConstStr t31=g31/(R*T);
ConstStr t32=g32/(R*T);
Function y=exp((G21*x2+G31*x3)*(t21*G21*x2+t31*G31*x3)/(x1+x2*G21+x3*G31)^2+(t13*G13*x3*x3+G13*G23*x2*x3*(t13-t23))/(x3+x1*G13+x2*G23)^2+(t12*G12*x2*x2+G12*G32*x2*x3*(t12-t32))/(x2+x1*G12+x3*G32)^2);
Data;
//T,x1,y=
293.15       0.0418       3.8504
298.15       0.052       3.5225
303.15       0.0644       3.2231
308.15       0.0792       2.9579
313.15       0.0986       2.6711
318.15       0.1196       2.4667
323.15       0.1468       2.2432
328.15       0.1792       2.0443
333.15       0.2205       1.8423
338.15       0.2708       1.6584
343.15       0.3334       1.4849

此时计算结果：
相关系数(R): 0.999976628442786
相关系数之平方(R^2): 0.999953257431802
决定系数(DC): 0.999953257381529
卡方系数(Chi-Square): 5.37442244394238E-5
F统计(F-Statistic): 192534.977511446

参数最佳估算
---------- -------------
g12         -221394.823719256
g13         10692.2673284175
g21         -337.136660727573
g23         10533.6794712838
g31         -338.589449050954
g32         -253770.053103296

如果小写g12和大写G12是不同的参数，那按四楼那样改。

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

6楼2016-04-18 13:58:53

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68128.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

我在2楼中回复说“代码格式书写没有问题”，现在看来不是的，原代码在变量的大小写书写区分上存在问题的。

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

7楼2016-04-18 14:00:29

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68128.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

如变量不区分大小写，还有更好的结果，决定系数可达5个“9”：
相关系数(R): 0.999996633225546
相关系数之平方(R^2): 0.999993266462426
决定系数(DC): 0.999993266385066
卡方系数(Chi-Square): 8.55335536774635E-6
F统计(F-Statistic): 1336584.12088942

参数最佳估算
---------- -------------
g12         776.942748582934
g13         -415.294727118021
g21         -0.118066444579496
g23         -2966.48681543665
g31         -0.118070641087448
g32         47.1205980648154

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

8楼2016-04-18 14:04:51

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

lipenggg

专家顾问 (正式写手)

专家经验: +65
应助: 122 (高中生)
金币: 4841.3
红花: 43
帖子: 722
在线: 435.7小时
虫号: 3929201
注册: 2015-06-17
性别: GG
专业: 流体力学
管辖: 计算模拟

内容已删除

赞一下(1人)

回复此楼

一团和气

9楼2016-04-18 16:52:48

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

普通回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68128.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

代码格式、公式书写应该都没问题，应该是1stOpt版本问题了。

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2016-04-18 11:53:49

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

zxw1982

木虫 (正式写手)

应助: 2 (幼儿园)
金币: 3396.8
帖子: 342
在线: 127.9小时
虫号: 260464
注册: 2006-06-17
性别: GG
专业: 能源化工

引用回帖:

2楼: Originally posted by 月只蓝 at 2016-04-18 11:53:49
代码格式、公式书写应该都没问题，应该是1stOpt版本问题了。

我怎么检查也搞不明白是什么问题，换成Wilson方程，也是拟合六个参数，就能运行。很奇怪

赞一下

回复此楼

3楼2016-04-18 12:12:13

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

zxw1982

木虫 (正式写手)

应助: 2 (幼儿园)
金币: 3396.8
帖子: 342
在线: 127.9小时
虫号: 260464
注册: 2006-06-17
性别: GG
专业: 能源化工

引用回帖:

6楼: Originally posted by 月只蓝 at 2016-04-18 13:58:53
小写g12和大写G12是一个参数吗，还是不同的参数？
须注意1stOpt是不区分字母的大小写的。
如果是同一个参数，改成如下，即可运行：
Parameters g12,g13,g21,g23,g31,g32;
Variable T,x1,y;
Constant R=8.314;
...

g12和G12在公式中是有不同含义的。我不知道这个软件不分大小写。非常感谢您指出这个问题

赞一下

回复此楼

10楼2016-04-18 18:20:23

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

相关版块跳转我要订阅楼主 zxw1982 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 270分求调剂 +5	maxjxbsk 2026-04-01	5/250	2026-04-06 15:32 by lin-da
[考研] 求调剂 +4	wos666 2026-04-03	5/250	2026-04-06 15:22 by wos666
[考研] 0857大类环境工程B区求调剂 +3	龚禹铭 2026-04-05	3/150	2026-04-06 10:22 by 蓝云思雨
[考研] 086000生物与医药298调剂求助 +9	元元青青 2026-03-31	12/600	2026-04-05 21:03 by 学员8dgXkO
[考研] 求调剂 +7	张.1 2026-04-05	7/350	2026-04-05 20:40 by 啵啵啵0119
[考研] 0854电子信息319求调剂（接受跨专业调剂） +3	星星不眨眼喽 2026-04-05	3/150	2026-04-05 20:20 by 啵啵啵0119
[考研] 调剂求助 +10	想换手机不想解� 2026-04-02	13/650	2026-04-05 09:41 by sam3303
[考研] 一志愿北京2，材料与化工308求调剂 +10	熊二想上岸 2026-04-04	10/500	2026-04-05 05:20 by houyaoxu
[考研] 一志愿北京化工大学，初试成绩350求调剂 +9	沿岸?贝壳 2026-04-04	14/700	2026-04-05 01:09 by 沿岸?贝壳
[考研] 283分求调剂 +7	小聂爱学习 2026-04-03	7/350	2026-04-04 21:51 by hemengdong
[考研] 283求调剂 +4	mcbbc 2026-04-03	5/250	2026-04-04 20:51 by imissbao
[考研] 292求调剂 +11	2022080213 2026-04-04	13/650	2026-04-04 18:38 by macy2011
[考研] 280求调剂 +21	咕噜晓晓 2026-04-02	22/1100	2026-04-04 11:12 by 猪会飞
[考研] 321求调剂 +17	y-yh 2026-04-01	20/1000	2026-04-03 12:57 by y-yh
[考研] 求调剂22408 288分 +5	new382 2026-04-02	5/250	2026-04-03 09:13 by 醉在风里
[考研] 一志愿陕西师范大学生物学317分 +5	1563日。 2026-04-02	5/250	2026-04-03 06:58 by ilovexiaobin
[考研] 材料化工340求调剂 +5	jhx777 2026-03-30	5/250	2026-04-02 12:45 by smileboy2006
[考研] 11408 321分求调剂 +3	huchun12138 2026-03-30	4/200	2026-04-01 22:48 by guanxin1001
[考研] 材料调剂 +10	Eujd1 2026-03-31	11/550	2026-04-01 11:23 by ivanqyq
[考研] 340求调剂 +4	希望如此i 2026-03-31	4/200	2026-03-31 16:40 by 690616278

24小时热门版块排行榜

[求助] 求助高手：1stopt软件拟合三元体系的NRTL 方程，为什么一运行软件就死？ 已有1人参与

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

【答案】应助回帖

» 本帖已获得的红花（最新10朵）

【答案】应助回帖

[求助] 求助高手：1stopt软件拟合三元体系的NRTL 方程，为什么一运行软件就死？已有1人参与