1stopt拟合结果的疑惑

li343 · 发表于 2010-4-20 14:46

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

以下三段代码中，K的值都相同，区别近在于代码1中的Kd等于代码2和3中的Hd/H。
疑问1：代码1和2拟合出来的方程系数为什么不一样？
疑问2：为什么将代码2变形为代码3时，拟合精度出现大幅度提高？这不是个别现象，所有如此处理的代码（Function Hd/H = a+b*K+c*K^2+d*K^3;）都出现这样的精度。

代码1：

Parameters a, b;
Variable K, Kd;
Function Kd = a+b*K+c*K^2+d*K^3;
Data;
0.661 0.2904
0.6613 0.3158
0.6593 0.3445
0.6436 0.3801
0.6498 0.3344
0.6336 0.3216
0.609 0.3378
0.6319 0.2943
0.6576 0.2745
0.6936 0.2425
0.6783 0.2627
0.6443 0.3019

复制代码

结果：
迭代数: 23
计算用时(时:分:秒:微秒): 00:00:00:797
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0264070555225598
残差平方和(SSE): 0.00836799097645868
相关系数(R): 0.698912402182278
相关系数之平方(R^2): 0.488478545924203
决定系数(DC): 0.488478545924203
卡方系数(Chi-Square): 0.01295675163964
F统计(F-Statistic): 9.54952215653931
参数最佳估算
---------- -------------
a          40.5545457048544
b          -196.016763020328
c          318.623624273836
d          -172.757524617675

代码2：

Parameters a, b;
Variable K, H, Hd;
Function Hd = H*(a+b*K+c*K^2+d*K^3);
Data;
0.661 9.2609 2.6898
0.6613 12.8403 4.0546
0.6593 17.4916 6.0254
0.6436 21.9374 8.3378
0.6498 25.6834 8.5873
0.6336 26.4722 8.5128
0.609 24.7251 8.3533
0.6319 22.8412 6.7211
0.6576 19.2163 5.2744
0.6936 14.8967 3.6131
0.6783 10.3519 2.7197
0.6443 8.0884 2.4419

复制代码

结果：
迭代数: 37
计算用时(时:分:秒:微秒): 00:00:01:203
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.516997264963694
残差平方和(SSE): 3.20743406375928
相关系数(R): 0.976585316853286
相关系数之平方(R^2): 0.953718881093433
决定系数(DC): 0.952200623736211
卡方系数(Chi-Square): 0.253781282259213
F统计(F-Statistic): 206.070834851427
参数最佳估算
---------- -------------
a          112.208456538171
b          -529.399962761009
c          834.989096449471
d          -438.986066191764

代码3：

Parameters a, b;
Variable K, H, Hd;
Function Hd/H = a+b*K+c*K^2+d*K^3;
Data;
0.661 9.2609 2.6898
0.6613 12.8403 4.0546
0.6593 17.4916 6.0254
0.6436 21.9374 8.3378
0.6498 25.6834 8.5873
0.6336 26.4722 8.5128
0.609 24.7251 8.3533
0.6319 22.8412 6.7211
0.6576 19.2163 5.2744
0.6936 14.8967 3.6131
0.6783 10.3519 2.7197
0.6443 8.0884 2.4419

复制代码

结果：
迭代数: 20
计算用时(时:分:秒:微秒): 00:00:00:859
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0264023110937081
残差平方和(SSE): 0.00836498437306733
相关系数(R): 0.999949966509266
相关系数之平方(R^2): 0.999899935521882
决定系数(DC): 0.999875339281326
卡方系数(Chi-Square): 0.000729076935108954
F统计(F-Statistic): 99925.5634293224
参数最佳估算
---------- -------------
a          40.7897484457353
b          -197.078763281855
c          320.219075739479
d          -173.554993489736

[ 本帖最后由 li343 于 2010-4-20 14:53 编辑 ]

账号		自动登录	找回密码
密码			我要加入

[1stopt] 1stopt拟合结果的疑惑

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。