ratfitlo.pro

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

 (  47   -8.3927e-01    1.5283e+00)  (  48   -8.3181e-01    1.5230e+00)  
 (  49   -8.2438e-01    1.5178e+00)  (  50   -8.1697e-01    1.5126e+00)  
 (  51   -8.0958e-01    1.5074e+00)  (  52   -8.0222e-01    1.5023e+00)  
 (  53   -7.9488e-01    1.4973e+00)  (  54   -7.8757e-01    1.4923e+00)  
 (  55   -7.8027e-01    1.4873e+00)  (  56   -7.7300e-01    1.4824e+00)  
 (  57   -7.6576e-01    1.4775e+00)  (  58   -7.5853e-01    1.4726e+00)  
 (  59   -7.5133e-01    1.4678e+00)  (  60   -7.4415e-01    1.4631e+00)  
 (  61   -7.3699e-01    1.4583e+00)  (  62   -7.2985e-01    1.4536e+00)  
 (  63   -7.2274e-01    1.4490e+00)  (  64   -7.1564e-01    1.4444e+00)  
 (  65   -7.0857e-01    1.4398e+00)  (  66   -7.0151e-01    1.4352e+00)  
 (  67   -6.9448e-01    1.4307e+00)  (  68   -6.8747e-01    1.4262e+00)  
 (  69   -6.8048e-01    1.4218e+00)  (  70   -6.7351e-01    1.4174e+00)  
 (  71   -6.6656e-01    1.4130e+00)  (  72   -6.5963e-01    1.4087e+00)  
 (  73   -6.5272e-01    1.4043e+00)  (  74   -6.4582e-01    1.4001e+00)  
 (  75   -6.3895e-01    1.3958e+00)  (  76   -6.3210e-01    1.3916e+00)  
 (  77   -6.2527e-01    1.3874e+00)  (  78   -6.1845e-01    1.3832e+00)  
 (  79   -6.1166e-01    1.3791e+00)  (  80   -6.0488e-01    1.3750e+00)  
 (  81   -5.9812e-01    1.3709e+00)  (  82   -5.9138e-01    1.3668e+00)  
 (  83   -5.8466e-01    1.3628e+00)  (  84   -5.7796e-01    1.3588e+00)  
 (  85   -5.7128e-01    1.3548e+00)  (  86   -5.6461e-01    1.3509e+00)  
 (  87   -5.5796e-01    1.3470e+00)  (  88   -5.5133e-01    1.3431e+00)  
 (  89   -5.4472e-01    1.3392e+00)  (  90   -5.3813e-01    1.3353e+00)  
 (  91   -5.3155e-01    1.3315e+00)  (  92   -5.2499e-01    1.3277e+00)  
 (  93   -5.1845e-01    1.3239e+00)  (  94   -5.1192e-01    1.3202e+00)  
 (  95   -5.0541e-01    1.3164e+00)  (  96   -4.9892e-01    1.3127e+00)  
 (  97   -4.9245e-01    1.3090e+00)  (  98   -4.8599e-01    1.3053e+00)  
 (  99   -4.7955e-01    1.3017e+00)  ( 100   -4.7313e-01    1.2981e+00)  
 ( 101   -4.6672e-01    1.2944e+00)  ( 102   -4.6033e-01    1.2909e+00)  
 ( 103   -4.5396e-01    1.2873e+00)  ( 104   -4.4760e-01    1.2837e+00)  
 ( 105   -4.4126e-01    1.2802e+00)  ( 106   -4.3493e-01    1.2767e+00)  
 ( 107   -4.2862e-01    1.2732e+00)  ( 108   -4.2233e-01    1.2697e+00)  
 ( 109   -4.1606e-01    1.2662e+00)  ( 110   -4.0979e-01    1.2628e+00)  
 ( 111   -4.0355e-01    1.2594e+00)  ( 112   -3.9732e-01    1.2560e+00)  
 ( 113   -3.9111e-01    1.2526e+00)  ( 114   -3.8491e-01    1.2492e+00)  
 ( 115   -3.7873e-01    1.2458e+00)  ( 116   -3.7256e-01    1.2425e+00)  
 ( 117   -3.6641e-01    1.2392e+00)  ( 118   -3.6027e-01    1.2358e+00)  
 ( 119   -3.5415e-01    1.2325e+00)  ( 120   -3.4804e-01    1.2293e+00)  
 ( 121   -3.4195e-01    1.2260e+00)  ( 122   -3.3588e-01    1.2227e+00)  
 ( 123   -3.2982e-01    1.2195e+00)  ( 124   -3.2377e-01    1.2163e+00)  
 ( 125   -3.1774e-01    1.2131e+00)  ( 126   -3.1173e-01    1.2099e+00)  
 ( 127   -3.0573e-01    1.2067e+00)  ( 128   -2.9974e-01    1.2035e+00)  
 ( 129   -2.9377e-01    1.2004e+00)  ( 130   -2.8781e-01    1.1972e+00)  
 ( 131   -2.8187e-01    1.1941e+00)  ( 132   -2.7595e-01    1.1910e+00)  
 ( 133   -2.7003e-01    1.1879e+00)  ( 134   -2.6413e-01    1.1848e+00)  
 ( 135   -2.5825e-01    1.1817e+00)  ( 136   -2.5238e-01    1.1787e+00)  
 ( 137   -2.4653e-01    1.1756e+00)  ( 138   -2.4069e-01    1.1726e+00)  
 ( 139   -2.3486e-01    1.1695e+00)  ( 140   -2.2905e-01    1.1665e+00)  
 ( 141   -2.2325e-01    1.1635e+00)  ( 142   -2.1746e-01    1.1605e+00)  
 ( 143   -2.1170e-01    1.1575e+00)  ( 144   -2.0594e-01    1.1546e+00)  
 ( 145   -2.0020e-01    1.1516e+00)  ( 146   -1.9447e-01    1.1486e+00)  
 ( 147   -1.8876e-01    1.1457e+00)  ( 148   -1.8306e-01    1.1428e+00)  
 ( 149   -1.7737e-01    1.1399e+00)  ( 150   -1.7170e-01    1.1370e+00)  
 ( 151   -1.6604e-01    1.1341e+00)  ( 152   -1.6039e-01    1.1312e+00)  
 ( 153   -1.5476e-01    1.1283e+00)  ( 154   -1.4914e-01    1.1254e+00)  
 ( 155   -1.4354e-01    1.1226e+00)  ( 156   -1.3795e-01    1.1197e+00)  
 ( 157   -1.3237e-01    1.1169e+00)  ( 158   -1.2681e-01    1.1140e+00)  
 ( 159   -1.2126e-01    1.1112e+00)  ( 160   -1.1572e-01    1.1084e+00)  
 ( 161   -1.1020e-01    1.1056e+00)  ( 162   -1.0469e-01    1.1028e+00)  
 ( 163   -9.9197e-02    1.1000e+00)  ( 164   -9.3714e-02    1.0972e+00)  
 ( 165   -8.8245e-02    1.0945e+00)  ( 166   -8.2788e-02    1.0917e+00)  
 ( 167   -7.7346e-02    1.0890e+00)  ( 168   -7.1916e-02    1.0862e+00)  
 ( 169   -6.6500e-02    1.0835e+00)  ( 170   -6.1096e-02    1.0808e+00)  
 ( 171   -5.5706e-02    1.0781e+00)  ( 172   -5.0329e-02    1.0753e+00)  
 ( 173   -4.4965e-02    1.0726e+00)  ( 174   -3.9615e-02    1.0700e+00)  
 ( 175   -3.4277e-02    1.0673e+00)  ( 176   -2.8952e-02    1.0646e+00)  
 ( 177   -2.3640e-02    1.0619e+00)  ( 178   -1.8342e-02    1.0593e+00)  
 ( 179   -1.3056e-02    1.0566e+00)  ( 180   -7.7830e-03    1.0540e+00)  
 ( 181   -2.5230e-03    1.0513e+00)  ( 182    2.7242e-03    1.0487e+00)  
 ( 183    7.9584e-03    1.0461e+00)  
 exit from full SQP
            termination reason    1e+00
          final value of tauqp   1.000e+00
      sum norm of slack vector   0.000e+00


  phase=  0  scf0=  1.0000e+00
  d = 
  -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

  multipliers: first estimate
  u =
    1   3.5772e-01    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
   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   66   0.0000e+00   67   0.0000e+00   68   0.0000e+00
   69   0.0000e+00   70   0.0000e+00   71   0.0000e+00   72   0.0000e+00
   73   0.0000e+00   74   0.0000e+00   75   0.0000e+00   76   0.0000e+00
   77   0.0000e+00   78   0.0000e+00   79   0.0000e+00   80   0.0000e+00
   81   0.0000e+00   82   0.0000e+00   83   0.0000e+00   84   0.0000e+00
   85   0.0000e+00   86   0.0000e+00   87   0.0000e+00   88   0.0000e+00
   89   0.0000e+00   90   0.0000e+00   91   0.0000e+00   92   0.0000e+00
   93   0.0000e+00   94   0.0000e+00   95   0.0000e+00   96   0.0000e+00
   97   0.0000e+00   98   0.0000e+00   99   0.0000e+00  100   0.0000e+00
  101   0.0000e+00  102   0.0000e+00  103   0.0000e+00  104   0.0000e+00
  105   0.0000e+00  106   0.0000e+00  107   0.0000e+00  108   0.0000e+00
  109   0.0000e+00  110   0.0000e+00  111   0.0000e+00  112   0.0000e+00
  113   0.0000e+00  114   0.0000e+00  115   0.0000e+00  116   0.0000e+00
  117   0.0000e+00  118   0.0000e+00  119   0.0000e+00  120   0.0000e+00
  121   0.0000e+00  122   0.0000e+00  123   0.0000e+00  124   0.0000e+00
  125   0.0000e+00  126   0.0000e+00  127   0.0000e+00  128   0.0000e+00
  129   0.0000e+00  130   0.0000e+00  131   0.0000e+00  132   0.0000e+00
  133   0.0000e+00  134   0.0000e+00  135   0.0000e+00  136   0.0000e+00
  137   0.0000e+00  138   0.0000e+00  139   0.0000e+00  140   0.0000e+00
  141   0.0000e+00  142   0.0000e+00  143   0.0000e+00  144   0.0000e+00
  145   0.0000e+00  146   0.0000e+00  147   0.0000e+00  148   0.0000e+00
  149   0.0000e+00  150   0.0000e+00  151   0.0000e+00  152   0.0000e+00
  153   0.0000e+00  154   0.0000e+00  155   0.0000e+00  156   0.0000e+00
  157   0.0000e+00  158   0.0000e+00  159   0.0000e+00  160   0.0000e+00
  161   0.0000e+00  162   0.0000e+00  163   0.0000e+00  164   0.0000e+00
  165   0.0000e+00  166   0.0000e+00  167   0.0000e+00  168   0.0000e+00
  169   0.0000e+00  170   0.0000e+00  171   0.0000e+00  172   0.0000e+00
  173   0.0000e+00  174   0.0000e+00  175   0.0000e+00  176   0.0000e+00
  177   0.0000e+00  178   0.0000e+00  179   0.0000e+00  180   0.0000e+00
  181   0.0000e+00  182   0.0000e+00  183   0.0000e+00  184   0.0000e+00
  185   0.0000e+00  186   0.0000e+00  187   0.0000e+00  188   0.0000e+00
  189   0.0000e+00  190   0.0000e+00  191   0.0000e+00  192   0.0000e+00
  193   0.0000e+00  194   0.0000e+00  195   0.0000e+00  196   0.0000e+00
  197   0.0000e+00  198   0.0000e+00  199   0.0000e+00  200   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
  401   0.0000e+00  402   0.0000e+00  403   0.0000e+00  404   0.0000e+00
  405   0.0000e+00  406   0.0000e+00  407   0.0000e+00  408   0.0000e+00
  409   0.0000e+00  410   0.0000e+00  411   0.0000e+00  412   0.0000e+00
  413   0.0000e+00  414   0.0000e+00  415   0.0000e+00  416   0.0000e+00
  417   0.0000e+00  418   0.0000e+00  419   0.0000e+00  420   0.0000e+00
  421   0.0000e+00  422   0.0000e+00


  phase=  0  scf0=  1.0000e+00
  d =