| 查看: 348 | 回复: 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 |
回复此楼» 猜你喜欢
一志愿重庆大学085700资源与环境总分308求调剂
已经有3人回复
-
南昌大学材料专硕311分求调剂
已经有6人回复
-
316求调剂
已经有6人回复
-
346求调剂[0856]
已经有7人回复
-
一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！
已经有3人回复
-
310求调剂
已经有3人回复
-
机械专硕299求调剂至材料
已经有4人回复
-
070300化学319求调剂
已经有7人回复
-
08工科 320总分 求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有12人回复
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
