24小时热门版块排行榜    

查看: 1992  |  回复: 9

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
回复此楼
春风再美比不上你的笑...没见过你的人不会明了
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金中了 +6 laoda193707 2026-08-06 6/300 2026-08-06 23:34 by dragonxp
[基金申请] 听说今天filecode变了 +17 布布和一二 2026-08-06 31/1550 2026-08-06 22:34 by qusibaya
[基金申请] filecode +10 等待解的谜 2026-08-06 14/700 2026-08-06 21:16 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +3 archvillain 2026-08-06 7/350 2026-08-06 21:05 by archvillain
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 6/300 2026-08-06 20:19 by 苏知砚
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 12/600 2026-08-06 15:11 by eddyin
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[教师之家] 咨询面上基金 +4 李长云 2026-07-31 7/350 2026-08-06 11:20 by 李长云
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +22 低垂的野花 2026-07-31 30/1500 2026-08-05 18:03 by 低垂的野花
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[论文投稿] 十年后又回来了,论文投稿求助 +3 哈哈114477 2026-08-01 3/150 2026-08-04 15:40 by tegsgjy20
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见