基于快速全局K-means谱聚类的城市交通控制子区的划分方法技术

本发明专利技术公开了一种基于快速全局K‑means谱聚类的城市交通控制子区的划分方法，涉及城市交通子区划分领域，其基本思路是根据路网内相邻信号交叉口之间的交通关联度作为划分依据，1）选定将要研究的路网，获取路网的道路拓扑结构、交通流数据，建立路网无向图；2）计算路网密度拉普拉斯矩阵L及其前k个最小的特征值对应的特征向量构建聚类特征矩阵；3）对聚类特征矩阵X实施快速全局K‑means聚类算法，将同一聚类里的交叉口划分到同一子区；4）判断进行聚类后子区内的交叉口数目｜V

Partition Method of Urban Traffic Control Subareas Based on Fast Global K-means Spectrum Clustering

The invention discloses a method for dividing urban traffic control sub-areas based on fast global K_means spectral clustering, which relates to the field of urban traffic sub-areas dividing. The basic idea of the method is to select the road network to be studied according to the traffic correlation degree between adjacent signal intersections in the road network. 1) To obtain the road topological structure and traffic flow data of the road network, and to establish the undirected road network. Fig. 2) Establish clustering feature matrix by calculating the Laplacian matrix L of road network density and the eigenvectors corresponding to its first k minimum eigenvalues; 3) implement fast global K means clustering algorithm for clustering feature matrix X, divide intersections in the same cluster into the same sub-region; 4) Judge the number of intersections in the sub-region after clustering.

