VRP问题 matlab源码 车辆路径

function P=make_it_better(center,customer,ff) [~,jieyue]=distance(center,customer); A=jieyue; %节约矩阵 B=[customer,ff]; %客户点及需求量 i=0; %记录有几条线路 examination=1; %控制每一条线路,使它的需求量小于maxiumA while(sum(sum(A>-2))>=1) %剩余三个点及以上的时候继续循环找下一条线路. %% 每一条线路的探索初始化 i=i+1; [m,w]=max(A); %w记录最大元所在的行数 [~,W]=max(m); %W记录列数 h=W;l=w(W);p{i}=[l,h]; %% 爬取该线路的所有点,限制条件为不超过maxium while(examination<=160) A(:,p{i})=-2; P{1,i}=p{i}; %记录当前线路 if max(max(A(p{i},:)))==-2 break end % 在不超过运输量的情况下选择节约最大的点并入线路 [w,~]=max(A(p{i},:)); %行 T=find(w>-2); %除去-2 [~,l]=sort(w(T),'descend'); %排序 q=find(B(T(l),3)<=(160-examination)); if isempty(q) break else t=T(l(q(1))); end p{i}=[P{1,i},t]; %更新线路 examination=sum(B(p{i},3)); end A(P{1,i},:)=-2; examination=1; %% 只剩一个点时候的处理 [tt0,~]=max(A);tt3=find(tt0>-2); if length(tt3)==1 i=i+1;P{1,i}=tt3; break end end