三维装箱和VRP结合的MATLAB算法

三维装箱和VRP结合的MATLAB算法,挺好的，需要输入数据后使用

%物体的长宽高
hi=[7 0 0; 15 0 0; 16 6 0; 17 0 0; 15 8 0; 12 6 16; 10 10 0; 12 7 9; 10 0 0; 14 0 0; 15 14 10; 9 14 7; 16 17 11; 14 12 8; 18 9 9];
wi=[5 0 0; 8 0 0; 15 5 0; 15 0 0; 7 10 0; 14 11 9; 7 7 0; 14 8 6; 7 0 0; 12 0 0; 15 13 13; 8 7 7; 11 13 10; 13 11 10; 12 12 6];
li=[30 0 0; 29 0 0; 33 36 0; 15 0 0; 13 15 0; 12 27 20; 23 21 0; 15 27 31; 24 0 0; 25 0 0; 31 19 16; 31 21 29; 34 26 28; 27 33 17; 22 23 34];
%客户的物品需求数
ni=[1 1 2 1 2 3 2 3 1 1 3 3 3 3 3];
ys=[1 0 0; 1 0 0; 0 0 0;0 0 0; 0 0 0; 0 0 1; 0 1 0; 0 0 1; 1 0 0; 0 0 0; 0 0 0; 0 0 0; 0 0 0; 0 0 0;0 1 1];
W=25;  
H=30;
L=60;%   车辆的长宽高
si=wi.*li; %每个货物的底面积
vi=si.*hi;  %每个货物的体积
ysa=0.7; %支撑面积约束
YS=0;
LIFO=0;
R=10000; %一个比较大的数
XP=[0 0 0 0 0 0;
    4 6 7 8 3 1;
    14 13 15 12 0 0;
    5 10 9 2 11 0;
    0 0 0 0 0 0];
XX=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
    6 0 5 1 0 2 3 4 0 0 0 0 0 0 0;
    0 0 0 0 0 0 0 0 0 0 0 4 2 1 3;
    0 4 0 0 1 0 0 0 3 2 5 0 0 0 0;
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
k=5;

km1=max（XX[]2）; %每一行的最大值，即车辆服务客户数
tic
itmax=50;
 pwz={};   %总的装箱位置
 pjg={};   %总的装箱结果
 zxcy=zeros（1itmax）;    %总的装箱惩罚值
 csj={};
 pjbwz1={};
 jbpar=XX;   %局部搜索的粒子就是XP 说明每个车辆服务的客户有哪些  罗列每辆车服务的客户是哪些
 
 for iter=1:itmax
     if iter==1
     jbcsj2=jbcsh（nihiviliwijbparkLIFOYSysXPkm1）; %求物品摆放的序列，jbcsj2是一个元胞
     [zxwz1zxjg1jbwz1cy1]=jbfitnesszx（LWHRhiwilinikjbcsj2XPysYSLIFO）;
     pwz{1iter}=zxwz1;
     pjg{1iter}=zxjg1;
    pjbwz1{1iter}=jbwz1;
    csj{1iter}=jbcsj2;
    bestwz=zxwz1;
    bestjg=zxjg1;
    bestfit=cy1;
    zxcy（1iter）=cy1;
     else
         
         jbcsj2=jbxh（csjpjbwz1iterk）
         [zxwz1zxjg1jbwz1cy1]=jbfitnesszx（LWHRhiwilinikjbcsj2XPysYSLIFO）
          pwz{1iter}=zxwz1;
          pjg{1iter}=zxjg1;
          pjbwz1{1iter}=jbwz1;
          csj{1iter}=jbcsj2;
          zxcy（1iter）=cy1;
          if cy1<=bestfit
              bestwz=zxwz1;
              bestjg=zxjg1;
              bestfit=cy1;
          end
     end
     if bestfit==0
         break
     end
 end
bestjg
bestwz
 bestfit
 toc

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       6707  2015-08-18 09:54  jbcsh.m

     文件       8638  2015-08-18 15:14  jbfitnesszx.asv

     文件       8597  2015-08-18 11:04  jbfitnesszx.m

     文件       2245  2015-08-18 09:56  cezhuangx.m

     文件       3478  2015-08-17 17:59  fitness1.m

     文件       6704  2015-08-18 09:51  jbcsh.asv

----------- ---------  ---------- -----  ----

                36369                    6