CN111055849A - Intersection intelligent driving method and system based on support vector machine - Google Patents
Intersection intelligent driving method and system based on support vector machine Download PDFInfo
- Publication number
- CN111055849A CN111055849A CN201811206747.7A CN201811206747A CN111055849A CN 111055849 A CN111055849 A CN 111055849A CN 201811206747 A CN201811206747 A CN 201811206747A CN 111055849 A CN111055849 A CN 111055849A
- Authority
- CN
- China
- Prior art keywords
- support vector
- vector machine
- sequence
- vehicle
- expanded
- 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
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
- B60W2050/0001—Details of the control system
- B60W2050/0019—Control system elements or transfer functions
- B60W2050/0028—Mathematical models, e.g. for simulation
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Automation & Control Theory (AREA)
- Human Computer Interaction (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Traffic Control Systems (AREA)
Abstract
Description
技术领域technical field
本发明是有关于一种路口智能驾驶方法及其系统,且尤其是有关一种基于支持向量机的路口智能驾驶方法及其系统。The present invention relates to an intelligent driving method at an intersection and a system thereof, and in particular to an intelligent driving method and system at an intersection based on a support vector machine.
背景技术Background technique
一般而言,十字路口或道路交会处具有多方向的车辆转弯或直行交会,因此在通过路口时,需要依驾驶的判断来进行加速、减速或定速行驶,一旦驾驶判断错误,交通事故便会发生。根据美国统计局的统计,2008年于十字路口或道路交会处发生交通事故的比例高达40%;而根据德国联邦统计局的统计,2013于十字路口或道路交会处发生交通事故的比例高达47.5%,在部分国家,发生交通事故的比例甚至高达98%。Generally speaking, there are multi-directional vehicle turns or straight intersections at intersections or road intersections. Therefore, when passing the intersection, it is necessary to accelerate, decelerate or drive at a constant speed according to the driving judgment. Once the driving judgment is wrong, traffic accidents will occur. occur. According to statistics from the US Bureau of Statistics, the proportion of traffic accidents at intersections or road junctions in 2008 was as high as 40%; according to the German Federal Statistics Office, the proportion of traffic accidents at intersections or road junctions in 2013 was as high as 47.5% , In some countries, the proportion of traffic accidents is even as high as 98%.
为了辅助驾驶通过路口时的决策判断,有业/学者发展出高度自动化车辆(highlyautomated vehicle,HAV),其包含人工智能方式进行机器学习,以辅助驾驶决策,而支持向量机便是其中一种机器学习方式,其透过建构模型进行预测或估计,便可进行决策。例如通过路口时,何时该加速、减速或定速。In order to assist in decision-making and judgment when driving through intersections, professionals/scholars have developed a highly automated vehicle (HAV), which includes artificial intelligence for machine learning to assist driving decision-making, and support vector machine is one of them. A learning approach that enables decision-making by constructing models to predict or estimate. For example, when to accelerate, decelerate, or stabilize when passing through an intersection.
在实际的情况下,驾驶会依据过去几个单位时间内的周遭信息,做出适合当下的驾驶行为,也就是说,真实的驾驶情况是具有时间相依关系的。然而,在习知的支持向量机的训练过程中,用来训练的数据虽然是依照连续的时间点来录制,但每一个时间所观察的数据皆被当作一笔独立的数据,而未考虑变数有时间相依性,因此判断出来的决策准确度仍有改善空间。In practical situations, driving will make appropriate driving behaviors based on the surrounding information in the past several units of time, that is to say, the real driving situation is time-dependent. However, in the training process of the conventional support vector machine, although the data used for training are recorded according to continuous time points, the data observed at each time is regarded as an independent data without considering The variables are time-dependent, so there is still room for improvement in the judgment accuracy.
有鉴于此,如何有效的提升支持向量机的决策准确度,遂成相关业者努力的目标。In view of this, how to effectively improve the decision-making accuracy of support vector machines has become the goal of the relevant industry.
发明内容SUMMARY OF THE INVENTION
本发明提供一种基于支持向量机的路口智能驾驶方法及其系统,其透过维度降低处理及时间补值处理后,可有效提升支持向量机的决策准确度。The present invention provides an intelligent driving method and system at an intersection based on a support vector machine, which can effectively improve the decision-making accuracy of the support vector machine through dimension reduction processing and time compensation processing.
依据本发明的一态样的一实施方式提供一种基于支持向量机的路口智能驾驶方法,其应用于一车辆且包含一支持向量机提供步骤、一数据处理步骤以及一决策步骤。于支持向量机提供步骤中,提供一支持向量机,所提供的支持向量机预先经过一训练过程,在训练过程中,提供一训练数据予支持向量机,训练数据由一原始数据经过一维度降低模块及一时间补值模块处理后获得,其中,原始数据包含多个训练样本,各训练样本包含通过一路口的一时间总值,及在时间总值内的多个取样时点中每一取样时点所对应的p个特征及一当下决策;维度降低模块将p个特征整合为k个新特征,时间补值模块提供一预设时间,时间补值模块将任一训练样本中任一取样时点及在前述任一取样时点之前的其他取样时点所分别对应的新特征视为一待扩展数列,当待扩展数列的长度小于预设时间内所具有的取样时点的数目时,于待扩展数列补入一预估值后形成一新待扩展数列,其中以待扩展数列的联合分配及所有新特征值的数据服从高斯分配求出预估值的条件分配,且时间补值模块将新待扩展数列重整为一扩展数列,扩展数列的长度等于预设时间内所具有的取样时点的数目,其中,p、k为正整数,且p>k;于数据处理步骤中,将一环境感测单元所获取的p个特征经由维度降低模块及时间补值模块处理后,提供予支持向量机进行分类;于决策步骤中,以支持向量机的分类结果决定车辆的驾驶行为。An embodiment according to an aspect of the present invention provides a support vector machine-based intelligent driving method at an intersection, which is applied to a vehicle and includes a support vector machine providing step, a data processing step, and a decision-making step. In the step of providing the support vector machine, a support vector machine is provided, and the provided support vector machine has undergone a training process in advance. During the training process, a training data is provided to the support vector machine, and the training data is reduced by one dimension from an original data module and a time compensation module are processed and obtained, wherein the original data includes a plurality of training samples, each training sample includes a time total value passing through an intersection, and each sample in a plurality of sampling time points within the time total value The p features corresponding to the time point and the next decision; the dimension reduction module integrates the p features into k new features, the time compensation module provides a preset time, and the time compensation module samples any one of the training samples. The time point and the new features corresponding to other sampling time points before any of the aforementioned sampling time points are regarded as a sequence to be expanded. When the length of the sequence to be expanded is less than the number of sampling time points in the preset time, A new sequence to be expanded is formed after an estimated value is added to the sequence to be expanded, wherein the conditional distribution of the predicted value is obtained by the joint distribution of the sequence to be expanded and the data of all new eigenvalues obeying Gaussian distribution, and the time compensation module Reshape the new sequence to be expanded into an expanded sequence, and the length of the expanded sequence is equal to the number of sampling time points in a preset time, wherein p and k are positive integers, and p>k; in the data processing step, The p features obtained by an environment sensing unit are processed by the dimension reduction module and the time compensation module, and then provided to the support vector machine for classification; in the decision-making step, the driving behavior of the vehicle is determined by the classification result of the support vector machine.
借此,训练数据及驾驶当下所获取的特征在经过维度降低模块及时间补值模块处理后,会具有时间相依关系,而能提升预测结果的准确度。In this way, after the training data and the features obtained while driving are processed by the dimension reduction module and the time compensation module, there will be a time dependency, which can improve the accuracy of the prediction results.
依据前述的基于支持向量机的路口智能驾驶方法的多个实施例,其中维度降低模块可采用一主成分分析法。或时间补值模块可采用一均匀缩放法。或预设时间可等于时间总值中的最大者。According to the foregoing embodiments of the support vector machine-based intelligent driving method at an intersection, the dimension reduction module may employ a principal component analysis method. Or the time compensation module can use a uniform scaling method. Or the preset time may be equal to the largest of the total time values.
依据前述的基于支持向量机的路口智能驾驶方法的多个实施例,其中p个特征可包含车辆相对一来车的一横向速度、车辆相对来车的一横向加速度、车辆相对来车的一纵向速度、车辆相对来车的一纵向加速度、车辆与来车的一距离、车辆与路口的一距离,及来车的一速度。原始数据中的多个特征可由环境感测单元取得,环境感测单元包含一雷达、一摄影机及一GPS定位装置中至少其中之一。According to various embodiments of the aforementioned SVM-based intelligent driving method at intersections, the p features may include a lateral speed of the vehicle relative to an oncoming vehicle, a lateral acceleration of the vehicle relative to an oncoming vehicle, and a longitudinal direction of the vehicle relative to an oncoming vehicle. Speed, a longitudinal acceleration of the vehicle relative to the oncoming vehicle, a distance between the vehicle and the oncoming vehicle, a distance between the vehicle and the intersection, and a speed of the oncoming vehicle. A plurality of features in the raw data can be obtained by an environment sensing unit, and the environment sensing unit includes at least one of a radar, a camera and a GPS positioning device.
依据本发明的一态样的另一实施方式提供一种基于支持向量机的路口智能驾驶方法,其应用于一车辆且包含一支持向量机提供步骤、一数据处理步骤及一决策步骤。于支持向量机提供步骤中,提供一支持向量机,所提供的支持向量机预先经过一训练过程,在训练过程中,提供一训练数据予支持向量机,训练数据由一原始数据经过一维度降低模块及一时间补值模块处理后获得,其中,原始数据包含多个训练样本,各训练样本包含通过一路口的一时间总值,及在时间总值内的多个取样时点中每一取样时点所对应的p个特征及一当下决策;维度降低模块将p个特征整合为k个新特征,时间补值模块提供一预设时间,时间补值模块将任一训练样本中任一取样时点及在前述任一取样时点之前的其他取样时点所分别对应的新特征视为一待扩展数列,且当待扩展数列的长度小于预设时间内所具有的取样时点的数目时,于待扩展数列补入一预估值后形成一新待扩展数列,且时间补值模块将新待扩展数列重整为一扩展数列,扩展数列的长度等于预设时间内所具有的取样时点的数目,其中,p、k为正整数,且p>k;于数据处理步骤中,将一环境感测单元所获取的p个特征经由维度降低模块及时间补值模块处理后,提供予支持向量机进行分类;于决策步骤中,以支持向量机的分类结果决定车辆的行为。Another embodiment according to an aspect of the present invention provides a support vector machine-based intelligent driving method at an intersection, which is applied to a vehicle and includes a support vector machine providing step, a data processing step, and a decision-making step. In the step of providing the support vector machine, a support vector machine is provided, and the provided support vector machine has undergone a training process in advance. During the training process, a training data is provided to the support vector machine, and the training data is reduced by one dimension from an original data module and a time compensation module are processed and obtained, wherein the original data includes a plurality of training samples, each training sample includes a time total value passing through an intersection, and each sample in a plurality of sampling time points within the time total value The p features corresponding to the time point and the next decision; the dimension reduction module integrates the p features into k new features, the time compensation module provides a preset time, and the time compensation module samples any one of the training samples. The time point and the new features corresponding to other sampling time points before any of the foregoing sampling time points are regarded as a sequence to be expanded, and when the length of the sequence to be expanded is less than the number of sampling time points in the preset time , a new sequence to be expanded is formed after the sequence to be expanded is filled with an estimated value, and the time compensation module reshapes the new sequence to be expanded into an expanded sequence, and the length of the expanded sequence is equal to the sampling time in the preset time. The number of points, wherein p and k are positive integers, and p>k; in the data processing step, the p features obtained by an environment sensing unit are processed by the dimension reduction module and the time compensation module, and then provided to The support vector machine is used for classification; in the decision-making step, the behavior of the vehicle is determined by the classification result of the support vector machine.
依据前述的基于支持向量机的路口智能驾驶方法的多个实施例,其中,时间补值模块可使用一主成分分析法。或时间补值模块可使用一均匀缩放法。According to various embodiments of the aforementioned SVM-based intelligent driving method at an intersection, the time compensation module can use a principal component analysis method. Or the time complement module can use a uniform scaling method.
依据本发明的另一态样的一实施方式提供一种基于支持向量机的路口智能驾驶系统,其应用于一车辆,基于支持向量机的路口智能驾驶系统包含一处理单元以及一环境感测单元;处理单元设置于车辆且包含一维度降低模块、一时间补值模块及一支持向量机。维度降低模块将多个取样时点中每一取样时点所对应的p个特征整合为k个新特征,其中,p、k为正整数,且p>k;时间补值模块提供一预设时间,时间补值模块将任一取样时点及在前述任一取样时点之前的其他取样时点所分别对应的新特征视为一待扩展数列,且当待扩展数列的长度小于预设时间内所具有的取样时点的数目时,于待扩展数列补入一预估值后形成一新待扩展数列,且时间补值模块将新待扩展数列重整为一扩展数列,扩展数列的长度等于预设时间内所具有的取样时点的数目;支持向量机经过一训练数据训练,训练数据由一原始数据经过维度降低模块及时间补值模块处理后获得,原始数据包含多个训练样本,各训练样本包含通过一路口的一时间总值,及在时间总值内的各取样时点所对应的p个特征及一当下决策。环境感测单元设置于车辆且信号连接处理单元,环境感测单元用以取得p个特征;其中,环境感测单元所获取的p个特征经由处理单元的维度降低模块及时间补值模块处理后,提供予支持向量机进行分类,支持向量机的分类结果用以决定车辆的驾驶行为。An embodiment according to another aspect of the present invention provides a support vector machine-based intersection intelligent driving system, which is applied to a vehicle. The support vector machine-based intersection intelligent driving system includes a processing unit and an environment sensing unit ; The processing unit is arranged in the vehicle and includes a dimension reduction module, a time compensation module and a support vector machine. The dimension reduction module integrates the p features corresponding to each sampling time point in the plurality of sampling time points into k new features, wherein p and k are positive integers, and p>k; the time compensation module provides a preset Time, the time complement value module regards any sampling time point and the new features corresponding to other sampling time points before any of the aforementioned sampling time points as a sequence to be expanded, and when the length of the sequence to be expanded is less than the preset time When there is the number of sampling time points in the to-be-expanded sequence, a new sequence to be expanded is formed after an estimated value is added to the sequence to be expanded, and the time complement value module reshapes the new sequence to be expanded into an expanded sequence, the length of the expanded sequence is equal to the number of sampling time points in the preset time; the support vector machine is trained by a training data, and the training data is obtained from a raw data after being processed by a dimension reduction module and a time compensation module, and the original data includes multiple training samples, Each training sample includes a time total for passing through an intersection, and p features and a current decision corresponding to each sampling time point within the time total. The environment sensing unit is disposed on the vehicle and is connected to the processing unit with a signal, and the environment sensing unit is used to obtain p features; wherein, the p features obtained by the environment sensing unit are processed by the dimension reduction module and the time compensation module of the processing unit. , provided to the support vector machine for classification, and the classification result of the support vector machine is used to determine the driving behavior of the vehicle.
依据前述的基于支持向量机的路口智能驾驶系统的多个实施例,其中p个特征包含:车辆相对一来车的一横向速度、车辆相对来车的一横向加速度、车辆相对来车的一纵向速度、车辆相对来车的一纵向加速度、车辆与来车的一距离、车辆与路口的一距离及来车的一速度。或环境感测单元可包含一雷达、一摄影机及一GPS定位装置中至少其中之一。或当下决策可包含加速、减速或定速中至少其中之一。According to various embodiments of the aforementioned SVM-based intersection intelligent driving system, the p features include: a lateral speed of the vehicle relative to an oncoming vehicle, a lateral acceleration of the vehicle relative to an oncoming vehicle, and a longitudinal direction of the vehicle relative to an oncoming vehicle. Speed, a longitudinal acceleration of the vehicle relative to the oncoming vehicle, a distance between the vehicle and the oncoming vehicle, a distance between the vehicle and the intersection and a speed of the oncoming vehicle. Or the environment sensing unit may include at least one of a radar, a camera and a GPS positioning device. Or the current decision may include at least one of acceleration, deceleration or constant speed.
附图说明Description of drawings
图1绘示依照本发明第一实施例的一种基于支持向量机的路口智能驾驶方法的流程图;FIG. 1 is a flowchart illustrating a method for intelligent driving at intersections based on support vector machines according to the first embodiment of the present invention;
图2绘示依照图1的基于支持向量机的路口智能驾驶方法的一第一模拟训练;2 illustrates a first simulation training of the SVM-based intelligent driving method at intersections according to FIG. 1;
图3绘示依照图1的基于支持向量机的路口智能驾驶方法的一第二模拟训练;3 illustrates a second simulation training of the SVM-based intelligent driving method at intersections according to FIG. 1;
图4绘示依照图1的基于支持向量机的路口智能驾驶方法的一第三模拟训练;FIG. 4 shows a third simulation training of the intelligent driving method based on the support vector machine at intersection according to FIG. 1;
图5绘示图2的第一模拟训练的第一累积率;FIG. 5 illustrates the first accumulation rate of the first simulated training of FIG. 2;
图6绘示图3的第二模拟训练的第一累积率;FIG. 6 illustrates the first accumulation rate of the second simulation training of FIG. 3;
图7绘示图4的第三模拟训练的第一累积率;以及FIG. 7 illustrates the first accumulation rate of the third simulation training of FIG. 4; and
图8绘示依照本发明第三实施例的一种基于支持向量机的路口智能驾驶系统的方块图。FIG. 8 is a block diagram illustrating an intelligent driving system at an intersection based on a support vector machine according to a third embodiment of the present invention.
具体实施方式Detailed ways
以下将参照附图说明本发明的实施例。为明确说明起见,许多实务上的细节将在以下叙述中一并说明。然而,阅读者应了解到,这些实务上的细节不应用以限制本发明。也就是说,在本发明部分实施例中,这些实务上的细节是非必要的。此外,为简化附图起见,一些习知惯用的结构与元件在附图中将以简单示意的方式绘示;并且重复的元件将可能使用相同的编号或类似的编号表示。Embodiments of the present invention will be described below with reference to the accompanying drawings. For the sake of clarity, many practical details are set forth in the following description. The reader should understand, however, that these practical details should not be used to limit the invention. That is, in some embodiments of the present invention, these practical details are unnecessary. Furthermore, for the purpose of simplifying the drawings, some well-known and conventional structures and elements will be shown in the drawings in a simplified and schematic manner; and repeated elements will possibly be designated by the same or similar numerals.
请参阅图1,其中图1绘示依照本发明第一实施例的一种基于支持向量机的路口智能驾驶方法100的流程图。基于支持向量机的路口智能驾驶方法100应用于一车辆且包含一支持向量机提供步骤110、一数据处理步骤120以及一决策步骤130。Please refer to FIG. 1 , wherein FIG. 1 shows a flowchart of a
于支持向量机提供步骤110中,提供一支持向量机,支持向量机预先经过一训练过程,在训练过程中,提供一训练数据予支持向量机,训练数据由一原始数据经过一维度降低模块及一时间补值模块处理后获得,其中,原始数据包含多个训练样本,各训练样本包含通过一路口的一时间总值,及在时间总值内的多个取样时点中每一取样时点所对应的p个特征及一当下决策,维度降低模块将p个特征整合为k个新特征,时间补值模块提供一预设时间,时间补值模块视任一训练样本中任一取样时点及在前述任一取样时点之前的其他取样时点所分别对应的新特征为一待扩展数列,且时间补值模块将待扩展数列重整为一扩展数列,扩展数列的长度等于预设时间内所具有的取样时点的数目,其中,p、k为正整数,且p>k。In the support vector
于数据处理步骤120中,将一环境感测单元所获取的p个特征经由维度降低模块及时间补值模块处理后,提供予支持向量机进行分类。In the
于决策步骤130中,以支持向量机的分类结果决定车辆的驾驶行为。In
借此,训练数据及驾驶当下所获取的特征在经过维度降低模块及时间补值模块处理后,会具有时间相依关系,而能提升预测结果的准确度。后面将详述基于支持向量机的路口智能驾驶方法100的细节。In this way, after the training data and the features obtained while driving are processed by the dimension reduction module and the time compensation module, there will be a time dependency, which can improve the accuracy of the prediction results. The details of the
支持向量机是一种监督式机器学习的分类器,其可以用来辅助判定车辆的行为。而在支持向量机提供步骤110中,所提供的支持向量机是经过训练过程训练,而能判定车辆在经过路口时的减速、加速或定速行为。A support vector machine is a supervised machine learning classifier that can be used to assist in determining vehicle behavior. In the support vector
在训练过程中,可以利用模拟方式来模拟车辆经过路口的状况,而形成多个训练样本。在第一实施例中,模拟的平台可以是Tass international公司开发的PreScan高级驾驶员辅助系统(Advanced Driver Assistance Systems;ADAS),其可以建置相关的路口信息以进行车辆经过路口的模拟。在其他实施例中,亦可以在实际的道路上取得多个训练样本,或用其他模拟软件,不以此为限。In the training process, a simulation method can be used to simulate the situation of the vehicle passing through the intersection, so as to form a plurality of training samples. In the first embodiment, the simulated platform may be the PreScan advanced driver assistance system (Advanced Driver Assistance Systems; ADAS) developed by Tass International, which can build relevant intersection information to simulate vehicles passing through the intersection. In other embodiments, a plurality of training samples may be obtained on the actual road, or other simulation software may be used, which is not limited thereto.
车辆经过路口一次,其所取得的数据可视为一个训练样本。也就是说,当车辆经过路口十次,可取得十个训练样本。每一个训练样本中,包含通过路口的时间总值,及在时间总值内的多个取样时点中每一取样时点所对应的p个特征及一当下决策。举例而言,假设在第1个训练样本中,通过路口的时间总值为2秒,而每间隔0.4秒取样一次,则会有5个取样时点,而每个取样时点均会搜集p个特征及一个当下决策,当下决策可以是加速、减速或定速,而p个特征可包含车辆相对一来车的一横向速度、车辆相对来车的一横向加速度、车辆相对来车的一纵向速度、车辆相对来车的一纵向加速度、车辆与来车的一距离、车辆与路口的一距离,及来车的一速度。When the vehicle passes through the intersection once, the data obtained can be regarded as a training sample. That is to say, when the vehicle passes through the intersection ten times, ten training samples can be obtained. Each training sample includes the total time of passing through the intersection, and the p features and a next decision corresponding to each sampling time point among a plurality of sampling time points within the time total value. For example, assuming that in the first training sample, the total time for passing the intersection is 2 seconds, and sampling is performed every 0.4 seconds, there will be 5 sampling time points, and each sampling time point will collect p features and a current decision, the current decision can be acceleration, deceleration or constant speed, and p features can include a lateral speed of the vehicle relative to an oncoming car, a lateral acceleration of the vehicle relative to the oncoming car, and a longitudinal direction of the vehicle relative to the oncoming car. Speed, a longitudinal acceleration of the vehicle relative to the oncoming vehicle, a distance between the vehicle and the oncoming vehicle, a distance between the vehicle and the intersection, and a speed of the oncoming vehicle.
第1个训练样本的数据可如表1所示。其中,单一个训练样本中,q个取样时点所取得的p个特征可形成一个原特征矩阵X,X=(x1,…,xp),其中xi=(xi1,..,xiq)T,而阅读者应该了解到,当具有n个训练样本时,就会具有n个原特征矩阵Xl对应不同的取样时点数目ql,n、q为正整数,l为1到n的正整数,i为1到p的正整数。n个训练样本中所有的当下决策可形成一个当下决策矩阵ZZ,ZZ=(z1,…,zn)。下文中,Tlw代表第l个训练样本中第w个取样时点,w为1到ql的正整数,xlwi代表在取样时点Tlw所获得的第i个特征,zlw代表在取样时点Tlw所获得的当下决策。因此,表1中的T11为第1个训练样本中的第1个取样时点,在第一实施例中即为0.4秒,T12为第1个训练样本中的第2个取样时点,在第一实施例中即为0.8秒,x111代表在第1个训练样本中第1个取样时点T11所获得的第1个特征,x122代表在第1个训练样本中第2个取样时点T12所获得的第2个特征,z13代表在第1个训练样本中第3个取样时点T13所获得的当下决策,以此类推,不再赘述。The data of the first training sample can be shown in Table 1. Among them, in a single training sample, the p features obtained at q sampling time points can form an original feature matrix X, X=(x 1 ,...,x p ), where x i =(x i1 ,..., x iq ) T , and the reader should know that when there are n training samples, there will be n original feature matrices X l corresponding to different sampling time points q l , where n and q are positive integers, and l is 1 A positive integer to n, where i is a positive integer from 1 to p. All the current decisions in the n training samples can form a current decision matrix ZZ, ZZ=(z 1 ,...,z n ). Hereinafter, T lw represents the w-th sampling time point in the l-th training sample, w is a positive integer from 1 to q l , x lwi represents the i-th feature obtained at the sampling time point T lw , and z lw represents the The current decision obtained at the sampling time point T lw . Therefore, T 11 in Table 1 is the first sampling time point in the first training sample, which is 0.4 seconds in the first embodiment, and T 12 is the second sampling time point in the first training sample , which is 0.8 seconds in the first embodiment, x 111 represents the first feature obtained at the first sampling time point T 11 in the first training sample, and x 122 represents the second feature in the first training sample The second feature obtained at the sampling time point T 12 , z 13 represents the current decision obtained at the third sampling time point T 13 in the first training sample, and so on, and will not be repeated.
表1、第1个训练样本Table 1. The first training sample
又,假设在第2个训练样本中,通过路口的时间总值为2.4秒,而每间隔0.4秒取样一次,共有6个取样时点。则第2个训练样本的数据可如表2所示。Also, suppose that in the second training sample, the total time for passing through the intersection is 2.4 seconds, and sampling is performed every 0.4 seconds, with a total of 6 sampling time points. Then the data of the second training sample can be shown in Table 2.
表2、第2个训练样本Table 2. The second training sample
假设只有2个训练样本,则原始数据即包含表1及表2的数据。Assuming that there are only 2 training samples, the original data includes the data in Table 1 and Table 2.
上述的原始数据会经过维度降低模块及时间补值模块处理后转为训练数据。而维度降低模块可使用主成分分析法(Principal Component Analysis;PCA)、偏最小平方回归法(Partial Least Squares Regression;PLSR)、多维标度法(MultidimensionalScaling;MDS)、投影寻踪法(Projection Pursuit method)、主成分回归法(PrincipalComponent Regression;PCR)、二次判别分析法(Quadratic Discriminant Analysis;QDA)、正规化判别分析法(Regularized Discriminant Analysis;RDA)及线性判别分析法(Linear Discriminant Analysis;LDA)等。较佳地,时间补值模块使用主成分分析法,主成分分析法的相关公式如式(1)、式(2)及式(3)所示。The above-mentioned raw data will be converted into training data after being processed by the dimension reduction module and the time compensation module. The dimension reduction module can use Principal Component Analysis (PCA), Partial Least Squares Regression (PLSR), Multidimensional Scaling (MDS), Projection Pursuit method ), Principal Component Regression (PCR), Quadratic Discriminant Analysis (QDA), Regularized Discriminant Analysis (RDA) and Linear Discriminant Analysis (LDA) Wait. Preferably, the time compensation module uses a principal component analysis method, and the relevant formulas of the principal component analysis method are as shown in formula (1), formula (2) and formula (3).
Y=aTX (1)。Y=a T X (1).
上式是基于一个训练样本中一个取样时点的公式,所以未包含代表训练样本数目的变数l及代表取样时点的变数w。其中,Y代表经整合后的新特征矩阵,其包含k个新特征,即Y=(y1,….,yk),其中yj表示第j个新特征,j为1到k的正整数。而当考虑n个训练样本及其所对应的取样时点数目ql时,Yl=(y11,….,ylk),a为系数矩阵,aji表示第i个特征xi所对应的系数。而第一实施例中,维度降低模块使各训练样本中的p个特征整合成1个新特征,也就是说,上述的k为1。因此,经维度降低模块处理后的第1个训练样本如表3所示,经维度降低模块处理后的第2个训练样本如表4所示。其中,yljw代表经重整后对应取样时点Tlw的第j个新特征,经维度降低模块处理后的数据为(Yl,zl)。The above formula is a formula based on one sampling time point in one training sample, so the variable l representing the number of training samples and the variable w representing the sampling time point are not included. Among them, Y represents the integrated new feature matrix, which contains k new features, namely Y=(y 1 ,....,y k ), where y j represents the jth new feature, and j is a positive value from 1 to k Integer. And when considering n training samples and their corresponding sampling time points q l , Y l =(y 11 ,....,y lk ), a is the coefficient matrix, and a ji represents the coefficient corresponding to the i-th feature x i . In the first embodiment, the dimension reduction module integrates p features in each training sample into a new feature, that is, the above k is 1. Therefore, the first training sample processed by the dimension reduction module is shown in Table 3, and the second training sample processed by the dimension reduction module is shown in Table 4. Wherein, y ljw represents the j-th new feature corresponding to the sampling time point T lw after reformation, and the data processed by the dimension reduction module is (Y l , z l ).
表3、经维度降低模块处理后的第1个训练样本Table 3. The first training sample processed by the dimension reduction module
表4、经维度降低模块处理后的第2个训练样本Table 4. The second training sample processed by the dimension reduction module
接着,上述数据会再经过时间补值模块处理以进行补值。由于监督式分类器输入的分类数据须为同长度且为数列数据,因此通过时间补值模块可将每一单位时间的累积数据拉成相同时间长度。Then, the above data will be processed by the time compensation module for compensation. Since the classification data input by the supervised classifier must be of the same length and sequence data, the accumulated data of each unit time can be pulled into the same time length through the time compensation module.
时间补值模块的补值方法可采用动态时间校正法(dynamic time warping;DTW)或均匀缩放法(Uniform scaling)。较佳地,时间补值模块采用均匀缩放法。The compensation method of the time compensation module may adopt a dynamic time warping (DTW) method or a uniform scaling method (Uniform scaling). Preferably, the time compensation module adopts a uniform scaling method.
在使用均匀缩放法时,时间补值模块可提供预设时间,其中,预设时间可等于时间总值中的最大者。也就是说,在第一实施例中,第1个训练样本的时间总值为2秒,第2个训练样本的时间总值为2.4秒,最大值为2.4,故预设时间可定为2.4秒,预设时间内所具有的取样时点的数目为6。When the uniform scaling method is used, the time compensation module can provide a preset time, wherein the preset time can be equal to the largest of the total time values. That is to say, in the first embodiment, the total time value of the first training sample is 2 seconds, the total time value of the second training sample is 2.4 seconds, and the maximum value is 2.4, so the preset time can be set to 2.4 seconds, the number of sampling time points in the preset time is 6.
在进行补值之前,时间补值模块视任一训练样本中任一取样时点及在前述任一取样时点之前的其他取样时点所分别对应的新特征为一待扩展数列。表5为待扩展数列表,其中Lljw表示待扩展数列,Lljw=(ylj1,…,yljw)。若经维度降低模块后使各训练样本中的p个特征整合成1个新特征,则表7中的j均为1。Before performing the complement value, the time complement value module regards the new features corresponding to any sampling time point in any training sample and other sampling time points before the aforementioned any sampling time point as a sequence to be expanded. Table 5 is a list of numbers to be expanded, wherein L ljw represents a sequence of numbers to be expanded, and L ljw =(y lj1 , . . . , y ljw ). If the p features in each training sample are integrated into a new feature after the dimension reduction module, then j in Table 7 is all 1.
表5、待扩展数列表Table 5. List of numbers to be expanded
更详细地说,在表5中,待扩展数Lljw包含第l个训练样本中第1个至第w个取样时点Tl1~Tlw对应的所有第j个新特征ylj1~yljw。举例而言,对扩展数列L111而言,w=1,j=1,因此扩展数列L111具有第1个取样时点T11的1个第1个新特征y111,且待扩展数列L111的长度为1,也就是其由一个数值组成。对待扩展数列L213而言,w=3,j=1,因此扩展数列L213具有第2个训练样本中第1个至第3个取样时点T21~T23所对应的3个第1个新特征y211~y213,且待扩展数列L213的长度为3,也就是其由三个数值组成,其他以此类推。因此,透过补值可使待扩展数列重整为一扩展数列,且扩展数列的长度等于预设时间内所具有的取样时点的数目。因为在第一实施例中,预设时间内所具有的取样时点的数目为6,故扩展数列会由6个数值组成。是以,经时间补值模块补值后,所有扩展数列均由6个数值组成,如表6所示,其中L* ljw表示扩展数列。若经维度降低模块后使各训练样本中的p个特征整合成1个新特征,则表7中的j均为1。In more detail, in Table 5, the number to be expanded L ljw includes all the j-th new features y lj1 ˜y ljw corresponding to the first to w-th sampling time points T l1 ˜T lw in the l-th training sample . For example, for the extended sequence L 111 , w=1, j=1, so the extended sequence L 111 has a first new feature y 111 at the first sampling time point T 11 , and the to-be-extended sequence L 111 has a length of 1, that is, it consists of a single value. For the to-be-expanded sequence L 213 , w=3, j=1, so the expanded sequence L 213 has three first to third sampling time points T 21 to T 23 in the second training sample. There are new features y 211 -y 213 , and the length of the sequence L 213 to be extended is 3, that is, it consists of three numerical values, and so on. Therefore, the to-be-expanded sequence can be reformed into an expanded sequence through the complementary value, and the length of the expanded sequence is equal to the number of sampling time points in the preset time. Since in the first embodiment, the number of sampling time points in the preset time is 6, the extended sequence consists of 6 values. Therefore, after the time compensation module complements the value, all the extended arrays are composed of 6 values, as shown in Table 6, where L * ljw represents the extended array. If the p features in each training sample are integrated into a new feature after the dimension reduction module, then j in Table 7 is all 1.
表6、扩展数列表Table 6. List of expansion numbers
表6中的扩展数列L* ljw=(L* ljw1,…,L* ljwr),其为待扩展数列Lljw经由式(4)及式(5)转换而来。The extended sequence L * ljw =(L * ljw1 , . . . , L * ljwr ) in Table 6 is the conversion of the to-be-extended sequence L1jw through equations (4) and (5).
其中r=1,...,qs,qs表示预设时间内所具有的取样时点的数目,qs小于等于ql,当预设时间等于时间总值中的最大者时,qs=max(ql)。表示地板公式,且Z为整数集,也就是说,将r×w/qs的结果无条件舍去仅保留整数。 where r = 1 , . s = max(q l ). represents the floor formula, and Z is a set of integers, that is to say, the result of r×w/q s is unconditionally discarded and only integers are retained.
举例而言,L113=(y111,y112,y113)扩展为L* 113=(L* 1131,L* 1132,L* 1133,L* 1134,L* 1135,L* 1136),而L* 1131=y111,L* 1132=y111(因为2x3/6=1,所以选择将位于L113第1个位置的y111补入L* 113的第二个位置),L* 1133=y111(因为3x3/6=1.5,无条件舍去仅保留整数为1,所以选择将位于L113第1个位置的y111补入L* 113的第三个位置),L* 1134=y112(因为4x3/6=2,所以选择将位于L113第2个位置的y112补入L* 113的第四个位置),L* 1135=Y112(因为5x3/6=2.5,无条件舍去仅保留整数为2,所以选择将位于L113第2个位置y112的补入L* 113的第五个位置),L* 1136=y113(因为6x3/6=3,所以选择将位于L113第3个位置y113的补入L* 113的第六个位置)。在此要特别说明的是,当上述的待扩展数列的长度大于预设时间内所具有的取样时点的数目,也就是说,待扩展数列的长度大于预计重整后的扩展数列的长度时,则仍可依式(4)及式(5)进行缩值,而达到本发明的目的。For example, L 113 = (y 111 , y 112 , y 113 ) expands to L * 113 = (L * 1131 , L * 1132 , L * 1133 , L * 1134 , L * 1135 , L * 1136 ), and L * 1131 = y 111 , L * 1132 = y 111 (because 2x3/6 = 1, so choose to fill y 111 at the first position of L 113 into the second position of L * 113 ), L * 1133 = y 111 (because 3x3/6=1.5, unconditional rounding only keeps the integer as 1, so choose to fill the y 111 at the first position of L 113 into the third position of L * 113 ), L * 1134 = y 112 (Because 4x3/6=2, so choose to fill y 112 at the second position of L 113 into the fourth position of L * 113 ), L * 1135 = Y 112 (because 5x3/6=2.5, unconditionally discarded Only the integer is kept as 2, so the selection will fill the fifth position of L * 113 with the 2nd position y 112 of L 113 ), L * 1136 = y 113 (since 6x3/6 = 3, the selection will be at L 113 3rd position y 113 's complement into L * 113 's sixth position). It should be noted here that, when the length of the above-mentioned sequence to be expanded is greater than the number of sampling time points in the preset time, that is, the length of the sequence to be expanded is greater than the length of the expanded sequence after expected reformation , the value can still be reduced according to formula (4) and formula (5), so as to achieve the purpose of the present invention.
原始数据经过维度降低模块处理后的数据为(Yl,zl),(Yl,zl)再经过时间补值模块处理后转为训练数据(L* l,zl),L* l=(L*l1,…,L* lk), 训练数据如表7所示,其中最后一列因其不需补值而可以维持原数值,但为了方便表示其与支持向量机的关系,仍以L* ljwr表示。若经维度降低模块后使各训练样本中的p个特征整合成1个新特征,则表7中的j均为1。The original data processed by the dimension reduction module is (Y l , z l ), and (Y l , z l ) are processed by the time compensation module and then converted into training data (L * l , z l ), L * l =(L * l1,...,L * lk ), The training data is shown in Table 7, in which the last column can maintain the original value because it does not need to be supplemented, but it is still represented by L * ljwr for the convenience of expressing its relationship with the support vector machine. If the p features in each training sample are integrated into a new feature after the dimension reduction module, then j in Table 7 is all 1.
表7、训练数据Table 7. Training data
上述的训练数据即可提供给支持向量机使用,以找出超平面,支持向量机的相关公式如式(6)至式(9)所示,其中支持向量机的开始目标如式(6)所示,式(6)代入式(7),并根据微积分及若且唯若原则将式(6)、(7)改写为式(8)所示。The above training data can be provided to the support vector machine to find out the hyperplane. The relevant formulas of the support vector machine are shown in equations (6) to (9), and the starting target of the support vector machine is shown in equation (6) As shown, formula (6) is substituted into formula (7), and formulas (6) and (7) are rewritten as formula (8) according to calculus and the if-and-only principle.
其中,C为惩罚系数(cost variable)且大于0,W为元素系数(entriesparameter)、ξl为松弛参数(slack variable)、b为截距项参数(intercept term)、αd、αe为拉格朗日乘数(lagrange multiplier)、Φ为径向基函数核(radial biasfunction),其将超平面拓展至非线性切割。L* l、L* d及L* e代表上述的扩展值L* ljwr,其为了简单示意而省略其他变数,d、e均为变数。Among them, C is the penalty coefficient (cost variable) and is greater than 0, W is the element coefficient (entries parameter), ξ l is the slack variable (slack variable), b is the intercept term parameter (intercept term), α d , α e are the tension parameters The lagrange multiplier and Φ are radial bias function kernels, which extend the hyperplane to nonlinear cuts. L * l , L * d , and L * e represent the above-mentioned extended value L * ljwr , other variables are omitted for simplicity, and d and e are both variables.
当支持向量机经过此训练数据后,即可找到超平面以辅助判断决策。After the support vector machine passes through this training data, it can find the hyperplane to assist in the decision-making.
于数据处理步骤120中,车辆在行经路口时,会由环境感测单元即时搜集p个特征,环境感测单元可包含雷达、摄影机、GPS定位装置等多个感测装置,感测装置能用来侦测距离、车速等,而能感测到p个特征,然感测装置的种类及数量不限于此。车辆在每一取样时点所搜集的p个特征均会进入维度降低模块及时间补值模块进行处理,处理方式如上所述。接着,于决策步骤130中,经处理的数据进入支持向量机,由于支持向量机已事先透过训练过程并找到超平面,因此,当即时取得的p个特征经处理后进入支持向量机,即可产生分类结果,并以分类结果决定车辆的驾驶行为,如减速、加速或定速。In the
请参阅图2、图3及图4,其中图2绘示依照图1的基于支持向量机的路口智能驾驶方法100的一第一模拟训练,图3绘示依照图1的基于支持向量机的路口智能驾驶方法100的一第二模拟训练,图4绘示依照图1的基于支持向量机的路口智能驾驶方法100的一第三模拟训练。在第一模拟训练、第二模拟训练及第三模拟训练中,车辆V1经过一T形路口,而在第一模拟训练时,横向道路R1上未有任何来车;在第二模拟训练时,横向道路R1的左侧有来车V2;在第三模拟训练时,横向道路R1的右侧有来车V2。Please refer to FIG. 2 , FIG. 3 and FIG. 4 , wherein FIG. 2 shows a first simulation training of the SVM-based
车辆V1的车速为每小时40公里,来车V2的车速界于每小时15公里至每小时40公里之间。所搜集的7个特征为车辆V1相对来车V2的一横向速度、车辆V1相对来车V2的一横向加速度、车辆V1相对来车V2的一纵向速度、车辆V1相对来车V2的一纵向加速度、车辆V1与来车V2的一距离、车辆V1与路口的一距离,及来车V2的一速度。且第一模拟训练、第二模拟训练及第三模拟训练各别包含20个训练样本,当下决策包含减速、定速及加速等行为,而第一模拟训练的预设时间为16.9秒,第二模拟训练的预设时间为28.8秒,第三模拟训练的预设时间为21.7秒。The speed of the vehicle V1 is 40 kilometers per hour, and the speed of the oncoming vehicle V2 is between 15 kilometers per hour and 40 kilometers per hour. The collected 7 features are a lateral velocity of the vehicle V1 relative to the oncoming vehicle V2, a lateral acceleration of the vehicle V1 relative to the oncoming vehicle V2, a longitudinal velocity of the vehicle V1 relative to the oncoming vehicle V2, and a longitudinal acceleration of the vehicle V1 relative to the oncoming vehicle V2. , a distance between the vehicle V1 and the oncoming vehicle V2, a distance between the vehicle V1 and the intersection, and a speed of the oncoming vehicle V2. In addition, the first simulation training, the second simulation training and the third simulation training each contain 20 training samples. The current decision-making includes deceleration, constant speed, and acceleration. The preset time of the first simulation training is 16.9 seconds, and the second simulation training The preset time for the simulation training is 28.8 seconds, and the preset time for the third simulation training is 21.7 seconds.
请参阅图5、图6及图7,其中图5绘示图2的第一模拟训练的第一累积率,图6绘示图3的第二模拟训练的第一累积率,图7绘示图4的第三模拟训练的第一累积率。第一模拟训练的平均第一累积率为0.9619,第二模拟训练的平均第一累积率为0.7588,第三模拟训练的平均第一累积率为0.8014,上述的第一累积率均在0.7以上,表示经维度降低处理后的信息已尽可能地解释原数据,而符合要求。Please refer to FIGS. 5 , 6 and 7 , wherein FIG. 5 shows the first accumulation rate of the first simulated training in FIG. 2 , FIG. 6 shows the first accumulation rate of the second simulated training in FIG. 3 , and FIG. 7 shows The first accumulation rate of the third simulation training of FIG. 4 . The average first accumulation rate of the first simulation training was 0.9619, the average first accumulation rate of the second simulation training was 0.7588, and the average first accumulation rate of the third simulation training was 0.8014. The above-mentioned first accumulation rates were all above 0.7. Indicates that the information after dimensionality reduction has explained the original data as much as possible and meets the requirements.
表8表示在与第一模拟训练相同的情况下,第一比较例与使用本案的基于支持向量机的路口智能驾驶方法100的决策结果准确度(AC)比较,由表8可看出,使用基于支持向量机的路口智能驾驶方法100的决策准确度较高。其中,第一比较例亦是使用支持向量机进行分类,但不同点在于,第一比较例的支持向量机只有经过原始数据训练。Table 8 shows the comparison between the decision result accuracy (AC) of the first comparative example and the
表8、第一模拟训练情境下的决策准确度比较Table 8. Comparison of decision accuracy in the first simulated training situation
表9表示在与第二模拟训练相同的情况下,第二比较例与使用本案的基于支持向量机的路口智能驾驶方法100的决策结果准确度(AC)比较,由表9可看出,使用基于支持向量机的路口智能驾驶方法100的决策准确度较高。其中,第二比较例亦是使用支持向量机进行分类,但不同点在于,第二比较例的支持向量机只有经过原始数据训练。Table 9 shows the comparison between the decision result accuracy (AC) of the second comparative example and the
表9、第二模拟训练情境下的决策准确度比较Table 9. Comparison of decision-making accuracy in the second simulated training situation
表10表示在与第三模拟训练相同的情况下,第三比较例与使用本案的基于支持向量机的路口智能驾驶方法100的决策结果准确度(AC)比较,由表10可看出,使用基于支持向量机的路口智能驾驶方法100的决策准确度较高。其中,第三比较例亦是使用支持向量机进行分类,但不同点在于,第三比较例的支持向量机只有经过原始数据训练。Table 10 shows the comparison of the decision result accuracy (AC) between the third comparative example and the
表10、第三模拟训练情境下的决策准确度比较Table 10. Comparison of decision-making accuracy in the third simulated training situation
在此要特别说明的是,上述的所有测试均为PreScan的模拟结果,然亦可以使用实际道路进行测试。It should be noted here that all the above-mentioned tests are the simulation results of PreScan, but can also be tested using actual roads.
在本发明的第二实施例中,是于支持向量机提供步骤中,让时间补值模块提供一预设时间,时间补值模块视任一训练样本中任一取样时点及在前述任一取样时点之前的其他取样时点所分别对应的新特征为一待扩展数列,且当待扩展数列的长度小于预设时间内所具有的取样时点的数目时,于待扩展数列补入前述取样时点的下一个取样时点的一预估值后形成一新待扩展数列,且时间补值模块将新待扩展数列重整为一扩展数列,扩展数列的长度等于预设时间内所具有的取样时点的数目。In the second embodiment of the present invention, in the step of providing the support vector machine, the time compensation module provides a preset time, and the time compensation module depends on any sampling time point in any training sample and at any one of the aforementioned The new features corresponding to other sampling time points before the sampling time point are a sequence to be expanded, and when the length of the sequence to be expanded is less than the number of sampling time points in the preset time, the sequence to be expanded is filled with the aforementioned A new sequence to be expanded is formed after an estimated value of the next sampling time point at the sampling time point, and the time complement value module reshapes the new sequence to be expanded into an expanded sequence, and the length of the expanded sequence is equal to the length of the sequence in the preset time. the number of sampling time points.
更详细地说,假设原始数据包含表1及表2的数据,预设时间为2.4秒,待扩展数列如表5所示,因为表5中累积取样时点0~T11时对应的待扩展数列长度为1,其小于预设时间内所具有的取样时点的数目(等于6),且在此累积取样时点0~T11中尚未取得下一个取样时点T12的数据,因此可以补入一预估值y'1j2对应取样时点T12。类似地,在表5中,所有累积取样时点对应的待扩展数列长度均小于6,因此需分别补入预估值,则新待扩展数列L'ljw表如表11所示,其中y'ljw表示预估值。在此需特别说明的是,在表11的新待扩展数列中,累积取样时点0~T11、0~T25对应的待扩展数列长度为6,而不需再进行补值,但为了清楚的表现出新待扩展数列与待扩展数列的差异,仍将其列出,而不以此限制本发明。第二实施例中,将新待扩展数列重整为待扩展数列以及后面的支持向量机训练方式均与第一实施例相同,而不再赘述。In more detail, assuming that the original data includes the data in Table 1 and Table 2, the preset time is 2.4 seconds, and the number sequence to be expanded is shown in Table 5, because the accumulated
表11、新待扩展数列表Table 11. List of new to-be-expanded numbers
第二实施例中,使用待扩展数列的联合分配及所有新特征yljw的数据服从高斯分配可求出预估值的条件分配,最终可求得预估值。由于上述经维度降低后的新特征yljw可配适数据服从高斯分配(或称高斯随机过程,gaussian random process,marginaldistribution),即yljw~GP(μjw,Σjw),其中(μjw,∑jw)可由式(10)及式(11)估算。此外,可以由式(12)求得预估值前的其他新特征yljw(相当于待扩展数列)的联合分配(jointdistribution),其中式(12)中的w大于2,式(14)中的c指第c个训练样本,m指第m个训练样本。In the second embodiment, the conditional assignment of the estimated value can be obtained by using the joint assignment of the sequence to be expanded and the data of all new features y ljw obeying the Gaussian assignment, and finally the estimated value can be obtained. Since the above-mentioned new features y ljw after dimensionality reduction can be adapted to the data subject to Gaussian distribution (or Gaussian random process, marginal distribution), that is, y ljw ~GP(μ jw , Σ jw ), where (μ jw , Σ jw ) can be estimated from equations (10) and (11). In addition, the joint distribution (jointdistribution) of other new features y ljw (equivalent to the sequence to be expanded) before the estimated value can be obtained from equation (12), where w in equation (12) is greater than 2, and in equation (14) where c refers to the c-th training sample, and m refers to the m-th training sample.
ylj[1:(w-1)]≡(ylj1,…,ylj(w-1))~GP(μj[1:(w-1)],Σj[1:(w-1)]) (12)。y lj [1:(w-1)]≡(y lj1 ,...,y lj(w-1) )~GP(μ j[1:(w-1)] ,Σj [1:(w-1 ) )] ) (12).
(μj[1:(w-1)],∑j[1:(w-1)])可由式(13)、式(14)及式(15)估算。(μ j[1:(w-1)] , Σ j[1:(w-1)] ) can be estimated by Equation (13), Equation (14) and Equation (15).
最终,可求得在过去时间点下预测下一单位时间的条件分配(conditionaldistribution),也就是说,在累积取样时点0~T11预测取样时点T12的条件分配,条件分配如式(16)、式(17)所示。Finally, the conditional distribution of predicting the next unit time at the past time point can be obtained, that is, the conditional distribution of the predicted sampling time point T12 at the cumulative
而可由估算,有了条件分配,可以根据条件分配预测下一取样时点的预估值(其可根据Monte Carlo概念生成预估值)。and by Estimation, with conditional assignments, an estimate of the next sampling point in time can be predicted based on the conditional assignment (which can generate an estimate according to the Monte Carlo concept).
其中,为yljw,ylj[1:w-1]的变异数矩阵(covariance matrix),GP(·,·)为高斯随机过程的缩写,μ为高斯随机过程平均数,Σ为高斯随机过程变异数矩阵(Covariance-variance矩阵),为变异数矩阵内的元素(Element)估计量。in, is y ljw , the covariance matrix of y lj[1:w-1] , GP(·,·) is the abbreviation of Gaussian random process, μ is the mean of Gaussian random process, and Σ is the variance of Gaussian random process matrix (Covariance-variance matrix), is an estimator of the elements in the variance matrix.
表12表示在加入预估值后的第一模拟训练情境下的决策结果准确度,表13表示在加入预估值后的第二模拟训练情境下的决策结果准确度,表14表示在加入预估值后的第三模拟训练情境下的决策结果准确度。由表12、表13及表14的结果可知,加入预估值可使决策准确度更提升。Table 12 shows the accuracy of decision-making results in the first simulated training situation after adding the estimated value, Table 13 shows the accuracy of the decision-making result in the second simulated training situation after adding the estimated value, and Table 14 shows the accuracy of the decision-making result in the second simulated training situation after adding the estimated value. The accuracy of the decision results in the third simulated training situation after evaluation. From the results in Table 12, Table 13 and Table 14, it can be seen that adding the estimated value can improve the decision accuracy.
表12、第一模拟训练情境下的决策准确度Table 12. Decision Accuracy in the First Simulation Training Scenario
表13、第二模拟训练情境下的决策准确度Table 13. Decision Accuracy in the Second Simulation Training Scenario
表14、第三模拟训练情境下的决策准确度Table 14. Decision Accuracy in the Third Simulation Training Scenario
请参阅图8,其中图8绘示依照本发明第三实施例的一种基于支持向量机的路口智能驾驶系统200的方块图。基于支持向量机的路口智能驾驶系统200应用于一车辆且包含一处理单元210以及一环境感测单元220,处理单元包含一维度降低模块211、一时间补值模块212及一支持向量机213,维度降低模块211将多个取样时点中每一取样时点所对应的p个特征整合为k个新特征,其中,p、k为正整数,且p>k;时间补值模块212提供一预设时间,时间补值模块212将任一取样时点及在前述任一取样时点之前的其他取样时点所分别对应的新特征视为一待扩展数列,且当待扩展数列的长度小于预设时间内所具有的取样时点的数目时,于待扩展数列补入一预估值后形成一新待扩展数列,且时间补值模块212将新待扩展数列重整为一扩展数列,扩展数列的长度等于预设时间内所具有的取样时点的数目;支持向量机213经过一训练数据训练,训练数据由一原始数据经过维度降低模块211及时间补值模块212处理后获得,原始数据包含多个训练样本,各训练样本包含通过一路口的一时间总值,及在时间总值内的各取样时点所对应的p个特征及一当下决策。Please refer to FIG. 8 , wherein FIG. 8 is a block diagram of an
环境感测单元220设置于车辆且信号连接处理单元210,环境感测单元220用以取得p个特征;其中,环境感测单元220所获取的p个特征经由处理单元210的维度降低模块211及时间补值模块212处理后,提供予支持向量机213进行分类,支持向量机213的分类结果用以决定车辆的驾驶行为。The
借此,可以辅助车辆决定经过路口时的加速、减速或定速行为。维度降低模块211、时间补值模块212的处理细节及其与支持向量机213之间的关系如上面所述,在此不再赘述。而p个特征可包含车辆相对一来车的一横向速度、车辆相对来车的一横向加速度、车辆相对来车的一纵向速度、车辆相对来车的一纵向加速度、车辆与来车的一距离、车辆与路口的一距离及来车的一速度。环境感测单元220可包含一雷达、一摄影机及一GPS定位装置中至少其中之一。或当下决策可包含加速、减速或定速中至少其中之一。In this way, the vehicle can be assisted in determining the acceleration, deceleration or constant speed behavior when passing through the intersection. The processing details of the
虽然本发明已以实施例揭露如上,然其并非用以限定本发明,任何熟悉此技艺者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰,因此本发明的保护范围当视所附的权利要求书所界定的范围为准。Although the present invention has been disclosed by the above embodiments, it is not intended to limit the present invention. Anyone skilled in the art can make various changes and modifications without departing from the spirit and scope of the present invention. Therefore, the present invention The scope of protection shall be subject to the scope defined by the appended claims.
Claims (13)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811206747.7A CN111055849B (en) | 2018-10-17 | 2018-10-17 | Intersection intelligent driving method and system based on support vector machine |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811206747.7A CN111055849B (en) | 2018-10-17 | 2018-10-17 | Intersection intelligent driving method and system based on support vector machine |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111055849A true CN111055849A (en) | 2020-04-24 |
| CN111055849B CN111055849B (en) | 2021-04-06 |
Family
ID=70296833
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811206747.7A Active CN111055849B (en) | 2018-10-17 | 2018-10-17 | Intersection intelligent driving method and system based on support vector machine |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111055849B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI765720B (en) * | 2021-05-24 | 2022-05-21 | 中國鋼鐵股份有限公司 | Method and system for detecting a temperature of conveyor system |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105260712A (en) * | 2015-10-03 | 2016-01-20 | 上海大学 | Method and system for detecting pedestrians in front of a vehicle |
| CN104494600B (en) * | 2014-12-16 | 2016-11-02 | 电子科技大学 | A Driver Intention Recognition Method Based on SVM Algorithm |
| CN106650644A (en) * | 2016-12-07 | 2017-05-10 | 上海交通大学 | Identification method and system for dangerous behaviors of driver |
| CN107284452A (en) * | 2017-07-18 | 2017-10-24 | 吉林大学 | Merge the following operating mode forecasting system of hybrid vehicle of intelligent communication information |
| CN107491720A (en) * | 2017-04-01 | 2017-12-19 | 江苏移动信息系统集成有限公司 | A kind of model recognizing method based on modified convolutional neural networks |
| JP2018063708A (en) * | 2016-10-11 | 2018-04-19 | ピーエルケー テクノロジーズ カンパニー リミテッドPLK Technologies Co., Ltd. | Moving object collision warning device and method for large vehicle |
-
2018
- 2018-10-17 CN CN201811206747.7A patent/CN111055849B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104494600B (en) * | 2014-12-16 | 2016-11-02 | 电子科技大学 | A Driver Intention Recognition Method Based on SVM Algorithm |
| CN105260712A (en) * | 2015-10-03 | 2016-01-20 | 上海大学 | Method and system for detecting pedestrians in front of a vehicle |
| JP2018063708A (en) * | 2016-10-11 | 2018-04-19 | ピーエルケー テクノロジーズ カンパニー リミテッドPLK Technologies Co., Ltd. | Moving object collision warning device and method for large vehicle |
| CN106650644A (en) * | 2016-12-07 | 2017-05-10 | 上海交通大学 | Identification method and system for dangerous behaviors of driver |
| CN107491720A (en) * | 2017-04-01 | 2017-12-19 | 江苏移动信息系统集成有限公司 | A kind of model recognizing method based on modified convolutional neural networks |
| CN107284452A (en) * | 2017-07-18 | 2017-10-24 | 吉林大学 | Merge the following operating mode forecasting system of hybrid vehicle of intelligent communication information |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI765720B (en) * | 2021-05-24 | 2022-05-21 | 中國鋼鐵股份有限公司 | Method and system for detecting a temperature of conveyor system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111055849B (en) | 2021-04-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111886603B (en) | Neural networks for object detection and representation | |
| Aoude et al. | Behavior classification algorithms at intersections and validation using naturalistic data | |
| Oh et al. | Estimation of rear-end crash potential using vehicle trajectory data | |
| EP3940487B1 (en) | Estimation of probability of collision with increasing severity level for autonomous vehicles | |
| US12049230B2 (en) | Method and apparatus for determining information related to a lane change of target vehicle, method and apparatus for determining a vehicle comfort metric for a prediction of a driving maneuver of a target vehicle and computer program | |
| Lefèvre et al. | Comparison of parametric and non-parametric approaches for vehicle speed prediction | |
| US9129519B2 (en) | System and method for providing driver behavior classification at intersections and validation on large naturalistic data sets | |
| Sander et al. | The potential of clustering methods to define intersection test scenarios: Assessing real-life performance of AEB | |
| CN110188797B (en) | Intelligent automobile rapid test method based on Bayesian optimization | |
| US11577750B2 (en) | Method and apparatus for determining a vehicle comfort metric for a prediction of a driving maneuver of a target vehicle | |
| JP6045846B2 (en) | Traffic accident occurrence prediction device, method and program | |
| CN115195713A (en) | Method for determining the trajectory of an at least partially assisted motor vehicle | |
| TWI690440B (en) | Intelligent driving method for passing intersections based on support vector machine and intelligent driving system thereof | |
| Hu et al. | A framework for probabilistic generic traffic scene prediction | |
| Hyeon et al. | Short-term speed forecasting using vehicle wireless communications | |
| CN113424209B (en) | Trajectory prediction using deep learning multi-predictor fusion and Bayesian optimization | |
| CN114730494A (en) | Method for estimating the coverage of a space of a traffic scene | |
| Oh et al. | In-depth understanding of lane changing interactions for in-vehicle driving assistance systems | |
| EP4264436B1 (en) | Generating unknown-unsafe scenarios, improving automated vehicles, computer system | |
| Platho et al. | Predicting velocity profiles of road users at intersections using configurations | |
| CN111055849B (en) | Intersection intelligent driving method and system based on support vector machine | |
| Jia et al. | Connected multi-vehicle crash risk assessment considering probability and intensity | |
| CN115257801B (en) | Track planning method and device, server and computer readable storage medium | |
| US20240403630A1 (en) | Adaptive driving style | |
| US20230195977A1 (en) | Method and system for classifying scenarios of a virtual test, and training method |
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 |