| 查看: 412 | 回复: 1 | |||
[交流]
固液平衡的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 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有18人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有18人回复
-
2027博士申请-临床检验诊断学、病原生物学、微生物方向
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有22人回复
-
越来越担心硕士毕业后没有竞争力了
已经有14人回复
-
论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章
已经有3人回复
-
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请
已经有3人回复
-
闲聊
已经有6人回复
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+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