Thermolucky 发表于 2011-9-29 18:36

求高版本1stopt帮忙运行一下,非常感谢!

Variable x(4),y;
Parameter k(3);
Function y=k1+k2*abs(x1-x2)^-0.45+k3*(x3/x4);
Data;
0.3052        0.2769        0.268434668        0.242913212        0.0016
0.2769        0.32684        0.242913212        0.26004133        -0.0034
0.2615        0.2769        0.255023851        0.242913212        0.0141
0.2769        0.27521        0.242913212        0.252335262        0.0209
0.2769        0.357        0.242913212        0.268491571        0.0164
0.2769        0.3794        0.242913212        0.275783762        -0.0181
0.2769        0.37721        0.242913212        0.266644198        -0.0079
0.3052        0.32684        0.268434668        0.26004133        -0.0014
0.3052        0.2615        0.268434668        0.255023851        -0.0116
0.3052        0.27521        0.268434668        0.252335262        -0.0147
0.3052        0.3794        0.268434668        0.275783762        0.0022
0.3052        0.37721        0.268434668        0.266644198        -0.0007
0.2615        0.32684        0.255023851        0.26004133        -0.0044
0.2615        0.27521        0.255023851        0.252335262        0.0076
0.2615        0.37721        0.255023851        0.266644198        -0.0143
0.32684        0.357        0.26004133        0.268491571        0.0087
0.32684        0.37721        0.26004133        0.266644198        -0.0031
0.32684        0.3776        0.26004133        0.267014537        0.0003
0.357        0.27521        0.268491571        0.252335262        -0.0159
0.263        0.1521        0.258217158        0.276460319        0.199
0.2769        0.1521        0.242913212        0.276460319        0.189
0.3052        0.1521        0.268434668        0.276460319        0.147
0.1521        0.32684        0.276460319        0.26004133        0.166
0.2615        0.1521        0.255023851        0.276460319        0.125
0.1521        0.27521        0.276460319        0.252335262        0.091
0.1521        0.357        0.276460319        0.268491571        0.132
0.1521        0.3794        0.276460319        0.275783762        0.149
0.1521        0.37721        0.276460319        0.266644198        0.149
0.263        0.201        0.258217158        0.273762537        0.198
0.2769        0.201        0.242913212        0.273762537        0.201
0.3052        0.201        0.268434668        0.273762537        0.149
0.32684        0.201        0.26004133        0.273762537        0.163
0.2615        0.201        0.255023851        0.273762537        0.126
0.27521        0.201        0.252335262        0.273762537        0.125
0.357        0.201        0.268491571        0.273762537        0.137
0.263        0.184        0.258217158        0.275858877        0.190
0.2769        0.184        0.242913212        0.275858877        0.202
0.3052        0.184        0.268434668        0.275858877        0.146
0.32684        0.184        0.26004133        0.275858877        0.160
0.2615        0.184        0.255023851        0.275858877        0.127
0.27521        0.184        0.252335262        0.275858877        0.122
0.357        0.184        0.268491571        0.275858877        0.132
0.37721        0.184        0.266644198        0.275858877        0.147
0.184        0.3776        0.275858877        0.267014537        0.153

dingd 发表于 2011-9-30 09:49

公式应写成:y=k1+k2*abs(x1-x2)^(-0.45)+k3*(x3/x4);

低版本也可运行啊,只是结果不好:

均方差(RMSE): 0.0707251457861792
残差平方和(SSE): 0.220090034844957
相关系数(R): 0.462620456244601
相关系数之平方(R^2): 0.214017686535962
决定系数(DC): 0.214017686535962
卡方系数(Chi-Square): 1.75277010876058
F统计(F-Statistic): 26.0820118326212

参数        最佳估算
----------        -------------
k1        0.49740821075128
k2        -0.0113040959580996
k3        -0.380394539349622

Thermolucky 发表于 2011-9-30 13:08

谢谢您的回复!我用的低版本运行时候是错误!我现在把我运行的结果给你看一下!
迭代数: 16
计算用时(时:分:秒:微秒): 00:00:04:625
优化算法: 准牛顿法(BFGS) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.08264348490418
残差平方和(SSE): 0.300517606272727
相关系数(R): 1.23652799814693E-15
相关系数之平方(R^2): 1.52900149020126E-30
决定系数(DC): 1.29302926187346E-15
卡方系数(Chi-Square): 0.570015345271548
F统计(F-Statistic): 6.42180625884528E-29

参数        最佳估算
----------        -------------
k1                 0.895038572458298
k2                 0.153751890756918
x1                 0.0978615316175077
x2                 0.811180888273234
k3                 2.78782801069012
x3                 -0.210231859075656
x4                 0.728042932880228

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

No        实测值y        计算值y
1        0.2769        0.2690127
2        0.32684        0.2690127
3        0.2769        0.2690127
4        0.27521        0.2690127
5        0.357        0.2690127
6        0.3794        0.2690127
7        0.37721        0.2690127
8        0.32684        0.2690127
9        0.2615        0.2690127
10        0.27521        0.2690127
11        0.3794        0.2690127
12        0.37721        0.2690127
13        0.32684        0.2690127
14        0.27521        0.2690127
15        0.37721        0.2690127
16        0.357        0.2690127
17        0.37721        0.2690127
18        0.3776        0.2690127
19        0.27521        0.2690127
20        0.1521        0.2690127
21        0.1521        0.2690127
页: [1]
查看完整版本: 求高版本1stopt帮忙运行一下,非常感谢!