应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 利用二元液液平衡实验数据模拟NRTL中的能量参数如何编程
31/1返回列表
查看: 412  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小指甲123

新虫 (小有名气)

[求助] 利用二元液液平衡实验数据模拟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;
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-10-13 19:21:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhchh008

金虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
http://muchong.com/bbs/viewthread.php?tid=7999270&authorid=99456,,你自己的帖子，，还没有明白？
赞一下 回复此楼
2楼2014-10-13 21:56:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小指甲123

新虫 (小有名气)
引用回帖:
2楼: Originally posted by zhchh008 at 2014-10-13 21:56:38
http://muchong.com/bbs/viewthread.php?tid=7999270&authorid=99456,,你自己的帖子，，还没有明白？

不是 的，这两个不一样，之前那个帖子是已知参数，解方程组问题，这一个是来回归参数的，感觉编程好难
赞一下 回复此楼
3楼2014-10-13 22:24:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小指甲123 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