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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）
101/1返回列表
查看: 1605  |  回复: 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的回帖

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的回帖

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的回帖

Iamcat

铁杆木虫 (著名写手)

【答案】应助回帖

★
fegg7502: 金币+1, 应助指数+1, 3ks 2013-01-29 09:19:41
第二组数据拟合不好，应该是需要调整参数。你看一下，根据你的要求调整一下我再给你跑一下
赞一下 回复此楼
若以色见我，以音声求我，是人行邪道，不能见如来
6楼2013-01-28 17:24:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

applemyj

木虫 (小有名气)

抽疯狂
引用回帖:
6楼: Originally posted by Iamcat at 2013-01-28 17:24:46
第二组数据拟合不好，应该是需要调整参数。你看一下，根据你的要求调整一下我再给你跑一下

其实，这个就是我不好调的地方，因为方程我只能一个一个调，然后对比参数，经常y2调好y1不行~同学要不你也帮我调调看，a的个数可以增减的，数值也可以在以上的数值中任意取。我明天把我原来单个调试觉得可以的参数发给你，你再看看
赞一下 回复此楼
春风再美比不上你的笑...没见过你的人不会明了
7楼2013-01-28 22:59:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

applemyj

木虫 (小有名气)

抽疯狂
★
fegg7502: 金币+1, 鼓励交流 2013-01-30 08:55:47
引用回帖:
6楼: Originally posted by Iamcat at 2013-01-28 17:24:46
第二组数据拟合不好，应该是需要调整参数。你看一下，根据你的要求调整一下我再给你跑一下

同学，我把第二批的数据调整了一下，能帮我跑一下下面这些数据吗？？？
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
赞一下 回复此楼
春风再美比不上你的笑...没见过你的人不会明了
8楼2013-01-29 16:20:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Iamcat

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fegg7502: 金币+3, 应助指数+1, 3ks 2013-01-30 08:55:59
applemyj: 金币+15, ★★★很有帮助, 我好像之前忘了哈 2013-03-21 09:39:34
内容已删除
赞一下 回复此楼
若以色见我，以音声求我，是人行邪道，不能见如来
9楼2013-01-29 18:12:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

applemyj

木虫 (小有名气)

抽疯狂
★
fegg7502: 金币+1, 鼓励交流 2013-01-30 08:56:41
引用回帖:
9楼: Originally posted by Iamcat at 2013-01-29 18:12:31
实际上这次的拟合结果仍然不好，你可以看一下。另外，a1到a8是10倍的递增，这是你的数据要求的？

迭代数: 75
计算用时(时:分:秒:微秒): 00:00:26:553
优化算法: 通用全局优化法(UGO1)
计算结束原因: 达到收敛 ...

不是的。其实我本来只是知道y1,y2和x的，是要拟合出来一系列的a和g。

看到有一篇文献是先把一组a值假设出来，比如说假设成一系列10倍的递增，然后利用SharedModel拟合出对应的一系列g值，所以我才这样做的，而且也是估计着取了8组值。

因为第一次接触这个软件，，，所以，难道可以不进行假设的前提，直接利用什么函数拟合出来多组a和g嘛？？？指点一下吧，都要硬着头皮现在学matlab了
赞一下 回复此楼
春风再美比不上你的笑...没见过你的人不会明了
10楼2013-01-29 22:04:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 applemyj 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