应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求大神跑下1stOpt，一个sharemodel，谢谢！
81/1返回列表
查看: 1767  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-06-25 11:08:11
已阅   关注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朵）

jiaqiong424
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 的主题更新
81/1返回列表
信息提示
关闭
请填处理意见
关闭