24小时热门版块排行榜    

查看: 778  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 17:45 by hZiFeudZyoGR
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 17:33 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 17:25 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:14 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +8 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 17:11 by hZiFeudZyoGR
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 16:34 by xiangy672
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 13:54 by GkCVyODlx4mR
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:03 by l0VvVHGBGRLv
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:39 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:17 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 02:40 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见