当前位置: 首页 > 计算模拟 >用1stOpt软件模拟,求解结果为什么有x2=-0.043??? x2不是应该是个变值么????

用1stOpt软件模拟,求解结果为什么有x2=-0.043??? x2不是应该是个变值么????

作者 11石宝
来源: 小木虫 300 6 举报帖子
+关注

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 返回小木虫查看更多

今日热帖
  • 精华评论
  • dingd

    把:
    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);

  • 11石宝

    引用回帖:
    2楼: Originally posted by dingd at 2017-02-14 17:18:11
    把:
    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

  • dingd

    那就是你的版本问题了,这不一起OK

  • 11石宝

    引用回帖:
    4楼: Originally posted by dingd at 2017-02-15 09:55:09
    那就是你的版本问题了,这不一起OK

    那你把你的计算结果发我看看呗

  • lipenggg

    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];

  • lipenggg

    迭代数: 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,

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