【毕业学位论文】（Word原稿）基于出租车GPS数据分析深圳交通拥堵情况-统计教育学

基于出租车 析 深圳交通 拥堵情况 解放军信息工程大学 摘要： 本文基于出租车 据，从区域交通量 布角度出发， 结合深圳市交通地图， 研究 深圳市 交通小区的划分，建立载客出租车和人们出行的 析筛选出拥堵的路段时段以及拥堵的路口时段。 本文特点在于把经过处理的数据显示在地里信息系统软件上，并结合 马克威分析系统完成小区的划分以及交通拥堵模型的检验。 首先， 采用随机抽样的方法获取数据样本，并 利用C+编程对 数据进行优化 处理 ， 得到 出租车载客 起讫点 时刻的 据，基于聚类与交通小区划分的相 似性 原理 ， 利用马克威分析系统中的快速 聚类 方 法进行交通小区的划分 ，共求得 200个小区以及各小区的坐标 。 其次 ， 利用合理假设，建立交通量 布概率模型，采用 法 对模型的参数进行估计、对 模型的整体拟合效果进行了检验 ，利用马克威分析系统处理数据得出小区交通量，进而求解出交通量的 布。 再次 ， 根据各区域总体居民出行时耗 (出行距离 )与出行比例的总体分布，将交通小区发生量离散为不同出行时耗 (出行距离 )区间上的出行量，并将相应区间内的出行量按区间所包含的交通吸引小区的面积及吸引量大小进行分配，得到了人们出行时耗 特征图。 最后 ， 根据出租车载客后的行驶数据，编写程序筛选出有用的数据，并且把这些数据显示在 找车辆比较拥堵的路段以及拥堵的路口时段。通过马克威分析系统软件处理 据，得到各个路段的车辆数据。 关键字： 交通小区、 布概率模型 、聚类分析、马克威分析系统 一、 问题的提出 近年来，随着国民经济的快速增长，人流、物流、信息流以前所未有的密度涌向大中城市并向周边辐射，城市化进程明显加快，城市规模不断扩大，人口不断集中。此种状态在带动城市交通需求高速增长， 机动车辆快速增加的同 时，也促使 城市道路负荷加重，交通拥堵现象日益加剧。交通拥堵使交通延误增大，行车速度降低、时间损失、燃料费用增加、排污量增大、城市环境恶化、并诱发交通事故，直接影响人们的工作效率和身体健康。 出租车是城市客运交通的重要组成部分，是常规公共交通的重要补充。随着出租车数量的不断增长，出租车交通量已成为城市道路交通总量中的重要组成部分，成为影响城市道路交通分布及分配预测准确性的重要因素之一。由于出租车运行的随机性，没有固定的起讫点和运行线路，给定量分析带来一定的困难。 各大城市 越来越多的出租车 的安装了 端，这 些终端能够每隔 1分钟向出租车管理中心发送本车的位置、速度和方向等信息，是车辆 始数据主要保存出租车上装配的 些数据包括序号，车牌号码， 度，纬度，车辆状态 (空车、重车 )，车辆速度，车辆方向(8个方向 )等信息。 这些 现阶段在交通拥堵方面主要集中在研究某个路口的情况，并且用缓冲区方法 ，而没有对整个区域有个整体的描述，而 对于整个区域 的研究则是基于交通需求建立模型，数据采集的精度不高。本文基于 出租车 立 深圳整个区域的交通拥堵模型。 二、 数据的采集和预处理 我们从网络 （什么网） 上下载了深圳市出租车的 于数据量较大 ，硬件设备无法对这些数据全部进行处理 ，我们采取了随机抽样的方式，从一万多辆出租车中随机抽取了 100 辆出租车的 据进行建模分析。然而，原有的一些数据的经纬度超出深圳市范围，还有些车辆的状态明显错误） ，我们用 C+语言编程对数据进行了筛选处理，得到了有效的 用这些数据我们进行了下面的数学建模分析。 数据格式见表 1。 表格 1 数据格式 字段 示例 备注 车牌号 日期 时间 经度 纬度 运营状态 速度 方位角 粤 011/04/18 08:44:27 53 0 0 为空载； 1 为重载 该点车辆的瞬时速度 行车方向 三、 问题的分析 与模型的建立 划分小区 理数据 根据出租车载客的起讫点，结合深圳市的交通地图，恰当的划分交通小区，并选择小区中的某一点，用其经纬数值作为该小区的坐标 。 租车的状态为 0表 示载人，为 1表示空车。在间隔几分钟后得到的下一数据多和上一数据的状态相同。要得到出租车载客的起讫点，我们需要 对数据进行处理，使得到的数据中相邻的出租车运营状态 不相同。当运营状态由 0变为 1时表示出租车由空车到有客人上车，当运营状态由 1变为 0时表示出租车由载客到客人下车。处理前后的数据见表 2、表 3： 表格 2 原始数据 jd wd v 011/04/24 18:27:32 0 4 粤 011/04/24 18:28:27 6 0 粤 011/04/24 18:29:19 37 7 粤 011/04/24 18:29:58 23 0 粤 011/04/24 18:30:53 0 4 粤 011/04/24 18:31:48 8 3 粤 011/04/24 18:32:44 17 3 粤 011/04/24 18:33:38 22 4 粤 011/04/24 18:34:33 12 4 表格 3 优化后的数据 jd wd v 011/04/24 18:29:19 37 7 粤 011/04/24 18:29:58 23 0 粤 011/04/24 18:33:38 22 4 型建立 类算法 经过上一步的数据优化，即可得到由起讫点数据组成的数据集 。 现在需要对这些数据进行聚类 运算：一段时间内的起点或者讫点被划分为若干个区域，每个区域中的点分布紧凑，区域与区域之间自然分开 本文采用了 类算法 。是数据挖掘技术中一种经典的基于划分的聚类算法 。 其目的在于从大量数据点中找出具有代表性的数据点，即中心点，然后再根据这些中心点进行后续的处理 。 类过程如下： 首先从 过运算其它点与这些聚类中心点的相似度（距离），将其分别分配给与其相似度最高（距离最 近）的中心点所在的聚类；然后对划分好后的聚类重新运算聚类中心这一过程不断重复直到标准测度函数开始收敛 。 设 X=, 21 R 为待聚类样本的全体 (称为论域 )， 即道路网络 ， , 21 为观测样本 各个路段 )的特征矢量或模式矢量 ， 对应特征空间中的一个对象 ， 设 c 为聚类数 ， 为样本数 ， 聚类中心点集 , 21 ,且, 为硬划分矩阵。若按照最近距离进行样本划分，则样本硬划分矩阵计算如下： ,2,1;,2,10,m i ，其它， 式中，以类内平方误差和最小化为聚类目标函数，则聚类目标函数表示 为： m (用马克威分析系统聚类分析 利用马克威分析系统，通过聚类分析得到了 200个交通小区 的 坐标 ， 如图表1， 各个小区内的交通量和各小区之间的欧氏距离。 在 图表 2。 图表 1 交通小区坐标分布图 14 2 小区的大致区域 根据小区划分和出租车 据，给出载客出租车的 ：某时刻从坐标 (, ) ( , )、 ( , ) 的出租车有多少辆 。 域交通量 布概率模型 假设 某一小区交通发生量分布到不同小区的 与这些小区的吸引量成正比。 这是 许多分布模型的基础，不再解释。 整个深圳区域共划分为 m 个小区，则区域交通量 布矩阵为： . . . . . . . . . 示从小区 i 到小区 j 的交通 量。设区域交通为 , 则有 : ,总 并记：总i,j=1,2, ,m 式中 示从 i 小区到 j 小区的交通量与整个研究区域的交通量之比。 取随机变量 x 、 y 分别表示个交通的起点、终点，利用参数 ( x 、 y ) 的联合分布。 表格 4 联合分布 , 0 i=1,2, ,m 上式为 ( x , y) 率分布。 型求解 把条件概率 | 简记为 P( j / I ）用来表示区域 i 的交通生成量与分布到其它各 小区之间的关系。 由于 |,即： 要求解 得先求出 P x=i 与 pj/i。 为讨论方便 , 将 记为 P , 取得小区 i 的交通生成量为 行全样本抽样 , 设 小区 i 的交通生成量 , 则 然是随机变量 , 且服从二项分布 P ) , P =P( j/ i) i , j = 1 , 2 , , . 若仅仅考虑机会均等 , P 的估计量为 p=X/ Q , 这仅是一个均值 ,而对 j 区城和其它各小区的区别未能反映 ,过去的信息和运行状况未起作用 , 通过对 P 采用 计 , x y 1 , 2 , , , m x=i 1 2 m y=j , , , ,pm p1 就能克服上述缺点。 法认为参数 P 不是一个确定的常 t , 而被看作是一个随机变量 P 的一个值 , 因而符合客观实际。 型检验 根据已获得的关于 P 的信息 , 可建立准确程度不同的 P 先验分布 , 作为对 P 进行修正的后脸分布的依据。 如果所研究的各个小区的社会经济地理 特点没什么大的区别 , 各区的相对位置近似。基于这样的情况 , 如果对各区间的交通分布情况了解不多 , 可假设 P 的先验分布为均匀分布的 , 密度函数为 : 其他当,010,1)( P 为交通量， f(p) 为密度函数。 由概率和前述可知 , x 的概率函数也服从二项式分布 ,则有 : )1()/( 式中： X 是随机变量 ( X=0,1,2， ， ； 小区 I 的交通生成量。 所以 X , P 的联合分布为： 1)1()/()(),( X 的边缘分布为： )1(),()(10P 的后验分布为 : 10 11)(),( P 的 计量就是后验分布的均值 : 10 2111101 上式中的 通量后验分布均值 ;为变量 ;为交通生成量 。 用马克威分析系统处理数据 根据第一题有关划 分小区的情况，以及出租车 用马克威分析系统软件可以得出区域交通分布情况。 如图表 4、 5、 6。 表格 5 部分小区交通量以及所占的交通比例 小区 小区 1 小区 2 小区 3 小区 4 小区 5 小区 6 交通量 270 31 18 42 10 6 交通比例 表 3 部分小区交通量 表格 6 小区交通分布情况 j i 1 , 2 , , , m 总量 1 2 m , , , 求解交通量 布 对于某一固定小区来说 ,设 0 为 平均值 , 02为 方差 , 参数 P 显然服从正态分布 N ( 0 , 02 ) , 故 P 的先验概率密度函数为 : 20200 2ex )( 这样可做到综合考虑交通分布情况 , 避免随机因素的影响 。 最后转化为具体的交通量 布 , 由 P = 得 从而求得整个区域的交通量 布 。如表格 7。 表格 7 部分小区交通量 分布表 小区 1 小区 2 小区 3 小区 4 小区 5 小区 6 小区 7 小区 1 区 2 区 3 区 4 区 5 区 6 区 7 合理的假设条件下， 推断 人们出行的 空分布 。 型假设 本文提出的出行 阵估算方法中假定：城市一定区域内交通小区的出行时耗 (出行距离 )特征具有相似性，区域总体出行特征可近似代表各小区征，并将城市区域划 分为中心区、过渡区、外围区等 (视城市规模和空间布局可进一步细化为多个区域，区域划分越多，计算精度越高，但计算量也越大 )。 就假定条件而言，城市不同的区域其出行特征存在一定差异，但具有相同区位和特征的区域内的居民出行特征具有一定相似性，而且区域总体的出行特征正是各交通小区出性特征的综合， 图 7为 深圳市 调查结果，从图中可以看出各小区的出行比例与出行时耗关系曲线与区域总体关系曲线基本一致，并以 2x 检验证明了假设条件的成立。 型算法 本方法的基本思路：根据各区域 总体居民出行时耗 (出行距离 )与出行比例的总体分布，将交通小区发生量离散为不同出行时耗 (出行距离 )区间上的出行量，并将相应区间内的出行量按区间所包含的交通吸引小区的面积及吸引量大小进行分配。 具体计算思路如下图 4、 5所示，图 4中曲线 f(x)为小区 区 上的出行比例为图 5阴影面积 ,设小区 对应的区间 上出行量P；图 4为城市规划小区的空间分布，阴影为空间上小区 所覆盖的区域，然后根据阴影区域内各吸引小区的大致面积和吸引量将区间出行量P按比例进行分配。 图表 4 小区 i 时耗分布图 图表 5 时耗覆盖区域示例 按照上述思路，主要计算公式 如下： （ 1） 式中，i 到交通小区 j 的出行量 ; k 个时耗 (距离 )区间交通小区 j 的出行量；交通小区 i 所在区域出行分布特征中第 k 个时耗 (距离 )区间的出行比例， 且1k 交通小区 j 隶属于交通小区 i 的第 k 个出行区间的隶属系数，可由下式计算。 0,2m a x,2m 其它 2,2 （ 2） 其中，距离 )； 距离 )； 时耗 (距 离 )区间 限。 式 (2)实际反映了交通小区 j 在交通小区 i 出行时耗 (距离 )轴上覆盖区间的情况，由图 4 示例可以理解，图中交通小区 j 相对于交通小区 i 在出行时耗 (距离 )区间 15， 20、 20， 25、 25， 30、 30， 35上的隶属系数 别为： 0，易知式 (2)能保证1k 由计算过程可知，根据式 (1)计算得到的 阵能够满足矩阵发生总量约束ik j ，但不一定能够满足吸引总量约束 i q。因此上述结果还需要进行修正，可根据交通小区吸引量进行 公式类同式 (1)，如下： （ 3） 式中，j 所在区域出行分布特征中第 k 各时耗（距离）区间的出行比例； 交通发生区 i 隶属交通吸引区 j 的第 k 个出行时耗区间的隶属系数，其它符号意义同上。公式（ 2）变 为： 0,2m m 其它 2,2 （ 4）其中， 距离 )直径；其它符号意义同上。 按照式 (1)和式 (3)计算，分别得到满足出行总量约束和吸引总量约束的 其元素分别为多 最终计算 2 （ 5） 型修正 实例分 析验证，利用上式 (5)得到的 阵，基本能满足约束条件的精度要求，如果仍不能满足，利用下式 (6)进行迭代运算，逐步修正结果： 22)1()()1()()1()( （ 6） 式中， 1 k=1 时 05)所求得 阵； 阵所对应的交通 小区 1 q， )0( j 已求得的第 阵所对应的交通小区 j 的出行吸引量，即 1 q，且 0j 的出行吸引量 。根据式 (6)进行迭代运算，直到满足精度要求)1()()1()(近于 1为止。 图表 6 计算公式图解 图表 7 出行时耗特征 由图 7 可以看出人们出行的交通量随出行距离的增加先增长然后呈下降趋势，这与实际情况基本吻合。 根据出租车载客后的行驶数据，筛选出拥堵的路段时段以及拥堵的路口时段。 通拥堵分析 交通拥堵 既是一个经济问题，又是一个社会问题，已经成为现代城市经济发展过程中的瓶颈，备受人们的广泛关注。交通拥堵程度的判断，广大市民往往是凭经验和直觉进行的，如，“ 人车挤成一团 ”、“ 前看不到头，后看不到尾 ”等等，缺乏量化手段，具有一定的主观性和片面性。因此，对交通拥堵程度的量化及评价的研究具有重要意义和理论价值。 通拥堵的定义 交通拥堵的定义：某路段在某个时段平均行驶速度小于 15 公里 /小时为拥堵。根据出租车载客后的行驶数据，编写程序筛选出有用的数据，并且把这些数据显示在 找车辆比较拥堵的路段 以及拥堵的路口时段。 图表 8 红色表示拥堵的路段 通过马克威分析系统软件处理 到各个路段的车辆数据。 表格 8 深圳拥堵路段 拥堵路段 1 布吉关 2 沙湾关 3 南头关（瓶颈路段） 4 华强北 5 东门中路 6 清水河 7 深南路华侨城路段 8 深盐路盐田港区 9 盐坝高速入口 10 布吉路路段 11 深惠路路段 12 坂田五和大道 13 罗沙路 路段 14 107 国 道 设平均速度