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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

mushixilin

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5.5
帖子: 43
在線: 29.2小時
蟲號: 3099019
注冊: 2014-03-29

[求助] 如何使用1stOpt進(jìn)行三自變量三因變量的非線性擬合已有2人參與

正在學(xué)習(xí)使用1stOpt，見到的示例都是只有單個因變量的非線性擬合。想要使用1stopt對三元液液相平衡的NRTL熱力學(xué)數(shù)據(jù)進(jìn)行擬合，由于出現(xiàn)三個因變量，不太確定如何編程。編寫的程序如下，如不使用sharedmodel，可以得到結(jié)果，但不理想；如使用sharedmodel，軟件無法計算。還望前輩可以給予一定指導(dǎo)。

Title "NRTL simulation";
//sharedmodel;
Parameters alpha[0,1],tau12,tau13,tau21,tau23,tau31,tau32;
Variables x1,x2,x3,y1,y2,y3;

ConstStr F12=exp(-alpha*tau12);
ConstStr F13=exp(-alpha*tau13);
ConstStr F21=exp(-alpha*tau21);
ConstStr F23=exp(-alpha*tau23);
ConstStr F31=exp(-alpha*tau31);
ConstStr F32=exp(-alpha*tau32);

ConstStr lngammax1=(x1*0*1+x2*tau21*F21+x3*tau31*F31)/(x1*1+x2*F21+x3*F31)+
      (x1*1/(x1*1+x2*F21+x3*F31)*(0-(x1*0*1+x2*tau21*F21+x3*tau31*F31)/(x1*1+x2*F21+x3*F31)))+
      (x2*F12/(x1*F12+x2*1+x3*F32)*(tau12-(x1*tau12*F12+x2*0*1+x3*tau32*F32)/(x1*F12+x2*0+x3*F32)))+
      (x3*F13/(x1*F12+x2*F23+x3*1)*(tau13-(x1*tau13*F13+x2*tau23*F23+x3*0*1)/(x1*F13+x2*F23+x3*1)));
ConstStr lngammax2=(x1*tau12*F12+x2*0*1+x3*tau32*F32)/(x1*F12+x2*0+x3*F32)+
      (x1*F21/(x1*1+x2*F21+x3*F31)*(tau21-(x1*0*1+x2*tau21*F21+x3*tau31*F31)/(x1*1+x2*F21+x3*F31)))+
      (x2*1/(x1*F12+x2*1+x3*F32)*(0-(x1*tau12*F12+x2*0*1+x3*tau32*F32)/(x1*F12+x2*0+x3*F32)))+
      (x3*F23/(x1*F12+x2*F23+x3*1)*(tau23-(x1*tau13*F13+x2*tau23*F23+x3*0*1)/(x1*F13+x2*F23+x3*1)));
ConstStr lngammax3=(x1*tau13*F13+x2*tau23*F23+x3*0*1)/(x1*F13+x2*F23+x3*1)+
      (x1*F31/(x1*1+x2*F21+x3*F31)*(tau31-(x1*0*1+x2*tau21*F21+x3*tau31*F31)/(x1*1+x2*F21+x3*F31)))+
      (x2*F32/(x1*F12+x2*1+x3*F32)*(tau32-(x1*tau12*F12+x2*0*1+x3*tau32*F32)/(x1*F12+x2*0+x3*F32)))+
      (x3*1/(x1*F12+x2*F23+x3*1)*(0-(x1*tau13*F13+x2*tau23*F23+x3*0*1)/(x1*F13+x2*F23+x3*1)));

ConstStr gammax1=exp(lngammax1);
ConstStr gammax2=exp(lngammax2);
ConstStr gammax3=exp(lngammax3);

ConstStr lngammay1=(y1*0*1+y2*tau21*F21+y3*tau31*F31)/(y1*1+y2*F21+y3*F31)+
      (y1*1/(y1*1+y2*F21+y3*F31)*(0-(y1*0*1+y2*tau21*F21+y3*tau31*F31)/(y1*1+y2*F21+y3*F31)))+
      (y2*F12/(y1*F12+y2*1+y3*F32)*(tau12-(y1*tau12*F12+y2*0*1+y3*tau32*F32)/(y1*F12+y2*0+y3*F32)))+
      (y3*F13/(y1*F12+y2*F23+y3*1)*(tau13-(y1*tau13*F13+y2*tau23*F23+y3*0*1)/(y1*F13+y2*F23+y3*1)));
