24小时热门版块排行榜    

查看: 772  |  回复: 3

小指甲123

新虫 (小有名气)

[求助] 会用1stopt的老师 帮忙看一下对于一环套一环的参数怎样回归 已有1人参与

Parameters g12,g21;
Variable T,x11,x12,y11,y12;
SharedModel;
Constant  alfa=0.2;
                R=8.314;
Conststr  
                tao12=g12/R/T;
                tao21=g21/R/T;
               G12=exp(-alfa*tao12);
              G21=exp(-alfa*tao21);
              x21=1-x11;
              x22=1-x12;
             r11=exp(x21^2*(tao21*(G21/(x11+x21*G21))^2+tao12*G12/(x21+x11*G12)^2));
            r12=exp(x22^2*(tao21*(G21/(x12+x22*G21))^2+tao12*G12/(x22+x12*G12)^2));
            r21=exp(x11^2*(tao12*(G12/(x21+x11*G12))^2+tao21*G21/(x11+x21*G21)^2));
            r22=exp(x12^2*(tao12*(G12/(x22+x12*G12))^2+tao21*G21/(x12+x22*G21)^2));(这些量都是g12 g21 以及变量组合而成的函数)
Function
          y11=(r12*r22)/(r11*r22 - r12*r21) - (r12*r21)/(r11*r22 - r12*r21);
          y12=(r11*r22)/(r11*r22 - r12*r21) - (r11*r21)/(r11*r22 - r12*r21);
Data;
T         x11        x12            y11        y12
303.15      0.7766      0.02522       0.7766      0.02522
313.15      0.7702      0.02548       0.7702      0.02548
323.15      0.7605      0.02683       0.7605      0.02683
333.15      0.7534      0.02740       0.7534      0.02740
其中y11 y12 是因变量
能不能请您帮忙用1stopt跑一下,我不知问题出在哪,跑不出来
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
小指甲123: 金币+10, ★★★★★最佳答案 2014-10-22 10:14:20
注意:1stOpt不区分大小写:
CODE:
Parameters g12,g21;
Variable T,x11,x12,y11,y12;
SharedModel;
Constant  alfa=0.2,
                R=8.314;
Conststr
                tao12=g12/R/T,
                tao21=g21/R/T,
               G12_=exp(-alfa*tao12),
              G21_=exp(-alfa*tao21),
              x21=1-x11,
              x22=1-x12,
             r11=exp(x21^2*(tao21*(G21_/(x11+x21*G21_))^2+tao12*G12_/(x21+x11*G12_)^2)),
            r12=exp(x22^2*(tao21*(G21_/(x12+x22*G21_))^2+tao12*G12_/(x22+x12*G12_)^2)),
            r21=exp(x11^2*(tao12*(G12_/(x21+x11*G12_))^2+tao21*G21_/(x11+x21*G21_)^2)),
            r22=exp(x12^2*(tao12*(G12_/(x22+x12*G12_))^2+tao21*G21_/(x12+x22*G21_)^2));//(这些量都是g12 g21 以及变量组合而成的函数)
Function
          y11=(r12*r22)/(r11*r22 - r12*r21) - (r12*r21)/(r11*r22 - r12*r21);
          y12=(r11*r22)/(r11*r22 - r12*r21) - (r11*r21)/(r11*r22 - r12*r21);
Data;
T         x11        x12            y11        y12
303.15      0.7766      0.02522       0.7766      0.02522
313.15      0.7702      0.02548       0.7702      0.02548
323.15      0.7605      0.02683       0.7605      0.02683
333.15      0.7534      0.02740       0.7534      0.02740

均方差(RMSE):0.00161223431942831
残差平方和(SSE):2.07943960059398E-5
相关系数(R): 0.988044683040149
相关系数之平方(R^2): 0.976232295683909
确定系数(DC): -1.93184123950086
F统计(F-Statistic): 454.59049493873

参数                  最佳估算
--------------------        -------------
g12        385.853819730494
g21        9803.43568253252
2楼2014-10-21 22:30:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小指甲123

新虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2014-10-21 22:30:07
注意:1stOpt不区分大小写:

Parameters g12,g21;
Variable T,x11,x12,y11,y12;
SharedModel;
Constant  alfa=0.2,
                R=8.314;
Conststr
                tao12=g12/R/T,
                 ...

非常感谢你的点拨,但是在我电脑上运行不出来,显示Fatal: Line too long (more than 1023 characters),是不是版本太低的缘故呢
3楼2014-10-22 10:18:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

应该是版本问题吧!
4楼2014-10-22 10:46:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小指甲123 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 05:48 by winsaint
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 02:40 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:14 by ZPa0EcMwuECS
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:12 by ZPa0EcMwuECS
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:03 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 00:52 by ZPa0EcMwuECS
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见