24小时热门版块排行榜    

查看: 538  |  回复: 2

明雨天地

新虫 (初入文坛)

[求助] 求大虾们利用高版本1stopt帮我跑下程序 已有1人参与

自己的软件版本太低,盼好心人、高版本1stopt出现~~~
非常感谢,一共有两批数据,都是同样的拟合过程,用的是通用全局优化法,拟合出来的g值要>0,所以如果不满足要求的话,可能Constant a的值需要调整,麻烦诸位大虾了
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;
第一组数据
714.954         218.475         0.100
790.545         205.354         0.158
843.715         191.499         0.251
880.668         172.557         0.398
916.241         179.116         0.631
956.146         175.205         1.000
995.362         179.390         1.585
1039.920         184.729         2.512
1082.480         192.036         3.981
1127.350         199.892         6.310
1173.930         208.977         10.000
1222.630         219.333         15.849
1273.370         231.457         25.119
1325.440         245.833         39.811
1376.670         263.674         63.096
1421.090         284.422         100.000
第二组数据
41.624         69.531         0.100
42.107         72.187         0.158
42.964         77.657         0.251
45.814         84.821         0.398
49.226         94.419         0.631
54.841         107.288         1.000
63.455         123.353         1.585
76.759         142.650         2.512
97.070         162.804         3.981
124.786         178.922         6.310
155.685         187.642         10.000
183.444         189.678         15.849
205.171         189.985         25.119
220.368         192.472         39.811
227.868         199.387         63.096
229.082         207.257         100.000
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
第一组:
均方差(RMSE):17.1145813010252
残差平方和(SSE):9373.08457950085
相关系数(R): 0.97371354912747
相关系数之平方(R^2): 0.948118075754413
确定系数(DC): 0.782365064327089
F统计(F-Statistic): 337.676757625548

参数                  最佳估算
--------------------        -------------
g1        16367.42794862
g2        9.61611855404897E-12
g3        243.997605139196
g4        269.059869287778
g5        188.681545840358
g6        298.558728261664
g7        560.35225311584
g8        4.41044543631369E-12

第二组:
均方差(RMSE):43.9464923913789
残差平方和(SSE):61801.4141921767
相关系数(R): 0.949451290628643
相关系数之平方(R^2): 0.901457753276396
确定系数(DC): 0.292328991114505
F统计(F-Statistic): 9.56015330185864

参数                  最佳估算
--------------------        -------------
g1        25105.6865445908
g2        4.29055528438205E-12
g3        4.64851133630867E-13
g4        196.527125439921
g5        54.3067052152756
g6        18.7330481903885
g7        1.72507414162872E-15
g8        8.505768725931E-21
2楼2014-09-04 13:23:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

明雨天地

新虫 (初入文坛)

谢谢大侠哈  还真巧两次都得到了您的帮助  这是上次数据的修改版,数据匹配性挺好的  谢谢了  
3楼2014-09-04 15:34:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 明雨天地 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
请填处理意见