function?[G_kG_best_routeG_best_lengthKbest_routebest_lengthlength_ave]=VRP（CNC_maxmAlphaBetaRhoQWt）?%#ok
%%
%--------------------------------------------------------------------------
%%參數說明
%%G_K:各代最優車輛數目
%%G_best_route：各代最優路線
%%G_best_length：各代最優路線長度
%%K：最優車輛數目
%%best_route：最優路線
%%best_length：最優路線長度
%%length_ave：各代平均長度
%%C：DC和工廠的坐標
%%NC_max：最大迭代次數
%%m：螞蟻數目
%%Alpha：重要度系數
%%Beta：能見度系數
%%Rho：揮發度系數
%%Q：信息更新參數
%%W:車輛載重量
%%
%%第一步?初始化變量和系數

m=60;Alpha=1;Beta=1;gama=2;Rho=0.15;NC_max=50;Q=15;W=9;qq=0.05;?%#ok
C=[0?0
????0?-1
????0?3
????-2?-2
????-3?-3
????3?-1
????-4?0
????-4?-1
????1?-2
????1?-1
????1?3
????3?4
????-3?0
????2?0
????1?-3
????2?-1
????2?1
????1?-4
????-3?2
????-1?-1];?%#ok
t=[0?1.5?1.8?2?0.8?1.5?1.0?2.5?3.0?1.7?0.6?0.2?2.4?1.9?2.0?0.7?0.5?2.2?3.1?0.1];

%構造倉庫和工廠等之間的距離矩陣

n=size（C1）;%n表示問題的規模（城市個數）?
D=zeros（nn）;%D表示完全圖的賦權鄰接矩陣?
for?i=1:n?
for?j=1:n?
if?i~=j?
D（ij）=（（C（i1）-C（j1））^2+（C（i2）-C（j2））^2）^0.5;?
else?
D（ij）=eps;?
end?
D（ji）=D（ij）;?
end?
end?

%構造節省量矩陣
U=zeros（nn）;%U表示工廠之間的連接和其與倉庫之間連接能節省的距離
for?i=1:n?
for?j=1:n?
if?i~=j?
U（ij）=D（i1）+D（j1）-D（ij）;
else?
U（ij）=eps;?
end?
U（ji）=U（ij）;?
end?
end?
load_w=0;
Eta=1./D;%Eta為啟發因子，這里設為距離的倒數?
Tau=ones（nn）;%Tau為信息素矩陣?
Tabu=zeros（mn+20）;%存儲并記錄路徑的生成?
NC=1;%迭代計數器?
G_best_route=[NC_maxn+20];%各代最佳路線?
G_best_length=inf.*ones（NC_max1）;%各代最佳路線的長度?
length_ave=zeros（NC_max1）;%各代路線的平均長度

%%第二步，把螞蟻放到DC內
while?NC<=NC_max%停止條件之一：達到最大迭代次數?
????Tabu（:1）=randint（m1[11]）;
%%第三步，m只螞蟻按照要求的方法選擇工廠，并完成周游?
for?i=1:m
????visited=Tabu（i:）;
????visited=visited（visited>0）;
????to_visit=setdiff（1:nvisited）;
????c_temp=length（to_visit）;
????j=1;
???????while?j<=n
??????????if?~isempty（to_visit）
????%visit_thisant=[];
????%while?to_visit~=[1]
????%for?ee=1:（2*c_temp）?
???%while?c_temp>=2
???%if?length（to_visit）>1
????%visited=Tabu（i:）;
????%to_visit=setdiff（1:nvisited）;
????%to_visit=[1to_visit];
????%c_temp=length（to_visit）;
????%visit_thisant=[];
???%?if?c_temp~=1
????????%for?b=1:c_temp
????????????%if?（load_w+t（to_visit（b）））<=W
????????????????%visit_thisant=[visit_thisantto_visit（b）];
????????????%end
????????%end
????????%jj=visited（visited>0）;
%%?按照規則選下一個工廠或者是回到倉庫
%for?k=1:length（visit_thisant）
????for?k=1:length（to_visit）
??x（k）=（Tau（visited（end）to_visit（k））^Alpha）*（Eta（visited（end）to_visit（k））^Beta）*（U（visited（end）to_visit（k））^gama）;
????end???
????????ww=rand;
if?ww????Select=find（max（x））;
????%Tabu（ilength（visited（visited>0））+1）=to_visit（Select（1））;?
else
x=x/（sum（x））;?
%按概率原則選取下一個城市?
xcum=cumsum（x）;?
Select=find（xcum>=rand）;
%Tabu（ilength（visited（visited>0））+1）=to_visi（Select（1））;?
end
if?isempty（Select）
????Select=1;
????load_w=load_w+t（Select）;
else
load_w=load_w+t（to_visit（Select（1）））;
end
if?load_w>W
????Select=1;
???????j=j-1;
????load_w=0;
????Tabu（ilength（visited）+1）=Select（1）;
else
Tab

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????215??2019-07-14?21:46??b?（2）\VVRPe\hslogic.txt

?????文件???????5592??2008-12-19?14:33??b?（2）\VVRPe\VRP.m

?????目錄??????????0??2019-07-24?19:24??b?（2）\VVRPe

?????目錄??????????0??2018-08-20?18:00??b?（2）

-----------?---------??----------?-----??----

?????????????????5807????????????????????4