本发明涉及一种基于时空特性分析的短时交通流预测方法，属于城市智能交通管理技术领域。
背景技术：
随着城市化进程加快，智能交通管理技术愈发广泛的应用于城市道路交通治理和车辆运输疏导、居民出行规划等实际工程当中，使得交通运输和居民出行逐步走向智能化，动态化和信息化。而该技术的核心基础部分便是短时交通流量的动态分析与预测，此项技术在交通工程领域发挥了不可或缺的作用。实时交通流量的精准预测为交通管理和指导提供了很大的便益，使智能交通系统充分发挥出其高效，方便，快捷的现实作用，智能选择合适路径，规避交通拥塞，减少居民用户出行损耗，最大程度减缓交通环境压力，合理分配道路资源。
短时交通流预测技术是要预测目标时段，路段或路径上未来的交通状况，时间间隔短，要求预测值与实际流量差值尽可能小，便于车辆管理和调控。此项技术话题一直十分活跃却难以取得突破性的进展，原因在于交通流过程是时变的，不同的空间位置、道路环境、时间段下，其状态特征存在巨大差异，难以准确把握。所以此项技术仅凭单一理论预测方法无法成功，需在把握交通流量内在联系与特性的基础上，建立高效的预测模型，方能取得不错的实际预测效果。
迄今为止，世界各国学者致力于研究短时交通流量预测技术，通过努力，斩获众多理论成果，并将之运用于实际工程，不断推陈出新。包括如kalman滤波法、时间序列法、神经网络法、支持向量机等。由于基于单一模型的预测方法都存在缺陷，国内外学者研究了很多基于组合理论的预测方法，如灰色神经网络预测方法、基于小波包和最小二乘支持向量机的预测方法、基于蚁群优化支持向量机的预测方法等。然而，交通流预测质量主要取决于两个问题：一、数据来源的质量问题，二、交通流预测模型的缺陷问题。
技术实现要素：
本发明所要解决的技术问题是提供一种基于时空特性分析的短时交通流预测方法，通过提高实验数据质量和优化参数改善预测模型，不仅提高预测模型的可行性和有效性，而且提高短时交通流预测的精确度，实现智能交通系统实时准确的交通流预测。
本发明为了解决上述技术问题采用以下技术方案：本发明设计了一种基于时空特性分析的短时交通流预测方法，用于实现目标道路的交通流量预测，包括如下步骤：
步骤a.获得目标道路指定历史时间区间内的实时交通流量，然后进入步骤b；
步骤b.针对目标道路指定历史时间区间内实时交通流量，进行预处理操作，获得目标道路指定历史时间区间内各个时间节点对应的交通流量，然后进入步骤c；
步骤c.将目标道路所对应指定历史时间区间内各个时间节点的交通流量，作为各个数据样本点，并针对各个数据样本点，进行维度空间坐标转换，获得各个数据样本点分别所对应的维度空间坐标，然后进入步骤d；
步骤d.基于各个数据样本点分别所对应的维度空间坐标，针对各个数据样本点进行聚类处理，获得各个聚类簇，然后进入步骤e；
步骤e.根据各个聚类簇，获得目标路段分别对应各聚类时间节点的交通流量，然后进入步骤f；
步骤f.基于目标路段分别对应各聚类时间节点的交通流量，结合预设预测时长跨度，获得各组训练样本(xk,yk)，然后进入步骤g；其中，k∈{1,…,k}，k表示训练样本的组数，xk、yk均属于目标路段分别对应各聚类时间节点的交通流量，xk所对应聚类时间节点与yk所对应聚类时间节点之间的差值等于预测时长跨度，且按时序、xk所对应聚类时间节点在yk所对应聚类时间节点之前；
步骤g.采用各组训练样本(xk,yk)，基于xk为输入，yk为输出，针对预设分类器进行训练，获得对应预测时长跨度的目标路段短时交通流预测模型，然后进入步骤h；
步骤h.基于所获目标路段短时交通流预测模型，针对目标路段，实现预测时长跨度的交通流量预测。
作为本发明的一种优选技术方案：所述步骤b中，采用历史平均法、相邻数据平均法或加权平均法中的一种或几种任意组合，针对目标道路指定历史时间区间内实时交通流量进行预处理操作。
作为本发明的一种优选技术方案：所述步骤b中，针对目标道路指定历史时间区间内实时交通流量，进行预处理操作，包括丢失值处理、错误值处理和异常值处理。
作为本发明的一种优选技术方案：所述步骤d中，基于各个数据样本点分别所对应的维度空间坐标，针对各个数据样本点进行optics聚类处理，获得各个聚类簇。
作为本发明的一种优选技术方案：所述步骤d中，初始化用于存放本次循环中核心点、及其直接密度可达点的集合ⅰ，以及包含各核心点分别所对应簇的集合ⅱ，然后基于各个数据样本点分别所对应的维度空间坐标，执行如下步骤,针对各个数据样本点进行聚类处理，获得各个聚类簇；
步骤d1.判断是否存在非核心点、以及未被划簇的数据样本点，是则进入步骤d2；否则进入步骤d10；
步骤d2.针对非核心点、以及未被划簇的各个数据样本点，任选一个数据样本点，作为本次循环的核心点，并进入步骤d3；
步骤d3.将该核心点划分至集合ⅱ当中，并获得该核心点所对应的所有密度可达的数据样本点，然后进入步骤d4；
步骤d4.针对该核心点所对应的所有密度可达的数据样本点，判断其中直接密度可达的数据样本点是否在集合ⅱ当中，是则进入步骤d5；否则步骤d6；
步骤d5.分别针对该核心点所对应直接密度可达的各个数据样本点，判断数据样本点对应的可达距离是否小于预设半径，是则将该数据样本点划分至该核心点所对应的簇中，并更新集合ⅰ中该核心点所对应的簇，且返回步骤d1；否则将数据样本点直接删除，并返回步骤d1；
步骤d6.针对该核心点所对应的所有密度可达的数据样本点，按可达距离升序，针对各个数据样本点进行排序，划入该核心点所对应的簇，并加入集合ⅰ，然后进入步骤d7；
步骤d7.由该核心点所对应簇中，剪切出第一个数据样本点，作为本次循环的待判断数据样本点，并进入步骤d8；
步骤d8.判断该待判断数据样本点是否属于集合ⅱ，是则将该待判断数据样本点，作为核心点，并返回步骤d3；否则进入步骤d9；
步骤d9.判断集合ⅰ是否为空，是则返回步骤d1；否则返回步骤d7；
步骤d10.即完成针对各个数据样本点的聚类处理，获得各个聚类簇。
作为本发明的一种优选技术方案：所述步骤e中，根据各个聚类簇，获得目标路段分别对应各聚类时间节点的交通流量，并针对各个交通流量进行数据归一化处理，更新目标路段分别对应各聚类时间节点的交通流量，然后进入步骤f。
作为本发明的一种优选技术方案：所述步骤g中，采用各组训练样本(xk,yk)，基于xk为输入，yk为输出，针对svm分类器进行训练，获得对应预测时长跨度的目标路段短时交通流预测模型。
作为本发明的一种优选技术方案：还包括步骤gh，执行完步骤g之后，进入步骤gh，执行完步骤gh之后，进入步骤h；
步骤gh.采用细菌觅食算法，针对目标路段短时交通流预测模型中的参数进行优化。
作为本发明的一种优选技术方案，所述步骤gh，包括如下步骤：
步骤gh-1.针对目标路段短时交通流预测模型，将其中核函数系数、惩罚因子作为待优化参数输入细菌觅食模型当中；
步骤gh-2.根据优性原则，设置细菌觅食算法的主要参数，细菌个数s，单向运动的最大步数n0趋向性循环操作次数n1，复制性循环操作次数n2，迁徙性循环操作次数n3，迁徙概率p；
步骤gh-3.根据细菌觅食算法，针对核函数系数、惩罚因子进行寻优，优化更新核函数系数、惩罚因子，进而更新目标路段短时交通流预测模型。
本发明所述一种基于时空特性分析的短时交通流预测方法采用以上技术方案与现有技术相比，具有以下技术效果：本发明设计的基于时空特性分析的短时交通流预测方法，以对交通流的内在联系与特性进行分析，找寻数据的内在规律，改善模型输入数据的质量；并运用得当方法对预测模型的重要参数进行优化，改善预测模型效果，从而提高预测模型的可行性和有效性，提高短时交通流预测精确度，实现智能交通系统实时准确的交通流预测。
附图说明
图1是本发明所设计基于时空特性分析的短时交通流预测方法的流程示意图；
图2是本发明所设计基于时空特性分析的短时交通流预测方法中聚类方法流程示意图；
图3是本发明所设计基于时空特性分析的短时交通流预测方法中参数优化方法流程示意图。
具体实施方式
下面结合说明书附图对本发明的具体实施方式作进一步详细的说明。
时空相关性分析：
(1)、周期性：由于人们的生活习惯，可以通过数据观察出，同一路段一周内的交通流数据呈现周期性变化，周期近似为一天。
(2)、时延性：不同路段间交通流数据存在相关性，即时延性，通过前一路段的车辆有部分同样也通过下一路段，然而这种交通流变化的反应需要一段时间可体现在数据当中。
(3)、混沌性：由于路况天气的不确定和多变性，导致交通流数据可能出现短期内的不稳定的非线性变化以及异常值的诞生。
如图1所示，本发明设计了一种基于时空特性分析的短时交通流预测方法，用于实现目标道路的交通流量预测，包括如下步骤：
步骤a.获得目标道路指定历史时间区间内的实时交通流量，然后进入步骤b。
步骤b.采用历史平均法、相邻数据平均法或加权平均法中的一种或几种任意组合，针对目标道路指定历史时间区间内实时交通流量，进行预处理操作，包括丢失值处理、错误值处理和异常值处理，获得目标道路指定历史时间区间内各个时间节点对应的交通流量，然后进入步骤c。
对于上述历史平均法、相邻数据平均法或加权平均法，具体如下：
(1)、历史平均法：将不同日期相同时刻的数据的均值作为缺失值的估计值(如公式1-1)。
其中，y(t)为缺失值的修正值，yi(t)为第i天、第t时刻的值，n为天数。
(2)、相邻数据平均法：将相邻时刻的均值作为估计值(如公式1-2)。
y(t)＝(y(t-1)+y(t+1))/2(1-2)
(3)、加权平均法：利用相邻历史数据和相邻实际测量值的加权平均值作为估计值如公式(1-3)。
y(t)＝γy(t-1)+(1-γ)yi(t)(1-3)
其中，γ为加权系数。
步骤c.将目标道路所对应指定历史时间区间内各个时间节点的交通流量，作为各个数据样本点，并针对各个数据样本点，进行维度空间坐标转换，获得各个数据样本点分别所对应的维度空间坐标，然后进入步骤d。
步骤d.基于各个数据样本点分别所对应的维度空间坐标，针对各个数据样本点进行optics聚类处理，获得各个聚类簇，然后进入步骤e。
如图2所示，上述步骤d中，初始化用于存放本次循环中核心点、及其直接密度可达点的集合ⅰ，以及包含各核心点分别所对应簇的集合ⅱ，然后基于各个数据样本点分别所对应的维度空间坐标，执行如下步骤,针对各个数据样本点进行聚类处理，获得各个聚类簇；
步骤d1.判断是否存在非核心点、以及未被划簇的数据样本点，是则进入步骤d2；否则进入步骤d10；
步骤d2.针对非核心点、以及未被划簇的各个数据样本点，任选一个数据样本点，作为本次循环的核心点，并进入步骤d3；
步骤d3.将该核心点划分至集合ⅱ当中，并获得该核心点所对应的所有密度可达的数据样本点，然后进入步骤d4；
步骤d4.针对该核心点所对应的所有密度可达的数据样本点，判断其中直接密度可达的数据样本点是否在集合ⅱ当中，是则进入步骤d5；否则步骤d6；
步骤d5.分别针对该核心点所对应直接密度可达的各个数据样本点，判断数据样本点对应的可达距离是否小于预设半径，是则将该数据样本点划分至该核心点所对应的簇中，并更新集合ⅰ中该核心点所对应的簇，且返回步骤d1；否则将数据样本点直接删除，并返回步骤d1；
步骤d6.针对该核心点所对应的所有密度可达的数据样本点，按可达距离升序，针对各个数据样本点进行排序，划入该核心点所对应的簇，并加入集合ⅰ，然后进入步骤d7；
步骤d7.由该核心点所对应簇中，剪切出第一个数据样本点，作为本次循环的待判断数据样本点，并进入步骤d8；
步骤d8.判断该待判断数据样本点是否属于集合ⅱ，是则将该待判断数据样本点，作为核心点，并返回步骤d3；否则进入步骤d9；
步骤d9.判断集合ⅰ是否为空，是则返回步骤d1；否则返回步骤d7；
步骤d10.即完成针对各个数据样本点的聚类处理，获得各个聚类簇。
步骤e.根据各个聚类簇，获得目标路段分别对应各聚类时间节点的交通流量，并针对各个交通流量进行数据归一化处理，提高模型精准度和降低计算复杂度，即更新目标路段分别对应各聚类时间节点的交通流量，然后进入步骤f。
步骤f.基于目标路段分别对应各聚类时间节点的交通流量，结合预设预测时长跨度，获得各组训练样本(xk,yk)，然后进入步骤g；其中，k∈{1,…,k}，k表示训练样本的组数，xk、yk均属于目标路段分别对应各聚类时间节点的交通流量，xk所对应聚类时间节点与yk所对应聚类时间节点之间的差值等于预测时长跨度，且按时序、xk所对应聚类时间节点在yk所对应聚类时间节点之前。
步骤g.采用各组训练样本(xk,yk)，基于xk为输入，yk为输出，针对svm分类器进行训练，获得对应预测时长跨度的目标路段短时交通流预测模型，然后进入步骤gh。
上述步骤g具体如下：
(1)、将聚类簇作为svm训练子样本集作为输入，即采用各组训练样本(xk,yk)，基于xk为输入，yk为输出。
(2)、本模型选择核函数，要求具有较强的学习性能，即误差较小的径向基(rbf)核函数，如公式(3-1)：
(3)、通过剖析已知数据，求出核函数系数ε、惩罚因子c等重要参数。
惩罚因子c的选择公式如下：
由于k(xi,x)≤1，可得|f(x)|≤cn
其中
(4)、构建核函数优化问题求解，如公式(3-4)：
(5)、决策函数由求得的最优解建立，预测值即公式(3-5)
步骤gh.采用细菌觅食算法(bfoa)，针对目标路段短时交通流预测模型中的参数进行优化，然后进入步骤h。如此通过细菌间的趋向性操作，复制操作和迁徙操作，利用群体间的竞争与协作实现核心参数寻优，获得最优解，优化预测模型性能。
如图3所示，上述步骤gh，包括如下步骤：
步骤gh-1.针对目标路段短时交通流预测模型，将其中核函数系数、惩罚因子作为待优化参数输入细菌觅食模型当中；
步骤gh-2.根据优性原则，设置细菌觅食算法(bfoa)的主要参数，细菌个数s，单向运动的最大步数n0趋向性循环操作次数n1，复制性循环操作次数n2，迁徙性循环操作次数n3，迁徙概率p；
步骤gh-3.根据细菌觅食算法(bfoa)，针对核函数系数、惩罚因子进行寻优，优化更新核函数系数、惩罚因子，进而更新目标路段短时交通流预测模型。
具体来说，采用细菌觅食算法(bfoa)进行参数优化，具体如下：
(1)确定影响bfoa寻优能力和优化速度的参数：细菌个数s，单向运动的最大步数n0趋向性循环操作次数n1，复制性循环操作次数n2，迁徙性循环操作次数n3，迁徙概率p。
(2)趋向性操作
即细菌向有利于自身环境的区域移动，包括翻转运动(任意方向单位步长)和游动运动(延上一步移动方向单位步长)。第i个细菌进行翻转运动时根据式(1)更新位置。
其中，pi(l,m,n)表示第i个细菌的第1次趋向操作，第m次复制操作，以及第n次迁徙操作后的位置；c(i)表示按选定方向游动的单位步长向量；
如果在pi(l+1,m,n)处的适应度优于pi(l,m,n)处，则保持
(3)复制操作
细菌达到设定的趋向性循环次数时，细菌进行自我复制，产生新的个体。自我复制前，对细菌种群的适应度降序排列，较大的前一半保留进行分裂，后一半淘汰。保证种群大小不变，复制后保留母细菌特性。
(4)迁徙操作
细菌完成固定周期内的复制操作后，以概率p进行迁徙操作。如果细菌满足迁徙发生概率则该个体死亡，同时解空间任意位置重新生成一个新个体，相当于把原来个体移到一个新的位置，保持其他细菌位置不变。迁徙操作随机生成的细菌个体更接近全局最优解，从而跳出局部最优解的可能性。
步骤h.基于所获目标路段短时交通流预测模型，针对目标路段，实现预测时长跨度的交通流量预测。
上面结合附图对本发明的实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。