买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：江苏大学

摘要：本发明涉及一种基于集成LSSVR模型的路网中短时交通流的预测方法。该预测方法在传统方法的基础上，根据不同路段的交通流量对不同滞后时间的历史交通流数据的依赖程度的不同，先是利用不同的滞后时间，搭建多个不同的非线性LSSVR模型，学习到不同路段的交通流对不同滞后时间的历史数据的依赖，然后通过一个线性LSSVR模型，对多个非线性LSSVR的输出进行整合,进行最终预测。同时，针对模型中参数的调节问题，利用和声搜索算法予以解决。与传统的预测方法相比，本发明可以达到更好的预测效果。

主权项：1.一种基于集成LSSVR模型的路网中短时交通流的预测方法，其特征在于，包括如下步骤：1采集路网中所有监测器的交通流数据，其中，路网包含多个站点，站点集合为S1，每一个站点包括多个监测器；2通过监测器的交通流数据在站点集合S1中筛选站点，组成集合S2，并计算S2的交通流数据；3构建路网短时交通流量值矩阵D2；4构造LSSVR模型的训练集和验证集；5构造单个LSSVR模型；6构建集成LSSVR模型；7利用智能优化算法对集成LSSVR模型参数进行优化；8将待预测的数据集，按步骤2到步骤4处理成符合集成LSSVR的模型的输入形式，输入到模型中进行预测；所述步骤6中，构建集成LSSVR模型方法如下：6.1定义一个计数器t5，令t5＝1；6.2依照步骤4和步骤5，构造之后时间为t的非线性LSSVR模型，及其对应的训练集和测试集，其预测输出记为6.3如果t5＝P，则t5＝t5+1并转至步骤6.4，否则转至步骤6.2；6.4构建线性LSSVR模型的输入D″T，构建方法如下：D″T＝[f1x...fPx]6.5依照步骤5构建线性LSSVR模型，其训练集和测试集的输入矩阵按照步骤6.4构造，输入矩阵与非线性LSSVR模型相同，其输出即为集成LSSVR模型的预测。

