【路徑規(guī)劃】基于遺傳算法求解帶時間窗帶容量公交車車輛路徑規(guī)劃問題（CVRPTW）matlab源碼
2021-10-21 23:05:43

%

%

clear

clc

close all

tic

%% 用importdata這個函數來讀取文件

c101=importdata('.datac101.txt');

cap=50;                                                       %車輛最大裝客量

%% 提取數據信息

E=c101(1,5);                                                   %發(fā)車中心時間窗開始時間

L=c101(1,6);                                                   %發(fā)車中心時間窗結束時間

vertexs=c101(:,2:3);                                           %所有點的坐標x和y

customer=vertexs(2:end,:);                                       %站點坐標

cusnum=size(customer,1);                                         %站點數

v_num=21;                                                       %車輛最多使用數目

demands=c101(2:end,4)./2;                                         %每個站點乘客量

a=c101(2:end,5);                                               %站點時間窗開始時間[a[i],b[i]]

b=c101(2:end,6);                                               %站點時間窗結束時間[a[i],b[i]]

s=c101(2:end,7);                                               %站點的服務時間

h=pdist(vertexs);

dist=squareform(h);                                             %距離矩陣，滿足三角關系，暫用距離表示花費c[i][j]=dist[i][j]

alpha=10;                                                       %違反的容量約束的懲罰函數系數

belta=100;                                                     %違反時間窗約束的懲罰函數系數

MAXGEN=100;                                                     %迭代次數

%% 遺傳算法

NIND=100;                                                       %種群大小

Pc=0.9;                                                         %交叉概率

Pm=0.05;                                                       %變異概率

GGAP=0.9;                                                       %代溝(Generation gap)

N=cusnum+v_num-1;                                               %染色體長度=站點數目+車輛最多使用數目-1

addpath('.GA')

[FF_GA,bestVC_GA]=ga(cusnum,a,b,L,s,dist,demands,cap,NIND,N,alpha,belta,GGAP,Pm,Pc,MAXGEN);

rmpath('.GA')

%%

%% 禁忌搜索設置參數

lamda=0.015;

delta=0.5;

vehicles_customer=cell(cusnum,1);                           %每輛車所經過的站點

TbLength=20; %禁忌長度

addpath('.TW')

[FF_TW,bestVC_TW]=TW(cusnum,a,b,L,s,dist,demands,cap,alpha,belta,MAXGEN,TbLength,delta,lamda);

rmpath('.TW')

disp('遺傳算法最優(yōu)解:')

draw_Best(bestVC_GA,vertexs);

title('遺傳算法最優(yōu)派車方案路線圖')

disp('禁忌搜索算法最優(yōu)解:')

draw_Best(bestVC_TW,vertexs);

title('禁忌搜索算法最優(yōu)派車方案路線圖')

figure(3);

hold on;box on

xlim([0,MAXGEN])

title('迭代曲線')

xlabel('代數')

ylabel('最優(yōu)值')

plot(1:MAXGEN,FF_GA,'b-',1:MAXGEN,[FF_TW,FF_TW(end)],'r-')

legend('遺傳算法','禁忌搜索算法')

toc