24小时热门版块排行榜    

查看: 829  |  回复: 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

金虫 (正式写手)

瓜籽皮儿

等啊等。。等结果。。。也用matlab做了下。。。。
4楼2014-06-29 12:02:01
已阅   关注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

金虫 (正式写手)

瓜籽皮儿

引用回帖:
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 22:18 by xxniao123
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +8 Kittylucky 2026-08-27 8/400 2026-08-27 18:33 by raiy123
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 基金未中,这种答复是模板吗? +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
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见