第五章-轨道交通系统客流预测(2).pptVIP

重力模型法（Gravity Method） 模拟物理学中的牛顿的万有引力定律。 基本假定：交通区i到交通区j的交通分布量与交通区i的交通发生量、交通区j的交通吸引量成正比，与交通区i和j之间的交通阻抗成反比。 无约束重力模型 其中: Tij——从分区i到分区j的出行分布量 Pi ——分区i的交通发生量 Aj ——分区j的交通吸引量 K ——平衡系数 tij——分区i，j之间的出行阻抗（距离、行程时间、广义费用） f(tij)——阻抗函数 α、?——模型系数 常见的交通阻抗函数有以下几种形式： 幂函数 指数函数 组合函数 最常用的形式 上式中，K,α,β,γ为模型系数，在现状OD表已知的条件下，Pi, Aj, Tij和tij已知，K,α,?,?可以用多元线性回归方法求得。 模型简化，一般可事先取下列经验值之一： （1） α=β （2） α=β=1.0 （3） α=β=0.5 若令α=β=1.0，则有： 无约束重力模型预测过程 ①利用现状OD数据，以及两小区之间的交通阻抗，用最小二乘法标定模型； ②利用建立的模型，根据交通生成和吸引预测得出未来个交通小区的产生量和吸引量，预测未来各交通小区的分布量； ③进行收敛性的计算，通常是利用平均增长系数法； 无约束重力模型计算出来的交通分布量，不满足上述约束的任何一个，通常要进行收敛迭代。 交通守恒约束： D O 1 2 3 合计 1 17 7 4 28 2 7 38 6 51 3 4 5 17 26 合计 28 50 27 105 表1 现状OD表 D O 1 2 3 合计 1 38.6 2 91.9 3 36 合计 39.3 90.3 36.9 166.5 表2 将来的交通发生吸引量 例4：表1和表2给出的现状OD表和将来发生与吸引交通量，以及表3和表4给出的现状和将来行驶时间，试利用重力模型和平均增长系数法，求出将来OD表。设定收敛标准为 1%。 D O 1 2 3 1 7 17 22 2 17 15 23 3 22 23 7 表3 现状行驶时间 D O 1 2 3 1 4 9 11 2 9 8 12 3 11 12 4 表4 将来行驶时间 解：（1）用下面的无约束重力模型： 两边取对数，得 已知数据 待标定参数 令： 则： 通过表1和表3获取9个样本数据 样本点 Tij Pi Aj PiAj tij y lnTij x1 ln(PiAj) x2 lntij i=1,j=1 17 28 28 784 7 2.8332 6.6644 1.9459 i=1,j=2 7 28 50 1400 17 1.9459 7.2442 2.8332 i=1,j=3 4 28 27 756 22 1.3863 6.6280 3.0910 i=2,j=1 7 51 28 1428 17 1.9459 7.2640 2.8332 i=2,j=2 38 51 50 2550 15 3.6376 7.8438 2.7081 i=2,j=3 6 51 27 1377 23 1.7918 7.2277 3.1355 i=3,j=1 4 26 28 728 22 1.3863 6.5903 3.0910 i=3,j=2 5 26 50 1300 23 1.6094 7.1701 3.1355 i=3,j=3 17 26 27 702 7 2.8332 6.5539 1.9459 采用最小二乘法对这9个样本数据进行标定，得出 =-2.084 =1.173 =-1.455 标定的重力模型为 (2) 第一次计算得到的OD表 D O 1 2 3 合计 1 88.862 72.458 18.940 180.260 2 75.542 237.912 46.164 359.619 3 18.791 43.932 76.048 138.771 合计 183.195 354.302 141.152 D O 1 2 3 合计 1 38.6 2 91.9 3 36 合计 39.3 90.3 36.9 166.5 表2 将来的交通发生吸引量 （3）通过无约束重力模型计算得到的OD表不满足出行分布的约束条件，因此还要用其它方法继续进行迭代，这里采用平均增长系