四阶段法-出行分布计算 底特律法

底特律法计算未来出行分布

% trans_planning 计算未来出行量分布
function main()
clc;clear;
confirm = 0;
while confirm  == 0
    inner = 0; % 是否使用算例标识
    e = input('确定本次计算使用误差精度（推荐0.03）：');
    OD = input('请输入当前交通量方阵（输入1使用算例1，输入2使用算例2）：');
    if OD == 1
        inner = 1;confirm = 1;
        Gi1 = [2 1 4 3];% 产生量增长系数
        Gj1 = Gi1;       % 吸引量增长系数
        OD1 = [ 0 25 10 5;
        0 0 15 10;
        0 0 0 5;
        0 0 0 0];  
        OD1 = OD1 + OD1'; % 初始出行分布
        OD = OD1;
        Gi = Gi1;Gj = Gj1;
    elseif OD ==2
        inner = 1;confirm = 1;
        OD2 = [17 7 4;
        7 38 6;
        4 5 17];
        P_0 = sum(OD2,1);
        A_0 = sum(OD2,2)';
        P_f = [38.6 91.9 36];
        A_f = [39.3 90.3 36.9];
        Gi2 = P_f./P_0;
        Gj2 = A_f./A_0;