CN113254866A - 一种交通短时流量智能预测方法 - Google Patents
一种交通短时流量智能预测方法 Download PDFInfo
- Publication number
- CN113254866A CN113254866A CN202110555809.0A CN202110555809A CN113254866A CN 113254866 A CN113254866 A CN 113254866A CN 202110555809 A CN202110555809 A CN 202110555809A CN 113254866 A CN113254866 A CN 113254866A
- Authority
- CN
- China
- Prior art keywords
- model
- formula
- following formula
- subsequences
- kernel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
- G06F17/156—Correlation function computation including computation of convolution operations using a domain transform, e.g. Fourier transform, polynomial transform, number theoretic transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明提供一种交通短时流量智能预测方法,包括以下步骤,利用第一模型分别提取原始交通流序列中的非线性特征和高频噪声,将原始交通流序列分解为几个中心频率不同的子序列;分别建立第二模型和第三模型,并对不同中心频率的子序列进行预测,采用第二模型模型预测低中心频率子序列,采用第三模型预测高中心频率子序列;将各子序列的预测值相加。本发明具有较高的精度和稳定性。
Description
技术领域
本发明涉及交通领域,具体涉及一种交通短时流量智能预测方法。
背景技术
随着ITS的发展,短期交通流预测一直是研究热点,其预测精度对ITS性能的提高起着关键作用。因此,发展高精度的预测模型是必要的。目前,各种传统的预测模型已经被建立起来,其实验结果证明是有效的,但仍存在缺陷需要改进。因此,迫切需要探索更可靠、更精确的短期交通流预测方法。
发明内容
本发明的发明目的是,提供一种更可靠、更精确的短期交通流预测方法。
本发明提供一种交通短时流量智能预测方法,包括以下步骤,
利用第一模型分别提取原始交通流序列中的非线性特征和高频噪声,将原始交通流序列分解为几个中心频率不同的子序列;
分别建立第二模型和第三模型,并对不同中心频率的子序列进行预测,采用第二模型模型预测低中心频率子序列,采用第三模型预测高中心频率子序列;
将各子序列的预测值相加。
进一步的,
所述第一模型包括,
式中,s(t)表示信号中包含高频噪声;s0(t)表示过滤后的信号;||·||2表示函数的第二范数,f(t)表示待求信号,t表示时间序列;
将上述公式通过傅立叶变换方法得到以下公式:
式中,s0(ω)、s(ω)分别表示s0(t)和s(t)的傅立叶变换;α是过滤因子;
根据下述公式将s(t)分为k个子系列,
式中,si(t)表示分解的子系列,k表示子系列的个数。
进一步的,
所述第二模型,在通过变分模态分解后,为低频子系列建立第二模型;
设有样本空间{(xi,yi),i=1,2,...n},xi∈Rd,yi∈R,R表示样本空间,d表示自变量维数,n表示样本的个数。xi,yi表示样本,存在一个映射f:Rd→R,以最小化回归值f(xi)和对象值yi之间的均方根误差,
然后引入第二模型:
式中,v表示惩罚因子,μ表示灵敏度,e表示预测值和对象值之间的误差,ei表示第i个样本预测值和对象值之间的误差,T表示矩阵转置;
在上述公式中应用拉格朗日乘数法,重写为以下公式。
式中,Ω∈Rn×n由ψ(xi,xj)组成,ψ(xi,xj)是xi和xj、V=[1,1,...1]T∈Rn;γ=ν/μ;υ=[y1,y2,...yn]∈Rn的核函数,高斯核函数因如下述公式所示;I∈Rn×n为单位矩阵;
K(xi,xj)=exp[-||xi-xj||2/(2σ2)],σ>0
式中,σ代表内核参数。然后,通过在样本空间中进行训练,生成下述公式,
式中,xi是样本空间的每个向量。
进一步的,
所述第二模型模型包括,
设获得的时间序列样本为{x1,x2,...xn},将时间序列样本重构为下述公式,
式中,sti(i=1,2,...d)和std+1分别代表{xi,xi+1,...xi+n-d-1}(i=1,2,...d)和{xd+1,xd+2,...xn}的标准差;
核密度函数是高斯核密度函数采用以下公式所示,
将联合分布核函数公式改写为下述公式
本发明的有益效果是,本发明不仅可以有效地将短期交通流划分为子序列,而且还可以根据各自的中心频率,采用两种合适的预测方法对分解后的子序列进行处理。更重要的是,它补充了几种传统模式的优缺点。因此,本发明非常适合于短期交通流预测。
附图说明
图1为本发明流程图。
具体实施方式
短时流量预测是智能交通***(ITS)的重要组成部分。然而,由于交通流数据存在很强的非线性和随机性,一些传统的预测方法难以提供令人满意的结果。为此,本发明提出了一种基于变分模态分解、最小二乘支持向量和条件核密度估计的混合模型(VMD-LSSVM-CKDE)。其中,第一模型首先将隐含的非线性和高频噪声分解为若干个子序列来处理。然后,根据不同的中心频率建立各子序列的第二模型和第三模型,并将相应的预测结果叠加得到最终的预测结果。最后,以一组主干道交叉口实测数据为例,对该方法的性能进行了评价。实验结果表明,该模型的性能优于其他模型。例如,与传统模型相比,该方法在绝对误差、平均相对百分比误差、均方根误差、均方根相对误差和等系数等指标上的改进分别为0.9920、2.25%、1.3457、4.30%和0.9917。
下面多本发明中的缩略语进行说明
AI,人工智能;AMIRA,自回归综合移动平均模型;KF,卡尔曼滤波器;SVM,支持向量;SARIMA,季节自回归综合移动平均模型;ANN,人工神经网络;CNN,卷积神经网络;LSTM,长短期记忆;ELM,极限学***均绝对误差;MRPE,平均相对百分比误差;RMSR,均方根误差;RMSRE,均方根相对误差;EC,等系数。
本发明的发明构思是,首先利用第一模型分别提取原始交通流序列中的非线性特征和高频噪声,将原始交通流序列分解为几个中心频率不同的子序列;其次,分别建立第二模型和第三模型,并对不同中心频率的子序列进行预测,第二模型预测低中心频率子序列,第三模型预测高中心频率子序列。第三,将各子序列的预测值相加,得到最终预测值。最后,基于一组交通量数据的案例研究,对其性能进行了评价。实验结果表明,该模型优于常规模型。与传统预测模型相比,本发明的有益效果是:(1)第一模型成功地应用于短期交通流预测领域。同时,得到了满意的结果。就我们所知,新发展的分解模型可以解决传统方法的问题,如模式混合。(2)根据子序列的中心频率不同,分别提出了不同的预测方法。特别的是,与常见的单一预测模型不同,该模型能够适应子序列的中心频率,并选择合适的模型(LSSVM或CKDE)对分解后的子序列进行预测。(3)利用一组城市道路实测数据集作为实验数据,对所提模型的性能进行了全面的研究。为了更好地展示其性能,我们使用了四个参考模型进行比较。结果表明,该方法具有较高的精度和稳定性。
短期交通流数据具有低频,非线性和随机性等特点,因此很难通过单一模型进行有效预测,本发明建立了一个新的组合预测模型。
下面对第一模型进行说明。
将短期交通流数据分解为具有不同中心频率的几个子系列后,第一模型可以区分非线性和随机性的组成部分。
作为一种新的分解方法,第一模型能够自适应地消除高频噪声,表达式如公式(1)所示:
式中,s(t)表示信号中包含高频噪声;s0(t)表示过滤后的信号;||·||2表示函数的第二范数,f(t)表示待求信号,t表示时间序列;
将公式(1)通过傅立叶变换方法得到以下公式(2):
式中,s0(ω)、s(ω)分别表示s0(t)和s(t)的傅立叶变换;α是过滤因子。显然,可以通过在降低真实信号损失的条件下转动滤波器系数来滤除高频噪声。
可以根据公式(3)根据下述公式将s(t)分为k个子系列,
式中,si(t)表示分解的子系列,k表示子系列的个数。
在继续下一个步骤之前,需要对子系列执行几个过程:
1)为了获得子序列的单边频谱,使用希尔伯特变换来处理子序列。
2)为了避免子系列之间的频域重叠,每个子系列都应与指数混合,以便将其中心频率转换为新的基带;
3)每个子序列的带宽可以通过解调信号的高斯平滑度来估计。
因此分解问题改写为(4):
式中,si(t)和ωi表示分解后的子序列和相应的中心频率;δ(t)是狄拉克函数;j是虚单位;*表示卷积运算。
为了消除约束项,可以通过添加惩罚项和Lagrange项将公式(4)转换为(5)。
式中,<·>表示内部产品;λ(t)表示拉格朗日乘数;α表示过滤因子,i为范围在1到n的自然数。
通过变分原理和ADMM([20]~[22])解决了优化问题(5)。ADMM也称为乘法器的交替方向方法。以这种方式,通过逐步迭代来获得的解。详细的迭代过程如公式(6)、(7)、(8)所示。
当满足公式(9)时,整个迭代过程将停止。
式中,ξ表示收敛系数,它决定了原始交通流与重建交通流之间的误差。
下面对第二模型进行说明。
在通过第一模型分解后,为低频子系列建立第二模型。它通常用于解决回归问题,并在描述交通流数据中存在的强非线性方面表现出色。
假设有一个样本空间{(xi,yi),i=1,2,...n},xi∈Rd,yi∈R,R表示样本空间,n表示样本的个数,xi,yi表示样本,d为样本自变量的维数,存在一个映射f:Rd→R,以最小化回归值f(xi)和对象值yi之间的均方根误差。
然后在公式中引入第二模型(12)
式中,v表示惩罚因子,μ表示灵敏度,e表示预测值和对象值之间的误差,ei表示第i个样本预测值和对象值之间的误差,T表示矩阵转置。
在公式(12)中应用了拉格朗日乘数法,因此公式(12)可以如下重写。
式中,L(·)表示拉格朗日函数;βi表示拉格朗日乘数。和b的解由KKT条件和矩阵公式(14)获得,KKT条件是非线性规划(nonlinear programming)最佳解的必要条件。KKT条件将Lagrange乘数法(Lagrange multipliers)所处理涉及等式的约束优化问题推广至不等式。
式中,Ω∈Rn×n由ψ(xi,xj)组成,ψ(xi,xj)是xi和xj、V=[1,1,...1]T∈Rn;γ=ν/μ;υ=[y1,y2,...yn]∈Rn的核函数。通常,高斯核函数因简单实用被广泛运用,如公式(15)所示。
K(xi,xj)=exp[-||xi-xj||2/(2σ2)],σ>0 (15)
式中,σ代表内核参数。然后,通过在样本空间中进行训练,将公式(11)改写为公式(16)。
式中,xi是样本空间的每个向量。
下面对第三模型进行说明。
第三模型是一种相对于核密度估计的改进方法,它针对更高频率的子序列而开发,有关核密度估计。假设获得的时间序列样本为{x1,x2,...xn}。考虑到时间序列之间存在复杂的关系。因此可以将时间序列样本重构为公式(17)。
式中,sti(i=1,2,...d)和std+1分别代表{xi,xi+1,...xi+n-d-1}(i=1,2,...d)和{xd+1,xd+2,...xn}的标准差。
另一方面,最广泛使用的核密度函数是高斯核密度函数,如公式(21)所示。
结合公式(19)(21),可以将公式(18)改写为公式(22)
本发明进行短期交通流预测的主要步骤:1)选择合适的过滤因子对原始交通流数据进行过滤;2)将过滤后的交通流数据分解为若干子序列;3)采用第二模型或第三模型对不同中心频率的子序列进行预测;4)将几个子序列的预测值相加,得到最终的预测值。
根据公式(2),原始交通流数据的损失与过滤因子的变化呈正相关,其中过滤因子的选择不要太大。本发明采用的滤波因子的取值范围在0到100之间,同时满足以下准则:1)残差必须服从高斯分布;2)残差分布与高斯分布的拟合误差均方误差最小;
通过第一模型分解得到不同中心频率的子级数,中心频率越高的子级数振荡越剧烈,反之亦然。因此,利用第二模型自身特有的拟合特征对低中心频率进行预测。另一方面,作为一种非参数估计方法,第三模型是一种估计低振幅、强振荡时间序列的有效方法。因此利用第三模型预测高中心频率子序列。众所周知,第二模型是在假定训练部分残差必须服从均值为零的高斯分布的情况下建立的,这可以作为预测结果是否有效的判断标准。也就是说,随着中心频率的增加,第二模型对子序列训练部分的残差不太可能服从高斯分布,就这样选择临界中心频率。具体选择过程如表1所示。
表1 预测方法的选择
本发明所提出的预测模型的具体过程如图1所示。
本发明针对传统预测方法的不足,提出了一种新的预测方法显著提高了短期交通流预测的准确性,具有优越性和可靠性。
本发明具有更高的预测精度,在短期交通流中的应用是有效的。其次,从与单一模型比较的角度来看:与ARIMA模型相比,本发明模型具有解释时间序列非线性分量的能力;
综上所述,本发明不仅可以有效地将短期交通流划分为子序列,而且还可以根据各自的中心频率,采用两种合适的预测方法对分解后的子序列进行处理。更重要的是,它补充了几种传统模式的优缺点。因此,本发明非常适合于短期交通流预测。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种交通短时流量智能预测方法,其特征在于,包括以下步骤,
利用第一模型分别提取原始交通流序列中的非线性特征和高频噪声,将原始交通流序列分解为几个中心频率不同的子序列;
分别建立第二模型和第三模型,并对不同中心频率的子序列进行预测,采用第二模型模型预测低中心频率子序列,采用第三模型预测高中心频率子序列;
将各子序列的预测值相加。
3.如权利要求1所述的一种交通短时流量智能预测方法,其特征在于,
所述第二模型,在通过变分模态分解后,为低频子系列建立第二模型;
设有样本空间{(xi,yi),i=1,2,...m},xi∈Rd,yi∈R,R表示样本空间,d表示自变量维数,n表示样本的个数,xi,yi表示样本,存在一个映射f:Rd→R,以最小化回归值f(xi)和对象值yi之间的均方根误差,
然后引入第二模型:
式中,v表示惩罚因子,μ表示灵敏度,e表示预测值和对象值之间的误差,ei表示第i个样本预测值和对象值之间的误差,T表示矩阵转置;
在上述公式中应用拉格朗日乘数法,重写为以下公式,
式中,L(·)表示拉格朗日函数,表示拉格朗日乘数,和b的解由KKT条件和下述公式获得
式中,Ω∈Rm×m由ψ(xi,xj)组成,ψ(xi,xj)是xi和xj、V=[1,1,...1]T∈Rm;γ=ν/μ;I∈Rn ×n为单位矩阵;υ=[y1,y2,...ym]∈Rm的核函数,高斯核函数如下述公式所示,
K(xi,xj)=exp[-||xi-xj||2/(2σ2)],σ>0
式中,σ代表内核参数,然后,通过在样本空间中进行训练,生成下述公式,
式中,xi是样本空间的每个向量。
4.如权利要求1所述的一种交通短时流量智能预测方法,其特征在于,所述第二模型包括,
设获得的时间序列样本为{x1,x2,...xn},将时间序列样本重构为下述公式,
式中,sti(i=1,2,...d)和std+1分别代表{xi,xi+1,...xi+n-d-1}(i=1,2,...d)和{xd+1,xd+2,...xn}的标准差;
核密度函数是高斯核密度函数采用以下公式所示,
将联合分布核函数公式改写为下述公式
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110555809.0A CN113254866B (zh) | 2021-05-21 | 2021-05-21 | 一种交通短时流量智能预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110555809.0A CN113254866B (zh) | 2021-05-21 | 2021-05-21 | 一种交通短时流量智能预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113254866A true CN113254866A (zh) | 2021-08-13 |
CN113254866B CN113254866B (zh) | 2022-02-01 |
Family
ID=77183521
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110555809.0A Active CN113254866B (zh) | 2021-05-21 | 2021-05-21 | 一种交通短时流量智能预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113254866B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103730006A (zh) * | 2014-01-26 | 2014-04-16 | 吉林大学 | 一种短时交通流量的组合预测方法 |
US20180137412A1 (en) * | 2016-11-16 | 2018-05-17 | Cisco Technology, Inc. | Network traffic prediction using long short term memory neural networks |
CN112801386A (zh) * | 2021-02-04 | 2021-05-14 | 北京工业大学 | 一种基于三阶段模型的短时交通流预测方法 |
-
2021
- 2021-05-21 CN CN202110555809.0A patent/CN113254866B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103730006A (zh) * | 2014-01-26 | 2014-04-16 | 吉林大学 | 一种短时交通流量的组合预测方法 |
US20180137412A1 (en) * | 2016-11-16 | 2018-05-17 | Cisco Technology, Inc. | Network traffic prediction using long short term memory neural networks |
CN112801386A (zh) * | 2021-02-04 | 2021-05-14 | 北京工业大学 | 一种基于三阶段模型的短时交通流预测方法 |
Non-Patent Citations (1)
Title |
---|
李淑庆、王延鹏、俞艇、赵磊娜: "基于韦布尔分布的城市快速路路段通行能力可靠性分布", 《重庆交通大学学报(自然科学版)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113254866B (zh) | 2022-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Madan et al. | Predicting computer network traffic: a time series forecasting approach using DWT, ARIMA and RNN | |
Li et al. | Short-term load-forecasting method based on wavelet decomposition with second-order gray neural network model combined with ADF test | |
Ratto et al. | State dependent parameter metamodelling and sensitivity analysis | |
CN110334875B (zh) | 计及评估指标冲突的风电功率组合概率预测方法 | |
CN109886464B (zh) | 基于优化奇异值分解生成特征集的低信息损失短期风速预测方法 | |
CN108647807B (zh) | 河流流量的预测方法 | |
CN109359404B (zh) | 基于经验小波消噪和神经网络融合的中长期径流预报方法 | |
KR20160019897A (ko) | 시계열의 고속 그룹화 기법 | |
CN111950759A (zh) | 一种基于两级分解、lstm和at的短期风速预测方法 | |
WO2024037393A1 (zh) | 拓扑识别模型的压缩方法、装置、电子设备及介质 | |
Ji et al. | Strategies for multi-step-ahead available parking spaces forecasting based on wavelet transform | |
CN111415008B (zh) | 一种基于vmd-foa-grnn的船舶流量预测方法 | |
CN111553513A (zh) | 一种基于二次分解与回声状态网络的中长期径流预测方法 | |
Chou et al. | Application of wavelet‐based multi‐model Kalman filters to real‐time flood forecasting | |
Chou | Applying multi-resolution analysis to differential hydrological grey models with dual series | |
CN111598313A (zh) | 一种风功率预测方法 | |
Lian | Runoff forecasting model based on CEEMD and combination model: a case study in the Manasi River, China | |
Garai et al. | An MRA Based MLR Model for Forecasting Indian Annual Rainfall Using Large Scale Climate Indices | |
CN112381279B (zh) | 一种基于vmd和bls组合模型的风电功率预测方法 | |
CN113254866B (zh) | 一种交通短时流量智能预测方法 | |
Shabri | A hybrid model for stream flow forecasting using wavelet and least Squares support vector machines | |
CN117155806A (zh) | 一种通信基站流量预测方法与装置 | |
Cai et al. | Phasednn-a parallel phase shift deep neural network for adaptive wideband learning | |
CN116192665B (zh) | 数据处理方法、装置、计算机设备及存储介质 | |
CN116822742A (zh) | 一种基于动态分解-重构集成处理的电力负荷预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |