testxxxx.pro

sqp程序包。用sqp算法实现非线性约束的优化求解

  -3.8305e-03   3.8361e-03  -4.6390e-03


 start unimin

    phi=  4.4005e-04   dphi= -8.0864e-07    psi=  1.1796e-16 tau0/2=  5.0000e-01
     fx=  4.4005e-04  dscal=  1.0000e+00    scf=  1.0000e+00   upsi=  1.1796e-16
    sig=  1.0000e+00     fx=  4.3937e-04    psi=  2.2898e-16   upsi=  2.2898e-16
    sig=  2.0000e+00     fx=  4.3897e-04    psi=  2.2204e-16   upsi=  2.2204e-16
    sig=  4.0000e+00     fx=  4.3898e-04    psi=  3.6082e-16   upsi=  3.6082e-16


 end unimin

 sig=  2.0000e+00  num. f-evaluations 3
 list of inactive hit constraints
none
BFGS-update as in Pantoja and Mayne
  tk  =   2.193202e-04
 xsik =   0.000000e+00



================================================================================
             5-th iteration step
   scf=  1.0000e+00 psist=  1.0000e-00   psi=  2.2204e-16  upsi=  2.2204e-16
  fxst=  1.5154e-15    fx=  4.3897e-04
  x=
   2.4576e-01   8.9207e-02   8.0012e-02   3.1515e-02   1.9064e-01   8.3766e-02
   6.2624e-02   1.6426e-01   5.2222e-02
 valid permutation of x

  9   2   3   4   5   6   7   8   1 

         quasi-Newton-matrix full update

row/column        1              2              3              4      
      1     6.514574e-01  -5.150521e-03   1.004969e-03   3.130357e-02 
      2     0.000000e+00   9.138795e-01  -1.521636e-01  -1.544834e-03 
      3     0.000000e+00   0.000000e+00   9.165990e-01  -1.353084e-01 
      4     0.000000e+00   0.000000e+00   0.000000e+00   4.536531e-01 
      5     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      6     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      7     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      8     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      9     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 

row/column        5              6              7              8      
      1    -5.661833e-01   1.503683e-02   1.322174e-01  -4.319814e-01 
      2    -7.591856e-02  -1.776131e-01  -1.995842e-01  -9.320904e-02 
      3    -6.699217e-02  -1.838898e-01  -2.078852e-01  -9.218905e-02 
      4     9.414972e-02  -1.222751e-01  -1.526177e-01  -1.195199e-02 
      5     6.310726e-01  -9.844229e-02   7.064668e-02  -7.421376e-01 
      6     0.000000e+00   8.651274e-01  -3.132077e-01  -2.051329e-01 
      7     0.000000e+00   0.000000e+00   7.445617e-01  -5.830602e-02 
      8     0.000000e+00   0.000000e+00   0.000000e+00   1.476141e-01 
      9     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 