ConstStr lngammay2=(y1*tau12*F12+y2*0*1+y3*tau32*F32)/(y1*F12+y2*0+y3*F32)+
      (y1*F21/(y1*1+y2*F21+y3*F31)*(tau21-(y1*0*1+y2*tau21*F21+y3*tau31*F31)/(y1*1+y2*F21+y3*F31)))+
      (y2*1/(y1*F12+y2*1+y3*F32)*(0-(y1*tau12*F12+y2*0*1+y3*tau32*F32)/(y1*F12+y2*0+y3*F32)))+
      (y3*F23/(y1*F12+y2*F23+y3*1)*(tau23-(y1*tau13*F13+y2*tau23*F23+y3*0*1)/(y1*F13+y2*F23+y3*1)));
ConstStr lngammay3=(y1*tau13*F13+y2*tau23*F23+y3*0*1)/(y1*F13+y2*F23+y3*1)+
      (y1*F31/(y1*1+y2*F21+y3*F31)*(tau31-(y1*0*1+y2*tau21*F21+y3*tau31*F31)/(y1*1+y2*F21+y3*F31)))+
      (y2*F32/(y1*F12+y2*1+y3*F32)*(tau32-(y1*tau12*F12+y2*0*1+y3*tau32*F32)/(y1*F12+y2*0+y3*F32)))+
      (y3*1/(y1*F12+y2*F23+y3*1)*(0-(y1*tau13*F13+y2*tau23*F23+y3*0*1)/(y1*F13+y2*F23+y3*1)));

ConstStr gammay1=exp(lngammay1);
ConstStr gammay2=exp(lngammay2);
ConstStr gammay3=exp(lngammay3);

Function y1=x1*gammax1/gammay1;
            y2=x2*gammax2/gammay2;
            y3=x3*gammax3/gammay3;

Data;
0.9561 0.0289 0.0150 0.0010 0.0011 0.9979
0.9101 0.0739 0.0160 0.0013 0.0028 0.9959
0.7982 0.1675 0.0343 0.0010 0.0071 0.9919
0.6980 0.2558 0.0462 0.0026 0.0164 0.9810
0.6461 0.2941 0.0598 0.0018 0.0253 0.9729
0.5686 0.3488 0.0826 0.0011 0.0404 0.9585
0.5194 0.3921 0.0885 0.0037 0.0525 0.9438
0.4919 0.4101 0.0980 0.0045 0.0583 0.9372
0.4664 0.4319 0.1017 0.0079 0.0688 0.9233
0.4497 0.4410 0.1093 0.0111 0.0736 0.9153
0.4374 0.4498 0.1128 0.0125 0.0884 0.8991
0.4189 0.4515 0.1296 0.0186 0.0942 0.8872

回復(fù)此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有184人回復(fù)
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)

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

利用Excel數(shù)據(jù)前處理，結(jié)合1stopt數(shù)據(jù)自定義擬合初值，美美origin數(shù)據(jù)作圖分析！已經(jīng)有72人回復(fù)

1樓 2017-02-17 15:15:08

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

dingd

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

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
mushixilin: 金幣+15, ★★★很有幫助 2017-02-18 09:50:46

sharedmodel;可以實(shí)現(xiàn)多因子參數(shù)共享擬合，1.5的估計不行。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

mushixilin

2樓2017-02-17 17:15:32

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

mushixilin

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5.5
帖子: 43
在線: 29.2小時
蟲號: 3099019
注冊: 2014-03-29

送紅花一朵

引用回帖:

2樓: Originally posted by dingd at 2017-02-17 17:15:32
sharedmodel;可以實(shí)現(xiàn)多因子參數(shù)共享擬合，1.5的估計不行。

謝謝，得到大神的解答。在小木蟲上看了很多類似的帖子，有提到1.5不可以使用sharedmodel，也有蟲友說lib文件夾可以用，下載后使用了一下，發(fā)現(xiàn)1.5可以使用sharedmodel功能，可以計算較為簡單的sharedmodel問題。但是對于我這個問題，可能由于代碼太長，無法計算，提示如下：Fatal: Line too long (more than 1023 characters)，Compile failed, check your program codes please!。

