MATLAB下料问题原创程序

tic
clear all; clc
l = [1743	1680	1532	1477	1313	1285	1232	1217	1180	1177  1105	1055	1046	1032	1030	975	893	882	847	845 ...
    830	795	766	745	730	719	714	695	645	630 610	600	590	588	582	578	540	488	355	343 320	515	414	411	405	328	313	290	275	265 ...
    255	184	155 ];
n = [ 4	216	104	38	4	60	4	8	6	10  8	2	4	8	8	2	8	301	6	38 30	8	4	4	34	18	4	4	90	30 ...
     30	212	108	482	196	8	32	4	52	42  8	8	8	60	136	4	68	286	602	286   192	57	24];
% [lpai weizhi ] = sort（l‘descend‘）;
% npai = n（weizhi）;
% qing = [l ; n];
TotalLeft = 1:53; NLeft0 = n;  FourDay = [  5791215182025 283648]; SixDay = [ 41124293238404650 ]; YiYoLiao = 0; 
ChangeCostLv = .08*.01; FourDayLeft =  FourDay;  SixDayLeft =  SixDay;  NFourDayLeft = n（FourDay）; NSixDayLeft = n（SixDay）; 
WasteLiao = 0;   NumMax = 1; 

for jj = 1:NumMax
    TotalLeft = 1:53; NLeft = NLeft0（TotalLeft）; ii = 0;
    ShengLiao = 3000; 
    Others = TotalLeft;
    ShengLiaoTotal = zeros（1length（Others））;
    for i = 1:length（Others）
          lDangQian = l（ Others（i） ）;  NumberOfPresent = 0; 
          ShengLiaoo = ShengLiao; 
          while ShengLiaoo >= lDangQian 
               if   ShengLiaoo < lDangQian + 5
                        ShengLiaoo = 0;
               else
                       ShengLiaoo = ShengLiaoo - lDangQian - 5;
                end
          end
               ShengLiaoTotal（i） = ShengLiaoo;
                  
      end
           
      [ShengLiaoTotal weizhi] = min（ShengLiaoTotal）;
      ShengLiaoTotal = [];
      DangQian = Others（weizhi）;  lDangQian = l（DangQian）;  
      
     while  ~isempty（TotalLeft） 
          MoShi = [];
          ii = ii + 1; tt = 0; ShengLiao = 3000;
%           DangQian = TotalLeft（1）;  lDangQian = l（DangQian）;
          while ShengLiao >= lDangQian  
               if   ShengLiao < lDangQian + 5  
                    tt = tt + 1;
%                     NLeft（weizhi） = NLeft（weizhi） -1;
                     NLeft（weizhi） = NLeft（weizhi） -1;
                    ShengLiao = 0;
                    MoShi（tt） = DangQian;
               else 
                    tt = tt + 1;
%                     NLeft（weizhi） = NLeft（weizhi） -1;
                     NLeft（weizhi） = NLeft（weizhi） -1;
                    ShengLiao = ShengLiao - lDangQian - 5;
                    MoShi（tt） = DangQian;
                          
               end
               TotalLeft（NLeft==0） = [];
               NLeft（NLeft==0） = [];
               if  isempty（TotalLeft）
                   break;
               end
               
                Others = TotalLeft;
                for i = 1:length（Others）
                   lDangQian = l（ Others（i） ）;  
                   ShengLiaoo = ShengLiao; 
                   while ShengLiaoo >= lDangQian 
                        if   ShengLiaoo < lDangQian + 5
                             ShengLiaoo = 0;
                        else
                            ShengLia