24小时热门版块排行榜    

查看: 830  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

jiaqiong424

金虫 (正式写手)

瓜籽皮儿

[求助] 请用1stopt帮我运行下这个修改后的sharemodel谢谢!!

请用1stopt帮我运行下这个修改后的sharemodel谢谢!!
1.如果方便的话,帮我贴下图,谢谢了!
2.附加问个小问题,x1,y1 是30组数,而x2,y2是22组数,两组数不一样数目,会不会影响1stopt的拟合共享模型的结果?
希望出来的a大概在12左右,x0是大概是个百位数,想先运行下看看吧,谢谢!

Parameters a, b, c, d, x0;
Variable x1,x2,y1,y2;
SharedModel;
Function y1 = a-b/((x1-x0-c*293)+Sqrt ((x1-x0-c*293)^2+d*x1*293));
         y2 = a-b/((x2-x0-c*313.85)+Sqrt((x2-x0-c*313.85)^2+d*x2*313.85));
Data;
// x1 ,   x2 ,   y1   ,    y2
1.495726496        1.49573           6.893169877        8.25219
10.47008547        121.15385        6.711033275        6.73905
15.70512821        131.62393        6.626970228        6.59895
20.94017094        142.09402        6.528896673        6.38879
30.66239316        161.53846        6.34676007        6.02452
31.41025641        180.98291        6.290718039        5.66025
45.61965812        191.45299        5.996497373        5.45009
59.82905983        202.67094        5.6882662        5.23993
61.32478632        210.89744        5.646234676        5.04378
75.53418803        221.36752        5.338003503        4.81961
76.28205128        231.83761        5.295971979        4.60946
90.49145299        241.55983        4.973730298        4.37128
91.23931624        252.02991        4.931698774        4.1331
105.4487179        263.24786        4.567425569        3.83888
110.6837607        271.47436        4.385288967        3.62872
121.1538462        281.19658        4.119089317        3.36252
130.8760684        291.66667        3.796847636        3.05429
135.3632479        301.38889        3.684763573        2.78809
150.3205128        310.36325        3.222416813        2.54991
151.0683761        322.32906        3.124343257        2.18564
165.2777778        330.55556        2.704028021        1.90543
170.5128205        339.52991        2.465849387        1.63923
179.4871795                              2.185639229
189.957265                                1.821366025
190.7051282                               1.737302977
204.1666667                               1.246935201
210.1495726                               1.022767075
219.8717949                               0.65849387
230.3418803                               0.210157618
234.8290598                              -0.014010508
回复此楼

» 本帖@通知

» 猜你喜欢

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

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

jiaqiong424

金虫 (正式写手)

瓜籽皮儿

引用回帖:
2楼: Originally posted by dingd at 2014-06-27 14:35:21
参考下:
均方差(RMSE):1.34870665973305
残差平方和(SSE):94.5885020084306
相关系数(R): 0.939737634027215
相关系数之平方(R^2): 0.883106820807067
确定系数(DC): 0.481504930047079
F统计(F-Statistic): ...

谢谢!dingd 您好,我又追加了金币,可否再帮我再跑一次下面的修改后的程序,谢谢!
Parameters a, b, c, d, x0;
Variable x1,x2,y1,y2;
SharedModel;
Function y1 = a-b/((293-x0-c*x1)+Sqrt ((293-x0-c*x1)^2+d*x1*293));
         y2 = a-b/((313.85-x0-c*x2)+Sqrt((313.85-x0-c*x2)^2+d*x2*313.85));
Data;
// x1 ,   x2 ,   y1   ,    y2
1.495726496        1.49573           6.893169877        8.25219
10.47008547        121.15385        6.711033275        6.73905
15.70512821        131.62393        6.626970228        6.59895
20.94017094        142.09402        6.528896673        6.38879
30.66239316        161.53846        6.34676007        6.02452
31.41025641        180.98291        6.290718039        5.66025
45.61965812        191.45299        5.996497373        5.45009
59.82905983        202.67094        5.6882662        5.23993
61.32478632        210.89744        5.646234676        5.04378
75.53418803        221.36752        5.338003503        4.81961
76.28205128        231.83761        5.295971979        4.60946
90.49145299        241.55983        4.973730298        4.37128
91.23931624        252.02991        4.931698774        4.1331
105.4487179        263.24786        4.567425569        3.83888
110.6837607        271.47436        4.385288967        3.62872
121.1538462        281.19658        4.119089317        3.36252
130.8760684        291.66667        3.796847636        3.05429
135.3632479        301.38889        3.684763573        2.78809
150.3205128        310.36325        3.222416813        2.54991
151.0683761        322.32906        3.124343257        2.18564
165.2777778        330.55556        2.704028021        1.90543
170.5128205        339.52991        2.465849387        1.63923
179.4871795                              2.185639229
189.957265                                1.821366025
190.7051282                               1.737302977
204.1666667                               1.246935201
210.1495726                               1.022767075
219.8717949                               0.65849387
230.3418803                               0.210157618
234.8290598                              -0.014010508
3楼2014-06-27 20:22:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
jiaqiong424: 金币+30, ★★★很有帮助 2014-06-27 20:24:11
参考下:
均方差(RMSE):1.34870665973305
残差平方和(SSE):94.5885020084306
相关系数(R): 0.939737634027215
相关系数之平方(R^2): 0.883106820807067
确定系数(DC): 0.481504930047079
F统计(F-Statistic): 5.76310866738768

参数                  最佳估算
--------------------        -------------
a        -0.429670166828949
b        -2276.35836923505
x0        39.4882598299995
c        -0.598329856341573
d        -0.958630288119618
2楼2014-06-27 14:35:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jiaqiong424

金虫 (正式写手)

瓜籽皮儿

等啊等。。等结果。。。也用matlab做了下。。。。
4楼2014-06-29 12:02:01
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 22:18 by xxniao123
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见