row/column        9      
      1     1.942527e-01 
      2    -2.140850e-01 
      3    -2.214208e-01 
      4    -1.426742e-01 
      5     1.536949e-01 
      6    -3.254034e-01 
      7    -6.564027e-01 
      8     1.009825e-01 
      9     1.076767e-01 

  multipliers: first estimate
  u =
    1   8.7848e-04


  del=  1.19048e-03  b2n0=  1.36736e-03   b2n=  5.92458e-04   gfn=  2.67945e-03


 values of restrictions
 (   1    2.2204e-16    3.0000e+00)  
 

   diag[r]=
  -1.0000e+00


 gradient of restriction nr.   1
   1.0000e+00    1.0000e+00    1.0000e+00    1.0000e+00    1.0000e+00  
   1.0000e+00    1.0000e+00    1.0000e+00    1.0000e+00  
 
  multipliers: second estimate
  u =
    1   8.7848e-04
  condition number of r     1.000000000000000e+00
  condition number of a     7.246277668182425e+01


  current scaling,  scf =   1.0000e+00 clow =            1 eta =   1.7133e-20


  scalres=
   1.0000e+00   1.2000e+00   1.2000e+00   1.0000e+00   1.4000e+00   1.2000e+00
   1.2000e+00   1.2000e+00   1.2000e+00   1.2000e+00   1.2000e+00   1.2000e+00
   1.2000e+00   1.2000e+00   1.0000e+00   1.0000e+00   1.0000e+00   1.0000e+00
   1.0000e+00   1.0000e+00   1.0000e+00   1.0000e+00   1.0000e+00


  phase=  2  scf0=  1.0000e+00
  d =
   1.3754e-02  -3.0931e-03  -3.2466e-03  -4.0279e-03   7.6737e-03  -3.4882e-03
  -5.7370e-03   5.0886e-03  -6.9239e-03


  dd=
   4.9295e-17   1.7089e-17   1.5189e-17   5.1675e-18   3.7956e-17   1.5968e-17
   1.1618e-17   3.2527e-17   9.4792e-18


 start unimin

    phi=  4.3897e-04   dphi= -1.8697e-06    psi=  2.2204e-16 tau0/2=  5.0000e-01
     fx=  4.3897e-04  dscal=  1.0000e+00    scf=  1.0000e+00   upsi=  2.2204e-16
    sig=  1.0000e+00     fx=  4.3724e-04    psi=  2.5674e-16   upsi=  2.5674e-16
    sig=  2.0000e+00     fx=  4.3578e-04    psi=  5.0654e-16   upsi=  5.0654e-16
    sig=  4.0000e+00     fx=  4.3367e-04    psi=  2.3419e-15   upsi=  2.3419e-15
    sig=  8.0000e+00     fx=  4.3598e-04    psi=  3.8773e-03   upsi=  3.8773e-03


 end unimin

 sig=  4.0000e+00  num. f-evaluations 4
 list of inactive hit constraints
none
BFGS-update as in Pantoja and Mayne
  tk  =   4.033173e-04
 xsik =   0.000000e+00



================================================================================
             6-th iteration step
   scf=  1.0000e+00 psist=  1.0000e-00   psi=  2.3419e-15  upsi=  2.3419e-15
  fxst=  1.5154e-15    fx=  4.3367e-04
  x=
   3.0078e-01   7.6834e-02   6.7026e-02   1.5403e-02   2.2133e-01   6.9813e-02
   3.9676e-02   1.8461e-01   2.4527e-02
 valid permutation of x

  9   2   3   4   5   6   7   8   1 

         quasi-Newton-matrix full update

row/column        1              2              3              4      
      1     6.476311e-01   4.813609e-03   2.696282e-03  -1.038047e-02 
      2     0.000000e+00   9.092624e-01  -1.545685e-01   3.686758e-02 
      3     0.000000e+00   0.000000e+00   9.162038e-01  -1.206537e-01 
      4     0.000000e+00   0.000000e+00   0.000000e+00   2.649865e-01 
      5     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      6     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      7     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      8     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      9     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 

row/column        5              6              7              8      
      1    -5.693311e-01   2.360137e-02   1.435724e-01  -4.341768e-01 
      2    -7.025364e-02  -1.865911e-01  -2.119186e-01  -8.926953e-02 
      3    -6.531163e-02  -1.872898e-01  -2.124701e-01  -9.105034e-02 
      4     8.786577e-02  -5.742864e-02  -5.843072e-02  -6.495910e-02 
      5     6.299883e-01  -1.005611e-01   6.724047e-02  -7.386552e-01 
      6     0.000000e+00   8.653064e-01  -3.125718e-01  -2.053169e-01 
      7     0.000000e+00   0.000000e+00   7.453536e-01  -5.858391e-02 
      8     0.000000e+00   0.000000e+00   0.000000e+00   1.475717e-01 
      9     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 

