ratfitlo.pro

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

  -5.4510e-01   2.7255e-01  -1.3628e-01   6.8138e-02  -3.4069e-02   1.7034e-02
  -1.6915e-01   8.4577e-02  -4.2289e-02   2.1144e-02  -1.0572e-02


 start unimin

    phi=  1.2088e+02   dphi= -1.2088e+02    psi=  1.2088e+02 tau0/2=  5.0000e+03
     fx=  0.0000e+00  dscal=  1.0000e+00    scf=  1.0000e+00   upsi=  1.0074e+02
    sig=  1.0000e+00     fx=  0.0000e+00    psi=  7.8701e+01   upsi=  7.8065e+01


 end unimin

 sig=  1.0000e+00  num. f-evaluations 1
 list of inactive hit constraints
 201   202   203   204   205   206   207   208   209   210   211   212   213  
 214   215   216   217   218   219   220   221   222   223   224   225   226  
 227   228   229   230   231   232   233   234   235   236   237   238   239  
 240   241   242   243   244   245   246   247   248   249   250   251   252  
 253   254   255   256   257   258   259   260   261   262   263   264   265  
 266   267   268   269   270   271   272   273   274   275   276   277   278  
 279   280   281   282   283   284   285   286   287   288   289   290   291  
 292   293   294   295   296   297   298   299   300   301   302   303   304  
 305   306   307   308   309   310   311   312   313   314   315   316   317  
 318   319   320   321   322   323   324   325   326   327   328   329   330  
 331   332   333   334   335   336   337   338   339   340   341   342   343  
 344   345   346   347   348   349   350   351   352   353   354   355   356  
 357   358   359   360   361   362   363   364   365   366   367   368   369  
 370   371   372   373   374   375   376   377   378   379   380   381   382  
 383   384   385   386   387   388   389   390   391   392   393   394   395  
 396   397   398   399   400  
BFGS-update as in Pantoja and Mayne
  tk  =   4.341973e-01
 xsik =   0.000000e+00



