24小时热门版块排行榜    

CyRhmU.jpeg
查看: 336  |  回复: 1

scut huang

新虫 (初入文坛)

[交流] 固液平衡的NRTL方程的程序怎么编写,学了好久还是一团糟啊?已有1人参与

function f= mynrtl(q)
syms a t12 t21 R g12 g21 G12 G21 i ;
T=[298.15 318.15 333.15 343.15];
x1=[0.5552 0.852 0.99688 0.999598];
x2=1-x1;
a=0.30;
R=8.314;
H=47700;
Tm=20;
beta0=[1,1];
A=nlinfit(x1,x2,r1,@mynrtl,beta0);
  for i=1:4
      t12=g12/R*T(i);
      t21=g21/R*T(i);
      G12=exp(-a*t12);
      G21=exp(-a*t21);
      r1(i)=exp((x2(i)^2)*(t21*(G21/(x1(i)+x2(i)*G21))^2+t12*G12/(x2(i)+x1(i)*G12)^2));
      xcal(i)=exp(H*(1/Tm-1/T))/r1(i);
      F=sum(((x1(i)-xcal(i))^2);
  end
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็


小木虫: 金币+0.5, 给个红包,谢谢回帖
NRTL 不是 Non-random two-liquid model

你用来做固液平衡?
https://en.wikipedia.org/wiki/Non-random_two-liquid_model

建议分解自己的问题为一个个小问题,
然后逐个问题地、自己用搜索引擎找答案。
ส็็็็็็็็็็็็็็็็็็็็
2楼2016-04-19 10:20:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 scut huang 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见