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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » 关于二元混合溶剂溶解度数据利用wilson、NRTL模拟的1stOPT代码
111/2返回列表12下一页
查看: 1739  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

嗨皮你我

新虫 (初入文坛)

[求助] 关于二元混合溶剂溶解度数据利用wilson、NRTL模拟的1stOPT代码 已有2人参与

各位大神，我测得了不同比例的某物质在二元混合溶剂中的溶解度随温度的变化，欲选用wilson、NRTL方程对其回归求得参数，下面是我的数据，但是1stOp软件t程序不大会写，总是提示function有错，请求高手指教
x2=0       
283.3        3.02194E-05
288.31        0.000157704
293.3        0.0004041
298.25        0.000686379
303.15        0.001026437
308.18        0.001405125
313.05        0.001810204
318.1        0.002726428
323.05        0.003769643
x2=0.05       
283.13        0.000349201
288.15        0.00046566
293.05        0.000981802
298.1        0.001500779
303.23        0.002474932
308.3        0.003373367
313.35        0.004337507
318.35        0.005850682
323.3        0.007716984
x2=0.1       
283.2        0.000828269
288.23        0.001597956
293.16        0.002383646
298.24        0.003548518
303.2        0.004900833
308.15        0.006296129
313.15        0.008064994
318.2        0.010216362
323.2        0.013075324
x2=0.15       
283.15        0.002970424
288.15        0.003905376
293.06        0.005280804
298.15        0.006616207
303.15        0.008243561
308.15        0.010303864
313.2        0.013119118
318.15        0.016480102
323.07        0.020616259
x2=0.2       
283.1        0.005204799
288.15        0.006279544
293.05        0.007921776
298.05        0.00986987
303.15        0.012566941
308.15        0.015421326
313.18        0.019498661
318.25        0.024651302
323.15        0.030872327
x2=0.25       
283.22        0.008048279
288.27        0.009556376
293.3        0.011981522
298.2        0.014520227
303.2        0.01789219
308.1        0.022389299
313.05        0.027783872
318.07        0.035092411
323.1        0.044652622
x2=0.3       
283.05        0.011149853
288.05        0.013651819
293.06        0.016540268
298.07        0.020332495
303.1        0.024782096
308.23        0.030544405
313.25        0.037978064
318.35        0.048406476
323.2        0.061434365
x2=0.35       
283.15        0.015106036
288.22        0.018196111
293.26        0.021838514
298.2        0.026852153
303.07        0.032794157
308.15        0.040592546
313        0.050666982
318.1        0.064596043
323.1        0.083519503
x2=0.4       
283.1        0.020213298
288.15        0.024960264
293.25        0.029919577
298.3        0.035835999
303.2        0.042800335
308.25        0.052897495
313.17        0.065920884
318.21        0.082267938
323.2        0.106310795
x2=0.45       
283.09        0.022163167
288.1        0.028017074
293.1        0.034843674
298.25        0.0420929
303.15        0.051962363
308.15        0.065005896
313.2        0.08215546
318.25        0.104250344
323.1        0.137575799
x2=0.5       
283.15        0.027972969
288.13        0.03399462
293.05        0.040605327
298.01        0.048112851
303.1        0.058477267
308.05        0.072571638
313.2        0.09047995
318.25        0.114083364
323.15        0.147697283

下面是我的程序
VarConstant x2=;
Constant R=8.314,V1=140.118,V2=97.9,V3=108.1053,H=23100,Tm=349.8;
Parameters g12,g21,g31,g13,g32,g23;
Variable T,x1;
constStr a21=V1/V2*exp(-g21/R/T),a12=V2/V1*exp(-g12/R/T),a31=V1/V3*exp(-g31/R/T),a13=V3/V1*exp(-g13/R/T),a23=V3/V2*exp(-g23/R/T),a32=V2/V3*exp(-g32/R/T);
constStr lnr=1-ln(x1+a12*x2+a13*(1-x1-x2))-x1/(x1+a12*x2+a13*(1-x1-x2))-a21*x2/(x2+a21*x1+a23(1-x1-x2))-a31*(1-x1-x2)/((1-x1-x2)+a31*x1+a32*x2);
constStr x1+x2+x3=1;
Function x1=exp(-H/R*(1/T-1/Tm)-lnr);
Data;
//T,x1
283.3        3.02194E-05
288.31        0.000157704
293.3        0.0004041
298.25        0.000686379
303.15        0.001026437
308.18        0.001405125
313.05        0.001810204
318.1        0.002726428
323.05        0.003769643
data;
283.13        0.000349201
288.15        0.00046566
293.05        0.000981802
298.1        0.001500779
303.23        0.002474932
308.3        0.003373367
313.35        0.004337507
318.35        0.005850682
323.3        0.007716984
data;
283.2        0.000828269
288.23        0.001597956
293.16        0.002383646
298.24        0.003548518
303.2        0.004900833
308.15        0.006296129
313.15        0.008064994
318.2        0.010216362
323.2        0.013075324
data;
283.15        0.002970424
288.15        0.003905376
293.06        0.005280804
298.15        0.006616207
303.15        0.008243561
308.15        0.010303864
313.2        0.013119118
318.15        0.016480102
323.07        0.020616259
data;
283.1        0.005204799
288.15        0.006279544
293.05        0.007921776
298.05        0.00986987
303.15        0.012566941
308.15        0.015421326
313.18        0.019498661
318.25        0.024651302
323.15        0.030872327
data;
283.22        0.008048279
288.27        0.009556376
293.3        0.011981522
298.2        0.014520227
303.2        0.01789219
308.1        0.022389299
313.05        0.027783872
318.07        0.035092411
323.1        0.044652622
data;
283.05        0.011149853
288.05        0.013651819
293.06        0.016540268
298.07        0.020332495
303.1        0.024782096
308.23        0.030544405
313.25        0.037978064
318.35        0.048406476
323.2        0.061434365
data;
283.15        0.015106036
288.22        0.018196111
293.26        0.021838514
298.2        0.026852153
303.07        0.032794157
308.15        0.040592546
313        0.050666982
318.1        0.064596043
323.1        0.083519503
data;
283.1        0.020213298
288.15        0.024960264
293.25        0.029919577
298.3        0.035835999
303.2        0.042800335
308.25        0.052897495
313.17        0.065920884
318.21        0.082267938
323.2        0.106310795
data;
283.09        0.022163167
288.1        0.028017074
293.1        0.034843674
298.25        0.0420929
303.15        0.051962363
308.15        0.065005896
313.2        0.08215546
318.25        0.104250344
323.1        0.137575799
data;
283.15        0.027972969
288.13        0.03399462
293.05        0.040605327
298.01        0.048112851
303.1        0.058477267
308.05        0.072571638
313.2        0.09047995
318.25        0.114083364
323.15        0.147697283