【技术实现步骤摘要】
基于快速全局K-means谱聚类的城市交通控制子区的划分方法
本专利技术涉及城市交通子区划分领域，具体为一种基于快速全局K-means谱聚类的城市交通控制子区的划分方法。
技术介绍
随着社会经济的发展，城市机动车保有量在逐年增加，造成高峰时段城市路网中某些信号交叉口时常出现交通拥堵的现象，导致相应地平均出行时间上升，同时也可能会影响与之相邻的路段或者交叉口出现拥堵现象导致交通瘫痪，极大地影响了人们的出行和生活。城市路网因其不同的地理位置和特殊的交通特征，呈现出数量庞大、结构复杂等特征，对整个区域的路网进行信号协调控制其实现难度大且不具备合理性，已有的相关研究表明，将城市路网划分成可以进行协调控制的交通子区，能够充分发挥交通信号控制方案的作用，有效提高城市路网的运行效率。目前在城市交通子区划分领域，谱聚类方法的应用比较广泛，谱聚类后得到的交通子区内部交通流特征具有较高的相似性，而各子区之间的交通流特征具有较大的差异性；但谱聚类一般是与传统K-means方法结合使用，其聚类结果的准确性依赖于初始聚类中心的选择，且最终的聚类结果不是全局最优解；本专利技术将谱聚类与快速全局K-means聚类算法结合起来进行交通控制子区划分，快速全局K-means算法是采用一种增量式的方法确定聚类中心，不仅可以解决传统K-means聚类因随机选择原始聚类中心而导致聚类结果不准确的问题，而且也可以避免将噪音点作为聚类中心，减少计算时间复杂度，实现交通子区的快速划分，便于针对不同子区信号控制方案的实施，有效地缓解城市道路拥堵问题。
技术实现思路
本专利技术为了解决现有城市交通控制子区划分效率低下，不能实现最佳聚类结果的问题，提供了一种基于快速全局K-means谱聚类的城市交通控制子区的划分方法。本专利技术是通过如下技术方案来实现的：一种基于快速全局K-means谱聚类的城市交通控制子区的划分方法，包括以下步骤：1)选择将要研究的路网，建立路网无向图G＝(V,E)；1.1)选择待研究路网G，获取该路网的道路拓扑结构，在交叉口的进口道和出口道位置设计安装检测器，获取相邻交叉口i驶入交叉口j的车流分支数h，相邻交叉口i到交叉口j的平均行程时间T，相邻交叉口i到交叉口j的车流量Q，相邻交叉口i到交叉口j的车道数N；1.2)基于获取的道路拓扑结构以及实际交通流数据，计算相邻交叉口之间的交通关联度，相邻交叉口i和交叉口j之间的交通关联度计算公式如下：其中，wij是相邻交叉口i与交叉口j之间的交通关联度；If为路段流量不均衡系数；h为来自上游交叉口的车流驶入的分支数；为到达下游交叉口的车流量总和；Qmax为来自上游交叉口车流驶入分支中最大的车流量，为Qs中的最大值；T为车辆在两个相邻交叉口之间的平均行程时间；N为由上游交叉口驶向下游交叉口的车道数；1.3)将路网G内的信号交叉口用顶点V表示，相邻交叉口之间的路段用边E表示，则所研究的路网G可以用有权重的无向图G＝(V,E)表示，其中，V＝{V1,V2,...,Vn}是路网G中所有信号交叉口的集合，Vi表示第i个交叉口；E＝{E12,E23,...,Eij}是连接路网G中两个相邻信号交叉口的所有路段的集合，Eij表示连接相邻交叉口i与交叉口j之间的路段，其上的权重wij表示相邻交叉口i与交叉口j之间的关联度；wij的取值越大表示两个相邻交叉口之间的交通关联度越大；2)根据步骤1)中得到的相邻信号交叉口之间的关联度信息，计算路网G的密度拉普拉斯矩阵L，构建聚类特征矩阵X：2.1)对于路网G中任意一个信号交叉口i，交叉口i的加权度值di定义为与交叉口i相邻的所有交叉口关联度之和，即：用加权度值di构造一个对角矩阵D，对角矩阵D为路网G的度矩阵，其对角线上的元素就是各交叉口对应的加权度值di；2.2)拉普拉斯矩阵L可表示为L＝D-W，其中，矩阵W＝{wij}n×n为权值矩阵，其元素wij为相邻交叉口i与交叉口j之间的关联度；若两个信号交叉口是不相邻的，则其对应位置的元素为零；2.3)将拉普拉斯矩阵L标准化后得矩阵计算标准化后拉普拉斯矩阵L’的特征值λ＝{λ1,λ2,...,λn}，将特征值λi按正序排列，计算其前k个特征值对应的特征向量x1,x2,...,xk，构建聚类特征矩阵X＝{x1,x2,...,xk}；3)对聚类特征矩阵X＝{x1,x2,...,xk}进行快速全局K-means聚类，将属于同一聚类的信号交叉口划分到同一个控制子区内：3.1)将聚类特征矩阵X＝{x1,x2,...,xk}的行向量进行单位化得到矩阵Y；3.2)先将所有数据聚为一类，即K＝1时，计算矩阵Y中所有样本数据的均值作为第一类的初始聚类中心：并设置聚类簇的个数q＝1；3.3)计算第q(q≥2)类的聚类中心，以c1,c2,…,cq-1作为前q-1个初始聚类中心，计算下一个聚类中心cq，引入其中，表示yj到相距最近的聚类中心的距离，计算使得mn最大的样本点yn作为下一簇的初始聚类中心cq；3.4)结合3.3)中得到的q个聚类中心c1,c2,…,cq，将其作为K-means算法的初始聚类中心进行聚类分析，计算每个数据点yi和每个聚类中心ci之间的欧式距离，按照距离就近原则将样本yi划分到与其距离最小的聚类中，更新聚类中心直至误差平方准则函数收敛，得到聚类结果，将属于同一聚类的信号交叉口划分到同一个控制子区内；然后用数字1、2、…、n对划分后的交通子区进行唯一编号；3.5)q＝q+1，若聚类簇的个数q>K时，输出聚类划分结果；否则，继续转到3.3)；4)判断划分后的交通子区内信号交叉口的数目是否满足约束条件Mmin≤|Vi|≤Mmax：4.1)交通控制小区内部的交叉口太少会导致交通控制小区数量过多，交通控制小区内部的交叉口过多会降低交通信号协调优化的有效性和精度；根据划分后的情况统计当前交通子区内的交叉口总数|Vi|，判断子区划分后的信号交叉口数目|Vi|是否满足约束条件：Mmin≤|Vi|≤Mmax，其中，Mmin、Mmax为交通子区内所允许的最小、最大信号交叉口数目；4.2)若交通子区内的交叉口数目满足该约束条件时，则保存聚类划分结果；4.3)若不满足时，则要调整路网G内相邻控制子区划分的边界，即若控制子区内的交叉口数目|Vi|≤Mmin时，则将该子区随机合并至相邻的控制子区中；若控制子区内的交叉口数目|Vi|≥Mmax时，则使用规范分割准则将该子区划分为子区A与子区B，其目标函数为：其中，cut(A,B)表示子区A与子区B之间的割；表示子区A中的信号交叉口和整个路网G中的信号交叉口之间的权重值之和，每一次边界调整都要重新执行步骤4)，并更新边界调整后的控制子区编号，直至子区内交叉口数目|Vi|满足约束条件Mmin≤|Vi|≤Mmax为止。本专利技术所叙述的基于快速全局K-means谱聚类的城市交通控制子区的划分方法的基本思路是根据路网内相邻信号交叉口之间的交通关联度作为划分依据，采用快速全局K-means谱聚类的方法对路网G进行交通控制子区的划分，保证划分后的子区内交通流特征相似度高，不同子区间交通流特征差异大，减少了计算时间复杂度，提高了子区划分效率，解决目前交通控制子区划分效率低下的问题，便于针对不同交通子区信号控制方案的实施。本文档来自技高网...

