NRTL方程擬合斂到參數(shù)全為0
各位大佬們,我在1stopt上擬合固液相平衡的NRTL方程,為什么收斂到參數(shù)全為0啊。。。有沒有大佬幫忙看看方程哪里出了問題。。。謝謝了!
數(shù)據(jù)是找的文獻的數(shù)據(jù)。
拿這套數(shù)據(jù)擬合其他模型都沒什么問題。。
parameters a,b,c,r;
variable y,x;
function ln(r)=exp(((1-x)^2)*((b/(8.314*y)*exp(c*a/(8.314*y))/((x+(1-x)*exp(c*a/(8.314*y)))^2))+((a/(8.314*y))*exp(c*b/(8.314*y)))/(((1-x)+(1-x)*exp(c*b/(8.314*y)))^2)));
data;
//y,x
278.15 0.001541
283.15 0.002004
288.15 0.002704
293.15 0.003392
298.15 0.004415
303.15 0.005561
308.15 0.007136
313.15 0.009151
318.15 0.01113
圖為NRTL方程原形式。
急急急 求求大家?guī)兔Ψ治龇治觯?br />

PWN3F1)AM5{5W4{)S}0X1OB.png
返回小木蟲查看更多
今日熱帖
京公網(wǎng)安備 11010802022153號
方程寫成y=……
因變量在方程左邊,變量,數(shù)據(jù)也要對應
棒
上邊寫錯了。。是function r=exp(((1-x)^2)*((b/(8.314*y)*exp(c*a/(8.314*y))/((x+(1-x)*exp(c*a/(8.314*y)))^2))+((a/(8.314*y))*exp(c*b/(8.314*y)))/(((1-x)+(1-x)*exp(c*b/(8.314*y)))^2)));
利用matlab的lsqcurvefit函數(shù),可以得到
a= 0.088241733690097801 ,b= -0.052836750140336099 ,c= -0.024187951865466500 ,r=0
123a.jpg
,
計算結(jié)果
123a.jpg