================================================================================
             2-th iteration step
   scf=  1.0000e+00 psist=  1.2088e+02   psi=  7.8701e+01  upsi=  7.8065e+01
  fxst=  0.0000e+00    fx=  0.0000e+00
  x=
  -5.4510e-01  -7.2745e-01  -6.3628e-01  -2.6520e-01  -2.8407e-01  -1.8297e-01
   8.3085e-01   1.0846e+00   9.5771e-01   1.0211e+00   9.8943e-01
 valid permutation of x

  1   2   3   4  11   6   7   8   9  10   5 

  multipliers: first estimate
  u =
    1   0.0000e+00    2   0.0000e+00    3   0.0000e+00    4   0.0000e+00
    5   0.0000e+00    6   0.0000e+00    7   0.0000e+00    8   0.0000e+00
    9   0.0000e+00   10   0.0000e+00   11   0.0000e+00   12   0.0000e+00
   13   0.0000e+00   14   0.0000e+00   15   0.0000e+00   16   0.0000e+00
   17   0.0000e+00   18   0.0000e+00   19   0.0000e+00   20   0.0000e+00
   21   0.0000e+00   22   0.0000e+00   23   0.0000e+00   24   0.0000e+00
   25   0.0000e+00   26   0.0000e+00   27   0.0000e+00   28   0.0000e+00
   29   0.0000e+00   30   0.0000e+00   31   0.0000e+00   32   0.0000e+00
   33   0.0000e+00   34   0.0000e+00   35   0.0000e+00   36   0.0000e+00
  201   0.0000e+00  202   0.0000e+00  203   0.0000e+00  204   0.0000e+00
  205   0.0000e+00  206   0.0000e+00  207   0.0000e+00  208   0.0000e+00
  209   0.0000e+00  210   0.0000e+00  211   0.0000e+00  212   0.0000e+00
  213   0.0000e+00  214   0.0000e+00  215   0.0000e+00  216   0.0000e+00
  217   0.0000e+00  218   0.0000e+00  219   0.0000e+00  220   0.0000e+00
  221   0.0000e+00  222   0.0000e+00  223   0.0000e+00  224   0.0000e+00
  225   0.0000e+00  226   0.0000e+00  227   0.0000e+00  228   0.0000e+00
  229   0.0000e+00  230   0.0000e+00  231   0.0000e+00  232   0.0000e+00
  233   0.0000e+00  234   0.0000e+00  235   0.0000e+00  236   0.0000e+00
  237   0.0000e+00  238   0.0000e+00  239   0.0000e+00  240   0.0000e+00
  241   0.0000e+00  242   0.0000e+00  243   0.0000e+00  244   0.0000e+00
  245   0.0000e+00  246   0.0000e+00  247   0.0000e+00  248   0.0000e+00
  249   0.0000e+00  250   0.0000e+00  251   0.0000e+00  252   0.0000e+00
  253   0.0000e+00  254   0.0000e+00  255   0.0000e+00  256   0.0000e+00
  257   0.0000e+00  258   0.0000e+00  259   0.0000e+00  260   0.0000e+00
  261   0.0000e+00  262   0.0000e+00  263   0.0000e+00  264   0.0000e+00
  265   0.0000e+00  266   0.0000e+00  267   0.0000e+00  268   0.0000e+00
  269   0.0000e+00  270   0.0000e+00  271   0.0000e+00  272  -0.0000e+00
  273   0.0000e+00  274   0.0000e+00  275   0.0000e+00  276   0.0000e+00
  277   0.0000e+00  278   0.0000e+00  279   0.0000e+00  280   0.0000e+00
  281   0.0000e+00  282   0.0000e+00  283   0.0000e+00  284   0.0000e+00
  285   0.0000e+00  286   0.0000e+00  287   0.0000e+00  288   0.0000e+00
  289   0.0000e+00  290   0.0000e+00  291   0.0000e+00  292   0.0000e+00
  293   0.0000e+00  294   0.0000e+00  295   0.0000e+00  296   0.0000e+00
  297   0.0000e+00  298   0.0000e+00  299   0.0000e+00  300   0.0000e+00
  301   0.0000e+00  302   0.0000e+00  303   0.0000e+00  304   0.0000e+00
  305   0.0000e+00  306   0.0000e+00  307   0.0000e+00  308   0.0000e+00
  309   0.0000e+00  310   0.0000e+00  311   0.0000e+00  312  -0.0000e+00
  313   0.0000e+00  314   0.0000e+00  315   0.0000e+00  316   0.0000e+00
  317   0.0000e+00  318   0.0000e+00  319   0.0000e+00  320   0.0000e+00
  321   0.0000e+00  322   0.0000e+00  323   0.0000e+00  324   0.0000e+00
  325   0.0000e+00  326   0.0000e+00  327   0.0000e+00  328   0.0000e+00
  329   0.0000e+00  330   0.0000e+00  331   0.0000e+00  332   0.0000e+00
  333   0.0000e+00  334   0.0000e+00  335   0.0000e+00  336   0.0000e+00
  337   0.0000e+00  338   0.0000e+00  339   0.0000e+00  340   0.0000e+00
  341   0.0000e+00  342   0.0000e+00  343   0.0000e+00  344   0.0000e+00
  345   0.0000e+00  346   0.0000e+00  347   0.0000e+00  348   0.0000e+00
  349   0.0000e+00  350   0.0000e+00  351   0.0000e+00  352   0.0000e+00
  353   0.0000e+00  354   0.0000e+00  355   0.0000e+00  356   0.0000e+00
  357   0.0000e+00  358   0.0000e+00  359   0.0000e+00  360   0.0000e+00
  361   0.0000e+00  362   0.0000e+00  363   0.0000e+00  364   0.0000e+00
  365   0.0000e+00  366   0.0000e+00  367   0.0000e+00  368   0.0000e+00
  369   0.0000e+00  370   0.0000e+00  371  -0.0000e+00  372   0.0000e+00
  373   0.0000e+00  374   0.0000e+00  375   0.0000e+00  376   0.0000e+00
  377   0.0000e+00  378   0.0000e+00  379   0.0000e+00  380   0.0000e+00
  381   0.0000e+00  382   0.0000e+00  383   0.0000e+00  384   0.0000e+00
  385   0.0000e+00  386   0.0000e+00  387   0.0000e+00  388   0.0000e+00
  389   0.0000e+00  390   0.0000e+00  391   0.0000e+00  392   0.0000e+00
  393   0.0000e+00  394   0.0000e+00  395   0.0000e+00  396   0.0000e+00
  397   0.0000e+00  398  -0.0000e+00  399   0.0000e+00  400   0.0000e+00
   37   0.0000e+00   38   0.0000e+00   39   0.0000e+00   40   0.0000e+00
   41   0.0000e+00   42   0.0000e+00   43   0.0000e+00   44   0.0000e+00
   45   0.0000e+00   46   0.0000e+00   47   0.0000e+00   48   0.0000e+00
   49   0.0000e+00   50   0.0000e+00   51   0.0000e+00   52   0.0000e+00
   53   0.0000e+00   54   0.0000e+00   55   0.0000e+00   56   0.0000e+00
   57   0.0000e+00   58   0.0000e+00   59   0.0000e+00   60   0.0000e+00
   61   0.0000e+00   62   0.0000e+00   63   0.0000e+00   64   0.0000e+00
   65   0.0000e+00


  del=  1.00000e-01  b2n0=  0.00000e+00   b2n=  0.00000e+00   gfn=  0.00000e+00


 values of restrictions
 (   1   -1.7784e-01    1.5336e+00)  (   2   -1.7254e-01    1.5295e+00)  
 (   3   -1.6725e-01    1.5254e+00)  (   4   -1.6199e-01    1.5213e+00)  
 (   5   -1.5674e-01    1.5173e+00)  (   6   -1.5150e-01    1.5134e+00)  
 (   7   -1.4629e-01    1.5094e+00)  (   8   -1.4109e-01    1.5055e+00)  
 (   9   -1.3591e-01    1.5017e+00)  (  10   -1.3075e-01    1.4978e+00)  
 (  11   -1.2560e-01    1.4940e+00)  (  12   -1.2047e-01    1.4902e+00)  
 (  13   -1.1535e-01    1.4865e+00)  (  14   -1.1025e-01    1.4827e+00)  
 (  15   -1.0517e-01    1.4790e+00)  (  16   -1.0010e-01    1.4754e+00)  
 (  17   -9.5048e-02    1.4717e+00)  (  18   -9.0012e-02    1.4681e+00)  
 (  19   -8.4991e-02    1.4645e+00)  (  20   -7.9985e-02    1.4610e+00)  
 (  21   -7.4995e-02    1.4574e+00)  (  22   -7.0019e-02    1.4539e+00)  
 (  23   -6.5058e-02    1.4504e+00)  (  24   -6.0112e-02    1.4469e+00)  
 (  25   -5.5181e-02    1.4435e+00)  (  26   -5.0264e-02    1.4401e+00)  
 (  27   -4.5362e-02    1.4367e+00)  (  28   -4.0473e-02    1.4333e+00)  
 (  29   -3.5599e-02    1.4299e+00)  (  30   -3.0739e-02    1.4266e+00)  
 (  31   -2.5893e-02    1.4233e+00)  (  32   -2.1061e-02    1.4200e+00)  
 (  33   -1.6242e-02    1.4167e+00)  (  34   -1.1437e-02    1.4135e+00)  
 (  35   -6.6456e-03    1.4103e+00)  (  36   -1.8675e-03    1.4070e+00)  
 ( 201   -4.1531e-01    1.5336e+00)  ( 202   -3.1562e-01    1.5295e+00)  
 ( 203   -3.1594e-01    1.5254e+00)  ( 204   -3.1627e-01    1.5213e+00)  
 ( 205   -3.1661e-01    1.5173e+00)  ( 206   -3.1695e-01    1.5134e+00)  
 ( 207   -3.1730e-01    1.5094e+00)  ( 208   -3.1765e-01    1.5055e+00)  
 ( 209   -3.1802e-01    1.5017e+00)  ( 210   -3.1839e-01    1.4978e+00)  
 ( 211   -3.1876e-01    1.4940e+00)  ( 212   -3.1914e-01    1.4902e+00)  
 ( 213   -3.1953e-01    1.4865e+00)  ( 214   -3.1992e-01    1.4827e+00)  
 ( 215   -3.2032e-01    1.4790e+00)  ( 216   -3.2073e-01    1.4754e+00)  
 ( 217   -3.2114e-01    1.4717e+00)  ( 218   -3.2156e-01    1.4681e+00)  
 ( 219   -3.2198e-01    1.4645e+00)  ( 220   -3.2241e-01    1.4610e+00)  
 ( 221   -3.2284e-01    1.4574e+00)  ( 222   -3.2328e-01    1.4539e+00)  
 ( 223   -3.2373e-01    1.4504e+00)  ( 224   -3.2418e-01    1.4469e+00)  
 ( 225   -3.2464e-01    1.4435e+00)  ( 226   -3.2510e-01    1.4401e+00)  
 ( 227   -3.2557e-01    1.4367e+00)  ( 228   -3.2604e-01    1.4333e+00)  
 ( 229   -3.2652e-01    1.4299e+00)  ( 230   -3.2700e-01    1.4266e+00)  
 ( 231   -3.2749e-01    1.4233e+00)  ( 232   -3.2799e-01    1.4200e+00)  
 ( 233   -3.2849e-01    1.4167e+00)  ( 234   -3.2899e-01    1.4135e+00)  
 ( 235   -3.2950e-01    1.4103e+00)  ( 236   -3.3001e-01    1.4070e+00)  
 ( 237   -3.3053e-01    1.4038e+00)  ( 238   -3.3105e-01    1.4007e+00)  
 ( 239   -3.3158e-01    1.3975e+00)  ( 240   -3.3211e-01    1.3944e+00)  
 ( 241   -3.3265e-01    1.3912e+00)  ( 242   -3.3319e-01    1.3881e+00)  
 ( 243   -3.3374e-01    1.3850e+00)  ( 244   -3.3429e-01    1.3820e+00)  
 ( 245   -3.3485e-01    1.3789e+00)  ( 246   -3.3540e-01    1.3759e+00)  
 ( 247   -3.3597e-01    1.3728e+00)  ( 248   -3.3654e-01    1.3698e+00)  
 ( 249   -3.3711e-01    1.3668e+00)  ( 250   -3.3768e-01    1.3638e+00)  
 ( 251   -3.3826e-01    1.3609e+00)  ( 252   -3.3885e-01    1.3579e+00)  
 ( 253   -3.3944e-01    1.3550e+00)  ( 254   -3.4003e-01    1.3520e+00)  
 ( 255   -3.4063e-01    1.3491e+00)  ( 256   -3.4123e-01    1.3462e+00)  
 ( 257   -3.4183e-01    1.3433e+00)  ( 258   -3.4244e-01    1.3405e+00)  
 ( 259   -3.4305e-01    1.3376e+00)  ( 260   -3.4366e-01    1.3347e+00)  
 ( 261   -3.4428e-01    1.3319e+00)  ( 262   -3.4490e-01    1.3291e+00)  
 ( 263   -3.4553e-01    1.3263e+00)  ( 264   -3.4616e-01    1.3235e+00)  
 ( 265   -3.4679e-01    1.3207e+00)  ( 266   -3.4742e-01    1.3179e+00)  
 ( 267   -3.4806e-01    1.3151e+00)  ( 268   -3.4870e-01    1.3123e+00)  
 ( 269   -3.4935e-01    1.3096e+00)  ( 270   -3.5000e-01    1.3069e+00)  
 ( 271   -3.5065e-01    1.3041e+00)  ( 272   -3.5130e-01    1.3014e+00)  
 ( 273   -3.5196e-01    1.2987e+00)  ( 274   -3.5262e-01    1.2960e+00)  
 ( 275   -3.5328e-01    1.2933e+00)  ( 276   -3.5395e-01    1.2906e+00)  
 ( 277   -3.5462e-01    1.2879e+00)  ( 278   -3.5529e-01    1.2853e+00)  
 ( 279   -3.5596e-01    1.2826e+00)  ( 280   -3.5664e-01    1.2800e+00)  
 ( 281   -3.5732e-01    1.2773e+00)  ( 282   -3.5800e-01    1.2747e+00)  
 ( 283   -3.5868e-01    1.2721e+00)  ( 284   -3.5937e-01    1.2695e+00)  
 ( 285   -3.6006e-01    1.2669e+00)  ( 286   -3.6075e-01    1.2643e+00)  
 ( 287   -3.6145e-01    1.2617e+00)  ( 288   -3.6214e-01    1.2591e+00)  
 ( 289   -3.6284e-01    1.2565e+00)  ( 290   -3.6354e-01    1.2539e+00)  
 ( 291   -3.6424e-01    1.2514e+00)  ( 292   -3.6495e-01    1.2488e+00)  
 ( 293   -3.6566e-01    1.2463e+00)  ( 294   -3.6636e-01    1.2437e+00)  
 ( 295   -3.6708e-01    1.2412e+00)  ( 296   -3.6779e-01    1.2387e+00)  
 ( 297   -3.6850e-01    1.2362e+00)  ( 298   -3.6922e-01    1.2337e+00)  
 ( 299   -3.6994e-01    1.2311e+00)  ( 300   -3.7066e-01    1.2286e+00)  
 ( 301   -3.7138e-01    1.2262e+00)  ( 302   -3.7210e-01    1.2237e+00)