结合幅度信息的扩展目标随机有限集跟踪方法

柳超 孙进平 陈小龙 张志国

柳超, 孙进平, 陈小龙, 等. 结合幅度信息的扩展目标随机有限集跟踪方法[J]. 雷达学报, 待出版. doi:  10.12000/JR19071
引用本文: 柳超, 孙进平, 陈小龙, 等. 结合幅度信息的扩展目标随机有限集跟踪方法[J]. 雷达学报, 待出版. doi:  10.12000/JR19071
LIU Chao, SUN Jinping, CHEN Xiaolong, et al. Random finite set-based extended target tracking method with amplitude information[J]. Journal of Radars, in press. doi:  10.12000/JR19071
Citation: LIU Chao, SUN Jinping, CHEN Xiaolong, et al. Random finite set-based extended target tracking method with amplitude information[J]. Journal of Radars, in press. doi:  10.12000/JR19071

结合幅度信息的扩展目标随机有限集跟踪方法

doi: 10.12000/JR19071
基金项目: 国家自然科学基金资助项目(61471019, U1633122)
详细信息
    作者简介:

    柳 超(1984–),男,山东宁阳人。北京航空航天大学博士生,研究方向为雷达数据处理。E-mail: LC2016@buaa.edu.cn

    孙进平(1975–),男,甘肃天水人,北京航空航天大学教授,博士生导师,主要研究方向为高分辨率雷达信号处理,数据处理,稀疏微波成像。E-mail: sunjinping@buaa.edu.cn

    陈小龙(1985–),男,山东烟台人,海军航空大学副教授,主要研究方向为雷达动目标检测、海杂波抑制、雷达信号精细化处理等。E-mail: cxlcxl1203@163.com

    张志国(1995–),男,山东聊城人,北京航空航天大学博士生,研究方向为雷达数据处理。E-mail: zzguo2016@163.com

    通讯作者:

    孙进平 sunjinping@buaa.edu.cn

  • 中图分类号: TP391.41

Random Finite Set-based Extended Target Tracking Method With Amplitude Information

Funds: The National Natural Science Foundation of China (61471019, U1633122)
More Information
图(12)
计量
  • 文章访问数:  546
  • HTML全文浏览量:  221
  • PDF下载量:  30
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-07-25
  • 修回日期:  2019-10-26
  • 网络出版日期:  2019-11-26

结合幅度信息的扩展目标随机有限集跟踪方法

doi: 10.12000/JR19071
    基金项目:  国家自然科学基金资助项目(61471019, U1633122)
    作者简介:

    柳 超(1984–),男,山东宁阳人。北京航空航天大学博士生,研究方向为雷达数据处理。E-mail: LC2016@buaa.edu.cn

    孙进平(1975–),男,甘肃天水人,北京航空航天大学教授,博士生导师,主要研究方向为高分辨率雷达信号处理,数据处理,稀疏微波成像。E-mail: sunjinping@buaa.edu.cn

    陈小龙(1985–),男,山东烟台人,海军航空大学副教授,主要研究方向为雷达动目标检测、海杂波抑制、雷达信号精细化处理等。E-mail: cxlcxl1203@163.com

    张志国(1995–),男,山东聊城人,北京航空航天大学博士生,研究方向为雷达数据处理。E-mail: zzguo2016@163.com

    通讯作者: 孙进平 sunjinping@buaa.edu.cn
  • 中图分类号: TP391.41

摘要: 基于随机有限集的扩展目标跟踪方法通常根据量测的空间信息进行量测划分,在杂波密集环境下有可能将杂波量测划入目标单元,从而造成跟踪性能的下降。为此,该文将目标和杂波的幅度信息引入高斯逆威沙特概率假设密度(GIW-PHD)滤波器,通过计算量测子集的幅度似然寻找最优的量测划分方法。此外,计算量测单元的中心时,采用幅度加权的方法计算量测单元的质量中心,以取代目前广泛使用的几何中心,从而进一步降低杂波对滤波器的干扰。在信杂比分别为13 dB和6 dB的条件下,通过对Rayleigh杂波中Swerling 1型起伏目标的跟踪结果证明了所提方法相比高斯逆威沙特概率假设密度滤波器具有更优的势估计和状态估计性能。

English Abstract

