CN107561503A - 一种基于多重渐消因子的自适应目标跟踪滤波方法 - Google Patents
一种基于多重渐消因子的自适应目标跟踪滤波方法 Download PDFInfo
- Publication number
- CN107561503A CN107561503A CN201710752493.8A CN201710752493A CN107561503A CN 107561503 A CN107561503 A CN 107561503A CN 201710752493 A CN201710752493 A CN 201710752493A CN 107561503 A CN107561503 A CN 107561503A
- Authority
- CN
- China
- Prior art keywords
- target
- measurement
- moment
- target tracking
- covariance
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于多重渐消因子的自适应目标跟踪滤波方法,属于目标跟踪领域,具体涉及一种自适应目标跟踪滤波方法。本发明包括:步骤一,建立运动目标的状态模型和量测模型;步骤二,初始化运动状态向量和状态误差协方差矩阵,计算目标运动状态向量预测值、量测向量预测值、新息、新息协方差和新息协方差的估计值;步骤三,计算指数加权因子和多重渐消因子;步骤四,计算状态预测协方差矩阵、卡尔曼增益、滤波值和滤波协方差;步骤五,重复执行步骤二~步骤四,直到目标跟踪结束。本发明渐消因子的计算简单,可以在系统模型未知或噪声统计信息不准确的情况下防止滤波器发散,提高了目标跟踪精度。本发明可运用于雷达目标跟踪。
Description
技术领域
本发明涉及一种自适应目标跟踪滤波方法,属于目标跟踪领域。
背景技术
在目标跟踪领域,卡尔曼滤波是应用最为广泛的跟踪滤波算法。卡尔曼滤波器通过对运动目标运动过程和传感器测量过程进行建模,获得在最小均方误差意义下对目标运动状态的最优估计。然而,卡尔曼滤波器只有在其数学模型确知的条件下才是最佳的。卡尔曼滤波器数学模型的建立包括建立状态方程、量测方程和确定初始状态估计、初始协方差、确定过程噪声和量测噪声的统计特性。但是在实际应用中,要精确建立滤波模型是相当困难的,而且运动目标的噪声统计特性也有可能随时间发生变化。在滤波模型与实际过程的数学模型不匹配时,会使滤波器精度下降,严重时会导致滤波器发散。为了防止滤波器发散,必须对卡尔曼滤波器进行调整。
因此,将卡尔曼滤波器应用于实际问题时,主要的工作就是建立滤波数学模型和寻求适用的调整算法。目前常用的方法有衰减记忆滤波、限定记忆滤波、自适应滤波以及上述方法的一系列组合。其中,自适应渐消记忆滤波器采用渐消因子来抑制滤波器的记忆长度,以便充分利用现时的观测数据,减小陈旧量测值的影响,因此渐消滤波的关键在于渐消因子的确定。然而,目前的算法中渐消因子的自适应估计公式较为复杂,而且算法通常引入的渐消因子是标量的单渐消因子,P(k+1|k)的各个通道仅具有相同的调节能力,目标跟踪精度较差。另外,在少数采用多重渐消因子的方案中,只将P(k+1|k)中与量测向量相关的通道乘以渐消因子,其他通道简单的置1不做调整,而没有充分利用目标运动状态向量各个方向元素之间的关系,即P(k+1|k)各个通道间的相关性,无法同时提高目标在各个维度的跟踪精度。
发明内容
本发明为解决现有技术存在的渐消因子的自适应估计较为复杂、目标跟踪精度较差的问题,提供了一种基于多重渐消因子的自适应目标跟踪滤波方法。
本发明所述一种基于多重渐消因子的自适应目标跟踪滤波方法,通过以下技术方案实现:
步骤一,建立运动目标的状态模型和量测模型,将状态模型的运动状态标记为N维的运动状态向量X,量测模型的量测值标记为M维的量测向量Z,根据运动模型获得状态误差协方差矩阵P和过程噪声协方差矩阵Q,根据量测模型获得量测噪声协方差矩阵R;
步骤二,获得第k时刻的运动状态向量X(k)、状态误差协方差矩阵P(k),按照卡尔曼递推公式计算运动目标在第k+1时刻的运动状态向量预测值量测向量预测值并计算得到第k+1时刻的目标的新息协方差S(k+1);根据第k+1时刻目标的量测向量Z(k+1)计算第k+1时刻新息v(k+1),并由v(k+1)得到第k+1时刻新息协方差的估计值
步骤三,计算第k+1时刻的指数加权因子根据和S(k+1)计算得到第k+1时刻多重渐消因子λ(k+1);
步骤四,通过第k+1时刻多重渐消因子λ(k+1)计算第k+1时刻的状态预测协方差矩阵P(k+1|k),并得到第k+1时刻的卡尔曼增益K(k+1),最终得到目标在第k+1时刻的滤波值和滤波协方差P(k+1|k+1);
步骤五,使k=k+1,判断目标跟踪是否结束,如果没有结束那么返回执行步骤二~步骤四,如果结束则完成所述目标跟踪滤波方法。
本发明与现有技术相比较,最为突出的特点和显著的有益效果是:
第一,根据目标运动状态向量和观测向量的关系,利用新息计算得到改进的多重渐消因子,渐消因子的计算简单,适应性强;
第二,引入可调节的指数加权因子,调整权重,减少陈旧量测值的影响,防止滤波器发散;
第三,根据目标运动状态向量各个方向元素之间的关系,即状态预测协方差矩阵P(k+1|k)各个通道间的相关性,引入改进的多重渐消因子,对P(k+1|k)的各个通道分别进行加权修正,因此目标跟踪精度较高。
选取均方根误差(RMSE)作为性能的量度,进行100次蒙特卡罗仿真实验;图2是本发明和匀速运动(CV)模型的卡尔曼滤波方法的跟踪位置对比图,由于模型失配,卡尔曼滤波算法的滤波轨迹逐渐偏离了目标的实际位置,而本发明的滤波轨迹结果始终没有发散,一直与目标运动趋势一致,目标跟踪精度提升10倍以上;图3是100次蒙特卡洛仿真实验后本发明和匀速运动(CV)模型的卡尔曼滤波方法、单渐消因子方法和多渐消因子方法进行比较目标位置的跟踪估计效果比较图,整体趋势上本发明在距离维的跟踪误差都要明显优于其他三种算法。
附图说明
图1是本发明的流程图;
图2是目标跟踪位置的对比图;
图3是目标的位置均方根误差随时间变化曲线;
图4是目标的速度均方根误差随时间变化曲线。
具体实施方式
具体实施方式一:如图1所示,本实施方式给出的一种基于多重渐消因子的自适应目标跟踪滤波方法,具体是按照以下步骤进行的:
步骤一,建立运动目标的状态模型和量测模型,将状态模型的运动状态标记为N维的运动状态向量X,量测模型的量测值标记为M维的量测向量Z,根据运动模型获得状态误差协方差矩阵P和过程噪声协方差矩阵Q,根据量测模型获得量测噪声协方差矩阵R;
步骤二,获得第k时刻的运动状态向量X(k)、状态误差协方差矩阵P(k),按照卡尔曼递推公式计算运动目标在第k+1时刻的运动状态向量预测值量测向量预测值并计算得到第k+1时刻的目标的新息协方差S(k+1);根据第k+1时刻目标的量测向量Z(k+1)计算第k+1时刻新息v(k+1),并由v(k+1)得到第k+1时刻新息协方差的估计值
步骤三,计算第k+1时刻的指数加权因子根据 和S(k+1)计算得到第k+1时刻多重渐消因子λ(k+1);
步骤四,通过第k+1时刻多重渐消因子λ(k+1)计算第k+1时刻的状态预测协方差矩阵P(k+1|k),并得到第k+1时刻的卡尔曼增益K(k+1),最终得到目标在第k+1时刻的滤波值和滤波协方差P(k+1|k+1);
步骤五,使k=k+1,判断目标跟踪是否结束,如果没有结束那么返回执行步骤二~步骤四,如果结束则完成所述目标跟踪滤波方法。
具体实施方式二:本实施方式与具体实施方式一不同的是:步骤一中目标的运动模型为:
X(k+1)=F(k)X(k)+G(k)V(k)
其中,F(k)是N×N阶的状态转移矩阵,G(k)为N×M维的过程噪声分布矩阵,V(k)为M维过程噪声,是具有零均值的高斯白噪声,过程噪声协方差矩阵为Q(k)。
具体实施方式三:本实施方式与具体实施方式二不同的是:步骤一中目标的量测模型为:
Z(k)=H(k)X(k)+W(k)
其中,H(k)是M×N阶的量测矩阵,W(k)为M维的量测噪声,是具有零均值的高斯白噪声,量测噪声协方差矩阵为R(k)。
具体实施方式四:本实施方式与具体实施方式三不同的是:步骤二具体为:当k=1时,采用传感器获得的量测向量Z(0)、Z(1),并初始化运动状态向量X(1)和状态误差协方差矩阵P(1);当k>1时,用k时刻的滤波值和滤波协方差P(k|k)初始化运动状态向量X(k)和状态误差协方差矩阵P(k),P(k)=P(k|k);第k+1时刻运动状态向量预测值量测向量预测值新息协方差S(k+1)、新息v(k+1)、新息协方差的估计值的计算方法具体为:
S(k+1)=H(k)(F(k)P(k|k)FT(k)+G(k)Q(k)GT(k))HT(k)+R(k+1)
其中,η(k)=diag(λ(1)1(k)/(1+λ(1)1(k)),…,λ(j)j(k)/(1+λ(j)j(k)),…,λ(M)M(k)/(1+λ(M)M(k))),Sk(j,j),代表S(k)和在第j行第j列的元素值。
其它步骤及参数与具体实施方式三相同。
具体实施方式五:本实施方式与具体实施方式四不同的是:步骤三的计算方法具体为:
第k+1时刻指数加权因子其中,b是调节因子,0≤b<1;
第k+1时刻多重渐消因子
其中,λj(k+1)代表第k+1时刻第j个通道的渐消因子,满足1≤j≤M,M是通道个数;是阶数为Mj×Mj的单位阵,Mj是每个子通道的维数,Mj=N/M,Sk+1(j,j),代表S(k+1)和在第j行第j列的元素值。
具体实施方式六:
本实施方式与具体实施方式五不同的是:步骤四的计算方法具体为:
P(k+1|k)=λ(k+1)(F(k)P(k|k)FT(k)+G(k)Q(k)GT(k))
K(k+1)=P(k+1|k)HT(k+1)(HT(k+1)P(k+1|k)HT(k+1)+R(k+1))-1
P(k+1|k+1)=P(k+1|k)-K(k+1)H(k+1)P(k+1|k)。
实施例
为展示本发明的实施方式的效果,设置一组仿真环境进行验证:
如图1所示,首先,在步骤一,建立运动目标的状态模型和量测模型。仿真目标在x-y平面作匀加速运动,目标的初始位置为(1000m,1000m),初始速度为(100m/s,50m/s),目标加速度为(-1m/s2,-3m/s2),对目标进行100s的观测,用匀速运动(CV)模型对其建模,则目标的运动模型为:
X(k+1)=F(k)X(k)+G(k)V(k)
式中,目标的运动状态向量:
将目标初始状态设置为:
X(0)=[1000 100 1000 50]T
系统的状态转移矩阵:
设雷达采样间隔T=1s;
过程噪声分布矩阵:
过程噪声协方差矩阵:
过程噪声分量设置为q1=q2=0.1;
从观测雷达处获得目标在在x-y方向的位置,雷达的测距误差为σx、σy,则量测模型为:
Z(k)=H(k)X(k)+W(k)
式中,量测向量:
Z(k)=[Z1(k),Z2(k)]T=[x,y]T
量测矩阵:
量测噪声协方差矩阵为:
其中,测距误差σx=50m,σy=50m;
接着,在步骤二,初始化第k时刻的运动状态向量X(k),状态误差协方差矩阵P(k),
当k=1时
其中,Z1(0)、Z2(0)为Z(0)中的元素,Z1(1)、Z2(1)为Z(1)的元素;
当k>1时
P(k)=P(k|k)
计算运动目标在第k+1个时刻的运动状态向量预测值
量测向量预测值
新息协方差S(k+1):
S(k+1)=H(k)(F(k)P(k|k)FT(k)+G(k)Q(k)GT(k))HT(k)+R(k+1)
新息v(k+1):
新息协方差的估计值
其中,η(k)=diag(λ(1)1(k)/(1+λ(1)1(k)),…,λ(j)j(k)/(1+λ(j)j(k)),…,λ(M)M(k)/(1+λ(M)M(k)));
然后,在步骤三,计算指数加权因子
其中,调节因子b=0.1
在本实施例中,运动状态向量中x,y方向的元素个数为2,且距离和速度之间为导数关系,可以乘以相同的渐消因子,因此Mj=2,多重渐消因子可以表示为:
其中,
之后,在步骤四,计算状态预测协方差矩阵P(k+1|k):
P(k+1|k)=λ(k+1)(F(k)P(k|k)FT(k)+G(k)Q(k)GT(k))
卡尔曼增益K(k+1):
K(k+1)=P(k+1|k)HT(k+1)(HT(k+1)P(k+1|k)HT(k+1)+R(k+1))-1
第k+1时刻的滤波值
第k+1时刻的滤波值和滤波协方差P(k+1|k+1):
P(k+1|k+1)=P(k+1|k)-K(k+1)H(k+1)P(k+1|k)
最后,在步骤五,返回执行步骤二,并使k=k+1,直到目标跟踪结束。
选取均方根误差(RMSE)作为性能的量度,进行100次蒙特卡罗仿真实验。
将本发明提出的基于多重渐消因子的自适应目标跟踪滤波方法与匀速运动(CV)模型的卡尔曼滤波方法、单渐消因子方法和多渐消因子方法进行比较。
仿真结果如图2到图4所示
图2是本发明和匀速运动(CV)模型的卡尔曼滤波方法的跟踪位置对比图,由于模型失配,卡尔曼滤波算法的滤波轨迹逐渐偏离了目标的实际位置,而本发明的滤波轨迹结果始终没有发散,一直与目标运动趋势一致。
图3是100次蒙特卡洛仿真实验后本发明和其他三种算法目标位置的跟踪估计效果比较图。从图中可看出,卡尔曼滤波方法由于模型失配,导致滤波结果发散,因此该算法在距离维的跟踪误差最大;单渐消因子方法由于采用标量的单渐消因子,状态预测协方差矩阵的各个通道仅具有相同的调节能力,虽然最终滤波结果没有发散,但距离维的跟踪误差较大;多渐消因子方法采用多重的渐消因子,跟踪性能有所提升,但是由于其只在距离维乘以了渐消因子,而在速度维不作调整,没有充分利用目标距离和速度的关系,而本发明方法在距离维和速度维都进行了调整,因此多渐消因子方法在距离维的跟踪误差要大于本发明方法。综上所述,整体趋势上本发明在距离维的跟踪误差都要明显优于其他三种算法。
图4是100次蒙特卡洛仿真实验后本发明和其他三种算法目标速度的跟踪估计效果对比图。同样地,从图中可看出,卡尔曼滤波方法由于模型失配,导致滤波结果发散,因此该算法在速度维的跟踪误差最大;单渐消因子方法由于采用标量的单渐消因子,状态预测协方差矩阵的各个通道仅具有相同的调节能力,虽然最终滤波结果没有发散,但速度维的跟踪误差同样较大;多渐消因子方法由于只在距离维乘以了渐消因子,而在速度维不作调整,没有充分利用目标距离和速度的关系,所以虽然该方法在距离维跟踪效果较好,但是速度维的跟踪效果较差,误差甚至大于单渐消因子方法;本发明方法在距离维和速度维都进行了调整,因子整体趋势上本发明在速度维的跟踪误差都要明显优于其他三种算法。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (6)
1.一种基于多重渐消因子的自适应目标跟踪滤波方法,其特征在于,所述方法包括如下步骤:
步骤一,建立运动目标的状态模型和量测模型,将状态模型的运动状态标记为N维的运动状态向量X,量测模型的量测值标记为M维的量测向量Z,根据运动模型获得状态误差协方差矩阵P和过程噪声协方差矩阵Q,根据量测模型获得量测噪声协方差矩阵R;
步骤二,获得第k时刻的运动状态向量X(k)、状态误差协方差矩阵P(k),按照卡尔曼递推公式计算运动目标在第k+1时刻的运动状态向量预测值量测向量预测值并计算得到第k+1时刻的目标的新息协方差S(k+1);根据第k+1时刻目标的量测向量Z(k+1)计算第k+1时刻新息v(k+1),并由v(k+1)得到第k+1时刻新息协方差的估计值
步骤三,计算第k+1时刻的指数加权因子根据和S(k+1)计算得到第k+1时刻多重渐消因子λ(k+1);
步骤四,通过第k+1时刻多重渐消因子λ(k+1)计算第k+1时刻的状态预测协方差矩阵P(k+1|k),并得到第k+1时刻的卡尔曼增益K(k+1),最终得到目标在第k+1时刻的滤波值和滤波协方差P(k+1|k+1);
步骤五,使k=k+1,判断目标跟踪是否结束,如果没有结束那么返回执行步骤二~步骤四,如果结束则完成所述目标跟踪滤波方法。
2.根据权利要求1所述的一种基于多重渐消因子的自适应目标跟踪滤波方法,其特征在于,所述步骤一中目标的运动模型为:
X(k+1)=F(k)X(k)+G(k)V(k)
其中,F(k)是N×N阶的状态转移矩阵,G(k)为N×M维的过程噪声分布矩阵,V(k)为M维过程噪声,是具有零均值的高斯白噪声,过程噪声协方差矩阵为Q(k)。
3.根据权利要求2所述的一种基于多重渐消因子的自适应目标跟踪滤波方法,其特征在于,所述步骤一中目标的量测模型为:
Z(k)=H(k)X(k)+W(k)
其中,H(k)是M×N阶的量测矩阵,W(k)为M维的量测噪声,是具有零均值的高斯白噪声,量测噪声协方差矩阵为R(k)。
4.根据权利要求3所述的一种基于多重渐消因子的自适应目标跟踪滤波方法,其特征在于,所述步骤二具体为:当k=1时,采用传感器获得的量测向量Z(0)、Z(1),并初始化运动状态向量X(1)和状态误差协方差矩阵P(1);当k>1时,用k时刻的滤波值和滤波协方差P(k|k)初始化运动状态向量X(k)和状态误差协方差矩阵P(k),P(k)=P(k|k);第k+1时刻运动状态向量预测值量测向量预测值新息协方差S(k+1)、新息v(k+1)、新息协方差的估计值的计算方法具体为:
<mrow>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>Z</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
S(k+1)=H(k)(F(k)P(k|k)FT(k)+G(k)Q(k)GT(k))HT(k)+R(k+1)
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>Z</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>S</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mi>v</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mi>v</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>k</mi>
<mo>&GreaterEqual;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,η(k)=diag(λ(1)1(k)/(1+λ(1)1(k)),…,λ(j)j(k)/(1+λ(j)j(k)),…,λ(M)M(k)/(1+λ(M)M(k))),Sk(j,j),代表S(k)和在第j行第j列的元素值。
5.根据权利要求4所述的一种基于多重渐消因子的自适应目标跟踪滤波方法,其特征在于,所述步骤三的计算方法具体为:
第k+1时刻指数加权因子其中,b是调节因子,0≤b<1;
第k+1时刻多重渐消因子
其中,λj(k+1)代表第k+1时刻第j个通道的渐消因子,满足1≤j≤M,M是通道个数; 是阶数为Mj×Mj的单位阵,Mj是每个子通道的维数,Mj=N/M,Sk+1(j,j),代表S(k+1)和在第j行第j列的元素值。
6.根据权利要求5所述的一种基于多重渐消因子的自适应目标跟踪滤波方法,其特征在于,所述步骤四的计算方法具体为:
P(k+1|k)=λ(k+1)(F(k)P(k|k)FT(k)+G(k)Q(k)GT(k))
K(k+1)=P(k+1|k)HT(k+1)(HT(k+1)P(k+1|k)HT(k+1)+R(k+1))-1
<mrow>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
P(k+1|k+1)=P(k+1|k)-K(k+1)H(k+1)P(k+1|k)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710752493.8A CN107561503B (zh) | 2017-08-28 | 2017-08-28 | 一种基于多重渐消因子的自适应目标跟踪滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710752493.8A CN107561503B (zh) | 2017-08-28 | 2017-08-28 | 一种基于多重渐消因子的自适应目标跟踪滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107561503A true CN107561503A (zh) | 2018-01-09 |
CN107561503B CN107561503B (zh) | 2020-08-14 |
Family
ID=60977568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710752493.8A Active CN107561503B (zh) | 2017-08-28 | 2017-08-28 | 一种基于多重渐消因子的自适应目标跟踪滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107561503B (zh) |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108646237A (zh) * | 2018-05-09 | 2018-10-12 | 西安电子科技大学 | 基于当前统计模型的雷达机动目标跟踪优化方法 |
CN109163720A (zh) * | 2018-08-27 | 2019-01-08 | 广西科技大学 | 基于渐消记忆指数加权的卡尔曼滤波跟踪方法 |
CN109188422A (zh) * | 2018-08-08 | 2019-01-11 | 中国航空工业集团公司雷华电子技术研究所 | 一种基于lu分解的卡尔曼滤波目标跟踪方法 |
CN109459033A (zh) * | 2018-12-21 | 2019-03-12 | 哈尔滨工程大学 | 一种多重渐消因子的机器人无迹快速同步定位与建图方法 |
CN109798888A (zh) * | 2019-03-15 | 2019-05-24 | 京东方科技集团股份有限公司 | 移动设备的姿态确定装置、方法和视觉里程计 |
CN110061716A (zh) * | 2019-01-15 | 2019-07-26 | 河海大学 | 一种基于最小二乘和多重渐消因子的改进kalman滤波方法 |
CN110824420A (zh) * | 2019-11-14 | 2020-02-21 | 中国人民解放军火箭军工程大学 | 一种分布式目标跟踪方法 |
CN110989856A (zh) * | 2019-11-07 | 2020-04-10 | 北京集创北方科技股份有限公司 | 一种坐标预测方法、装置、设备及可存储介质 |
CN111044054A (zh) * | 2020-01-06 | 2020-04-21 | 哈尔滨工业大学 | 一种基于单目slam算法利用卡尔曼滤波进行位姿优化的方法 |
CN111077518A (zh) * | 2019-12-20 | 2020-04-28 | 哈尔滨工业大学 | 一种基于距离-多普勒量测的跟踪滤波方法及装置 |
CN111865268A (zh) * | 2020-07-10 | 2020-10-30 | 南京理工大学 | 一种简易快速的卡尔曼滤波方法 |
CN111884555A (zh) * | 2020-07-29 | 2020-11-03 | 江南大学 | 一种永磁同步电机转子转速与位置的滤波估计方法 |
CN112284414A (zh) * | 2020-10-10 | 2021-01-29 | 东南大学 | 一种基于多渐消因子自适应动机座旋转调制精对准方法 |
CN112835027A (zh) * | 2021-01-21 | 2021-05-25 | 同济大学 | 一种基于强跟踪滤波的雷达目标跟踪方法 |
CN112912759A (zh) * | 2018-10-24 | 2021-06-04 | 株式会社电装 | 物体跟踪装置 |
CN113030959A (zh) * | 2021-04-01 | 2021-06-25 | 成都汇蓉国科微系统技术有限公司 | 一种雷达和光电联动反低慢小目标3d地图画面消抖方法 |
CN113075652A (zh) * | 2021-02-25 | 2021-07-06 | 中国人民解放军空军预警学院 | 一种高超声速飞行器三维跟踪方法 |
CN113156383A (zh) * | 2021-04-21 | 2021-07-23 | 杭州加速科技有限公司 | 基于国产fpga的雷达信号抗干扰方法、系统及数据处理设备 |
CN113155122A (zh) * | 2021-04-01 | 2021-07-23 | 广州大学 | 基于自适应滤波的机动目标跟踪方法 |
CN113253615A (zh) * | 2021-06-22 | 2021-08-13 | 季华实验室 | 基于分布式电动底盘的运动状态观测方法及系统 |
CN115685128A (zh) * | 2022-11-14 | 2023-02-03 | 中国人民解放军空军预警学院 | 一种机动目标场景下的雷达目标跟踪算法及电子设备 |
CN115964603A (zh) * | 2023-02-10 | 2023-04-14 | 成都理工大学 | 一种基于改进Kalman滤波的机动目标跟踪方法 |
CN116736287A (zh) * | 2023-06-19 | 2023-09-12 | 哈尔滨工业大学 | 基于多探测模式的地波多径目标跟踪方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101252531B1 (ko) * | 2013-01-22 | 2013-04-09 | 국방과학연구소 | 적응적 칼만필터를 이용한 목표 위치추적장치 |
CN104408744A (zh) * | 2014-11-17 | 2015-03-11 | 电子科技大学 | 一种用于目标跟踪的强跟踪容积卡尔曼滤波方法 |
CN107045125A (zh) * | 2017-03-17 | 2017-08-15 | 电子科技大学 | 一种基于预测值量测转换的交互多模型雷达目标跟踪方法 |
-
2017
- 2017-08-28 CN CN201710752493.8A patent/CN107561503B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101252531B1 (ko) * | 2013-01-22 | 2013-04-09 | 국방과학연구소 | 적응적 칼만필터를 이용한 목표 위치추적장치 |
CN104408744A (zh) * | 2014-11-17 | 2015-03-11 | 电子科技大学 | 一种用于目标跟踪的强跟踪容积卡尔曼滤波方法 |
CN107045125A (zh) * | 2017-03-17 | 2017-08-15 | 电子科技大学 | 一种基于预测值量测转换的交互多模型雷达目标跟踪方法 |
Non-Patent Citations (2)
Title |
---|
徐小良等: "基于量化新息的容积粒子滤波融合目标跟踪算法", 《自动化学报》 * |
薛海建等: "基于自适应多重渐消因子卡尔曼滤波的SINS初始对准方法", 《系统工程与电子技术》 * |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108646237A (zh) * | 2018-05-09 | 2018-10-12 | 西安电子科技大学 | 基于当前统计模型的雷达机动目标跟踪优化方法 |
CN109188422A (zh) * | 2018-08-08 | 2019-01-11 | 中国航空工业集团公司雷华电子技术研究所 | 一种基于lu分解的卡尔曼滤波目标跟踪方法 |
CN109163720A (zh) * | 2018-08-27 | 2019-01-08 | 广西科技大学 | 基于渐消记忆指数加权的卡尔曼滤波跟踪方法 |
CN112912759A (zh) * | 2018-10-24 | 2021-06-04 | 株式会社电装 | 物体跟踪装置 |
CN109459033A (zh) * | 2018-12-21 | 2019-03-12 | 哈尔滨工程大学 | 一种多重渐消因子的机器人无迹快速同步定位与建图方法 |
CN110061716A (zh) * | 2019-01-15 | 2019-07-26 | 河海大学 | 一种基于最小二乘和多重渐消因子的改进kalman滤波方法 |
CN109798888B (zh) * | 2019-03-15 | 2021-09-17 | 京东方科技集团股份有限公司 | 移动设备的姿态确定装置、方法和视觉里程计 |
CN109798888A (zh) * | 2019-03-15 | 2019-05-24 | 京东方科技集团股份有限公司 | 移动设备的姿态确定装置、方法和视觉里程计 |
CN110989856A (zh) * | 2019-11-07 | 2020-04-10 | 北京集创北方科技股份有限公司 | 一种坐标预测方法、装置、设备及可存储介质 |
CN110989856B (zh) * | 2019-11-07 | 2024-03-22 | 北京集创北方科技股份有限公司 | 一种坐标预测方法、装置、设备及可存储介质 |
CN110824420A (zh) * | 2019-11-14 | 2020-02-21 | 中国人民解放军火箭军工程大学 | 一种分布式目标跟踪方法 |
CN111077518A (zh) * | 2019-12-20 | 2020-04-28 | 哈尔滨工业大学 | 一种基于距离-多普勒量测的跟踪滤波方法及装置 |
CN111077518B (zh) * | 2019-12-20 | 2020-11-06 | 哈尔滨工业大学 | 一种基于距离-多普勒量测的跟踪滤波方法及装置 |
CN111044054A (zh) * | 2020-01-06 | 2020-04-21 | 哈尔滨工业大学 | 一种基于单目slam算法利用卡尔曼滤波进行位姿优化的方法 |
CN111865268A (zh) * | 2020-07-10 | 2020-10-30 | 南京理工大学 | 一种简易快速的卡尔曼滤波方法 |
CN111884555A (zh) * | 2020-07-29 | 2020-11-03 | 江南大学 | 一种永磁同步电机转子转速与位置的滤波估计方法 |
CN111884555B (zh) * | 2020-07-29 | 2021-11-12 | 山东宇恒智能动力科技有限公司 | 一种永磁同步电机转子转速与位置的滤波估计方法 |
CN112284414A (zh) * | 2020-10-10 | 2021-01-29 | 东南大学 | 一种基于多渐消因子自适应动机座旋转调制精对准方法 |
WO2022073376A1 (zh) * | 2020-10-10 | 2022-04-14 | 东南大学 | 一种基于多渐消因子自适应动机座旋转调制精对准方法 |
CN112835027A (zh) * | 2021-01-21 | 2021-05-25 | 同济大学 | 一种基于强跟踪滤波的雷达目标跟踪方法 |
CN113075652A (zh) * | 2021-02-25 | 2021-07-06 | 中国人民解放军空军预警学院 | 一种高超声速飞行器三维跟踪方法 |
CN113030959A (zh) * | 2021-04-01 | 2021-06-25 | 成都汇蓉国科微系统技术有限公司 | 一种雷达和光电联动反低慢小目标3d地图画面消抖方法 |
CN113155122A (zh) * | 2021-04-01 | 2021-07-23 | 广州大学 | 基于自适应滤波的机动目标跟踪方法 |
CN113030959B (zh) * | 2021-04-01 | 2023-07-04 | 成都汇蓉国科微系统技术有限公司 | 一种雷达和光电联动反低慢小目标3d地图画面消抖方法 |
CN113156383B (zh) * | 2021-04-21 | 2024-02-20 | 杭州加速科技有限公司 | 基于fpga的雷达信号抗干扰方法、系统及数据处理设备 |
CN113156383A (zh) * | 2021-04-21 | 2021-07-23 | 杭州加速科技有限公司 | 基于国产fpga的雷达信号抗干扰方法、系统及数据处理设备 |
CN113253615A (zh) * | 2021-06-22 | 2021-08-13 | 季华实验室 | 基于分布式电动底盘的运动状态观测方法及系统 |
CN115685128A (zh) * | 2022-11-14 | 2023-02-03 | 中国人民解放军空军预警学院 | 一种机动目标场景下的雷达目标跟踪算法及电子设备 |
CN115685128B (zh) * | 2022-11-14 | 2023-04-11 | 中国人民解放军空军预警学院 | 一种机动目标场景下的雷达目标跟踪算法及电子设备 |
CN115964603A (zh) * | 2023-02-10 | 2023-04-14 | 成都理工大学 | 一种基于改进Kalman滤波的机动目标跟踪方法 |
CN115964603B (zh) * | 2023-02-10 | 2023-05-30 | 成都理工大学 | 一种基于改进Kalman滤波的机动目标跟踪方法 |
CN116736287A (zh) * | 2023-06-19 | 2023-09-12 | 哈尔滨工业大学 | 基于多探测模式的地波多径目标跟踪方法 |
CN116736287B (zh) * | 2023-06-19 | 2024-04-05 | 哈尔滨工业大学 | 基于多探测模式的地波多径目标跟踪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107561503B (zh) | 2020-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107561503B (zh) | 一种基于多重渐消因子的自适应目标跟踪滤波方法 | |
CN111985093B (zh) | 一种带噪声估计器的自适应无迹卡尔曼滤波状态估计方法 | |
CN107728138B (zh) | 一种基于当前统计模型的机动目标跟踪方法 | |
CN110503071B (zh) | 基于变分贝叶斯标签多伯努利叠加模型的多目标跟踪方法 | |
CN108470089B (zh) | 一种基于最小二乘样本拟合的复信号时延估计方法 | |
CN108682023A (zh) | 基于Elman神经网络的紧耦合无迹卡尔曼跟踪滤波算法 | |
CN105894033B (zh) | 一种海杂波背景下的弱目标检测方法及系统 | |
CN110673089B (zh) | 未知视距和非视距分布情况下基于到达时间的定位方法 | |
CN109581281B (zh) | 基于到达时间差和到达频率差的移动目标定位方法 | |
CN106980104B (zh) | 用于传感器阵列的信号波达方向自校正方法 | |
CN109671100B (zh) | 一种分布式可变扩散组合系数粒子滤波直接跟踪方法 | |
CN109471061B (zh) | 一种鲁棒处理模型参数误差的接收信号强度差定位方法 | |
CN103237320B (zh) | 无线传感器网络基于混合量化卡尔曼融合的目标跟踪方法 | |
CN112881979A (zh) | 一种基于ekf滤波的初始状态自适应融合定位方法 | |
CN110716203A (zh) | 一种被动声纳目标的时频分析与跟踪方法 | |
Hasan et al. | Adaptive α-β-filter for target tracking using real time genetic algorithm | |
CN107966676A (zh) | 复杂噪声环境下阵列天线角度和信源个数联合估计方法 | |
CN112446010B (zh) | 自适应弱敏秩卡尔曼滤波方法及其应用 | |
CN109239681B (zh) | 一种运动干扰稳健滤波方法 | |
CN115469314A (zh) | 一种均匀圆环阵稳健水下目标方位跟踪方法及系统 | |
CN107544050A (zh) | 白噪声背景下一种构造自适应阈值估计信号源数目的方法 | |
CN105022025B (zh) | 基于稀疏处理的信号波达方向估计方法 | |
CN109375160B (zh) | 纯方位无源定位中一种测角误差估计方法 | |
CN109188420B (zh) | 基于深度长短期记忆网络的窄带雷达目标跟踪方法 | |
Jin et al. | Iterative doubly constrained robust Capon beamformer |
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 |