24小时热门版块排行榜    

查看: 600  |  回复: 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 的主题更新
信息提示
请填处理意见