1stopt中sharedmodel在低版本运行不了，求助大侠帮运行
Title "818";Parameters A1,A2,A3,B,D1,D2,D3;
Variable x1,y1,x2,y2,x3,y3;
Sharedmodel;
Function y1=2*((1-x1)*(1-x1)-1/(1-x1))*(A1+A2/(1-x1)+2*A3/((1-x1)*(1-x1)*(1-x1)))*((1/(1-x1))^(-B*log(1))+((1/(1-x1))^(-D1))*(D2/(1+D3*log(1))));
y2=2*((1-x2)*(1-x2)-1/(1-x2))*(A1+A2/(1-x2)+2*A3/((1-x2)*(1-x2)*(1-x2)))*((1/(1-x2))^(-B*log(0.828125))+((1/(1-x2))^(-D1))*(D2/(1+D3*log(0.828125))));
y3=2*((1-x3)*(1-x3)-1/(1-x3))*(A1+A2/(1-x3)+2*A3/((1-x3)*(1-x3)*(1-x3)))*((1/(1-x3))^(-B*log(0.4296875))+((1/(1-x3))^(-D1))*(D2/(1+D3*log(0.4296875))));
Date;//x1,y1,x2,y2,x3,y3
5.58E-05 0.00642774 0 0.0132314 0 0.0176991
0.0100369 0.246055 0.00894166 0.212315 0.00528559 0.13766
0.0199765 0.377522 0.0180299 0.300283 0.0105712 0.185841
0.0300602 0.440154 0.0269716 0.335968 0.015942 0.20354
0.0399629 0.47526 0.0360598 0.373282 0.0212276 0.212389
0.0498655 0.510366 0.0450015 0.417137 0.0265132 0.221239
0.0599387 0.545468 0.0539431 0.467527 0.0317988 0.238938
0.0700156 0.590401 0.0630314 0.526083 0.0370844 0.274336
0.079928 0.651073 0.071973 0.592814 0.04237 0.318584
0.0900185 0.731404 0.0810613 0.65137 0.0477408 0.380531
0.100118 0.835334 0.0900029 0.724637 0.0530264 0.446411
0.110045 0.93337 0.0989446 0.797903 0.058312 0.509341
0.119978 1.0491 0.108033 0.884237 0.0635976 0.557522
0.130075 1.14517 0.116974 0.950967 0.0688832 0.588004
0.139994 1.22354 0.126063 1.00952 0.0741688 0.610619
0.149907 1.28618 0.135004 1.05338 0.0795396 0.637168
0.159994 1.35668 0.143946 1.11194 0.0848252 0.672566
0.170074 1.40948 0.153034 1.1705 0.0901108 0.711898
0.17999 1.47998 0.161976 1.22906 0.0953964 0.761062
0.189906 1.55048 0.171064 1.29578 0.100682 0.80531
0.2 1.63868 0.180006 1.34617 0.105968 0.831858
0.209912 1.69935 0.188948 1.39003 0.111338 0.849558
0.22 1.77182 0.198036 1.43388 0.116624 0.867257
0.23008 1.82462 0.206977 1.47773 0.12191 0.898722
0.239992 1.88529 0.216066 1.52158 0.127195 0.938053
0.249902 1.93809 0.225007 1.57361 0.132481 0.973451
0.259988 2.00859 0.233949 1.60929 0.137766 1.00393
0.270072 2.07123 0.243037 1.65968 0.143137 1.02655
0.279994 2.15746 0.251979 1.70353 0.148423 1.0354
0.290078 2.22009 0.261067 1.74738 0.153708 1.04818
0.299994 2.2906 0.270009 1.79287 0.158994 1.0708
0.310075 2.34537 0.27895 1.82856 0.16428 1.10226
0.31999 2.4139 0.288039 1.87894 0.169565 1.12783
0.330074 2.47654 0.29698 1.93097 0.174936 1.14651
0.339993 2.55491 0.306069 1.98952 0.180222 1.15929
0.34991 2.62541 0.31501 2.05625 0.185507 1.17699
0.360003 2.71361 0.323952 2.11481 0.190793 1.19862
0.370094 2.79394 0.33304 2.1652 0.196078 1.23009
0.380013 2.87231 0.341982 2.23193 0.201364 1.25664
0.389933 2.95265 0.35107 2.30356 0.206735 1.28319
0.400026 3.04084 0.360012 2.385 0.21202 1.3058
0.409956 3.14675 0.368953 2.4648 0.217306 1.32743
0.419892 3.26838 0.378042 2.54623 0.222592 1.35398
0.429993 3.37624 0.386983 2.61133 0.227877 1.38053
0.439919 3.47231 0.396072 2.69276 0.233163 1.40315
0.450013 3.56051 0.405013 2.76603 0.238534 1.42478
0.460103 3.64084 0.413955 2.83276 0.243819 1.43363
0.470019 3.71135 0.423043 2.88314 0.249105 1.44248
0.479932 3.77202 0.431985 2.91229 0.25439 1.44248
0.490005 3.80712 0.441073 2.92673 0.259676 1.44248
0.499894 3.80683 0.450015 2.92647 0.264962 1.44248
0.50994 3.77113 0.458956 2.89679 0.270332 1.42478
0.52014 3.69217 0.468045 2.8377 0.275618 1.39823
0.529988 3.58568 0.47577 2.805668 0.280904 1.37168 均方差(RMSE): 0.0539457521231941
残差平方和(SSE): 0.471443355886211
相关系数(R): 0.999859920672014
相关系数之平方(R^2): 0.999719860966246
决定系数(DC): 0.999721350870139
F统计(F-Statistic): 2736.6102196876
参数 最佳估算
-------------------- -------------
a1 -25.0661126112369
a2 20.9932537276705
a3 -2.57913035151389
b -0.185201671868749
d1 -0.231980800917566
d2 -0.813317107378621
d3 0.0435096561920393 回复 2 # dingd 的帖子
结果很棒，在此不胜感激，大侠能否再帮我运行一下，我之前忘了设置一下B和D3的范围，谢谢！
Title "818";
Parameters A1,A2,A3,B[-9,-2],D1,D2,D3[-15,-5];
Variable x1,y1,x2,y2,x3,y3;
Sharedmodel;
Function y1=2*((1-x1)*(1-x1)-1/(1-x1))*(A1+A2/(1-x1)+2*A3/((1-x1)*(1-x1)*(1-x1)))*((1/(1-x1))^(-B*log(1))+((1/(1-x1))^(-D1))*(D2/(1+D3*log(1))));
y2=2*((1-x2)*(1-x2)-1/(1-x2))*(A1+A2/(1-x2)+2*A3/((1-x2)*(1-x2)*(1-x2)))*((1/(1-x2))^(-B*log(0.828125))+((1/(1-x2))^(-D1))*(D2/(1+D3*log(0.828125))));
y3=2*((1-x3)*(1-x3)-1/(1-x3))*(A1+A2/(1-x3)+2*A3/((1-x3)*(1-x3)*(1-x3)))*((1/(1-x3))^(-B*log(0.4296875))+((1/(1-x3))^(-D1))*(D2/(1+D3*log(0.4296875))));
Date;//x1,y1,x2,y2,x3,y3
5.58E-05 0.00642774 0 0.0132314 0 0.0176991
0.0100369 0.246055 0.00894166 0.212315 0.00528559 0.13766
0.0199765 0.377522 0.0180299 0.300283 0.0105712 0.185841
0.0300602 0.440154 0.0269716 0.335968 0.015942 0.20354
0.0399629 0.47526 0.0360598 0.373282 0.0212276 0.212389
0.0498655 0.510366 0.0450015 0.417137 0.0265132 0.221239
0.0599387 0.545468 0.0539431 0.467527 0.0317988 0.238938
0.0700156 0.590401 0.0630314 0.526083 0.0370844 0.274336
0.079928 0.651073 0.071973 0.592814 0.04237 0.318584
0.0900185 0.731404 0.0810613 0.65137 0.0477408 0.380531
0.100118 0.835334 0.0900029 0.724637 0.0530264 0.446411
0.110045 0.93337 0.0989446 0.797903 0.058312 0.509341
0.119978 1.0491 0.108033 0.884237 0.0635976 0.557522
0.130075 1.14517 0.116974 0.950967 0.0688832 0.588004
0.139994 1.22354 0.126063 1.00952 0.0741688 0.610619
0.149907 1.28618 0.135004 1.05338 0.0795396 0.637168
0.159994 1.35668 0.143946 1.11194 0.0848252 0.672566
0.170074 1.40948 0.153034 1.1705 0.0901108 0.711898
0.17999 1.47998 0.161976 1.22906 0.0953964 0.761062
0.189906 1.55048 0.171064 1.29578 0.100682 0.80531
0.2 1.63868 0.180006 1.34617 0.105968 0.831858
0.209912 1.69935 0.188948 1.39003 0.111338 0.849558
0.22 1.77182 0.198036 1.43388 0.116624 0.867257
0.23008 1.82462 0.206977 1.47773 0.12191 0.898722
0.239992 1.88529 0.216066 1.52158 0.127195 0.938053
0.249902 1.93809 0.225007 1.57361 0.132481 0.973451
0.259988 2.00859 0.233949 1.60929 0.137766 1.00393
0.270072 2.07123 0.243037 1.65968 0.143137 1.02655
0.279994 2.15746 0.251979 1.70353 0.148423 1.0354
0.290078 2.22009 0.261067 1.74738 0.153708 1.04818
0.299994 2.2906 0.270009 1.79287 0.158994 1.0708
0.310075 2.34537 0.27895 1.82856 0.16428 1.10226
0.31999 2.4139 0.288039 1.87894 0.169565 1.12783
0.330074 2.47654 0.29698 1.93097 0.174936 1.14651
0.339993 2.55491 0.306069 1.98952 0.180222 1.15929
0.34991 2.62541 0.31501 2.05625 0.185507 1.17699
0.360003 2.71361 0.323952 2.11481 0.190793 1.19862
0.370094 2.79394 0.33304 2.1652 0.196078 1.23009
0.380013 2.87231 0.341982 2.23193 0.201364 1.25664
0.389933 2.95265 0.35107 2.30356 0.206735 1.28319
0.400026 3.04084 0.360012 2.385 0.21202 1.3058
0.409956 3.14675 0.368953 2.4648 0.217306 1.32743
0.419892 3.26838 0.378042 2.54623 0.222592 1.35398
0.429993 3.37624 0.386983 2.61133 0.227877 1.38053
0.439919 3.47231 0.396072 2.69276 0.233163 1.40315
0.450013 3.56051 0.405013 2.76603 0.238534 1.42478
0.460103 3.64084 0.413955 2.83276 0.243819 1.43363
0.470019 3.71135 0.423043 2.88314 0.249105 1.44248
0.479932 3.77202 0.431985 2.91229 0.25439 1.44248
0.490005 3.80712 0.441073 2.92673 0.259676 1.44248
0.499894 3.80683 0.450015 2.92647 0.264962 1.44248
0.50994 3.77113 0.458956 2.89679 0.270332 1.42478
0.52014 3.69217 0.468045 2.8377 0.275618 1.39823
0.529988 3.58568 0.47577 2.805668 0.280904 1.37168 均方差(RMSE): 0.0615732825983757
残差平方和(SSE): 0.614185599050189
相关系数(R): 0.99981391103778
相关系数之平方(R^2): 0.999627856704662
决定系数(DC): 0.999636982299969
F统计(F-Statistic): 3012.78392386749
参数 最佳估算
-------------------- -------------
a1 -2.7262374042727
a2 1.30194358655782
a3 -0.0135132525181409
b -2
d1 -6.0835837994044
d2 0.0330995783682136
d3 -5 回复 4 # dingd 的帖子
太感激了{:3_53:} 这里边这么多数码眼前好多小星星 date应该为data
页:
[1]