用1stOpt软件模拟,求解结果为什么有x2=-0.043??? x2不是应该是个变值么????
Parameters a[35,50],b[0.02,0.1],c[0,1],d[0,100],e[0,10],f[0,10];
Variable x1,x2,y;
Function a*x1/(b+x1+c*x2+x1*x2*d+x1*x1/e+x2*x2/f);
Data;
//x1,x2,y;
0.03276245 0.040560018 11.32441
0.037514287 0.038002633 10.7704
0.02777479 0.027442326 10.22368
0.035865345 0.031787511 11.11952
0.033207492 0.041754224 13.16454
0.039172673 0.052949472 14.32423
0.058890928 0.069234367 22.9482
0.06598988 0.076544576 24.10042
0.063845034 0.068018698 22.79831
0.070846079 0.072672711 20.17602
0.060763037 0.076786118 23.25512
0.083675767 0.083558949 22.65642
0.148866584 0.187125943 28.94045
0.153289562 0.178016219 29.99842
0.150439553 0.170146036 28.94679
0.134706648 0.17161249 29.09359
0.14891361 0.186164469 28.96675
0.164153345 0.173520858 26.9409
0.353145247 0.355085849 14.83138
0.313251484 0.334727787 19.4787
0.373761142 0.374815352 10.70323
0.360901466 0.357844418 9.11525
0.353697346 0.41236388 11.29219
0.361937093 0.426641669 9.5724
迭代数: 45
计算用时(时:分:秒:微秒): 00:00:04:57
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 2.21610773764785
残差平方和(SSE): 117.867204116704
相关系数(R): 0.954540798172775
相关系数之平方(R^2): 0.911148135376319
决定系数(DC): 0.911106451585493
卡方系数(Chi-Square): 3.63514471177025
F统计(F-Statistic): 225.603132395452
参数 最佳估算
---------- -------------
a 35.0008402251607
b 0.091965993506086
c 1.41822435291068E-19
d 46.7982996861949
e 0.105217875549044
f 0.0512551223671831
x2 -0.0438194107686669
====== 结果输出 =====
No 实测值y 计算值y
1 11.32441 10.8995663
2 10.7704 12.6999846
3 10.22368 9.0366163
4 11.11952 12.0737297
5 13.16454 11.0674608
6 14.32423 13.3302251
7 22.9482 20.5067568
8 24.10042 22.7597855
9 22.79831 22.1056664
10 20.17602 24.1460145
11 23.25512 21.1246104
12 22.65642 27.1035900
13 28.94045 28.3730903
14 29.99842 27.9885139
15 28.94679 28.2389663
16 29.09359 29.4024200
17 28.96675 28.3691260
18 26.9409 26.9669868
19 14.83138 13.0983733
20 19.4787 14.9597131
21 10.70323 12.2901755
22 9.11525 12.7834027
23 11.29219 13.0754968
24 9.5724 12.7423713 返回小木虫查看更多
把:
Function a*x1/(b+x1+c*x2+x1*x2*d+x1*x1/e+x2*x2/f);
改成:
Function y=a*x1/(b+x1+c*x2+x1*x2*d+x1*x1/e+x2*x2/f);
一样的结果,参数最佳估算还是会出现x2
那就是你的版本问题了,这不一起OK
那你把你的计算结果发我看看呗
Parameters a[35,50],b[0.02,0.1],c[0,1],d[0,100],e[0,10],f[0,10];
第一行最末位分号格式不对,应该是下面的格式。
Parameter a[35,50],b[0.02,0.1],c[0,1],d[0,100],e[0,10],f[0,10];
迭代数: 70
计算用时(时:分:秒:微秒): 00:00:02:510
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 4.75781852697606
残差平方和(SSE): 543.284091255279
相关系数(R): 0.945028578963777
相关系数之平方(R^2): 0.893079015058297
决定系数(DC): 0.590263881876634
卡方系数(Chi-Square): 15.709566322621
F统计(F-Statistic): 183.759421426908
参数 最佳估算
---------- -------------
a 50
b 0.0688742464559206
c 7.22240286699123E-13
d 6.01350655206507E-19
e 0.209102801713676
f 9.99999999979151
====== 结果输出 =====
No 实测值y 计算值y
1 11.32441 15.3189388
2 10.7704 16.5606630
3 10.22368 13.8301898
4 11.11952 16.1566866
5 13.16454 15.4410730
6 14.32423 16.9335632
7 22.9482 20.3310061
8 24.10042 21.1133153
9 22.79831 20.9087167
10 20.17602 21.5663017
11 23.25512 20.5441576
12 22.65642 22.4052627
13 28.94045 22.7468198
14 29.99842 22.6956683
15 28.94679 22.7633210
16 29.09359 22.9635244
17 28.96675 22.7485756
18 26.9409 22.4926189
19 14.83138 17.1256682
20 19.4787 18.1573185
21 10.70323 16.6150917
22 9.11525 16.9361024
23 11.29219 17.0398251
24 9.5724 16.8265762,