24小时热门版块排行榜    

查看: 529  |  回复: 2

小指甲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;
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 YXM1FDCSE8Bn 2026-11-01 6/300 2026-11-02 09:57 by dcc2yOIRGGTT
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 YXM1FDCSE8Bn 2026-10-31 5/250 2026-11-02 09:42 by dcc2yOIRGGTT
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 YXM1FDCSE8Bn 2026-10-31 6/300 2026-11-02 09:36 by dcc2yOIRGGTT
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 09:31 by dcc2yOIRGGTT
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 06:29 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 05:54 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:34 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:24 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
请填处理意见