全文数据：基于集成LSSVR模型的路网中短时交通流的预测方法技术领域本发明涉及智能交通技术领域，具体为一种基于集成LSSVR模型的路网中短时交通流的预测方法。背景技术交通流量预测作为智能交通系统ITS的基础，较为准确的短时预测是实现科学规范的交通规划、交通诱导和交通控制的关键技术。是实现交通控制与诱导的前提，也是使ITS从“被动式反应”转变到“主动式动作”的关键。一些经典的模型包括历史平均模型，卡尔曼状态空间滤波模型，自回归移动平均模型ARIMA季节性自回归移动平均模型SARIMA，K近邻等已被用于短时交通流的预测。这些方法主要关注了交通流在特定位置的时间相关性并在相对稳定的交通流数据上表现良好，此外，这些模型假定交通流是线性的，但在真实的交通中，不同位置的交通流在时空上相互影响，是非线性的。因此，这些方法在短时交通流预测上不能取得理想的效果。人工神经网络ANN是另一种流行的数据驱动的交通流预测方法。其中，得到广泛认可的出色ANN模型有多层感知机MLP，卷积神经网络CNN，长短时记忆LSTM网络，自编码器AE等。近期，多隐层的ANN被用于交通流预测，并取得了较一些传统方法更好的结果。尽管ANN在构建非线性模型上是如此的强大，它仍然存在一些问题。例如，其学习过程缺乏可解释性并且在网络模型很复杂的情况下，训练ANN将需要大量的数据和算力。支持向量机SVM，是一种有监督的机器学习模型，已被成功的应用在许多实际问题中。通过引入不敏感损失函数ε，SVM可以转化为支持向量回归SVR，解决非线性回归问题。在此基础上，最小二乘支持向量回归LSSVR通过最小二乘损失函数替换掉不敏感损失函数ε的方式，简化了SVR。在这种方式下，LSSVR的解可以通过求解线性方程组而不是复杂的二次规划问题得到，并同时不牺牲泛化性能。SVR和LSSVR都已被用于交通流预测，并得到了可观的成果。但是，已有的基于SVR的预测模型多是基于某个确定的滞后时间下的单一模型，这种方式没有充分考虑到路网交通中不同路段对滞后时间的依赖性不同的问题，其一侧性能仍有提升的空间。发明内容针对上述存在的问题，本发明提供了一种基于集成LSSVR模型的路网中短时交通流的预测方法，以提高预测效果。为了实现上述目的，本发明具体技术方案如下：一种基于集成LSSVR模型的路网中短时交通流的预测方法，包括如下步骤：1采集路网中所有监测器的交通流数据，其中，路网包含多个站点，站点集合为S1，每一个站点包括多个监测器；2通过监测器的交通流数据在站点集合S1中筛选站点，组成集合S2，并计算S2的交通流数据；3构建路网短时交通流量值矩阵D2；4构造LSSVR模型的训练集和验证集；5构造单个LSSVR模型；6构建集成LSSVR模型；7利用智能优化算法对集成LSSVR模型参数进行优化；8将待预测的数据集，按步骤2到步骤4处理成符合集成LSSVR的模型的输入形式，输入到模型中进行预测。进一步地，上述步骤2包括如下步骤：2.1统计各监测器数据的缺失率，删除缺失率大于1％的监测器所对应的站点，剩下的所有站点组成集合S2；2.2利用相邻数据平均法对S2中每个站点的每个监测器的缺失值进行补全；2.3对集合S2中的所有站点，将同一站点内的不同监测器同一时段上的交通流数据相加，得到S2的交通流数据，站点S2的交通流数据为N阶数列。进一步地，上述步骤3中，路网短时交通流量值矩阵的构建方法如下：3.1将S2中每个站点的交通流数据作为一个向量，将所有站点所对应的向量叠加得到M*N阶的样本数据矩阵D1，其中M为S2中的站点个数，N为每个站点采集的交通流数据的个数；3.2归一化处理样本数据矩阵D1，得到路网短时交通流量值矩阵D2，记为D2的计算公式如下：式中D2为路网短时交通流量值矩阵，xmax为D1中所有元素的最大值。进一步地，上述步骤4中，构造LSSVR模型的训练集和验证集包括如下步骤：4.1构建滞后时间为T的LSSVR模型，其中，预测的目标站点为Mt，集成LSSVR模型中的非线性LSSVR模型的个数为P，非线性LSSVR模型对应的滞后时间T的取值为1至P；4.2定义计数器变量t，令t＝1；4.3利用D2中的元素，构造一个M*T阶的矩阵dt，构造方法如下：4.4将dt转化为一个1*M*T的矩阵d′t，转换方法如下：4.5计算对应的期望输出yt，计算方法如下：4.6判断t是否小于N-T，若是则t＝t+1转步骤4.3，否则转步骤4.7；4.7构造样本集的输入矩阵D′T，构建方法如下：4.8取Ntext作为测试集的样本个数，训练集的样本个数Ntrain为：Ntrain＝N-P-Ntest4.9从样本集的输入矩阵D′T中分割出训练集的输入矩阵D′Ttrain和测试集的输入矩阵D′Ttest，方法如下：4.10构造样本集的期望输出矩阵YT，构造方法如下：4.11从样本集的期望输出矩阵YT中分割出训练集的期望输出矩阵YTtrain和测试集的期望输出矩阵YTtest，方法如下：进一步地，上述步骤5中，构造单个LSSVR模型的方法如下：5.1计算LSSVR模型的α和b的值，计算方法如下：α＝H-1y-H-1Ob其中K是一个M×M的核矩阵，其元素为Kij＝Kxi，xj＝xiTxj，O是一个全1向量，I是一个具有合适维数的单位矩阵，y是训练集的期望输出Ntrain*1，γ是平衡参数；5.2利用核函数求解LSSVR模型的核矩阵K，其中，线性LSSVR模型的核矩阵K的计算方法如下：Kxi，xj＝xiTxj式中xi为测试集的输入矩阵的第i行的向量，xj为训练集的输入矩阵的第j行的向量，Kxi，xj为一个数值，是核矩阵K第i行j列的元素；非线性LSSVR模型的核矩阵K的计算方法如下：式中σ为核宽度参数；5.3计算LSSVR模型的预测结果fx，计算方法如下：进一步地，上述步骤6中，构建集成LSSVR模型方法如下：6.1定义一个计数器t5，令t5＝1；6.2依照步骤4和步骤5，构造之后时间为t的非线性LSSVR模型，及其对应的训练集和测试集。其预测输出记为6.3如果t5＝P，则t5＝t5+1并转至步骤6.4，否则转至步骤6.2；6.4构建线性LSSVR模型的输入D″T，构建方法如下：D″T＝[f1x…fPx]6.5依照步骤5构建线性LSSVR模型，其训练集和测试集的输入矩阵按照步骤6.4构造，输入矩阵与非线性LSSVR模型相同，其输出即为集成LSSVR模型的预测。与传统的预测方法相比，本发明提高了预测精度和预测速度。附图说明图1是短时交通流的预测方法流程图。图2是本发明的集成LSSVR模型图。具体实施方式下面结合附图对本发明作进一步描述，一种基于集成LSSVR模型的路网中短时交通流的预测方法流程如图1所示，包括如下步骤：1采集路网中所有监测器的交通流数据，监测器监测到的交通流数据是由监测器每隔时间Δt采集的交通流观测值组成的N阶数列；其中，路网包含多个站点，站点集合为S1，每一个站点包括多个监测器；2筛选站点并计算站点的交通流数据，包括如下步骤：2.1统计各监测器数据的缺失率，删除缺失率大于1％的监测器所对应的站点，剩下的所有站点组成集合S2；2.2利用相邻数据平均法对S2中每个站点的每个监测器的缺失值进行补全；2.3对S2中的所有站点，将同一站点内的不同监测器同一时段上的交通流数据相加，得到站点的交通流数据，站点的交通流数据为N阶数列；3构建路网短时交通流量值矩阵，包括如下步骤：3.1将S2中每个站点的交通流数据作为一个向量，将所有站点所对应的向量叠加得到M*N阶的样本数据矩阵D1，其中M为S2中的站点个数，N为每个站点采集的交通流数据的个数；本具体实施例中，M＝24,采集的天数为10,Δt为15min，那么每个站点一天的交通流数据个数为6015*24＝96，则N＝96*10＝960；3.2归一化处理样本数据矩阵D1，得到路网短时交通流量值矩阵D2，记为D2的计算公式如下：式中D2为路网短时交通流量值矩阵，xmax为D1中所有元素的最大值；4构造LSSVR模型的训练集和验证集，包括如下步骤：4.1构建滞后时间为T的LSSVR模型，其中，预测的目标站点为Mt，集成LSSVR模型中的非线性LSSVR模型的个数为P，非线性LSSVR模型对应的滞后时间T的取值为1至P；4.2定义计数器变量t，令t＝1；4.3利用D2中的元素，构造一个M*T阶的矩阵dt，构造方法如下：4.4将dt转化为一个1*M*T的矩阵d′t即为一个输入样本，转换方法如下：4.5计算对应的期望输出yt，计算方法如下：4.6判断t是否小于N-T，若是则t＝t+1转步骤4.3，否则转步骤4.7；4.7构造样本集的输入矩阵D′T，构建方法如下：4.8取Ntest作为测试集的样本个数，则训练集的样本个数Ntrain为：Ntrain＝N-P-Ntest4.9从样本集的输入矩阵D′T中分割出训练集的输入矩阵D′Ttrain和测试集的输入矩阵D′Ttest，方法如下：4.10构造样本集的期望输出矩阵YT，构造方法如下：4.11从样本集的期望输出矩阵YT中分割出训练集的期望输出矩阵YTtrain和测试集的期望输出矩阵YTtest，方法如下：5构造单个LSSVR模型5.1计算LSSVR模型的α和b的值，计算方法如下α＝H-1y-H-1Ob其中K是一个M×M的核矩阵，其元素为Kij＝Kxi，xj＝xiTxj，O是一个全1向量，I是一个具有合适维数的单位矩阵，y是训练集的期望输出Ntrain*1，γ是平衡参数；5.2利用核函数求解LSSVR模型的核矩阵K，其中，线性LSSVR模型的核矩阵K的计算方法如下：Kxi，xj＝xiTxj式中xi为测试集的输入矩阵的第i行的向量，xj为训练集的输入矩阵的第j行的向量，Kxi，xj为一个数值，是核矩阵K第i行j列的元素；非线性LSSVR模型的核矩阵K的计算方法如下：式中σ为核宽度参数；5.3计算LSSVR模型的预测结果fx，计算方法如下：6构建集成LSSVR模型，集成LSSVR模型包括两层，如图2所示。6.1定义一个计数器t5，令t5＝1；6.2依照步骤4和步骤5，构造之后时间为t的非线性LSSVR模型，及其对应的训练集和测试集。其预测输出记为6.3如果t5＝P，则t5＝t5+1并转至步骤6.4，否则转至步骤6.2；6.4构建线性LSSVR模型的输入D″T，构建方法如下：D″T＝[f1x...fPx]6.5依照步骤5构建线性LSSVR模型，其训练集和测试集的输入矩阵按照步骤6.4构造，输入矩阵与非线性LSSVR模型相同，其输出即为集成LSSVR模型的预测。7利用智能优化算法对集成LSSVR模型参数进行优化，本具体实施例选用的是和声搜索算法；8将待预测的数据集，按步骤2到步骤4处理成符合该集成LSSVR的模型的输入形式后，输入到模型中进行预测。

