重力模型演示幻灯片.ppt
重力模型重力模型1重力模型法重力模型法（GravityMethod）模拟物理学中的模拟物理学中的牛顿的万有引力定律牛顿的万有引力定律基本假定：基本假定：交通区交通区i到交通区到交通区j的交通分布量的交通分布量与交通区与交通区i的交通量、交通区的交通量、交通区j的交通吸引量的交通吸引量成正比，与交通区成正比，与交通区i和和j之间的交通阻抗参数，之间的交通阻抗参数，如两区中心间交通的距离、时间或费用等成如两区中心间交通的距离、时间或费用等成反比。反比。2无约束重力模型无约束重力模型Casey在1955年提出了如下重力模型，该模型也是最早出现的重力模型：分别表示i小区和j小区的人口 表示i，j小区之间的距离 表示参数 模型本身不满足交模型本身不满足交通守恒约束条件：通守恒约束条件：(用出行费用函数 来表示)(用出行人数代替了总人数)3改进的重力模型可表示为：改进的重力模型可表示为：常见的交通阻抗函数有以下几种形式：常见的交通阻抗函数有以下几种形式：幂函数：幂函数：指数函数：指数函数：组合函数：组合函数：为参数，根据现状为参数，根据现状OD调查资料，利用最小二乘法调查资料，利用最小二乘法确定。确定。4例：例：按例按例3中表中表3和表和表4给出的现状给出的现状OD表和将来发生与吸引交通表和将来发生与吸引交通量，以及表量，以及表5和表和表6给出的现状和将来行驶时间，试利用重力模给出的现状和将来行驶时间，试利用重力模型和平均增长系数法，求出将来型和平均增长系数法，求出将来OD表。设定收敛标准为表。设定收敛标准为O/D123合计合计O/D123合计合计117.07.04.028.0138.627.038.06.051.0291.934.05.017.026.0336.0合计合计28.050.027.0105.0合计合计39.390.336.9166.5表3 现状OD表（单位：万次）表4 将来的发生与吸引交通量 512312317.017.022.014.09.011.0217.015.023.029.08.012.0322.023.07.0311.012.04.0表表5 5 现状行驶时间现状行驶时间 表表6 6 将来行驶时间将来行驶时间6解：（解：（1）用下面的无约束重力模型：）用下面的无约束重力模型：两边取对数，得两边取对数，得已知数据已知数据待标定参数待标定参数令：令：则：则：a0,a1,a2为待定系数7样本点样本点（）()()i=1,j=117282878472.83326.66441.9459i=1,j=2728501400171.94597.24422.8332i=1,j=342827756221.38636.62803.0910i=2,j=1751281428171.94597.26402.8332i=2,j=23851502550153.63767.84382.7081i=2,j=3651271377231.79187.22773.1355i=3,j=142628728221.38636.59033.0910i=3,j=2526501300231.60947.17013.1355i=3,j=317262770272.83326.55391.9459通过表通过表3和表和表5获取获取9个样本数据个样本数据8采用最小二乘法对这采用最小二乘法对这9个样本数据进行标定，得出个样本数据进行标定，得出=-2.084=1.173=-1.455 标定的重力模型为标定的重力模型为9O/D123合计合计188.86272.45818.940180.260275.542237.91246.164359.619318.79143.93276.048138.771合计合计183.195354.302141.152678.650(2)第一次计算得到的OD表 10（3）通过无约束重力模型计算得到的OD表不满足出行分布的约束条件，因此还要用其它方法继续进行迭代，这里采用平均增长系数法进行迭代计算。重新计算 和 计算结果如下面表所示 Ui为表4最后一列的值；Vj为表4最后一行的值Oi为每次计算得到的OD表每一行的合计值；Dj为每次计算得到的OD表 每一列的合计值11O/D123合计合计增长系数增长系数119.04616.9924.50440.5410.9521217.75560.71711.93390.4051.016534.45311.29719.80435.5541.0125合计合计41.25489.00536.241166.500增长系数增长系数0.95261.01451.0182用平均增长系数法第一次迭代计算用平均增长系数法第一次迭代计算OD表表12用平均增长系数法第三次迭代计算用平均增长系数法第三次迭代计算OD表表O/D123合计合计增长系数增长系数117.82316.6844.43838.9460.9911217.12762.31812.29191.7361.001834.27611.54420.31036.1300.9964合计合计39.22690.54637.040166.812增长系数增长系数1.00190.99730.996213（1 1）根据现状根据现状OD调查资料，利用最小二乘法确定参数，将确调查资料，利用最小二乘法确定参数，将确定的参数代入模型，得到标定的重力模型定的参数代入模型，得到标定的重力模型参数标定。（还有参数标定。（还有很多其他参数标定的方法）。很多其他参数标定的方法）。（2 2）利用）利用标定的重力模型计算得到表。标定的重力模型计算得到表。（3 3）无约束重力模型计算得到的无约束重力模型计算得到的OD表不满足出行分布的表不满足出行分布的约束条件，因此还要用其它方法继续进行迭代约束条件，因此还要用其它方法继续进行迭代。（例如：增长（例如：增长系数法等）系数法等）（4 4）迭代完成后得到最终的表。迭代完成后得到最终的表。重力计算步骤重力计算步骤14修正重力模型修正重力模型1.乌尔希斯重力模型乌尔希斯重力模型为交通阻抗函数,一般形式：待定系数待定系数根据现状根据现状OD调查资料拟和确定，一般可采用调查资料拟和确定，一般可采用试算法等数值方式，以某一指标作为控制目标，通过用模试算法等数值方式，以某一指标作为控制目标，通过用模型计算和实际调查所得指标的误差比较确定。型计算和实际调查所得指标的误差比较确定。的计算过的计算过程如下：程如下：15先假定一个 值，利用现状OD统计资料所得的 ，以及 代入模型中进行计算，所得出的计算交通分布称为GM分布。GM分布的平均行程时间采用下式计算：GM分布与现状分布的每次运行的平均行程时间之间的相对误差为 。当交通按GM分布与按实际分布每次运行的平均相对误差不大于某一限定值（常用3%）时，计算即可结束；当误差超过限定值时需改动待定系数 ，进行下一轮计算。调整方法为：如果GM分布的 大于现状分布 ，可增大 值；反之，则减小 值。162.美国公路局重力模型（B.P.R.模型）式中，为调整系数(也叫地域间结合度)，其计算公式为：其中，表示i小区到j小区的实际分布交通量与计算分布交通量之比；表示i小区到j小区的实际分布交通量与i小区的出行发生量之比。17的计算方法为：首先令 =1，根据现状OD表标定模型，计算 。将现状数据代入模型，计算出OD分布。根据上面的公式计算 。假定 的值在将来不发生变化，预测时不做任何修改而直接使用。标定 的方法与乌尔希斯重力模型 相同。18这两种模型均能满足出行产生约束条件，即：，因此都称为单约束重力模型。用上述两种重力模型进行交通分布预测时，首先是将预测的交通产生量和吸引量以及将来的交通阻抗参数带入模型进行计算。通常计算出的交通吸引量与给定的交通吸引量并不相同，因此需要进行进一步迭代计算。19为如下形式：为如下形式：双约束重力模型双约束重力模型2021以幂指数交通阻抗函数 为例介绍其计算方法：第1步：令m=0，m为迭代次数；第2步：给出 （可以用最小二乘法求出）；第3步：令 ，求出 ()；第4步：求出 和 ；第5步：收敛判定。若满足收敛条件，则结束计算；反之，令m+1=m，返回第2步重新计算。22优点优点（1 1）直观上容易理解；）直观上容易理解；（2 2）考虑路网的变化和土地利用对人们的出行产生的影响；）考虑路网的变化和土地利用对人们的出行产生的影响；（3 3）特定交通小区之间的）特定交通小区之间的ODOD交通量为零时，也能预测；交通量为零时，也能预测；（4 4）能比较敏感地反映交通小区之间行驶时间变化的情况。）能比较敏感地反映交通小区之间行驶时间变化的情况。缺点缺点（1 1）缺乏对人的出行行为的分析；）缺乏对人的出行行为的分析；（2 2）将出行费用视为定值；）将出行费用视为定值；（3 3）重力模型使用了同一时间段；）重力模型使用了同一时间段；（4 4）求内内交通量时的行驶时间难以给出；）求内内交通量时的行驶时间难以给出；（5 5）交通小区之间的距离小时，有夸大预测的可能性；）交通小区之间的距离小时，有夸大预测的可能性；（6 6）必须借助于其它方法进行收敛计算。）必须借助于其它方法进行收敛计算。重力模型的特点重力模型的特点23