贊一下

回復(fù)此樓

3樓2017-02-17 22:53:23

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

mushixilin

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5.5
帖子: 43
在線: 29.2小時
蟲號: 3099019
注冊: 2014-03-29

將代碼進(jìn)行了一定處理，仍然提示Fatal: Line too long (more than 1023 characters)，Compile failed, check your program codes please!。是因?yàn)樽詈蟮墓竭^于復(fù)雜嗎？

Title "NRTL simulation";
sharedmodel;
Parameters alpha[0,1],t12,t13,t21,t23,t31,t32;
Variables x1,x2,x3,y1,y2,y3;
ConstStr G12=exp(-alpha*t12),
      G13=exp(-alpha*t13),
      G21=exp(-alpha*t21),
      G23=exp(-alpha*t23),
      G31=exp(-alpha*t31),
      G32=exp(-alpha*t32),
      a=(x2*t21*G21+x3*t31*G31)/(x1+x2*G21+x3*G31),
      b=(x1*t12*G12+x3*t32*G32)/(x1*G12+x2+x3*G32),
      c=(x1*t13*G13+x2*t23*G23)/(x1*G13+x2*G23+x3),
      d=x1/(x1+x2*G21+x3*G31),
      e=x2/(x1*G12+x2+x3*G32),
      f=x3/(x1*G13+x2*G23+x3),
      j=(y2*t21*G21+y3*t31*G31)/(y1+y2*G21+y3*G31),
      k=(y1*t12*G12+y3*t32*G32)/(y1*G12+y2+y3*G32),
      l=(y1*t13*G13+y2*t23*G23)/(y1*G13+y2*G23+y3),
      m=y1/(y1+y2*G21+y3*G31),
      n=y2/(y1*G12+y2+y3*G32),
      o=y3/(y1*G13+y2*G23+y3);
Conststr lngammax1=a-d*a+e*G12*(t12-b)+f*G13*(t13-c),
      lngammax2=b+d*G12*(t21-a)-e*b+f*G23*(t23-c),
      lngammax3=c+d*G31*(t31-a)+e*G32*(t32-b)-f*c,
      lngammay1=j-m*j+n*G12*(t12-k)+o*G13*(t13-l),
      lngammay2=k+m*G12*(t21-j)-n*k+o*G23*(t23-l),
      lngammay3=l+m*G31*(t31-j)+n*G32*(t32-k)-o*l;
Conststr gammax1=exp(lngammax1),
      gammax2=exp(lngammax2),
      gammax3=exp(lngammax3),
      gammay1=exp(lngammay1),
      gammay2=exp(lngammay2),
      gammay3=exp(lngammay3);
Function y1=x1*gammax1/gammay1;
      y2=x2*gammax2/gammay2;
      y3=x3*gammax3/gammay3;
Data;
0.9561 0.0289 0.0150 0.0010 0.0011 0.9979
0.9101 0.0739 0.0160 0.0013 0.0028 0.9959
0.7982 0.1675 0.0343 0.0010 0.0071 0.9919
0.6980 0.2558 0.0462 0.0026 0.0164 0.9810
0.6461 0.2941 0.0598 0.0018 0.0253 0.9729
0.5686 0.3488 0.0826 0.0011 0.0404 0.9585
0.5194 0.3921 0.0885 0.0037 0.0525 0.9438
0.4919 0.4101 0.0980 0.0045 0.0583 0.9372
0.4664 0.4319 0.1017 0.0079 0.0688 0.9233
0.4497 0.4410 0.1093 0.0111 0.0736 0.9153
0.4374 0.4498 0.1128 0.0125 0.0884 0.8991
0.4189 0.4515 0.1296 0.0186 0.0942 0.8872

贊一下

回復(fù)此樓

4樓2017-02-18 16:20:57

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

lipenggg

專家顧問 (正式寫手)

專家經(jīng)驗(yàn): +65
應(yīng)助: 122 (高中生)
金幣: 4841.3
紅花: 43
帖子: 722
在線: 435.7小時
蟲號: 3929201
注冊: 2015-06-17
性別: GG
專業(yè): 流體力學(xué)
管轄: 計算模擬

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
mushixilin: 金幣+5, ★★★很有幫助 2017-02-18 22:05:54

