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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 1stOpt回归模型参数的问题
51/1返回列表
查看: 1141  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

最爱明日香

铁虫 (初入文坛)

[求助] 1stOpt回归模型参数的问题 已有1人参与

用NRTL方程关联三元液液相平衡数据时，利用1stOpt和实验测的相平衡数据来回归模型的最佳参数，总是出现如下错误提示：在”DatSet“中有错误数据。但是如果只选用前五个实验点来关联就可以运行，只要加入了后面五个实验点就出现错误提示。请教诸位大神如何解决这个问题，感激不尽！
以下是1stOpt程序:
Parameters g12,g21,g13,g31,g32,g23;
DataSet;
T,xa,xb,xc,xd,xe,xf,ag,ah,ai=
303.15        0.9797        0                0.0203        0.0210        0                0.9790     0.3  0.3  0.3
303.15        0.9744        0.0052        0.0204        0.0163        0.0036        0.9801        0.3  0.3  0.3
303.15        0.9728        0.0075        0.0197        0.0114        0.0056        0.9831        0.3  0.3  0.3
303.15        0.9697        0.0085        0.0217        0.0124        0.0086        0.9790        0.3  0.3  0.3
303.15        0.9697        0.0107        0.0196        0.0134        0.0108        0.9759        0.3  0.3  0.3
303.15     0.9662        0.0129        0.0208        0.0107        0.0125        0.9768     0.3  0.3  0.3
303.15     0.9649        0.0143        0.0207        0.0126        0.0154        0.9720     0.3  0.3  0.3
303.15     0.9631        0.0163        0.0205        0.0137        0.0178        0.9686     0.3  0.3  0.3
303.15     0.9620        0.0172        0.0208        0.0142        0.0219        0.9639     0.3  0.3  0.3
303.15     0.9613        0.0184        0.0203        0.0126        0.0251        0.9623     0.3  0.3  0.3
EndDataSet;
MinFunction Sum(i=1:10)((xe-xb*exp(((xa*exp(-ag*g12/(8.314*T))+xc*exp(-ai*g32/(8.314*T)))*(g12*xa*exp(-ag*g12/(8.314*T))/(8.314*T)+g32*xc*exp(-ai*g32/(8.314*T))/(8.314*T))/(xb+xa*exp(-ag*g12/(8.314*T))+xc*exp(-ai*g32/(8.314*T)))^2+(g23*xc^2*exp(-ai*g23/(8.314*T))/(8.314*T)+xa*xc*(g23-g13)*exp(-ah*g13/(8.314*T))*exp(-ai*g23/(8.314*T))/(8.314*T))/(xa*exp(-ah*g13/(8.314*T))+xc+xb*exp(-ai*g23/(8.314*T)))^2+(xa^2*g21*exp(-ag*g21/(8.314*T))/(8.314*T)+xa*xc*(g21-g31)*exp(-ag*g21/(8.314*T))*exp(-ah*g31/(8.314*T))/(8.314*T))/(xb*exp(-ag*g21/(8.314*T))+xc*exp(-ah*g31/(8.314*T))+xa)^2)-
                         ((xd*exp(-ag*g12/(8.314*T))+(1-xd-xe)*exp(-ai*g32/(8.314*T)))*(g12*xd*exp(-ag*g12/(8.314*T))/(8.314*T)+g32*(1-xd-xe)*exp(-ai*g32/(8.314*T))/(8.314*T))/(xe+xd*exp(-ag*g12/(8.314*T))+(1-xd-xe)*exp(-ai*g32/(8.314*T)))^2+(g23*(1-xd-xe)^2*exp(-ai*g23/(8.314*T))/(8.314*T)+xd*(1-xd-xe)*(g23-g13)*exp(-ah*g13/(8.314*T))*exp(-ai*g23/(8.314*T))/(8.314*T))/(xd*exp(-ah*g13/(8.314*T))+(1-xd-xe)+xe*exp(-ai*g23/(8.314*T)))^2+(xd^2*g21*exp(-ag*g21/(8.314*T))/(8.314*T)+xd*(1-xd-xe)*(g21-g31)*exp(-ag*g21/(8.314*T))*exp(-ah*g31/(8.314*T))/(8.314*T))/(xe*exp(-ag*g21/(8.314*T))+(1-xd-xe)*exp(-ah*g31/(8.314*T))+xd)^2)))^2+
                       (xd-xa*exp(((xb*exp(-ag*g21/(8.314*T))+xc*exp(-ai*g31/(8.314*T)))*(g21*xb*exp(-ag*g21/(8.314*T))/(8.314*T)+g31*xc*exp(-ah*g31/(8.314*T))/(8.314*T))/(xa+xb*exp(-ag*g21/(8.314*T))+xc*exp(-ah*g31/(8.314*T)))^2+(g12*xb^2*exp(-ag*g12/(8.314*T))/(8.314*T)+xb*xc*(g12-g32)*exp(-ag*g12/(8.314*T))*exp(-ai*g32/(8.314*T))/(8.314*T))/(xa*exp(-ag*g12/(8.314*T))+xb+xc*exp(-ai*g32/(8.314*T)))^2+(xc^2*g13*exp(-ah*g13/(8.314*T))/(8.314*T)+xb*xc*(g13-g23)*exp(-ah*g13/(8.314*T))*exp(-ai*g23/(8.314*T))/(8.314*T))/(xa*exp(-ah*g13/(8.314*T))+xb*exp(-ai*g23/(8.314*T))+xc)^2)-
                         ((xe*exp(-ag*g21/(8.314*T))+(1-xd-xe)*exp(-ai*g31/(8.314*T)))*(g21*xe*exp(-ag*g21/(8.314*T))/(8.314*T)+g31*(1-xd-xe)*exp(-ah*g31/(8.314*T))/(8.314*T))/(xd+xe*exp(-ag*g21/(8.314*T))+(1-xd-xe)*exp(-ah*g31/(8.314*T)))^2+(g12*xe^2*exp(-ag*g12/(8.314*T))/(8.314*T)+xe*(1-xd-xe)*(g12-g32)*exp(-ag*g12/(8.314*T))*exp(-ai*g32/(8.314*T))/(8.314*T))/(xd*exp(-ag*g12/(8.314*T))+xe+(1-xd-xe)*exp(-ai*g32/(8.314*T)))^2+((1-xd-xe)^2*g13*exp(-ah*g13/(8.314*T))/(8.314*T)+xe*(1-xd-xe)*(g13-g23)*exp(-ah*g13/(8.314*T))*exp(-ai*g23/(8.314*T))/(8.314*T))/(xd*exp(-ah*g13/(8.314*T))+xe*exp(-ai*g23/(8.314*T))+(1-xd-xe))^2)))^2+
                       ((xf)-xc*(exp((xa*exp(-ah*g13/(8.314*T))+xb*exp(-ai*g23/(8.314*T)))*(g13*xa*exp(-ah*g13/(8.314*T))/(8.314*T)+g23*xb*exp(-ai*g23/(8.314*T))/(8.314*T))/(xc+xb*exp(-ai*g23/(8.314*T))+xa*exp(-ah*g13/(8.314*T)))^2+(g31*xa^2*exp(-ah*g31/(8.314*T))/(8.314*T)+xb*xa*(g31-g21)*exp(-ag*g21/(8.314*T))*exp(-ah*g31/(8.314*T))/(8.314*T))/(xb*exp(-ag*g21/(8.314*T))+xa+xc*exp(-ah*g31/(8.314*T)))^2+(xb^2*g32*exp(-ai*g32/(8.314*T))/(8.314*T)+xa*xb*(g32-g12)*exp(-ag*g12/(8.314*T))*exp(-ai*g32/(8.314*T))/(8.314*T))/(xa*exp(-ag*g12/(8.314*T))+xc*exp(-ai*g32/(8.314*T))+xb)^2)-
                         ((xd*exp(-ah*g13/(8.314*T))+xe*exp(-ai*g23/(8.314*T)))*(g13*xd*exp(-ah*g13/(8.314*T))/(8.314*T)+g23*xe*exp(-ai*g23/(8.314*T))/(8.314*T))/((1-xd-xe)+xe*exp(-ai*g23/(8.314*T))+xd*exp(-ah*g13/(8.314*T)))^2+(g31*xd^2*exp(-ah*g31/(8.314*T))/(8.314*T)+xe*xd*(g31-g21)*exp(-ag*g21/(8.314*T))*exp(-ah*g31/(8.314*T))/(8.314*T))/(xe*exp(-ag*g21/(8.314*T))+xd+(1-xd-xe)*exp(-ah*g31/(8.314*T)))^2+(xe^2*g32*exp(-ai*g32/(8.314*T))/(8.314*T)+xd*xe*(g32-g12)*exp(-ag*g12/(8.314*T))*exp(-ai*g32/(8.314*T))/(8.314*T))/(xd*exp(-ag*g12/(8.314*T))+(1-xd-xe)*exp(-ai*g32/(8.314*T))+xe)^2)))^2);
