CN106487358B - 一种机动目标转弯跟踪方法 - Google Patents
一种机动目标转弯跟踪方法 Download PDFInfo
- Publication number
- CN106487358B CN106487358B CN201610871781.0A CN201610871781A CN106487358B CN 106487358 B CN106487358 B CN 106487358B CN 201610871781 A CN201610871781 A CN 201610871781A CN 106487358 B CN106487358 B CN 106487358B
- Authority
- CN
- China
- Prior art keywords
- state
- calculating
- covariance matrix
- prediction
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 239000011159 matrix material Substances 0.000 claims abstract description 56
- 238000000034 method Methods 0.000 claims abstract description 47
- 238000012417 linear regression Methods 0.000 claims abstract description 11
- 238000001914 filtration Methods 0.000 claims abstract description 9
- 238000005070 sampling Methods 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000002123 temporal effect Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 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
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0202—Two or more dimensional filters; Filters for complex signals
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0202—Two or more dimensional filters; Filters for complex signals
- H03H2017/0205—Kalman filters
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Feedback Control In General (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于信号处理技术领域,提供了一种机动目标转弯跟踪方法,通过雷达测得机动目标的距离和角度,并结合目标物的转弯特性构建状态空间模型经过初始化滤波方法、时间预测、定点迭代和更新误差协方差矩阵等步骤后,能够对机动目标转弯跟踪的非线性系统状态进行估计。该方法在最大相关熵卡尔曼滤波器基础上,利用统计学中的统计线性回归方法,将最大相关熵代价函数融入容积卡尔曼滤波器的框架中,从而应用到非线性系统状态估计上。该方法在做非线性系统状态估计的时候较传统的容积卡尔曼方法,有明显的鲁邦性和和抗噪能力的提升。
Description
技术领域
本发明涉及信号处理技术领域,具体涉及一种机动目标转弯跟踪方法。
背景技术
系统的状态估计是信号处理中的一个重要问题。卡尔曼滤波方法是解决线性系统的状态估计的主要方法之一,其充分利用了系统的状态空间模型和观测数据,通过求解优化问题,使得状态估计的误差最小,从而来获得系统状态的最优估计。非线性卡尔曼滤波器包括无迹卡尔曼滤波器,容积卡尔曼滤波器,以及稀疏网格求积分滤波器,这三种滤波器都是解决非线性系统的状态估计的有效工具。
状态空间模型为:
其中f(·)、h(·)分别表示状态向量、观测向量、状态转移函数和观测函数,m表示状态向量的维数,n表示观测向量的维数,vk和wk为随机噪声项,i、j表示变量。
估计效果的好坏可以用RMSE指标来衡量,其定义为:
其中,xi,k表示真实的状态,表示对状态的估计,下标k表示第k时刻,下标i表示第i次蒙特卡罗,各个时刻RMSE的数值越小表示估计效果越好。
容积卡尔曼滤波器,作为状态估计中的一种有效方法,广泛应用于跟踪、导航。然而从贝叶斯滤波框架中得出的容积卡尔曼滤波,为了应用数值积分方法,需要假设噪声项为高斯分布,这严重影响了容积卡尔曼滤波在实际应用中的抗噪能力。同时也由于系统内部的变化不定,目前考虑了鲁邦性的滤波方法都是对特定应用定制特定的方法。
发明内容
针对现有技术中的缺陷,本发明提供一种机动目标转弯跟踪方法,适用于非线性系统,提高了系统的鲁棒性和抗噪能力。
一种机动目标转弯跟踪方法,基于统计线性回归的最大相关熵容积卡尔曼滤波方法实现,包括以下步骤:
S1:通过雷达测得机动目标的距离和角度,并结合目标物的转弯特性构建如下的状态空间模型:
其中,xk和yk是位置,和代表x和y方向上的速度,rk和θk分别表示雷达测得的距离和角度,ωk是转弯速度;T=1,(sx,sy)=(0,0),vk和wk为0;
S2:初始化滤波:获取从高斯分布中随机产生的初始状态估计;
S3:时间预测:根据上一个时刻的状态估计和状态估计误差协方差矩阵计算第一采样点,并根据第一采样点计算经过状态转移函数变换后得到的第一变换采样点、状态预测以及预测误差协方差矩阵;
S4:定点迭代:基于时间预测得到的状态预测以及预测误差协方差矩阵计算第二采样点,并根据第二采样点计算经过观测函数变换后得到的第二变换采样点、观测预测值和观测估计误差协方差矩阵;利用统计线性回归和定点迭代的方法计算状态估计;
S5:更新状态估计误差协方差矩阵,并在获得带噪观测值后返回步骤S3循环。
优选地,所述步骤S2具体为:从均值为协方差矩阵为P0的高斯分布中随机产生初始状态估计令k←1,熵计算所用核函数取高斯核函数其中σ为高斯核函数的核宽,k为离散时刻。
优选地,所述步骤S3包括:
S3a:获取上一个时刻的状态估计和状态估计误差协方差矩阵Pk-1,并计算第一采样点Xi,k-1:其中Sk-1是Pk-1的平方根矩阵,ej是第j个元素为1的单位列向量;n为k时刻下的状态向量xk的维数,i取值1~2n;
S3b:根据下式计算经过状态转移函数xk=f(xk-1)变换后的第一变换采样点和状态预测
S3c:计算状态预测对应的预测误差协方差矩阵Pk|k-1:
其中,Qk-1为状态转移噪声的协方差矩阵。
优选地,所述步骤S4包括:
S4a:基于时间预测得到的状态预测和对应的预测误差协方差矩阵Pk|k-1重新计算第二采样点Xi,k|k-1,其中Sk|k-1是Pk|k-1的平方根矩阵;
S4b:根据下式计算经过观测函数zk=h(xk)变换后的第二变换采样点和观测预测值
S4c:计算观测估计误差协方差矩阵Pxz;
S4d:计算统计线性回归系数矩阵并定义以下关系:
dk=Wkxk+ek;
其中rk为观测噪声,zk为k时刻的带噪观测值;I为单位矩阵;另外φk=diag([φp,k|k-1,φr]),φp,k|k-1和φr分别表示Pk|k-1和观测噪声的协方差矩阵Rk的平方根矩阵;
S4e:令j←1,采用定点迭代的方法计算状态估计:
其中,
上标j表示第j次迭代,对角系数矩阵和中对角元素由误差元素的第l个元素通过核函数κσ计算到,其中l=1~(n+m),
S4f:计算终止条件判断是否小于常数ε,如果是,记录迭代次数为Nk,并终止迭代,执行步骤S5;否则继续S4e。
优选地,所述步骤S5包括:
S5a:更新状态估计误差协方差矩阵
S5b:令k=k+1,返回步骤S3。
由上述技术方案可知,本发明提供的机动目标转弯跟踪方法,在最大相关熵卡尔曼滤波器基础上,利用统计学中的统计线性回归方法,将最大相关熵代价函数融入容积卡尔曼滤波器的框架中,从而应用到机动目标转弯跟踪的非线性系统状态估计上。该方法在做非线性系统状态估计的时候较传统的容积卡尔曼方法,有明显的鲁邦性和和抗噪能力的提升。
附图说明
图1为实施例1中CKF和MCCKF-SLR对转弯机动目标跟踪效果图。
图2为实施例1中CKF和MCCKF-SLR对位置估计误差的对比图。
图3为实施例2中高斯噪声项增加了新噪声项后仿真得到的时间序列的估计误差对比图。
图4为实施例2中高斯噪声项增加了噪声项后仿真得到的时间序列的估计误差对比图。
具体实施方式
下面将结合附图对本发明技术方案的实施例进行详细的描述。以下实施例仅用于更加清楚地说明本发明的技术方案,因此只作为示例,而不能以此来限制本发明的保护范围。
需要注意的是,除非另有说明,本申请使用的技术术语或者科学术语应当为本发明所属领域技术人员所理解的通常意义。
由于现有卡尔曼滤波器存在抗噪能力低、鲁邦性差等缺陷,因此考虑替换卡尔曼滤波器求解优化过程的代价函数,是一个很好的解决这两个问题的办法。最近基于最大相关熵代价函数提出的最大相关熵卡尔曼滤波器,作为新的线性滤波器较传统的卡尔曼滤波器在鲁邦性和抗噪能力方面提高明显,而迄今还没有人在适用于线性系统的最大相关熵卡卡尔曼方法基础上,提出适用于非线性系统的最大相关熵卡尔曼方法。
一种机动目标转弯跟踪方法,包括以下步骤:
S1:通过雷达测得机动目标的距离和角度,并结合目标物的转弯特性构建如下的状态空间模型:
其中,xk和yk是位置,和代表x和y方向上的速度,rk和θk分别表示雷达测得的距离和角度,ωk是转弯速度;T=1,(sx,sy)=(0,0),vk和wk为0;
S2:初始化滤波:从均值为协方差矩阵为P0的高斯分布中随机产生初始状态估计令k←1,熵计算所用核函数取高斯核函数其中σ为高斯核函数的核宽,k为离散时刻。
S3:时间预测:
S3a:获取上一个时刻的状态估计和状态估计误差协方差矩阵Pk-1,并计算第一采样点Xi,k-1:其中Sk-1是Pk-1的平方根矩阵,ej是第j个元素为1的单位列向量;n为k时刻下的状态向量xk的维数,i取值1~2n;
S3b:根据下式计算经过状态转移函数xk=f(xk-1)变换后的第一变换采样点和状态预测
S3c:计算状态预测对应的预测误差协方差矩阵Pk|k-1:
其中,Qk-1为状态转移噪声的协方差矩阵;
S4:定点迭代:
S4a:基于时间预测得到的状态预测和对应的预测误差协方差矩阵Pk|k-1重新计算第二采样点Xi,k|k-1,其中Sk|k-1是Pk|k-1的平方根矩阵;
S4b:根据下式计算经过观测函数zk=h(xk)变换后的第二变换采样点和观测预测值
S4c:计算观测估计误差协方差矩阵Pxz;
S4d:计算统计线性回归系数矩阵并定义以下关系:
dk=Wkxk+ek;
其中rk为观测噪声,zk为k时刻的带噪观测值;I为单位矩阵;另外φk=diag([φp,k|k-1,φr]),φp,k|k-1和φr分别表示Pk|k-1和观测噪声的协方差矩阵Rk的平方根矩阵;
S4e:令j←1,采用定点迭代的方法计算状态估计:
其中,
上标j表示第j次迭代,对角系数矩阵和中对角元素由误差元素
的第l个元素通过核函数κσ计算到,其中l=1~(n+m),
S4f:计算终止条件判断是否小于常数ε,如果是,记录迭代次数为Nk,并终止迭代,执行步骤S5;否则继续S4e。
S5:更新状态估计误差协方差矩阵:
S5a:更新状态估计误差协方差矩阵
S5b:令k=k+1,返回步骤S3。
该方法在最大相关熵卡尔曼滤波器基础上,利用统计学中的统计线性回归方法,将最大相关熵代价函数融入容积卡尔曼滤波器的框架中,从而应用到非线性系统状态估计上。能够在状态估计过程中提供的很好的鲁邦性,而且在非高斯噪声环境中也表现明显的优势,性能较传统的容积卡尔曼算法有很大的提升,可广泛用于非线性系统,尤其系统的噪声项是在非高斯情况下。
该方法利用了统计学中的统计线性回归(SLR)做线性化,将最大相关熵代价函数融入了卡尔曼滤波器(CKF)的框架中,因此该算法继承CKF的优点,可以直接应用于非线性系统的状态估计问题中。
算法分为时间预测和观测更新两部分,观测更新部分选取了定点迭代的方法来多步计算状态更新,直到满足迭代收敛的条件。定点迭代过程,选取预测更新的结果作为观测更新部分的定点迭代算法的初值,当满足终止条件终止迭代;经验表明,一般3到4步迭代即可。
在实际情况中,很多突变是突发的,而且噪声的统计特性很难被精确建模,因而我们的发明相比传统的容积卡尔曼滤波器有更好实际应用价值。
实施例1:
如图1-2所示,离散状态空间模型如下:
其中xk和yk是位置,和代表x和y方向上的速度,rk和θk分别表示雷达测得的距离和角度,转弯速度ωk仿真中设置为:
仿真中的其他参数设置为T=1,(sx,sy)=(0,0)。vk和wk是均值为0,协方差矩阵为Q和R的高斯噪声,其中,
而
首先经过初始化滤波算法。状态实际初值和协方差矩阵初值P0分别设置为P0=diag([100,10,100,10,1×10-4]),状态估计的初值从高斯分布中随机选取。高斯核函数中的核宽σ设置为5。然后依次经过时间预测、定点迭代和观测更新处理后,得到图1、2的仿真图,本发明提供的方法的仿真结果见图中曲线MCCKF-SLR,可见较传统CKF,性能优势明显。其中常数ε=1×10-6。
实施例2:
更换仿真模型为单变量非平稳增长模型。状态空间模型如下:
其中qk=1。根据上述方法处理后进行仿真,即可得到图3和4。状态的初始值和方差分别设置为P0=1。和实施例1相比,实施例2改变了模型中的噪声rk的统计特性。在图3的仿真中,以0.5的概率在原来的模型的噪声基础上增加了新的噪声项σr=1。用相同方法在图4的仿真中加的噪声项为从图3和4中可以看出,随着混合高斯噪声的非高斯程度的增加,该算法的性能较CKF有越来越大的提升。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围,其均应涵盖在本发明的权利要求和说明书的范围当中。
Claims (5)
1.一种机动目标转弯跟踪方法,其特征在于,基于统计线性回归的最大相关熵容积卡尔曼滤波方法实现,包括以下步骤:
S1:通过雷达测得机动目标的距离和角度,并结合目标物的转弯特性构建如下的状态空间模型:
其中,xk和yk是位置,和代表x和y方向上的速度,rk和θk分别表示雷达测得的距离和角度,ωk是转弯速度;T=1,(sx,sy)=(0,0),vk和wk为0;
S2:初始化滤波:获取从高斯分布中随机产生的初始状态估计;
S3:时间预测:根据上一个时刻的状态估计和状态估计误差协方差矩阵计算第一采样点,并根据第一采样点计算经过状态转移函数变换后得到的第一变换采样点、状态预测以及预测误差协方差矩阵;
S4:定点迭代:基于时间预测得到的状态预测以及预测误差协方差矩阵计算第二采样点,并根据第二采样点计算经过观测函数变换后得到的第二变换采样点、观测预测值和观测估计误差协方差矩阵;利用统计线性回归和定点迭代的方法计算状态估计;
S5:更新状态估计误差协方差矩阵,并在获得带噪观测值后返回步骤S3循环。
2.根据权利要求1所述的机动目标转弯跟踪方法,其特征在于,所述步骤S2具体为:从均值为协方差矩阵为P0的高斯分布中随机产生初始状态估计令k←1,熵计算所用核函数取高斯核函数其中σ为高斯核函数的核宽,k为离散时刻。
3.根据权利要求2所述的机动目标转弯跟踪方法,其特征在于,所述步骤S3包括:
S3a:获取上一个时刻的状态估计和状态估计误差协方差矩阵Pk-1,并计算第一采样点Xi,k-1:其中Sk-1是Pk-1的平方根矩阵,ej是第j个元素为1的单位列向量;n为k时刻下的状态向量xk的维数,i取值1~2n;
S3b:根据下式计算经过状态转移函数xk=f(xk-1)变换后的第一变换采样点和状态预测
S3c:计算状态预测对应的预测误差协方差矩阵Pk|k-1:
其中,Qk-1为状态转移噪声的协方差矩阵。
4.根据权利要求3所述的机动目标转弯跟踪方法,其特征在于,所述步骤S4包括:
S4a:基于时间预测得到的状态预测和对应的预测误差协方差矩阵Pk|k-1重新计算第二采样点Xi,k|k-1,其中Sk|k-1是Pk|k-1的平方根矩阵;
S4b:根据下式计算经过观测函数zk=h(xk)变换后的第二变换采样点和观测预测值
S4c:计算观测估计误差协方差矩阵Pxz;
S4d:计算统计线性回归系数矩阵并定义以下关系:
dk=Wkxk+ek;
其中rk为观测噪声,zk为k时刻的带噪观测值;I为单位矩阵;另外φk=diag([φp,k|k-1,φr]),φp,k|k-1和φr分别表示Pk|k-1和观测噪声的协方差矩阵Rk的平方根矩阵;
S4e:令j←1,采用定点迭代的方法计算状态估计:
其中,
上标j表示第j次迭代,对角系数矩阵和中对角元素由误差元素的第l个元素通过核函数κσ计算到,其中l=1~(n+m),
S4f:计算终止条件判断是否小于常数ε,如果是,记录迭代次数为Nk,并终止迭代,执行步骤S5;否则继续S4e。
5.根据权利要求4所述的机动目标转弯跟踪方法,其特征在于,所述步骤S5包括:
S5a:更新状态估计误差协方差矩阵
S5b:令k=k+1,返回步骤S3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610871781.0A CN106487358B (zh) | 2016-09-30 | 2016-09-30 | 一种机动目标转弯跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610871781.0A CN106487358B (zh) | 2016-09-30 | 2016-09-30 | 一种机动目标转弯跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106487358A CN106487358A (zh) | 2017-03-08 |
CN106487358B true CN106487358B (zh) | 2019-05-10 |
Family
ID=58268301
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610871781.0A Expired - Fee Related CN106487358B (zh) | 2016-09-30 | 2016-09-30 | 一种机动目标转弯跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106487358B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107368849B (zh) * | 2017-06-30 | 2020-04-28 | 西安交通大学 | 一种基于互相关熵的共用空间模式空域特征提取方法 |
CN108562290B (zh) * | 2018-07-13 | 2020-10-27 | 深圳市戴升智能科技有限公司 | 导航数据的滤波方法、装置、计算机设备及存储介质 |
CN110034746B (zh) * | 2019-04-02 | 2023-10-13 | 汕头大学 | 一种基于最大协同熵卡尔曼滤波方法 |
CN110135012B (zh) * | 2019-04-24 | 2023-12-22 | 北京启迪清云智慧能源有限公司 | 一种系统线性回归预测模型的回归系数确定方法 |
CN110532517A (zh) * | 2019-08-28 | 2019-12-03 | 北京理工大学 | 基于改进的arukf的燃气管道参数估计方法 |
CN111010145A (zh) * | 2019-12-10 | 2020-04-14 | 西南大学 | 基于范数正则化离散线性系统的滤波方法、离散线性系统 |
CN111256873B (zh) * | 2020-02-27 | 2021-05-25 | 南昌航空大学 | 基于卡尔曼滤波和迭代学习的分布式光纤温度预测方法 |
CN111756353B (zh) * | 2020-06-12 | 2024-04-16 | 杭州电子科技大学 | 一种基于非线性融合滤波的液位仪噪声优化方法 |
CN112836354B (zh) * | 2021-01-12 | 2022-08-30 | 中南大学 | 一种目标跟踪定位方法、系统、装置及可读存储介质 |
CN113032988B (zh) * | 2021-03-18 | 2024-04-12 | 杭州电子科技大学 | 基于最大相关熵的高阶扩展卡尔曼滤波器设计方法 |
CN113378469B (zh) * | 2021-06-21 | 2022-04-08 | 大庆油田有限责任公司 | 一种基于卡尔曼滤波与支持向量机的测井曲线预测方法 |
CN113449384A (zh) * | 2021-07-07 | 2021-09-28 | 中国人民解放军军事科学院国防科技创新研究院 | 一种基于中心误差熵准则扩展卡尔曼滤波的姿态确定方法 |
CN113792412B (zh) * | 2021-08-13 | 2022-10-11 | 中国人民解放军军事科学院国防科技创新研究院 | 一种基于中心误差熵准则容积卡尔曼滤波的航天器姿态确定方法 |
CN114609912B (zh) * | 2022-03-18 | 2023-10-03 | 电子科技大学 | 基于伪线性最大相关熵卡尔曼滤波的仅测角目标追踪方法 |
CN117353705A (zh) * | 2023-10-08 | 2024-01-05 | 哈尔滨工业大学 | 一种带有未知延迟概率的一步时延跟踪滤波方法及系统 |
CN117216482B (zh) * | 2023-11-07 | 2024-03-01 | 南开大学 | 借助迭代策略增强卡尔曼滤波器抗野值干扰能力的方法 |
CN117713750B (zh) * | 2023-12-14 | 2024-05-17 | 河海大学 | 一种基于分数次幂的一致性卡尔曼滤波状态估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217175A (zh) * | 2013-04-10 | 2013-07-24 | 哈尔滨工程大学 | 一种自适应容积卡尔曼滤波方法 |
CN104168005A (zh) * | 2014-08-25 | 2014-11-26 | 北京理工大学 | 带有未知观测噪声协方差阵递推估计的卡尔曼滤波方法 |
CN104408744A (zh) * | 2014-11-17 | 2015-03-11 | 电子科技大学 | 一种用于目标跟踪的强跟踪容积卡尔曼滤波方法 |
CN105737828A (zh) * | 2016-05-09 | 2016-07-06 | 郑州航空工业管理学院 | 一种基于强跟踪的相关熵扩展卡尔曼滤波的组合导航方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8483432B2 (en) * | 2009-12-23 | 2013-07-09 | General Electric Company | Methods for automatic segmentation and temporal tracking |
-
2016
- 2016-09-30 CN CN201610871781.0A patent/CN106487358B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217175A (zh) * | 2013-04-10 | 2013-07-24 | 哈尔滨工程大学 | 一种自适应容积卡尔曼滤波方法 |
CN104168005A (zh) * | 2014-08-25 | 2014-11-26 | 北京理工大学 | 带有未知观测噪声协方差阵递推估计的卡尔曼滤波方法 |
CN104408744A (zh) * | 2014-11-17 | 2015-03-11 | 电子科技大学 | 一种用于目标跟踪的强跟踪容积卡尔曼滤波方法 |
CN105737828A (zh) * | 2016-05-09 | 2016-07-06 | 郑州航空工业管理学院 | 一种基于强跟踪的相关熵扩展卡尔曼滤波的组合导航方法 |
Non-Patent Citations (2)
Title |
---|
几种面向弹道目标跟踪算法的性能评估;王品 等;《深圳大学学报理工版》;20120930;第29卷(第5期);正文第392-398页 |
基于高斯和均方根容积卡尔曼滤波的姿态角辅助目标跟踪算法;单甘霖 等;《电子与信息学报》;20140731;第36卷(第7期);正文第1579-1584页 |
Also Published As
Publication number | Publication date |
---|---|
CN106487358A (zh) | 2017-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106487358B (zh) | 一种机动目标转弯跟踪方法 | |
CN105205313B (zh) | 模糊高斯和粒子滤波方法、装置及目标跟踪方法、装置 | |
CN110503071B (zh) | 基于变分贝叶斯标签多伯努利叠加模型的多目标跟踪方法 | |
Huang et al. | Robust student’st based nonlinear filter and smoother | |
CN103902812B (zh) | 一种粒子滤波方法、装置及目标跟踪方法、装置 | |
Wang et al. | Efficient structural reliability analysis based on adaptive Bayesian support vector regression | |
Zhao et al. | Design of adaptive robust square-root cubature Kalman filter with noise statistic estimator | |
CN104121907B (zh) | 一种基于平方根容积卡尔曼滤波器的飞行器姿态估计方法 | |
CN105354860B (zh) | 基于箱粒子滤波的扩展目标CBMeMBer跟踪方法 | |
CN103955892B (zh) | 一种目标跟踪方法及扩展截断无迹卡尔曼滤波方法、装置 | |
CN111127523B (zh) | 基于量测迭代更新的多传感器gmphd自适应融合方法 | |
Wang et al. | Two stage particle filter for nonlinear Bayesian estimation | |
CN111798494B (zh) | 广义相关熵准则下的机动目标鲁棒跟踪方法 | |
Miao et al. | Neural network-aided variational Bayesian adaptive cubature Kalman filtering for nonlinear state estimation | |
CN115981148A (zh) | 无人机地面移动目标跟踪方法 | |
CN112415468A (zh) | 一种基于多伯努利滤波的doa跟踪方法 | |
CN106772354B (zh) | 基于并行模糊高斯和粒子滤波的目标跟踪方法及装置 | |
CN113452349B (zh) | 一种基于贝叶斯序贯重要性积分的卡尔曼滤波方法 | |
Mukherjee et al. | Likelihood function modeling of particle filter in presence of non-stationary non-gaussian measurement noise | |
Hu et al. | IMM fusion estimation with multiple asynchronous sensors | |
He et al. | An EM algorithm for target tracking with an unknown correlation coefficient of measurement noise | |
Huang et al. | A robust and efficient system identification method for a state-space model with heavy-tailed process and measurement noises | |
CN112379350A (zh) | 智能车辆毫米波雷达多目标跟踪方法、装置及设备 | |
CN109188422B (zh) | 一种基于lu分解的卡尔曼滤波目标跟踪方法 | |
CN109375160B (zh) | 纯方位无源定位中一种测角误差估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190510 Termination date: 20210930 |
|
CF01 | Termination of patent right due to non-payment of annual fee |