1.5版本太低了，即使能計算共享參數(shù)擬合，你的代碼有點(diǎn)復(fù)雜，軟件功能不夠。
需要正版高版本軟件來計算。

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

mushixilin

一團(tuán)和氣

5樓2017-02-18 18:09:04

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

mushixilin

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5.5
帖子: 43
在線: 29.2小時
蟲號: 3099019
注冊: 2014-03-29

送紅花一朵

引用回帖:

5樓: Originally posted by lipenggg at 2017-02-18 18:09:04
1.5版本太低了，即使能計算共享參數(shù)擬合，你的代碼有點(diǎn)復(fù)雜，軟件功能不夠。
需要正版高版本軟件來計算。

謝謝前輩解答。
試了試，字符數(shù)多但擬合公式簡單的sharedmodel問題可以解決，我這個公式確實(shí)有點(diǎn)兒復(fù)雜，算不了。
之前在用matlab，可惜對初值要求太大，專程學(xué)了1stOpt，沒想到還是解不了。我再想想辦法。

贊一下

回復(fù)此樓

6樓2017-02-18 22:11:33

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

小靈魚兒

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1343.9
散金: 285
紅花: 1
帖子: 147
在線: 138.9小時
蟲號: 1722780
注冊: 2012-03-28
性別: MM
專業(yè): 有機(jī)化工

引用回帖:

6樓: Originally posted by mushixilin at 2017-02-18 22:11:33
謝謝前輩解答。
試了試，字符數(shù)多但擬合公式簡單的sharedmodel問題可以解決，我這個公式確實(shí)有點(diǎn)兒復(fù)雜，算不了。
之前在用matlab，可惜對初值要求太大，專程學(xué)了1stOpt，沒想到還是解不了。我再想想辦法。...

請問你的問題解決了嗎？

贊一下

回復(fù)此樓

不打無準(zhǔn)備之戰(zhàn)

7樓2017-03-19 17:47:48

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

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

返回列表

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

[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8	AZMK 2026-03-25	12/600	2026-03-28 18:15 by AZMK
[考研] 321求調(diào)劑 +6	璞玉~~ 2026-03-25	7/350	2026-03-28 17:48 by 璞玉~~
[考研] 0703化學(xué)求調(diào)劑 +9	奶油草莓. 2026-03-22	10/500	2026-03-28 13:30 by 唐沐兒
[考研] 070300求調(diào)劑306分 +4	26要上岸 2026-03-27	4/200	2026-03-28 13:06 by 唐沐兒
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +6	腳滑的守法公民 2026-03-27	6/300	2026-03-28 11:02 by gjlllb
[考研] 311求調(diào)劑 +3	希望上岸阿小楊 2026-03-23	3/150	2026-03-28 07:57 by 熱情沙漠
[考研] 275求調(diào)劑 +10	jjjjjjjjjjl 2026-03-27	10/500	2026-03-27 23:47 by barnett0632
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。�！ +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[有機(jī)交流] 高溫高壓反應(yīng)求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3	zzffylgg 2026-03-24	3/150	2026-03-27 10:42 by shangxh
[考研] 一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑 +4	舉個栗子oi 2026-03-24	5/250	2026-03-26 23:15 by 不吃魚的貓
[考研] 271求調(diào)劑 +6	生如夏花… 2026-03-22	6/300	2026-03-26 16:48 by 張凱十八號
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4	gdlf9999 2026-03-24	4/200	2026-03-25 23:01 by boxking200
[考研] 0854AI CV方向招收調(diào)劑 +4	章小魚567 2026-03-23	4/200	2026-03-25 17:04 by CoderLoser
[考研] 318求調(diào)劑 +5	plum李子 2026-03-21	8/400	2026-03-25 09:26 by aa331100
[考研] B區(qū)考研調(diào)劑 +4	yqdszhdap－ 2026-03-22	5/250	2026-03-25 08:51 by baoball
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3	an10101 2026-03-24	7/350	2026-03-25 00:00 by 山鬼0-
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106

24小時熱門版塊排行榜

[求助] 如何使用1stOpt進(jìn)行三自變量三因變量的非線性擬合 已有2人參與

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

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

» 本帖已獲得的紅花（最新10朵）

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

» 本帖已獲得的紅花（最新10朵）

[求助] 如何使用1stOpt進(jìn)行三自變量三因變量的非線性擬合已有2人參與

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