1.一种基于大数据的公交线路时刻表协同优化方法，其特征在于，步骤如下：
一.公交线路运行数据和客流数据获取
(1)公交线路运行数据获取
(1.1)判断公交车辆运行方向
公交线路运行区分上行和下行两个方向，公交车辆按照排班顺序在两个方向循环往复；公交线路的上行和下行为独立的运行线路、停靠站点和运行时段，公交车辆运行状态与客流也是相对独立的，分为两个运行方向的数据；
(1.2)车辆到站数据获取
建立车辆到站数据表，包含站点编号S_ID、车辆到站时间S_TIME、车辆的车牌号NUMBER和车辆运行方向DIRECTION共四个字段，见表1；
表1 车辆到站数据表
在车辆到站数据表中记录所有上行方向的车辆到达站点S的数据，步骤如下：
(1)提取公交GPS数据中车辆运行方向为“上行”的所有数据；
(2)先按照车辆的车牌号对步骤(1)中提取的公交GPS数据进行排序，再根据车辆到站时间对同一车辆的车牌号进一步排序；在连续同一车辆的车牌号中，前一条公交GPS数据匹配到路段L1，后一条公交GPS数据匹配到路段L2，则L1和L2的中间值作为车辆到达站点S的时间；
(3)将站点S编号、车辆到达站点S的时间、车辆的车牌号以及车辆运行方向写入对应的车辆到站数据表中，直至处理完所有上行方向的公交GPS数据，得到车辆到站数据表；
上行方向和下行方向的车辆到站数据获取的方式相同；
(1.3)车辆行程时间获取
建立车辆行程时间数据表，包括车辆的车牌号NUMBER、运行方向DIRECTION、首站到站时间O_TIME、末站到站时间E_TIME和行程时间INTERVAL五个字段，见表2；
表2 车辆行程时间数据表
获取上行线路车辆行程时间，步骤如下：
(1)在车辆到站数据表中提取车辆运行方向为“上行”的所有数据，先按照车辆的车牌号进行排序，再根据同一车辆的车牌号的车辆到站时间进一步排序；
(2)针对每一条车辆到站时间数据，若站点编号与车辆运行方向的首站一致，则该条数据为首站发车数据，若站点编号与车辆运行方向的末站一致，则该条数据为末站到站数据；
将车辆的车牌号、车辆运行方向、首站到站时间和末站到站时间分别写入对应的车辆行程时间数据表中；
(3)计算同一车辆的车牌号的首站到站时间与末站到站时间的间隔，并将其写入车辆行程时间数据表的行程时间中，直至遍历所有上行方向的车辆到站时间数据，得到车辆行程时间数据表；
上行方向和下行方向的车辆行程时间获取的方式相同；
(2)客流数据获取
(2.1)乘客上车站点获取
轨道交通刷卡数据中记录乘客进站和出站的站点信息，直接判断乘客轨道交通出行的上车站点和下车站点；
公交刷卡数据不记录乘客上车站点信息，判断乘客公交出行的上车站点时间，对公交刷卡数据中的刷卡时间和车辆的车牌号字段与车辆到站数据表中的车辆到站时间和车辆的车牌号字段进行匹配融合；
建立Mix_IC数据表，包括IC卡记录编码ID、线路名称LINE、车辆的车牌号NUMBER、线路运行方向DIRECTION、上车站点编号U_ID、上车站点名称U_NAME、上车时间U_TIME、下车站点编号D_ID、下车站点名称D_NAME和下车时间D_TIME共十个字段，见表3；
表3 Mix_IC数据表
车辆到站数据与公交刷卡数据匹配融合，步骤如下：
(1)提取任一条公交刷卡数据，将公交刷卡数据中ID、LINE和NUMBER三个字段的数据对应写入Mix_IC数据表的ID、LINE和NUMBER字段中；
(2)在车辆到站数据表中匹配与步骤(1)提取的公交刷卡数据NUMBER字段一致，且车辆到站时间与刷卡时间字段最接近的一条记录，将此条记录的站点编号和车辆到站时间字段对应写入Mix_IC数据表的上车站点编号和上车时间字段中，直至遍历所有公交刷卡数据；
同一次轨道交通出行的数据分别记录在“进站”和“出站”两条记录中，经处理后存入Mix_IC数据表中，步骤如下：
(1)先按照IC卡记录编码对轨道交通刷卡数据进行排列，再对同一个IC卡记录编码的记录时间进一步排列；
(2)提取步骤(1)两条连续的轨道交通刷卡数据，其IC卡记录编码和线路名称字段都相同，前一条数据的线路运行方向为“进站”，后一条数据的线路运行方向 为“出站”；将IC卡记录编码和线路名称的数据对应写入Mix_IC数据表的IC卡记录编码和线路名称字段中；
(3)将前一条数据的进站站点名称和进站刷卡时间对应写入Mix_IC数据表的上车站点名称和上车时间字段中，将后一条数据的出站站点名称和出站刷卡时间对应写入Mix_IC数据表的下车站点名称和下车时间字段中；
(4)按照站点编号、上车站点名称和下车站点名称，匹配进出站站点编号，并写入Mix_IC数据表的上车站点编号和下车站点编号字段中，直至遍历所有轨道交通刷卡数据；
(2.2)乘客下车站点获取
通过分析居民的乘车行为，基于家庭的公交出行人数占总出行人数的比例较高，即大部分乘客一日出行是以家庭为起点，又以家庭为终点，其出行活动链构成一个封闭的环形；
对于同一天使用公交IC卡超过2次的乘客，其乘坐地面公交的下车站点根据出行活动链进行判断；
使用Mix_IC数据表记录的数据进行下车站点判断，并将获取的下车站点数据写入Mix_IC数据表的下车站点编号字段，乘客下车站点的判断，步骤如下：
(1)在Mix_IC数据表中提取同一ID的数据，并按上车时间进行排列；
(2)判断同一ID的前后两条数据是否为同一条线路：若判断为同一条线路，则后一次刷卡上车的站点为前一次刷卡下车的站点；若判断为不同线路，则根据后一次刷卡上车站点匹配前一次刷卡线路中距离最近的换乘站点，作为前一条数据的下车站点；
(3)同一ID首次刷卡的上车地点为最后一次刷卡的下车地点，直至遍历Mix_IC数据表中所有数据；
(2.3)乘客换乘行为判断
判断乘客连续两次采取公交出行是否为换乘的依据：前一次乘车下车时间与后一次乘车上车时间的间隔是否超过设定的换乘时间上限，设定的换乘时间上限为30分钟；
首先对乘客下车时间进行匹配，步骤如下：
(1)提取Mix_IC数据表中的一条数据，以Mix_IC数据表中车辆的车牌号、线路运行方向和下车站点编号字段与车辆到站数据表中车辆的车牌号、线路运行方向和站点编号字段进行匹配，得到该条件下的车辆到站时间；一天中车辆在上下行方向多次往返运行，车辆到站数据表在该下车站点名称应有多个车辆到站时间；
(2)以Mix_IC数据表中上车时间字段，选择大于该上车时间且时间间隔最短的车辆到站时间，作为下车时间；将匹配的下车时间写入Mix_IC数据表的下车时间字段，直至遍历Mix_IC数据表中所有数据；
判断乘客连续两次刷卡记录是否属于换乘，步骤如下：
(1)先按照IC卡记录编码对Mix_IC数据表进行排列，再根据上车时间字段对同一IC卡记录编码进一步排序；
(2)提取相同ID的两条连续记录，将后一条数据的上车刷卡时间与前一条数据的下车刷卡时间作差，若时间间隔小于30分钟，则认为是换乘；遍历Mix_IC数据表中所有数据；
二.模型建立
模型假设
(1)将地面公交和轨道交通的上下行线路看做两条不同的线路，即同一条线路仅有一个运行方向，假设所有车辆严格按照公交时刻表从首站发出，忽略道路交通条件对车辆运行的影响，车辆在各站间行程时间分布已知；同一时段车辆的发车间隔相同，同一条线路的公交车辆数量固定，车辆型号均相同；公交车辆按照发车顺序运行，不存在串车和滞站现象；非换乘乘客到达站点的过程服从均匀分布，且同一时段乘客到达率相同；乘客在换乘站点的换乘步行时间已知，且换乘步行时间小于地面公交和轨道交通的发车间隔；当车辆到达站点后乘客即上车，忽略由于车辆拥挤原因继续候车的情况；线路运营车辆均为全程车，忽略区间车调度方案和大站快车调度方案；
(2)参数定义
S——地面公交站点序列，s∈S；
——换乘站点序列，
——换乘站点数
i
B——在研究时段发车且车次为i的地面公交；
Rj——在研究时段发车且车次为j的轨道交通；
eB——研究时段内地面公交线路B的发车数量；
eR——研究时段内轨道线路R的发车数量；
a——研究时段时长，a＝Tend-Tstart；
hB——研究时段内地面公交线路B的发车间隔；
hR——研究时段内轨道线路R的发车间隔；
hB,min，hB,max——地面公交的最小和最大发车间隔；
hR,min，hR,max——轨道交通的最小和最大发车间隔；
——研究时段内在站点s乘坐地面公交线路B的乘客到达率；
——在换乘站点 从地面公交站点到轨道交通站点的步行时间；  ——地面公交Bi到达换乘站点 的时刻；
——轨道交通Rj到达换乘站点 的时刻；
——研究时段内车次序列为1的地面公交车辆从首站开出的时刻；
——研究时段内车次序列为1的轨道交通车辆从首站开出的时刻；
——地面公交线路车辆从首站开至换乘站点 的行程时间；
——轨道线路车辆在研究区间从首站开至换乘站点 的行程时间；
tB——研究时段内地面公交车辆行程时间；
tR——研究时段内轨道线路车辆在研究区间的行程时间；
——地面公交线路首站与换乘站点 的站距；
——轨道线路首站与换乘站点 的站距；
LB——地面公交线路的全程长度；
LR——轨道线路在研究区间的长度；
Uc——单位乘客候车时间的时间价值，元/分钟；
Up——地面公交线路单位发车次数的运营成本，元/公里；
CB——研究时段内无换乘的乘客乘坐地面公交的候车时间成本；
CBR——研究时段内乘客从地面公交换乘至轨道交通的换乘候车时间成本；
CRB——研究时段内乘客从轨道交通换乘至地面公交的换乘候车时间成本；
CP——研究时段内公交运营企业的成本；
CS——研究时段内系统总成本；
(3)约束条件
根据模型假设，乘客在换乘站点的换乘步行时间小于地面公交和轨道交通的发车间隔，则对所有
地面公交与轨道交通的发车间隔应不小于最小发车间隔，同时不大于最大发车间隔，则有：
hB,min≤hB≤hB,max   (3.2)
hR,min≤hR≤hR,max   (3.3)
(4)目标函数
研究时段内地面公交发车数量eB与发车间隔hB的关系为：
研究时段内轨道交通发车数量eR与发车间隔hR的关系为：
1)非换乘的乘客乘坐地面公交的候车时间成本
乘客到达站点s乘坐在研究时段内发车且车次为i的地面公交Bi，其候车成本可表示为：
因此非换乘的乘客的候车时间成本可表示为：
2)乘客从地面公交换乘至轨道交通的换乘候车时间成本
在研究时段内，车次为i的地面公交Bi到达换乘站点 的时刻 可表示为：
其中地面公交线路车辆从首站开至换乘站点 的行程时间 可用式4.6计算：
车次为j轨道交通Rj到达换乘站点 的时刻可表示为：
其中轨道线路车辆从首站开至换乘站点 的行程时间 可用式4.8计算：
一位乘客在换乘站点 从地面公交Bi换乘至轨道交通Rj需要的步行时间为  其换乘候车时间可表示为：
其中 为布尔变量，
乘坐在研究时段内发车的地面公交并在换乘站点 下车的乘客总人数为 在换乘站点 由地面公交换乘至轨道交通的换乘率为 则乘坐在研究时段内发车的地面公交在换乘站点 下车换乘至轨道交通的乘客总人数可按式4.11计算
将在换乘站点 从地面公交换乘至轨道交通的乘客平均分配至每辆地面公交车次，则乘坐地面公交Bi且在换乘站点 换乘至轨道交通的乘客数量可表示为：
由于存在换乘步行时间 不超过地面公交发车间隔的假设，所以不会出现从k时段发车的地面公交换乘至k+1时段发车的轨道交通的情况；因此在换乘站点 由地面公交换乘至轨道交通的乘客换乘候车的时间成本为：
在研究时段在所有换乘站点乘客从地面公交换乘至轨道交通的换乘候车时间成本为：
其中
3)乘客从轨道换乘至地面公交的换乘候车时间成本
一位乘客在换乘站点 从轨道交通Rj换乘至地面公交Bi需要的步行时间为  其换乘候车时间可表示为：
其中 为布尔变量，
在换乘站点 上车乘坐在研究时段内发车的地面公交的乘客总人数为 在换乘站点由轨道交通换乘至地面公交的换乘率为 则在换乘站点 上车换乘至在研究时段内发车的地面公交的乘客总人数可按式4.18计算：
将在换乘站点 从轨道交通换乘至地面公交的乘客平均分配至每辆地面公交车次，则在换乘站点 由轨道交通换乘至地面公交Bi的乘客数量可表示为：
由于存在换乘步行时间 不超过轨道交通发车间隔的假设，所以不会出现从k时段发车的轨道交通换乘至k+1时段发车的地面公交的情况；因此在k时段在换乘站点 由轨道交通换乘至地面公交的乘客换乘候车的时间成本为：
在研究时段在所有换乘站点乘客从轨道交通换乘至地面公交的换乘候车时间成本为：
其中：
4)公交企业的运营成本
公交企业的运营成本与线路全程长度和研究时段内的发车数量有关，可用式4.23表示：
CP＝UPLBeB  (4.23)
5)系统总成本
研究时段系统总成本包括：非换乘的乘客乘坐地面公交的候车时间成本、乘客从地面公交换乘至轨道交通的换乘候车时间成本、乘客从轨道交通换乘至地面公交的换乘候车时间成本和公交企业的运营成本，即：
CS＝CB+CBR+CRB+CP  (4.24)
目标函数为系统总成本最小，即minCS；决策变量为地面公交的发车间隔hB。