clc;clear?all;
format?long;
M=11;
N=12;
zq=100;
%c=[ceil（rand*M）?ceil（rand*M）?ceil（rand*M）?ceil（rand*M）?ceil（rand*M）?ceil（rand*M）?ceil（rand*M）?ceil（rand*M）?ceil（rand*M）?ceil（rand*M）?ceil（rand*M）?ceil（rand*M）];%定義染色體編碼，各個貨物對應(yīng)的車1*12
c=ceil（rand（zqN）*M）;
L=[3.86?6.9?4.87?7.15?7.15?7.15?3.9?7.15?7.15?7.15?7.994?7.12];%各個12貨物的長度1*12
L=L+0.2;
K=13;%箱子的長度
e=zeros（1N+1）;
e（N+1）=100000;

for?jhcs=1:10
s=zeros（zqM）;
%初始化種群
for?zqs=1:zq
for?xh=1:40
????for?j=1:M
???????s（zqsj）=0;?
????end
A=zeros（MN）;
for?j=1:M%遍歷各個車
????t=1;
????for?k=1:N%遍歷各個貨物
????????if（c（zqsk）==j）
????????????s（zqsj）=s（zqsj）+L（k）;%計算每輛箱子上貨物的長度s[]
????????????A（jt）=k;
????????????t=t+1;
????????end
????end
end
for?j=1:M??%遍歷各個車
????if（（s（zqsj）==0）&&（j????????for?i=1:N
????????????if?（c（zqsi）==max（c（zqs:）））
????????????????c（zqsi）=j;
????????????end
????????end
????????%break;
????end
end

for?j=1:M??%遍歷各個車
????d=10;
????T=0;
????if（s（zqsj）>=K）%第j輛車貨物超過箱子長度
????????for?t=1:N%建立此箱子對應(yīng)的貨物長度矩陣
????????????if?（（A（jt）~=0）&&（d>L（A（jt））））
????????????????d=L（A（jt））;???????????%尋找最短長度
????????????????p=j;q=t;???????????????????%保存序號
????????????end
????????end
????????%找到最短的車序號之后再放到每個車上看看長度夠不夠，夠就插入
????????for?j=1:M
????????????if?（s（zqsj）+d???????????????%將此貨物放入此車中?
???????????????c（zqsA（pq））=j;
???????????????T=1;
???????????????break;
????????????end
????????end
????????if?（T==0）%表明所有的車上貨物都不能容納新貨物
????????????c（zqsA（pq））=max（c（zqs:））+1;
????????????break;
????????end
????end
end
end
end

num=zeros（1zq）;
f=zeros（1zq）;
for?zqs=1:zq
????%計算每個種群中所用箱子的個數(shù)，存于num中
????for?j=1:M
????if?（s（zqsj）>0）
????????num（zqs）=num（zqs）+1;
????end
????end
????f（zqs）=num（zqs）*K*（K*num（zqs）-sum（s（zqs:）））;
end

d=[cf‘];
d=sortrows（dN+1）;
if（d（N+1）????e=d;
end
d（:N+1）=[];
c=d;

%[sfindex]=sort（f）;
%for?zqs=1:zq
?%???if?（zqs??%??b=c（index（zqs）:）;
???%?c（index（zqs）: