24小时热门版块排行榜    

查看: 634  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298求调剂 +5 axyz3 2026-02-28 5/250 2026-03-01 06:45 by 刘兵
[考研] 272求调剂 +4 田智友 2026-02-28 4/200 2026-03-01 06:43 by 刘兵
[考研] 285求调剂 +6 满头大汗的学生 2026-02-28 6/300 2026-03-01 06:29 by Trying]
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 化工专硕348,一志愿985求调剂 +4 弗格个 2026-02-28 6/300 2026-02-28 22:00 by wang_dand
[考研] 292求调剂 +3 yhk_819 2026-02-28 3/150 2026-02-28 21:57 by gaoxiaoniuma
[考研] 295求调剂 +5 19171856320 2026-02-28 5/250 2026-02-28 21:39 by gaoxiaoniuma
[考博] 26申博 +4 想申博! 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 求调剂 +4 repeatt?t 2026-02-28 4/200 2026-02-28 21:16 by gaoxiaoniuma
[考研] 高分子化学与物理调剂 +4 好好好1233 2026-02-28 7/350 2026-02-28 20:42 by 好好好1233
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[考研] 0856材料求调剂 +10 hyf hyf hyf 2026-02-28 11/550 2026-02-28 18:50 by 无际的草原
[考博] 博士自荐 +3 kkluvs 2026-02-28 3/150 2026-02-28 16:59 by StarAura
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见