本发明提供一种专用道公交干线协调子区划分方法，该方法先通过统计周期内的公交车辆通过下游交叉口的一次性通过率，并建立一次性通过率矩阵；再将一次性通过率矩阵进行欧几里得距离计算，获取其欧氏距离矩阵；利用k‑means聚类算法对欧式距离矩阵进行分类，根据分类结果进行交叉口的分区，最后依照分区结果，对分区进行相位差调整。
1.一种专用道公交干线协调子区划分方法，其特征在于，包括如下步骤：(1)计算公交车辆到达下游交叉口的预测时刻，若所述预测时刻处于下游路口绿灯时段内，则判定该公交车辆能够一次性通过下游交叉口，否则判定为不能一次性通过下游交叉口；
(2)统计一周期T内的公交车辆通过下游交叉口的一次性通过率，并建立一次性通过率矩阵λ＝[λ1 λ2 … λi],其中λi＝Ni/ni，其中，Ni为一周期内公交车辆在路段i的下游交叉口一次性通过的车辆数，ni为一周期内路段i的专用道公交车流量，矩阵λ为干线的交叉口一次性通过率矩阵；
(3)建立一周期内公交车辆通过所有交叉口的一次性通过矩阵λper为其中k＝T/0.5；
(4)将上述一次性通过矩阵λper进行伪F统计量分析，计算最佳分类类别数为PFS；
(5)将一次性通过率矩阵λ进行欧几里得距离计算，获取其欧几里得距离矩阵，具体如下：得到相应的欧几里得距离矩阵
(6)利用k-means聚类算法对欧几里得距离矩阵进行分类，共分为PFS类，所得到的分类结果如下：L＝[L1 L2 … Lm]
其中，L1 L2 … Lm为利用聚类算法对欧几里得距离矩阵进行分类后的分类项；
其中，S1 S2 … Sm为根据分类结果L得到的一次性通过率分类的分类项；
(7)选取一次性通过率 其中 判断 是否成立；若成立，则两个交叉口属于同一分区；否则执行(8)操作；
(8)判断 是否成立，若成立，则两个交叉口属于同一分区；否则两交叉口不在同一分区，其中Ckj为第j个交叉口的周期；
(9)判断此分类是否判断结束：j+1≤length(Sm)，若满足上式则依次选取下一元素进行(7)操作，否则，进行(10)操作；
(10)判断分类循环是否结束：m+1≤length(S)，若满足上式则依次选择下一类进行(7)操作，否则，停止循环；
(11)整理所得分区，依照分区结果，对一次性通过率高的分区进行相位差调整实现分区内干线协调，对于一次性通过率较低的分区通过采用公交主动或被动优先策略增加公交优先权。
一种专用道公交干线协调子区划分方法
技术领域
[0001] 本发明涉及一种专用道公交干线协调子区划分方法，具体涉及一种基于专用道公交一次性通过率的聚类分区方法。
背景技术
[0002] 交通控制子区划分是以路网供给和需求空间分布不均为研究，用以降低管理和控制的复杂程度，提高通行效率和系统稳定性。作为现代交通工程研究的重要内容之一，在城市交通规划和控制领域的研究中均有着较多的成果和应用。与社会交通相比，公交车辆受到运营调度、信号控制、公交站点等多重影响，因此，公交优先必然会对控制子区的划分提出新的要求。同样，公交优先进行干线协调的前提，需以公交干线协调的子区划分为基础，但是，目前交通子区划分基本都是针对社会车辆的控制子区，鲜有存在公交干线协调分区的研究。与社会交通相比，公交车辆受到运营调度、信号控制、公交站点等多重影响，因此，公交优先必然会对控制子区的划分提出新的要求。
发明内容
[0003] 为了解决上述问题，本发明提出了一种专用道公交干线协调子区划分方法，具体采用如下步骤：
[0004] (1)计算公交车辆到达下游交叉口的预测时刻，若所述预测时刻处于下游路口绿灯时段内，则判定该公交车辆能够一次性通过下游交叉口，否则判定为不能一次性通过下游交叉口；
[0005] (2)统计一周期T内的公交车辆通过下游交叉口的一次性通过率，并建立一次性通过率矩阵λ＝[λ1 λ2…λi],其中λi＝Ni/ni，其中，Ni为一周期内专用道公交车辆在路段i的下游交叉口一次性通过的车辆数，ni为一周期内路段i的专用道公交车流量，矩阵λ为干线的交叉口一次性通过率矩阵；
[0006] (3)建立一周期内干线公交车辆通过所有交叉口的一次性通过矩阵λper为[0007]
[0008] (4)将上述一次性通过矩阵λper进行伪F统计量分析，计算最佳分类类别数为PFS；
[0009] (5)将一次性通过率矩阵λ进行欧几里得距离计算，获取其欧氏距离矩阵，具体如下：
[0010]
[0011] 得到相应的欧氏距离矩阵
[0012] D＝[d1 d2…di-1]
[0013] (6)利用k-means聚类算法对欧式距离矩阵进行分类，共分为PFS类，所得到的分类结果如下：
[0014] L＝[L1 L2…Lm]
[0015] 依照此分类得到的一次性通过率分类如下：
[0016] S＝[S1 S2…Sm]
[0017] (7)选取一次性通过率 其中 判断 是否成立；若成立，则两个交叉口属于同一分区；否则执行(8)操作。
[0018] (8)判断 是否成立，若成立，则两个交叉口属于同一分区；否则两交叉口不在同一分区，其中Ckj为第j个交叉口的周期；
[0019] (9)判断此分类是否判断结束：j+1≤length(Sm)，若满足上式则依次选取下一元素进行(7)操作，否则，进行(10)操作；
[0020] (10)判断分类循环是否结束：m+1≤length(S)，若满足上式则依次选择下一类进行(7)操作，否则，停止循环；
[0021] (11)整理所得分区，依照分区结果，对一次性通过率高的分区进行相位差调整实现分区内干线协调，对于一次性通过率较低的分区通过采用公交主动或被动优先策略增加公交优先权。
附图说明
[0022] 图1是行程时间推算示意图。
[0023] 图2是分区判断流程图。
[0024] 图3是k-means聚类算法流程图。
具体实施方式
[0025] 为获取专用道公交车辆到达下游路口的具体时间，就需要对其路段行程时间进行分析判断。结合专用道公交车辆同社会车辆之间的行驶特性，两者的行程时间差异主要区别于公交车辆包含驻站时间一项而社会车辆不需要驻站，具体模型分析如下：
[0026] 如图1所示，通过分析每辆车的班次n、路段i到达时刻tiarr、到站时刻tiarrs、离站时刻tidep、到达下游路口停止线时刻tiarn
[0027]
[0028] 其中titp为路段i的站前行程时间，titb为路段i的站后行程时间。
[0029] 通过班次n的站前行程时间和站后行程时间统计，获得该班次在路段i的平均站前行程时间 和平均站后行程时间
[0030]
[0031] 其中m为该路段该班次的统计时段车辆数。
[0032] 通过分析班次n的前后车到站时间间隔tpar和该班次的后车驻站时间tdwell之间的回归关系，
[0033]
[0034] 经回归验证发现其回归参数形式如下：tdwell＝α6tpar6+α5tpar5+…+α1tpar+α0[0035] 因此可以获得班次n的路段i的预测总行程时间
[0036] 据此获得该车辆的到达下游交叉口的预测时刻
[0037] 根据预测结果连同下游路口绿灯时段[tgmin tgmax]判断，若 则判定该车能够不停车通过下游交叉口，否则判定为不能一次性通过下游交叉口。
[0038] 鉴于专用道公交车辆在实际运行过程中由于流量基数偏小和分布随机性大等问题，造成车辆在周期内分布不均匀，因此单周期的不停车通过率不具备代表性，因此选用30分钟为统计周期，其一次性通过率描述如下：λ＝[λ1 λ2…λi],其中λi＝Ni/ni。
[0039] 其中，Ni为专用道公交车辆在路段i的下游交叉口一次性通过的车辆数，ni为路段i特定时间内的专用道公交车流量，矩阵λ为干线的交叉口一次性通过率矩阵。
[0040] 因特定时段内的信号控制策略是相对不变的，所以假定时段时长为时间T(30min的整数倍)，则该时段的干线所有交叉口的一次性通过矩阵为
[0041] 其中k＝T/0.5
[0042] 为了确定上述一次性通过率最佳分类类别数目，本专利以一次性通过矩阵为输入，通过伪F统计量算法衡量几种分类情况下的类间及类内关系，进而通过其最佳FPS值确定最佳分类类别数，其具体步骤如下：
[0043] 将上述一次性通过矩阵λper进行伪F统计量分析(FPS)，其i个样本的离差平方和为[0044] 其中，Wt为分类t下的样本离差平方和， 为分类t下的样本均值；
[0045] 因此所有k个分类下各自样本的离差平方和的和为：
[0046] 其类间偏差平方和为：
[0047] 故其伪F统计量，即最佳分类类别数为： 其中n为样本总量，k为类别数。
[0048] 为加快聚类分析算法的运算速度，需要经原有的将一次性通过率矩阵进行欧几里得距离计算，获取其欧氏距离矩阵，具体如下：
[0049]
[0050] 得到相应的欧氏距离矩阵
[0051] D＝[d1 d2…di-1]
[0052] 利用k-means聚类算法对欧式距离矩阵进行分类，共分为PFS类，所得到的分类结果如下：
[0053] L＝[L1 L2…Lm]
[0054] 依照此分类得到的一次性通过率分类如下：
[0055] S＝[S1 S2…Sm]
[0056] 其中k-means聚类算法流程见附图3。
[0057] a.选取一次性通过率 其中 判断 是否成立；若成立，则两个交叉口属于同一分区；否则执行下步操作。
[0058] b.判断 是否成立，若成立，则两个交叉口属于同一分区；否则两交叉口不在同一分区。其中Ckj为第j个交叉口的周期。
[0059] c.判断此分类是否判断结束，即j+1≤length(Sm)，若满足则依次选取下一元素进行a操作，否则，进行下步操作。
[0060] 判断分类循环是否结束，即m+1≤length(S)，若满足则依次选择下一类进行a操作，否则，停止循环。
[0061] 整理所得分区，依照分区结果，对一次性通过率高的分区进行相位差调整等策略实现分区内干线协调，对于一次性通过率较低的分区通过采用公交主动/被动优先策略增加公交优先权。
