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

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

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

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

x

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

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

[1stopt] 求高版本1stopt帮忙运行一下，非常感谢！

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