24小时热门版块排行榜    

Znn3bq.jpeg
查看: 660  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂求收留 +27 果然有我 2026-04-10 28/1400 2026-04-12 15:32 by 89436494
[考研] 286求调剂 +26 Faune 2026-04-06 26/1300 2026-04-12 10:24 by lhj2009
[考研] 求调剂 +11 月@163.com 2026-04-07 13/650 2026-04-11 22:55 by BruceLiu320
[考研] 291分调剂 +5 上岸小莹加油 2026-04-09 6/300 2026-04-11 21:06 by 逆水乘风
[考研] 303求调剂 +14 SereinQ 2026-04-10 15/750 2026-04-11 20:43 by 蓝云思雨
[考研] 还有化工二轮调剂的学校吗 5+14 化工人999 2026-04-09 48/2400 2026-04-11 10:27 by 89436494
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 362求调剂 +10 我要考大 2026-04-06 14/700 2026-04-10 17:00 by luoyongfeng
[考研] 297求调剂 +27 GENJIOW 2026-04-07 30/1500 2026-04-09 23:20 by wolf97
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 材料与化工专硕306分找合适调剂 +27 沧海轻舟e 2026-04-06 28/1400 2026-04-08 22:06 by wdyheheeh
[考研] 304求调剂 +16 c297914 2026-04-05 17/850 2026-04-08 13:00 by grayjzr
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 11408 325分 +3 jgtxuxgkx 2026-04-07 3/150 2026-04-07 23:10 by lbsjt
[考研] 318求调剂 +5 李青山山山 2026-04-07 5/250 2026-04-07 18:24 by 蓝云思雨
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +6 白易辰 2026-04-06 7/350 2026-04-06 22:46 by 875465
[考研] 一志愿武汉理工大学080200机械工程308分,求调剂 +4 终不似从前 2026-04-05 4/200 2026-04-06 11:46 by 考研学校招点人
信息提示
请填处理意见