机场延误 遗传算法

各个飞机由于其飞行情况，乘坐旅客类型，重要程度，旅客对航班延误的耐受度的不同，从而对航班延误的敏感性不同。即有一些航班非常重要，一旦延误一点时间就会造成巨大损失，有一些航班则不那么重要，就算延误一点时间也无伤大雅。由于现有的飞行空间有限，可以通过调整各个航班飞机的起飞次序可以使得航班延误造成的损失（时间和经济）最少。很显然这是可以通过数学上的最优化模型来进行优化，而这种优化很大程度上是对航班起飞计划的优化而很少对硬件设施进行优化，所以花费的成本较低，但产生的效果较好（运输延误和起飞延误的权重相差不大）。所以我们主要考虑对流量控制进行优化。

%function y = Aviation_Lag（FSFStzsCgfCafTpdfsTprfsQfsMgfsMafsFTLagYr）

%LagY = 10;
%r = 1.1;

FS = [];

F = 100;
T = 40;

FSu = ones（1F）;
FSd = 2*FSu;

for i = 1:F
    z = rand（1）;
    if z > 0.7
        FSt（i） = 2;
    else
        FSt（i） = 1;
    end
end

Fdashs = [];

for i = 1:F
    z = rand（1）;
    if z > 0.8
        zs（i） = 1;
    else
        zs（i） = 0;
    end
end

Cgf = 7;
Caf = 10;

for i = 1:F
Tpdfs（i） = ceil（rand（1）.*10）;
Tprfs（i） = ceil（T - 10.*rand（1））;
end

Mgfs = T.*ones（1F）;
Mafs = Mgfs;

Qfs = FSt.*100;

% 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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Create TdfsM/TrfsM Matrix Using standard GA
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
alpha = 1/3;   % Set a threshold for GA process
Cthres = sum（Qfs）.*alpha;   
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

StdTime = Tprfs - Tpdfs;
% Inner GA starts

NIND = 20;            % num of individuals
GGAP=1;               % generation gap
XOVR=0.7;               % rate of recombine
MUTR=1/4;          % rate of mutation
MAXGEN=10;           % Max number of generation
INSR=0.8;             % probabily of insertion
SUBPOP=1;             % number of subpop
MIGR=0.2;             % rate of migration
MIGGEN=20; % every 20 gen miggration
kk=1;

baseVT = （T-StdTime-1）.*ones（1F）;   % 计时时间段数从1开始，最大到T - StdTime因航班而异，而这是起飞的时间段向量范围；
ChromT = crtbp（NINDbaseVT）; % 随机构造标准的起飞/降落时间表，严格满足要求：先起飞后降落，且都在预定的整个时间段内完成飞行任务；
ChromT1 = ChromT + 1;
baseVL = StdTime + Mgfs + Mafs;       % 这是降落的时间段向量范围，要叠加在起飞时间段的基础上
BX = max（ChromT1）;
TempB = T.*ones（1F） - BX;
for i = 1:F
    if baseVL（i） > TempB（i）
        baseVL（i） = TempB（i）;
    end
end
ChromL = crtbp（NINDbaseVL）;
ChromL = ChromL + ChromT1;
ChromStd = [ChromT1ChromL];
%zox = ChromStd（:1+F:2*F） - ChromStd（:1:F）

gen = 0;

TdfsM = zeros（FT）;
TrfsM = zeros（FT）;

for i = 1:NIND
    
% TdfsM/TrfsM--计划航班序列的实际起飞/降落时间（为两个矩阵，大小都为F行T列，每一行代表一个航班的起（降）情况，且每行只有一个1，代表航班在该时刻位置起（降））;
    for mm = 1:F
        for nn = 1:T
            if ChromT（imm） == nn
                TdfsM（mmnn） = 1;
            end
            if ChromL（imm） == nn
                TrfsM（mmnn） = 1;
     

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       6693  2007-02-24 17:30  GA\GA\Aviation_Lag.asv

     文件       6781  2007-02-25 17:35  GA\GA\Aviation_Lag.m

     文件       3977  2007-03-04 00:34  GA\GA\Avi_Lag.m

     文件       4124  2007-02-28 22:21  GA\GA\Avi_Lag2.asv

     文件       4231  2007-03-04 00:34  GA\GA\Avi_Lag2.m

     文件       4517  2007-03-10 16:36  GA\GA\Avi_Lag2_C.asv

     文件       4988  2007-03-10 17:13  GA\GA\Avi_Lag2_C.m

     文件       3983  2015-04-25 22:13  GA\GA\Avi_Laggai.m

     文件       4720  2007-03-23 20:55  GA\GA\Avi_Lag_C.m

     文件       2296  2007-02-27 10:13  GA\GA\Avi_Viol.m

     文件       3217  2007-02-26 01:47  GA\GA\bs2rv.m

     文件        630  2007-02-15 00:42  GA\GA\CapacityCon.asv

     文件        750  2007-02-15 00:52  GA\GA\CapacityCon.m

     文件        750  2007-03-06 21:49  GA\GA\CapacityCon_C.m

     文件       7989  2015-04-25 21:09  GA\GA\CoGA.m

     文件       3674  2006-08-30 17:23  GA\GA\CoGA2forT_P2.m

     文件      18471  2007-03-11 14:59  GA\GA\CoGA2_C.m

     文件       9133  2015-05-01 15:20  GA\GA\CoGAgai.m

     文件       7766  2007-03-02 20:21  GA\GA\CoGA_2.m

     文件      18748  2007-03-23 21:04  GA\GA\CoGA_C.m

     文件      18748  2015-05-01 15:22  GA\GA\CoGA_Cgai.m

     文件       7348  2015-04-25 21:33  GA\GA\CoGA_V.m

     文件        128  2007-03-04 11:12  GA\GA\ContinousConF1.asv

     文件        400  2015-04-25 22:18  GA\GA\create.m

     文件       1168  2007-02-24 16:53  GA\GA\crtbase.m

     文件       2187  2007-02-24 16:53  GA\GA\crtbp.m

     文件       2093  2007-02-24 16:53  GA\GA\CRTRP.m

     文件      13511  2007-03-03 21:56  GA\GA\HangBan1.mat

     文件      15912  2007-02-25 19:29  GA\GA\hs_err_pid3112.log

     文件        491  2007-03-02 17:30  GA\GA\LandCon.m

............此处省略39个文件信息