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

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

mushixilin

銅蟲 (初入文壇)

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

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

正在學(xué)習使用1stOpt，見到的示例都是只有單個因變量的非線性擬合。想要使用1stopt對三元液液相平衡的NRTL熱力學(xué)數(shù)據(jù)進行擬合，由于出現(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論文潤色/翻譯怎么收費? 已經(jīng)有158人回復(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ù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團隊招收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è)作家)

計算強帖: 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;可以實現(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;可以實現(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

將代碼進行了一定處理，仍然提示Fatal: Line too long (more than 1023 characters)，Compile failed, check your program codes please!。是因為最后的公式過于復(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)驗: +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ù)擬合，你的代碼有點復(fù)雜，軟件功能不夠。
需要正版高版本軟件來計算。

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

贊一下

回復(fù)此樓

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

mushixilin

一團和氣

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ù)擬合，你的代碼有點復(fù)雜，軟件功能不夠。
需要正版高版本軟件來計算。

謝謝前輩解答。
試了試，字符數(shù)多但擬合公式簡單的sharedmodel問題可以解決，我這個公式確實有點兒復(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è): 有機化工

引用回帖:

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

請問你的問題解決了嗎？

贊一下

回復(fù)此樓

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

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

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

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

返回列表

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

[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +3	AZMK 2026-03-27	5/250	2026-03-28 16:19 by xxxsssccc
[考研] 求化學(xué)調(diào)劑 +4	wulanna 2026-03-28	4/200	2026-03-28 13:37 by 唐沐兒
[考研] 070300求調(diào)劑306分 +4	26要上岸 2026-03-27	4/200	2026-03-28 13:06 by 唐沐兒
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +5	dongfang59 2026-03-27	5/250	2026-03-28 04:53 by wxiongid
[考研] 303求調(diào)劑 +7	安憶靈 2026-03-22	8/400	2026-03-27 11:46 by sanrepian
[論文投稿] Journal of Mechanical Science and Technology +3	Russ_ss 2026-03-25	5/250	2026-03-27 10:49 by 陸小果畫大餅
[考研] 284求調(diào)劑 +11	junqihahaha 2026-03-26	12/600	2026-03-27 04:37 by wxiongid
[考研] 071000生物學(xué)求調(diào)劑，初試成績343 +6	小小甜面團 2026-03-25	6/300	2026-03-26 23:01 by 不吃魚的貓
[考研] 349求調(diào)劑 +4	李木子啊哈哈 2026-03-25	4/200	2026-03-26 22:49 by fmesaito
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5	星沉uu 2026-03-26	6/300	2026-03-26 19:02 by macy2011
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4	JKSOIID 2026-03-26	4/200	2026-03-26 15:58 by 不吃魚的貓
[考研] 考研調(diào)劑 +6	來好運來來來 2026-03-21	7/350	2026-03-25 22:43 by 418490947
[考研] 0854人工智能方向招收調(diào)劑 +4	章小魚567 2026-03-24	4/200	2026-03-25 13:29 by 2177681040
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +5	WW.' 2026-03-23	7/350	2026-03-24 14:50 by sprinining
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎？ +3	lishucheng1 2026-03-22	5/250	2026-03-24 08:22 by gltch
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5	vv迷 2026-03-21	8/400	2026-03-22 14:20 by ColorlessPI

24小時熱門版塊排行榜

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

» 猜你喜歡

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

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

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

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

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

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

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