ratfitlo.pro

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

 ( 303   -3.7283e-01    1.2212e+00)  ( 304   -3.7356e-01    1.2187e+00)  
 ( 305   -3.7428e-01    1.2162e+00)  ( 306   -3.7501e-01    1.2138e+00)  
 ( 307   -3.7574e-01    1.2113e+00)  ( 308   -3.7648e-01    1.2089e+00)  
 ( 309   -3.7721e-01    1.2064e+00)  ( 310   -3.7795e-01    1.2040e+00)  
 ( 311   -3.7868e-01    1.2015e+00)  ( 312   -3.7942e-01    1.1991e+00)  
 ( 313   -3.8016e-01    1.1967e+00)  ( 314   -3.8090e-01    1.1943e+00)  
 ( 315   -3.8164e-01    1.1919e+00)  ( 316   -3.8238e-01    1.1894e+00)  
 ( 317   -3.8313e-01    1.1870e+00)  ( 318   -3.8387e-01    1.1846e+00)  
 ( 319   -3.8462e-01    1.1823e+00)  ( 320   -3.8536e-01    1.1799e+00)  
 ( 321   -3.8611e-01    1.1775e+00)  ( 322   -3.8686e-01    1.1751e+00)  
 ( 323   -3.8760e-01    1.1727e+00)  ( 324   -3.8835e-01    1.1704e+00)  
 ( 325   -3.8910e-01    1.1680e+00)  ( 326   -3.8985e-01    1.1656e+00)  
 ( 327   -3.9060e-01    1.1633e+00)  ( 328   -3.9136e-01    1.1609e+00)  
 ( 329   -3.9211e-01    1.1586e+00)  ( 330   -3.9286e-01    1.1562e+00)  
 ( 331   -3.9361e-01    1.1539e+00)  ( 332   -3.9437e-01    1.1516e+00)  
 ( 333   -3.9512e-01    1.1493e+00)  ( 334   -3.9588e-01    1.1469e+00)  
 ( 335   -3.9663e-01    1.1446e+00)  ( 336   -3.9739e-01    1.1423e+00)  
 ( 337   -3.9814e-01    1.1400e+00)  ( 338   -3.9890e-01    1.1377e+00)  
 ( 339   -3.9965e-01    1.1354e+00)  ( 340   -4.0041e-01    1.1331e+00)  
 ( 341   -4.0116e-01    1.1308e+00)  ( 342   -4.0192e-01    1.1285e+00)  
 ( 343   -4.0267e-01    1.1262e+00)  ( 344   -4.0343e-01    1.1239e+00)  
 ( 345   -4.0418e-01    1.1217e+00)  ( 346   -4.0494e-01    1.1194e+00)  
 ( 347   -4.0569e-01    1.1171e+00)  ( 348   -4.0645e-01    1.1148e+00)  
 ( 349   -4.0720e-01    1.1126e+00)  ( 350   -4.0796e-01    1.1103e+00)  
 ( 351   -4.0871e-01    1.1081e+00)  ( 352   -4.0946e-01    1.1058e+00)  
 ( 353   -4.1022e-01    1.1036e+00)  ( 354   -4.1097e-01    1.1013e+00)  
 ( 355   -4.1172e-01    1.0991e+00)  ( 356   -4.1247e-01    1.0969e+00)  
 ( 357   -4.1322e-01    1.0946e+00)  ( 358   -4.1397e-01    1.0924e+00)  
 ( 359   -4.1472e-01    1.0902e+00)  ( 360   -4.1547e-01    1.0880e+00)  
 ( 361   -4.1622e-01    1.0857e+00)  ( 362   -4.1697e-01    1.0835e+00)  
 ( 363   -4.1771e-01    1.0813e+00)  ( 364   -4.1846e-01    1.0791e+00)  
 ( 365   -4.1920e-01    1.0769e+00)  ( 366   -4.1995e-01    1.0747e+00)  
 ( 367   -4.2069e-01    1.0725e+00)  ( 368   -4.2143e-01    1.0703e+00)  
 ( 369   -4.2217e-01    1.0682e+00)  ( 370   -4.2291e-01    1.0660e+00)  
 ( 371   -4.2365e-01    1.0638e+00)  ( 372   -4.2439e-01    1.0616e+00)  
 ( 373   -4.2512e-01    1.0594e+00)  ( 374   -4.2586e-01    1.0573e+00)  
 ( 375   -4.2659e-01    1.0551e+00)  ( 376   -4.2732e-01    1.0530e+00)  
 ( 377   -4.2806e-01    1.0508e+00)  ( 378   -4.2879e-01    1.0486e+00)  
 ( 379   -4.2951e-01    1.0465e+00)  ( 380   -4.3024e-01    1.0443e+00)  
 ( 381   -4.3097e-01    1.0422e+00)  ( 382   -4.3169e-01    1.0401e+00)  
 ( 383   -4.3241e-01    1.0379e+00)  ( 384   -4.3313e-01    1.0358e+00)  
 ( 385   -4.3385e-01    1.0336e+00)  ( 386   -4.3457e-01    1.0315e+00)  
 ( 387   -4.3528e-01    1.0294e+00)  ( 388   -4.3600e-01    1.0273e+00)  
 ( 389   -4.3671e-01    1.0252e+00)  ( 390   -4.3742e-01    1.0230e+00)  
 ( 391   -4.3813e-01    1.0209e+00)  ( 392   -4.3884e-01    1.0188e+00)  
 ( 393   -4.3954e-01    1.0167e+00)  ( 394   -4.4024e-01    1.0146e+00)  
 ( 395   -4.4094e-01    1.0125e+00)  ( 396   -4.4164e-01    1.0104e+00)  
 ( 397   -4.4234e-01    1.0083e+00)  ( 398   -4.4303e-01    1.0062e+00)  
 ( 399   -4.4372e-01    1.0042e+00)  ( 400   -4.4441e-01    1.0021e+00)  
 (  37    2.8973e-03    1.4038e+00)  (  38    7.6489e-03    1.4007e+00)  
 (  39    1.2387e-02    1.3975e+00)  (  40    1.7113e-02    1.3944e+00)  
 (  41    2.1826e-02    1.3912e+00)  (  42    2.6525e-02    1.3881e+00)  
 (  43    3.1212e-02    1.3850e+00)  (  44    3.5887e-02    1.3820e+00)  
 (  45    4.0549e-02    1.3789e+00)  (  46    4.5198e-02    1.3759e+00)  
 (  47    4.9835e-02    1.3728e+00)  (  48    5.4460e-02    1.3698e+00)  
 (  49    5.9073e-02    1.3668e+00)  (  50    6.3673e-02    1.3638e+00)  
 (  51    6.8261e-02    1.3609e+00)  (  52    7.2838e-02    1.3579e+00)  
 (  53    7.7402e-02    1.3550e+00)  (  54    8.1955e-02    1.3520e+00)  
 (  55    8.6495e-02    1.3491e+00)  (  56    9.1025e-02    1.3462e+00)  
 (  57    9.5542e-02    1.3433e+00)  (  58    1.0005e-01    1.3405e+00)  
 (  59    1.0454e-01    1.3376e+00)  (  60    1.0903e-01    1.3347e+00)  
 (  61    1.1350e-01    1.3319e+00)  (  62    1.1796e-01    1.3291e+00)  
 (  63    1.2241e-01    1.3263e+00)  (  64    1.2685e-01    1.3235e+00)  
 (  65    1.3127e-01    1.3207e+00)  
 

   diag[r]=
   1.0000e+00  -1.3425e-01  -3.1981e-02   2.4580e-03   9.6653e-04   3.6442e-05
  -5.1135e-06  -1.3134e-06   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00   0.0000e+00
   0.0000e+00



 singular case : full regularized SQP
  del = 1.000000000000000e-01


  scalres=
   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.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.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.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.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.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.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.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.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.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.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.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.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.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.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.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.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.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.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.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.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
   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   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
   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   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
   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   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
   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   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
   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   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
   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   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
   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   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
   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   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
   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   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
   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   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
   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   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
   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   1.0000e+00   1.0000e+00   1.0000e+00
   1.0000e+00   1.0000e+