不知道哪里出现问题了，请求指教！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-07-11 16:14:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

嗨皮你我

新虫 (初入文坛)
自己顶一下，大神们请现身呐，急用！
赞一下 回复此楼
2楼2014-07-14 15:08:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
2楼: Originally posted by 嗨皮你我 at 2014-07-14 15:08:53
自己顶一下，大神们请现身呐，急用！

这方面的1stopt代码已经很多了吧，你去专家dingd的主页，看他的回帖。
想更深入的学习的话，小木虫里就有1stopt的帮助手册。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2014-07-14 15:33:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

x2, x3是什么？参数？起什么作用的？
赞一下 回复此楼
4楼2014-07-14 19:32:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

嗨皮你我

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by dingd at 2014-07-14 19:32:07
x2, x3是什么？参数？起什么作用的？

大神您好，跪求指导，是这样的：我做的是固体物质A在不同比例的两种混合溶剂（B和C）中的溶解度数据回归

x1是溶液中溶质的摩尔分率，x2、x3是溶液中两种溶剂B和C的摩尔分率，即有关系x1+x2+x3=1，其实我上边列出的那一大列数据中的x2应该改为xB，是溶剂B在混合溶剂B+C中所占比例，xB+xC=1，因为随着溶质A不断加入溶液中，x1、x2、x3都在变化，不是固定值。

比如这一组数据，
xB=0.05
    T                 xA
283.13        0.000349201
288.15        0.00046566
293.05        0.000981802
298.1        0.001500779
303.23        0.002474932
308.3        0.003373367
313.35        0.004337507
318.35        0.005850682
323.3        0.007716984
是指在混合溶剂中B摩尔分率为0.05时，A在混合溶剂B+C中摩尔溶解度x1随温度变化。

单对此组数据写程序还好，但是这样的话每个组成下都会拟合出6个参数，上面的10组数据就会出现10组参数，没法统一，我想整体拟合，最终只有6个参数，请求指导
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 溶解度模型.JPG
    • 2014-07-15 11:18:03, 46.48 K
5楼2014-07-15 11:15:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

嗨皮你我

新虫 (初入文坛)
引用回帖:
5楼: Originally posted by 嗨皮你我 at 2014-07-15 11:15:29
大神您好，跪求指导，是这样的：我做的是固体物质A在不同比例的两种混合溶剂（B和C）中的溶解度数据回归

x1是溶液中溶质的摩尔分率，x2、x3是溶液中两种溶剂B和C的摩尔分率，即有关系x1+x2+x3=1，其实我上边列出 ...

附件图片是一楼问题中所说的溶解度模型，我编的程序根本就不对，但是式子什么的应该都没错，请大神指教！！
赞一下 回复此楼
6楼2014-07-15 11:18:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

不懂你的专业，还是就1#代码而言，x3起什么作用，除了：x1+x2+x3=1外，其它任何地方都没有出现x3，x3是常数还是要求的参数？
赞一下 回复此楼
7楼2014-07-15 21:16:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

嗨皮你我

新虫 (初入文坛)
引用回帖:
7楼: Originally posted by dingd at 2014-07-15 21:16:47
不懂你的专业，还是就1#代码而言，x3起什么作用，除了：x1+x2+x3=1外，其它任何地方都没有出现x3，x3是常数还是要求的参数？

我把x3 用1-x1-x2代替了，您看附件图片中的关系式就行，那里面写得很清楚，要求的参数有g12,g21,g31,g13,g32,g23六个，而方程中T是自变量，x1是函数，知道的初值就是在每一个x2下：x1随T的变化
赞一下 回复此楼
8楼2014-07-16 14:14:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

嗨皮你我

新虫 (初入文坛)
大神现身呐，求助！
回复此楼
9楼2014-07-18 09:52:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

萧之夏

金虫 (小有名气)
你能给我解释下这是什么方程吗?
Function x1=exp(-H/R*(1/T-1/Tm)-lnr);
方程中H代表什么？
赞一下 回复此楼
10楼2016-03-29 20:20:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 嗨皮你我 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