row/column        9      
      1     2.079024e-01 
      2    -2.289682e-01 
      3    -2.269301e-01 
      4    -3.677145e-03 
      5     1.464112e-01 
      6    -3.260102e-01 
      7    -6.566579e-01 
      8     1.037954e-01 
      9     7.123440e-02 

  multipliers: first estimate
  u =
    1   8.6746e-04


  del=  2.59921e-03  b2n0=  3.58830e-03   b2n=  7.34633e-04   gfn=  2.64264e-03


 values of restrictions
 (   1    2.3419e-15    3.0000e+00)  
 

   diag[r]=
  -1.0000e+00


 gradient of restriction nr.   1
   1.0000e+00    1.0000e+00    1.0000e+00    1.0000e+00    1.0000e+00  
   1.0000e+00    1.0000e+00    1.0000e+00    1.0000e+00  
 
  multipliers: second estimate
  u =
    1   8.6746e-04
  condition number of r     1.000000000000000e+00
  condition number of a     1.654263398044409e+02
additional increase of eta due to large clow


  current scaling,  scf =   1.0000e+00 clow =            2 eta =   2.2273e-20


  scalres=
   5.5087e-01   1.2000e+00   1.2000e+00   1.0000e+00   1.4000e+00   1.2000e+00
   1.2000e+00   1.2000e+00   1.2000e+00   1.2000e+00   1.2000e+00   1.2000e+00
   1.2000e+00   1.2000e+00   1.0000e+00   1.0000e+00   1.0000e+00   1.0000e+00
   1.0000e+00   1.0000e+00   1.0000e+00   1.0000e+00   1.0000e+00


  phase=  2  scf0=  1.0000e+00
  d =
   6.3409e-02  -1.6716e-02  -1.5159e-02  -6.8013e-03   3.4085e-02  -1.8074e-02
  -2.8698e-02   2.2482e-02  -3.4527e-02


  dd=
  -1.9749e-17  -4.1118e-18  -3.4816e-18  -1.6751e-19  -1.4192e-17  -3.6339e-18
  -1.5313e-18  -1.1644e-17  -4.6902e-19


 start unimin

    phi=  4.3367e-04   dphi= -1.2876e-05    psi=  1.2901e-15 tau0/2=  5.0000e-01
     fx=  4.3367e-04  dscal=  1.0000e+00    scf=  1.0000e+00   upsi=  2.3419e-15
    sig=  1.0000e+00     fx=  4.3219e-04    psi=  6.7672e-03   upsi=  1.1049e-02
    sig=  5.0000e-01     fx=  4.2755e-04    psi=  6.7657e-16   upsi=  1.2282e-15


 end unimin

 sig=  5.0000e-01  num. f-evaluations 2
 list of inactive hit constraints
   2  
BFGS-update as in Pantoja and Mayne
  tk  =   8.585775e-03
 xsik =   0.000000e+00



================================================================================
             7-th iteration step
   scf=  1.0000e+00 psist=  5.5087e-01   psi=  6.7657e-16  upsi=  1.2282e-15
  fxst=  1.5154e-15    fx=  4.2755e-04
  x=
   3.3248e-01   6.8476e-02   5.9446e-02   1.2003e-02   2.3837e-01   6.0776e-02
   2.5328e-02   1.9585e-01   7.2629e-03
 valid permutation of x

  8   2   3   4   5   6   7   9   1 

         quasi-Newton-matrix full update

row/column        1              2              3              4      
      1     6.472574e-01   6.149216e-03   2.807240e-03  -1.654468e-02 
      2     0.000000e+00   9.085171e-01  -1.544371e-01   4.540556e-02 
      3     0.000000e+00   0.000000e+00   9.165875e-01  -1.173645e-01 
      4     0.000000e+00   0.000000e+00   0.000000e+00   1.973352e-01 
      5     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      6     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      7     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      8     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 
      9     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 