【技术保护点】
1.一种基于快速全局K‑means谱聚类的城市交通控制子区的划分方法，其特征在于：包括以下步骤：1)选择将要研究的路网，建立路网无向图G＝(V,E)；1.1)选择待研究路网G，获取该路网的道路拓扑结构，在交叉口的进口道和出口道位置设计安装检测器，获取相邻交叉口i驶入交叉口j的车流分支数h，相邻交叉口i到交叉口j的平均行程时间T，相邻交叉口i到交叉口j的车流量Q，相邻交叉口i到交叉口j的车道数N；1.2)基于获取的道路拓扑结构以及实际交通流数据，计算相邻交叉口之间的交通关联度，相邻交叉口i和交叉口j之间的交通关联度计算公式如下：

【技术特征摘要】
1.一种基于快速全局K-means谱聚类的城市交通控制子区的划分方法，其特征在于：包括以下步骤：1)选择将要研究的路网，建立路网无向图G＝(V,E)；1.1)选择待研究路网G，获取该路网的道路拓扑结构，在交叉口的进口道和出口道位置设计安装检测器，获取相邻交叉口i驶入交叉口j的车流分支数h，相邻交叉口i到交叉口j的平均行程时间T，相邻交叉口i到交叉口j的车流量Q，相邻交叉口i到交叉口j的车道数N；1.2)基于获取的道路拓扑结构以及实际交通流数据，计算相邻交叉口之间的交通关联度，相邻交叉口i和交叉口j之间的交通关联度计算公式如下：其中，wij是相邻交叉口i与交叉口j之间的交通关联度；If为路段流量不均衡系数；h为来自上游交叉口的车流驶入的分支数；为到达下游交叉口的车流量总和；Qmax为来自上游交叉口车流驶入分支中最大的车流量，为Qs中的最大值；T为车辆在两个相邻交叉口之间的平均行程时间；N为由上游交叉口驶向下游交叉口的车道数；1.3)将路网G内的信号交叉口用顶点V表示，相邻交叉口之间的路段用边E表示，则所研究的路网G可以用有权重的无向图G＝(V,E)表示，其中，V＝{V1,V2,...,Vn}是路网G中所有信号交叉口的集合，Vi表示第i个交叉口；E＝{E12,E23,...,Eij}是连接路网G中两个相邻信号交叉口的所有路段的集合，Eij表示连接相邻交叉口i与交叉口j之间的路段，其上的权重wij表示相邻交叉口i与交叉口j之间的关联度；2)根据步骤1)中得到的相邻信号交叉口之间的关联度信息，计算路网G的密度拉普拉斯矩阵L，构建聚类特征矩阵X：2.1)对于路网G中任意一个信号交叉口i，交叉口i的加权度值di定义为与交叉口i相邻的所有交叉口关联度之和，即：用加权度值di构造一个对角矩阵D，对角矩阵D为路网G的度矩阵，其对角线上的元素就是各交叉口对应的加权度值di；2.2)拉普拉斯矩阵L可表示为L＝D-W，其中，矩阵W＝{wij}n×n为权值矩阵，其元素wij为相邻交叉口i与交叉口j之间的关联度；若两个信号交叉口是不相邻的，则其对应位置的元素为零；2.3)将拉普拉斯矩阵L标准化后得矩阵计算标准化后拉普拉斯矩阵L’的特征值λ＝{λ1,λ2,...,λn}，将特征值λi按正序排列，计算其前k个特征值对应的特征向量x1,x2,...,xk，构建聚类特征矩阵X＝{x1,x2,...,xk}；3)对聚类特征矩阵X...

【专利技术属性】
技术研发人员：闫飞，王昆，阎高伟，张曼，丁洁，
申请(专利权)人：太原理工大学，
类型：发明
国别省市：山西,14