这个问题的全局最优差不多就是85.26604,通过reformulate,可以快速得到一个误差还可以接受的全局最优解,如下:
Global optimal solution found.
Objective value: 85.2660408856
Objective bound: 85.2660408856
Infeasibilities: 0.00000000000
Extended solver steps: 1
Total solver iterations: 2174
Elapsed runtime seconds: 0.33
Model Class: NLP
Total variables: 17
Nonlinear variables: 16
Integer variables: 0
Total constraints: 10
Nonlinear constraints: 1
Total nonzeros: 42
Nonlinear nonzeros: 16
Variable Value Reduced Cost
PX( A) 7.25000000337 -0.195495486466E-06
PX( B) 3.25488261883 0.00000000000
PY( A) 7.74999998429 0.267527632047E-06
PY( B) 5.65233243070 0.00000000000
我用另外的工具得到这样的结果:
>>> Objective value is: 85.266040800894487006190
>>> Corresponding solution vector is:
>>> Variable no. Value
>>> 1 7.2499999999996180832795
>>> 2 3.2548829854934648864173
>>> 3 7.7499999999998916422328
>>> 4 5.6523323932727604912429
Global optimal solution found.
Objective value: 85.2660408856
Objective bound: 85.2660408856
Infeasibilities: 0.00000000000
Extended solver steps: 1
Total solver iterations: 2174
Elapsed runtime seconds: 0.33
Model Class: NLP
Total variables: 17
Nonlinear variables: 16
Integer variables: 0
Total constraints: 10
Nonlinear constraints: 1
Total nonzeros: 42
Nonlinear nonzeros: 16
Variable Value Reduced Cost
PX( A) 7.25000000337 -0.195495486466E-06
PX( B) 3.25488261883 0.00000000000
PY( A) 7.74999998429 0.267527632047E-06
PY( B) 5.65233243070 0.00000000000
我用另外的工具得到这样的结果:
>>> Objective value is: 85.266040800894487006190
>>> Corresponding solution vector is:
>>> Variable no. Value
>>> 1 7.2499999999996180832795
>>> 2 3.2548829854934648864173
>>> 3 7.7499999999998916422328
>>> 4 5.6523323932727604912429