data1.asv

一个用拟牛顿算法求解优化问题的程序

mswp=[34  54  35;  
10  22  11   ;
15  24  16   ;
12  37  13   ;
27  67  28   ;
23   34  24   ;
4   7   5    ;
1   4   2     ;
29  49  30     ;
29  52  30     ;
52  97  53    ;
28  32  29     ;
27  4974  26   ;  
 35  74  36     ;
 53  80  54    ;
 55  91  56     ;
86  152   87     ;
246  636  247     ;
388  1141  389    ;
29   52   30     ;
173   294   174   ;  
383  561   384    ;
312  914   313     ;
6   14  7    ;
10  79  11    ;
12  27   13    ;
44  48   45   ;
48  51  49     ;
 6  14   7     ;
 18   39   19   ;
 7   11   8    ;
12   15   13    ;
12   15   13  ;
12    15  13   ;
13    16   14    ;
13    45  14    ;
63   340   64   ;
112   637  113    ;
216   1223   217  ;  
1   3   2     ;
1   3   2    ;
 1  3    2    ;
1   3   2    ;
2   10    3  ;   
 14   4988  12;    
2   11  3   ]
  

asll=[1   2   2     30  97   31     ;
2   2   2     9   36   10     ;
5   2   3     14  53   15    ;
6   2   2     12  24   13   ;
7   3   3     28  95  29    ;
8   3   15     17  69  18   ;
9   3   15     4   7   5    ;
11   3   10     1   4   2     ;
12   3   10     27   106   28   ;
13   4   4     29   52   30    ;
14   4   6     52   97   53   ;
15   4   11     28   32  29   ;
16   4   20     22   102  23    ;
18   6   13     34   161  35   ;
19   11   65     56   142  57    ;
20   20   31     55   91  56     ;
21   8   8     77   178  78    ;
21   50   50     241  629  242   ;
21   100   100     379  1188  380    ;
22   4   4     29  52  30    ;
23   2   2     183   309   184   ;
24   4   8     383   561  384   ;
24   50   100     312  914  313  ;  ;
25   2   4     6  14  7;
25   50   52     21  85  22    ;
26   3   3     12  27  13    ;
26   50   50     44  48  45    ;
26   100   100     48  51   49     ;
28   3   3     6  14  7     ;
28   10   10     18  39  19    ;
29   3   3     7  11   8    ;
29   50   50     12  15   13     ;
29   100   100     12   15   13   ;
29   200   200     12  15  13   ;
29   500   500     13   16  14   ;
30   3   3     13    45   14    ;
30   50   50     63  340  64    ;
30   100   100     112  637  113    ;
30   200   200     216  1223  217    ;
32   2   3     1  3   2     ;
32   50   100     1   3   2     ;
32   500   1000     1  3   2 ;
32   1000   1001     1   3  2    ;
33   2   3     2   10   3    ;
33   10   11     15   4987  13  ;
34   4   5     2   11  3 ]   
aa=[];
aa1=[];
for i=1:46
    aa(i)=mswp(i,2)+mswp(i,3)*5;
    aa1(i)=asll(i,5)+asll(i,6)*5;
end
for i=1:46
    gg=aa1(i)/aa(i);
end
for j=1:46