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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
51/1返回列表
查看: 1613  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

applemyj

木虫 (小有名气)

抽疯狂

[求助] 求大哥大姐帮跑一下1stopt~新手啊，版本低啊

自己的软件版本太低，盼好心人、高版本1stopt出现~~~
非常感谢，俺也知道大家的时间很宝贵，所以尽量多给金币了~~~
一共有两批数据，都是同样的拟合过程

Constant a1=0.0001,a2=0.001,a3=0.01,a4=0.1,a5=1,a6=10,a7=100,a8=1000;
Parameter g(1:8);
Variables y1,y2,x;
SharedModel;
Function y1=g1*a1^2*x^2/(1+a1^2*x^2)+g2*a2^2*x^2/(1+a2^2*x^2)+g3*a3^2*x^2/(1+a3^2*x^2)+g4*a4^2*x^2/(1+a4^2*x^2)+g5*a5^2*x^2/(1+a5^2*x^2)+g6*a6^2*x^2/(1+a6^2*x^2)+g7*a7^2*x^2/(1+a7^2*x^2)+g8*a8^2*x^2/(1+a8^2*x^2);
         y2=g1*a1*x/(1+a1^2*x^2)+g2*a2*x/(1+a2^2*x^2)+g3*a3*x/(1+a3^2*x^2)+g4*a4*x/(1+a4^2*x^2)+g5*a5*x/(1+a5^2*x^2)+g6*a6*x/(1+a6^2*x^2)+g7*a7*x/(1+a7^2*x^2)+g8*a8*x/(1+a8^2*x^2);
Data;
第一组数据：
42.75        9.227        0.06283
45.16        9.471        0.0791
46.97        9.602        0.09958
48.7         9.642        0.1254
50.3                9.692        0.1578
51.83        9.666        0.1987
53.31        9.599        0.2501
54.69        9.584        0.3149
56.06        9.458        0.3964
57.36        9.305        0.4991
58.65        9.237        0.6283
59.89        9.147        0.791
61.02        9.061        0.9958
62.11        8.981        1.254
63.24        9.023        1.578
64.27        9.073        1.987
65.35        9.146        2.501
66.45        9.306        3.149
67.48        9.622        3.964
68.61        9.99               4.991
69.82        10.46        6.283
71.04        10.95        7.91
72.26        11.61        9.958
73.69        12.34        12.54
75.2               13.19        15.78
76.84        14.14        19.87
78.54        15.22        25.01
80.45        16.42        31.49
82.49        17.76        39.64
84.62        19.21        49.91
86.93        21.01        62.83
86.72        20.87        62.83
第二组数据：
3.69         5.22         0.06283
3.856        4.242        0.0791
3.687        4.166        0.09958
3.544        4.172        0.1254
3.426        4.203        0.1578
3.319        4.264        0.1987
3.214        4.348        0.2501
3.13               4.466        0.3149
3.055        4.608        0.3964
2.981        4.742        0.4991
2.887        4.897        0.6283
2.805        5.119        0.791
2.766        5.346        0.9958
2.691        5.602        1.254
2.65               5.876        1.578
2.629        6.192        1.987
2.588        6.546        2.501
2.571        6.909        3.149
2.565        7.35               3.964
2.583        7.832        4.991
2.6           8.381        6.283
2.659        8.999        7.91
2.776        9.695        9.958
2.973        10.49        12.54
3.283        11.35        15.78
3.715        12.31        19.87
4.242        13.42        25.01
4.874        14.71        31.49
5.636        16.25        39.64
6.635        18.02        49.91
8.058        19.71        62.83
8.086        19.55        62.83
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
春风再美比不上你的笑...没见过你的人不会明了
1楼2013-01-27 16:47:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Iamcat

铁杆木虫 (著名写手)
★ ★
fegg7502: 金币+2, 鼓励交流 2013-01-29 09:19:31
上面是第一组数据的拟合，但是y1的拟合结果很好，y2的不算好。

1.jpg



1-y2.jpg
赞一下 回复此楼
若以色见我，以音声求我，是人行邪道，不能见如来
5楼2013-01-28 17:23:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

applemyj

木虫 (小有名气)