row/column        5              6              7              8      
      1    -5.694394e-01   2.471190e-02   1.449934e-01  -4.342573e-01 
      2    -6.969559e-02  -1.879095e-01  -2.137168e-01  -8.879895e-02 
      3    -6.502777e-02  -1.870320e-01  -2.122724e-01  -9.066328e-02 
      4     1.110876e-01  -3.140569e-02  -2.052741e-02  -8.908771e-02 
      5     6.262819e-01  -1.027276e-01   6.421073e-02  -7.362782e-01 
      6     0.000000e+00   8.657012e-01  -3.114033e-01  -2.051051e-01 
      7     0.000000e+00   0.000000e+00   7.465895e-01  -5.781255e-02 
      8     0.000000e+00   0.000000e+00   0.000000e+00   1.478632e-01 
      9     0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00 

row/column        9      
      1     2.096072e-01 
      2    -2.311887e-01 
      3    -2.268264e-01 
      4     6.273427e-02 
      5     1.369227e-01 
      6    -3.265085e-01 
      7    -6.556372e-01 
      8     1.016133e-01 
      9     4.256082e-02 

  multipliers: first estimate
  u =
    1   4.0354e-03    2   3.2002e-03


  del=  4.16260e-03  b2n0=  9.90286e-03   b2n=  8.34206e-04   gfn=  2.61482e-03


 values of restrictions
 (   1    1.2282e-15    3.0000e+00)  (   2    1.1084e-02    2.9291e+00)  
 

   diag[r]=
  -1.0000e+00   8.5771e-02


 gradient of restriction nr.   1
   1.0000e+00    1.0000e+00    1.0000e+00    1.0000e+00    1.0000e+00  
   1.0000e+00    1.0000e+00    1.0000e+00    1.0000e+00  
 

 gradient of restriction nr.   2
  -1.0000e+00   -1.0000e+00   -8.8000e-01   -1.0000e+00   -1.0000e+00  
  -9.6000e-01   -9.4000e-01   -1.0000e+00   -1.0000e+00  
 
  multipliers: second estimate
  u =
    1   4.0354e-03    2   3.2002e-03
  condition number of r     1.165900200484237e+01
  condition number of a     4.637968608929334e+02


  current scaling,  scf =   1.0000e+00 clow =            2 eta =   2.2273e-20


  scalres=
   5.5087e-01   1.2000e+00   1.2000e+00   1.0000e+00   1.4000e+00   1.2000e+00
   1.2000e+00   1.2000e+00   1.2000e+00   1.2000e+00   1.2000e+00   1.2000e+00
   1.2000e+00   1.2000e+00   1.0000e+00   1.0000e+00   1.0000e+00   1.0000e+00
   1.0000e+00   1.0000e+00   1.0000e+00   1.0000e+00   1.0000e+00


  phase=  0  scf0=  1.0000e+00
  d =
   1.5707e-01  -4.7109e-02  -3.7820e-02   1.0124e-02   8.1442e-02  -4.9333e-02
  -7.6209e-02   5.3418e-02  -9.1585e-02


 start unimin

    phi=  4.2755e-04   dphi= -5.1264e-05    psi=  6.7657e-16 tau0/2=  5.0000e-01
     fx=  4.2755e-04  dscal=  1.0000e+00    scf=  1.0000e+00   upsi=  1.2282e-15
    sig=  1.0000e+00     fx=  5.1853e-04    psi=  2.3306e-01   upsi=  2.6735e-01
    sig=  5.0000e-01     fx=  4.5469e-04    psi=  8.2260e-02   upsi=  9.6304e-02
    sig=  1.2500e-01     fx=  4.2545e-04    psi=  2.3055e-03   upsi=  4.1852e-03
    sig=  1.2500e-02     fx=  4.2691e-04    psi=  7.4585e-16   upsi=  1.3540e-15


 end unimin

 sig=  1.2500e-02  num. f-evaluations 4
 list of inactive hit constraints
none
BFGS-update as in Pantoja and Mayne
  tk  =   5.453890e-02
 xsik =   0.000000e+00



================================================================================