权利要求：1.一种基于集成LSSVR模型的路网中短时交通流的预测方法，其特征在于,包括如下步骤：1采集路网中所有监测器的交通流数据，其中，路网包含多个站点，站点集合为S1，每一个站点包括多个监测器；2通过监测器的交通流数据在站点集合S1中筛选站点，组成集合S2，并计算S2的交通流数据；3构建路网短时交通流量值矩阵D2；4构造LSSVR模型的训练集和验证集；5构造单个LSSVR模型；6构建集成LSSVR模型；7利用智能优化算法对集成LSSVR模型参数进行优化；8将待预测的数据集，按步骤2到步骤4处理成符合集成LSSVR的模型的输入形式，输入到模型中进行预测。2.如权利要求1所述的基于集成LSSVR模型的路网中短时交通流的预测方法，其特征在于所述步骤2包括如下步骤：2.1统计各监测器数据的缺失率，删除缺失率大于1％的监测器所对应的站点，剩下的所有站点组成集合S2；2.2利用相邻数据平均法对S2中每个站点的每个监测器的缺失值进行补全；2.3对集合S2中的所有站点，将同一站点内的不同监测器同一时段上的交通流数据相加，得到S2的交通流数据，站点S2的交通流数据为N阶数列。3.如权利要求1所述的基于集成LSSVR模型的路网中短时交通流的预测方法，其特征在于所述步骤3中，路网短时交通流量值矩阵的构建方法如下：3.1将S2中每个站点的交通流数据作为一个向量，将所有站点所对应的向量叠加得到M*N阶的样本数据矩阵D1，其中M为S2中的站点个数，N为每个站点采集的交通流数据的个数；3.2归一化处理样本数据矩阵D1，得到路网短时交通流量值矩阵D2，记为D2的计算公式如下：式中D2为路网短时交通流量值矩阵，xmax为D1中所有元素的最大值。4.如权利要求1所述的基于集成LSSVR模型的路网中短时交通流的预测方法，其特征在于所述步骤4中，构造LSSVR模型的训练集和验证集包括如下步骤：4.1构建滞后时间为T的LSSVR模型，其中，预测的目标站点为Mt，集成LSSVR模型中的非线性LSSVR模型的个数为P，非线性LSSVR模型对应的滞后时间T的取值为1至P；4.2定义计数器变量t，令t＝1；4.3利用D2中的元素，构造一个M*T阶的矩阵dt，构造方法如下：4.4将dt转化为一个1*M*T的矩阵d′t，转换方法如下：4.5计算对应的期望输出yt，计算方法如下：4.6判断t是否小于N-T，若是则t＝t+1转步骤4.3，否则转步骤4.7；4.7构造样本集的输入矩阵D′T，构建方法如下：4.8取Ntest作为测试集的样本个数，训练集的样本个数Ntrain为：Ntrain＝N-P-Ntest4.9从样本集的输入矩阵D′T中分割出训练集的输入矩阵D′Ttrain和测试集的输入矩阵D′Ttest，方法如下：4.10构造样本集的期望输出矩阵YT，构造方法如下：4.11从样本集的期望输出矩阵YT中分割出训练集的期望输出矩阵YTtrain和测试集的期望输出矩阵YTtest，方法如下：5.如权利要求1所述的基于集成LSSVR模型的路网中短时交通流的预测方法，其特征在于所述步骤5中，构造单个LSSVR模型的方法如下：5.1计算LSSVR模型的α和b的值，计算方法如下：α＝H-1y-H-1Ob其中K是一个M×M的核矩阵，其元素为Kij＝Kxi，xj＝xiTxj，O是一个全1向量，I是一个具有合适维数的单位矩阵，y是训练集的期望输出Ntrain*1，γ是平衡参数；5.2利用核函数求解LSSVR模型的核矩阵K，其中，线性LSSVR模型的核矩阵K的计算方法如下：Kxi，xj＝xiTxj式中xi为测试集的输入矩阵的第i行的向量，xj为训练集的输入矩阵的第j行的向量，Kxi，xj一个数值，是核矩阵K第i行j列的元素；非线性LSSVR模型的核矩阵K的计算方法如下：式中σ为核宽度参数；5.3计算LSSVR模型的预测结果fx，计算方法如下：6.如权利要求1所述的基于集成LSSVR模型的路网中短时交通流的预测方法，其特征在于所述步骤6中，构建集成LSSVR模型方法如下：6.1定义一个计数器t5，令t5＝1；6.2依照步骤4和步骤5，构造之后时间为t的非线性LSSVR模型，及其对应的训练集和测试集，其预测输出记为6.3如果t5＝P，则t5＝t5+1并转至步骤6.4，否则转至步骤6.2；6.4构建线性LSSVR模型的输入D″T，构建方法如下：D″T＝[f1x...fPx]6.5依照步骤5构建线性LSSVR模型，其训练集和测试集的输入矩阵按照步骤6.4构造，输入矩阵与非线性LSSVR模型相同，其输出即为集成LSSVR模型的预测。

百度查询： 江苏大学 基于集成LSSVR模型的路网中短时交通流的预测方法