24小时热门版块排行榜    

查看: 2048  |  回复: 1

li343

至尊木虫 (职业作家)


[交流] 1stopt数据拟合问题—麦夸特法拟合结果不能用

CODE:
Parameters a, b, c, d, e;
Variable x1, x2, x3, x4, y;
Function y = a+b*x1+c*x2+d*x3+e*x4;
Data;
0.3431        0.3558         13.6000         0.7200         0.5913
0.2646        0.2274         14.6000         0.7800         0.6521
0.2306        0.1697         17.9000         0.8200         0.702
0.2332        0.1730         22.1000         0.8400         0.7761
0.2747        0.2558         25.5000         0.8400         0.7458
0.3162        0.3487         27.6000         0.8400         0.6718
0.3842        0.4909         28.6000         0.8200         0.5539
0.3633        0.4378         28.4000         0.8200         0.6168
0.3929        0.4678         27.2000         0.7800         0.5804
0.4336        0.5123         24.2000         0.7200         0.5084
0.4592        0.5293         19.6000         0.6600         0.4544
0.4375        0.5051         15.3000         0.6600         0.4801

麦夸特法:
迭代数: 30
计算用时(时:分:秒:微秒): 00:00:00:875
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0174734705867846
残差平方和(SSE): 0.0036638660921667
相关系数(R): 0.984066837057457
相关系数之平方(R^2): 0.968387539796267
决定系数(DC): 0.968387539796266
卡方系数(Chi-Square): 0.00305221476593393
F统计(F-Statistic): 306.330963662834

参数        最佳估算
----------        -------------
a                  111.20677645553
b                  -198.27959261801
c                  74.8838404196295
d                  0.582893244596483
e                  -107.190616584768

粒子群算法:
迭代数: 5001
计算用时(时:分:秒:微秒): 00:00:00:391
优化算法: 粒子群算法CFPSO 1
计算结束原因: 达到最大迭代数
均方差(RMSE): 0.0203361413325662
残差平方和(SSE): 0.00496270373157731
相关系数(R): 0.978444362793864
相关系数之平方(R^2): 0.957353371083091
决定系数(DC): 0.957180947591719
卡方系数(Chi-Square): 0.00381491899822074
F统计(F-Statistic): 224.48512236415

参数        最佳估算
----------        -------------
a                  1.61463343624522
b                  -2.32255537453418
c                  0.37776143548061
d                  0.00858171702132069
e                  -0.688766293128575

用另外一组数据对上述拟合结果进行验证发现麦夸特法误差太大,基本失效,粒子群算法可以通过检验。但问题是粒子群算法的结果不唯一,且从计算得到的相关系数(R)来看精度较麦夸特法低。请问如何才能使麦夸特法的拟合结果得到有效的结果(不限定a,b,c,d,e的范围)?
验证用数据为:
x1        x2        x3        x4        y
0.3546         0.3498         14.2125         0.6825         0.6061
0.2976         0.2506         16.1125         0.7363         0.6641
0.2450         0.1798         18.8625         0.7563         0.7279
0.2585         0.1785         23.0750         0.7913         0.7651
0.3225         0.2740         26.5125         0.7675         0.6866
0.3096         0.2780         28.0000         0.7975         0.6592
0.3834         0.4219         29.4250         0.7500         0.6039
0.3802         0.4043         29.0875         0.7588         0.6315
0.4108         0.4420         27.7250         0.7238         0.5876
0.4506         0.4979         25.1500         0.6588         0.5615
0.4614         0.5157         20.5750         0.6325         0.4860
0.4391         0.4983         15.7625         0.6275         0.4999
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木材

金虫 (小有名气)



小木虫(金币+0.5):给个红包,谢谢回帖
内容已删除
2楼2011-11-18 22:48:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 li343 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金不中,共勉 +5 eulota 2026-08-26 5/250 2026-08-26 23:33 by young90
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见