CN105353353B - 多重搜索粒子概率假设密度滤波的多目标跟踪方法 - Google Patents
多重搜索粒子概率假设密度滤波的多目标跟踪方法 Download PDFInfo
- Publication number
- CN105353353B CN105353353B CN201510791334.XA CN201510791334A CN105353353B CN 105353353 B CN105353353 B CN 105353353B CN 201510791334 A CN201510791334 A CN 201510791334A CN 105353353 B CN105353353 B CN 105353353B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- msubsup
- particle
- 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
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
-
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
- G01S13/56—Discriminating between fixed and moving objects or between objects moving at different speeds for presence detection
-
- 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
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)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种多重搜索粒子概率假设密度滤波的多目标跟踪方法,属于雷达数据处理领域。基于粒子概率假设密度滤波的多目标跟踪方法存在一个明显的缺陷,即当目标出现漏检时,重采样会造成粒子多样性的迅速退化,进而造成目标丢失的现象,因此该算法难以适应目标检测概率较低时的多目标跟踪。本发明提出的多重搜索粒子概率假设密度滤波即立足于解决此类问题。本发明具有结构简单,计算快速,易于硬件实现,同时克服了基于一般的粒子概率假设密度滤波方法应用的局限性,对非线性非高斯系统具有较强的适应性,因此具有较强的工程应用价值和推广前景。
Description
技术领域
本发明涉及一种雷达数据处理方法,特别是涉及一种低检测概率情况下雷达对多目标的跟踪方法。
背景技术
粒子概率假设密度滤波(Particle probability hypothesis density filter,PPHDF)是一种密集杂波环境下对多目标进行跟踪的有效方法。PPHDF通过将量测和目标状态建模为随机集,可以非常方便的从个数时变的量测中估计出个数时变且未知的目标状态,即可同时对目标个数和目标状态进行估计;同时,PPHDF可以避免目标和量测之间的关联问题,极大的降低了多目标跟踪算法的复杂性和计算量。因此,PPHDF在多目标跟踪领域得到了广泛的关注和研究。该方法主要通过以下步骤实现:
(1)初始化,得到初始粒子集;
(2)对已存在的粒子集进行一步预测得到预测粒子集,并生成用于搜索新目标的搜索粒子集,将预测粒子集和搜索粒子集合成新的预测粒子集;
(3)利用新量测对预测粒子集粒子权重进行更新;
(4)对权重更新后的粒子集进行重采样,并得到目标个数和各目标状态估计。
基于PPHDF的多目标跟踪方法存在一个明显的缺陷,即当目标出现漏检时,重采样会造成粒子多样性的迅速退化,进而造成目标丢失的现象,因此该算法难以适应目标检测概率较低时的多目标跟踪。
发明内容
本发明的目的是提出一种多重搜索粒子概率假设密度滤波(MS-PPHDF)的多目标跟踪方法,解决一般的PPHDF方法在目标检测概率较低的情况下容易出现目标丢失的问题。
本发明提出的MS-PPHDF方法的技术方案包括以下步骤:
步骤1:变量初始化
K是雷达关机时刻;
T是雷达扫描周期;
L0为代表1个目标的粒子数;
D0为目标出现的初始分布;
Jk为搜索新目标的粒子数;
Sk表示搜索1个消失目标的粒子数;
Lk为k时刻滤波器采用的粒子总数;
γk为平均目标出现概率;
PD为目标检测概率;
λk为平均每帧的杂波个数;
Fk为状态转移矩阵;
Gk为过程噪声分布矩阵;
Qk为过程噪声协方差;
Rk为量测噪声协方差;
步骤2:令k=0,初始化粒子集
对任意p∈{1,2,…,L0},从初始分布D0中采样粒子并赋予该粒子权重得到初始粒子集其中表示粒子代表的目标状态,包含了目标的位置和速度信息;
步骤3:令k=k+1,获得k时刻的雷达量测
将雷达接收到的信号进行A/D变换,得到k时刻的雷达量测集送雷达数据处理计算机,其中表示k时刻雷达得到的第q个量测,包含了目标的距离量测多普勒速度量测以及方位量测等信息,而Mk则表示k时刻雷达得到的量测个数;
步骤4:预测
(1)若k≤2,令Ik=0,直接转(4),否则定义集合
其中,表示k-2时刻存在而k-1时刻消失的第n个目标的状态,Ndis,k-1表示消失的目标数;
(2)若令Ik=0,直接转(4),否则令Ik=Ndis,k-1Sk,对任意根据状态转移方程
进行一步预测,其中
(3)对任意n∈{1,2,…,Ndis,k-1}和任意p∈{Lk-1+(n-1)Sk+1,…,Lk-1+nSk},根据
采样粒子并赋予该粒子权重其中
为过程噪声分布矩阵,vk为过程噪声,其噪声协方差为Qk;
(4)对任意p∈{Lk-1+Ik+1,…,Lk-1+Ik+Jk},根据初始分布D0采样“新生”粒子并赋予该粒子权重
(5)对任意p∈{1,2,…,Lk-1},根据
采样粒子并赋予该粒子权重
步骤5:更新
对任意p∈{1,2,…,Lk-1+Ik+Jk},利用量测集Zk对粒子权重进行更新
其中
gk(z|x)为量测似然函数;
步骤6:重采样
(1)计算所有粒子的权重和
(2)对粒子集进行重采样,得到新的粒子集其中为k时刻估计的目标数,Round(x)表示取与x最接近的整数;步骤7:目标状态估计
若直接转步骤8,否则采用K-均值聚类分析的方法将粒子集划分为个类其中Lk,n表示第n个类包含的粒子数,满足则估计的目标状态为
步骤8:重复步骤3~步骤8,直至雷达关机。
和背景技术相比,本发明的有益效果说明:
本发明提出的MS-PPHDF多目标跟踪方法,可以解决一般的PPHDF方法在目标检测概率较低的情况下容易出现目标丢失的问题,提高了PPHDF算法的适应范围。
附图说明
附图1是本发明提出的MS-PPHDF多目标跟踪方法的整体流程图;
附图2是本发明实施例中MS-PPHDF方法跟踪多目标的效果展示,附图中星号“*”表示目标的真实位置,圆圈“o”表示估计的目标状态;
附图3是本发明实施例中一般的PPHDF方法跟踪多目标的效果展示,附图中星号“*”表示目标的真实位置,圆圈“o”表示估计的目标状态。
具体实施方式
下面结合附图对本发明提出的MS-PPHDF方法进行详细描述。
不失一般性,设置一个二维的仿真场景,监测区域S=[-60km,60km]×[-60km,60km],总仿真时间K=50s。假设目标可在监测区域内随机出现和消失,平均目标出现概率γk=0.2,目标出现的初始分布D0服从均值为x0和协方差为Qb的正态分布,这里取x0=[30km0.2km/s 30km-0.1km/s]T和Qb=diag([1km 0.5km/s 1km 0.5km/s]),x方向和y方向过程噪声的标准差均为0.01km,目标持续存在的概率为ek|k-1=0.95且与目标状态无关。雷达位于点(0km,-10km),检测概率为PD=0.75,平均每帧的杂波个数为λk=4,且在整个监测区域内均匀分布,距离量测误差标准差、多普勒速度量测误差标准差和方位角量测误差标准差分别为0.2km、0.04km/s和0.0087rad,雷达采样间隔T=1s;代表1个目标的粒子数L0=3000,搜索新目标的粒子数Jk=4000,搜索1个消失目标的粒子数Sk=1000。
其步骤如附图1所示。
(1)根据以上仿真条件进行变量初始化
由以上仿真条件可知
雷达量测方程为
其中Wk为相互独立的零均值高斯白噪声,其协方差为
(2)按发明内容部分步骤2所述的方法进行粒子集初始化,得到初始粒子集;
(3)按发明内容部分步骤3所述的方法获得当前时刻雷达量测;
(4)按发明内容部分步骤4所述的方法对粒子集进行一步预测,得到预测的粒子集;
(5)按发明内容部分步骤5所述的方法对预测的粒子集的权重进行更新;
(6)按发明内容部分步骤6所述的方法对权重更新后的粒子集进行重采样;
(7)按发明内容部分步骤7所述的方法得到目标的状态估计;
(8)循环执行发明内容部分步骤3~步骤8,直至雷达关机。
实施例条件中,在目标检测概率PD=0.75的情况下,本发明提出的MS-PPHDF方法仍能实现杂波环境下对多目标的有效跟踪(见附图2),而一般的PPHDF方法在目标出现漏检时出现了目标的丢失的情况(见附图3),因此本发明方法可以有效的改善一般PPHDF方法在目标检测概率较低情况下的多目标跟踪性能(附图2和附图3对比)。
Claims (1)
1.多重搜索粒子概率假设密度滤波的多目标跟踪方法,其特征在于,包括以下步骤:
步骤1,变量初始化:
K是雷达关机时刻;
T是雷达扫描周期;
L0为代表1个目标的粒子数;
D0为目标出现的初始分布;
Jk为搜索新目标的粒子数;
Sk表示搜索1个消失目标的粒子数;
Lk为k时刻滤波器采用的粒子总数;
γk为平均目标出现概率;
PD为目标检测概率;
λk为平均每帧的杂波个数;
Fk为状态转移矩阵;
Gk为过程噪声分布矩阵;
Qk为过程噪声协方差;
Rk为量测噪声协方差;
步骤2,令k=0,初始化粒子集:
对任意p∈{1,2,…,L0},从初始分布D0中采样粒子并赋予该粒子权重得到初始粒子集其中表示粒子代表的目标状态,包含了目标的位置和速度信息;
步骤3,令k=k+1,获得k时刻的雷达量测:
将雷达接收到的信号进行A/D变换,得到k时刻的雷达量测集送雷达数据处理计算机,其中表示k时刻雷达得到的第q个量测,包含了目标的距离量测多普勒速度量测以及方位量测信息,而Mk则表示k时刻雷达得到的量测个数;
步骤4,预测:
(1)若k≤2,令Ik=0,直接转(4),否则定义集合
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>d</mi>
<mi>i</mi>
<mi>s</mi>
<mo>,</mo>
<mi>k</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<mo>{</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>d</mi>
<mi>i</mi>
<mi>s</mi>
<mo>,</mo>
<mi>k</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo>|</mo>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>d</mi>
<mi>i</mi>
<mi>s</mi>
<mo>,</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>}</mo>
</mrow>
其中,表示k-2时刻存在而k-1时刻消失的第n个目标的状态,Ndis,k-1表示消失的目标数;
(2)若令Ik=0,直接转(4),否则令Ik=Ndis,k-1Sk,对任意根据状态转移方程
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
<mo>,</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>x</mi>
<mrow>
<mi>d</mi>
<mi>i</mi>
<mi>s</mi>
<mo>,</mo>
<mi>k</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
<mi>n</mi>
</msubsup>
</mrow>
进行一步预测,其中
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>T</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>T</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
(3)对任意n∈{1,2,...,Ndis,k-1}和任意p∈{Lk-1+(n-1)Sk+1,Lk-1+(n-1)Sk+2,…,Lk-1+nSk},根据
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>F</mi>
<mi>k</mi>
</msub>
<msubsup>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
<mo>,</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>G</mi>
<mi>k</mi>
</msub>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
</mrow>
采样粒子并赋予该粒子权重其中
<mrow>
<msub>
<mi>G</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>T</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mi>T</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
为过程噪声分布矩阵,vk为过程噪声,其噪声协方差为Qk;
(4)对任意p∈{Lk-1+Ik+1,Lk-1+Ik+2,…,Lk-1+Ik+Jk},根据初始分布D0采样“新生”粒子并赋予该粒子权重
(5)对任意p∈{1,2,…,Lk-1},根据
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>F</mi>
<mi>k</mi>
</msub>
<msubsup>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>G</mi>
<mi>k</mi>
</msub>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
</mrow>
采样粒子并赋予该粒子权重
步骤5,更新:
对任意p∈{1,2,…,Lk-1+Ik+Jk},利用量测集Zk对粒子权重进行更新
<mrow>
<msubsup>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>z</mi>
<mo>&Element;</mo>
<msub>
<mi>Z</mi>
<mi>k</mi>
</msub>
</mrow>
</munder>
<mfrac>
<mrow>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>z</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&lambda;</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<msubsup>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
</mrow>
其中
<mrow>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>z</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
<msub>
<mi>g</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>|</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>I</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>J</mi>
<mi>k</mi>
</msub>
</mrow>
</munderover>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>z</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
<msubsup>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
</mrow>
gk(z|x)为量测似然函数;
步骤6,重采样:
(1)计算所有粒子的权重和
<mrow>
<msub>
<mover>
<mi>N</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>I</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>J</mi>
<mi>k</mi>
</msub>
</mrow>
</msubsup>
<msubsup>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mi>p</mi>
</msubsup>
</mrow>
(2)对粒子集进行重采样,得到新的粒子集其中为k时刻估计的目标数,Round(x)表示取与x最接近的整数;步骤7,目标状态估计:
若直接转步骤8,否则采用K-均值聚类分析的方法将粒子集划分为个类其中Lk,n表示第n个类包含的粒子数,满足则估计的目标状态为
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>L</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
</munderover>
<msubsup>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
<mi>p</mi>
</msubsup>
<mo>,</mo>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mover>
<mi>N</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
</msub>
</mrow>
步骤8,重复步骤3~步骤8,直至雷达关机。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510791334.XA CN105353353B (zh) | 2015-11-17 | 2015-11-17 | 多重搜索粒子概率假设密度滤波的多目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510791334.XA CN105353353B (zh) | 2015-11-17 | 2015-11-17 | 多重搜索粒子概率假设密度滤波的多目标跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105353353A CN105353353A (zh) | 2016-02-24 |
CN105353353B true CN105353353B (zh) | 2017-08-18 |
Family
ID=55329354
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510791334.XA Active CN105353353B (zh) | 2015-11-17 | 2015-11-17 | 多重搜索粒子概率假设密度滤波的多目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105353353B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107037423A (zh) * | 2016-11-09 | 2017-08-11 | 谭顺成 | 结合幅值信息的phd滤波多目标跟踪方法 |
CN106772357B (zh) * | 2016-11-23 | 2019-01-29 | 谭顺成 | 信噪比未知条件下的ai-phd滤波器多目标跟踪方法 |
CN107656265B (zh) * | 2017-09-19 | 2021-03-30 | 电子科技大学 | 针对多帧检测前跟踪短航迹的粒子滤波融合方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103310115B (zh) * | 2013-06-27 | 2015-12-23 | 西安电子科技大学 | 一种多目标跟踪的杂波估计方法 |
CN104019816A (zh) * | 2014-05-30 | 2014-09-03 | 西北工业大学 | 基于全局时空相关的概率假设密度滤波器航迹提取方法 |
CN104518756B (zh) * | 2014-12-16 | 2018-01-16 | 中国人民解放军海军航空工程学院 | 基于概率假设密度phd的修正粒子滤波器 |
-
2015
- 2015-11-17 CN CN201510791334.XA patent/CN105353353B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105353353A (zh) | 2016-02-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101975575B (zh) | 基于粒子滤波的被动传感器多目标跟踪方法 | |
CN105301584B (zh) | 同时解距离模糊的ipphdf机动多目标跟踪方法 | |
CN109633590B (zh) | 基于gp-vsmm-jpda的扩展目标跟踪方法 | |
CN103885057B (zh) | 自适应变滑窗多目标跟踪方法 | |
CN104931934B (zh) | 一种基于pam聚类分析的雷达点迹凝聚方法 | |
CN105372659A (zh) | 道路交通监控多目标探测跟踪方法及跟踪系统 | |
CN105205313A (zh) | 模糊高斯和粒子滤波方法、装置及目标跟踪方法、装置 | |
CN106054169A (zh) | 基于跟踪信息的多站雷达信号融合检测方法 | |
CN104715154B (zh) | 基于kmdl准则判据的核k‑均值航迹关联方法 | |
CN104730537A (zh) | 基于多尺度模型的红外/激光雷达数据融合目标跟踪方法 | |
CN105093198A (zh) | 一种分布式外辐射源雷达组网探测的航迹融合方法 | |
CN109031229B (zh) | 一种杂波环境下目标跟踪的概率假设密度方法 | |
CN105353353B (zh) | 多重搜索粒子概率假设密度滤波的多目标跟踪方法 | |
CN107703496A (zh) | 一种交互式多模伯努利滤波的机动弱目标检测前跟踪方法 | |
Kim et al. | DeepTC: ConvLSTM network for trajectory prediction of tropical cyclone using spatiotemporal atmospheric simulation data | |
CN105549004B (zh) | 解距离测量模糊的impm‑pphdf方法 | |
Yang et al. | Aircraft tracking based on fully conventional network and Kalman filter | |
CN107037423A (zh) | 结合幅值信息的phd滤波多目标跟踪方法 | |
CN106019253A (zh) | 基于箱粒子cphd的多扩展目标跟踪方法 | |
CN106707272A (zh) | 一种基于随机集理论的多目标跟踪方法 | |
CN105353352B (zh) | 改进搜索策略的mm‑pphdf机动多目标跟踪方法 | |
CN104050641A (zh) | 基于形状方位描述符的集中式多传感器编队目标粒子滤波算法 | |
CN116047495B (zh) | 一种用于三坐标雷达的状态变换融合滤波跟踪方法 | |
Ebert et al. | Deep radar sensor models for accurate and robust object tracking | |
CN106772357A (zh) | 信噪比未知条件下的ai‑phd滤波器 |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20190627 Address after: 264001 scientific research office, two Road 188, Zhifu District, Yantai, Shandong Patentee after: Naval Aviation University of PLA Address before: 264001 scientific research department, 188 road two, Zhifu District, Yantai, Shandong Patentee before: Naval Aeronautical Engineering Institute PLA |