回复此楼

» 猜你喜欢

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

最爱明日香

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by dingd at 2014-05-05 16:00:12
改为如下(6.0运行没问题，其它版本不敢保证)：

Constant T=303.15, ag=0.3, ah=0.3, ai=0.3;
Parameters g12,g21,g13,g31,g32,g23;
DataSet;
xa,xb,xc,xd,xe,xf=
        0.9797        0        0.0203        0.0210        0        0.9790
...

可以运行了，谢谢楼上。只是运行的结果和你给的很不一样，请问你用的是什么算法。还有为什么每次重新运行的结果也不一样呢，再次感谢~
赞一下 回复此楼
3楼2014-05-06 10:16:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
改为如下(6.0运行没问题，其它版本不敢保证)：
CODE:
Constant T=303.15, ag=0.3, ah=0.3, ai=0.3;
Parameters g12,g21,g13,g31,g32,g23;
DataSet;
xa,xb,xc,xd,xe,xf=
        0.9797        0        0.0203        0.0210        0        0.9790
        0.9744        0.0052        0.0204        0.0163        0.0036        0.9801
        0.9728        0.0075        0.0197        0.0114        0.0056        0.9831
        0.9697        0.0085        0.0217        0.0124        0.0086        0.9790
        0.9697        0.0107        0.0196        0.0134        0.0108        0.9759
        0.9662        0.0129        0.0208        0.0107        0.0125        0.9768
        0.9649        0.0143        0.0207        0.0126        0.0154        0.9720
        0.9631        0.0163        0.0205        0.0137        0.0178        0.9686
        0.9620        0.0172        0.0208        0.0142        0.0219        0.9639
        0.9613        0.0184        0.0203        0.0126        0.0251        0.9623
