一种结合重力模型与Fratar模型的交通分布预测方法与流程

技术特征：

1.一种结合重力模型与Fratar模型的交通分布预测方法，其特征在于，包括以下步骤：

步骤10)采集各小区的现状发生吸引交通量、现状OD分布以及相关基础数据，其中，所述各小区的现状发生吸引交通量包括小区i的出行发生量O_i和小区j的出行吸引量D_j，i＝1、2...n，j＝1、2...n，n表示小区个数；

步骤20)标定无约束重力模型，包含步骤201)至步骤204)：

步骤201)确定无约束重力模型形式：q_ij表示小区i、j之间的交通量，c_ij表示小区i、j之间的阻抗，α，β，γ为无约束重力模型待标定参数；

步骤202)对两边取对数，得ln(q_ij)＝lnα+βln(O_iD_j)-γln(c_ij)；

步骤203)令Y＝ln(q_ij)，a₀＝lnα，a₁＝β，a₂＝-γ，X₁＝ln(O_iD_j)，X₂＝ln(c_ij)，则Y＝a₀+a₁X₁+a₂X₂，其中a₀，a₁，a₂为待标定系数，X₁，X₂，Y为包含样本数据集的向量；

步骤204)通过各小区现状发生吸引交通量及现状OD分布确定样本数据集X₁，X₂，Y，采用最小二乘法对样本数据进行标定，确定无约束重力模型参数；

步骤30)以现状发生吸引交通量为基础预测各小区未来年的发生吸引交通量，所述各小区未来年的发生吸引交通量包括第i个小区未来年的出行发生量P_i和第i个小区未来年的出行吸引量A_i；

步骤40)将步骤30)得到各小区未来年的发生吸引交通量代入步骤20)中既标定的无约束重力模型，计算得到未来年OD分布q_ij；

步骤50)：以步骤40)中应用无约束重力模型得到的OD分布q_ij作为初始OD，运行一次Fratar模型得到新的未来年预测OD分布q_ij¹，包含步骤501)至步骤503)；

步骤501)：依据无约束重力模型求解得到小区i、j之间OD分布q_ij，汇总得到各小区发生吸引交通量：

步骤502)：计算Fratar模型相关系数其中，F_Oi表示应用Fratar模型时第i个小区的出行发生收敛系数，F_Dj表示第j个小区的出行吸引收敛系数，L_ij表示第i个小区的出行发生调整系数，L_ji表示第j个小区的出行吸引调整系数；

步骤503)：求解新的未来年OD分布预测q_ij¹＝q_ij*F_Oi*F_Dj*(L_ii+L_jj)/2；

步骤60)：收敛性判断：根据步骤50)中Fratar模型一次循环得到的OD分布q_ij¹，汇总得到各小区新的未来年预测发生吸引交通量作为

$P_i^1=\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{q_{ij}}^1,A_i^1=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{q_{ij}}^1$

若各小区发生吸引交通量误差均在可接受范围以内，即：

$\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\[IF\left(0.99\<\frac{P_i}{P_i^1}\<1.01\right)\]+\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\[IF\left(0.99\<\frac{A_i}{A_i^1}\<1.01\right)\]==2n$

转入步骤70)，否则令转步骤40)；

步骤70)预测结果满足收敛条件，结束循环，此时的OD分布即为未来年各小区间OD分布预测结果。

2.根据权利要求1所述的结合重力模型与Fratar模型的交通分布预测方法，其特征在于：所述步骤10)中相关基础数据包括小区人口、面积、土地利用、区位因素数据。

3.根据权利要求1所述的结合重力模型与Fratar模型的交通分布预测方法，其特征在于：所述步骤30)中以现状发生吸引交通量为基础预测各小区未来年的发生吸引交通量的方法包括原单位法、增长率法、交叉分类法、函数法。