柳超, 孙进平, 陈小龙, 等. 结合幅度信息的扩展目标随机有限集跟踪方法[J]. 雷达学报, 待出版. doi:  10.12000/JR19071
引用本文: 柳超, 孙进平, 陈小龙, 等. 结合幅度信息的扩展目标随机有限集跟踪方法[J]. 雷达学报, 待出版. doi:  10.12000/JR19071
LIU Chao, SUN Jinping, CHEN Xiaolong, et al. Random finite set-based extended target tracking method with amplitude information[J]. Journal of Radars, in press. doi:  10.12000/JR19071
Citation: LIU Chao, SUN Jinping, CHEN Xiaolong, et al. Random finite set-based extended target tracking method with amplitude information[J]. Journal of Radars, in press. doi:  10.12000/JR19071
    • 多目标跟踪的目的是从一系列被噪声、漏检和虚警污染的传感器量测中估计出目标的数目和状态[1]。解决此类问题的一个重要环节是建立量测和目标状态之间的关系模型。传统跟踪方法[2]假设一个目标最多只能产生一个量测,所建立的模型称为标准量测模型,目标被称为点目标。这一模型显著降低了多目标跟踪的复杂度,因而获得了广泛应用。但是随着传感器技术的发展,标准量测模型在许多情况下已不再适用。一方面当传感器的分辨率比较高或目标与传感器之间的距离比较近时,目标可能会占据多个分辨单元,从而产生多个量测。例如,当采用车载雷达对车辆进行跟踪和采用激光测距传感器对行人进行跟踪时,都可能引起一个目标产生多个量测的现象,且量测的分布往往能反映出目标的形状特征。这种类型的目标称为扩展目标[3]。另一方面,海上的舰船或空中的飞机等一些军事目标经常以编队的形式出现,当这些目标与传感器之间的距离较远时可能会出现一个“群”目标产生多个量测的情形。由于雷达分辨率的制约,群内的目标往往是不可分辨的,此时合理的处理方式是将这一群目标视为一个独立的目标进行跟踪,因而这类目标被称为“群目标”。由于不满足标准量测模型的要求,扩展/群目标的跟踪需要采用非标准量测模型。

      随着信号处理以及计算机技术的发展,在过去的十多年间,对扩展/群目标跟踪问题的研究已获得越来越多的关注。由于群目标和扩展目标在量测过程上的相似性,许多文献把它们作为一类目标进行研究。Gilholm和Salmond[3]假设每个时刻接收的量测数为泊松分布,在此前提下提出一种扩展目标跟踪方法。Boers等人[4]采用类似的方法对具有1维扩展特性的“竿”形目标进行联合检测和跟踪。Mahler[5]在文献[3]所提模型的基础上推导了一种扩展目标概率假设密度(Extended Target Probability Hypothesis Density, ET-PHD)滤波器,并指出ET-PHD实现的前提是对量测集进行合理划分。Granström等人[6,7]利用高斯混合(Gaussian Mixture, GM)实现了Mahler提出的ET-PHD滤波器,称为ET-GM-PHD,并且提出了距离划分和子划分两种量测划分方法。但是在文献[6,7]中,为降低算法的复杂度,作者只估计了目标中心的运动特性,而忽略了对目标扩展的估计。其后,为实现扩展目标数目的更准确估计,Orguner[8,9]用高斯混合实现了带势估计的PHD滤波器。Koch[10]在假设扩展目标具有椭圆形状的前提下,提出了一种随机矩阵模型[11-14],对目标的运动状态采用高斯分布建模,扩展状态采用逆威沙特分布建模。Granström[11]提出了基于该模型的高斯逆威沙特PHD(Gaussian Inverse Wishart PHD, GIW-PHD)滤波器,弥补了ET-GM-PHD无法估计目标扩展状态的不足,并且针对大小不同且空间邻近的扩展目标分别提出了预测划分和最大期望划分方法。此外,一些学者还提出了其他量测划分方法[15-17]

      目前,有关扩展目标跟踪的文献均直接或间接采用空间信息划分量测,在杂波密集情况下一些杂波量测可能在某些划分中被划入目标单元,从而降低目标扩展状态和运动状态估计的准确度。为解决这一问题,将目标和杂波的幅度信息(Amplitude Information, AI)引入GIW-PHD滤波器,通过为量测单元叠加幅度似然,增大目标量测单元的权重,降低包含杂波的量测单元的权重,从而达到抑制杂波的目的。此外,为了降低杂波对计算中心量测的干扰,采用幅度加权的方法计算量测单元的质量中心,以取代文献[11]中的几何中心,从而可以获得更准确的中心量测估计,并进一步提升对扩展目标的状态估计准确度。

    • ${{\xi}} _k^{\left( i \right)}$表示$k$时刻第$i$个扩展目标的状态,并将$k$时刻扩展目标状态的集合[11]表示为

      $$ \begin{split} {{{X}}_k} \;&= \left\{ {{{\xi}} _k^{\left( i \right)}} \right\}_{i = 1}^{{N_{{{\xi}} ,k}}},\;{{\xi}} _k^{\left( i \right)} \triangleq \left( {\bar {{\xi}} _k^{\left( i \right)},\alpha _k^{\left( i \right)}} \right),\\ & \bar {{\xi}} _k^{\left( i \right)} \triangleq \left( {{{x}}_k^{\left( i \right)},\bar {{X}}_k^{\left( i \right)}} \right) \end{split} $$ (1)

      其中,${N_{{{\xi}} ,k}}$为目标的数量,$\bar{{ \xi}} _k^{\left( i \right)}$包含了第$i$个目标的运动状态${{x}}_k^{\left( i \right)}$以及扩展状态$\bar {{X}}_k^{\left( i \right)}$, $\alpha _k^{\left( i \right)}$为表征目标幅度的参数。令$\left| \cdot \right|$表示集合的势,则有$\left| {{{{X}}_k}} \right| = {N_{{{\xi}} ,k}}$

      假设每个扩展目标的运动状态服从线性高斯模型

      $${{x}}_{k + 1}^{\left( i \right)} = \left( {{{{F}}_{k + 1|k}} \otimes {{{I}}_d}} \right){{x}}_k^{\left( i \right)} + {{w}}_{k + 1}^{\left( i \right)}$$ (2)

      其中,${{w}}_{k + 1}^{\left( i \right)}$是协方差为${{\Delta}} _{k + 1|k}^{\left( i \right)} = {{{Q}}_{k + 1|k}} \otimes \bar {{X}}_{k + 1}^{\left( i \right)}$的零均值高斯过程噪声,$d$为目标扩展维度,$\bar {{X}}_k^{\left( i \right)}$$d \times d$维对称正定矩阵,${{{I}}_d}$$d$维单位矩阵,$ \otimes $为张量积运算符,${{{F}}_{k + 1|k}}$${{{Q}}_{k + 1|k}}$

      $$ {{{F}}_{k + 1|k}} = \left[ {\begin{array}{*{20}{c}} 1&{{T_{\rm{s}}}}&{{{T_{\rm{s}}^2} / 2}} \\ 0&1&{{T_{\rm{s}}}} \\ 0&0&{{{\rm e}^{ - {{{T_{\rm{s}}}} / \theta }}}} \end{array}} \right] \hspace{45pt} $$ (3)
      $${{{Q}}_{k + 1|k}} = {\Sigma ^2}\left( {1 - {{\rm{e}}^{ - 2{{{T_{\rm{s}}}} / \theta }}}} \right){\rm{diag}}\left( {\left[ {0\;0\;1} \right]} \right)$$ (4)

      其中,${T_{\rm{s}}}$为采样间隔,$\Sigma $为标量加速度标准差,$\theta $为机动相关时间。

      $k$时刻目标产生的量测集表示为

      $${{{Z}}_k} = \left\{ {{{z}}_k^{\left( j \right)}} \right\}_{j = 1}^{{N_{{{z}},k}}}$$ (5)

      其中,${{z}}_k^{\left( j \right)}: = \left[ {{{z}}_{{\rm{P}},k}^{\left( j \right)};a_k^{\left( j \right)}} \right]$, ${{z}}_{{\rm{P}},k}^{\left( j \right)}$为位置量测,$a_k^{\left( j \right)}$为幅度量测,${N_{{{z}},k}} = \left| {{{{Z}}_k}} \right|$为量测数。假设传感器具有线性高斯量测模型,则目标的位置量测与目标运动状态之间有如下关系[11]

      $${{z}}_{{\rm{P}},k}^{\left( j \right)} = \left( {{{{H}}_k} \otimes {{{I}}_d}} \right){{x}}_k^{\left( j \right)} + {{e}}_k^{\left( j \right)}$$ (6)

      其中,${{{H}}_k} = \left[ {1\;0\;0} \right]$, ${{e}}_k^{\left( j \right)}$为白色高斯噪声,其协方差由扩展矩阵$\bar {{X}}_k^{\left( j \right)}$给定。每个目标产生的量测数服从泊松分布,泊松率$\gamma \left( {{{{\xi}} _k}} \right)$为目标状态${{{\xi}} _k}$的函数。

      杂波量测数也采用泊松分布建模,单个时间单位体积内的泊松率为${\beta _{{\rm{FA}},k}}$,观测区域的体积为$V$,因此单位时间内杂波量测的均值为${\beta _{{\rm{FA}},k}}V$

    • GIW-PHD滤波器[11]采用高斯逆威沙特分布的混合来近似$k$时刻扩展目标的PHD

      $$ \begin{split} {D_{k|k}}\left( {{{{\xi}} _k}} \right) =& \sum\limits_{j = 1}^{{J_{k|k}}} {w_{k|k}^{\left( j \right)}{\cal N}\left( {{{{x}}_k}|{{m}}_{k|k}^{\left( j \right)},{{P}}_{k|k}^{\left( j \right)} \otimes {{\bar {{X}}}_k}} \right)} \\ & \times {\cal I}{\cal W}\left( {{{\bar X}_k};v_{k|k}^{\left( j \right)},{{V}}_{k|k}^{\left( j \right)}} \right) \end{split} $$ (7)

      其中,${J_{k|k}}$为高斯项数,$w_{k|k}^{\left( j \right)}$, ${{m}}_{k|k}^{\left( j \right)}$${{P}}_{k|k}^{\left( j \right)} \otimes {\bar {{X}}_k}$分别为第$j$个高斯项的权重,均值和协方差。${\cal I}{\cal W}\left( {\bar {{X}};v,{{V}}} \right)$为自由度为$v$且逆尺度矩阵为${{V}}$的逆威沙特分布。

      定义第$j$个GIW项为

      $${{\xi}} _{k|k}^{\left( j \right)} = \left( {{{m}}_{k|k}^{\left( j \right)},{{P}}_{k|k}^{\left( j \right)},v_{k|k}^{\left( j \right)},{{V}}_{k|k}^{\left( j \right)}} \right)$$ (8)

      运动状态的协方差和扩展状态的估计为

      $$\hat {{P}}_{k|k}^{\left( j \right)} = \frac{{{{P}}_{k|k}^{\left( j \right)} \otimes {{V}}_{k|k}^{\left( j \right)}}}{{v_{k|k}^{\left( j \right)} + s - sd - 2}}$$ (9)
      $$ {\hat{ \bar {{X}}}}_{k|k}^{\left( j \right)} = \frac{{{{V}}_{k|k}^{\left( j \right)}}}{{v_{k|k}^{\left( j \right)} - 2d - 2}} \hspace{15pt} $$ (10)

      其中,$s$为目标运动状态的维度。

      GIW-PHD滤波器的具体滤波过程见文献[11],为节省篇幅,这里仅给出更新部分。

      $k$时刻更新后的PHD为GIW的混合式

      $$ {D_{k|k}}\left( {{{\bar{{ \xi}} }_k}} \right) = D_{k|k}^{{\rm{ND}}}\left( {{{\bar {{\xi}} }_k}} \right) + \sum\limits_{{{P}}\angle {{{Z}}_k}} {\sum\limits_{{{W}} \in {{P}}} {D_{k|k}^{{\rm D}}\left( {{{\bar {{\xi}} }_k},{ W}} \right)} } $$ (11)

      其中,${{P}}$为一个划分,${{W}}$为划分后的量测子集,$D_{k|k}^{\rm{ND}}\left( \cdot \right)$为漏检的PHD

      $$ \begin{split} D_{k|k}^{\rm{ND}}\left( {{{\bar {{\xi}} }_k}} \right) =\;& \sum\limits_{j = 1}^{{J_{k|k - 1}}} {w_{k|k}^{\left( j \right)}{\cal N}\left( {{{{x}}_k}|{{m}}_{k|k}^{\left( j \right)},{{P}}_{k|k}^{\left( j \right)}} \right)} \\ & \times {\cal I}{\cal W}\left( {{{\bar {{X}}}_k};v_{k|k}^{\left( j \right)},{{V}}_{k|k}^{\left( j \right)}} \right) \end{split} $$ (12)
      $$ \tag{13a} w_{k|k}^{\left( j \right)} = \left( {1 - \left( {1 - {{\rm e}^{ - \gamma \left( j \right)}}} \right)p_{\rm{D}}^{\left( j \right)}} \right)w_{k|k - 1}^{\left( j \right)} \hspace{10pt} $$
      $$ \tag{13b} \bar {{\xi}} _{k|k}^{\left( j \right)} = \bar {{\xi}} _{k|k - 1}^{\left( j \right)} \hspace{114pt} $$

      $D_{k|k}^{\rm{D}}\left( {{{\bar {{\xi}} }_k},{ W}} \right)$为联合似然函数$\displaystyle\prod\limits_{{{{z}}_k} \in {{W}}} {\dfrac{{{\phi _{{{{z}}_{{\rm P},k}}}}\left( \cdot \right)}}{{{\lambda _k}{c_k}\left( {{{{z}}_{{\rm P},k}}} \right)}}} $与式(7)中的预测GIW项的乘积,可以写为

      $$ \begin{split} & \beta _{{\rm{FA}},k}^{ - \left| {W} \right|}{\cal L}_k^{\left( {j,{{ W}}} \right)}{\cal N}\left( {{{{x}}_k}|{{m}}_{k|k}^{\left( {j,{ W}} \right)},{{P}}_{k|k}^{\left( {j,{ W}} \right)} \otimes {{\bar {{X}}}_k}} \right) \\ &\qquad \times {\cal I}{\cal W}\left( {{{\bar {{X}}}_k};v_{k|k}^{\left( {j,{ W}} \right)},{{V}}_{k|k}^{\left( {j,{ W}} \right)}} \right) \end{split} $$ (14)

      其中

      $$ \tag{15a} {{m}}_{k|k}^{\left( {j,{{W}}} \right)} = {{m}}_{k|k - 1}^{\left( j \right)} + \left( {{{K}}_{k|k - 1}^{\left( {j,{{W}}} \right)} \otimes {{{I}}_d}} \right)\varepsilon _{k|k - 1}^{\left( {j,{{W}}} \right)} \hspace{15pt} $$
      $$ \tag{15b} {{P}}_{k|k}^{\left( {j,{{W}}} \right)} = {{P}}_{k|k - 1}^{\left( j \right)} - {{K}}_{k|k - 1}^{\left( {j,{{W}}} \right)}{{S}}_{k|k - 1}^{\left( {j,{{W}}} \right)}{\left( {{{K}}_{k|k - 1}^{\left( {j,{{W}}} \right)}} \right)^{\rm{T}}} $$
      $$ \tag{15c} v_{k|k}^{\left( {j,{{W}}} \right)} = v_{k|k - 1}^{\left( j \right)} + \left| {{W}} \right| \hspace{98pt} $$
      $$ \tag{15d} {{V}}_{k|k}^{\left( {j,{{W}}} \right)} = {{V}}_{k|k - 1}^{\left( j \right)} + {{N}}_{k|k - 1}^{\left( {j,{{W}}} \right)} + {{Z}}_{{\rm{P}}{\rm{,}}k}^{{W}} \hspace{46pt} $$
      $$ \tag{15e} \bar {{z}}_{{\rm{P}},k}^{ W} = \frac{1}{{\left| {{W}} \right|}}\sum\limits_{{{z}}_k^{\left( i \right)} \in {{W}}} {{{z}}_{{\rm{P}},k}^{\left( i \right)}} \hspace{93pt} $$
      $$ \tag{15f} {{Z}}_{{\rm{P}},k}^{{W}} = \sum\limits_{{{z}}_k^{\left( i \right)} \in {{W}}} {\left( {{{z}}_{{\rm{P}},k}^{\left( i \right)} - \bar {{z}}_{{\rm{P}},k}^{{W}}} \right){{\left( {{{z}}_{{\rm{P}},k}^{\left( i \right)} - \bar {{z}}_{{\rm{P}},k}^{{W}}} \right)}^{\rm{T}}}} \hspace{5pt} $$
      $$ \tag{15g} {{S}}_{k|k - 1}^{\left( {j,{{W}}} \right)} = {{{H}}_k}{{P}}_{k|k - 1}^{\left( j \right)}{{H}}_k^{\rm{T}} + \frac{1}{{\left| {{W}} \right|}} \hspace{56pt} $$
      $$ \tag{15h} {{K}}_{k|k - 1}^{\left( {j,{{W}}} \right)} = {{P}}_{k|k - 1}^{\left( j \right)}{{H}}_k^{\rm{T}}{\left( {{{S}}_{k|k - 1}^{\left( {j,{{W}}} \right)}} \right)^{ - 1}} \hspace{46pt} $$
      $$ \tag{15i} \varepsilon _{k|k - 1}^{\left( {j,{{W}}} \right)} = \bar {{z}}_{{\rm{P}},k}^{{W}} - \left( {{{{H}}_k} \otimes {{{I}}_d}} \right){{m}}_{k + 1|k}^{\left( j \right)} \hspace{41pt} $$
      $$ \tag{15j} {{N}}_{k|k - 1}^{\left( {j,{{W}}} \right)} = {\left( {{{S}}_{k|k - 1}^{\left( {j,{{W}}} \right)}} \right)^{ - 1}}\varepsilon _{k|k - 1}^{\left( {j,{{ W}}} \right)}{\left( {\varepsilon _{k|k - 1}^{\left( {j,{{W}}} \right)}} \right)^{\rm{T}}} \hspace{20pt} $$

      式(14)中的似然函数为

      $$ \begin{split} & {\cal L}_k^{\left( {j,{{W}}} \right)} \\ & = \frac{1}{{{{{\left( {{\pi ^{\left| {{W}} \right|}}\left| {{W}} \right|{{S}}_{k|k - 1}^{\left( {j,{{ W}}} \right)}} \right)}^{\frac{d}{2}}}}}}\frac{{{{\left| {{{V}}_{k|k - 1}^{\left( j \right)}} \right|}^{\frac{{v_{k|k - 1}^{\left( j \right)}}}{2}}}{\varGamma _d}\left( {\dfrac{{v_{k|k}^{\left( {j,{{W}}} \right)}}}{2}} \right)}}{{{{\left| {{{V}}_{k|k}^{\left( {j,{{W}}} \right)}} \right|}^{\frac{{v_{k|k}^{\left( {j,{{ W}}} \right)}}}{2}}}{\varGamma _d}\left( {\dfrac{{v_{k|k - 1}^{\left( j \right)}}}{2}} \right)}}\\ \end{split} $$ (16)

      其中,$\left| {{V}} \right|$表示${{V}}$的行列式,$\left| {{W}} \right|$表示${{W}}$中的量测数。更新后GIW项的权重为

      $$ w_{k|k}^{\left( {j,{{W}}} \right)} = \frac{{{w_{{P}}}}}{{{d_{{W}}}}}{\left( {\frac{{\gamma \left( j \right)}}{{{\beta _{{\rm{FA}},k}}}}} \right)^{\left| { W} \right|}}p_{\rm{D}}^{\left( j \right)}{\cal L}_k^{\left( {j,{{W}}} \right)}w_{k|k - 1}^{\left( j \right)} $$ (17)

      其中,${w_{{P}}}$为某个划分${{P}}$的权重,${d_{{W}}}$为某个量测单元${{W}}$的权重[11]。式(13)~式(17)中下标为$k|k - 1$的量表示从$k{\rm{ - }}1$时刻到$k$时刻的预测值。

    • 图1为扩展目标与量测的关系图。其中,T1和T2为两个大小不同的扩展目标,红色圆点为扩展目标的量测,C1~C3为杂波。由图中可见杂波C1距离目标量测较远,容易划入杂波单元;但是C2和C3距离目标量测较近,因而有可能在某种划分下被划入目标单元,导致目标扩展状态的估计误差,并影响运动状态的估计效果。这种情况下,仅利用空间信息已经很难区分目标和杂波的量测,一种可能的途径是利用目标和杂波的幅度信息。

      图  1  扩展目标量测

      Figure 1.  Measurements of the extended targets

      Lerro[18]建立了幅度信息模型。Clark[19]将这一模型引入了标准PHD滤波器,并且用仿真证明利用目标的幅度信息可以更好的区分目标和杂波。但他所采用的标准PHD滤波器本身仅适用于标准量测模型即点目标,并不适用于扩展目标跟踪场景[6]。本节给出结合幅度信息的GIW-PHD滤波器,称为AI-GIW-PHD滤波器。通过将目标和杂波的幅度信息引入GIW-PHD滤波器,可以增大只包含目标量测的单元的幅度似然,而降低包含杂波量测的单元的幅度似然。由于GIW-PHD在滤波时需要考虑所有可能的量测划分方法和量测子集,因此,这种操作的结果是估计目标扩展状态和运动状态时,仅包含目标量测的单元获得的权重增大,而包含杂波量测的单元的权重降低,从而抑制了杂波对状态估计的影响。

    • 假定目标的幅度和运动状态相互独立,则$k$时刻目标的量测似然函数$g\left( {{{z}}|{{\xi}} } \right)$和杂波的量测似然函数$c\left( {{z}} \right)$[19]可分别表示为

      $$ g\left( {{{{z}}_k}|{{{\xi}} _k}} \right){\rm{ = }}{g_{{{{z}}_{{\rm{P}},k}}}}({{{z}}_{{\rm{P}},k}}\left| {{{\bar {{\xi}} }_k}} \right.){g_{{a_k}}}({a_k}|{\alpha _k}),{a_k} \ge 0 $$ (18)
      $$ c({{{z}}_k}){\rm{ = }}{c_{{{{z}}_{{\rm{P}},k}}}}({{{z}}_{{\rm{P}},k}}){c_{{a_k}}}({a_k}),{a_k} \ge 0 \hspace{50pt} $$ (19)

      其中,${g_{{{{z}}_{{\rm{P}},k}}}}({{{z}}_{{\rm{P}},k}}\left| {{{\bar {{\xi}} }_k}} \right.)$${c_{{{{z}}_{{\rm{P}},k}}}}({{{z}}_{{\rm{P}},k}})$分别为目标和杂波关于运动状态的似然函数,${g_{{a_k}}}({a_k}|{\alpha _k})$${c_{{a_k}}}({a_k})$分别为目标和杂波关于幅度的似然函数。由此计算检测概率$p_{\rm{D}}^\tau $和虚警概率$p_{{\rm{FA}}}^\tau $[18]

      $$p_{\rm{D}}^\tau = \int\limits_\tau ^\infty {{g_a}(a|\alpha ){\rm{d}}a} $$ (20)
      $$ p_{{\rm{FA}}}^\tau = \int\limits_\tau ^\infty {{c_a}(a){\rm{d}}a} \hspace{5pt} $$ (21)

      其中,$\tau $为检测门限。则过门限检测后,目标和杂波的幅度似然函数分别为

      $$g_a^\tau (a|\alpha ) = \frac{1}{{p_D^\tau }}{g_a}(a|\alpha ),\;\;a > \tau $$ (22)
      $$ c_a^\tau (a) = \frac{1}{{p_{FA}^\tau }}{c_a}(a),\;\;a > \tau \hspace{15pt} $$ (23)
    • 当杂波为复高斯分布且幅度为Swerling 1型起伏时,杂波和目标的幅度为Rayleigh分布。杂波幅度的概率密度函数(Probability Density Function, PDF)为[19]

      $${c_a}(a) = \frac{{2a}}{{{\sigma _n}}}\exp \left( { - \frac{{{a^2}}}{{{\sigma _n}}}} \right), a > 0$$ (24)

      其中,${\sigma _n}$为复高斯分布的协方差。目标幅度${A_k}$的PDF为

      $$ {g_a}({A_k}) = \frac{{2{A_k}}}{{{\sigma _t}}}\exp \left( { - \frac{{{A_k}^{\!\! 2}}}{{{\sigma _t}}}} \right),\;{A_k} > 0 $$ (25)

      其中,${\sigma _t}$为目标幅度的均方值。则目标加杂波的幅度的PDF为

      $${g_a}(a) = \frac{{2a}}{{{\sigma _n} + {\sigma _t}}}\exp \left( { - \frac{{{a^2}}}{{{\sigma _n} + {\sigma _t}}}} \right), a > 0$$ (26)

      现在令$\alpha = {{{\sigma _t}} / {{\sigma _n}}}$,并令${\sigma _n} = 1$,则有$\alpha = {\sigma _t}$。将${\sigma _n}$$\alpha $代入式(24)和式(26)可得

      $$ {c_a}(a) = 2a\exp \left( { - {a^2}} \right),\;\;a > 0 \hspace{47pt} $$ (27)
      $${g_a}(a|\alpha ) = \frac{{2a}}{{1 + \alpha }}\exp \left( { - \frac{{{a^2}}}{{1 + \alpha }}} \right),\;\;a > 0$$ (28)

      由于信杂比的定义[19]

      $${\rm{SCR}} = 10\lg \left( {\frac{{{\sigma _n} + {\sigma _t}}}{{{\sigma _n}}}} \right)$$ (29)

      ${\sigma _n}$$\alpha $代入式(29)可得

      $${\rm{SCR}} = 10\lg \left( {1 + \alpha } \right)$$ (30)
    • 下面给出AI-GIW-PHD的具体算法。

      由于AI-GIW-PHD只改进了GIW-PHD的更新部分,而预测部分保持不变,因此这里只给出AI-GIW-PHD的更新步骤。

      $$ {D_{k|k}}\left( {{{\bar {{\xi}} }_k}} \right) = D_{k|k}^{\rm{ND}}\left( {{{\bar {{\xi}} }_k}} \right) + \sum\limits_{{{P}}\angle {{{W}}_k}} {\sum\limits_{{{W}} \in {{P}}} {D_{k|k}^{\rm{D}}\left( {{{\bar {{\xi}} }_k},{{W}}} \right)} } $$ (31)
      $$ \begin{split} D_{k|k}^{\rm{ND}}\left( {{{\bar {{\xi}} }_k}} \right) =\;& \sum\limits_{j = 1}^{{J_{k|k - 1}}} {w_{k|k}^{\left( j \right)}{\cal N}\left( {{{{x}}_k}|{{m}}_{k|k}^{\left( j \right)},{{P}}_{k|k}^{\left( j \right)}} \right)} \\ & \times {\cal I}{\cal W}\left( {{{\bar {{X}}}_k};v_{k|k}^{\left( j \right)},{{V}}_{k|k}^{\left( j \right)}} \right) \hspace{35pt} \end{split} $$ (32)
      $$ w_{k|k}^{\left( j \right)} = \left( {1 - \left( {1 - {{\rm{e}}^{ - \gamma \left( j \right)}}} \right)p_{\rm{D}}^{\tau ,\left( j \right)}} \right)w_{k|k - 1}^{\left( j \right)} \hspace{18pt} $$ (33)

      其中,$p_{\rm{D}}^{\tau ,\left( j \right)}$为根据式(20)计算出的检测概率。

      每个单元中所有量测的似然函数为

      $$ \prod\limits_{{{{z}}_k} \in { W}} {\frac{{{\phi _{{{{z}}_{{\rm{P}},k}}}}\left( \cdot \right)p_{\rm{D}}^\tau \left( \alpha \right)g_a^\tau (a|\alpha )}}{{{\lambda _k}{c_k}\left( {{{{z}}_{{\rm{P}},k}}} \right)p_{{\rm{FA}}}^\tau c_a^\tau (a)}}} $$ (34)

      其中,$p_{\rm{D}}^\tau \left( \alpha \right)$$p_{{\rm{FA}}}^\tau $分别由式(20)和式(21)给出,$g_a^\tau (a|\alpha )$$c_a^\tau (a)$分别由式(22)和式(23)给出。

      根据式(34),增加幅度信息后式(16)中的${\cal L}_k^{\left( {j,{{ W}}} \right)}$变为

      $$\tilde {\cal L}_k^{\left( {j,{{W}}} \right)} = {\cal L}_k^{\left( {j,{{W}}} \right)}\prod\limits_{{{{z}}_k} \in {{W}}} {\frac{{p_{\rm{D}}^\tau \left( \alpha \right)g_a^\tau (a|\alpha )}}{{p_{{\rm{FA}}}^\tau c_a^\tau (a)}}} $$ (35)

      将式(35)代入式(17)可得GIW项的新的权重为

      $$ w_{k|k}^{\left( {j,{{W}}} \right)} = \frac{{{w_{{P}}}}}{{{d_{{W}}}}}{\left( {\frac{{\gamma \left( j \right)}}{{{\beta _{{\rm{FA}},k}}}}} \right)^{\left| {{W}} \right|}}p_{\rm{D}}^{\tau ,\left( j \right)}\tilde {\cal L}_k^{\left( {j,{{W}}} \right)}w_{k|k - 1}^{\left( j \right)} $$ (36)
    • 根据式(15(e)),当计算某个单元的中心量测时,对该单元中的所有量测取几何平均值。然而,当该单元中混入了杂波量测时,采用这一方法将导致计算得到的中心量测偏离真实的中心量测,从而增大状态估计误差。为解决这一问题,计算中心量测时,采用幅度加权的方法来计算量测的质量中心,以取代式(15(e))中所用的几何中心

      $$ \bar {{z}}_{{\rm{P}},k}^{ W} = \frac{1}{{\left| {{W}} \right|}}\sum\limits_{{{z}}_k^{\left( i \right)} \in {{W}}} {{{\rm{\omega}} ^{\left( i \right)}}{{z}}_{{\rm{P}},k}^{\left( i \right)}} $$ (37)

      其中,幅度的加权因子${{\rm{\omega}} ^{\left( i \right)}}$

      $$ {{\rm{\omega}} ^{\left( i \right)}}{\rm{ = }}\frac{{{a^{\left( i \right)}}}}{{\displaystyle\sum\limits_{{{z}}_k^{\left( i \right)} \in {{W}}} {{a^{\left( i \right)}}} }} $$ (38)

      其中,${a^{\left( i \right)}}$为量测$i$的幅度。

      通过幅度加权,当量测单元中包含了杂波时,由于通常情况下杂波的幅度小于目标量测的幅度,因此能够更多地利用目标量测的信息,并减少杂波对中心量测的影响。

    • 为检验所提AI-GIW-PHD算法的性能,通过仿真实验将其与标准的GIW-PHD滤波器进行对比。

    • 雷达观测区域的面积为[–500, 500] × [–500, 500] m。在总计31 s的观测时间内有4个目标进入观测区域。目标的真实航迹如图2所示。图中每种颜色的航迹代表一个目标,可以看到有2个扩展目标为匀速直线运动,另外2个目标为匀速转弯运动,且目标的航迹之间出现了多次交叉。

      图  2  扩展目标真实航迹(SCR=13 dB)

      Figure 2.  Real trajectories of the extended targets (SCR=13 dB)

      目标的真实扩展为

      $$ \bar {{X}}_k^{\left( i \right)} = {{R}}_k^{\left( i \right)}{\rm{diag}}\left( {\begin{array}{*{20}{c}} {\upsilon _i^2}&{\iota _i^2} \end{array}} \right){\left( {{{R}}_k^{\left( i \right)}} \right)^{\rm{T}}} $$

      其中,${{R}}_k^{\left( i \right)}$为旋转矩阵,用于保证目标的主轴平行于目标的运动方向,${\upsilon _i}$${\iota _i}$分别为长轴和短轴的长度。本文仿真主要考察幅度对跟踪性能的影响,因此4个目标的长短轴均为$\left( {{\upsilon _i},{\iota _i}} \right) = \left( {20,5} \right)$。扩展目标产生量测的期望数与扩展目标的尺寸${\cal V}_k^{\left( i \right)} \triangleq \pi \sqrt {\left| {X_k^{\left( i \right)}} \right|} = \pi {\upsilon _i}{\iota _i}$有关,这里采用如下的模型

      $$ \gamma _k^{\left( i \right)} = \left\lfloor {\sqrt {\frac{4}{\pi }{\cal V}_k^{\left( i \right)}} + 0.5} \right\rfloor = \left\lfloor {2\sqrt {{\upsilon _i}{\iota _i}} + 0.5} \right\rfloor $$ (39)

      每个扩展目标产生的量测,在幅度上均独立的服从参数为${\sigma _t}$的Rayleigh分布,且在帧间独立的服从Swerling 1型起伏。

      杂波在每一帧的总数服从均值为1000的泊松分布,虚警率为0.05,则理论上每一帧的虚警数服从均值为50的泊松分布。

      运动模型的参数为:${T_{\rm{s}}} = 1 \; {\rm{s}}$, $\theta = 1 \; {\rm{s}}$, $\Sigma $=0.1 ${\rm{m}}/{{\rm{s}}^2}$, $\tau {\rm{ = 5}} \; {\rm{s}}$。新生目标的参数为${J_{b,k}}$=4, $w_{b,k}^{\left( j \right)}=0.1$,并且有$m_{b,k}^{\left( j \right)} = {\left[ {{{\left( {{{x}}_0^{\left( j \right)}} \right)}^{\rm{T}}}{\bf{0}}_4^{\rm{T}}} \right]^{\rm{T}}}$, $P_{b,k}^{\left( j \right)}$= diag([1002 252 252]), $v_{b,k}^{\left( j \right)} = 7$, ${{V}}_{b,k}^{\left( j \right)} = {\rm{diag}}\left( {\left[ {1{\rm{ 1}}} \right]} \right)$,即新生目标的起始位置为每个目标航迹的起始处。

      设置目标的信杂比分别为SCR=13 dB和SCR= 6 dB。本文通过计算每种算法的OSPA误差[20]来评价算法的跟踪性能。

    • 图3图4为SCR=13 dB时未利用幅度信息的GIW-PHD和利用幅度信息的AI-GIW-PHD在某次仿真中的航迹滤波结果。图3中,黑色“*”表示目标真实位置,红色“Ο”表示状态估计值。从该图可以看到,GIW-PHD出现了比较明显的航迹漏检现象。这是因为当杂波率较高时,目标量测与杂波量测有可能相距很近,状态估计容易受到杂波的干扰,导致一些GIW项权重降低,从而估计出的目标数目也会减少。图4中,蓝色“.”表示目标真实位置,黑色“Ο”表示状态估计值。可以看到,利用了幅度信息的AI-GIW-PHD可以非常准确地估计出目标航迹,基本上没有漏检。这是因为,由于采用了幅度信息,相距较近的目标和杂波的量测能够得到更好的区分,从而降低了杂波对目标GIW项的干扰。

      图  3  GIW-PHD单次航迹估计结果(SCR=13 dB)

      Figure 3.  Track estimation of the GIW-PHD in a run (SCR=13 dB)

      图  4  AI-GIW-PHD单次航迹估计结果(SCR=13 dB)

      Figure 4.  Track estimation of the AI-GIW-PHD in a run (SCR=13 dB)

      图5为GIW-PHD在该次仿真中的扩展状态估计结果。图5中绿色填充带表示目标的真实扩展状态,红色椭圆表示GIW-PHD估计出的目标扩展状态。可以看到,由于存在大量的航迹漏检,GIW-PHD对目标扩展状态的估计也出现了大量漏检。当能够估计出目标航迹时,GIW-PHD也能够比较准确地估计出目标的扩展状态,但在滤波开始时刻还存在比较明显的估计错误。

      图  5  GIW-PHD单次扩展状态估计结果(SCR=13 dB)

      Figure 5.  Extended state estimation of the GIW-PHD in a run (SCR=13 dB)

      图6为AI-GIW-PHD在该次仿真中的扩展状态估计结果。图6中蓝色填充带表示目标的真实扩展状态,黑色椭圆表示AI-GIW-PHD估计出的目标扩展状态。可以看到,AI-GIW-PHD始终能够比较准确地估计出目标的扩展状态,并且在滤波初始时刻也没有出现明显的错误估计。

      图  6  AI-GIW-PHD单次扩展状态估计结果(SCR=13 dB)

      Figure 6.  Extended state estimation of the AI-GIW-PHD in a run (SCR=13 dB)

      图7为SCR=13 dB时GIW-PHD和AI-GIW-PHD经过20次蒙特卡洛仿真后统计得到的平均OSPA位置误差。可以看到,在大部分时刻AI-GIW-PHD的误差比GIW-PHD稍大一些。但是,这并不能说明AI-GIW-PHD的滤波性能不如GIW-PHD。这是因为与单目标跟踪时位置误差的计算不同,多目标跟踪时OSPA位置误差的计算需要统计所有航迹的滤波误差。由于GIW-PHD出现了明显的航迹漏检,用于误差统计的航迹数少于AI-GIW-PHD,因而统计得到的OSPA位置误差小于AI-GIW-PHD。当GIW-PHD不发生漏检时,其位置误差会更大。这与文献[19]中的情况相似。虽然从图7无法直接判定GIW-PHD和AI-GIW-PHD的滤波性能,但是根据图3图4可以看到,当没有漏检时,这两种算法都能够准确地跟踪目标。

      图  7  平均OSPA位置误差(SCR=13 dB)

      Figure 7.  Averaged OSPA location error (SCR=13 dB)

      图8为SCR=13 dB时GIW-PHD和AI-GIW-PHD经过20次仿真后得到的平均目标数目估计值。可以看到,随着滤波的进行,GIW-PHD总体上会出现越来越明显的数目欠估计现象。这是因为杂波对GIW项的干扰越来越严重。相比之下,结合了幅度信息的AI-GIW-PHD始终保持了对目标数目的准确估计。这是因为,通过利用幅度信息,杂波与目标量测的区分会更加容易,从而使GIW项受到的杂波干扰显著下降。

      图  8  平均势估计结果(SCR=13 dB)

      Figure 8.  Averaged cardinality estimation (SCR=13 dB)

      图9为SCR=13 dB时GIW-PHD和AI-GIW-PHD经过20次仿真后得到的平均OSPA势误差。通过该图,可以更加清楚地看到在目标数目估计方面,AI-GIW-PHD显著优于GIW-PHD。这再次证明了幅度信息对扩展目标跟踪的有效性。

      图  9  平均OSPA势误差(SCR=13 dB)

      Figure 9.  Averaged OSPA cardinality error (SCR=13 dB)

      图10为SCR=6 dB时GIW-PHD和AI-GIW-PHD经过20次仿真后统计得到的平均OSPA位置误差。结合图10图7可以看到,随着SCR从13 dB降到6 dB,未采用幅度信息的GIW-PHD的跟踪误差明显增大,而采用了幅度信息的AI-GIW-PHD的跟踪误差变化不大。这是因为当SCR降低后目标的检测概率也会下降,因而GIW-PHD算法中GIW项受到的杂波干扰会大大增加,导致状态估计出现较大偏差。相比之下,AI-GIW-PHD由于利用了幅度信息从而能够较好地区分目标和杂波,减轻了GIW项受到的杂波干扰,保持了状态估计的准确度。

      图  10  平均OSPA位置误差(SCR=6 dB)

      Figure 10.  Averaged OSPA location error (SCR=6 dB)

      图11为SCR=6 dB时GIW-PHD和AI-GIW-PHD经过20次仿真后得到的平均目标数目估计值。可以看到,在低SCR条件下,GIW-PHD由于受到杂波的干扰导致目标数目估计始终存在很大的误差,而AI-GIW-PHD在大部分时刻保持了对目标数目的准确估计,仅在倒数第2帧出现了比较明显的偏差。

      图  11  平均势估计结果(SCR=6 dB)

      Figure 11.  Averaged cardinality estimation (SCR=6 dB)

      图12为SCR=6 dB时GIW-PHD和AI-GIW-PHD经过20次仿真后得到的平均OSPA势误差。可以看到,与图11显示的结果一致,AI-GIW–PHD在目标数目估计方面显著优于GIW-PHD。

      图  12  平均OSPA势误差(SCR=6 dB)

      Figure 12.  Averaged OSPA cardinality error (SCR=6 dB)

    • GIW-PHD扩展目标滤波器在进行量测划分和状态估计时都只利用了量测的空间信息,在杂波密集环境下状态估计和势估计性能不佳。为解决这一问题,提出了结合幅度信息的GIW-PHD滤波器,利用量测的幅度信息来区分相距较近的目标和杂波量测,赋予目标量测单元更高的似然值,从而提升目标GIW项的权重,降低杂波的影响。在计算量测单元的中心时采用幅度加权的质量中心取代仅包含空间信息的几何中心,从而进一步降低杂波的干扰。仿真结果表明,所提算法对扩展目标的势估计性能显著增强,在低信杂比环境下的状态估计性能明显优于GIW-PHD滤波器。但是,计算幅度似然会在一定程度上增加算法的运算成本。

参考文献 (20)

目录

    /

    返回文章
    返回