capacitycon.asv

遗传算法例程：通过MATLAB

function y = CapacityCon(gk,bk,Cthres,Cpdk,Cprk,t)

% gk为第k个机场容量曲线的系数，bk为第k个机场容量曲线的范围，二者大小都是1*4，见word文档
% Cthres为容量曲线拐点处的数值，大小1*3，(对应横轴变量Cpdk的几个奇异点数值）

% Cpdk和Cprk是第k个机场的预测数值,大小1*T（决策变量）

if Cpdk(t) <= Cthres(1)
    sum = bk(1) - Cprk(t) - gk(1).*Cpdk(t);
elseif Cpdk(t) > Cthres(1) & Cpdk(t) <= Cthres(2)
    sum = bk(2) - Cprk(t) - gk(2).*Cpdk(t);
elseif Cpdk(t) > Cthres(2) & Cpdk(t) <= Cthres(3)
    sum = bk(3) - Cprk(t) - gk(3).*Cpdk(t);
elseif Cpdk(t) > Cthres(3)
    sum = bk(4) - Cprk(t) - gk(4).*Cpdk(t);
end

y = sum;