fhn 发表于 2010-11-13 10:09

请dingd主任帮忙

本帖最后由 fhn 于 2010-11-13 10:19 编辑

主任,你好,我用的是1stopt1.5,但是sharedmodel好像不好用,请问是版本的问题吗?谢谢。如果是版本问题,您能帮我算下下面的程序吗?谢谢。
Title "eps";
Parameter P1,P2,P3,P4,P5,P6,P7,P8,P9,P10,P11;
Variable x,y1,y2 ;
SharedModel ;
Function y1=9.4-P1*10000*1.6*1.6/(0.22*9.11*8.854*(39.4784*x*x+P2*P2))-P3*(39.4784*x*x-P4*P4)/((39.4784*x*x-P4*P4)*(39.4784*x*x-P4*P4)+39.4784*x*x*P5*P5)-P6*(39.4784*x*x-P7*P7)/((39.4784*x*x-P7*P7)*(39.4784*x*x-P7*P7)+39.4784*x*x*P8*P8)-P9*(39.4784*x*x-P10*P10)/((39.4784*x*x-P10*P10)*(39.4784*x*x-P10*P10)+39.4784*x*x*P11*P11);
         y2=P1*10000*1.6*1.6*P2/(0.22*9.11*8.854*(248.05*x*x*x+6.283*x*P2*P2))+P3*6.283*x*P5/((39.4784*x*x-P4*P4)*(39.4784*x*x-P4*P4)+39.4784*x*x*P5*P5)+P6*6.283*x*P8/((39.4784*x*x-P7*P7)*(39.4784*x*x-P7*P7)+39.4784*x*x*P8*P8)+P9*6.283*x*P11/((39.4784*x*x-P10*P10)*(39.4784*x*x-P10*P10)+39.4784*x*x*P11*P11);
Data;
0.30732    -44.811    7.0532
0.32195    -45.63    8.1728
0.33658    -45.925    9.9178
0.35122    -45.6    11.737
0.36585    -44.456    13.553
0.38049    -42.839    15.406
0.39512    -41.166    16.912
0.40975    -39.403    17.786
0.42439    -37.577    18.25
0.43902    -36.099    18.436
0.45366    -35.11    18.094
0.46829    -34.275    17.236
0.48292    -33.44    16.313
0.49756    -32.838    15.579
0.51219    -32.544    14.928
0.52683    -32.384    14.326
0.54146    -32.278    13.842
0.55609    -32.165    13.403
0.57073    -31.836    12.999
0.58536    -31.252    12.833
0.6    -30.673    12.896
0.61463    -30.16    12.839
0.62927    -29.412    12.548
0.6439    -28.377    12.309
0.65853    -27.414    12.179
0.67317    -26.653    11.829
0.6878    -25.837    11.174
0.70244    -24.947    10.512
0.71707    -24.299    9.9076
0.7317    -23.945    9.1251
0.74634    -23.647    8.1958
0.76097    -23.427    7.3994
0.77561    -23.513    6.7239
0.79024    -23.815    5.9606
0.80487    -24.07    5.2
0.81951    -24.308    4.6718
0.83414    -24.681    4.316
0.84878    -25.092    3.9772
0.86341    -25.399    3.7151
0.87804    -25.663    3.5835
0.89268    -25.891    3.4735
0.90731    -25.979    3.3843
0.92195    -25.993    3.4116
0.93658    -26.076    3.42
0.95121    -26.074    3.2324
0.96585    -25.805    3.0472
0.98048    -25.523    3.1081
0.99512    -25.549    3.1731
1.0098    -25.668    2.8971


dingd 发表于 2010-11-14 18:42

杯具,沦为做题的了?

fhn 发表于 2010-11-14 19:10

回复 2 # dingd 的帖子

dingd主任,不是这样的,我想知道是否是因为版本问题,如果是的话,我正想办法购买4.0版本,让您帮着算算主要是我想看看1stopt和origin算的差别。

dingd 发表于 2010-11-14 22:13

在拟合方面,不论是易用性还是拟合效果方面,Origin都远不如1stOpt。
网上的1.5 PJ版似乎无法运行这种模式。

下面两张图式不考虑参数范围限制的拟合结果,整体均方差为RMSE: 0.401694。用任何其它软件可以说都难以达到这一效果。当然,你的参数范围是有限制的,自己再试试吧,计算量不小。

fhn 发表于 2010-11-14 23:58

回复 4 # dingd 的帖子

感谢dingd主任,可否告诉我拟合出来的参数数值,我好和origin有个比较。

lidaxue1987 发表于 2011-1-17 22:13

"eps"

迭代数: 41
计算用时(时:分:秒:微秒): 00:03:55:96
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.146220094596482
残差平方和(SSE): 1.0476354871264
相关系数(R): 0.999775060161714
相关系数之平方(R^2): 0.999550170921358
决定系数(DC): 0.999550168454576
卡方系数(Chi-Square): -0.0188600726713256
F统计(F-Statistic): 104437.13015431

参数        最佳估算
----------        -------------
P1                 0.182162725931055
P2                 0.000231106083485798
P3                 199.100630560448
P4                 3.38571700571544
P5                 2.15215415935754
P6                 461.10843277154
P7                 5.40874270283985
P8                 3.28572383389455
P9                 562.236402878342
P10              1.74631891174537
P11              2.19802342782235

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

No        实测值y1        计算值y1
1        -44.811        -44.7550437
2        -45.63        -45.7663633
3        -45.925        -45.9568601
4        -45.6        -45.4367139
5        -44.456        -44.3596515
6        -42.839        -42.8928111
7        -41.166        -41.2047689
8        -39.403        -39.4483770
9        -37.577        -37.7562817
10        -36.099        -36.2387942
11        -35.11        -34.9705970
12        -34.275        -33.9938573
13        -33.44        -33.3088290
14        -32.838        -32.8784290
15        -32.544        -32.6351151
16        -32.384        -32.4916616
17        -32.278        -32.3586545
18        -32.165        -32.1588538
19        -31.836        -31.8381059
20        -31.252        -31.3706123
21        -30.673        -30.7556379
22        -30.16        -30.0145929
23        -29.412        -29.1807300
24        -28.377        -28.2964223
25        -27.414        -27.4039832
26        -26.653        -26.5435411
27        -25.837        -25.7523563
28        -24.947        -25.0593190
29        -24.299        -24.4882925
30        -23.945        -24.0541346
31        -23.647        -23.7642735
32        -23.427        -23.6188634
33        -23.513        -23.6098977
34        -23.815        -23.7226527
35        -24.07        -23.9364559
36        -24.308        -24.2264578
37        -24.681        -24.5645024
38        -25.092        -24.9225255
39        -25.399        -25.2726974
40        -25.663        -25.5907467
41        -25.891        -25.8563748
42        -25.979        -26.0535485
43        -25.993        -26.1718886
44        -26.076        -26.2054142
45        -26.074        -26.1526920
46        -25.805        -26.0155871
47        -25.523        -25.7989811
48        -25.549        -25.5091773
49        -25.668        -25.1530234
页: [1]
查看完整版本: 请dingd主任帮忙