EndDataSet;
MinFunction Sum(i=1:10,xa,xb,xc,xd,xe,xf)((xe-xb*exp(((xa*exp(-ag*g12/(8.314*T))+xc*exp(-ai*g32/(8.314*T)))*(g12*xa*exp(-ag*g12/(8.314*T))/(8.314*T)+g32*xc*exp(-ai*g32/(8.314*T))/(8.314*T))/(xb+xa*exp(-ag*g12/(8.314*T))+xc*exp(-ai*g32/(8.314*T)))^2+(g23*xc^2*exp(-ai*g23/(8.314*T))/(8.314*T)+xa*xc*(g23-g13)*exp(-ah*g13/(8.314*T))*exp(-ai*g23/(8.314*T))/(8.314*T))/(xa*exp(-ah*g13/(8.314*T))+xc+xb*exp(-ai*g23/(8.314*T)))^2+(xa^2*g21*exp(-ag*g21/(8.314*T))/(8.314*T)+xa*xc*(g21-g31)*exp(-ag*g21/(8.314*T))*exp(-ah*g31/(8.314*T))/(8.314*T))/(xb*exp(-ag*g21/(8.314*T))+xc*exp(-ah*g31/(8.314*T))+xa)^2)-
                         ((xd*exp(-ag*g12/(8.314*T))+(1-xd-xe)*exp(-ai*g32/(8.314*T)))*(g12*xd*exp(-ag*g12/(8.314*T))/(8.314*T)+g32*(1-xd-xe)*exp(-ai*g32/(8.314*T))/(8.314*T))/(xe+xd*exp(-ag*g12/(8.314*T))+(1-xd-xe)*exp(-ai*g32/(8.314*T)))^2+(g23*(1-xd-xe)^2*exp(-ai*g23/(8.314*T))/(8.314*T)+xd*(1-xd-xe)*(g23-g13)*exp(-ah*g13/(8.314*T))*exp(-ai*g23/(8.314*T))/(8.314*T))/(xd*exp(-ah*g13/(8.314*T))+(1-xd-xe)+xe*exp(-ai*g23/(8.314*T)))^2+(xd^2*g21*exp(-ag*g21/(8.314*T))/(8.314*T)+xd*(1-xd-xe)*(g21-g31)*exp(-ag*g21/(8.314*T))*exp(-ah*g31/(8.314*T))/(8.314*T))/(xe*exp(-ag*g21/(8.314*T))+(1-xd-xe)*exp(-ah*g31/(8.314*T))+xd)^2)))^2+
                       (xd-xa*exp(((xb*exp(-ag*g21/(8.314*T))+xc*exp(-ai*g31/(8.314*T)))*(g21*xb*exp(-ag*g21/(8.314*T))/(8.314*T)+g31*xc*exp(-ah*g31/(8.314*T))/(8.314*T))/(xa+xb*exp(-ag*g21/(8.314*T))+xc*exp(-ah*g31/(8.314*T)))^2+(g12*xb^2*exp(-ag*g12/(8.314*T))/(8.314*T)+xb*xc*(g12-g32)*exp(-ag*g12/(8.314*T))*exp(-ai*g32/(8.314*T))/(8.314*T))/(xa*exp(-ag*g12/(8.314*T))+xb+xc*exp(-ai*g32/(8.314*T)))^2+(xc^2*g13*exp(-ah*g13/(8.314*T))/(8.314*T)+xb*xc*(g13-g23)*exp(-ah*g13/(8.314*T))*exp(-ai*g23/(8.314*T))/(8.314*T))/(xa*exp(-ah*g13/(8.314*T))+xb*exp(-ai*g23/(8.314*T))+xc)^2)-
                         ((xe*exp(-ag*g21/(8.314*T))+(1-xd-xe)*exp(-ai*g31/(8.314*T)))*(g21*xe*exp(-ag*g21/(8.314*T))/(8.314*T)+g31*(1-xd-xe)*exp(-ah*g31/(8.314*T))/(8.314*T))/(xd+xe*exp(-ag*g21/(8.314*T))+(1-xd-xe)*exp(-ah*g31/(8.314*T)))^2+(g12*xe^2*exp(-ag*g12/(8.314*T))/(8.314*T)+xe*(1-xd-xe)*(g12-g32)*exp(-ag*g12/(8.314*T))*exp(-ai*g32/(8.314*T))/(8.314*T))/(xd*exp(-ag*g12/(8.314*T))+xe+(1-xd-xe)*exp(-ai*g32/(8.314*T)))^2+((1-xd-xe)^2*g13*exp(-ah*g13/(8.314*T))/(8.314*T)+xe*(1-xd-xe)*(g13-g23)*exp(-ah*g13/(8.314*T))*exp(-ai*g23/(8.314*T))/(8.314*T))/(xd*exp(-ah*g13/(8.314*T))+xe*exp(-ai*g23/(8.314*T))+(1-xd-xe))^2)))^2+
                       ((xf)-xc*(exp((xa*exp(-ah*g13/(8.314*T))+xb*exp(-ai*g23/(8.314*T)))*(g13*xa*exp(-ah*g13/(8.314*T))/(8.314*T)+g23*xb*exp(-ai*g23/(8.314*T))/(8.314*T))/(xc+xb*exp(-ai*g23/(8.314*T))+xa*exp(-ah*g13/(8.314*T)))^2+(g31*xa^2*exp(-ah*g31/(8.314*T))/(8.314*T)+xb*xa*(g31-g21)*exp(-ag*g21/(8.314*T))*exp(-ah*g31/(8.314*T))/(8.314*T))/(xb*exp(-ag*g21/(8.314*T))+xa+xc*exp(-ah*g31/(8.314*T)))^2+(xb^2*g32*exp(-ai*g32/(8.314*T))/(8.314*T)+xa*xb*(g32-g12)*exp(-ag*g12/(8.314*T))*exp(-ai*g32/(8.314*T))/(8.314*T))/(xa*exp(-ag*g12/(8.314*T))+xc*exp(-ai*g32/(8.314*T))+xb)^2)-
                         ((xd*exp(-ah*g13/(8.314*T))+xe*exp(-ai*g23/(8.314*T)))*(g13*xd*exp(-ah*g13/(8.314*T))/(8.314*T)+g23*xe*exp(-ai*g23/(8.314*T))/(8.314*T))/((1-xd-xe)+xe*exp(-ai*g23/(8.314*T))+xd*exp(-ah*g13/(8.314*T)))^2+(g31*xd^2*exp(-ah*g31/(8.314*T))/(8.314*T)+xe*xd*(g31-g21)*exp(-ag*g21/(8.314*T))*exp(-ah*g31/(8.314*T))/(8.314*T))/(xe*exp(-ag*g21/(8.314*T))+xd+(1-xd-xe)*exp(-ah*g31/(8.314*T)))^2+(xe^2*g32*exp(-ai*g32/(8.314*T))/(8.314*T)+xd*xe*(g32-g12)*exp(-ag*g12/(8.314*T))*exp(-ai*g32/(8.314*T))/(8.314*T))/(xd*exp(-ag*g12/(8.314*T))+(1-xd-xe)*exp(-ai*g32/(8.314*T))+xe)^2)))^2);

目标函数值(最小): 0.00173687399605989
g12: 1120356.41156086
g21: 1438658.65231283
g13: 10744.5036700117
g31: 36020.8511722266
g32: 4738.45016333894
g23: -3612.93060735352
赞一下(2人) 回复此楼
2楼2014-05-05 16:00:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕，总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
信息提示
关闭
请填处理意见
关闭