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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 帮忙看下1stopt代码
51/1返回列表
查看: 726  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yarlyju

金虫 (小有名气)

[求助] 帮忙看下1stopt代码

nrtl方程拟合，怎么不行？


Parameters x(1),x(2),x(3);
Variable x1 ,x2 ,r,T;
ConstStr
tao21=(x(2)-x(1))/(8.314*T),
tao12=(x(1)-x(2))/(8.314*T),
g21=exp(-x(3)*tao21),
g12=exp(-x(3)*tao12)
Function r=exp(x2^2*(tao21*g21^2/(x1+x2*g21)^2+tao12*g12^2/(x2+x1*g12)^2));
data;
0.0317        0.9683        4.7509        366.29
0.0424        0.9576        4.717        364.32
0.0863        0.9137        3.5852        360.33
0.13        0.87        2.905        358.1
0.1666        0.8334        2.5357        357.16
0.2137        0.7863        2.1719        356.07
0.293        0.707        1.7506        354.97
0.3525        0.6475        1.5273        354.5
0.395        0.605        1.4167        353.99
0.4531        0.5469        1.3164        353.4
0.506        0.494        1.2182        353.01
0.5629        0.4371        1.1619        352.64
0.6142        0.3858        1.1126        352.33
0.6395        0.3605        1.0951        352.15
0.6794        0.3206        1.0715        351.95
0.724        0.276        1.0497        351.77
0.774        0.226        1.0267        351.57
0.8436        0.1564        1.0082        351.48
0.8612        0.1388        1.0058        351.44
0.902        0.098        1.0018        351.42
0.9464        0.0536        0.9941        351.48
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-08-03 11:05:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+2, 多谢应助 2012-08-03 11:48:05
yarlyju: 金币+3, ★★★★★最佳答案 2012-08-03 15:08:06
内容已删除
赞一下 回复此楼
2楼2012-08-03 11:32:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yarlyju

金虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2012-08-03 11:32:37
参数和变量定义重名了，另外：
tao21=(p2-p1)/(8.314*T),
tao12=(p1-p2)/(8.314*T),

可简化为：
tao21=p1/(8.314*T),
tao12=-tao21,

否则会出现参数不唯一现象

Parameters p1,p2,p3;
Variable x1 ,x2 ...

谢谢
回复此楼
3楼2012-08-03 15:08:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yarlyju

金虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2012-08-03 11:32:37
参数和变量定义重名了，另外：
tao21=(p2-p1)/(8.314*T),
tao12=(p1-p2)/(8.314*T),

可简化为：
tao21=p1/(8.314*T),
tao12=-tao21,

否则会出现参数不唯一现象

Parameters p1,p2,p3;
Variable x1 ,x2 ...

x1和x(1)会被识别成一个符号吗
赞一下 回复此楼
4楼2012-08-03 15:12:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
应该是的。
回复此楼
5楼2012-08-03 15:58:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yarlyju 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:22 by o3EMxHUAQDtS
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:12 by o3EMxHUAQDtS
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