| 查看: 413 | 回复: 2 | |||
[求助]
利用二元液液平衡实验数据模拟NRTL中的能量参数如何编程 已有1人参与
|
|
MATLAB太弱，还请帮帮忙 Data=[303.15 0.7766 0.02522]; T=Data(:,1); x11=Data(:,2); %x11=组分1在第1相中的摩尔分数(实验值） x12=Data(:,3); %x12=组分1在第2相中的摩尔分数(实验值） x21=1-x11; x22=1-x12; alfa=0.2; R=8.314; tao12=g12./R./T; tao21=g21./R./T; G12=exp(-alfa*tao12); G21=exp(-alfa*tao21); x21cal=1-x11cal; %x11cal=组分1在第1相中的摩尔分数（计算值） x22cal=1-x12cal; %x11cal=组分2在第1相中的摩尔分数（计算值） garma11=exp(x21cal.^2.*(tao21.*(G21./(x11cal+x21cal.*G21)).^2+tao12.*G12./(x21cal+x11cal.*G12).^2)); garma12=exp(x22cal.^2.*(tao21.*(G21./(x12cal+x22cal.*G21)).^2+tao12.*G12./(x22cal+x12cal.*G12).^2)); garma21=exp(x11cal.^2.*(tao12.*(G12./(x21cal+x11cal.*G12)).^2+tao21.*G21./(x11cal+x21cal.*G21).^2)); garma22=exp(x12cal.^2.*(tao12.*(G12./(x22cal+x12cal.*G12)).^2+tao21.*G21./(x12cal+x22cal.*G21).^2)); garma11.*x11cal=garma12.*x12cal; garma21*x21cal=garma22*x22cal 目标函数是fmin sigma=(x11-x11cal).^2+(x12-x12cal).^2+(x21-x21cal).^2+(x22-x22cal).^2; |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于二元混合溶剂溶解度数据利用wilson、NRTL模拟的1stOPT代码
已经有10人回复
- 药动实验中，血浆样品用二氯甲烷液液萃取的多吗？
已经有15人回复
- 液液萃取多环芳烃
已经有9人回复
- matlab动力学模型编程求助
已经有6人回复
- aspen 数据回归求助
已经有3人回复
- 用matlab如何解方程？
已经有6人回复
- HSCCC溶剂系统
已经有6人回复
- 请问matlab拟合二元方程的问题
已经有3人回复
- 用Chemcad怎样回归二元交互作用参数？？
已经有4人回复
- 汽液平衡数据怎么算
已经有15人回复
- 二元幂函数参数求解
已经有10人回复
- 请问abaqus做90度ecap挤压模拟分析步应该怎么选择
已经有4人回复
- opls下L-J参数转换成amber下L-J参数的问题
已经有9人回复
- 请问哪位仁兄知道用IAST模型计算二元混合气体吸附的选择性？
已经有3人回复
- ASPEN数据回归
已经有7人回复
- 重金加人民币求拟合wilson方程和NRTL方程参数的程序
已经有17人回复
- 【求助】求个能计算粘度等参数的水蒸气性质计算软件
已经有6人回复
- 【讨论】请教一个气-液-固三相流模拟可能性的问题
已经有6人回复
- 【求助】分散液液微萃取形成乳浊液的关键是什么
已经有10人回复
- 【求助】Aspen二元交互参数求助
已经有4人回复
zhchh008
金虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 1613.2
- 散金: 861
- 红花: 13
- 帖子: 833
- 在线: 202.6小时
- 虫号: 99456
- 注册: 2005-11-11
- 性别: GG
- 专业: 资源化工
【答案】应助回帖
感谢参与，应助指数 +1
| http://muchong.com/bbs/viewthread.php?tid=7999270&authorid=99456,,你自己的帖子，，还没有明白？ |
2楼2014-10-13 21:56:38
3楼2014-10-13 22:24:15