| 5 | 1/1 | 返回列表 |
| 查看: 1812 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
關(guān)于二元混合溶劑溶解度數(shù)據(jù)利用wilson、NRTL模擬的1stOPT代碼 已有2人參與
|
|||
|
各位大神,我測得了不同比例的某物質(zhì)在二元混合溶劑中的溶解度隨溫度的變化,欲選用wilson、NRTL方程對其回歸求得參數(shù),下面是我的數(shù)據(jù),但是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 不知道哪里出現(xiàn)問題了,請求指教! |
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

鐵桿木蟲 (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +5 | @taotao 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿北化求調(diào)劑 +4 | Jsman 2026-03-22 | 4/200 |
|
|
[考研] 375求調(diào)劑 +3 | 雨夏整夜 2026-03-29 | 3/150 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +11 | 北道巷 2026-03-24 | 11/550 |
|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 打過很多競賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑。。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|