?? coga2_c.m
字號:
% FS--計劃航班編號序列，大小1*F（函數中沒有出現，可以不用）,FSt--計劃航班類型序列，大小1*F(2表示大型機,1表示中型機,沒有小型機);
% Fdashs--連續飛行的航班序列（為F的一個子集,表示為一個2進制串（大小1*F'),第i個位置為1表示連續航班序列包括第i個航班); --本函數沒有用到
% FSu--計劃航班序列的各個起飛機場,;FSd--計劃航班序列的各個降落機場; --本函數沒有用到，二者大小都是1*F；
% zs--各個航班是否取消的向量,為一個二進制串,大小1*F，第i個位置為1表示第i個航班取消,為0表示第i個航班不取消(這是一個決策向量);
% Cgf--地面延誤的單位時間成本(對大型機); Caf--空中延誤的單位時間成本(對大型機)，這兩個都是單一數字
% Tpdfs/Tprfs--這是計劃航班序列的計劃起飛和降落的時間段(都是1*F向量);
% Mgfs/Mafs--計劃航班序列的依次最大地面/空中的等待時間數(都是1*F向量);--本函數沒有用到
% TdfsM/TrfsM--計劃航班序列的實際起飛/降落時間(為兩個矩陣，大小都為F行T列，每一行代表一個航班的起（降）情況，且每行只有一個1，代表航班在該時刻位置起（降）);
% Qfs--計劃航班序列的取消成本，為1*F向量
% Only zs plays the role of j_v in this function
% F -- Num of total flights
% T -- Num of total time periods
% LagY & r -- Lagrangian Multiplexer
%%%%%%%%%%%%%%%%%%%%%%% Generating Test Data
FS = [];
T = 40;
% 只考慮單一機場的容量
FSu = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2;];
FSd = [2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1;];
F = length(FSu);
% 不連續航班的起降機場表，FSu表示起飛機場序列，FSu中的1表示考慮容量限制的機場，其他不考慮者一律用2表示；
% FSd表示降落的機場序列，和上面類似；
FSt = [2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,1,1,1,1,1,2,1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1;];
% 飛機的類型序列，1表示中型飛機，2表示大型飛機
% Generate zs according to TdfsM & TrfsM
Cgf = 7;
Caf = 10;
Tpdfs = [13,13,13,14,14,15,15,15,16,16,16,16,16,16,16,17,17,17,18,18,18,19,19,19,19,19,20,20,20,20,20,20,20,8,3,3,8,8,7,7,6,1,5,3,10,5,9,1,10,9,3,5,11,9,6,13,13,11,11,12,12,8,6,13,13,12,2,7,16,13,12;];
Tprfs = [20,20,23,26,22,21,23,27,24,31,31,27,27,24,25,21,28,29,25,29,23,25,33,31,25,33,24,27,36,25,32,32,27,13,13,13,13,13,13,13,14,15,15,15,15,15,15,15,16,16,16,16,16,16,16,16,17,18,18,19,19,19,20,20,20,20,20,20,20,20,20;];
% 以上是不連續航班的標準起降時間表
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
FS1u = [ones(1,10),2*ones(1,4)]; % 連續航班集合1（前面取消不會影響后面的航班序列）
FS1d = 3 - FS1u;
FS1t = [1,2,1,1,1,2,1,2,1,1,1,2,1,2]; % 航班類型序列
Tpdf1s = [14,15,16,16,17,17,18,19,20,20,3,3,7,12]; % F1序列的標準起降時間表
Tprf1s = [22,27,26,23,25,29,28,29,24,25,13,14,16,20];
F1con = [24,29,28,25,23,31,30,31,26,27,15,16,18,22];% F1后續航班起飛的標準時間表
F1 = length(FS1u);
StdTime1 = Tprf1s - Tpdf1s;
Mgf1s = 4.*ones(1,F1); % 地面最多等待4個時間段
Maf1s = Mgf1s./2; % 空中最多等待2個時間段
%Qf1s = FS1t.*20;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
FS2u = [ones(1,11),2*ones(1,10)]; % 連續航班集合2（前面取消會影響后面的航班序列）
FS2d = 3 - FS2u;
FS2t = [2,1,1,1,2,2,1,1,2,1,1,2,2,1,1,2,2,1,2,2,1];
Tpdf2s = [13,13,15,15,16,17,17,17,19,19,20,5,3,3,5,6,9,11,5,7,16]; % F2序列的標準起降時間表
Tprf2s = [21,21,23,27,30,25,24,29,33,32,25,13,14,14,15,16,16,17,18,18,20];
F2con = [23,23,25,29,32,27,27,31,35,34,29,15,16,16,17,18,18,19,20,21,22]; % F2后續航班起飛的標準時間表
F2 = length(FS2u);
StdTime2 = Tprf2s - Tpdf2s;
Mgf2s = 4.*ones(1,F2); % 地面最多等待4個時間段
Maf2s = Mgf2s./2; % 空中最多等待2個時間段
%Qf2s = FS2t.*20;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
tfsdash = 1; %松弛時間
%%%%%%%%%%%%%%%%%%%%%%%%%
Mgfs = 4.*ones(1,F); % 地面最多等待4個時間段
Mafs = Mgfs./2; % 空中最多等待2個時間段
StdTime = Tprfs - Tpdfs; % 標準飛行時間，實際飛行時間不能比這個短
Qfs = FSt.*20; % 取消費用，大概相當于單位時間損失的2倍
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 生成并修補非連續航班的染色體
NIND1=100;
Chrom11 = crtbp(NIND1,T.*ones(1,F));
for ii = 1:NIND1
for jj = 1:F
if Chrom11(ii,jj) <= Tpdfs(jj)
Chrom11(ii,jj) = Tpdfs(jj);
end
end
end % Repair chrom so that all the requirments are satisfied
Chrom12 = ceil(5.*rand(NIND1,F)) + Chrom11;
for ii = 1:NIND1
Chrom12(ii,:) = Chrom12(ii,:) + StdTime;
end
Chrom1 = [Chrom11,Chrom12];
Chromzs = zeros(NIND1,F); % 0 -- 不取消；1 -- 取消；
for ii = 1:NIND1
for jj = 1:F
if (Chrom11(ii,jj) > Tpdfs(jj) + Mgfs(jj))|(Chrom12(ii,jj) - Chrom11(ii,jj) > StdTime(jj) + Mafs(jj))
Chromzs(ii,jj) = 1; % 取消航班的情況（延時過長）
end
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 生成并修補連續航班序列F1(前面取消不影響后面)的染色體
ChromF11 = crtbp(NIND1,T.*ones(1,F1));
for ii = 1:NIND1
for jj = 1:F1
if ChromF11(ii,jj) <= Tpdf1s(jj)
ChromF11(ii,jj) = Tpdf1s(jj);
end
end
end
ChromF12 = ceil(5.*rand(NIND1,F1)) + ChromF11;
for ii = 1:NIND1
ChromF12(ii,:) = ChromF12(ii,:) + Tprf1s - Tpdf1s;
end
ChromzsF1 = zeros(NIND1,F1);
for ii = 1:NIND1
for jj = 1:F1
if (ChromF11(ii,jj) > Tpdf1s(jj) + Mgf1s(jj))|(ChromF12(ii,jj) - ChromF11(ii,jj) > Tprf1s(jj) - Tpdf1s(jj) + Maf1s(jj))
ChromzsF1(ii,jj) = 1; % 取消航班的情況（延時過長）
end
end
end
for ii = 1:NIND1
for jj = 1:F1
ChromzsF1_con(ii,jj) = round(rand(1));
ChromF1_con(ii,jj) = F1con(jj);
end
end
ChromF1 = [ChromF11,ChromF12];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 生成并修補連續航班序列F2(前面取消影響后面)的染色體
ChromF21 = crtbp(NIND1,T.*ones(1,F2));
for ii = 1:NIND1
for jj = 1:F2
if ChromF21(ii,jj) <= Tpdf2s(jj)
ChromF21(ii,jj) = Tpdf2s(jj);
end
end
end
ChromF22 = ceil(5.*rand(NIND1,F2)) + ChromF21;
for ii = 1:NIND1
ChromF22(ii,:) = ChromF21(ii,:) + Tprf2s - Tpdf2s;
end
ChromzsF2 = zeros(NIND1,F2);
for ii = 1:NIND1
for jj = 1:F2
if (ChromF21(ii,jj) > Tpdf2s(jj) + Mgf2s(jj))|(ChromF22(ii,jj) - ChromF21(ii,jj) > Tprf2s(jj) - Tpdf2s(jj) + Maf2s(jj))
ChromzsF2(ii,jj) = 1; % 取消航班的情況（延時過長）
end
end
end
for ii = 1:NIND1
for jj = 1:F2
if ChromzsF2(ii,jj) == 1
ChromzsF2_con(ii,jj) = 1;
ChromF2_con(ii,jj) = F2con(jj); % 這里的數值無關緊要
else
temp = ChromF22(ii,jj) - Tprf2s(jj); %%%%%%%%%%%%%%%%%%% WRONG,needs modification %%%%%%%%%%%%%%%%%%%%
if temp > tfsdash + 4
ChromzsF2_con(ii,jj) = 1; %%% 后續航班取消
ChromF2_con(ii,jj) = F2con(jj);
else
ChromzsF2_con(ii,jj) = 0; %%% 后續航班不取消
ChromF2_con(ii,jj) = temp + F2con(jj);
end
end
end
end
ChromF2 = [ChromF21,ChromF22];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
NIND2=20;
PRECI=10;
FieldD2=[rep([PRECI],[1,1]); rep([0;10],[1,1]);rep([1;0;1;1],[1,1])];
Chrom2=crtbp(NIND2,PRECI);
Chrom1u=Chrom1(5,:);
Chrom2u=Chrom2(1,:);
r=1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
MAXGEN=150;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:NIND1
% reqmin(i)=T_P2Lag(Chrom1(i,:),bs2rv(Chrom2u,FieldD2),r);
reqmin(i) = Avi_Lag2_C(FS,FSt,Chromzs(i,:),Chrom11(i,:),Chrom12(i,:),Cgf,Caf,Tpdfs,Tprfs,Qfs,Mgfs,Mafs,FSu,FSd,F,T,bs2rv(Chrom2u,FieldD2),r,FS1t,ChromF11(i,:),ChromF12(i,:),ChromzsF1(i,:),ChromzsF1_con(i,:),ChromF1_con(i,:),F1con,Tpdf1s,Tprf1s,F1,FS1u,FS1d,FS2t,ChromF21(i,:),ChromF22(i,:),ChromzsF2(i,:),ChromzsF2_con(i,:),ChromF2_con(i,:),F2con,Tpdf2s,Tprf2s,F2,FS2u,FS2d);
end
[R,S]=size(reqmin);
if R==1
reqmin=reqmin';
end
%disp(reqmin);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%disp('part1 done');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for j=1:NIND2
reqmax(j)=Avi_Lag2_C(FS,FSt,Chromzs(5,:),Chrom11(5,:),Chrom12(5,:),Cgf,Caf,Tpdfs,Tprfs,Qfs,Mgfs,Mafs,FSu,FSd,F,T,bs2rv(Chrom2(j,:),FieldD2),r,FS1t,ChromF11(5,:),ChromF12(5,:),ChromzsF1(5,:),ChromzsF1_con(5,:),ChromF1_con(5,:),F1con,Tpdf1s,Tprf1s,F1,FS1u,FS1d,FS2t,ChromF21(5,:),ChromF22(5,:),ChromzsF2(5,:),ChromzsF2_con(5,:),ChromF2_con(5,:),F2con,Tpdf2s,Tprf2s,F2,FS2u,FS2d);
end
[M,N]=size(reqmax);
if M==1
reqmax=reqmax';
end
%disp(Qfs);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%disp('part2 done');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Field1 = rep([1;T],[1,F]);
FieldF1 = rep([1;T],[1,F1]);
FieldF2 = rep([1;T],[1,F2]);
Field = [Field1,FieldF1,FieldF2];
Chrom = [Chrom1,ChromF1,ChromF2];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% CoGA Starts
for gen=1:MAXGEN
%tic;
for gen1=1:20
% GA begin
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -