24小时热门版块排行榜    

查看: 1813  |  回复: 7

jiaqiong424

金虫 (正式写手)

瓜籽皮儿

[求助] 求大神跑下1stOpt,一个sharemodel,谢谢!

请大神代跑下,一个sharemodel,麻烦给看下结果,谢谢!
IntParameter 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
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : jia
  • 2014-06-25 11:08:32, 1.44 K

» 猜你喜欢

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

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

jiaqiong424

金虫 (正式写手)

瓜籽皮儿

没人啊,自己顶下,先用matlab试试
2楼2014-06-25 15:01:20
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jiaqiong424

金虫 (正式写手)

瓜籽皮儿

单独用origin画了下两个函数的拟合图,如果都拟合的好的时候,两个函数的参数值 不太一样。
3楼2014-06-25 15:04:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
jiaqiong424: 金币+40, ★★★很有帮助, 十分感谢 2014-06-25 16:59:48
楼主的公式是否可以改动下,比如:“Sqrt(x-x0-c*313.85)^2”与“abs(x-x0-c*313.85)”效果上有区别吗?
前者有可能出现负数开方导致出错的现象,而后者则不会。改动后计算结果:

均方差(RMSE):1.34753484873093
残差平方和(SSE):94.4242087643035
相关系数(R): 0.934546401041351
相关系数之平方(R^2): 0.873376975699341
确定系数(DC): 0.486778559920511
F统计(F-Statistic): 6.02609707437939

参数                  最佳估算
--------------------        -------------
a        -2.87467666113097
b        31.397228617914
x0        2.2915940326294E-16
c        0.00497831152708185
d        -0.00640191796000243
4楼2014-06-25 16:08:52
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jiaqiong424

金虫 (正式写手)

瓜籽皮儿

引用回帖:
4楼: Originally posted by dingd at 2014-06-25 16:08:52
楼主的公式是否可以改动下,比如:“Sqrt(x-x0-c*313.85)^2”与“abs(x-x0-c*313.85)”效果上有区别吗?
前者有可能出现负数开方导致出错的现象,而后者则不会。改动后计算结果:

均方差(RMSE):1.3475348487309 ...

终于等到大神回复了,十分感谢!!!!!!!!!!!
我修改了下,谢谢您提的问题,之前的好像括号被我弄错了,那个开方还包括后面的d*x1*293部分,能麻烦您再帮我跑一次吗,并且把拟合的图形贴给贴一下,谢谢了!!
IntParameter 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

» 本帖已获得的红花(最新10朵)

5楼2014-06-25 16:58:27
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jiaqiong424

金虫 (正式写手)

瓜籽皮儿

送红花一朵
引用回帖:
5楼: Originally posted by jiaqiong424 at 2014-06-25 16:58:27
终于等到大神回复了,十分感谢!!!!!!!!!!!
我修改了下,谢谢您提的问题,之前的好像括号被我弄错了,那个开方还包括后面的d*x1*293部分,能麻烦您再帮我跑一次吗,并且把拟合的图形贴给贴一下,谢谢了 ...

在线等回复,谢谢
6楼2014-06-25 17:00:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jiaqiong424: 金币+20, ★★★很有帮助, 非常感谢,但是请问这个结果是否还有可能优化 2014-06-26 13:11:31
均方差(RMSE):0.0641941124282293
残差平方和(SSE):0.214285971663304
相关系数(R): 0.999476341637478
相关系数之平方(R^2): 0.998952957493036
确定系数(DC): 0.998952307048273
F统计(F-Statistic): 5493.33798212482

参数                  最佳估算
--------------------        -------------
a        -83.6026915531719
b        -102616.099056213
c        0.326376754921213
d        -6.69667656305467
x0        -661.550785752498
7楼2014-06-25 23:00:07
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jiaqiong424

金虫 (正式写手)

瓜籽皮儿

引用回帖:
7楼: Originally posted by dingd at 2014-06-25 23:00:07
均方差(RMSE):0.0641941124282293
残差平方和(SSE):0.214285971663304
相关系数(R): 0.999476341637478
相关系数之平方(R^2): 0.998952957493036
确定系数(DC): 0.998952307048273
F统计(F-Statistic): 5493.3 ...

想请问下,我想要的参数是大于0的(程序中IntParameter a, b, c, d, x0;
),是否可以模拟出大于0的参数结果?
8楼2014-06-26 13:30:20
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jiaqiong424 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +5 冬十三 2026-03-15 5/250 2026-03-15 18:38 by 无际的草原
[考研] 0703化学调剂 ,六级已过,有科研经历 +4 曦熙兮 2026-03-15 4/200 2026-03-15 18:01 by JourneyLucky
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 材料与化工(0856)304求B区调剂 +7 邱gl 2026-03-10 11/550 2026-03-14 12:18 by 邱gl
[考研] 337一志愿华南理工材料求调剂(有希望2吗?) +3 mysdl 2026-03-09 3/150 2026-03-14 02:53 by JourneyLucky
[考研] 求调剂,药学 +3 归零lbm 2026-03-09 5/250 2026-03-14 02:21 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 304求调剂(085602一志愿985) +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 327分求调剂086 +4 西红柿?小帅 2026-03-09 7/350 2026-03-10 14:47 by ruiyingmiao
信息提示
请填处理意见