抽疯狂
补充：拟合出来的g值要>0，所以如果不满足要求的话，可能Constant a的值需要调整，麻烦诸位大虾了
赞一下 回复此楼
春风再美比不上你的笑...没见过你的人不会明了
2楼2013-01-27 16:50:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Iamcat

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
applemyj: 金币+20, ★★★★★最佳答案 2013-01-28 17:16:05
迭代数: 50
计算用时(时:分:秒:微秒): 00:02:14:58
优化算法: 通用全局优化法(UGO1)
计算结束原因: 达到收敛判定标准
均方差(RMSE): 0.457554430452535
残差平方和(SSE): 13.3987876369116
相关系数(R): 0.9999538427357
相关系数之平方(R^2): 0.999907687601893
决定系数(DC): 0.999907687601894
F统计(F-Statistic): 3562.34348909742

参数                  最佳估算
--------------------        -------------
g1        927.857807878991
g2        9.57785631266135E-13
g3        30.3737389532088
g4        12.5751749127888
g5        11.2164978073129
g6        17.0073506738556
g7        1.69984592611317E-12
g8        38.0722834931112

====== 输出结果 =====

文件：数据文件 - 1
No        目标y1        计算y1        目标y2        计算y2
1        42.75        42.9208655882332        9.227        9.07301175786781
2        45.16        44.6824214427367        9.471        9.76896519890456
3        46.97        46.647739107148        9.602        10.1565588445128
4        48.7        48.6417149047618        9.642        10.1861644022169
5        50.3        50.4807769676214        9.692        9.91884968686892
6        51.83        52.0726793695168        9.666        9.49374507297026
7        53.31        53.4031433386455        9.599        9.06871699358951
8        54.69        54.5459257412626        9.584        8.76090705808797
9        56.06        55.6051504740333        9.458        8.62711010015922
10        57.36        56.6919283086835        9.305        8.65808444172752
11        58.65        57.884623561838        9.237        8.78943502508435
12        59.89        59.2090543844829        9.147        8.92412069327581
13        61.02        60.6209358191493        9.061        8.97205706822433
14        62.11        62.0279852602239        8.981        8.89533200269829
15        63.24        63.3273798375738        9.023        8.73074605029889
16        64.27        64.4760123961251        9.073        8.56849146194841
17        65.35        65.4822607395558        9.146        8.51193568323656
18        66.45        66.4161485371833        9.306        8.63763832848576
19        67.48        67.3696455632238        9.622        8.97648277669657
20        68.61        68.4399021188289        9.99        9.50877077753916
21        69.82        69.6936753378756        10.46        10.1663962825769
22        71.04        71.1462540772053        10.95        10.8557003983832
23        72.26        72.7427110475973        11.61        11.4961422946902
24        73.69        74.3827653652001        12.34        12.0707900142452
25        75.2        75.9629103469361        13.19        12.6445462252856
26        76.84        77.4584805102356        14.14        13.3499721258169
27        78.54        78.9136599671919        15.22        14.3221416933818
28        80.45        80.4572601371014        16.42        15.6622129496805
29        82.49        82.2508432131028        17.76        17.392306553707
30        84.62        84.4617666451246        19.21        19.4489901385217
31        86.93        87.1910850455587        21.01        21.6701667656268
32        86.72        87.1910850455587        20.87        21.6701667656268
赞一下 回复此楼
若以色见我，以音声求我，是人行邪道，不能见如来
3楼2013-01-28 17:11:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

applemyj

木虫 (小有名气)

抽疯狂
引用回帖:
3楼: Originally posted by Iamcat at 2013-01-28 17:11:40
迭代数: 50
计算用时(时:分:秒:微秒): 00:02:14:58
优化算法: 通用全局优化法(UGO1)
计算结束原因: 达到收敛判定标准
均方差(RMSE): 0.457554430452535
残差平方和(SSE): 13.3987876369116
相关系数(R): 0.99 ...

嘻嘻，还有一组数据呢，，谢过好心人~~~
话说大哥你真凑巧，我前一秒刚追加了金币，再刷新的时候你的回复就出来了
赞一下 回复此楼
春风再美比不上你的笑...没见过你的人不会明了
4楼2013-01-28 17:14:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