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

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

cheerfulchen

鐵桿木蟲 (小有名氣)

應(yīng)助: 6 (幼兒園)
金幣: 4934.6
紅花: 2
帖子: 127
在線: 62.4小時
蟲號: 946644
注冊: 2010-01-22
性別: GG
專業(yè): 分離過程

[求助] 用lsqnonlin擬合 NRTL三元方程回歸參數(shù) 出現(xiàn)錯誤求解脫

function NRTLBT

a0=[1 1 1 1 1 1 1 1 1 ];
x=Lsqnonlin(@csjs,a0)
function z=csjs
x1 =[0.0022 0.0358 0.0134 0.0097 0.0073 0.0035];
x2 =[0.0224 0.339 0.1227 0.0857 0.0721 0.032];
x3=1-x1-x2;
y1e = [0.0098 0.0379 0.0329 0.037 0.0316 0.0168];
y2e = [0.0098 0.0379 0.0329 0.037 0.0316 0.0168];
y3e=1-y1e-y2e;
p =[2561 3112 3454 3600 3600 3710];
ps1 = [2530.011324 3731.190276 3731.190276 3731.190276 3731.190276 3731.190276];
ps2 = [38.27422323 91.99677979 91.99677979 91.99677979 91.99677979 91.99677979];
ps3 = [19.78428868 47.08467815 47.08467815 47.08467815 47.08467815 47.08467815];
T=[333 353 353 353 353 353];
R=8.314;
for i = 1:6
%G12=exp(-a(1)*a(4)/R/T(i));tao12=a(4)/R/T(i);%a12=a(1),a13=a(2),a23=a(3),g12-g22=a(4),g13-g33=a(5),g23-g33=a(6)
%G13=exp(-a(2)*a(5)/R/T(i));tao13=a(5)/R/T(i);
%G23=exp(-a(3)*a(6)/R/T(i));tao23=a(6)/R/T(i);
a=zero[];
tao12=a(1)/R/T(i);tao21=a(2)/R/T(i);
tao13=a(3)/R/T(i);tao31=a(4)/R/T(i);
tao23=a(5)/R/T(i);tao32=a(6)/R/T(i);
G12=exp(-a(7)*a(1)/R/T(i));G21=exp(-a(7)*a(2)/R/T(i));
G13=exp(-a(8)*a(3)/R/T(i));G31=exp(-a(8)*a(4)/R/T(i));
G23=exp(-a(9)*a(5)/R/T(i));G32=exp(-a(9)*a(6)/R/T(i));
y1c(i) = x1(i)*ps1(i)/p(i)*exp((G21*x2(i)+G31*x3(i))*(tao21*G21*x2(i)+tao31*G31*x3(i))/(x1(i)+x2(i)*G21+x3(i)*G31)^2)...
+(tao12*G12*x2(i)^2+G12*G32*x2(i)*x3(i)*(tao12-tao32))/(x1(i)*G12+x2(i)+x3(i)*G32)^2+(tao13*G13*x3(i)^2+G13*G23*x2(i)*x3(i)*(tao13-tao23))/(x1(i)*G13+x2(i)*G23+x3(i))^2;
y2c(i) = x2(i)*ps2(i)/p(i)*exp((G12*x1(i)+G32*x3(i))*(tao12*G12*x1(i)+tao32*G32*x3(i))/(x1(i)*G12+x2(i)+x3(i)*G32)^2)...
+(tao23*G23*x3(i)^2+G13*G23*x1(i)*x3(i)*(tao23-tao13))/(x1(i)*G13+x2(i)*G23+x3(i))^2+(tao21*G21*x1(i)^2+G21*G31*x1(i)*x3(i)*(tao21-tao31))/(x1(i)+x2(i)*G21+x3(i)*G31)^2;
y3c(i) = x3(i)*ps3(i)/p(i)*exp((G13*x1(i)+G23*x2(i))*(tao13*G13*x1(i)+tao23*G23*x2(i))/(x1(i)*G13+x2(i)*G23+x3(i))^2)...
+(tao21*G31*x1(i)^2+G21*G31*x1(i)*x2(i)*(tao31-tao21))/(x1(i)+x2(i)*G21+x3(i)*G31)^2+(tao32*G32*x2(i)^2+G12*G32*x1(i)*x2(i)*(tao32-tao12))/(x1(i)*G12+x3(i)*G32+x2(i))^2;

z(i)=(y1e(i)-y1c(i))^2 + (y2e(i) - y2c(i))^2+(y3e(i) - y3c(i))^2;

end

%Matlab函數(shù)：lsqnonlin('NRTL',x0)

[ Last edited by cheerfulchen on 2012-4-6 at 09:43 ]

回復(fù)此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有51人回復(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ù)試實施細(xì)則已出。歡迎加入光子學(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)價值貼推薦，對您同樣有幫助:

ABAQUS安裝成功后，運行出現(xiàn)問題，請大家?guī)兔匆幌?/a> 已經(jīng)有14人回復(fù)

添加聚乙烯醇后出現(xiàn) 團聚，求教原因已經(jīng)有10人回復(fù)
材料內(nèi)出現(xiàn)納米晶已經(jīng)有6人回復(fù)
辰華-恒電流計時電位法測催化劑耐久性如何設(shè)參數(shù)？已經(jīng)有7人回復(fù)
求高手看下這個DNA單鏈的結(jié)構(gòu)錯誤何在？已經(jīng)有13人回復(fù)
氘代水中能出現(xiàn)活潑氫么？已經(jīng)有5人回復(fù)
這是什么錯誤？已經(jīng)有3人回復(fù)
R語言圖書-求解差分方程和并行計算！已經(jīng)有43人回復(fù)
化藥6類四種介質(zhì)溶出度擬合時用做方法學(xué)嗎？已經(jīng)有8人回復(fù)
求助這個非線性方程改怎么解��？？已經(jīng)有4人回復(fù)
求助：origin高斯擬合已經(jīng)有10人回復(fù)
擬合過過彈性模量的沒？已經(jīng)有43人回復(fù)
求繩子的曲線方程已經(jīng)有6人回復(fù)
求matlab 解方程的問題已經(jīng)有5人回復(fù)
對數(shù)據(jù)進行曲線擬合后如何得到擬合誤差已經(jīng)有30人回復(fù)
脈沖微分方程內(nèi)涵已經(jīng)有3人回復(fù)
重金加人民幣求擬合wilson方程和NRTL方程參數(shù)的程序已經(jīng)有17人回復(fù)

1樓 2012-04-05 20:57:21

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

chaiqing123

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 679.7
紅花: 1
帖子: 71
在線: 65.2小時
蟲號: 1704607
注冊: 2012-03-20
性別: GG
專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)

目標(biāo)函數(shù)錯了本身就包括求和累加

贊一下

回復(fù)此樓

2樓2012-04-26 09:01:58

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

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