CN107346460B - 基于智能网联系统下前车运行信息的未来工况预测方法 - Google Patents
基于智能网联系统下前车运行信息的未来工况预测方法 Download PDFInfo
- Publication number
- CN107346460B CN107346460B CN201710585795.0A CN201710585795A CN107346460B CN 107346460 B CN107346460 B CN 107346460B CN 201710585795 A CN201710585795 A CN 201710585795A CN 107346460 B CN107346460 B CN 107346460B
- Authority
- CN
- China
- Prior art keywords
- vehicle
- prediction
- working condition
- front vehicle
- model
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000012937 correction Methods 0.000 claims abstract description 55
- 238000012843 least square support vector machine Methods 0.000 claims description 34
- 238000004891 communication Methods 0.000 claims description 25
- 238000012549 training Methods 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000012216 screening Methods 0.000 claims description 5
- 239000002245 particle Substances 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 3
- 238000004904 shortening Methods 0.000 claims description 3
- 238000012706 support-vector machine Methods 0.000 claims description 2
- 230000006870 function Effects 0.000 description 11
- 238000013528 artificial neural network Methods 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman problem
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L67/00—Network arrangements or protocols for supporting network services or applications
- H04L67/01—Protocols
- H04L67/12—Protocols specially adapted for proprietary or special-purpose networking environments, e.g. medical networks, sensor networks, networks in vehicles or remote metering networks
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computing Systems (AREA)
- Evolutionary Biology (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Educational Administration (AREA)
- Algebra (AREA)
- Bioinformatics & Computational Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Medical Informatics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于智能网联系统下前车运行信息的未来工况预测方法,属于智能交通技术领域,包括前车行驶工况数据获取、前车行驶工况划分、未来工况预测模型建立和本车未来工况在线预测,该工况预测方法能准确获取与本车未来工况信息最为接近的前车工况信息,并以此为基础通过建立最小二乘支持向量机与自回归滑动平均误差修正相结合带有预测模型精度判断的未来工况预测模型进行本车未来工况预测,预测结果具有实时性好、准确率高、可参考性强的优点,特别适合应用在固定线路运行的混合动力车辆。
Description
技术领域
本发明提供一种汽车行驶工况预测方法,特别涉及一种基于智能网联系统下前车运行信息的未来工况预测方法,属于智能交通技术领域。
背景技术
行驶工况是混合动力汽车能量管理策略设计考虑的重要因素之一,对提高整车燃油经济性有着至关重要的作用。开发对未来控制时域内工况进行合理精确预测的方法,进而结合预测能量管理算法实现混合动力系统实时最优控制,已成为混动汽车智能能量管理策略的有效方法。目前行驶工况预测的主要研究是根据车辆自身行驶一定周期后的工况数据后,结合建立的预测模型作出对车辆未来工况的预测。由于依赖车辆运行一定周期的历史数据积累之后作出的工况预测,未来工况预测结果存在滞后性、准确率低,参考性差等问题。如2013年8月14日申请公布的发明专利:申请公布号:CN 103246943A,基于马尔可夫链的汽车运行工况多尺度预测方法,该方法建立汽车运行工况的马尔可夫链预测模型,根据汽车运行工况的历史信息,通过极大似然估计计算出状态转移矩阵;运用马尔可夫链蒙特卡洛模拟方法,根据获得的状态转移矩阵进行不同时间尺度的汽车运行工况预测;再将不同尺度预测结果在原数据频率下融合,获得汽车运行工况多尺度预测结果。该方法基于汽车自身运行工况的历史信息,通过建立马尔可夫链预测模型完成对汽车工况的多尺度预测,由于汽车自身运行工况的历史信息的对未来工况预测的结果实时性差,且工况预测模型预测结果缺少对误差的深度分析的,不能保证预测模型和预测结果两者预测精度,因此未来工况预测结果存在滞后性、准确率低,参考性差等问题。
发明内容
本发明的目的是提供一种能准确的获取与本车未来工况信息最为接近的前车工况信息,并以此为基础进行未来工况预测,其预测结果具有实时性好、准确性高,可参考性强的基于智能网联系统下前车运行信息的未来工况预测方法,其技术内容为:
基于智能网联系统下前车运行信息的未来工况预测方法,包括前车行驶工况数据获取、前车行驶工况划分、未来工况预测模型建立和本车未来工况在线预测,其特征在于:
第一步,前车行驶工况数据获取:首先通过V2V车车通讯系统获取周围车辆行驶信息,通过V2I车路通讯系统获取交通路况信息,以及通过车载定位系统获取车辆当前位置信息及行驶路径道路信息,进行工况数据信息的采集;
再确定工况数据获取的最佳途径,根据V2V车车通讯系统获取前车的车距、行驶方向及道路信息进行分析,当本车与前车同向同路且距离小于S时,首先判断前车类型与本车类型是否为同类型车辆,若前车类型与本车类型不属于同类型车辆,则通过对前车类型进行筛选,匹配与本车类型最相近的车辆类型的前车作为目标前车,若前车与本车为同类型车辆且数量M=1时,直接将前车作为目标对象,若前车与本车为同类型车辆且数量M>1时,通过前车此前与本车在相同路段内的工况信息比对确定前车筛选的优先级,优先选定工况信息相似度较高的前车作为目标前车;然后对本车当前行驶时交通流与目标前车之前相同路段时交通流对比,若两者差值在允许值范围内,利用V2V车车通讯系统获取前车的工况信息;若两者差值超过允许值,由V2V车车通讯系统获取前车工况信息,同时结合V2I车路通讯系统获取实时交通流信息,对前车车速工况信息加以修正;当不满足本车与前车同向同路且距离小于S时,利用V2I车路通讯系统获取由远程监控平台系统提供的上一时刻前车工况数据;
所述的同类型车辆为同款同型号完全相同的车型,所述的最相近的车辆类型是指整车质量、动力源部件功率及滚阻都近似相同的车型;
第二步,前车行驶工况的划分:将获取的前车行驶工况信息,按照行驶工况路网划分为道路交叉路口和连接交叉路口的路段,再结合V2I车路通讯系统获取的数据信息分别确定连接路段和道路交叉口划分时间窗口长度;
所述的连接路段划分时间窗口长度的确定,先根据连接路段内道路路面等级划分,在相同的道路路面等级下再按照交通流密度等级及通行速度等级的一致性进行划分,逐步缩短连接路段划分时间窗口长度,再结合连接路段划分时间窗口长度限值,进而确定连接路段划分时间窗口长度;所述的道路交叉口划分时间窗口长度的确定,将车辆通过交叉路口的通行状态分为车辆加速通过交叉口和车辆停车起步通过交叉口两类,分别结合此时交通流密度等级及交叉口划分时间窗口长度限值,进而确定交叉口划分时间窗口长度;
第三步,未来工况预测模型建立:依据划分的前车行驶工况,收集划分时间窗口内的前车实时运行工况信息,以及本车预测时间步长Δt内的未来实际工况信息数据,进行滤波和参数化归一处理后,将当前接收到得前车的工况信息数据作为输入,本车预测时间步长Δt内的未来实际工况信息数据作为输出,构造训练样本集G;
建立基于最小二乘支持向量机模型(LS-SVM)与自回归滑动平均误差修正模型(ARMA)相结合的工况预测模型:先确定最小二乘支持向量机模型(LS-SVM)函数中核函数,并利用前车及本车工况数据的训练样本集作为粒子,利用粒子群优化算法确定核宽度σ和惩罚参数C;再利用训练样本集G对最小二乘支持向量机模型(LS-SVM)进行训练,确定拉格朗日乘子α和偏置值b,最终得到最小二乘支持向量机模型(LS-SVM)函数;利用最小二乘支持向量机模型(LS-SVM)预测的工况数据与实际工况数据两者差值,组成误差序列样本,根据误差序列样本自相关系数的截尾性和偏自相关系数拖尾性,初步确定自回归滑动平均误差修正模型(ARMA)阶数,再用最小二乘估计法确定模型中的自回归参数、滑动平均参数,最终得到自回归滑动平均误差修正模型(ARMA);
预测模型精度判断:将利用最小二乘支持向量机模型(LS-SVM)得到的工况预测值与自回归滑动平均误差修正模型(ARMA)得到的误差修正值Ei两者之和作为预测调整值与工况实际值yi进行比较,当满足预测精度要求时,将预测调整值作为最终预测值;当预测精度不满足要求,进一步调整自回归滑动平均误差修正模型(ARMA)得到新的误差修正值Ei(k),进而得到新的预测调整值再与工况实际值yi比较并进行预测精度分析,直到满足预测精度要求,最终得到最小二乘支持向量机模型(LS-SVM)与自回归滑动平均误差修正模型(ARMA)相结合带有预测模型精度判断的未来工况预测模型;
第四步,本车未来工况在线预测:将实时获取的划分时间窗口内的前车行驶工况数据作为输入,利用建立的未来工况预测模型进行工况预测,得到本车未来工况预测结果,再将本车预测结果和本车未来工况的实际结果,输入残差计算模型并进行预测精度判断,当不满足预测精度要求时,通过再次建立的自回归滑动平均误差修正模型(ARMA),在下一段的工况预测时,优先通过自回归滑动平均误差修正模型(ARMA)对所得的在线预测结果进行误差修正后,作为下一段工况的预测结果;当满足预测精度要求时,在下一段的工况预测时不进行误差修正,直接将当前预测结果输出;
经过误差修正模型对在线预测结果进行误差修正后,再进行残差计算及预测精度判断时,以误差修正后的预测结果为基础进行分析,重复上述精度判断后的处理方法。
本发明与现有技术相比,有益效果如下:
(1)相比于依赖车辆自身运行一定周期的历史数据积累之后作出的工况预测,该工况预测方法采用的工况数据获取途径,能较准确地获取与本车未来工况最为接近的前车工况信息,保证了预测模型输入数据的实时性和准确性,以此为基础进行工况预测,结果更加接近未来实际工况;
(2)该工况预测方法采用了基于行驶工况路网与交通信息相融合的行驶工况划分方法,针对连接路段和道路交叉分别确定口划分时间窗口长度,保证预测准确度的同时降低了预测的计算量;
(3)该工况预测方法还从其他两方面保证了预测结果的精确性:①建立了最小二乘支持向量机模型(LS-SVM)与自回归滑动平均误差修正模型(ARMA)相结合带有预测模型精度判断的未来工况预测模型,预测模型中加入了误差修正模型和预测模型精度判断,较好的保证了预测模型精度;②未来工况在线预测时进一步对当前在线预的测结果进行预精度分析,通过误差预测修正模型对下时间段的工况预测进行修正和调整,进一步保证预测结果准确性;
(4)由于固定线路运行的混合动力车辆在利用本方法进行前车行驶工况数据获取、工况划分时更为准确,因此该方法特别适合应用在固定线路运行的混合动力车辆,如固定线路运行的营运及作业的混合动力车辆,如采用混合动力系统的公交车辆、固定线路物流车、清扫作业车辆等。
附图说明
图1是本发明实施例的工况预测方法整体流程图。
图2是本发明实施例行驶工况数据获取途径的流程图。
图3是本发明实施例前车此前Ty段与本车当前Tn段两相同路段的示意图。
图4是本发明实施例行驶工况划分时间窗口长度的流程图。
具体实施方式
下面结合附图对本发明作进一步说明:
如图1所示,基于智能网联系统下前车运行信息的未来工况预测方法,包括前车行驶工况数据获取、前车行驶工况划分、未来工况预测模型建立和本车未来工况在线预测,其特征在于:
第一步,前车行驶工况数据获取:首先通过V2V车车通讯系统获取周围车辆行驶信息,通过V2I车路通讯系统获取交通路况信息,以及通过车载定位系统获取车辆当前位置信息及行驶路径道路信息,进行工况数据信息的采集;
再确定工况数据获取的最佳途径,如图2所示,根据V2V车车通讯系统获取前车的车距、行驶方向及道路信息进行分析,当本车与前车同向同路且距离小于S时,首先判断前车类型与本车类型是否为同类型车辆,若前车类型与本车类型不属于同类型车辆,则通过对前车类型进行筛选,匹配与本车类型最相近的车辆类型的前车作为目标前车,若前车与本车为同类型车辆且数量M=1时,直接将前车作为目标对象,若前车与本车为同类型车辆且数量M>1时,通过前车此前与本车在相同路段内的工况信息比对确定前车筛选的优先级,优先选定工况信息相似度较高的前车作为目标前车;然后通过V2I车路通讯系统获取目标前车此前Ty段与本车当前Tn段两相同路段的交通流信息,如图3所示,包括目标前车此前Ty段的交通流密度Ky、交通流流量Qy,本车当前Tn段的交通流密度Kn、交通流流量Qn,并计算两者交通流密度的差值Ek、交通流流量差值Eq,分别与设定交通流密度变化情况临界值k、交通流流量变化情况临界值q进行比较;若交通流变化较小即:Ek≤k且Eq≤q,利用V2V车车通讯系统获取前车的工况信息;若交通流变化较大即:Ek>k或Eq>q,由V2V车车通讯系统获取前车工况信息,同时结合V2I车路通讯系统获取实时交通流信息,建立基于交通流变化的车速修正模型,对获取的前车车速工况信息加以修正;
当不满足本车与前车同向同路且距离小于S时,利用V2I车路通讯系统获取由远程监控平台系统提供的上一时刻前车工况数据;
所述的同类型车辆为同款同型号完全相同的车型,所述的最相近的车辆类型是指整车质量、动力源部件功率及滚阻都近似相同的车型;
所述的相同路段内的工况信息对比,是在相同路段内对车辆瞬时速度vt、平均车速vave、最大速度vmax、车速变化频率f、路面坡度i、路面等级g、最大加速度amax、加速度均值am进行的比对;
所述的基于交通流变化的车速修正模型,可通过构建基于交通流变化的RBF神经网络车速修正模型得到,RBF神经网络车速修正模型的构建包括:(1)确定RBF神经网络车速修正模型输入参数矢量和输出参数矢量,输入参数矢量为当前接收到得前车工况信息包括的瞬时速度vy_t,平均车速vy_ave、最大速度vy_max、车速变化频率fy,以及本车当前路段交通流信息包括交通流密度Kn、交通流流量Qn,即{vy_t,vy_ave,vy_max,fy,Kn,Qn},输出参数矢量为本车所对应的未来实际工况的车速信息,包括本车瞬时速度vn_t,平均车速vn_ave、最大速度vn_max、车速变化频率fn,即{vn_t,vn_ave,vn_max,fn};(2)将输入参数矢量和输出参数矢量作为训练样本,输入到RBF神经网络模型中进行离线训练,选用自组织选取中心的RBF神经网络学习方法,求解确定隐含层基函数中心、奇函数的方差和隐含层单元输出单元权值,最终建立的基于交通流变化的RBF神经网络车速修正模型;
第二步,前车行驶工况的划分:将获取的前车行驶工况信息,按照行驶工况路网划分为道路交叉路口和连接交叉路口的连接路段,再结合V2I车路通讯系统获取的数据信息分别确定连接路段和道路交叉口划分时间窗口长度,如图4所示;
所述的连接路段划分时间窗口长度的确定,先根据连接路段内道路路面等级划分,在相同的道路路面等级下再按照交通流密度等级及通行速度等级的一致性进行划分,逐步缩短连接路段划分时间窗口长度,再结合连接路段划分时间窗口长度限值,进而确定连接路段划分时间窗口长度;所述的道路交叉口划分时间窗口长度的确定,将车辆通过交叉路口的通行状态分为车辆加速通过交叉口和车辆停车起步通过交叉口两类,分别结合此时交通流密度等级及交叉口划分时间窗口长度限值,进而确定交叉口划分时间窗口长度;
所述的连接路段和交叉口划分时间窗口长度限值,可分别对时间窗口长度从t1~t2进行预测精确度仿真分析,综合考虑不同时间窗口长度对计算量的影响,分别确定路段和交叉口划分时间窗口长度的上限值和下限值;
第三步,未来工况预测模型建立:依据划分的前车行驶工况,收集某段划分时间窗口内的前车实时运行工况信息,以及本车预测时间步长Δt内的未来实际工况信息数据,进行滤波和参数化归一处理后,将当前接收到得前车的工况信息数据作为输入,本车预测时间步长Δt内的未来实际工况信息数据作为输出,构造训练样本集G,所述的工况信息包括车辆瞬时速度vt、平均车速vave、最大速度vmax、车速变化频率f、路面坡度i、路面等级g、交通流密度ρ,最大加速度amax、加速度均值am;
由此构造训练样本集为:G={(x1,y1),(x2,y2),……(xi,yi),……(xn,yn)},其中,当前接收到得前车工况信息为:本车预测时间步长Δt内的未来的实际工况信息:
建立基于最小二乘支持向量机模型(LS-SVM)与自回归滑动平均误差修正模型(ARMA)相结合的工况预测模型:
由最小二乘支持向量机理论,最优回归估计函数是在一定约束条件下的最小化泛函,即:
约束条件:
利用拉格朗日乘子法将问题转化到其对偶空间并进行优化求解,可得LS-SVM回归估计函数模型为:
其中,为LS-SVM的核函数,选定核函数为RBF径向基核函数,由此可得:
并利用前车及本车工况数据的训练样本集作为粒子,利用粒子群优化算法确定核宽度σ和惩罚参数C;再利用训练样本集G对最小二乘支持向量机模型(LS-SVM)进行训练,确定拉格朗日乘子α和偏置值b,最终得到最小二乘支持向量机模型(LS-SVM)函数;
根据自回归滑动平均模型理论,可得ARMA(p,q)模型表达式为
yt=φ1yt-1+φ2yt-2+…+φpyt-p+εt-γ1εt-1-γ2εt-2-…-γqεt-q
其中,(p,q)为自回归滑动平均模型的阶,φ1,φ2,…,φp,为自回归参数,γ1,γ2,…,γq为滑动平均参数,εt为高斯白噪声序列;
利用最小二乘支持向量机模型(LS-SVM)预测的工况数据与实际工况数据两者差值,组成误差序列样本,根据误差序列样本自相关系数的截尾性和偏自相关系数拖尾性,初步确定自回归滑动平均误差修正模型(ARMA)阶数(p,q),再用最小二乘估计法确定模型中的自回归参数、滑动平均参数,最终得到自回归滑动平均误差修正模型(ARMA);
预测模型精度判断:将利用最小二乘支持向量机模型(LS-SVM)得到的工况预测值与自回归滑动平均误差修正模型(ARMA)得到的误差修正值Ei两者之和作为预测调整值与工况实际值yi进行比较,当满足预测精度要求时,将预测调整值作为最终预测值;当预测精度不满足要求,进一步调整自回归滑动平均误差修正模型(ARMA)得到新的误差修正值Ei(k),进而得到新的预测调整值再与工况实际值yi,直到满足预测精度要求;
第四步,本车未来工况在线预测:将实时获取的划分时间窗口内的前车行驶工况数据作为输入,利用建立的未来工况预测模型进行工况预测,得到本车未来工况预测结果,再将本车预测结果和本车未来工况的实际结果,输入残差计算模型并进行预测精度判断,当不满足预测精度要求时,通过再次建立的自回归滑动平均误差修正模型(ARMA),在下一段的工况预测时,优先通过自回归滑动平均误差修正模型(ARMA)对所得的在线预测结果进行误差修正后,作为下一段工况的预测结果;当满足预测精度要求时,在下一段的工况预测时不进行误差修正,直接将当前预测结果输出;
经过误差修正模型对在线预测结果进行误差修正后,再进行残差计算及预测精度判断时,以误差修正后的预测结果为基础进行分析,重复上述精度判断后的处理方法。
Claims (1)
1.一种基于智能网联系统下前车运行信息的未来工况预测方法,包括前车行驶工况数据获取、前车行驶工况划分、未来工况预测模型建立和本车未来工况在线预测,其特征在于:
第一步,前车行驶工况数据获取:首先通过V2V车车通讯系统获取周围车辆行驶信息,通过V2I车路通讯系统获取交通路况信息,以及通过车载定位系统获取车辆当前位置信息及行驶路径道路信息,进行工况数据信息的采集;
再确定工况数据获取的最佳途径,根据V2V车车通讯系统获取前车的车距、行驶方向及道路信息进行分析,当本车与前车同向同路且距离小于S时,首先判断前车类型与本车类型是否为同类型车辆,若前车类型与本车类型不属于同类型车辆,则通过对前车类型进行筛选,匹配与本车类型最相近的车辆类型的前车作为目标前车,若前车与本车为同类型车辆且数量M=1时,直接将前车作为目标对象,若前车与本车为同类型车辆且数量M>1时,通过前车此前与本车在相同路段内的工况信息比对确定前车筛选的优先级,优先选定工况信息相似度较高的前车作为目标前车;然后对本车当前行驶时交通流与目标前车之前相同路段时交通流对比,若两者差值在允许值范围内,利用V2V车车通讯系统获取前车的工况信息;若两者差值超过允许值,由V2V车车通讯系统获取前车工况信息,同时结合V2I车路通讯系统获取实时交通流信息,对前车车速工况信息加以修正;当不满足本车与前车同向同路且距离小于S时,利用V2I车路通讯系统获取由远程监控平台系统提供的上一时刻前车工况数据;
所述的同类型车辆为同款同型号完全相同的车型,所述的最相近的车辆类型是指整车质量、动力源部件功率及滚阻都近似相同的车型;
第二步,前车行驶工况的划分:将获取的前车行驶工况信息,按照行驶工况路网划分为道路交叉路口和连接交叉路口的路段,再结合V2I车路通讯系统获取的数据信息分别确定连接路段和道路交叉口划分时间窗口长度;
所述的连接路段划分时间窗口长度的确定,先根据连接路段内道路路面等级划分,在相同的道路路面等级下再按照交通流密度等级及通行速度等级的一致性进行划分,逐步缩短连接路段划分时间窗口长度,再结合连接路段划分时间窗口长度限值,进而确定连接路段划分时间窗口长度;所述的道路交叉口划分时间窗口长度的确定,将车辆通过交叉路口的通行状态分为车辆加速通过交叉口和车辆停车起步通过交叉口两类,分别结合此时交通流密度等级及交叉口划分时间窗口长度限值,进而确定交叉口划分时间窗口长度;
第三步,未来工况预测模型建立:依据划分的前车行驶工况,收集划分时间窗口内的前车实时运行工况信息,以及本车预测时间步长Δt内的未来实际工况信息数据,进行滤波和参数化归一处理后,将当前接收到得前车的工况信息数据作为输入,本车预测时间步长Δt内的未来实际工况信息数据作为输出,构造训练样本集G;
建立基于最小二乘支持向量机模型(LS-SVM)与自回归滑动平均误差修正模型(ARMA)相结合的工况预测模型:先确定最小二乘支持向量机模型(LS-SVM)函数中核函数,并利用前车及本车工况数据的训练样本集作为粒子,利用粒子群优化算法确定核宽度σ和惩罚参数C;再利用训练样本集G对最小二乘支持向量机模型(LS-SVM)进行训练,确定拉格朗日乘子α和偏置值b,最终得到最小二乘支持向量机模型(LS-SVM)函数;利用最小二乘支持向量机模型(LS-SVM)预测的工况数据与实际工况数据两者差值,组成误差序列样本,根据误差序列样本自相关系数的截尾性和偏自相关系数拖尾性,初步确定自回归滑动平均误差修正模型(ARMA)阶数,再用最小二乘估计法确定模型中的自回归参数、滑动平均参数,最终得到自回归滑动平均误差修正模型(ARMA);
预测模型精度判断:将利用最小二乘支持向量机模型(LS-SVM)得到的工况预测值与自回归滑动平均误差修正模型(ARMA)得到的误差修正值Ei两者之和作为预测调整值与工况实际值yi进行比较,当满足预测精度要求时,将预测调整值作为最终预测值;当预测精度不满足要求,进一步调整自回归滑动平均误差修正模型(ARMA)得到新的误差修正值Ei(k),进而得到新的预测调整值再与工况实际值yi比较并进行预测精度分析,直到满足预测精度要求,最终得到最小二乘支持向量机模型(LS-SVM)与自回归滑动平均误差修正模型(ARMA)相结合带有预测模型精度判断的未来工况预测模型;
第四步,本车未来工况在线预测:将实时获取的划分时间窗口内的前车行驶工况数据作为输入,利用建立的未来工况预测模型进行工况预测,得到本车未来工况预测结果,再将本车预测结果和本车未来工况的实际结果,输入残差计算模型并进行预测精度判断,当不满足预测精度要求时,通过再次建立的自回归滑动平均误差修正模型(ARMA),在下一段的工况预测时,优先通过自回归滑动平均误差修正模型(ARMA)对所得的在线预测结果进行误差修正后,作为下一段工况的预测结果;当满足预测精度要求时,在下一段的工况预测时不进行误差修正,直接将当前预测结果输出;
经过误差修正模型对在线预测结果进行误差修正后,再进行残差计算及预测精度判断时,以误差修正后的预测结果为基础进行分析,重复上述精度判断后的处理方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710585795.0A CN107346460B (zh) | 2017-07-18 | 2017-07-18 | 基于智能网联系统下前车运行信息的未来工况预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710585795.0A CN107346460B (zh) | 2017-07-18 | 2017-07-18 | 基于智能网联系统下前车运行信息的未来工况预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107346460A CN107346460A (zh) | 2017-11-14 |
CN107346460B true CN107346460B (zh) | 2018-06-12 |
Family
ID=60257094
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710585795.0A Expired - Fee Related CN107346460B (zh) | 2017-07-18 | 2017-07-18 | 基于智能网联系统下前车运行信息的未来工况预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107346460B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107909118B (zh) * | 2017-12-11 | 2022-02-22 | 北京映翰通网络技术股份有限公司 | 一种基于深度神经网络的配电网工况录波分类方法 |
CN110286584A (zh) * | 2018-03-19 | 2019-09-27 | 罗伯特·博世有限公司 | 机动车降温控制系统和方法 |
CN109017809B (zh) * | 2018-08-27 | 2020-05-19 | 北京理工大学 | 一种基于越野工况预测的能量分配方法 |
CN109754175B (zh) * | 2018-12-28 | 2023-04-07 | 广州明动软件股份有限公司 | 用于对行政审批事项的办结时限进行压缩预测的计算模型及其应用 |
CN109946730B (zh) * | 2019-03-06 | 2022-05-13 | 东南大学 | 一种车路协同下基于超宽带的车辆高可靠融合定位方法 |
KR20240116961A (ko) * | 2019-04-09 | 2024-07-30 | 닛산 지도우샤 가부시키가이샤 | 정보 처리 장치, 정보 처리 방법, 및 서버 |
CN110379165A (zh) * | 2019-07-26 | 2019-10-25 | 中国第一汽车股份有限公司 | 一种道路类型预测方法、装置、设备及存储介质 |
CN110610611B (zh) * | 2019-09-26 | 2021-08-03 | 江苏大学 | 面向混行交通流中智能网联车辆的行车安全性评价方法 |
CN111785024B (zh) * | 2020-07-17 | 2022-03-18 | 陕西工业职业技术学院 | 一种分区域分时域的城市车辆工况构建方法 |
CN112101415B (zh) * | 2020-08-13 | 2024-07-09 | 联合汽车电子有限公司 | 累碳量预测方法、装置、汽车、云端服务器及计算机可读存储介质 |
CN112580251B (zh) * | 2020-11-16 | 2022-07-12 | 北京理工大学 | 基于交通信息和模型预测控制的混动汽车能量管理方法 |
CN112668779B (zh) * | 2020-12-28 | 2024-07-12 | 重庆大学 | 基于自适应高斯过程的前车运动状态预测方法 |
CN112966860B (zh) * | 2021-02-23 | 2022-09-16 | 武汉育知联信息科技有限公司 | 一种基于附预链式数据分析预测系统 |
CN116843368B (zh) * | 2023-07-17 | 2024-01-26 | 杭州火奴数据科技有限公司 | 基于arma模型的营销数据处理方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004021468A (ja) * | 2002-06-14 | 2004-01-22 | Matsushita Electric Ind Co Ltd | 自動車、及び、走行制御システム |
CN102831768A (zh) * | 2012-08-15 | 2012-12-19 | 大连理工大学 | 一种基于车联网的混合动力客车行驶工况预测方法 |
CN103246943A (zh) * | 2013-05-31 | 2013-08-14 | 吉林大学 | 基于马尔可夫链的汽车运行工况多尺度预测方法 |
-
2017
- 2017-07-18 CN CN201710585795.0A patent/CN107346460B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004021468A (ja) * | 2002-06-14 | 2004-01-22 | Matsushita Electric Ind Co Ltd | 自動車、及び、走行制御システム |
CN102831768A (zh) * | 2012-08-15 | 2012-12-19 | 大连理工大学 | 一种基于车联网的混合动力客车行驶工况预测方法 |
CN103246943A (zh) * | 2013-05-31 | 2013-08-14 | 吉林大学 | 基于马尔可夫链的汽车运行工况多尺度预测方法 |
Non-Patent Citations (3)
Title |
---|
"Real time vehicle speed prediction using a Neural Network Traffic Model";Jungme Park 等;《The 2011 International Joint Conference on Neural Networks》;20111003;第2991-2996页 * |
"基于马尔可夫链的混合动力汽车行驶工况预测研究";张昕 等;《汽车工程》;20141025(第10期);第1216-1236页 * |
"车辆道路行驶工况构建典型方法研究";王川;《赤峰学院学报(自然科学版)》;20160910(第17期);第35-36页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107346460A (zh) | 2017-11-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107346460B (zh) | 基于智能网联系统下前车运行信息的未来工况预测方法 | |
CN107284452B (zh) | 融合智能通讯信息的混合动力汽车未来工况预测系统 | |
CN107610464B (zh) | 一种基于高斯混合时间序列模型的轨迹预测方法 | |
CN109733390B (zh) | 一种基于驾驶人特性的自适应换道预警方法 | |
CN103903452B (zh) | 交通流短时预测方法 | |
CN109410587B (zh) | 一种城市快速路的宏观交通流参数估计方法 | |
CN113095558B (zh) | 一种智能网联汽车的迭代优化多尺度融合车速预测方法 | |
CN107909179A (zh) | 一种插电式混合动力车辆行驶工况的预测模型构建方法及车辆能量管理方法 | |
CN111806240A (zh) | 一种基于中国工况的电动汽车续驶里程预测方法 | |
CN110910531B (zh) | 一种基于车载obd信息的路面摩擦系数快速检测方法 | |
CN111582586B (zh) | 用于减少抖动的多车队驾驶危险预测系统与方法 | |
JP7504647B2 (ja) | 車両の走行曲線及び/又は走行経路を生成するための機械学習システムを訓練するための、コンピュータにより実行されるコンピュータ実装方法、車両の走行曲線及び/又は走行経路を生成するための方法、車両の走行曲線及び/又は走行経路を評価するための方法、車両の駆動システムの経路固有の排出量を特定するための方法、車両の駆動システムを適合させるための方法、並びに、コンピュータプログラム、機械可読記憶媒体、及び、コンピュータ実装式の機械学習システム | |
CN111951430A (zh) | 一种车辆驾驶性评价方法及其系统 | |
CN108831181A (zh) | 一种用于公交车辆行程时间预测的模型建立方法及系统 | |
CN111832743A (zh) | 用于确定行驶变化过程的方法 | |
CN115675098B (zh) | 基于用户特征差异的智能网联汽车动能回收系统、方法、设备及存储介质 | |
CN104599500A (zh) | 基于灰熵分析和改进贝叶斯融合的交通流预测方法 | |
CN113642768A (zh) | 一种基于工况重构的车辆行驶能耗预测方法 | |
CN106297296A (zh) | 一种基于稀疏轨迹点数据的细粒度旅行时间分配方法 | |
CN112489422B (zh) | 一种基于算法模型分析的车流量预测方法 | |
CN116151472A (zh) | 一种基于车联网数据挖掘的电动汽车短期能耗预测方法 | |
Feng et al. | Energy consumption prediction strategy for electric vehicle based on LSTM-transformer framework | |
CN114741806A (zh) | 一种悬架优化方法、系统、装置、设备及介质 | |
CN116340767B (zh) | 一种电动汽车行程能耗概率分布预测方法、系统及产品 | |
CN118152761B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180612 Termination date: 20200718 |
|
CF01 | Termination of patent right due to non-payment of annual fee |