?? coga_c.m
字號:
% FS--計劃航班編號序列，大小1*F（函數(shù)中沒有出現(xiàn)，可以不用）,FSt--計劃航班類型序列，大小1*F(2表示大型機(jī),1表示中型機(jī),沒有小型機(jī));
% Fdashs--連續(xù)飛行的航班序列（為F的一個子集,表示為一個2進(jìn)制串（大小1*F'),第i個位置為1表示連續(xù)航班序列包括第i個航班); --本函數(shù)沒有用到
% FSu--計劃航班序列的各個起飛機(jī)場,;FSd--計劃航班序列的各個降落機(jī)場; --本函數(shù)沒有用到，二者大小都是1*F；
% zs--各個航班是否取消的向量,為一個二進(jìn)制串,大小1*F，第i個位置為1表示第i個航班取消,為0表示第i個航班不取消(這是一個決策向量);
% Cgf--地面延誤的單位時間成本(對大型機(jī)); Caf--空中延誤的單位時間成本(對大型機(jī))，這兩個都是單一數(shù)字
% Tpdfs/Tprfs--這是計劃航班序列的計劃起飛和降落的時間段(都是1*F向量);
% Mgfs/Mafs--計劃航班序列的依次最大地面/空中的等待時間數(shù)(都是1*F向量);--本函數(shù)沒有用到
% 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 -- Active Penalty Function
T = 40; % 仿真的總時段數(shù)
Cgf = 7; % 大型飛機(jī)地面單位時間成本
Caf = 10; % 大型飛機(jī)空中單位時間成本
%%%%%%%%%%%%%%%%%%%%%%% Generating Test Data
FS = []; % 沒有用處,但不要刪除
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); % 不連續(xù)航班的總數(shù)
% 不連續(xù)航班的起降機(jī)場表，F(xiàn)Su表示起飛機(jī)場序列，F(xiàn)Su中的1表示考慮容量限制的機(jī)場，其他不考慮者一律用2表示；
% FSd表示降落的機(jī)場序列，和上面類似；
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;]; % 飛機(jī)的類型序列，1表示中型飛機(jī)，2表示大型飛機(jī)
% Generate zs according to TdfsM & TrfsM
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;];
% 以上是不連續(xù)航班的標(biāo)準(zhǔn)起飛/降落時間表,用時段表示
Mgfs = 4.*ones(1,F); % 地面最多等待4個時間段
Mafs = Mgfs./2; % 空中最多等待2個時間段
StdTime = Tprfs - Tpdfs; % 標(biāo)準(zhǔn)飛行時間，實際飛行時間不能比這個短
Qfs = FSt.*20; % 取消費(fèi)用，大概相當(dāng)于單位時間損失的2倍
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 連續(xù)航班集合1（前面取消不會影響后面的航班序列）
FS1u = [ones(1,10),2*ones(1,4)]; % 連續(xù)航班1起飛機(jī)場序列
FS1d = 3 - FS1u; % 連續(xù)航班1降落機(jī)場序列
FS1t = [1,2,1,1,1,2,1,2,1,1,1,2,1,2]; % 連續(xù)航班1類型序列
Tpdf1s = [14,15,16,16,17,17,18,19,20,20,3,3,7,12]; % F1序列的標(biāo)準(zhǔn)起飛時間表
Tprf1s = [22,27,26,23,25,29,28,29,24,25,13,14,16,20]; % F1序列的標(biāo)準(zhǔn)降落時間表
F1con = [24,29,28,25,23,31,30,31,26,27,15,16,18,22]; % F1后續(xù)航班起飛的標(biāo)準(zhǔn)時間表
F1 = length(FS1u); % 連續(xù)航班1的總數(shù)
StdTime1 = Tprf1s - Tpdf1s; % 標(biāo)準(zhǔn)飛行時間
Mgf1s = 4.*ones(1,F1); % 地面最多等待4個時間段
Maf1s = Mgf1s./2; % 空中最多等待2個時間段
Qf1s = FS1t.*20; % 取消費(fèi)用
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 連續(xù)航班集合2（前面取消會影響后面的航班序列）
FS2u = [ones(1,11),2*ones(1,10)]; % 連續(xù)航班2起飛機(jī)場序列
FS2d = 3 - FS2u; % 連續(xù)航班2降落機(jī)場序列
FS2t = [2,1,1,1,2,2,1,1,2,1,1,2,2,1,1,2,2,1,2,2,1]; % 連續(xù)航班2類型序列
Tpdf2s = [13,13,15,15,16,17,17,17,19,19,20,5,3,3,5,6,9,11,5,7,16]; % F2序列的標(biāo)準(zhǔn)起飛時間表
Tprf2s = [21,21,23,27,30,25,24,29,33,32,25,13,14,14,15,16,16,17,18,18,20]; % F2序列的標(biāo)準(zhǔn)降落時間表
F2con = [23,23,25,29,32,27,27,31,35,34,29,15,16,16,17,18,18,19,20,21,22]; % F2后續(xù)航班起飛的標(biāo)準(zhǔn)時間表
F2 = length(FS2u); % 連續(xù)航班2的總數(shù)
StdTime2 = Tprf2s - Tpdf2s; % 標(biāo)準(zhǔn)飛行時間
Mgf2s = 4.*ones(1,F2); % 地面最多等待4個時間段
Maf2s = Mgf2s./2; % 空中最多等待2個時間段
Qf2s = FS2t.*20; % 取消費(fèi)用
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
tfsdash = 1; %松弛時間
%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 生成并修補(bǔ)非連續(xù)航班的染色體
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
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 生成并修補(bǔ)連續(xù)航班序列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(7.*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];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 生成并修補(bǔ)連續(xù)航班序列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(7.*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);
if temp > tfsdash + Mgf2s(jj)
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=300;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:NIND1
% reqmin(i)=T_P2Lag(Chrom1(i,:),bs2rv(Chrom2u,FieldD2),r);
reqmin(i) = Avi_Lag_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,Qf1s,FS2t,ChromF21(i,:),ChromF22(i,:),ChromzsF2(i,:),ChromzsF2_con(i,:),ChromF2_con(i,:),F2con,Tpdf2s,Tprf2s,F2,FS2u,FS2d,Qf2s);
end
[R,S]=size(reqmin);
if R==1
reqmin=reqmin';
end
%disp(reqmin);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%disp('part1 done');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for j=1:NIND2
reqmax(j)=Avi_Lag_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,Qf1s,FS2t,ChromF21(5,:),ChromF22(5,:),ChromzsF2(5,:),ChromzsF2_con(5,:),ChromF2_con(5,:),F2con,Tpdf2s,Tprf2s,F2,FS2u,FS2d,Qf2s);
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
tic;
for gen=1:MAXGEN
%tic;
for gen1=1:20
% GA begin
FitnVmin=ranking(reqmin);
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -