CN109740795B - 基于sa-pso的短时交通流预测方法 - Google Patents
基于sa-pso的短时交通流预测方法 Download PDFInfo
- Publication number
- CN109740795B CN109740795B CN201811469397.3A CN201811469397A CN109740795B CN 109740795 B CN109740795 B CN 109740795B CN 201811469397 A CN201811469397 A CN 201811469397A CN 109740795 B CN109740795 B CN 109740795B
- Authority
- CN
- China
- Prior art keywords
- sequence
- model
- equidistant
- traffic flow
- parameter
- 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.)
- Active
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
公开了一种基于SA‑PSO的短时交通流预测方法,属于数据预测领域,该方法基于背景值和初始条件是影响灰色建模精度的两个重要因素,根据车联网的短时交通流数据建立非等间距GM(1,1)模型,利用复化梯形公式对非等间距GM(1,1)模型中的背景值进行了优化,并基于新信息优先原理,对模型的初始条件进行了重构,利用相对误差平方和最小准则,对初始条件对应的最优时间参数进行求解,并利用基于模拟退火算法(SA)和粒子群优化算法(PSO)相结合的智能算法作为优化工具对模型背景值和初始条件中的自适应参数进行了优化提取,有效克服了现有非等间距GM(1,1)模型的不足,提高了非等间距GM(1,1)模型的预测精度,获得短时模拟交通流数据。
Description
技术领域
本发明涉及数据预测领域,更具体地,涉及一种基于SA-PSO的短时交通流预测方法。
背景技术
灰色预测模型是灰色***理论的一个重要组成部分,因其具有建模过程简单、所需样本少等优点,目前已成功应用到了电力、农业、管理、经济、高新技术产业等领域。灰色预测模型种类较多,由于众多原因,实际中存在很多原始数据不完整、样本序列为非等间距的情况,针对非等间距数据序列,非等间距GM(1,1)模型是应用最广泛的、最重要的灰色预测模型。传统非等间距GM(1,1)模型的预测精度易受背景值、初始条件的影响,当背景值和初始条件构造不合理时会出现预测失真、振荡现象,从而导致预测模型的适用性降低甚至不可用。因此,有必要开发一种预测精度高、适用性强的基于SA-PSO的短时交通流预测方法。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种基于SA-PSO的短时交通流预测方法,其能够通过对背景值进行重构,计算初始条件对应的最优时间参数,确定自适应参数,增加了模型使用的稳定性,提高了模型构建的理论深度;获得的短时模拟交通流数据具有适用性强、预测精度高、相对误差小的优点。
所述方法可以包括如下过程:步骤1:获取车联网的短时交通流数据,建立非等间距GM(1,1)模型的原始数据序列X(0)=(x(0)(k1),x(0)(k2),…,x(0)(kn));步骤2:对所述原始数据序列X(0)作1-AGO处理,得到1-AGO序列X(1),其中,i、j为数据序号,Δkj=kj-kj-1,Δk1=1;步骤3:根据所述1-AGO序列X(1),生成背景值Z(1)、计算矩阵B和Y,获得所述非等间距GM(1,1)模型的表达式为x(0)(ki)Δki+az(1)(ki)=b,通过最小二乘参数估计获得参数序列为/>其中,a为发展系数,b为灰作用量,步骤4:通过模拟退火算法和粒子群优化算法,确定背景值和初始条件的最优自适应参数p*与δ*;步骤5:利用步骤4所得的最优自适应参数p*重构背景值,建立所述非等间距GM(1,1)模型,利用步骤4所得的最优自适应参数δ*重构初始条件,求解所述非等间距GM(1,1)模型的白化微分方程的时间响应函数,即为所述1-AGO序列X(1)的模拟预测序列/>通过累减还原获得所述原始数据序列X(0)的模拟预测序列/>即为短时模拟交通流数据。
优选地,步骤3中生成背景值Z(1)包括:
其中,q为计算参数,p为自适应参数,p为一个整数值,p≥1,/>为1-AGO序列X(1)的拟合函数,/> M=x(0)(k1)-N。
优选地,步骤4中的初始条件为:
其中,δ为自适应参数,θ*为时间参数的最优值,其计算方法为:G为计算参数,
优选地,所述步骤4进一步包括:设定误差阈值,以所述相对误差平方和指标小于所述误差阈值时的p和δ为最优自适应参数p*与δ*。
优选地,所述相对误差平方和指标为
优选地,所述1-AGO序列X(1)的模拟预测序列为:
优选地,所述原始数据序列X(0)的模拟预测序列为:
其有益效果在于:
(1)利用复化梯形公式对背景值进行了重构;基于新信息优先原理,提出了一种基于时间参数的变权值初始条件优化方法,利用相对误差平方和最小准则,给出了初始条件对应的最优时间参数求解方法;利用模拟退火算法(SA)和粒子群优化算法(PSO)相结合的智能算法自适应地确定模型中的自适应参数,本发明增加了模型使用的稳定性、提高了模型构建的理论深度;
(2)根据本发明预测的短时模拟交通流数据具有适用性强、预测精度高、相对误差小的优点。
本发明的方法具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的基于SA-PSO的短时交通流预测方法的步骤的流程图。
图2示出了根据本发明的一个实施例的基于SA-PSO算法确定自适应参数的流程图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的基于SA-PSO的短时交通流预测方法的步骤的流程图。
在该实施例中,根据本发明的基于SA-PSO的短时交通流预测方法可以包括如下过程:步骤1:获取车联网的短时交通流数据,建立非等间距GM(1,1)模型的原始数据序列X(0)=(x(0)(k1),x(0)(k2),…,x(0)(kn));步骤2:对原始数据序列X(0)作1-AGO处理,得到1-AGO序列X(1),其中,i、j为数据序号,Δkj=kj-kj-1,Δk1=1;步骤3:根据1-AGO序列X(1),生成背景值Z(1)、计算矩阵B和Y,获得非等间距GM(1,1)模型的表达式为x(0)(ki)Δki+az(1)(ki)=b,通过最小二乘参数估计获得参数序列为/>其中,a为发展系数,b为灰作用量,/>步骤4:通过模拟退火算法和粒子群优化算法,确定背景值和初始条件的最优自适应参数p*与δ*;步骤5:利用步骤4所得的最优自适应参数p*重构背景值,建立非等间距GM(1,1)模型,利用步骤4所得的最优自适应参数δ*重构初始条件,求解非等间距GM(1,1)模型的白化微分方程的时间响应函数,即为1-AGO序列X(1)的模拟预测序列/>通过累减还原获得原始数据序列X(0)的模拟预测序列即为短时模拟交通流数据。
在一个示例中,步骤3中生成背景值Z(1)包括:
其中,q为计算参数,p为自适应参数,p为一个整数值,p≥1,/>为1-AGO序列X(1)的拟合函数,/>
M=x(0)(k1)-N。
在一个示例中,步骤4中的初始条件为:
其中,δ为自适应参数,θ*为时间参数的最优值,其计算方法为:G为计算参数,
在一个示例中,步骤4进一步包括:设定误差阈值,以所述相对误差平方和指标小于所述误差阈值时的p和δ为最优自适应参数p*与δ*。
在一个示例中,相对误差平方和指标为
在一个示例中,1-AGO序列X(1)的模拟预测序列为:
在一个示例中,原始数据序列X(0)的模拟预测序列为:
具体地,根据本发明的基于SA-PSO的短时交通流预测方法可以包括:
步骤1:获取车联网的短时交通流数据,建立非等间距预测模型的原始非负数据序列X(0)=(x(0)(k1),x(0)(k2),…,x(0)(kn))。
步骤2:针对原始数据序列作1-AGO处理,得到1-AGO序列X(1)为Δkj=kj-kj-1,Δk1=1,Δkj不为常数。
步骤3:根据1-AGO序列,生成背景值为计算矩阵为B和Y,其中,/>非等间距预测模型的表达式为x(0)(ki)Δki+az(1)(ki)=b,参数序列为/>
步骤4:通过模拟退火算法和粒子群优化算法相结合的智能算法,确定使相对误差平方和指标最小的最优自适应参数p*与δ*,其原理为:步骤401:初始化PSO算法群体,初始化各个粒子的速度V和位置等参数,确定各粒子初始极值为pBest=X;步骤402:建立非等间距GM(1,1)模型,利用SSRE指标计算各粒子的适应度值大小;步骤403:更新每个粒子的最优位置pBest和群体最优全局极值gBest;步骤404:更新各粒子速度V和计算各粒子位置Xnew=X+V;步骤405:利用模拟退火算法的Metropolis准则更新各粒子位置;步骤406:判断是否满足算法结束条件,即相对误差平方和指标是否为小于设定的误差阈值)或达到最大迭代次数,若是,计算结束,输出gBest,若否,则返回步骤402,进行迭代。最终获得的gBest即为最优自适应参数p*与δ*。其中,计算相对误差平方和指标为
步骤5:利用步骤4所得的最优自适应参数p*重构背景值,建立非等间距GM(1,1)模型,利用步骤4所得的最优自适应参数δ*重构初始条件为求解非等间距GM(1,1)模型的白化微分方程的时间响应函数,即为1-AGO序列X(1)的模拟预测序列/>为/>通过累减还原获得原始数据序列X(0)的模拟预测序列/>为/>即为对应的短时模拟交通流数据,其中,白化微分方程为/>
本方法通过对背景值和初始条件进行重构,并利用SA-PSO确定背景值和初始条件中的自适应参数,并计算初始条件对应的最优时间参数,增加了模型使用的稳定性,提高了模型构建的理论深度;获得的短时模拟交通流数据具有适用性强、预测精度高、相对误差小的优点。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
以计算车联网的短时交通流预测为例,包括:
步骤1:以非等间距的交通流量为x(0)(ki),建立非等间距预测模型的原始数据序列X(0)=(x(0)(k1),x(0)(k2),…,x(0)(kn))。
步骤2:针对原始数据序列作1-AGO处理,得到1-AGO序列X(1)为
步骤3:根据1-AGO序列,生成背景值为计算矩阵为B和Y,其中,/>非等间距预测模型的表达式为x(0)(ki)Δki+az(1)(ki)=b,参数序列为/>
图2示出了根据本发明的一个实施例的基于SA-PSO算法确定最优自适应参数的流程图。
步骤4:通过模拟退火算法和粒子群优化算法,确定使相对误差平方和指标最小的最优自适应参数p*与δ*,流程图如图2所示。其中根据对应的原始数据序列的模拟预测序列,计算相对误差平方和指标为
步骤5:利用步骤4所得的最优自适应参数p*重构背景值,建立非等间距GM(1,1)模型,利用步骤4所得的最优自适应参数δ*重构初始条件为非等间距GM(1,1)模型的白化微分方程为/>求解白化微分方程的时间响应函数,即为1-AGO序列X(1)的模拟预测序列/>为/>通过累减还原获得原始数据序列X(0)的模拟预测序列/>为即为短时模拟交通流量。
综上所述,本发明通过对背景值进行重构,计算初始条件对应的最优时间参数,确定自适应参数,增加了模型使用的稳定性,提高了模型构建的理论深度;获得的短时模拟交通流数据具有适用性强、预测精度高、相对误差小的优点。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (2)
1.一种基于SA-PSO的短时交通流预测方法,其特征在于,包括如下过程:
步骤1:获取车联网的短时交通流数据,建立非等间距GM(1,1)模型的原始数据序列X(0)=(x(0)(k1),x(0)(k2),…,x(0)(kn));
步骤2:对所述原始数据序列X(0)作1-AGO处理,得到1-AGO序列X(1),其中,为数据序号,Δkj=kj-kj-1,Δk1=1;
步骤3:根据所述1-AGO序列X(1),生成背景值Z(1)、计算矩阵B和Y,获得所述非等间距GM(1,1)模型的表达式为x(0)(ki)Δki+az(1)(ki)=b,通过最小二乘参数估计获得参数序列为其中,a为发展系数,b为灰作用量,/>
步骤4:通过模拟退火算法和粒子群优化算法,确定背景值和初始条件的最优自适应参数p*与δ*;
步骤5:利用步骤4所得的最优自适应参数p*重构背景值,建立所述非等间距GM(1,1)模型,利用步骤4所得的最优自适应参数δ*重构初始条件,求解所述非等间距GM(1,1)模型的白化微分方程的时间响应函数,即为所述1-AGO序列X(1)的模拟预测序列通过累减还原获得所述原始数据序列X(0)的模拟预测序列/>即为短时模拟交通流数据;
其中,步骤3中生成背景值Z(1)包括:
其中,q为计算参数,p为自适应参数,p为一个整数值,p≥1,/>为1-AGO序列X(1)的拟合函数,/> M=x(0)(k1)-N;
其中,步骤4中的初始条件为:
其中,δ为自适应参数,θ*为时间参数的最优值,其计算方法为:G为计算参数,
所述步骤4进一步包括:
设定误差阈值,以相对误差平方和指标小于所述误差阈值时的p和δ为最优自适应参数p*与δ*;
所述相对误差平方和指标为
所述原始数据序列X(0)的模拟预测序列为:
2.根据权利要求1所述的基于SA-PSO的短时交通流预测方法,其特征在于,所述1-AGO序列X(1)的模拟预测序列为:/>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811469397.3A CN109740795B (zh) | 2018-11-28 | 2018-11-28 | 基于sa-pso的短时交通流预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811469397.3A CN109740795B (zh) | 2018-11-28 | 2018-11-28 | 基于sa-pso的短时交通流预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109740795A CN109740795A (zh) | 2019-05-10 |
CN109740795B true CN109740795B (zh) | 2023-09-29 |
Family
ID=66359125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811469397.3A Active CN109740795B (zh) | 2018-11-28 | 2018-11-28 | 基于sa-pso的短时交通流预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109740795B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103730006A (zh) * | 2014-01-26 | 2014-04-16 | 吉林大学 | 一种短时交通流量的组合预测方法 |
KR20140121740A (ko) * | 2013-04-08 | 2014-10-16 | 연세대학교 산학협력단 | 모관 흡수력 측정을 통한 비탈면의 안전도 모니터링 장치 및 방법 |
CN107862108A (zh) * | 2017-10-12 | 2018-03-30 | 成都阜特科技股份有限公司 | 一种工业机械设备健康状态分析与预测方法及其*** |
CN107918808A (zh) * | 2017-11-16 | 2018-04-17 | 淮阴工学院 | 基于反余弦函数与幂函数变换与背景值组合优化的灰色预测模型方法 |
-
2018
- 2018-11-28 CN CN201811469397.3A patent/CN109740795B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140121740A (ko) * | 2013-04-08 | 2014-10-16 | 연세대학교 산학협력단 | 모관 흡수력 측정을 통한 비탈면의 안전도 모니터링 장치 및 방법 |
CN103730006A (zh) * | 2014-01-26 | 2014-04-16 | 吉林大学 | 一种短时交通流量的组合预测方法 |
CN107862108A (zh) * | 2017-10-12 | 2018-03-30 | 成都阜特科技股份有限公司 | 一种工业机械设备健康状态分析与预测方法及其*** |
CN107918808A (zh) * | 2017-11-16 | 2018-04-17 | 淮阴工学院 | 基于反余弦函数与幂函数变换与背景值组合优化的灰色预测模型方法 |
Non-Patent Citations (6)
Title |
---|
Reliability prediction using an unequal interval grey model;Yuhong Wang et al;《2010 Proceedings - Annual Reliability and Maintainability Symposium (RAMS)》;1-2 * |
基于GA-SAPSO优化BP神经网络的光伏***发电量预测;黄超等;《宿州学院学报》;20160331;第31卷(第3期);110-114 * |
基于SA-PSO的配电网电容器优化配置;杨洪等;《电力***及其自动化学报》;20090630;第21卷(第3期);111-114 * |
基于广义灰色模型的极限承载力建模与预测研究;李军亮;《中国博士学位论文全文数据库(电子期刊)》;20100115;第2010年卷(第01期);第4章 * |
基于组合插值的GM(1,1)模型背景值的改进;李凯等;《计算机应用研究》;20181031;第35卷(第10期);第2.1节 * |
李军亮.基于广义灰色模型的极限承载力建模与预测研究.《中国博士学位论文全文数据库(电子期刊)》.2010,第2010年卷(第01期),第4章. * |
Also Published As
Publication number | Publication date |
---|---|
CN109740795A (zh) | 2019-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109243538B (zh) | 一种预测疾病与LncRNA关联关系的方法及*** | |
KR101958674B1 (ko) | 시퀀스 재귀 필터링 3차원 변분(3d-var) 기반의 실측 해양 환경 데이터 동화방법 | |
CN103730006B (zh) | 一种短时交通流量的组合预测方法 | |
CN108763167A (zh) | 一种变分贝叶斯的自适应滤波方法 | |
CN110210495A (zh) | 基于并行LSTM自编码器动态特征提取的XGBoost软测量建模方法 | |
CN106875047B (zh) | 水库流域径流预报方法及*** | |
CN110674996B (zh) | 一种城市交通噪声预测方法 | |
CN107832789B (zh) | 基于平均影响值数据变换的特征加权k近邻故障诊断方法 | |
CN110930996B (zh) | 模型训练方法、语音识别方法、装置、存储介质及设备 | |
CN108681689B (zh) | 基于生成对抗网络的帧率增强步态识别方法及装置 | |
CN112561881A (zh) | 基于评价模型的红外图像自适应数据增强方法 | |
Sobol et al. | Derivative based global sensitivity measures | |
CN115937693A (zh) | 一种基于遥感图像的道路识别方法及*** | |
CN109740795B (zh) | 基于sa-pso的短时交通流预测方法 | |
CN117332688B (zh) | 随机荷载作用下x80管道钢疲劳裂纹扩展预测方法及*** | |
Singh et al. | Modified mean square error algorithm with reduced cost of training and simulation time for character recognition in backpropagation neural network | |
CN116935964A (zh) | 基于深度学习模型预测抗体翻译后修饰位点的方法和*** | |
CN101872482B (zh) | 一种基于免疫遗传粒子滤波的智能视频跟踪方法 | |
CN110751400B (zh) | 一种风险评估方法及装置 | |
CN113360896B (zh) | 一种横向联邦学习架构下的Free Rider攻击检测方法 | |
Lalande et al. | A new framework to accurately quantify soil bacterial community diversity from DGGE | |
WO2019149376A1 (en) | Method and system for processing input data using a neural network and normalizations | |
CN114819398A (zh) | 基于灰色马尔可夫链的北斗卫星钟差序列组合预测方法 | |
CN114942480A (zh) | 基于信息感知注意力动态协同网络的海洋站风速预报方法 | |
CN113870341A (zh) | 基于rgb和激光特征融合的高炉烧结矿粒径检测方法及*** |
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 |