SimWe仿真论坛 » H20：MathCAD Maple » 求助好心人用sharedmodel拟合

ferret 发表于 2008-6-6 17:49

求助好心人用sharedmodel拟合

我用的1.5pro版，出不来结果：（
origin,sigmaplot都试不过了，很难得到结果，想知道1stOpt能帮我找到这些参数吗？请有2.0,2.5版的好心人帮忙运行一下吧。谢谢了！
Title MaxwellG;
Parameters etas[0,],eta1[0,],eta2[0,],eta3[0,],eta4[0,],eta5[0,],lamda1[0,],lamda2[0,],lamda3[0,],lamda4[0,],lamda5[0,];
Variable  x,y1,y2;
SharedModel;
Function y1=(eta1*lamda1*x*x)/(1+lamda1*lamda1*x*x)+(eta2*lamda2*x*x)/(1+lamda2*lamda2*x*x)+(eta3*lamda3*x*x)/(1+lamda3*lamda3*x*x)+(eta4*lamda4*x*x)/(1+lamda4*lamda4*x*x)+(eta5*lamda5*x*x)/(1+lamda5*lamda5*x*x); y2=etas*x+(eta1*x)/(1+lamda1*lamda1*x*x)+(eta2*x)/(1+lamda2*lamda2*x*x)+(eta3*x)/(1+lamda3*lamda3*x*x)+(eta4*x)/(1+lamda4*lamda4*x*x)+(eta5*x)/(1+lamda5*lamda5*x*x);
Data;
0.0501    1.08    70.1
0.0794    3.2    126
0.1    5.57    163
0.126    7.21    215
0.158    12.3    270
0.2    12.3    341
0.251    28.1    438
0.316    39.1    553
0.398    51.3    706
0.501    85    883
0.631    112    1090
0.794    160    1360
1    234    1700
1.26    327    2110
1.58    459    2610
2    629    3200
2.51    871    3910
3.16    1170    4740
3.98    1580    5730
5.01    2100    6900
6.31    2770    8280
7.94    3620    9850
10    4710    11700
12.6    6040    13800
15.8    7670    16100
20    9640    18700
25.1    12000    21500
31.6    14800    24600
39.8    18100    28000
50.1    21900    31600
63.1    26300    35500
79.4    31400    39600
100    37100    43800
122    42400    47800
149    48200    51800
182    54800    55600
222    61800    59700
271    68900    63800
330    76800    67900
403    84600    72300
492    92100    76500
600    99200    80800

[[i] 本帖最后由 ferret 于 2008-6-6 21:02 编辑 [/i]]

20wangz 发表于 2008-6-6 21:30

运行结果如下：

maxwellg

迭代数: 20
计算用时(时:分:秒:微秒): 00:00:50:462
优化算法: 全局麦夸特法(Global Levenberg-Marquardt)
计算结束原因: 达到收敛判定标准
均方差(RMSE): 73.5615662727591
残差平方和(SSE): 454549.538730129
相关系数(R): 0.999997687784651
相关系数之平方(R^2): 0.999995375574648
决定系数(DC): 0.999995375574648
F统计(F-Statistic): 485323.974770657

参数                  最佳估算
--------------------        -------------
etas               65.8811655408653
eta1               538.082393127142
eta2               130.779882136573
eta3               217.698368124466
eta4               499.078610147539
eta5               378.956396056874
lamda1             0.0684855883450606
lamda2             0.00207095580518927
lamda3             0.00630355004658354
lamda4             0.361290542310684
lamda5             0.0192988314116657

====== 输出结果 =====

文件：Data file - 1
No        实测y1        计算y1        实测y2        计算y2
1        1.08        0.567413615170923        70.1        91.6983627507655
2        3.2        1.42459895186636        126        145.305937941404
3        5.57        2.25882985905274        163        182.979947291453
4        7.21        3.58391365757886        215        230.504684549345
5        12.3        5.63012498559776        270        288.948672202868
6        12.3        9.00704847280799        341        365.555519686354
7        28.1        14.1524198806414        438        458.385651861062
8        39.1        22.3460019264833        553        576.317157157644
9        51.3        35.2352712444579        706        724.337439500424
10        85        55.3094610034391        883        908.800771079116
11        112        86.4573708555608        1090        1138.80178088462
12        160        133.841224124405        1360        1421.93374226017
13        234        205.126156191252        1700        1770.19103173049
14        327        309.408883050053        2110        2193.1438164067
15        459        452.774138148033        2610        2687.92030417586
16        629        654.315045938203        3200        3297.54334538935
17        871        905.178576510438        3910        3988.07320907453
18        1170        1222.34773798378        4740        4811.08763817996
19        1580        1615.61904558812        5730        5788.93986946348
20        2100        2108.9617595695        6900        6952.32687006301
21        2770        2747.11452027573        8280        8337.74095609119
22        3620        3578.58384509457        9850        9952.03419927146
23        4710        4660.44063178733        11700        11805.8946326838
24        6040        6026.91582347415        13800        13886.4930262515
25        7670        7655.06923811385        16100        16140.0831211759
26        9640        9675.04958642473        18700        18740.3152817271
27        12000        11980.5140114578        21500        21527.6679064511
28        14800        14755.7900046495        24600        24654.1066637291
29        18100        18056.83400047        28000        28065.4871839791
30        21900        21910.9985261656        31600        31699.7325025807
31        26300        26341.7790685499        35500        35547.9178982199
32        31400        31331.5477410088        39600        39596.9908118758
33        37100        36970.672153926        43800        43890.3393169348
34        42400        42385.6630927789        47800        47748.2051286283
35        48200        48357.5601988247        51800        51722.6917826525
36        54800        54832.0182916092        55600        55756.507635612
37        61800        61706.4417208358        59700        59807.8272953097
38        68900        69006.8312511329        63800        63921.3192790549
39        76800        76543.8285888151        67900        68011.893925865
40        84600        84374.1500732831        72300        72164.1772714461
41        92100        92101.6857641108        76500        76355.8958817583
42        99200        99330.1842137076        80800        80788.5379220239

还有两张图，如果需要，请发信到我的信箱，我发给你。我的信箱是[email]wjh027@gmail.com[/email]。

ferret 发表于 2008-6-6 21:47

非常感谢 20wangz！
我刚发现是我之前犯了个差错,Sigmaplot里也可以算sharedmodel类型的的参数，出来的结果很相似。
不过1stopt好很多，因为不用调初始值。

查看完整版本: 求助好心人用sharedmodel拟合