clear
%A點集合，第一個點表示配送中心。
%?A=[50?19?49?49?53?85?33?59?33?82?70;
%?????50?16?92?30?42?90?74?32?23?59?41;
%????0?1.15?1.34?0.53?1.79?1.47?0.01?0.95?0.26?1.21?0.69];?
%?A=[];
load?A.mat
rong=3500;?
%?rong=2100;%最大負載
m=size（A2）;%客戶數量
c=zeros（mm）;%距離表
for?j=1:m?????????????????????????????????????????????????????
????for?i=（j+1）:m
????????c（ij）=sqrt（（A（1i）-A（1j））^2+（A（2i）-A（2j））^2）;
????end
end
?
p=zeros（m-1m-1）;%節約里程表去掉第一個配送中心
for?j=2:（m-1）?????????????????????????????????????????????????
????for?i=（j+1）:m
????????p（i-1j-1）=c（i1）+c（j1）-c（ij）;
????end
end
?
s=p（:）;?
????????????????????????????????
[hswz]=sort（s1‘descend‘）;%降序排序，hs值，wz對應的索引
hs（find（hs==0））=[];?%去掉值為0的元素
for?i=1:size（hs）
????[sub（i1）sub（i2）]=ind2sub（size（p）wz（i））;???%將P矩陣各元素索引轉換成坐標，并存儲到sub矩陣
end;
svt=[hssub];%%得到節約里程的降序表，值，行號，列號。
route=zeros（m-1m-1）;
sv=0;
Zhuang=[];
for?j=1:m-1
?for?i=1:size（svt1）?%求從最大節約值開始，可一起配送的兩個客戶，作為初始解
??if?A（3（svt（i2）+1））+A（3（svt（i3）+1））<=rong?????????????
?????solut=[svt（i2）svt（i3）];
?????sv=sv+svt（i1）;%存累計節約值
?????zhuang=A（3（solut（11）+1））+A（3（solut（12）+1））;%存當前載重量
?????ii=i;
?????break
??end
?end

?for?i=（ii+1）:size（svt）?
????
??????if?????（svt（i2）==solut（11））&&（isempty（find（svt（i3）==solut））==1）&&（（A（3（svt（i3）+1））+zhuang）<=rong）;???%從最大的小于初始解對應的最大節約值對應的坐標判斷（左坐標等于最優解的左坐標，并且右坐標不等于最優解的右坐標，并且容量不超）
??????????????solut=[svt（i3）solut];??%如滿足條件，將右坐標加到路徑的左側
??????????????sv=sv+svt（i1）;
??????????????zhuang=A（3（svt（i3）+1））+zhuang;
??????elseif?（svt（i2）==solut（1length（solut）））&&（isempty（find（svt（i3）==solut））==1）&&（（A（3（svt（i3）+1））+zhuang）<=rong）;
??????????????solut=[solutsvt（i3）];%如滿足條件，將右坐標加到路徑的右側
??????????????sv=sv+svt（i1）;
??????????????zhuang=A（3（svt（i3）+1））+zhuang;
??????elseif?（svt（i3）==solut（11））&&（isempty（find（svt（i2）==solut））==1）&&（（A（3（svt（i2）+1））+zhuang）<=rong）;
??????????????solut=[svt（i2）solut];
??????????????sv=sv+svt（i1）;
??????????????zhuang=A（3（svt（i2）+1））+zhuang;
??????elseif??（svt（i3）==solut（1length（solut）））&&（isempty（find（svt（i2）==solut））==1）&&（（A（3（svt（i2）+1））+zhuang）<=rong）;
??????????????solut=[solutsvt（i2）];
??????????????sv=sv+svt（i1）;
??????????????zhuang=A（3（svt（i2）+1））+zhuang;
??????else
??????????continue
??????end
?end

?for?i=size（svt）:-1:1???%刪除已經選到路徑中的點
?????if?（（isempty（find（svt（i2）==solut））==0）||（isempty（find（svt（i3）==solut））==0））;
?????????svt（i:）=[];
?????else
?????????continue
?????end
?end

route（j（1:length（solut）））=solut;??%將確定好的某一條路徑存到route矩陣的一行中
Zhuang=[Zhuangzhuang];

if?isempty（svt）??%直到判斷svt無元素，退出
????break
end
end

for?i=1:m-1??%判斷route中沒有的單一用戶，增加一行，加入到route中
????if?（isempty（find（route（:）==i））==1）
????????route（j+11）=i;
????????j=j+1;
????end
end
opt=2*sum（c（:1））-sv;%行駛距離

plot（A（11）A（21）‘s‘）
%?text（A（11）A（21）[‘（‘num2str（A（11））‘‘num2str（A（21））‘）‘]）
hold?on;

for?i=2:m
????plot?（A（1i）A（2i）‘o‘）
%?????text（A（1i）A（2i）[‘（‘num2s

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1572??2020-02-16?15:53??節約里程法\A.mat
?????文件????????4276??2020-04-26?16:03??節約里程法\cw.m
?????文件???????10918??2020-02-16?19:48??節約里程法\data（1）.xlsx
?????文件????????1085??2020-02-16?16:35??節約里程法\dis2dot.m
?????文件?????2127855??2020-02-16?14:53??節約里程法\考慮碳排放的多車場多車型VRP模型及算法研究_劉麗姣.caj
?????目錄???????????0??2020-02-16?19:48??節約里程法\