遗传算法多车辆路径优化，粒子群算法多车辆路径优化（代码完整，数据齐全）

%D 距离矩阵 %NIND 为种群个数 %X 参数是中国34个城市的坐标(初始给定) %MAXGEN 为停止代数，遗传到第MAXGEN代时程序停止,MAXGEN的具体取值视问题的规模和耗费的时间而定 %m 为适值淘汰加速指数,最好取为1,2,3,4,不宜太大 %Pc 交叉概率 %Pm 变异概率 %输出： %R 为最短路径 %Rlength 为路径长度 clear clc close all %% 加载数据 % load CityPosition2 [a1,ax1,ay1] = xlsread('各点坐标.xlsx'); xynum0 = []; for ii = 1:20 xy = ay1{ii+2,4}; for jj = 1:length(xy) if strcmp(xy(jj),',') xynum0 = [xynum0 ;str2double(xy(2:jj-1)) str2double(xy(jj+1:end-1))]; end end end numd = zeros(20); for ii = 1:20 for jj = 1:20 numd(ii,jj) = jwdis(xynum0(ii,1),xynum0(ii,2),xynum0(jj,1),xynum0(jj,2)); end end xynum = []; for ii = 1:20 xy = ay1{ii+2,4}; for jj = 1:length(xy) if strcmp(xy(jj),',') xynum = [xynum ;str2double(xy(2:jj-1)) str2double(xy(jj+1:end-1))]; end end end X=xynum; a2 = xlsread('各点需求量.xlsx'); a2num= a2(:,3); a3 = xlsread('点与点之间的运输距离.xlsx'); a3num = a3(2:end,2:end); a3num = numd; D=Distanse(xynum); %生成距离矩阵 D =a3num; N=size(D,1); %城市个数 %% 遗传参数 NIND=1000; %种群大小 MAXGEN=50; %最大遗传代数 Pc=0.9; %交叉概率 Pm=0.05; %变异概率 GGAP=0.9; %代沟 %% 初始化种群 Chrom=InitPop(NIND,N); %% 画出随机解的路径图 % DrawPath(Chrom(1,:),X) pause(0.0001) % %% 输出随机解的路径和总距离 % disp('初始种群中的一个随机值:') % OutputPath(Chrom(1,:)); Rlength=PathLength(D,Chrom(1,:),a2num); disp(['总距离：',num2str(Rlength)]); disp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~') %% 优化 gen=0; figure; hold on;box on xlim([0,MAXGEN]) title('优化过程') xlabel('代数') ylabel('最优值') ObjV=PathLength(D,Chrom,a2num); %计算路径长度 preObjV=min(ObjV); while gen<MAXGEN %% 计算适应度 ObjV=PathLength(D,Chrom,a2num); %计算路径长度 % fprintf('%d %1.10f\n',gen,min(ObjV)) line([gen-1,gen],[preObjV,min(ObjV)]);pause(0.0001) preObjV=min(ObjV); FitnV=Fitness(ObjV); %% 选择 SelCh=Select(Chrom,FitnV,GGAP); %% 交叉操作 SelCh=Recombin(SelCh,Pc); %% 变异 SelCh=Mutate(SelCh,Pm); %% 逆转操作 SelCh=Reverse(SelCh,D,a2num); %% 重插入子代的新种群 Chrom=Reins(Chrom,SelCh,ObjV); %% 更新迭代次数 gen=gen+1 ; end %% 画出最优解的路径图 load maydata.mat ObjV=PathLength(D,Chrom,a2num); %计算路径长度 [minObjV,minInd]=min(ObjV); % DrawPath(Chrom(minInd(1),:),X) %% 输出最优解的路径和总距离 disp('最优解:') p=OutputPath(Chrom(minInd(1),:),a2num,xynum0); disp(['总距离：',num2str(ObjV(minInd(1)))]); disp('-------------------------------------------------------------')