应用重力模型进行交通分布的详细步骤
第三章 交通需求预测-重力模型
l 基本假设为：交通区i到交通区j的出行分布
1、无约束重力模型
l
万有引力模型
模型为：
F = K⋅
l
量与i区的出行发生量、j区的出行吸引量成正 比，与i区和j区之间的交通阻抗成反比。 根据对约束情况的不同分类，重力模型有三种 形式：无约束重力模型、单约束重力模型和双 约束重力模型。
X ij = k ⋅
可采用先考虑宏观因素预测其总比例，再考 虑微观因素预测各交通区间出行方式的分担率的 方法。 出行总比例预测 条件类方式：根据车辆拥有量进行预测。
l
竞争类方式 取决于需求的出行方式其总比例预测可根据有关的社 会经济发展目标，结合其发展实际状况，通过综合分 析求得。如出租车等。 取决于有关政策的出行方式其总比例预测可按照有关 的发展策略，根据已有基础进行规划确定，如公交车 出行等。 各交通区间的出行比例预测 各交通区间某种出行方式的出行比例取决于该出行方 式的总比例、出行目的结构和出行距离，通过前述模 型以及根据出行调查等资料统计分析拟合建立的现状 关系曲线进行预测。
其中, ｃ:汽车(car)；ｂ:公共汽车(bus)
l
∑∑ A
j m
其中Tijm——从交通区i到交通j，第 m种交通方 式的交通量；
4、 回归模型法——产生分担组合模型
l
二、交通方式的分类
l l l l l l l
该模型是通过建立交通方式分担率与其相关因素 间的回归方程，作为预测交通方式模型。
可分为：自由类、条件类和竞争类。 1、自由类交通方式 主要指步行交通，影响因素（内在因素）包括： 出行目的、出行距离、气候条件等 2、条件类交通方式 主要指单位小汽车、单位大客车、私人小汽车、摩托 车等交通方式 影响因素（外在因素）包括：有关政策、社会、经济 的发展水平。 影响因素（内在因素）包括：车辆拥有量、出行目的、 出行距离等。
重力模型标定方法及分析
以幂 函数双约束 重力模 型的标 定为例 ：
＾
＝
由于重力模型是从万有引力 的定律 抽象而来 ， 因此 用于交通 分析时必然存在 一定 的弊端 。首 先模 型的物 理意义是 揭示 人 的 活动的社会现象 ， 又没有 完全立 足于人 的 出行规 律。其次 ， 但 对
收 稿 日期 ：０ ２０ —０ ２ １ －２１ 作者 简 介 ： 季 凯 （９ ７ ， ， 士 ， 程 师 １７ ．） 男 硕 工
其中 ， Ａ 为小 区 的交通吸引总量 ； Ｉ Ｐ 为小 区 ｉ 的交通发生总
第 ０ １２年 １ 期 ８卷 １ ２３ 第 ４ 月
季
凯： 重力模型标定方法及分析
Ｋ Ｐ ／Ｃ 。 ； ； Ａ
Ｊ
Ｋ ：［ ｉ ∑ ＡＣ］ 。 ／ ￣～
＝
于阻抗 函数而言 ， 出行距 离的 系数 为 常数 的假 设不 符合 实际 ， 对 对于距离太小 的情 况 ， 有时 预测值 过高 ， 导致 与实 际值 的误 差 较 大 。此外 ， 小区内出行时 间较难 确定 ， 使得 小 区内交通 量预 测结
， ｃ ）＝ ｅｐ ｆ 。 （ ｘ （ Ｃ ） ｌ ４ 多项式 函数 ： ）
广义费进行 测度 。 常用的阻抗 函数形式包括 以下 ５类 ：
１ 幂 函数 ： ）
， Ｃ ）＝ｃ 。 （
１ １ 重力模 型分 类 ．
随着交通研究者 的不 断努 力开 发 ， 重力模 型在 表达 形式 、 参 数标定 与检验 方法上 已有多种 形式 。按 照表达形 式 的不 同可 分
第 ３ 卷 第 １ 期 ８ １
・
１ ・ ８
２０ １２ 年 ４ 月
山 西 建 筑
交通规划设计之重力模型法
)
C -n ij
exp（
Cij）
二、重力模型参数标定
在 现 状 OD 表 已 知 的 条 件 下 ， Oi, Dj, Cij 和 tij 已 知 ， k,α,,可以用最小二乘法求得。对（7.3.1）式取对数:
tij
k
Oi
D
j
Cij
（7.3.1）
ln tij ln k ln Oi ln Dj ln Cij
1,
1 bmj 1 /ቤተ መጻሕፍቲ ባይዱbmj 1
第三步，将求得参数代入，用现状OD值求现 状理论分布表｛ tˆij ｝
第四步，计算现状实际OD分布表的平均交通
阻抗 R 1 ti
tij Rij
j
，再计算理论分布表中的平均
交通阻抗： Rˆ 1 ti
tij Rij
j
，求两者相对误差。
第五步，如果 满足要求(<3%)，则接受γ 的求解，否则： ①若 Rˆ R ，则理论分布量小于实际， 应减 少γ的值，可令γ＝ γ／2； ②若 ，则理论分布量大于实际，应增 大γ的值，令γ＝2 γ； ③返回第一步，重新计算。
j
S.t.
t ij Oi
j
Kij为调整系数，采用试算法确定。当Kij =1时，
即为乌尔希斯重力模型。
Kij (1 Yij )ij /(1 Yijij )
• λij—i小区到j小区的实际分布交通量与计算分布交 通量之比，λij=tij/t’ij； • Yij—i小区到j小区的实际分布交通量与i小区的出 行产生量之比，Yij=tij/Oi。
R 1503 100 2 505 4003 1005 200 4 3.4 1000
Rˆ 147.63 95.7 2 56.75 402.43 104.35 193.3 4 1000
transcad四阶段法
六、Connect连接 Connect连接就是将小区中心点连接到路网上， 目的是将小区的属性与路网的属性进行融合， 使小区真正成为路网中的一部分，只有这样才 能对路网进行交通分配。 在Connect连接之前需在路网节点层(Endpoints) 添加一个字段(table)，取名为index，在路网层 也添加一个字段，也取名为index。这些都是为 下一步ID转换作准备的。
四、用重力模型进行交通分布预测
注意点：
Dataview中须 选小区中心点 层 Productions选 出行发生量(O) Attractions选 出行吸引量(D)
五、检验路网的连通性 该步骤用以检验路网中各路段与节点之间是否 真正连接。其操作过程为：路网层Tools → Map Editing → Check Line Layer Connectivity 在Threshold中填上检查的间隔。
二、生成小区图
新建交通区层（面层） New－Geographic File,选择Area Geographic File， options选择第二项，输入层名:area 添加域字段：交通区编码ID，现状base，现状吸引量A base，现状交通区容量POP base，未来发生量P fur， 未来吸引量A fur，未来交通区容量POP fur 画交通区：Tools－Map Editing－Toolbox
应用重力模型进行交通分布的详细步骤
应用重力模型进行交通分布的详细步骤第一步：求阻抗矩阵Rij（Impedance Matrix）交通阻抗可表示为：出行距离和行程时间的长短，以及出行费用的大小等。
为真实地反映交通阻抗，依托工程道路网规划采用通常使用的平均行程时间表示。
小区之间的阻抗——平均行程时间越小表示小区之间阻抗越小，越大表示小区之间阻抗越大，因此以平均行程时间为道路权值求各小区之间的最短路径（Shortest Path），其值即为小区之间的阻抗R ij。
1、数据准备（1）创建路网图1表示的是TransCAD创建路网的界面。
（2）做选择集。
在Endpoints层，于dataview中选择质心点，将其作为一个选择集。
（3）各路段平均行程时间（Travel time）其中，平均行程时间=Length/平均车速*602、操作过程Networks/Paths—Multiple paths调出其对话框如图2所示。
3、运行结果（即为阻抗矩阵），如图3所示。
第二步：重力模型标定（校准）（Gravity Mode Calibration）1、数据准备（1）基础OD矩阵。
（2）阻抗矩阵（Shortest Paths），如图3所示。
重力模型标定（校准）（Gravity Mode Calibration）数据准备：基年OD矩阵的索引（质心层质心ID）与最短路径矩阵的索引（路网节点层质心ID）不匹配，并且因为下面将在路网节点层上操作，因此必须使基年OD索引与最短路径矩阵的索引相一致，以使两表数据相对应（转换为“质心ID”）。
操作方法：按其对话框4示意操作。
2、操作过程按对话框（如图5）操作即可。
3、运行结果（1）标定参数结果（这里选用伽马函数）：a=2.6288，b=0.2361，c=0.0，如图6所示，不过大看show report 里面参数更准确。
（2）K－Factor Flow：如图7所示。
第三步：创建综合阻抗因子f (Rij) (Synthetic Friction Factors)1、数据准备（1）创建空矩阵“Friction Factor shell”；（2）已标定的a、b、c值；（3）阻抗（最短路径）矩阵，如图8所示。
重力模型的解释及系数计算方法
重力模型法(gravity model)是一种最常用的方法，它根据牛顿的万有引力定律，即两物体间的引力与两物体的质量之积成正比，而与它们之间距离的平方成反比类推而成。
下式为Casey(1955)提出的重力模型。
其中，：i，j小区的人口； d为i，j小区间的距离,α为系数。
上式的约束条件为：s.t.同时满足守恒条件的α是不存在的，因此，将重力模型修改如下：其中，为交通阻抗函数。
交通阻抗函数的几种形式：指数函数：（1）幂函数：（2）组合函数：（3）为参数。
单约束型B.P.R.模型其中，调整系数。
发生侧得到保证，即：以下以幂指数交通阻抗函数为例介绍其计算方法：第1步令m=0，m为计算次数。
第2步给出n（可以用最小二乘法求出）。
第3步令第4步求出第5步收敛判定。
若下式满足，则结束计算；反之，令m+1=m，返回第２步重复计算。
,作业：按上次作业给出的现状OD表和将来生成、发生与吸引交通量，利用下式重力模型和弗拉塔算法，求出将来OD表。
收敛标准。
重力模型:其中，，，。
读者也可以利用以前给出的现状分布交通量和表4-1示现状行驶时间，估计出这3个参数。
表4-1 现状行驶时间表4-2将来行驶时间解：利用重力模型求解分布交通量如下：同理，可以计算出其它各交通小区之间的交通量如下表所示。
重力模型的优点：a.直观上容易理解；b.能考虑路网的变化和土地利用对人们的出行产生的影响；c.特定交通小区之间的OD交通量为零时，也能预测；d.能比较敏感地反映交通小区之间行驶时间变化的情况。
重力模型的缺点：a.重力模型仅仅是将物理法则简单直观上容易理解；b.能考虑路网的变化和土地利用对地应用到社会现象，尽管有类似性，需要更加贴合人们出行的方法；c.一般，人们的出行距离分布在全区域并非为定值，而重力模型将其视为定值；d.交通小区之间的行驶时间因交通方式和时间段的不同而异，而重力模型使用了同一时间；e.求内内交通量时的行驶时间难以给出；f.交通小区之间的距离小时，有夸大预测的可能性；g.利用重力模型计算出的分布交通量必须借助于其它方法进行收敛计算。
trsc交通规划实例详细步骤
t r s c交通规划实例详细步骤内部编号：（YUUT-TBBY-MMUT-URRUY-UOOY-DBUYI-0128）T r a n s C A D 交通规划实例作者：applepin启动TransCAD软件插入软件狗，启动TransCAD软件。
建立路段层点击新建文件图标，出现点击建立地理文件，点击OK，出现选择文更改名更改名选中Line 层地理文件类型，并更改名字，点击OK，出现路段层属性对话框：点击Add Field，逐一增加字段属性点击OK，出现保存对话框：选择保存路径及文件名称点击Save,出现路段层地图界面添加路段在路段层上，点击Tools→map editing→toolbox出现地图编辑工具栏：点击，逐一增加路段应用保存得到7条线路、5个节点的路网。
输入路段属性数据打开路段层数据库：逐一输入路段属性数据(通行时间等于路段长度除以速度)建立小区层点击新建文件图标，出现点击建立地理文件，点击OK，出现点击OK，出现小区层属性对话框逐一添加小区的字段属性点击OK，出现保存对话框：起好名字，点击Save。
再次回到路段层地图界面。
画小区。
点击Tools→map editing→toolbox出现工具栏：画出3个小区，并保存。
输入小区属性数据打开小区层数据库输入各小区的属性数据将小区质心点连接到路网（目的是做ID转换）。
在节点层上，增加Index。
原先节点层上只有经纬度两个字段两个字段点击数据库菜单，修改数据库属性增加一个新字段Index在小区层上，点击Tools→Map Editing →Connect 调出对话框：点击Fill,完成如下设置：点出OK，路网显示出已经连接，出现小区质心节点打开点层数据库，发现新的变化：打开路段层数据库，发现新的变化：增加的三个节点的小区质心节增加了6,7，8三个填充连接后新增路段的填充连接后新增路段(质心连杆)的值。
将其通行能力设为无穷大（大数即可）的值，通行时间设为很小的值。
TransCAD运用重力模型预测详细步骤
TransCAD 运用重力模型进行预测详细步骤一、将CAD 创建的DXF 文件导入到TransCAD用TransCAD 软件打开.dxf 文件会弹出窗口，按照下图所示操作即可，Road 层和Area 层需要分两次导入。
导入路网层：点击OK ，命名为Road 并保存（Save ）。
导入小区层：点击点击OK ，命名为Area 并保存（Save ）。
Road 层和Area 层叠加：右键选择Layers 重命名图层：点击点击检查路段连接性：tools-map editing-check line layerconnectivity无问题，则显示为黑色点儿，继续下一步；有问题则显示为其他颜色点儿，需要进行调整。
二、分别给Road 层和Area 层建立相应字段并填充数据切换到Road 层，点击Dataview ，选择Modify Table ，弹出如下窗口，建立luming 、danxiangchedaoshu 、daoludengji 、danchedaonengli 、AB-T 、BA-T 、AB-V 、BA-V 、AB-C 、BA-C 和reallength字段。
切换到Endpoints 层，点击Dataview ，选择Modify Table ，弹出如下窗口，建立xiaoqu字段。
点击填充Road层数据表：其中luming 、danxiangchedaoshu 、daoludengji 根据调查结果进行填充，本例快速路AB-V 和BA-V 取60km/h ，主干路AB-V 和BA-V 取40km/h ，次干路AB-V 和BA-V 取30km/h ，支路AB-V 和BA-V 取20km/h ，不同等级道路单车道通行能力视隔离情况（机非隔离、中央隔离）进行取值，本例中快速路取值1100pcu/h ；主干路无隔离，则取950pcu/h ，有一种隔离措施，则取1000pcu/h ，有两种隔离措施，则取1050pcu/h ；次干路则依次取值为800pcu/h 、850pcu/h 、900pcu/h ；支路无机非隔离带则取450pcu/h ，有机非隔离带则取650pcu/h 。
论述交通量预测的增长率法和重力模型法
论述交通量预测的增长率法和重力模型法0、引言所谓的交通量分布就是区与区之间的交通流，现状的区与区之间的交通分布已经从OD表中体现出来了。
交通量分布预算的目的就是根据现状OD分布量及各区因经济增长、土地开发等形成的交通量增长来推算各区之间将来的交通分布。
交通量预测主要有增长率法和重力模型法两种方法。
1、增长率法预测分布交通量增长率法是从已知的现有OD调查表和发生、吸引交通量的增长率求出OD 分布交通量的近似值，其次对、就、进行收敛计算，从而求得将来的分布交通量。
增长率法包括平均增长率法、底特律法和弗雷特法。
1.1平均增长率法1.1.1平均增长率法计算步骤①根据计算公式计算将来出行量式中：—区到区的将来出行量；—区到区的现在出行量；—区出行发生的增长系数；—区出行吸引的增长系数②检验吸引量和发生量是否与推算的交通量相符合，是否满足（为判定值），如符合计算完成；如不符合需要在第一轮的基础上重新计算增长系数，并重复步骤①，直到满足上述要求为止。
1.1.2平均增长率法算例【例1】已知1、2、3区的出行、增长系数及现状分布，如表2.1-1所示，求将来的出行分布。
（取）表1.1-1 出行、增长系数及现状分布解：求间的交通量于是有，，，，，，，，得到第一轮计算结果，如表1.1-2所示。
表1.1-2 第一轮计算结果因第一轮计算结果中新的调整系数不能满足的要求，因此需要进行第二轮计算，直到满足要求为止。
本例共需要进行四轮计算，才可得到最终结果。
1.2底特律法此方法假定区到区间的交通量同和成比例增加。
1.2.1底特律法计算步骤①根据计算公式计算将来出行量式中：其中—未来发生量合计；—未来吸引量合计②检验吸引量和发生量是否与推算的交通量相符合，是否满足（为判定值），如符合计算完成；如不符合需要在第一轮的基础上重新计算增长系数，并重复步骤①，直到满足上述要求为止。
1.3弗雷特法该方法假设，小区之间OD交通量的增长系数不仅与小区的发生增长系数和小区的吸引增长系数有关，还与整个规划区域的其他交通小区的增长系数有关。
重力模型在交通量分布中的应用研究
２ 重 力模 型 的标 定 。
２ ．１ 理 论 依 据 重力模型的标定实际上是对分离系数
ｏｏ ， ４
Ｋ ： ２×ｏ ｏ ×ｌ ｏ ０ ．４ ＿ ．８
，路运输 象 鼹壤撬越 厶 ＼ 爻遍量分本 南 的程序来 ｊ 完成 ，重 力模 型实际上是增长系数法的进
３ １ 第 一 次 迭代 ， ．
分布 中的应用研究
舔 靠 ，轰 Ｉ 交通量 粳型 恒其 盘 茨 务毒 ＆多 于公 路运输 系统礴络 规 划不多 ，本ｉ 玲鳇重 誊
魏粳 壁 ＆ 吼 骨 交 通 量 分 夺 的 摊 峨 ｜ ｜
ＺＡ ，
。
／、 （） ＾ ３
： ：
２
维普资讯
瞧 惫 嗨 晕熹 毒
ｃ 氓
ｉ
憾
、
钧鸭
｜ 譬
由于冬季气温变化较大，降雪对道岔
的对锁闭影响较大 ，我们在冬季加大对道 岔密贴、开程、油压等各项指标进行调整
和 测 试 ，冬 季 对雪 为 令 ， 及 时 清 扫 积 雪 ，
１五 个 区域 如 图 １ ） ；
‘＝ １
（） ２
（ ） 中：Ｔ＿ 区域 ｉ Ｊ １式 １ 到 的交通量 。
Ｐ 一区域 ｉ 的交 通 发 生 量 。 Ａ 区域 ｊ 的交 通 吸 引 量 。
２两个区域是交通量发生区（ ） ①和②
３三 个 区域是 交通 量 吸 引 区（ ，④和 ） ③
⑤ ） ；
广 工 ［ 匝—正 ］ —］
Ｋ 一区域 ｉ ｊ 到 的社会经济调整系数。 Ｆ 一 区域 ｉ ｊ 分 离 系数 。 到 的
臣 １２巨０１ ｌ １ ０ ５ ２ ２ ０ ３｛ ． ６ ．８ ８ １ ８ ８
（） ２ 式中： 『 区域 ｉ 的阻抗系数 ， ｗ＿－ 到ｊ
重力模型的过程和原理教案
重力模型的过程和原理教案重力模型是一种经济地理学中常用的分析工具，用于研究地理空间内不同地区之间的贸易、人口流动、投资等现象。
它通过考虑地理距离和经济规模大小两个因素，揭示了地理空间的相互作用对区域之间的联系和互动产生的影响。
一、重力模型的基本原理重力模型的基本原理是基于物理学中的引力定律，即两个物体之间的引力与它们的质量成正比，与它们之间的距离成反比。
将这个物理定律应用到经济地理学中，我们可以认为两个地区之间的联系程度与它们的经济规模大小成正比，与它们之间的地理距离成反比。
重力模型的基本公式可以表示为：T = k * (M1 * M2) / D^a，其中T表示两个地区之间的贸易、人口流动或投资的强度，k是一个常数，M1和M2分别表示两个地区的经济规模，D表示两个地区之间的地理距离，a是一个指数，用于衡量地理距离对联系程度的影响。
二、重力模型的应用过程1. 数据收集：首先需要收集关于地区经济规模、地理距离以及贸易、人口流动或投资强度的数据。
这些数据可以来自于统计局、商业机构、调查问卷等渠道。
2. 变量定义：根据研究的具体对象和目的，将收集到的数据转化为模型中的变量。
一般来说，经济规模可以用GDP、人口数量或其他相关指标表示，地理距离可以用实际距离或交通时间等方式衡量，贸易、人口流动或投资强度可以用贸易额、人口流动量或投资金额等指标表示。
3. 模型估计：根据收集到的数据和变量定义，利用计量经济学中的方法对重力模型进行估计。
传统的估计方法包括普通最小二乘法（OLS）、仪器变量法（IV）等。
4. 参数解释和检验：根据估计结果，解释模型中的参数。
一般来说，经济规模的系数表示经济规模对联系程度的影响，指数a的值表示地理距离对联系程度的影响。
为了确保参数的统计显著性，还需要进行假设检验。
5. 模型拟合度检验：为了评估模型的拟合程度和预测能力，一般需要计算模型的拟合度指标，比如决定系数（R-squared）等。
重力模型在未来交通分布预测中的应用
体 规划 确 定 各 小 区 未 来 人 口数
（ ） 由 交 通 调 查 得 到 各 小 区 工 作 岗 位 数 ． 根 据 市 区 用 ２
地 规 划 确 定 各 小 区 未 来 工作 岗位 数
（ ） 由 调 查 得 到 各 小 区 机 动 车 数 ， 根 据 人 口和 工 作 岗 ３
位 数 确 定 各 小 区 未 来 机动 车 数 （ ）利 用 重力 模 型 ，对 居 民 出行 ０ ４ Ｄ表 ．机 动 车 出行 ［ ） Ｄ
＝
Ｋ
易 （ 吖 （） １ ）
是
合 ，但 结 果 并 不 完 全理 想 ．主要 原 因 是 似 合 参 数 较 少 ．信 息 损 失 大 ，因 此 我们 作 了较 大 改 进 。 ２ 重 力模 型 在 丹 东 市 交 通 规 划 中 的 应 用
其 中 ｉ ｊ ｌ ，… … ，ｍ ， ＝ ，２ ，ｍ为交 通 小 区 的 个 数 ，
Ｙ （）＝ ｔ ｏ￣ （，１ ｏｘ （，２）＋ｔ （，３ ｉ ｏｏ ｔ ｉ ）＋ ｔ ｉ ＋ ２ ｏｘ ｉ ３
表 和 自行 车 出行 ０ Ｄ表 进 行 数 据 拟 合 ， 预 测 中采 用 丁 四 个影 响 因素 。
要 使 用 重 力模 型 ，必 须 先 对 交 通 的 生 成 量 和 吸 引 量 进 行 预 测 ，在 丹 东 市 的交 通 规 划 中 ，我 们 在 分 析 现 状 交 通 调 查 数 据 基 础 上 ，结 合 丹 东 市 总 体 规 划 ，市 区土 地 利 用 规 划 ，人 口 规划 ，选 择 了三 个 相 关 因 素 建 立 了 交通 生 成 模 型 。
了如 下 改 进 ， 骤 如 下 ： 步 （ ） 由 交 通 调 查 得 到 各 小 区现 状 人 口数 ，根 据 市 区 总 １
《土地利用与交通》课程报告-重力模型理论和应用
——重力模型的理论推导及应用方法
学号：
姓名：
专业：
课程教师：
目录
一、重力模型的理论推导- 1 -
二、重力模型的应用实例- 2 -
一、
1948年信息论创始人Shannon（香农）借用热力学中熵的概念提出了“信息熵”的概念：设A={a1,a2,… ,an}是某个信息系统X的符号表，A中各个符号出现的概率分布为{p1,p2,… ,pn}，则信息系统X的熵为：
K’2=( )-1=1.044183155
K’3=( )-1=1.067353957
于是完成第一遍迭代计算。再将求得的K’j（j=1,2 , 3）代入（1-12）式进行Ki的第二遍迭代。计算结果为：
K1=0.002551223K2=0.002965K3=0.002495651
再将新得到的K1、K2、K3代入式（1-13）第二次迭代求K’1、K’2和K’3
H（X）= - （1-1）
1957年Jaynes提出了最大信息熵原理：在所有备选的概率分布中挑选这样的分布作为信息系统X的分布，它在某些约束条件下，是信息熵达到最大值的分布。这是因为信息熵取最大值时对应的那组概率分布出现的可能性占绝对优势。后人又将最大信息熵原理推广到一般熵原理，作为决定系统概率分布的标准，广泛应用于气象学、交通科学等领域。一般熵原理就是：
首先令列约束系数K’1= K’2= K’3=1代入式（1-12）求两个行约束系数：
K1=( )-1=0.00242915
K2=( )-1=0.002823529
K3=( )-1=0.002376238
再将求得的K1、K2、K3代入式（1-13）中求K’1、K’2和K’3
K’1=( )-1=0.95215125
第三章 交通需求预测-重力模型
α的标定： 1、利用OD表统计Ti、Uj,确定tij 2、假定一α，并假定所有的Bj的初始值为1， 计算Ai 3、将得出的Ai代入公式求Bj，再将Bj代回求 Ai，直至前后两次计算结果大致相同 两个约 4、将所得Ai及Bj代入公式求出Xij，并检验所 得的交通的平均交通时间与实际交通的平均交 通时间之间的误差是否满足要求。
j
交通调整系数 kij及α的确定：先令 kij =1，此时公路局 模型同乌氏模型，用其计算α及 Xij，再通过下式计算 kij 。
l
1、计算时，将预测的发生量、吸引量及阻抗 参数代入模型，计算Xij。 2、通过迭代计算,使得
k ij = rij
其中:
1 − Yij 1 − Yij rij
′ X ij
Uj 1 T ′ = X ij × i + X ij 2 ∑ X ij ∑ X ij j i
3-34
l
l
l 【例题3.3-6】
按例题3.3-2中表3-17和表3-18给出 的现状OD表和将来发生与吸引交通量，以及表 3-14和表3-15给出的现状和将来行驶时间，试利 用重力模型和平均增长系数法，求出将来OD表。 设定收敛标准为 ε = 1% 。
采用行程时间函数作为出行阻抗函数最为常 见，以行程时间函数为出行阻抗函数的形式有 多种，常用的有以下三种:
α f ( t ij ) = t ij
① 幂函数： ② 指数函数： ③Gamma函数：
f ( t ij ) = exp( β t ij )
α f (t ij ) = t ij ⋅ exp( β t ij )
各类交通方式的预测
l l
1、自由类交通方式预测 建立步行与出行目的和出行距离的关系，即可 进行步行方式预测。 2、条件类交通方式预测 根据影响因素，这类交通方式的预测可采用先 预测车辆的拥有量，再预测其出行总比例，最 后预测各交通区之间的出行比例的程序。
4-2交通分布预测-重力模型
1、通过拟合现状OD，用最小二乘法确定待定系 数k,αβγ。 2、将预测的发生量Oi和吸引量Dj及阻抗代入无约 束重力模型公式，计算qij 3、若qij不满足约束条件，则用增长率法进行迭代 计算，使其满足约束条件。如迭代时，可采用下 列公式：
Dj 1 Oi ` q
7.5 重力模型的优缺点
缺点：
（1）重力模型最主要的缺陷是难以准确预测小区内出行 分布量。采用幂型、指数型函数时，由于区内出行距离很 小，预测结果往往比实际偏高。 解决这一问题的方法有两种：一是采用更为灵活的阻抗 函数型式，如复合型和离散型阻抗函数，但这会带来模型 标定上的困难；另外一种方法是区内出行不参加重力模型 的运算，而采用其他方法（如增长系数法）来单独处理。
Oi 16 28 40 84
表2 各区之间的行程时间
1 2 3
1 2 4 4
2 4 1 2
3 4 2 2
美国联邦公路局重力模型 模型形式为：
qij Oi
D
j
D j f cij kij
j
f cij kij
交通调整系数kij及α的确定：先令kij =1，此时公路局模型同乌 氏模型，用其计算α及qij，再通过下式计算kij 。
交通分布预测模型的选择交通分布预测模型的选择城市居民出行及市内货运分布预测城市居民出行及市内货运分布预测重力模型重力模型城市流动人口出行分布预测城市流动人口出行分布预测城市对外及过境客货运交通分布预测城市对外及过境客货运交通分布预测模型或平均增长率模型模型或平均增长率模型区域交通分布预测区域交通分布预测重力模型重力模型补充
7.5 重力模型的优缺点
优点：
模型形式直观，可解释性强，易被规划人员理解和 接 受； 能比较敏感地反映交通设施变化对出行的影响，适 用于中长期需求预测； 不需要完整的基年OD矩阵，如果有可信赖的模型参 数，甚至不需要基年OD矩阵； 特定交通小区（如新开发区）之间的分布量为零时， 也能进行预测。 能比较敏感地反映交通小区之间行驶时间变化的情况。
客流 重力模型
客流重力模型
客流重力模型是一种用来分析城市交通网络中客流分布的模型。
它基
于物理学中的万有引力定律和牛顿引力定律，将客流视为物理学中的
质点，通过考虑各节点之间的距离和吸引力，推断出客流的分布情况。
在客流重力模型中，节点之间的距离越近，吸引力就越大，客流就越高。
而同一个节点内的客流之间并没有特定的关系，每个客流都被视
为一个独立的质点，受到其他节点的吸引力所影响。
这种模型的重点
在于分析节点之间距离和吸引力的关系，从而推断出客流的分布情况。
在实际应用中，客流重力模型主要应用于城市交通规划和管理中。
根
据客流重力模型的分析结果，可以确定一些关键节点的客流分布情况，从而可以制定更科学和合理的交通规划策略。
例如，在规划城市公共
交通线路时，可以根据客流重力模型的分析结果，选择建设在客流最
集中的节点附近，以提高线路的运营效率。
除此之外，客流重力模型还可以用来分析不同节点之间的交通流动情况、分析城市发展趋势、评估城市的发展水平等。
通过分析客流分布
情况，可以更准确地把握城市的潜力和发展方向，为城市的长远发展
提供有力的支撑。
综上所述，客流重力模型是一种非常实用的工具，它可以帮助我们更好地了解城市的客流分布情况，规划和管理城市的交通网络。
未来随着城市的发展和人口的增长，客流重力模型的应用也将越来越广泛。
第五章交通运输模型
3 集计Logit选择模型
8、单击Model Fields used from view TAZ右侧的replace将I NC代替为income 9、单击Run，Output File Setting对话框将打开。可以接收 缺省的输出文件名或输入自己的文件名。单击OK来运行模 型。 10、选择File－Close All关闭所有的窗口。
2 重力模型出行分布
2 重力模型出行分布
8、单击Friction Factors 键来显示Friction Factors页。 9、在Purposes滚动列表中加亮HBW并做以下修改： • 单击Gamma单选按钮 • 输入28507作为a的值，输入0.02作为b的值，0.123作为 c的值 • 从Matrix File下拉列表中选择Shortest Path，并从Matrix 下拉列表中选择TRAVEL TIME。
1.2 预估小区回归分析模型参数
1.2 预估小区回归分析模型参数
7、单击OK显示Save Model As对话框。 8、输入“MYREG1”作为文件名，并单击Save键。 9、TransCAD预测模型，生成一个.MOD文件（称为MY REG1 .MOD）包含所有的回归方程参数。 10、单击Show Report，可以在文件底部查看预测的结果。 当完成后关闭Notepad程序。 11、选择File－Close All关闭dataview。
2 重力模型出行分布
2 重力模型出行分布
12、单击Ok来显示Save As对话框。 13、输入“My—graveval”作为文件名，并单击Save键。 TransCAD显示结果对话框说明三个模型都收敛了。 14、单击Close键。TransCAD显示Results Summary对 话框。 15、单击Close键。TransCAD显示结果矩阵。 16、选择File－Close All关闭地图和dataview。
交通规划 用简单引力模型预测运输需求分布
交通与汽车工程学院实验报告课程名称: 交通运输系统规划课程代码: 6011029 年级/专业/班学生姓名:学号:实验总成绩: 任课教师:开课学院: 交通与汽车工程学院实验中心名称: 汽车交通实验中心西华大学实验报告开课学院及实验室：交通与汽车工程学院实验时间：2013年12月19日2、实验设备、仪器及材料3、实验内容3.1 一般实验（非上机实验）：3.2 上机实验：，程序调试过程中出现的问题及解决方法3.2.4 程序运行的结果注解：理工科实验需记录实验过程中的数据、图表、计算、现象观察等，实验过程中出现的问题；其它如在计算机上进行的编程、仿真性或模拟性实验需记录程序核心代码以及程序在调式过程中出现的问题及解决方法；记录程序执行的结果。
4、实验总结4.1实验结果分析及问题讨论4.2实验总结心得体会注解：实验总结的内容根据不同学科和类型实验要求不一样，一般理工科类的实验需要对实验结果进行分析，并且对实验过程中问题进行讨论；在计算机上进行的编程、仿真性或模拟性实验需要对上机实践结果进行分析，上机的心得体会及改进意见。
其它实验应总结实验过程写出心得体会及改进意见。
说明：各门实验课程实验报告的格式及内容要求，请按照实验指导书的要求手工书写。
实验一用简单引力模型预测运输需求分布一、实验目的和任务能够运用简单引力模型预测未来年交通运输需求分布。
二、实验仪器、设备及材料计算机三、实验内容表1是只有三个小区的现状OD矩阵，表2是未来目标年各交通小区的出行产生量、吸引量，表3是交通小区之间的阻抗，用简单引力模型预测未来年运输需求分布矩阵。
表3 小区间交通阻抗（1）根据已知，把数据输入到Excel中，如图1所示：图1 （2）用简单引力模型：qij = k OiDj/ Rijγ1-1（3）两边取对数，得：log qij – log OiDj= log k –γlog Rij1-2（4）参数标定：令Y = log q ij– log O i D j , a = log k , b = - γ, X = log R ij（5）公式1-2可以转换成：Y = a + b X 1-3 （6）公式1-3为二元一次线性回归方程，a，b为待定系数。
重力模型的解释及系数计算方法
重力模型法(gravity model)是一种最常用的方法，它根据牛顿的万有引力定律，即两物体间的引力与两物体的质量之积成正比，而与它们之间距离的平方成反比类推而成。
下式为Casey(1955)提出的重力模型。
其中，：i，j小区的人口； d为i，j小区间的距离,α为系数。
上式的约束条件为：s.t.同时满足守恒条件的α是不存在的，因此，将重力模型修改如下：其中，为交通阻抗函数。
交通阻抗函数的几种形式：指数函数：（1）幂函数：（2）组合函数：（3）为参数。
单约束型B.P.R.模型其中，调整系数。
发生侧得到保证，即：以下以幂指数交通阻抗函数为例介绍其计算方法：第1步令m=0，m为计算次数。
第2步给出n（可以用最小二乘法求出）。
第3步令第4步求出第5步收敛判定。
若下式满足，则结束计算；反之，令m+1=m，返回第２步重复计算。
,作业：按上次作业给出的现状OD表和将来生成、发生与吸引交通量，利用下式重力模型和弗拉塔算法，求出将来OD表。
收敛标准。
重力模型:其中，，，。
读者也可以利用以前给出的现状分布交通量和表4-1示现状行驶时间，估计出这3个参数。
表4-1 现状行驶时间表4-2将来行驶时间解：利用重力模型求解分布交通量如下：同理，可以计算出其它各交通小区之间的交通量如下表所示。
重力模型的优点：a.直观上容易理解；b.能考虑路网的变化和土地利用对人们的出行产生的影响；c.特定交通小区之间的OD交通量为零时，也能预测；d.能比较敏感地反映交通小区之间行驶时间变化的情况。
重力模型的缺点：a.重力模型仅仅是将物理法则简单直观上容易理解；b.能考虑路网的变化和土地利用对地应用到社会现象，尽管有类似性，需要更加贴合人们出行的方法；c.一般，人们的出行距离分布在全区域并非为定值，而重力模型将其视为定值；d.交通小区之间的行驶时间因交通方式和时间段的不同而异，而重力模型使用了同一时间；e.求内内交通量时的行驶时间难以给出；f.交通小区之间的距离小时，有夸大预测的可能性；g.利用重力模型计算出的分布交通量必须借助于其它方法进行收敛计算。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
应用重力模型进行交通分布的详细步骤
第一步：求阻抗矩阵Rij（Impedance Matrix）
交通阻抗可表示为：出行距离和行程时间的长短，以及出行费用的大小等。
为真实地反映交通阻抗，依托工程公交规划采用通常使用的平均行程时间表示。
小区之间的阻抗——平均行程时间越小表示小区之间阻抗越小，越大表示小区之间阻抗越大，因此以平均行程时间为路权值求各小区之间的最短路径（Shortest Path），其值即为小区之间的阻抗R ij。
1、数据准备
（1）创建路网
图1表示的是TransCAD创建路网的界面。
（2）做选择集。
在Endpoints层，于dataview中选择质心点，将其作为一个选择集。
（3）各路段平均行程时间（Travel time）
其中，平均行程时间=Length/平均车速
2、操作过程
Networks/Paths—Multiple paths调出其对话框如图2所示。
3、运行结果（即为阻抗矩阵），如图3所示。
第二步：重力模型标定（校准）（Gravity Mode Calibration）1、数据准备
（1）公交基础OD矩阵。
（2）阻抗矩阵（Shortest Paths），如图3所示。
重力模型标定（校准）（Gravity Mode Calibration）数据准备：
基年OD矩阵的索引（质心层质心ID）与最短路径矩阵的索引（路网节点层质心ID）不匹配，并且因为下面将在路网节点层上操作，因此必须使基年OD 索引与最短路径矩阵的索引相一致，以使两表数据相对应（转换为“质心ID”）。
操作方法：按其对话框4示意操作。
2、操作过程
按对话框（如图5）操作即可。
3、运行结果
（1）标定参数结果：a=2.6288，b=0.2361，c=0.0，如图6所示，不过大看show report 里面参数更准确。
（2）K－Factor Flow：如图7所示。
第三步：创建综合阻抗因子f (Rij) (Synthetic Friction Factors)
1、数据准备
（1）创建空矩阵“Friction Factor shell”；
（2）已标定的a、b、c值；
（3）阻抗（最短路径）矩阵。
如图8所示。
2、操作过程
详见图9对话框所示。
3、运行结果
如图10所示。
第四步：应用重力模型（得2010年公交OD分布矩阵）1、数据准备
（1）已标定的a，b，c值；
（2）综合阻抗因子矩阵；
（3）阻抗（最短路径）矩阵。
（4）K－Factor矩阵
因为采用重力模型分布时要用到规划年交通出行量，所以必须在小区层上操作，因此综合阻抗因子矩阵、阻抗（最短路径）矩阵、K－Factor矩阵索引值必须与2010PCU_P、2010PCU_A的ID相匹配。
即要将这三个矩阵的ID转换回
来。
ID转换后的数据见图11。
2、操作过程
详见图12所示。
3、运行结果
规划年（2010年）公交OD矩阵（ID转换为实际小区号），详见图13。
操作结束，希望对大家有用！
我的QQ号是250234329.欢迎一起交流。