SimWe仿真论坛's Archiver

COMSOL 2008年会圆满结束!

fzl9708 发表于 2007-1-10 14:31

1stopt非线性拟合时多次运行结果相差太大的疑问

本人采用1stopt进行参数拟合,代码如下
Parameters E1[1,50000],eta1[1,1000],E2[1,10000],eta2[1,2000];
Variable x,y(6:12);
//t                  6                7              8              9               10              11      12
function y6=1/E1+x/eta1+(1-exp(-E2/eta2*x))/E2 ;
Data;
0.047        1.53E-04        1.53E-04        1.47E-04        1.53E-04        1.54E-04        1.49E-04        1.52E-04
0.097        2.66E-04        2.66E-04        2.54E-04        2.66E-04        2.67E-04        2.59E-04        2.66E-04
0.149        3.72E-04        3.71E-04        3.58E-04        3.74E-04        3.75E-04        3.62E-04        3.72E-04
0.203        4.77E-04        4.74E-04        4.59E-04        4.79E-04        4.79E-04        4.63E-04        4.77E-04
0.26        5.83E-04        5.77E-04        5.60E-04        5.86E-04        5.83E-04        5.65E-04        5.83E-04
0.32        6.90E-04        6.81E-04        6.62E-04        6.94E-04        6.87E-04        6.69E-04        6.90E-04
0.382        8.00E-04        7.86E-04        7.69E-04        8.05E-04        7.94E-04        7.75E-04        7.99E-04
0.447        9.12E-04        8.93E-04        8.78E-04        9.18E-04        9.02E-04        8.84E-04        9.11E-04
0.516        1.03E-03        1.00E-03        9.91E-04        1.03E-03        1.01E-03        9.97E-04        1.03E-03
0.587        1.14E-03        1.12E-03        1.11E-03        1.27E-03        1.13E-03        1.11E-03        1.14E-03
0.662        1.26E-03        1.23E-03        1.23E-03        1.39E-03        1.24E-03        1.24E-03        1.26E-03
0.741        1.39E-03        1.35E-03        1.36E-03        1.52E-03        1.37E-03        1.36E-03        1.39E-03
0.823        1.52E-03        1.48E-03        1.49E-03        1.64E-03        1.49E-03        1.49E-03        1.52E-03
0.91        1.65E-03        1.61E-03        1.63E-03        1.78E-03        1.62E-03        1.63E-03        1.65E-03
1        1.79E-03        1.74E-03        1.77E-03        1.91E-03        1.76E-03        1.77E-03        1.78E-03

在采用同组数据进行拟和时,多次运行得到的结果有时差别比较大,但是,其相关系数、残差平方和等都表明回归具有很高的质量,这种情况下怎么判断哪一个结果更合理

另外采用不同的回归方法得到的参数结果相差就更大,如以上代码采用简面体爬山法与麦夸特法得到的回归结果差别更大

想请教如何选择合适的回归方法及对回归结果的合理性进行识别?

哪位有关于1stopt中几种回归方法介绍的资料啊,想看一看,希望热心人提供共享

fzl9708 发表于 2007-1-10 14:39

本人采用同一种方法回归曾得到一下两个结果

E1                 21817.99942
eta1               657.4972261
E2                 4359.326634
eta2               1166.419372
残差平方和: 6.56921262993701E-11       
相关系数(R^2): 0.99998241999163       


E1                 12750.91385
eta1               190.709309
E2                 954.2788789
eta2               388.8222842
残差平方和: 7.25872833499016E-11       
相关系数(R^2): 0.999981338901178       

但是可以看到二者差别很大,这时候怎么确定哪一个更合理,很令人头疼

shamohu 发表于 2007-1-10 15:24

用差分或爬山法可得较稳定得结果:

优化算法: 标准差分进化算法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 2.09248811116761E-6
残差平方和(SSE): 6.56775974306669E-11
相关系数(R): 0.999991211713161
相关系数之平方(R^2): 0.999982423503557
决定系数(DC): 0.999982423503557
卡方系数(Chi-Square): 5.15594463345163E-8
F统计(F-Statistic): 208611.91981562

参数        最佳估算
----------        -------------
e1                 21216.5578502987
eta1               657.299336720303
e2                 4386.64987777283
eta2               1182.8148458846

fzl9708 发表于 2007-1-10 16:57

本人经过对方程y6进行多次回归,请楼上的看看这个结果:
均方差(RMSE): 2.13976694540969E-6
残差平方和(SSE): 6.86790387100187E-11
相关系数((R): 0.999991177676713
相关系数之平方(R^2): 0.999982355431258
决定系数(DC): 0.999981620264339
卡方系数(Chi-Square): 5.89339541530299E-8
F统计(F-Statistic): 736757.62842478

参数        最佳估算
----------        -------------
E1                 21188.6705988168
eta1               719.437129927333
E2                 1322.65092706657
eta2               1355.87514879997

这个结果与楼上的结果差别就比较大。

楼上的有没有经历过回归的效果基本相同的情况下,得到的回归参数差别较大,
我关注的是这种情况下怎么判断哪个更合理

shamohu 发表于 2007-1-10 17:01

我得到得参数结果至少整数部分都是一样得,小数部分有微小差别。

fzl9708 发表于 2007-1-10 17:43

[quote]原帖由 [i]shamohu[/i] 于 2007-1-10 17:01 发表
我得到得参数结果至少整数部分都是一样得,小数部分有微小差别。 [/quote]

请问shamohu你用的是哪个版本啊,我用的是1.0的

页: [1]
 

Powered by Discuz! Archiver 6.1.0  © 2001-2007 Comsenz Inc.