CN106054170B - 一种约束条件下的机动目标跟踪方法 - Google Patents
一种约束条件下的机动目标跟踪方法 Download PDFInfo
- Publication number
- CN106054170B CN106054170B CN201610339347.8A CN201610339347A CN106054170B CN 106054170 B CN106054170 B CN 106054170B CN 201610339347 A CN201610339347 A CN 201610339347A CN 106054170 B CN106054170 B CN 106054170B
- Authority
- CN
- China
- Prior art keywords
- model
- filtering
- state
- calculating
- measurement
- 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
- 238000000034 method Methods 0.000 title claims abstract description 75
- 230000033001 locomotion Effects 0.000 claims abstract description 49
- 238000005259 measurement Methods 0.000 claims abstract description 48
- 230000003068 static effect Effects 0.000 claims abstract description 3
- 238000001914 filtration Methods 0.000 claims description 64
- 239000011159 matrix material Substances 0.000 claims description 23
- 230000004927 fusion Effects 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 14
- 239000013598 vector Substances 0.000 claims description 12
- 230000002452 interceptive effect Effects 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 7
- 230000007704 transition Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 4
- 230000028161 membrane depolarization Effects 0.000 claims description 2
- 238000007500 overflow downdraw method Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 description 8
- 230000001133 acceleration Effects 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- 238000000342 Monte Carlo simulation Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000002699 waste material Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000010485 coping Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 238000011423 initialization method Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
- G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
- G01S13/723—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种约束条件下的机动目标跟踪方法,所述方法包括以下步骤:A、获得机动目标的位置量测信息;B、将所述位置量测信息转换为笛卡尔坐标量测;C、在一维空间建立约束坐标静态模型,即建立机动目标在道路x,y方向上的位置关于里程数s的表达式;D、依据所述机动目标的运动特点设置运动模型,根据所建立的表达式建立状态方程和量测方程;E、分别针对各所述运动模型,对转换得到的笛卡尔坐标量测进行滤波,得到一维变量里程数s及其速度的滤波结果。根据本发明的实施方式,可以提高机动目标跟踪的性能。
Description
技术领域
本发明涉及约束条件下的机动目标跟踪,尤其涉及一种将二维约束转到一维空间下进行处理的机动目标跟踪方法。
背景技术
在雷达目标跟踪领域,存在一些场景,目标的运动轨迹发展不由其本身的运动速度决定,而是受到外界环境的严重影响甚至严格约束。例如,在公路上行驶的车辆受到地形和公路形状的约束,离开公路的可能性极小,这就是典型的约束目标跟踪问题。
然而,在处理约束目标跟踪问题的时候,常规目标跟踪方法存在明显的缺陷:一、约束轨迹形状复杂,传统运动模型失配严重;二、常规方法不考虑约束条件先验信息,造成信息浪费;三、跟踪结果难以满足约束条件,存在性能损失。因此,如何合理利用先验约束信息,建立准确的约束目标运动模型并提出对应的跟踪方法,对于提高约束目标跟踪精度具有重要的意义。
目前比较常见的约束目标跟踪方法包括:
一、模型降阶法,如W.Wen and H.F.Durrant-Whyte,“Model-based multi-sensordata fusion,”Proceedings of IEEE International Conference on Robotics andAutomation,Nice,France,pp.1720--1726,1992(基于模型的多传感器数据融合),提取***模型中状态分量之间的约束关系式,利用约束关系式将状态分量互相表示来减少状态向量的维数;R.J.Hewett,M.T.Health,M.D.Butala and F.Kamalabadi,“A robust nullspace method for linear equality constrained state estimation,”IEEETransactions on Signal Processing,vol.58,pp.3961--3971,2010(针对线性等式约束状态估计的鲁棒零空间方法),通过零空间分解获得降阶模型,这类方法简化了***并获得了简单模型约束状态的最优估计,但缺点是模型降阶后物理意义不清晰,且只适用于简单问题难以推广。
二、完美量测法,如L.S.Wang,Y.T.Chiang and F.R.Chang,“Filtering methodfor nonlinear systems with constraints,”IEEE Proceedings of Control Theoryand Applications,vol.149(6),pp.525--531,2002(包含约束条件的非线性***滤波方法),即将等式约束作为没有噪声的伪量测引入,把约束估计问题便转化为观测扩维的常规估计问题来解决。但由于观测噪声协方差矩阵是奇异矩阵,因此在卡尔曼滤波求解时可能引起数值计算问题。
三、估计后投影法,如D.Simon and T.L.Chia,“Kalman filtering with stateequality constraints,”IEEE Transactions on Aerospace and Electronic Systems,vol.38(1),pp.128--136,2002(状态等式约束下的卡尔曼滤波),利用约束优化技术将无约束估计向约束子空间投影,寻找约束子空间中距离无约束估计最近的一点作为最终估计。这样得到的估计结果必定满足约束条件,但不是约束状态的最优估计。
四、模型修正法,主要有两类方法,一种如S.Ko and R.Bitmead,“Stateestimation for linear systems with state equality constraints,”Automatica,vol.43(8),pp.1363-1368,2007(状态等式约束下的线性***状态估计),通过左乘投影矩阵实现无约束动态模型和约束条件的融合。另一种如Z.S.Duan,X.R.Li and J.F.Ru,“Design and Analysis of Linear Equality Constrained Dynamic Systems,”Proceedings of the 15th International Conference on Information Fusion,Singapore,pp.2537--2544,2012(线性等式约束动态***的设计与分析),采用直接消除法得到满足约束条件的修正模型。这类方法仍然无法解决常规模型与复杂约束之间的失配问题。
上述几种方法都是在常规无约束模型的基础上研究约束跟踪问题,并没有***的解决模型失配等问题,且大多是基于目标无机动的假设开展研究。然而实际场景中,目标可能会有变速机动。以道路场景为例,目标可能在转弯处减速,超车时加速等。因此,提出一种更加贴近物理实际的机动目标运动建模和状态估计方法具有重要的实践意义。
发明内容
为了克服上述现有技术存在的一项或更多的缺陷,本发明提供一种基于一维约束坐标的机动目标跟踪方法。
为实现上述目的,本发明提供了一种约束条件下的机动目标跟踪方法,该方法包括以下步骤:A、获得机动目标的位置量测信息;B、将所述位置量测信息转换为笛卡尔坐标量测;C、在一维空间建立约束坐标静态模型,即建立机动目标在道路x,y方向上的位置关于一维变量里程数s的表达式;D、依据所述机动目标的运动特点设置运动模型,根据所建立的表达式建立状态方程和量测方程;E、分别针对各所述运动模型,对转换得到的笛卡尔坐标量测进行滤波,得到一维变量里程数s及其速度的滤波结果。
根据本发明的一种实施方式,所述机动目标为道路上行驶的车辆,所述机动目标的位置量测信息为车辆相对观测雷达坐标系原点的距离和方位角,在所述步骤B中,利用去偏量测转换方法将所述位置量测信息转换为笛卡尔坐标量测。
根据本发明的一种实施方式,在所述步骤C中,建立机动目标在道路x,y方向上的位置关于一维变量里程数s的表达式:
C1:当道路段为直线段时,直线段上的任意点位置满足表达式:
其中,方向角(x1,y1)为直线段起始点,(x2,y2)为直线段终点;
C2:当道路段为圆弧段时,圆弧段上的任意点位置满足表达式:
其中,(x1,y1)为圆弧段起始点,(x2,y2)为圆弧段终点,(x0,y0)为圆心,圆弧起始角和曲率半径
根据本发明的一种实施方式,所述步骤E包括计算各运动模型对应的模型概率的步骤,以及按模型概率针对各所述运动的滤波结果进行加权融合的步骤。
根据本发明的一种实施方式,所述步骤E滤波包括如下的步骤:
E1:将k-1时刻滤波结果相互作用初值及其方差与转换量测值一起作为k时刻第j个模型的输入值,利用公式进行状态一步预测,其中,Φk-1为状态转移矩阵,
E2、计算状态一步预测协方差
其中Qk-1过程噪声协方差矩阵;
E3、计算无迹变换δ采样点状态预测
其中κ是尺度参数,其数值满足是均方根矩阵的第l行或第l列,n是状态向量的维数,对应对应
E4、计算δ采样点量测预测
其中hk为量测函数,由所述一维坐标约束建模给出;
E5、计算量测预测
其中是无迹变换时对均值加权计算使用的权值
E6、计算新息指量测值与量测预测值的差
E7、计算量测预测协方差即新息协方差
其中,是无迹变换时对均值加权计算使用的权值,为量测噪声协方差矩阵
E8、计算量测和状态向量的交互协方差
E9、确定滤波增益
E10、进行状态更新,获得k时刻j模型滤波结果
E11、进行协方差更新,获得k时刻j模型滤波误差协方差
E12、模型概率更新:
在得到各个模型k时刻的输出和后,更新该时刻各模型对应的概率计算对应模型j的可能性为
且对运动模型j的模型概率进行更新:
E13、k时刻交互式多模型滤波的融合输出为
根据本发明的一种实施方式,所述方法还包括基于k时刻针对各运动模型的滤波结果外推下一时刻预测结果,并根据k时刻模型概率对针对各运动模型的预测结果加权融合,得到k+1时刻融合预测结果的步骤。
本发明一些实施方式的一维坐标机动目标建模方法保证了过程噪声的自由度为1,即只在沿道路方向上变化,比起传统的在二维自由空间建模的方法更加贴近物理实际。
根据本发明的实施方式,能够提高跟踪性能。
附图说明
结合附图,可以更好地理解本发明,但是附图仅仅是示例性的,不是对本发明的限制。
图1示出了本发明的机动目标跟踪方法的一种实施方式的流程图;
图2示出了本发明的机动目标跟踪方法的滤波流程图;
图3示出了本发明的一个实施例中道路段参数示意图;
图4和图5示出了本发明线性约束仿真实验的200次蒙特卡洛仿真结果;
图6示出了本发明非线性约束仿真实验的单次轨迹仿真结果;
图7和图8示出了本发明非线性约束仿真实验的200次蒙特卡洛仿真结果。
具体实施方式
下面结合附图和实施例详细说明本发明的实施方式。附图和实施方式中的描述仅仅是示例性的,不是对本发明的保护范围的限制。
图1示出了依据本发明的一种实施方式的流程图。如图1所示,
首先在步骤S210,从观测雷达处获得雷达位置量测信息。在一种实施方式中,雷达位置量测信息为机动目标(该实施例以道路上行驶的车辆为例)相对观测雷达坐标系原点的距离和方位角。至于雷达如何获得这些信息、以及本发明实施方式的方法如何从雷达处获得这些信息,本领域技术人员可以采取各种方式来实现,无论采用哪种方法都在本发明的保护范围内。
接着,在步骤S220,将来自观测雷达的雷达位置量测信息(如上文所示,即机动目标相对观测雷达坐标系原点的距离和方位角)转换为笛卡尔坐标量测。可以采用多种方法进行这种转换,无论采用哪种方法都在本发明的范围内。在一种实施方式中,优选地,采用去偏量测转换方法。
去偏量测转换方法如下:
其中,代表转换量测向量,代表转换后x方向位置量测,代表转换后y方向位置量测,和是从雷达获取的距离和方位角量测,是基于目标量测位置得到的平均真实偏差。
对应的协方差矩阵为
其中
然后,在步骤S230,获得车辆的运动轨迹的约束条件,即一维坐标运动建模。这里的一维坐标是以目标所在道路为坐标轴,目标在道路上行驶的里程数为坐标。将道路近似为直线段和曲线段来考虑,参考图3。
对于直线段来说,必备的参数是起始点(x1,y1)和终点(x2,y2)。有了这两个点,道路直线段的长度S和方向角θ可以表示为:
被估计的状态变量是弧长或者称其为路程,包括它对应的速度和加速度,用s来表示。那么直线段上的任意点可以表示为:
对于圆弧段来说,除了起始点(x1,y1)和终点(x2,y2)这两个先验信息外,我们还需要圆心(x0,y0)来确定圆弧起始角θ0和曲率半径R,具体如下:
应该注意,在该式中,两个等号分别对应利用圆弧段起点和终点计算曲率半径的形式,不是递进关系而是并列关系。
同样以s作为变量,圆弧段上的任意点可以表示为:
由此得到目标x,y方向位置关于里程数s的表达式,该表达式可以用于量测方程构建和滤波结果转换中。
接着,在步骤S240,根据目标运动特点设置运动模型。针对不同应用场景采用不同运动模型来描述目标运动过程,运动模型主要由目标运动特性决定,影响状态方程的建立。在本发明的一些实施方式中,道路上行驶的车辆主要做近匀速或近匀加速运动,因此,运动模型分为近匀速模型NCV和近匀加速模型NCA。
之后,在步骤S250,进行模型滤波器初始化。
为进行滤波,需要建立状态方程和***观测方程。状态方程和***观测方程应该在滤波器初始化前进行,即建立方程后再初始化,然后进行滤波。
状态方程是目标运动规律的假设,根据选定的运动模型写出目标状态的递推表达式xk=Φkxk-1+Γkvk-1,其中Φk是状态转移矩阵,Γk是过程噪声分布矩阵,与运动模型有关。
或
分别对应近匀速NCV模型和近匀加速NCA模型,状态变量分别为和采用NCA模型时需要先将状态扩维,引入加速度分量,滤波后再降维,T为雷达扫描周期,vk-1为零均值高斯白噪声,其协方差为
其中qs是沿道路方向的过程噪声标准差;
***观测方程为
其中hk(xk)由所述S1中推导出的里程s与x,y坐标间关系给出,wk为去偏转换测量噪声。
尽管为了描述的方便,将模型滤波器的初始化设为单独的一个步骤,但是根据上下文,其也可以认为是后面滤波步骤的一部分。
尽管滤波初始化可以通过一维观测独立完成,但这样会造成信息浪费,初始化精度也不理想。根据本发明的一种实施方式,采用最小均方误差融合的方法将多维观测融合进行滤波初始化,既避免了信息浪费又提高了初始化精度。具体过程如下:
对于随机变量s=f(x,y),其方差与x,y方差之间的关系为:
其中和是随机变量x,y的方差。
若s可分别由x,y单独表示,即s=f1(x),s=f2(y),则可分别求出s对x,y的方差:
此时估计的融合公式为:
其中和是分别是采用x,y一维观测进行初始化的结果,一维观测初始化方法随轨迹形状不同而不同,以直线段道路轨迹为例,一维观测初始化方法为 为融合后的方差:
之后,在步骤S260进行滤波。
执行交互式多模型滤波,此时已知上一时刻各模型滤波器状态估计结果状态估计误差协方差上一时刻模型概率(上一时刻滤波结果),本时刻转换量测(雷达得到),过程噪声协方差矩阵Qk-1,量测噪声协方差矩阵状态转移矩阵Φk-1,量测函数hk(先验已知)。
首先是进行状态估计值之间的相互作用,目的是得到对应j模型滤波器的混合初值及其协方差
子滤波器模型Mj在k-1时刻的混合状态初值由模型Mi的输出相应的模型概率以及过渡概率ρij来计算,即
在k-1时刻,子滤波器模型Mj的状态误差协方差阵为
其中,为k-1时刻第i个子滤波器的状态估计,为相应的状态协方差矩阵,N为模型数目,
将状态向量及其方差与转换量测值一起作为k时刻第j个模型的输入值,由于状态变量与量测间存在非线性关系,选用无迹卡尔曼滤波器(UKF)进行滤波。
各模型滤波器单独滤波。由于状态向量和量测向量之间是非线性关系,需要非线性滤波方法来处理。可以采用多种非线性滤波方法,无论采用哪种方法都在本发明的范围内。在一种实施方式中,优选地,采用无迹卡尔曼滤波方法。该方法通过一系列精确选择的δ采样点经过非线性模型的映射来传递随机量的统计特性,是一种有效的非线性滤波方法。
单独滤波目的是得到对应j模型滤波器的滤波结果及其协方差滤波过程参考图2。
计算状态一步预测
其中,Φk-1为状态转移矩阵,为相互作用得到的混合初值。
计算协方差一步预测
其中,Qk-1为过程噪声协方差矩阵
计算无迹变换δ采样点状态预测
其中κ是尺度参数,可以为任何数值,只要(n+κ)≠0。是均方根矩阵的第l行或第l列,n是状态向量的维数,对应n=2;对应
计算δ采样点量测预测
其中hk为量测函数,由一维坐标约束建模给出
计算量测预测
其中是无迹变换时对均值加权计算使用的权值
计算新息指量测值与量测预测值的差
计算量测预测协方差即新息协方差
其中是无迹变换时对均值加权计算使用的权值,为量测噪声协方差矩阵
量测和状态向量的交互协方差
确定滤波增益
进行状态更新,获得k时刻j模型滤波结果
进行协方差更新,获得k时刻j模型滤波误差协方差
其中和分别是对均值和协方差加权计算时使用的权值
同样其中κ是尺度参数,可以为任何数值,只要(n+κ)≠0;n是状态向量的维数,对应对应a决定采样点在均值周围的散布,优选设定为比较小的正值,β用于引入变量的先验分布信息,优选高斯分布下设定β=2。
模型概率更新,目的是得到j模型对应的模型概率
在得到各个模型k时刻的输出和后,更新该时刻各模型对应的概率计算对应模型j的可能性为
模型j的模型更新为
状态估计融合,目的是得到融合滤波结果及其协方差Pk|k
k时刻交互式多模型滤波的融合输出为
于是便得到了对量测值的滤波结果及其协方差,接着对目标状态(包括位置和速度)进行外推预测。
最后,在步骤S270,将滤波结果转换到笛卡尔坐标系下
利用步骤S230中推导出的里程s与x,y坐标间关系将交互式多模型滤波结果转换到笛卡尔坐标系下,得到目标位置和速度滤波结果然后转第一步等待下一时刻k+1时刻雷达观测数据。
本发明所述的公路坐标运动建模方法和跟踪方法为约束目标跟踪提供了一条新途径,合理利用等式约束条件的同时有效避免了信息浪费。本发明的一些实施方式的公路坐标运动建模方法的一个优点在于保证了过程噪声的自由度为1,即只在沿道路方向上变化,比起传统的在二维自由空间建模的方法更加贴近物理实际。这种建模方法还使得目标动态模型与其运动轨迹独立,在保证滤波结果是满足约束条件的同时大大降低了复杂道路约束对滤波精度的影响,有效提高滤波性能。在跟踪方法中引入交互式多模型方法,提高了应对目标可能存在的变速机动的能力,多个模型同时作用有效匹配目标实际运动模型,减少模型失配导致的性能恶化。
为展示本发明的实施方式的效果,设置两组仿真场景分别对应线性等式约束和非线性等式约束。首先在线性等式约束场景中验证方法性能,跟踪一个在道路直线段上匀速行驶的目标。将本发明提出的约束坐标卡尔曼滤波方法(CCKF)与几种经典的约束估计算法进行比较,包括完美量测法(PMKF),估计后投影法(EPKF),线性等式约束卡尔曼滤波方法(LECKF)以及模型降维法(MRKF)。此外还选取一种经典的无约束算法,转换量测法(CMKF)作为对比。选取均方根误差(RMSE)作为性能的量度,进行200次蒙特卡洛仿真。
仿真结果参见图4、5所示,可以明显看出CCKF方法成功将约束条件引入跟踪***中,性能明显优于非约束的CMKF方法。与其他约束估计算法相比,CCKF也展现出一定的性能优势。这是由于其采用了不同的目标运动建模方法,运动模型更加贴近物理实际。
接着在非线性场景中验证交互式多模型约束坐标卡尔曼滤波方法(IMM-CCKF)跟踪机动目标的能力,跟踪一个在近似田径跑道形状道路段上变速运动的目标。将IMM-CCKF与包含匀速(CV)、匀加速(CA)和坐标转弯(CT)模型的交互式多模型算法(IMM-CVCACT、IMM-CVCA)对比。同样选取均方根误差(RMSE)作为性能的量度。单次轨迹仿真结果如图6所示,可以看出,滤波结果是满足约束条件的,证明IMM-CCKF可以有效引入约束条件并保证滤波结果满足约束条件。
200次蒙特卡洛仿真结果如图7、8所示,可以明显看出IMM-CCKF算法性能整体优于其他两种交互式多模型算法,尤其是在模型转换的时候,其他两种算法均出现了明显的滤波发散现象,而IMM-CCKF的曲线波动相对较为平缓,基本保持在相近的水平上。这是由于IMM-CCKF算法运动模型更加贴近物理实际,同时目标的动态模型与其运动轨迹独立,非线性约束对滤波过程的影响很小,因此得到了比较好的结果。
因而,本发明的一些实施方式将目标所在道路作为一维坐标系,于是道路上目标的运动便可以由其行驶里程数和速度随时间的变化来描述。再将整个公路分成多个道路段分别用直线段和圆弧段来近似,建立目标位置和行驶里程数之间的关系,在一维空间内对行驶里程数进行滤波。滤波跟踪过程利用无迹卡尔曼滤波(UKF)来处理一维状态变量和笛卡尔坐标量测之间存在的非线性关系,并结合交互式多模型(IMM)方法,引入近匀速(NCV)和近匀加速(NCA)模型应对目标存在机动的情况,得到一维滤波结果,最后将一维滤波结果转换到笛卡尔坐标系,得到目标位置的滤波结果。
另外,本发明所述的建模方法使得目标动态模型与其运动轨迹独立,在保证滤波结果是满足约束条件的同时大大降低了复杂道路约束对滤波精度的影响,有效提高滤波性能。
本发明所述的公路坐标机动目标跟踪方法中引入交互式多模型方法,提高了应对目标可能存在的变速机动的能力,多个模型同时作用有效匹配目标实际运动模型,减少模型失配导致的性能恶化。
应该注意,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种约束条件下的机动目标跟踪方法,其特征在于,所述方法包括以下步骤:
A、获得机动目标的位置量测信息;
B、将所述位置量测信息转换为笛卡尔坐标量测;
C、在一维空间建立约束坐标静态模型,即建立机动目标在道路x,y方向上的位置关于一维变量里程数s的表达式;
D、依据所述机动目标的运动特点设置运动模型,根据所建立的表达式建立状态方程和量测方程;该步骤中采用最小均方误差融合的方法将多维观测融合进行滤波初始化;
E、分别针对各所述运动模型,对转换得到的笛卡尔坐标量测进行滤波,得到一维变量里程数s及其速度的滤波结果;
在所述步骤C中,如下地建立机动目标在道路x,y方向上的位置关于里程数s的表达式:
C1:当道路段为直线段时,直线段上的任意点位置满足表达式:
其中,方向角(x1,y1)为直线段起始点,(x2,y2)为直线段终点;
C2:当道路段为圆弧段时,圆弧段上的任意点位置满足表达式:
其中,(x1,y1)为圆弧段起始点,(x2,y2)为圆弧段终点,(x0,y0)为圆心,圆弧起始角和曲率半径
所述步骤E包括计算各运动模型对应的模型概率的步骤,以及按模型概率针对各所述运动的滤波结果进行加权融合的步骤;
所述步骤E的滤波包括如下的步骤:
E1:将k-1时刻滤波结果相互作用初值及其方差与转换量测值一起作为k时刻第j个模型的输入值,利用公式进行状态一步预测,其中,Φk-1为状态转移矩阵,
E2、计算状态一步预测协方差
其中Qk-1过程噪声协方差矩阵;
E3、计算无迹变换δ采样点状态预测
其中κ是尺度参数,其数值满足是均方根矩阵的第l行或第l列,n是状态向量的维数,对应n=2;对应n=3;
E4、计算δ采样点量测预测
其中hk为量测函数,由所述一维坐标约束建模给出;
E5、计算量测预测
其中Wi m是无迹变换时对均值加权计算使用的权值
E6、计算新息指量测值与量测预测值的差
E7、计算量测预测协方差即新息协方差
其中,Wi c是无迹变换时对均值加权计算使用的权值,为量测噪声协方差矩阵
E8、计算量测和状态向量的交互协方差
E9、确定滤波增益
E10、进行状态更新,获得k时刻j模型滤波结果
E11、进行协方差更新,获得k时刻j模型滤波误差协方差
E12、模型概率更新:
在得到各个模型k时刻的输出和后,更新该时刻各模型对应的概率计算对应模型j的可能性为
且对运动模型j的模型概率进行更新:
E13、k时刻交互式多模型滤波的融合输出为
2.根据权利要求1所述的方法,其特征在于,所述机动目标为道路上行驶的车辆,所述机动目标的位置量测信息为车辆相对观测雷达坐标系原点的距离和方位角,在所述步骤B中,利用去偏量测转换方法将所述位置量测信息转换为笛卡尔坐标量测。
3.根据权利要求1所述的方法,其特征在于,还包括基于k时刻针对各运动模型的滤波结果外推下一时刻预测结果,并根据k时刻模型概率对针对各运动模型的预测结果加权融合,得到k+1时刻融合预测结果的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610339347.8A CN106054170B (zh) | 2016-05-19 | 2016-05-19 | 一种约束条件下的机动目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610339347.8A CN106054170B (zh) | 2016-05-19 | 2016-05-19 | 一种约束条件下的机动目标跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106054170A CN106054170A (zh) | 2016-10-26 |
CN106054170B true CN106054170B (zh) | 2017-07-25 |
Family
ID=57176647
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610339347.8A Active CN106054170B (zh) | 2016-05-19 | 2016-05-19 | 一种约束条件下的机动目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106054170B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109099927A (zh) * | 2018-09-26 | 2018-12-28 | 北京永安信通科技股份有限公司 | 对象定位方法、对象定位装置和电子设备 |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108020838B (zh) * | 2016-11-02 | 2021-08-20 | 惠州市德赛西威汽车电子股份有限公司 | 一种毫米波雷达信号在自适应巡航中的处理方法 |
CN106548475B (zh) * | 2016-11-18 | 2019-03-26 | 西北工业大学 | 一种适用于空间非合作自旋目标运动轨迹的预测方法 |
CN107045125B (zh) * | 2017-03-17 | 2020-01-14 | 电子科技大学 | 一种基于预测值量测转换的交互多模型雷达目标跟踪方法 |
CN108872975B (zh) * | 2017-05-15 | 2022-08-16 | 蔚来(安徽)控股有限公司 | 用于目标跟踪的车载毫米波雷达滤波估计方法、装置及存储介质 |
CN107387064A (zh) * | 2017-07-27 | 2017-11-24 | 河南科技学院 | 一种新的排爆机器人隧进定位方法 |
CN108037663A (zh) * | 2017-12-07 | 2018-05-15 | 电子科技大学 | 一种基于变维滤波算法的机动目标跟踪方法 |
CN108573271B (zh) * | 2017-12-15 | 2022-06-28 | 上海蔚来汽车有限公司 | 多传感器目标信息融合的优化方法及装置、计算机设备和记录介质 |
CN108572362B (zh) * | 2018-01-30 | 2020-07-31 | 哈尔滨工业大学 | 一种tws雷达空时联合关联跟踪方法及装置 |
CN109001699B (zh) * | 2018-01-30 | 2021-07-16 | 哈尔滨工业大学 | 基于带噪声目的地信息约束的跟踪方法 |
CN108507587B (zh) * | 2018-06-07 | 2021-09-17 | 南京邮电大学 | 一种基于坐标变换的三维车载定位导航方法 |
CN108871365B (zh) * | 2018-07-06 | 2020-10-20 | 哈尔滨工业大学 | 一种航向约束下的状态估计方法及*** |
CN108802721B (zh) * | 2018-08-22 | 2019-04-05 | 哈尔滨工业大学 | 一种任意直线约束下目标跟踪方法 |
CN109581302B (zh) * | 2018-12-12 | 2021-04-16 | 北京润科通用技术有限公司 | 一种车载雷达数据跟踪方法及*** |
CN109752002B (zh) * | 2018-12-26 | 2020-09-04 | 北京壹氢科技有限公司 | 一种多平台无源目标跟踪方法 |
CN110208792B (zh) * | 2019-06-26 | 2020-02-11 | 哈尔滨工业大学 | 同时估计目标状态和轨迹参数的任意直线约束跟踪方法 |
CN110501686A (zh) * | 2019-09-19 | 2019-11-26 | 哈尔滨工程大学 | 基于一种新型自适应高阶无迹卡尔曼滤波的状态估计方法 |
CN110927727B (zh) * | 2019-11-25 | 2021-09-03 | 深圳市智慧海洋科技有限公司 | 目标定位跟踪方法和装置 |
CN111798491B (zh) * | 2020-07-13 | 2022-09-06 | 哈尔滨工业大学 | 一种基于Elman神经网络的机动目标跟踪方法 |
CN112733907A (zh) * | 2020-12-31 | 2021-04-30 | 上海商汤临港智能科技有限公司 | 数据融合方法、装置、电子设备以及存储介质 |
CN114660546B (zh) * | 2022-05-23 | 2022-08-16 | 南京天朗防务科技有限公司 | 一维距离像目标真实尺寸的估计方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SG101444A1 (en) * | 2001-04-20 | 2004-01-30 | Singapore Tech Aerospace Ltd | A maneuvering target tracking method via modifying the interacting multiple model (imm) and the interacting acceleration compensation (iac) algorithms |
CN102568004A (zh) * | 2011-12-22 | 2012-07-11 | 南昌航空大学 | 一种高机动目标跟踪算法 |
-
2016
- 2016-05-19 CN CN201610339347.8A patent/CN106054170B/zh active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109099927A (zh) * | 2018-09-26 | 2018-12-28 | 北京永安信通科技股份有限公司 | 对象定位方法、对象定位装置和电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN106054170A (zh) | 2016-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106054170B (zh) | 一种约束条件下的机动目标跟踪方法 | |
CN109885883B (zh) | 一种基于gk聚类算法模型预测的无人车横向运动的控制方法 | |
CN111985093B (zh) | 一种带噪声估计器的自适应无迹卡尔曼滤波状态估计方法 | |
CN110849369B (zh) | 机器人跟踪方法、装置、设备及计算机可读存储介质 | |
CN106443661B (zh) | 基于无迹卡尔曼滤波的机动扩展目标跟踪方法 | |
CN110930495A (zh) | 基于多无人机协作的icp点云地图融合方法、***、装置及存储介质 | |
CN109284677B (zh) | 一种贝叶斯滤波目标跟踪算法 | |
CN111178385A (zh) | 一种鲁棒在线多传感器融合的目标跟踪方法 | |
CN105510882B (zh) | 基于目标机动参数估计的快速自适应采样周期跟踪方法 | |
CN108871365B (zh) | 一种航向约束下的状态估计方法及*** | |
CN111965618B (zh) | 一种融合多普勒量测的转换量测跟踪方法及*** | |
Zhao et al. | Adaptive non-linear joint probabilistic data association for vehicle target tracking | |
Li et al. | Joint localization based on split covariance intersection on the Lie group | |
Giefer et al. | Extended object tracking on the affine group aff (2) | |
Dahal et al. | Extended object tracking in curvilinear road coordinates for autonomous driving | |
Xu et al. | Dynamic vehicle pose estimation and tracking based on motion feedback for LiDARs | |
CN108507587B (zh) | 一种基于坐标变换的三维车载定位导航方法 | |
Huang et al. | Incremental non-Gaussian inference for SLAM using normalizing flows | |
CN116047495B (zh) | 一种用于三坐标雷达的状态变换融合滤波跟踪方法 | |
CN113030940A (zh) | 一种转弯机动下的多星凸型扩展目标跟踪方法 | |
CN112379350A (zh) | 智能车辆毫米波雷达多目标跟踪方法、装置及设备 | |
Zhu et al. | An adaptive interactive multiple-model algorithm based on end-to-end learning | |
Havangi | An adaptive particle filter based on PSO and fuzzy inference system for nonlinear state systems | |
CN114445456B (zh) | 基于部分模型的数据驱动智能机动目标跟踪方法及装置 | |
CN114815899A (zh) | 基于3d激光雷达传感器的无人机三维空间路径规划方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |