24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

嗨皮你我

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 438.5
帖子: 41
在线: 16.8小时
虫号: 1975097
注册: 2012-09-04
专业: 传递过程

[求助] 关于二元混合溶剂溶解度数据利用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

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

回复此楼