| 查看: 377 | 回复: 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文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有19人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有15人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有19人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有19人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有24人回复
-
面上提前没消息，有中的吗
已经有11人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有14人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有12人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有14人回复
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
50