CN103324835B - 概率假设密度滤波器目标信息的保持方法及信息保持*** - Google Patents
概率假设密度滤波器目标信息的保持方法及信息保持*** Download PDFInfo
- Publication number
- CN103324835B CN103324835B CN201310208275.XA CN201310208275A CN103324835B CN 103324835 B CN103324835 B CN 103324835B CN 201310208275 A CN201310208275 A CN 201310208275A CN 103324835 B CN103324835 B CN 103324835B
- Authority
- CN
- China
- Prior art keywords
- gaussian
- moment
- target
- posterior
- term
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000005259 measurement Methods 0.000 claims abstract description 45
- 230000004083 survival effect Effects 0.000 claims abstract description 22
- 239000000203 mixture Substances 0.000 claims abstract description 7
- 238000000605 extraction Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 45
- 238000001514 detection method Methods 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 5
- 239000002245 particle Substances 0.000 claims description 2
- 230000014759 maintenance of location Effects 0.000 claims 2
- 239000000284 extract Substances 0.000 claims 1
- 238000004321 preservation Methods 0.000 claims 1
- 230000004927 fusion Effects 0.000 abstract description 2
- 238000007689 inspection Methods 0.000 abstract 2
- 230000008901 benefit Effects 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明适用于多传感器信息融合领域,提供了概率假设密度滤波器目标信息的保持方法,包括:步骤1、根据前一时刻的后验矩和高斯项预测当前时刻的后验矩和高斯项;步骤2、根据预测的当前时刻后验矩和高斯项以及当前时刻的测量集对后验矩和高斯项进行更新;步骤3、对更新后的高斯项裁减与合并;步骤4、根据裁减与合并后的高斯项,提取权重的高斯项作为滤波器的输出,相应高斯项中的均值和方差分别为存活目标的状态估计和误差估计。所述方法通过修正概率假设密度滤波器更新方程将漏检目标的信息驻留在后验的更新矩中,避免了漏检目标的信息丢失,提高了目标数估计和目标状态提取的有效性,进而提高了高斯混合概率假设密度滤波器的多目标跟踪能力。
Description
技术领域
本发明属于多传感器信息融合技术领域,尤其涉及概率假设密度滤波器目标信息的保持方法及信息保持***。
背景技术
在存在虚警、漏检和目标数未知的情况下,Mahler提出的概率假设密度滤波器是解决目标检测和跟踪的新方法。概率假设密度滤波器避免了观测值与状态值间的直接关联,其最大优点是能从后验矩中估计出目标数。为解决概率假设密度滤波器的预测与更新方程中积分运算难以处理的问题,Vo等提出了粒子概率假设密度滤波器和高斯混合概率假设密度滤波器。目前,概率假设密度滤波器已在被动定位、被动雷达目标跟踪、视频跟踪、声纳图像中的目标跟踪和组目标跟踪等领域取得了比较广泛的应用。
然而,正如Mahler所指出的,概率假设密度滤波器是一个弱记忆力的滤波器,一旦某一时刻目标被漏检,即测量集中没有源于该目标的测量,漏检目标的信息会立刻从更新的后验矩中丢弃,从而导致漏检目标的信息丢失。目标信息的丢失会导致目标数估计的不稳定和漏检目标的状态难以从后验矩中提取出来。信息丢失问题是概率假设密度滤波器中需要探索和解决的一个关键技术问题。
发明内容
本发明所要解决的技术问题在于提供概率假设密度滤波器目标信息的保持方法及信息保持***,旨在解决目标信息丢失导致目标数估计的不稳定和漏检目标的状态难以从后验矩中提取出来的问题。
本发明是这样实现的,概率假设密度滤波器目标信息的保持方法,包括以下步骤:
步骤1、根据前一时刻的后验矩和高斯项以及表征目标在前一时刻是否被漏的信息标识预测当前时刻的后验矩和高斯项;
步骤2、根据预测的当前时刻的后验矩和高斯项以及当前时刻的测量集确定出当前时刻各目标被漏的信息标识;利用所述当前时刻的测量集、当前时刻目标被漏的信息标识以及预测的当前时刻的后验矩和高斯项求取当前时刻更新的高斯项;
步骤3、对更新后的高斯项进行裁减与合并,裁减与合并后的高斯项作为当前时刻的高斯项,裁减与合并后的高斯项加权和构成当前时刻的后验矩,当前时刻的高斯项和后验矩作为下一次滤波器递归的输入;
步骤4、根据裁减与合并后的高斯项,提取权重的高斯项作为滤波器的输出,相应高斯项中的均值和方差分别为存活目标的状态估计和误差估计。
进一步地,所述步骤1中,以k-1表示前一时刻,以k表示当前时刻;k-1时刻构成后验矩的高斯项集合表示为其中w、m和P分别表示高斯项的权重、均值和方差,α为表明一个目标是否被漏的标识,i为索引号,取值从1至Jk-1,Jk-1为k-1时刻高斯项的数目;
k-1时刻后验矩表示为则k时刻预测的后验矩表示为vk|k-1(x|Z(k-1))=vS,k|k-1(x|Z(k-1))+γk(x),其中,x表示目标状态,Z(k-1)表示1至k-1时刻所有测量构成的集合,γk(x)为新生目标矩,并且
和pS,k分别为目标i的幸存概率和目标的先验幸存概率,Fk-1为状态转移矩阵,Qk-1为过程噪声方差矩阵,Jγ,k、和均为已知的模型参数,分别表示为新生高斯项的数目、第j个高斯项的权重、状态的均值和方差,N表示高斯分布,x表示目标的状态;γk(x)表示新生目标矩;预测的高斯项由和组成,其中为幸存目标i的预测信息标识,为新生目标j的信息标识,j表示索引号,取值从1至Jγ,k。
进一步地,所述步骤2中,k时刻预测的后验矩表示为
k时刻构成预测后验矩的高斯项集合表示为i=1,2,…,Jk|k-1,其中Jk|k-1=Jk-1+Jγ,k,则k时刻更新后验矩vk(x|Z(k))表示为其中, 而Z(k)为1至k时刻所有测量构成的集合,Ki为滤波器增益矩阵,Hk为观测矩阵,Rk为观测噪声的方差矩阵,为目标i被漏的概率,为目标i的检测概率,为目标l的检测概率,λc(zj)为杂波密度,zj为测量集合Zk中的第j个测量,Zk为k时刻测量集合,I为单位矩阵,为观测矩阵的转置矩阵,l为索引号,取值从1至Jk|k-1;更新高斯项由和组成,其中,i=1,2,…,Jk|k-1,j=1,2,…,nk,nk为测量集Zk中测量的个数。
进一步地,所述步骤3中,k时刻构成更新后验矩的高斯项表示为i=1,2,…,Jk|k,其中Jk|k=Jk|k-1(nk+1),删除权重充分小的高斯项,即删除的高斯项,其中,τ为裁减门限;
将距离充分小的高斯项合并成一个,即将dij<U的高斯项合并成一个,其中U为合并门限,合并距离dij定义为多个高斯项的合并方法如下: 其中,L为合并高斯项上标形成的集合,上标T和b分别表示矩阵的转置和合并后高斯项的索引号,合并后高斯项标识取为合并前权重最大高斯项的标识。
进一步地,所述步骤1中,所述目标的幸存概率根据k-1时刻目标是否被漏的标识加以确定,当被漏标识α为1,表明所述目标在k-1时刻漏检,其幸存概率取为0.5;当α为0,表明所述目标在k-1时刻被检测到,其幸存概率取为先验幸存概率pS,k。
进一步地,所述步骤2中先确定目标是否被漏:当目标被漏时,所述标识α为1,并且其预测的后验矩将完全地驻留在更新的后验矩中。
本发明还提供概率假设密度滤波器目标信息的信息保持***,其中,包括:
预测模块,根据前一时刻的后验矩和高斯项以及表征目标在前一时刻是否被漏的信息标识预测当前时刻的后验矩和高斯项;
更新模块,与所述预测模块相连,根据预测的当前时刻的后验矩和高斯项以及当前时刻的测量集确定出当前时刻各目标被漏的信息标识;利用所述当前时刻的测量集、当前时刻目标被漏的信息标识以及预测的当前时刻的后验矩和高斯项求取当前时刻更新的高斯项;
裁减与合并模块,与所述更新模块相连,对更新后的高斯项进行裁减与合并,裁减与合并后的高斯项作为当前时刻的高斯项,裁减与合并后的高斯项加权和构成当前时刻的后验矩,当前时刻的高斯项和后验矩作为下一次滤波器递归的输入;
目标状态提取模块,与所述裁减与合并模块相连,根据裁减与合并后的高斯项,提取权重的高斯项作为滤波器的输出,相应高斯项中的均值和方差分别为存活目标的状态估计和误差估计。
进一步地,在所述预测模块中,以k-1表示前一时刻,以k表示当前时刻;k-1时刻构成后验矩的高斯项集合表示为其中w、m和P分别表示高斯项的权重、均值和方差,α为表明一个目标是否被漏的标识,i为索引号,取值从1至Jk-1,Jk-1为k-1时刻高斯项的数目;
k-1时刻后验矩表示为则k时刻预测的后验矩表示为vk|k-1(x|Z(k-1))=vS,k|k-1(x|Z(k-1))+γk(x),其中,x表示目标状态,Z(k-1)表示1至k-1时刻所有测量构成的集合,γk(x)为新生目标矩,并且
和pS,k分别为目标i的幸存概率和目标的先验幸存概率,Fk-1为状态转移矩阵,Qk-1为过程噪声方差矩阵,Jγ,k、和均为已知的模型参数,分别表示为新生高斯项的数目、第j个高斯项的权重、状态的均值和方差,N表示高斯分布,x表示目标的状态;γk(x)表示新生目标矩;
预测的高斯项由i=1,2,…,Jk-1和j=1,2,…,Jγ,k组成,其中为幸存目标i的预测信息标识, 为新生目标j的信息标识。
进一步地,在所述更新模块中,k时刻预测的后验矩表示为
k时刻构成预测后验矩的高斯项集合表示为i=1,2,…,Jk|k-1,其中Jk|k-1=Jk-1+Jγ,k,则k时刻更新后验矩vk(x|Z(k))表示为
其中, 而Z(k)为1至k时刻所有测量构成的集合,Ki为滤波器增益矩阵,Hk为观测矩阵,Rk为观测噪声的方差矩阵,为目标i被漏的概率,为目标i的检测概率,为目标l的检测概率,λc(zj)为杂波密度,zj为测量集合Zk中的第j个测量,Zk为k时刻测量集合,I为单位矩阵,为观测矩阵的转置矩阵,l为索引号,取值从1至Jk|k-1;更新高斯项由和组成,其中,i=1,2,…,Jk|k-1,j=1,2,…,nk,nk为测量集Zk中测量的个数。
进一步地,在所述裁减与合并模块中,k时刻构成更新后验矩的高斯项表示为i=1,2,…,Jk|k,其中Jk|k=Jk|k-1(nk+1),删除权重充分小的高斯项,即删除的高斯项,其中,τ为裁减门限;
将距离充分小的高斯项合并成一个,即将dij<U的高斯项合并成一个,其中U为合并门限,合并距离dij定义为多个高斯项的合并方法如下: 其中,L为合并高斯项上标形成的集合,上标T和b分别表示矩阵的转置和合并后高斯项的索引号,合并后高斯项标识取为合并前权重最大高斯项的标识。
与现有概率假设密度滤波器相比,本发明的优点是:本发明利用预测的后验矩和当前时刻的测量集确定存活的真目标是否被漏检了,如果目标确实被漏检,能将漏检目标的信息驻留在后验矩中,从而避免了漏检目标的信息丢失,进而提高了目标数估计的稳定性和目标状态提取的可靠性。
附图说明
图1是本发明提供的概率假设密度滤波器中目标信息的保持方法的流程图;
图2是检测率为0.85时本发明与现有方法得到的平均目标数;
图3是检测率为0.85时本发明与现有方法的目标数均方根误差;
图4是检测率为0.85时本发明与现有方法的平均OSPA距离。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明依据后验的预测矩和当前时刻的测量集判定存活的真目标在当前时刻是否被漏检,如果目标被漏检,通过修正概率假设密度滤波器更新方程将漏检目标的信息驻留在后验的更新矩中,以使漏检目标的状态能被滤波器提取出来作为其输出。
如图1所示,本发明提供的概率假设密度滤波器中目标信息的保持方法包括以下步骤:
步骤1、根据前一时刻的后验矩和高斯项以及表征目标在前一时刻是否被漏的信息标识预测当前时刻的后验矩和高斯项。
以k-1表示前一时刻,以k表示当前时刻;k-1时刻构成后验矩的高斯项集合表示为其中w、m和P分别表示高斯项的权重、均值和方差,α为表明一个目标是否被漏的标识,Jk-1为k-1时刻高斯项的数目,i为索引号,取值从1至Jk-1。所述后验矩的高斯项集合由多个高斯项构成。
k-1时刻的后验矩表示为则k时刻预测的后验矩表示为vk|k-1(x|Z(k-1))=vS,k|k-1(x|Z(k-1))+γk(x),其中,x表示目标状态,Z(k-1)表示1至k-1时刻所有测量构成的集合,γk(x)为新生目标矩,并且 和pS,k分别为目标i的幸存概率和目标的先验幸存概率,Fk-1为状态转移矩阵,Qk-1为过程噪声方差矩阵,Jγ,k、和均为已知的模型参数,分别表示为新生高斯项的数目、第j个高斯项的权重、状态的均值和方差,N表示高斯分布,x表示目标的状态,γk(x)表示新生目标矩。所述后验矩由所述高斯项构成,即所述后验矩为高斯项的加权和。
预测的高斯项由和组成,其中为幸存目标i的预测信息标识,为新生目标j的信息标识,j表示索引号,取值从1至Jγ,k。
步骤2、根据预测的当前时刻的后验矩和高斯项以及当前时刻的测量集确定出当前时刻各目标被漏的信息标识;利用所述当前时刻的测量集、当前时刻目标被漏的信息标识以及预测的当前时刻的后验矩和高斯项求取当前时刻更新的高斯项。
根据所述步骤1可得,k时刻预测的后验矩表示为
k时刻构成预测后验矩的高斯项集合表示为i=1,2,…,Jk|k-1,其中Jk|k-1=Jk-1+Jγ,k。
根据预测的后验矩和高斯项可得,k时刻更新后验矩vk(x|Z(k))表示为其中, 而Z(k)为1至k时刻所有测量构成的集合,Ki为滤波器增益矩阵,Hk为观测矩阵,Rk为观测噪声的方差矩阵,为目标i被漏的概率,为目标i的检测概率,为目标l的检测概率,λc(zj)为杂波密度,zj为测量集合Zk中的第j个测量,Zk为k时刻测量集合,I为单位矩阵,为观测矩阵的转置矩阵,l为索引号,取值从1至Jk|k-1;更新高斯项由和组成,其中,i=1,2,…,Jk|k-1,j=1,2,…,nk,nk为测量集Zk中测量的个数。
步骤3、对更新后的高斯项进行裁减与合并,裁减与合并后的高斯项作为当前时刻的高斯项,裁减与合并后的高斯项加权和构成当前时刻的后验矩。
根据步骤2可得,k时刻构成更新后验矩的高斯项表示为i=1,2,…,Jk|k,其中Jk|k=Jk|k-1(nk+1)。删除权重充分小的高斯项,即删除的高斯项,其中τ为裁减门限。
将距离充分小的高斯项合并成一个,即将dij<U的高斯项合并成一个,其中U为合并门限。合并距离dij定义为多个高斯项的合并方法如下: 其中,L为合并高斯项上标形成的集合,上标T和b分别表示矩阵的转置和合并后高斯项的索引号,合并后高斯项标识取为合并前权重最大高斯项的标识。
步骤4、根据裁减与合并后的高斯项,假设裁减与合并后余下的高斯项为i=1,2,…,Jk,提取权重的高斯项作为滤波器的输出,相应高斯项中的均值和方差分别为存活目标的状态估计和误差估计。裁减与合并步骤后,余下的高斯项为i=1,2,…,Jk,作为下一次滤波器递归的输入。
进一步地,在所述步骤1中,所述目标的幸存概率根据k-1时刻目标是否被漏的标识加以确定,当被漏标识α为1,表明所述目标在k-1时刻漏检,其幸存概率取为0.5;当α为0,表明所述目标在k-1时刻被检测到,其幸存概率取为先验幸存概率pS,k。
进一步地,所述步骤2中,确定目标是否被漏的所述标识α是由预测的后验矩和k时刻的测量集确定,当目标被漏时,所述标识α为1,并且其预测的后验矩将完全地驻留在更新的后验矩中。
进一步地,在所述步骤3和4中,所述目标是否被漏的标识α随同所述权重w、均值m和方差P一起传送。
本发明所述的概率假设密度滤波器中目标信息的保持方法通过确定存活的目标是否被漏检将漏检目标的信息在更新的后验矩中,从而提高了概率假设密度滤波器的记忆力,避免了漏检目标的信息丢失。在杂波密度λc=5×10-6m-2、目标检测概率为0.85、既有新目标出现和已存在目标消失的情况下,本发明与现有方法对10批目标跟踪时100次MonteCarlo实验得到的平均目标数、目标数均方根误差和平均OSPA(Optimal SubpatternAssignment,最优亚模式分配)距离分别如图2、图3和图4所示。从图2和图3中可看出,与现有方法相比,本发明能获得更加精确和稳定的目标数估计,平均目标数与真实目标数更为接近,目标数估计的均方根误差更小。从图4中可看出,本发明的多目标跟踪精度明显高于现有方法,其OSPA距离比现有方法得到的OSPA距离要小得多。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.概率假设密度滤波器目标信息的保持方法,其特征在于,包括以下步骤:
步骤1、根据前一时刻的后验矩和高斯项以及表征目标在前一时刻是否被漏的信息标识预测当前时刻的后验矩和高斯项;
以k-1表示前一时刻,以k表示当前时刻,所述前一时刻的高斯项集合表示为其中w、m和P分别表示高斯项的权重、均值和方差,α为表明一个目标是否被漏的信息标识,i为索引号,取值从1至Jk-1,Jk-1为k-1时刻高斯项的数目;
k-1时刻后验矩表示为则k时刻预测的后验矩表示为vk|k-1(x|Z(k-1))=vS,k|k-1(x|Z(k-1))+γk(x),其中,x表示目标状态,Z(k -1)表示1至k-1时刻所有测量构成的集合,γk(x)为新生目标矩,并且
和pS,k分别为目标i的幸存概率和目标的先验幸存概率,Fk-1为状态转移矩阵,Qk-1为过程噪声方差矩阵,Jγ,k、和均为已知的模型参数,分别表示为新生高斯项的数目、第j个高斯项的权重、状态的均值和方差,N表示高斯分布,x表示目标的状态,γk(x)表示新生目标矩;
预测的高斯项由和组成,其中 为幸存目标i的预测信息标识, 为新生目标j的信息标识,j表示索引号,取值从1至Jγ,k;
步骤2、根据预测的当前时刻的后验矩和高斯项以及当前时刻的测量集确定出当前时刻各目标被漏的信息标识;利用所述当前时刻的测量集、当前时刻目标被漏的信息标识以及预测的当前时刻的后验矩和高斯项求取当前时刻更新的高斯项;
步骤3、对更新后的高斯项进行裁减与合并,裁减与合并后的高斯项为当前时刻的高斯项,裁减与合并后的高斯项加权和构成当前时刻的后验矩,当前时刻的高斯项和后验矩作为下一次滤波器递归的输入;
步骤4、根据当前时刻的高斯项,提取权重的高斯项作为滤波器的输出,相应高斯项中的均值和方差分别为存活目标的状态估计和误差估计。
2.根据权利要求1所述的保持方法,其特征在于,所述步骤2中,k时刻预测的后验矩表示为k时刻构成预测后验矩的高斯项集合表示为i=1,2,…,Jk|k-1,则k时刻更新后验矩vk(x|Z(k))表示为
其中,
而Z(k)为1至k时刻所有测量构成的集合,Ki为滤波器增益矩阵,Hk为观测矩阵,Rk为观测噪声的方差矩阵,为目标i被漏的概率,为目标i的检测概率,为目标l的检测概率,λc(zj)为杂波密度,zj为测量集合Zk中的第j个测量,Zk为k时刻测量集合,I为单位矩阵,为观测矩阵的转置矩阵,l为索引号,取值从1至Jk|k-1;更新高斯项由和组成,其中,i=1,2,…,Jk|k-1,j=1,2,…,nk,nk为测量集Zk中测量的个数。
3.根据权利要求2所述的保持方法,其特征在于,所述步骤3中,k时刻构成更新后验矩的高斯项表示为i=1,2,…,Jk|k,删除权重充分小的高斯项,即删除的高斯项,其中,τ为裁减门限;
将距离充分小的高斯项合并成一个,即将dij<U的高斯项合并成一个,其中U为合并门限,合并距离dij定义为多个高斯项的合并方法如下: 其中,L为合并高斯项上标形成的集合,上标T和b分别表示矩阵的转置和合并后高斯项的索引号,合并后高斯项标识取为合并前权重最大高斯项的标识。
4.根据权利要求1所述的保持方法,其特征在于,所述步骤1中,所述目标的幸存概率根据k-1时刻目标是否被漏的标识加以确定,当被漏标识α为1,表明所述目标在k-1时刻漏检,其幸存概率取为0.5;当α为0,表明所述目标在k-1时刻被检测到,其幸存概率取为先验幸存概率pS,k。
5.根据权利要求2所述的保持方法,其特征在于,所述步骤2中先确定目标是否被漏:当目标被漏时,所述标识α为1,并且其预测的后验矩将完全地驻留在更新的后验矩中。
6.概率假设密度滤波器目标信息的信息保持***,其特征在于,包括:预测模块、更新模块、裁减与合并模块和目标状态提取模块;
所述预测模块根据前一时刻的后验矩和高斯项以及表征目标在前一时刻是否被漏的信息标识预测当前时刻的后验矩和高斯项;
以k-1表示前一时刻,以k表示当前时刻,所述前一时刻的高斯项集合表示为其中w、m和P分别表示高斯项的权重、均值和方差,α为表明一个目标是否被漏的信息标识,i为索引号,取值从1至Jk-1,Jk-1为k-1时刻高斯项的数目;
k-1时刻后验矩表示为则k时刻预测的后验矩表示为vk|k-1(x|Z(k-1))=vS,k|k-1(x|Z(k-1))+γk(x),其中,x表示目标状态,Z(k -1)表示1至k-1时刻所有测量构成的集合,γk(x)为新生目标矩,并且
和pS,k分别为目标i的幸存概率和目标的先验幸存概率,Fk-1为状态转移矩阵,Qk-1为过程噪声方差矩阵,Jγ,k、和均为已知的模型参数,分别表示为新生高斯项的数目、第j个高斯项的权重、状态的均值和方差,N表示高斯分布,x表示目标的状态;γk(x)表示新生目标矩;
预测的高斯项由和组成,其中 为幸存目标i的预测信息标识, 为新生目标j的信息标识,j表示索引号,取值从1至Jγ,k;
所述更新模块与所述预测模块相连,根据预测的当前时刻的后验矩和高斯项以及当前时刻的测量集确定出当前时刻各目标被漏的信息标识;利用所述当前时刻的测量集、当前时刻目标被漏的信息标识以及预测的当前时刻的后验矩和高斯项求取当前时刻更新的高斯项;
所述裁减与合并模块与所述更新模块相连,对更新后的高斯项进行裁减与合并,裁减与合并后的高斯项作为当前时刻的高斯项,裁减与合并后的高斯项加权和构成当前时刻的后验矩,当前时刻的高斯项和后验矩作为下一次滤波器递归的输入;
所述目标状态提取模块与所述裁减与合并模块相连,根据裁减与合并后的高斯项,提取权重的高斯项作为滤波器的输出,相应高斯项中的均值和方差分别为存活目标的状态估计和误差估计。
7.根据权利要求6所述的信息保持***,其特征在于,在所述更新模块中,k时刻预测的后验矩表示为k时刻构成预测后验矩的高斯项集合表示为i=1,2,…,Jk|k-1,则k时刻更新后验矩vk(x|Z(k))表示为
其中, 而Z(k)为1至k时刻所有测量构成的集合,Ki为滤波器增益矩阵,Hk为观测矩阵,Rk为观测噪声的方差矩阵,为目标i被漏的概率,为目标i的检测概率,为目标l的检测概率,λc(zj)为杂波密度,zj为测量集合Zk中的第j个测量,Zk为k时刻测量集合,I为单位矩阵,为观测矩阵的转置矩阵,l为索引号,取值从1至Jk|k-1;更新高斯项由和组成,其中,i=1,2,…,Jk|k-1,j=1,2,…,nk,nk为测量集Zk中测量的个数。
8.根据权利要求7所述的信息保持***,其特征在于,在所述裁减与合并模块中,k时刻构成更新后验矩的高斯项表示为i=1,2,…,Jk|k,删除权重充分小的高斯项,即删除的高斯项,其中,τ为裁减门限;
将距离充分小的高斯项合并成一个,即将dij<U的高斯项合并成一个,其中U为合并门限,合并距离dij定义为多个高斯项的合并方法如下: 其中,L为合并高斯项上标形成的集合,上标T和b分别表示矩阵的转置和合并后高斯项的索引号,合并后高斯项标识取为合并前权重最大高斯项的标识。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310208275.XA CN103324835B (zh) | 2013-05-30 | 2013-05-30 | 概率假设密度滤波器目标信息的保持方法及信息保持*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310208275.XA CN103324835B (zh) | 2013-05-30 | 2013-05-30 | 概率假设密度滤波器目标信息的保持方法及信息保持*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103324835A CN103324835A (zh) | 2013-09-25 |
CN103324835B true CN103324835B (zh) | 2016-09-28 |
Family
ID=49193574
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310208275.XA Expired - Fee Related CN103324835B (zh) | 2013-05-30 | 2013-05-30 | 概率假设密度滤波器目标信息的保持方法及信息保持*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103324835B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732564A (zh) * | 2018-03-13 | 2018-11-02 | 中国电子科技集团公司第二十八研究所 | 一种双雷达修正序贯高斯混合概率假设密度滤波方法 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103679753A (zh) * | 2013-12-16 | 2014-03-26 | 深圳大学 | 一种概率假设密度滤波器的轨迹标识方法及轨迹标识*** |
CN103701433A (zh) * | 2013-12-29 | 2014-04-02 | 哈尔滨理工大学 | 一种多重测量丢失情形下时变目标跟踪***的量化滤波方法 |
CN103902829B (zh) * | 2014-04-11 | 2017-02-15 | 深圳大学 | 传递边缘分布及存在概率的目标跟踪方法与目标跟踪*** |
CN104063615B (zh) * | 2014-07-03 | 2017-02-15 | 深圳大学 | 一种基于变系数α‑β滤波器的目标跟踪方法与跟踪*** |
CN104778358B (zh) * | 2015-04-09 | 2017-11-10 | 西安工程大学 | 多传感器存在监测区域部分重叠的扩展目标跟踪方法 |
CN105975772B (zh) * | 2016-05-04 | 2019-02-05 | 浙江大学 | 基于概率假设密度滤波的多目标检测前跟踪方法 |
CN108333569B (zh) * | 2018-01-19 | 2021-01-12 | 杭州电子科技大学 | 一种基于phd滤波的异步多传感器融合多目标跟踪方法 |
CN108344981B (zh) * | 2018-01-19 | 2020-08-25 | 杭州电子科技大学 | 面向杂波的多传感器异步检测tsbf多目标跟踪方法 |
CN109194305B (zh) * | 2018-08-20 | 2021-07-13 | 电子科技大学 | 基于密度估计的数字化仪均值滤波方法 |
CN111257865B (zh) * | 2020-02-07 | 2021-09-24 | 电子科技大学 | 一种基于线性伪量测模型的机动目标多帧检测跟踪方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102307041A (zh) * | 2011-04-29 | 2012-01-04 | 浙江大学 | 当前统计模型的概率假设密度粒子滤波器的设计及滤波器 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5324391B2 (ja) * | 2009-10-22 | 2013-10-23 | キヤノン株式会社 | 画像処理装置およびその制御方法 |
-
2013
- 2013-05-30 CN CN201310208275.XA patent/CN103324835B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102307041A (zh) * | 2011-04-29 | 2012-01-04 | 浙江大学 | 当前统计模型的概率假设密度粒子滤波器的设计及滤波器 |
Non-Patent Citations (4)
Title |
---|
一种非线性GM-PHD滤波新方法;王品等;《电子学报》;20120815;第40卷(第8期);全文 * |
基于FCM聚类的粒子滤波多目标跟踪算法;陈龙等;《仪器仪表学报》;20111115;第32卷(第11期);全文 * |
概率假设密度滤波算法及其在多目标跟踪中的应用;张鹤冰;《中国优秀博士论文全文数据库》;20120401;全文 * |
独立元子空间算法及其在故障检测上的应用;张沐光等;《化工学报》;20100215;第61卷(第2期);全文 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732564A (zh) * | 2018-03-13 | 2018-11-02 | 中国电子科技集团公司第二十八研究所 | 一种双雷达修正序贯高斯混合概率假设密度滤波方法 |
CN108732564B (zh) * | 2018-03-13 | 2020-04-17 | 中国电子科技集团公司第二十八研究所 | 一种双雷达修正序贯高斯混合概率假设密度滤波方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103324835A (zh) | 2013-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103324835B (zh) | 概率假设密度滤波器目标信息的保持方法及信息保持*** | |
Vo et al. | A random finite set conjugate prior and application to multi-target tracking | |
CN104076355B (zh) | 基于动态规划的强杂波环境中弱小目标检测前跟踪方法 | |
CN105761276B (zh) | 基于迭代ransac自适应新生目标强度估计的gm-phd多目标跟踪方法 | |
CN104766320B (zh) | 阈值化量测下的多伯努利滤波弱目标检测与跟踪方法 | |
US10935653B2 (en) | Multi-target tracking method and tracking system applicable to clutter environment | |
WO2012009947A1 (zh) | 特征迹驱动的多目标跟踪装置及方法 | |
CN106291533A (zh) | 一种基于amd的分布式多传感器融合算法 | |
CN108802722A (zh) | 一种基于虚拟谱的弱目标检测前跟踪方法 | |
CN103472445A (zh) | 一种针对多目标场景的检测跟踪一体化方法 | |
CN104793200B (zh) | 一种基于迭代处理的动态规划检测前跟踪方法 | |
CN102915545A (zh) | 一种基于OpenCV的视频目标跟踪算法 | |
CN109031229B (zh) | 一种杂波环境下目标跟踪的概率假设密度方法 | |
CN106934324A (zh) | 基于简化多假设算法的雷达数据关联方法 | |
Wang et al. | Multiple model particle filter track-before-detect for range ambiguous radar | |
CN110889862B (zh) | 一种网络传输攻击环境中多目标跟踪的组合测量方法 | |
CN104181524A (zh) | 一种自适应粒子数的多目标粒子滤波检测前跟踪方法 | |
CN104501812A (zh) | 基于自适应新生目标强度的滤波算法 | |
CN111830501B (zh) | Hrrp历史特征辅助的信号模糊数据关联方法及*** | |
CN110376581B (zh) | 基于高斯混合概率假设密度滤波器的显式多目标跟踪方法 | |
CN103424742B (zh) | 一种序贯处理测量数据的目标跟踪方法与目标跟踪*** | |
CN112213697B (zh) | 一种基于贝叶斯决策理论用于雷达欺骗干扰识别的特征融合方法 | |
Zhang et al. | Adaptive Gaussian mixture probability hypothesis density for tracking multiple targets | |
WO2016000487A1 (zh) | 一种基于变系数α-β滤波器的目标跟踪方法与跟踪*** | |
Du et al. | ADOS-CFAR algorithm for multibeam seafloor terrain detection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160928 Termination date: 20170530 |
|
CF01 | Termination of patent right due to non-payment of annual fee |