| 查看: 1832 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
关于二元混合溶剂溶解度数据利用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 不知道哪里出现问题了，请求指教！ |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有95人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有11人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 自己这一段时间过柱子的总结
已经有11人回复
9楼2014-07-18 09:52:45
2楼2014-07-14 15:08:53
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68129.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼2014-07-14 15:33:47
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2014-07-14 19:32:07
