本发明涉及一种公交线路车辆调度中的时间控制点设置方法，属于公交线路车辆调度领域。
背景技术：
优先发展公共交通是缓解城市交通拥堵的重要措施。合理的公交线路车辆调度方案对于提高公交车辆准点程度以及服务水平具有重要意义，是提高公交出行方式吸引力的关键。公交线路车辆调度方案一般包括时间控制点设置、发车频率确定、时刻表制定三个环节。时间控制点设置是指从公交线路的各个停靠站中挑选出若干个重要的站点，而时刻表编制是确定公交车在各个时间控制点的发车时刻(驶离时刻)。因此，时间控制点设置是制定公交线路调度方案的基础。
根据公交车辆调度方法，如果公交车早于时刻表规定时间到达时间控制点，则需要在站点停车等待；如果晚于时刻表规定时间到达时间控制点，则在服务乘客结束后立即驶离站点。受路网交通流随机因素的干扰，同一运营时段内运行的各辆公交车在相邻时间控制点之间的行程时间必然存在差异，即存在公交车无法准点到达时间控制点的情况。如果各辆公交车在相邻时间控制点之间的行程时间差异较大，则导致部分公交车在时间控制点的停车等待时间较长，使得公交车在线路上的整体运行时间增加，严重降低公交车运行效率。
然而，现有的时间控制点设置都是依靠工程经验，根据各个时间控制点需相距较远、客流量较大等原则，从公交线路途径站点中选择若干个作为时间控制点，并未考虑时间控制点设置对公交车辆准点率以及公交车在线路上的总运行时间的影响，无法生成最佳的时间控制点设置方案，严重降低了公交车辆准点率与运行效率。
技术实现要素：
本发明的目的是为了解决现有的时间控制点设置未考虑时间控制点设置对公交车辆准点率以及公交车在线路上的总运行时间的影响，无法生成最佳的时间控制点设置方案，严重降低了公交车辆准点率与运行效率的缺点，而提出一种公交线路车辆调度中的时间控制点设置方法。
一种公交线路车辆调度中的时间控制点设置方法，包括：
步骤一、确定一条公交线路的时间控制点总数量；
步骤二、将公交线路的始发站与终点站设置为时间控制点；
步骤三、根据时间控制点的总数量以及公交线路的站点数量，生成i种时间控制点备选方案；
步骤四、针对第i种时间备选方案，计算最近连续10个工作日的给定运营时段内所有公交车在所有相邻时间控制点之间的行程时间方差的平均值之和，将计算结果称为第i种时间备选方案的评价指标；
步骤五、寻找数值最小的评价指标，并将其对应的时间备选方案中的时间控制点选为给定运营时段内该公交线路的时间控制点。
本发明的有益效果为：以所有公交车在相邻时间控制点之间的行程时间方差最小为原则提出时间控制点的设置方法，一方面可以提高线路上整体公交车的准点率，另一方面可以减少公交车辆在时间控制点的停车等待时间，减少公交车在线路上的总运行时间，提高公交线路运行效率，可减少乘客约7％的出行时间。
附图说明
图1为本发明的公交线路车辆调度中的时间控制点设置方法的流程图。
具体实施方式
具体实施方式一：本实施方式的公交线路车辆调度中的时间控制点设置方法，如图1所示，包括如下步骤：
步骤一、确定一条公交线路的时间控制点总数量。
步骤二、将公交线路的始发站与终点站设置为时间控制点。
步骤三、根据时间控制点的总数量以及公交线路的站点数量，生成i种时间控制点备选方案。
步骤四、针对第i种时间备选方案，计算最近连续10个工作日的给定运营时段内所有公交车在所有相邻时间控制点之间的行程时间方差的平均值之和，将计算结果称为第i种时间备选方案的评价指标。运营时段可以为早高峰等时段。
本步骤的计算过程的一个实施例为：
对于第一个工作日，计算在某运营时段内所有公交车从第1个时间控制点到第2个时间控制点的行程时间的方差，然后计算所有公交车从第2个时间控制点到第3个时间控制点的行程时间的方差，以此类推，直至计算出所有公交车从第n-1个时间控制点到第n个时间控制点的行程时间的方差，然后求这些方差的平均值，作为第1个工作日的计算结果，再对每一个工作日如此计算，将每一个工作日的计算结果累加起来(本实施例中为10个工作日)得到最终结果。
步骤五、寻找数值最小的评价指标，并将其对应的时间备选方案中的时间控制点选为给定运营时段内该公交线路的时间控制点。
具体实施方式二：本实施方式与具体实施方式一不同的是：
步骤一中，时间控制点总数量n通过如下公式计算：
其中，floor[·]为向下取整函数；m为公交线路的站点数量；a为中间变量。式中的20％根据实际需要也可以替换为其他值。
其它步骤及参数与具体实施方式一相同。
具体实施方式三：本实施方式与具体实施方式一或二不同的是：
步骤三中，时间控制点备选方案的种类数i通过如下公式确定：
其中：c为数学中的组合数；且
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四：本实施方式与具体实施方式一至三之一不同的是：
步骤四中，第i种时间备选方案的评价指标di具体通过如下公式计算：
其中，dj,j+1为在给定运营时段内，各个工作日的所有公交车从时间控制点j行驶至时间控制点j+1的行程时间方差的平均值。
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五：本实施方式与具体实施方式一至四之一不同的是：
步骤五中，数值最小的评价指标dk具体为：
dk＝min{d1,d2,…,di,…di}。
其它步骤及参数与具体实施方式一至四之一相同。
本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。