版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

嗨皮你我

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 438.5
帖子: 41
在線: 16.8小時(shí)
蟲號(hào): 1975097
注冊(cè): 2012-09-04
專業(yè): 傳遞過程

[求助] 關(guān)于二元混合溶劑溶解度數(shù)據(jù)利用wilson、NRTL模擬的1stOPT代碼已有2人參與

各位大神，我測得了不同比例的某物質(zhì)在二元混合溶劑中的溶解度隨溫度的變化，欲選用wilson、NRTL方程對(duì)其回歸求得參數(shù)，下面是我的數(shù)據(jù)，但是1stOp軟件t程序不大會(huì)寫，總是提示function有錯(cuò)，請(qǐng)求高手指教
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ǐng)求指教！

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有84人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有29人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有53人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

自己這一段時(shí)間過柱子的總結(jié) 已經(jīng)有11人回復(fù)

1樓 2014-07-11 16:14:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

嗨皮你我

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 438.5
帖子: 41
在線: 16.8小時(shí)
蟲號(hào): 1975097
注冊(cè): 2012-09-04
專業(yè): 傳遞過程

自己頂一下，大神們請(qǐng)現(xiàn)身吶，急用！

贊一下

回復(fù)此樓

2樓2014-07-14 15:08:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

引用回帖:

2樓: Originally posted by 嗨皮你我 at 2014-07-14 15:08:53
自己頂一下，大神們請(qǐng)現(xiàn)身吶，急用！

這方面的1stopt代碼已經(jīng)很多了吧，你去專家dingd的主頁，看他的回帖。
想更深入的學(xué)習(xí)的話，小木蟲里就有1stopt的幫助手冊(cè)。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

3樓2014-07-14 15:33:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

x2, x3是什么？參數(shù)？起什么作用的？

贊一下

回復(fù)此樓

4樓2014-07-14 19:32:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

嗨皮你我

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 438.5
帖子: 41
在線: 16.8小時(shí)
蟲號(hào): 1975097
注冊(cè): 2012-09-04
專業(yè): 傳遞過程

引用回帖:

4樓: Originally posted by dingd at 2014-07-14 19:32:07
x2, x3是什么？參數(shù)？起什么作用的？

大神您好，跪求指導(dǎo)，是這樣的：我做的是固體物質(zhì)A在不同比例的兩種混合溶劑（B和C）中的溶解度數(shù)據(jù)回歸

x1是溶液中溶質(zhì)的摩爾分率，x2、x3是溶液中兩種溶劑B和C的摩爾分率，即有關(guān)系x1+x2+x3=1，其實(shí)我上邊列出的那一大列數(shù)據(jù)中的x2應(yīng)該改為xB，是溶劑B在混合溶劑B+C中所占比例，xB+xC=1，因?yàn)殡S著溶質(zhì)A不斷加入溶液中，x1、x2、x3都在變化，不是固定值。

比如這一組數(shù)據(jù)，
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時(shí)，A在混合溶劑B+C中摩爾溶解度x1隨溫度變化。

單對(duì)此組數(shù)據(jù)寫程序還好，但是這樣的話每個(gè)組成下都會(huì)擬合出6個(gè)參數(shù)，上面的10組數(shù)據(jù)就會(huì)出現(xiàn)10組參數(shù)，沒法統(tǒng)一，我想整體擬合，最終只有6個(gè)參數(shù)，請(qǐng)求指導(dǎo)

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 溶解度模型.JPG

2014-07-15 11:18:03, 46.48 K

5樓2014-07-15 11:15:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

嗨皮你我

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 438.5
帖子: 41
在線: 16.8小時(shí)
蟲號(hào): 1975097
注冊(cè): 2012-09-04
專業(yè): 傳遞過程

引用回帖:

5樓: Originally posted by 嗨皮你我 at 2014-07-15 11:15:29
大神您好，跪求指導(dǎo)，是這樣的：我做的是固體物質(zhì)A在不同比例的兩種混合溶劑（B和C）中的溶解度數(shù)據(jù)回歸

x1是溶液中溶質(zhì)的摩爾分率，x2、x3是溶液中兩種溶劑B和C的摩爾分率，即有關(guān)系x1+x2+x3=1，其實(shí)我上邊列出 ...

附件圖片是一樓問題中所說的溶解度模型，我編的程序根本就不對(duì)，但是式子什么的應(yīng)該都沒錯(cuò)，請(qǐng)大神指教��！

贊一下

回復(fù)此樓

6樓2014-07-15 11:18:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

不懂你的專業(yè)，還是就1#代碼而言，x3起什么作用，除了：x1+x2+x3=1外，其它任何地方都沒有出現(xiàn)x3，x3是常數(shù)還是要求的參數(shù)？

贊一下

回復(fù)此樓

7樓2014-07-15 21:16:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

嗨皮你我

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 438.5
帖子: 41
在線: 16.8小時(shí)
蟲號(hào): 1975097
注冊(cè): 2012-09-04
專業(yè): 傳遞過程

引用回帖:

7樓: Originally posted by dingd at 2014-07-15 21:16:47
不懂你的專業(yè)，還是就1#代碼而言，x3起什么作用，除了：x1+x2+x3=1外，其它任何地方都沒有出現(xiàn)x3，x3是常數(shù)還是要求的參數(shù)？

我把x3 用1-x1-x2代替了，您看附件圖片中的關(guān)系式就行，那里面寫得很清楚，要求的參數(shù)有g(shù)12,g21,g31,g13,g32,g23六個(gè)，而方程中T是自變量，x1是函數(shù)，知道的初值就是在每一個(gè)x2下：x1隨T的變化

贊一下

回復(fù)此樓

8樓2014-07-16 14:14:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

嗨皮你我

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 438.5
帖子: 41
在線: 16.8小時(shí)
蟲號(hào): 1975097
注冊(cè): 2012-09-04
專業(yè): 傳遞過程

大神現(xiàn)身吶，求助！

回復(fù)此樓

9樓2014-07-18 09:52:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

蕭之夏

金蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 796.3
散金: 217
帖子: 208
在線: 149.3小時(shí)
蟲號(hào): 2160502
注冊(cè): 2012-12-02
性別: GG
專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)

你能給我解釋下這是什么方程嗎?
Function x1=exp(-H/R*(1/T-1/Tm)-lnr);
方程中H代表什么？

贊一下

回復(fù)此樓

10樓2016-03-29 20:20:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主嗨皮你我的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +9	秋有木北 2026-03-14	9/450	2026-03-18 19:49 by macy2011
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4	1孫悟空 2026-03-17	4/200	2026-03-18 17:59 by fivewind
[考研] 314求調(diào)劑 +8	無懈可擊的巨人 2026-03-12	8/400	2026-03-18 14:50 by haxia
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +7	手機(jī)用戶 2026-03-17	7/350	2026-03-18 14:34 by vgtyfty
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考博] 26申博 +4	八旬速覽 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321

24小時(shí)熱門版塊排行榜

[求助] 關(guān)于二元混合溶劑溶解度數(shù)據(jù)利用wilson、NRTL模擬的1stOPT代碼 已有2人參與

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

» 本帖附件資源列表

【答案】應(yīng)助回帖

[求助] 關(guān)于二元混合溶劑溶解度數(shù)據(jù)利用wilson、NRTL模擬的1stOPT代碼已有2人參與