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

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

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)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有27人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有261人回復(fù)
UJN物理學(xué)專業(yè)調(diào)劑已經(jīng)有12人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名已經(jīng)有0人回復(fù)
山西大同大學(xué)物理學(xué)碩士研究生招收調(diào)劑生已經(jīng)有0人回復(fù)
[調(diào)劑信息]211智能人工感知方向國家青年特聘專家課題組招收調(diào)劑研究生已經(jīng)有0人回復(fù)
211智能人工感知方向國家青年特聘專家課題組招收調(diào)劑研究生已經(jīng)有0人回復(fù)
[調(diào)劑信息]211智能人工感知方向國家青年特聘專家課題組招收調(diào)劑研究生已經(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的回帖

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!。是因為最后的公式過于復(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的回帖

查看全部 7 個回答

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;可以實現(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的回帖

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ā)自手機(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的回帖

查看全部 7 個回答

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

[考研] 071000生物學(xué)調(diào)劑 +6	知昭蔓 2026-04-02	6/300	2026-04-02 15:28 by sophie2180
[考研] 279求調(diào)劑 +4	傅文秋 2026-04-02	4/200	2026-04-02 15:02 by cal0306
[考研] 材料調(diào)劑 +6	一樣YWY 2026-04-02	6/300	2026-04-02 14:56 by liu823948201
[考研] 0832食品科學(xué)與工程學(xué)碩282調(diào)劑 +4	魚在水中游a 2026-04-02	7/350	2026-04-02 14:12 by baoball
[考研] 初試301，代碼085701環(huán)境工程，本碩一致，四六級已過，有二區(qū)一作，共發(fā)表5篇論文 +6	axibli 2026-04-01	6/300	2026-04-02 13:42 by Ecowxq666！
[考研] 290求調(diào)劑 +3	1314捧花 2026-04-02	3/150	2026-04-02 13:34 by Ecowxq666！
[考研] 08工科求調(diào)劑290分 +5	1314捧花 2026-04-02	8/400	2026-04-02 13:16 by 喬噠噠噠
[考研] 292求調(diào)劑 +14	是妍子也是研子 2026-03-30	14/700	2026-04-02 10:17 by sanrepian
[考研] 070300化學(xué)求調(diào)劑 +14	小黃鴨寶 2026-03-30	14/700	2026-04-02 10:07 by oooqiao
[考研] 362求調(diào)劑 +13	西南交材料專碩3 2026-03-31	13/650	2026-04-01 17:38 by JYD2011
[考研] 土木304求調(diào)劑 +5	頂級擦擦 2026-03-31	5/250	2026-04-01 08:15 by fdcxdystjk￥
[考研] 材料科學(xué)與工程求調(diào)劑 +13	深V宿舍吧 2026-03-29	13/650	2026-03-31 19:50 by Dyhoer
[考研] 263求調(diào)劑 +3	DDDDuu 2026-03-27	3/150	2026-03-31 16:21 by 土木碩士招生
[考研] 一志愿中海洋材料357 +4	麥恩莉. 2026-03-30	4/200	2026-03-31 14:35 by 記事本2026
[考研] 一志愿浙江大學(xué)工科動力工程370,數(shù)一121,專業(yè)課135，現(xiàn)在能去哪里 +3	080700調(diào)劑 2026-03-30	4/200	2026-03-31 12:00 by KLMY666
[考研] 083000學(xué)碩274求調(diào)劑 +12	Li李魚 2026-03-26	12/600	2026-03-31 10:01 by cal0306
[考研] 08工科求調(diào)劑286 +5	tgs_001 2026-03-28	5/250	2026-03-31 08:18 by 一只好果子?
[考研] 求調(diào)劑 +10	家佳佳佳佳佳 2026-03-29	10/500	2026-03-30 18:34 by 544594351
[考研] 071010 323 分求調(diào)劑 +3	Baekzhy 2026-03-27	3/150	2026-03-30 14:24 by andresqi
[考研] 081200-314 +3	LILIQQ 2026-03-27	4/200	2026-03-28 09:41 by 保護(hù)地球你我做?/a>

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

mushixilin

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

» 猜你喜歡

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

mushixilin

dingd

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

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

mushixilin

lipenggg

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

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

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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

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