当前位置: 首页 > 计算模拟 >NRTL方程拟合敛到参数全为0

NRTL方程拟合敛到参数全为0

作者 09天岚
来源: 小木虫 300 6 举报帖子
+关注

各位大佬们,我在1stopt上拟合固液相平衡的NRTL方程,为什么收敛到参数全为0啊。。。有没有大佬帮忙看看方程哪里出了问题。。。谢谢了!
数据是找的文献的数据。
拿这套数据拟合其他模型都没什么问题。。
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方程原形式。
急急急   求求大家帮忙分析分析!

NRTL方程拟合敛到参数全为0
PWN3F1)AM5{5W4{)S}0X1OB.png 返回小木虫查看更多

今日热帖
  • 精华评论
  • 独孤神宇

    方程写成y=……
    因变量在方程左边,变量,数据也要对应

  • 09天岚

    上边写错了。。是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)));

  • hzlhm

    引用回帖:
    4楼: Originally posted by 09天岚 at 2021-05-26 12:31:26
    上边写错了。。是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函数,可以得到
    a= 0.088241733690097801 ,b= -0.052836750140336099 ,c= -0.024187951865466500 ,r=0
    NRTL方程拟合敛到参数全为0-1
    123a.jpg


  • hzlhm

    引用回帖:
    6楼: Originally posted by 09天岚 at 2021-05-26 22:22:16
    其实我可能表达不清楚  这个r是活度系数。。就是想要通过拟合得到一个可用的r用以计算溶解度。。。所以正常情况下r不可以为0
    ...

    计算结果
    NRTL方程拟合敛到参数全为0-2
    123a.jpg

猜你喜欢