Skip to main content

Gaussian mixture probability hypothesis density filter for multipath multitarget tracking in over-the-horizon radar

Abstract

Conventional multitarget tracking systems presume that each target can produce at most one measurement per scan. Due to the multiple ionospheric propagation paths in over-the-horizon radar (OTHR), this assumption is not valid. To solve this problem, this paper proposes a novel tracking algorithm based on the theory of finite set statistics (FISST) called the multipath probability hypothesis density (MP-PHD) filter in cluttered environments. First, the FISST is used to derive the update equation, and then Gaussian mixture (GM) is introduced to derive the closed-form solution of the MP-PHD filter. Moreover, the extended Kalman filter (EKF) is presented to deal with the nonlinear problem of the measurement model in OTHR. Eventually, the simulation results are provided to demonstrate the effectiveness of the proposed filter.

1 Introduction

Over-the-horizon radar (OTHR) exploits skywave propagation of high-frequency signals to detect and track targets, which are different from the conventional radar. It has received wide attention because of its wide area surveillance, long detection range, strong anti-stealth ability, the capability of the long early warning time, and so on. In OTHR, a significant problem is the effect of multipath propagation, which causes multiple detections via different propagation paths for a target with missed detections and false alarms at the receiver [16]. Nevertheless, the conventional tracking algorithms, such as probabilistic data association (PDA) [79], presume that a single-measurement per target, it may consider the other measurements of the same target as clutter, and multiple tracks are produced when a single target is present. Therefore, these methods cannot effectively solve the multipath propagation problem.

In order to solve the multipath propagation problem of OTHR, lots of algorithms have been proposed. Such as the multipath viterbi data association (MVDA) [10] and modified probabilistic data association (MPDA) [11]. However, these conventional algorithms involve data association problem and the complexity of the computation. In addition, it is noted that most of the existing algorithms in OTHR focus on a single target. Therefore, the multitarget tracking algorithm in OTHR is still an open problem.

Most of the conventional algorithms about multitarget tracking, such as multiple hypothesis tracker (MHT) [12], joint probabilistic data association (JPDA) [13, 14], and probability hypothesis density (PHD) filter [15], assume the following measurement model: (1) every target produces at most one measurement and (2) any measurement is produced by a target or clutter. In this paper, we consider the measurement model which satisfies these assumptions as standard measurement model. However, many measurement models in real-life target tracking scenarios do not satisfy these prerequisite of assumption, which are treated as nonstandard measurement model. Recently, the multiple detection joint probabilistic data association (MD-JPDA) filter [16] based on the JPDA framework was proposed to deal with the multiple detection targets, which can apply to OTHR. However, this algorithm involves complexity data association. To solve this problem of multitarget tracking, more effective methods should be taken into account to implement in OTHR.

In 2009, Mahler derived a series of “second-generation” PHD/cardinalized probability hypothesis density (CPHD) filters [17] which addressed the tracking problems for nonstandard multitarget measurement model, such as extended targets [1821], unresolved target [22], unknown clutter [23], and superpositional sensors [24]. As for the multipath propagation in OTHR, the receiver can obtain more than one measurement from the same target, which means that it violates the second assumption of the standard measurement model. Therefore, the multipath tracking model in OTHR is one of the nonstandard measurement models, such that conventional multitarget tracking algorithms would no longer be applicable.

To solve the multipath propagation problem of multitarget, we propose a novel tracking algorithm based on the theory of finite set statistics (FISST) called the multipath probability hypothesis density (MP-PHD) filter in OTHR. First, inspired by Mahler’s work in [17], the FISST is used to derive the update equation of the MP-PHD filter. Then, a method of Gaussian mixture (GM) is introduced to derive the closed-form solution of the MP-PHD filter. In addition, the extended Kalman filter (EKF) is presented to deal with the nonlinear problem of the measurement model in OTHR. The simulation results demonstrate that the MP-PHD filter can accurately estimate the target state and the target number for low detection probability under the multitarget tracking circumstance.

This paper is organized as follows. Section 2 presents the problem formulation for OTHR. Section 3 then formulates the random finite set (RFS) measurement model for the OTHR tracking problem and derives the MP-PHD filter in OTHR based on the FISST. A closed-form implementation of MP-PHD filter is presented in Section 4. In Section 5, the simulation results are presented. Finally, some conclusions and future works are discussed in Section 6.

2 Problem formulation

2.1 Dynamic model

The geometry of the target and OTHR system is depicted in Fig. 1. The target state vector of OTHR at time k is defined by \( {x}_k={\left[\rho (k),\overset{.}{\rho }(k),b(k),\overset{.}{b}(k)\right]}^{\prime },{x}_k\in \Im \), where ρ(k) = ρ 1(k), \( \overset{.}{\rho }(k) \), b(k), and \( \overset{.}{b}(k) \) are the ground range, range rate, bearing, and bearing rate and is the space of the ground coordinates, the dash denotes transposition.

Fig. 1
figure 1

Geometry of a planar OTHR measurement model with a collocated transmitter and receiver. h r and h t denote reflecting ionosphere heights

Since the distance between the receiver and the targets is large, we usually assume that the state equation of OTHR is linear and discrete-time. It could be modeled as

$$ {x}_k=F{x}_{k-1}+{u}_{k-1} $$
(1)

where u k − 1 is a zero mean, white Gaussian noise with covariance Q k − 1, and the state transition matrix F is given by

$$ F=\left[\begin{array}{cccc}\hfill 1\hfill & \hfill T\hfill & \hfill 0\hfill & \hfill 0\hfill \\ {}\hfill 0\hfill & \hfill 1\hfill & \hfill 0\hfill & \hfill 0\hfill \\ {}\hfill 0\hfill & \hfill 0\hfill & \hfill 1\hfill & \hfill T\hfill \\ {}\hfill 0\hfill & \hfill 0\hfill & \hfill 0\hfill & \hfill 1\hfill \end{array}\right] $$
(2)

where T is the sampling period.

2.2 Measurement model of OTHR

OTHR achieves beyond the line-of-sight horizon target tracking by exploiting skywave propagation of high-frequency signals via the ionospheric, and depending on the particular ionospheric propagation conditions, there may be several propagation paths. As shown in Fig. 1, the OTHR signal reflects from the transmitter to the target by an ionospheric at height h t, then the signal reflects from the target to the receiver by an ionospheric at height h r. As is well known, there are many ionospheric layers present in the atmosphere, and we suppose that there are two ionospheric layers E and F for simplicity, with constant heights h E and h F (this suppose is the same as other algorithms in OTHR). Therefore, there are four propagation modes: mode EE (transmit on E and receive on E); mode EF (transmit on E and receive on F); mode FE (transmit on F and receive on E); and mode FF (transmit on F and receive on F). It means that the receiver can obtain more than one measurement from the same target at one time.

The OTHR measurements consist of a slant range Rg = r 1 + r 2, Doppler f d , and Az = π/2 ‐ θ (θ is the azimuth as shown in Fig. 1) of the form z k  = [Rg(k), f d (k), Az(k)]′, \( {z}_k\in \mathcal{G} \), where \( \mathcal{G} \) is the space of the slant coordinates.

The measurement model of OTHR involves the mapping from the ground coordinates \( \left(\rho, \overset{.}{\rho },b,\overset{.}{b}\right) \) to the slant coordinates (Rg, f d , Az). From Fig. 1, it can be shown that the mapping is expressed as [2]

$$ \left\{\begin{array}{c}\hfill \mathrm{R}\mathrm{g}={r}_1+{r}_2\hfill \\ {}\hfill {f}_d=\frac{\overset{.}{\rho }}{4}\left\{\frac{\rho }{r_1}+\frac{\eta }{r_2}\right\}\hfill \\ {}\hfill \mathrm{A}\mathrm{z}={ \sin}^{-1}\left\{\rho \sin (b)/\left(2{r}_1\right)\right\}\hfill \end{array}\right. $$
(3)

where

$$ {r}_1=\sqrt{{\left(\rho /2\right)}^2+{h}_r^2} $$
(4)
$$ {r}_2=\sqrt{{\left(\rho /2\right)}^2-d\rho \sin (b)/2+{\left(d/2\right)}^2+{h}_t^2} $$
(5)
$$ \eta =\rho -d \sin (b) $$
(6)

Then, the measurement model of OTHR can be expressed as [11]

$$ {z}_k=\Big\{\kern1em \begin{array}{c}{h}_1\left({x}_k\right)+{w}_{k,1}\kern1em if\kern.3em mode\kern.5em EE\kern1em \\ {}\kern1em {h}_2\left({x}_k\right)+{w}_{k,2}\kern1em if\kern.3em mode\kern.5em EF\kern1em \\ {}\kern1em {h}_3\left({x}_k\right)+{w}_{k,3}\kern1em if\kern.3em mode\kern.5em FE\kern1em \\ {}\kern0.6em {h}_4\left({x}_k\right)+{w}_{k,4}\kern1em if\kern.3em mode\kern.5em FF\kern1em \\ {}\kern1em \mathrm{clutter}\kern3.5em otherwise\kern1em \end{array} $$
(7)

where x k is the state variable; w k,i is the zero mean, white Gaussian noise with known covariance R k,i , and h i () is the nonlinear measurement function of ith propagation mode, i = 1, 2, 3, 4.

$$ {h}_i(x)=\left[\begin{array}{c}\hfill \sqrt{{\left(\rho /2\right)}^2+{h}_r^2}+\sqrt{{\left(\rho /2\right)}^2-d\rho \sin (b)/2+{h}_t^2}\hfill \\ {}\hfill \frac{\overset{.}{\rho }}{4}\left(\frac{\rho }{\sqrt{{\left(\rho /2\right)}^2+{h}_r^2}}+\frac{\rho -d \sin (b)}{\sqrt{{\left(\rho /2\right)}^2-d\rho \sin (b)/2+{h}_t^2}}\right)\hfill \\ {}\hfill { \sin}^{-1}\left\{\rho \sin (b)/\left(2\sqrt{{\left(\rho /2\right)}^2+{h}_t^2}\right)\right\}\hfill \end{array}\right] $$
(8)

where, for a given propagation mode i, h t and h r must be replaced by the virtual ionosphere heights, where the details can be found in [11].

3 The RFS model and PHD filter implementation in OTHR

In the standard PHD filter, it is presumed that a target produces at most one measurement. Due to the multipath effect of OTHR, a single target may generate several measurements. Therefore, the measurement model of the MP-PHD filter proposed in this paper is different from the standard PHD filter. In this section, we first introduce the RFS measurement model for OTHR and then derive the update equation of the MP-PHD filter based on the theory of FISST.

3.1 The RFS measurement model for OTHR

At each time k, a finite set of measurements of OTHR, denoted by \( {Z}^k=\left\{{z}_{k,1},{z}_{k,2},\cdots {z}_{k,{N}_k}\right\} \), where \( {z}_{k,1},{z}_{k,2},\cdots {z}_{k,{N}_k} \) are the received measurements at time k and N k is the number of measurements. Since the measurements include several detections from targets and clutter, the collection of measurements of OTHR can be modeled by RFS

$$ {Z}^k={\varTheta}_{k,1}\left({x}_k\right)\cup {\varTheta}_{k,2}\left({x}_k\right)\cup {\varTheta}_{k,3}\left({x}_k\right)\cup {\varTheta}_{k,4}\left({x}_k\right)\cup {\varGamma}_k $$
(9)

where Θ k,i (x k ), i = 1, , 4 denotes the measurement originated from the ith propagation path and Γ k denotes the RFS of clutter. It is presumed that conditional on x k , Θ k,i (x k ), i = 1, , 4 and Γ k are independent RFSs.

3.2 The update equation of MP-PHD filter in OTHR

Note that both the MP-PHD filter and the standard PHD filter recursion require two steps: prediction and update. In the following subsection, we only derive the update equation according to the above RFS measurement model since the prediction step of the MP-PHD filter is identical to the standard PHD filter.

In ref. [17], the probability generating functional (PGFL) of update equation for the multitarget Bayes filter can be written as

$$ {G}_{k\Big|k}\left[h\right]=\frac{\frac{\delta F}{\delta {Z}_k}\left[0,h\right]}{\frac{\delta F}{\delta {Z}_k}\left[0,1\right]} $$
(10)

where F[g, h] is two-variable PGFL as

$$ F\left[g,h\right]={\displaystyle \int {h}^X}\cdot {G}_k\left[g\Big|X\right]\cdot {f}_{k\Big|k-1}\left(X\Big|{Z}^{\left(k-1\right)}\right)\delta X $$
(11)
$$ {G}_k\left[g\Big|X\right]={\displaystyle \int {g}^Z\cdot {f}_k\left(Z\Big|X\right)\delta Z} $$
(12)

where f k|k − 1(X|Z (k − 1)) is the predicted multitarget distribution, f k (Z|X) is the multitarget likelihood function and h X is defined by h X = 1 when X =  and h X = ∏ xX h(x) otherwise. The update PHD filter can be given by

$$ {D}_{k\Big|k}(x)=\frac{\delta {G}_{k\Big|k}}{\delta x}\left[1\right] $$
(13)

Therefore, the update equation for MP-PHD filter can be derived by the following procedures:

  1. 1)

    Derive the PGFL G k [g|X] by using the OTHR measurements and exploiting Eq. (11) to derive a closed form of F[g, h]

  2. 2)

    Derive the derivatives of F[g, h] by using Eq. (10) to derive a closed form of G k|k [h]

  3. 3)

    According to Eq. (11), derive a closed-form update for the MP-PHD filter in OTHR

In the following part of this subsection, we will present the details of the derivation. First, we derive the PGFL G k [g|X] according to the OTHR measurements in Eq. (9). Since Θ k (x) = Θ k,1(x) Θ k,2(x) Θ k,3(x) Θ k,4(x), the PGFL G k [g|x] of Θ k (x) is

$$ \begin{array}{l}{G}_k\left[g\Big|x\right]={\displaystyle \prod_{i=1}^4{G}_{k,i}\left[g\Big|x\right]}={\displaystyle \prod_{i=1}^4\left(1-{p}_{D,k}(x)+{p}_{D,k}(x){p}_{g,i}(x)\right)}\\ {}\begin{array}{ccc}\hfill \hfill & \hfill \hfill & \hfill \hfill \end{array}={\left(1-{p}_{D,k}(x)\right)}^4+{\left(1-{p}_{D,k}(x)\right)}^3{p}_{D,k}(x)\left({\displaystyle \sum_{i=1}^4{p}_{g,i}(x)}\right)+{\left(1-{p}_{D,k}(x)\right)}^2{p}_{D,k}^2(x)\\ {}\begin{array}{ccc}\hfill \hfill & \hfill \hfill & \hfill \hfill \end{array}\left({\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{p}_{g,i}(x){p}_{g,j}(x)}}\right)+\left(1-{p}_{D,k}(x)\right){p}_{D,k}^3(x)\left({\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{p}_{g,i}(x){p}_{g,j}(x)}}}{p}_{g,m}(x)\right)\\ {}\begin{array}{ccc}\hfill \hfill & \hfill \hfill & \hfill \hfill \end{array}+{p}_{D,k}^4(x)\left({\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{p}_{g,i}(x){p}_{g,j}(x){p}_{g,m}(x){p}_{g,n}(x)}}}}\right)\end{array} $$
(14)

where p g,i (x) = ∫g(z)g k,i (z|x)dz, p D,k (x) is the detection probability, and g k,i (z|x) is the likelihood of the ith propagation path. Consequently, if we abbreviate q D,k (x) = 1 − p D,k (x), Eq. (11) can be summarized as follows:

$$ \begin{array}{c}F\left[g,h\right]={e}^{\lambda c\left[g\right]-\lambda}\cdot {G}_{k\Big|k-1}\left[h\right({q}_{D,k}^4+{q}_{D,k}^3{p}_{D,k}\left({\displaystyle \sum_{i=1}^4{p}_{g,i}}\right)+{q}_{D,k}^2{p}_{D,k}^2\left({\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{p}_{g,i}{p}_{g,j}}}\right)\\ {}+{q}_{D,k}{p}_{D,k}^3\left({\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{p}_{g,i}{p}_{g,j}{p}_{g,m}}}}\right)+{p}_{D,k}^4\left({\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{p}_{g,i}{p}_{g,j}{p}_{g,m}{p}_{g,n}}}}}\right)\left)\right]\end{array} $$
(15)

where c[g] = ∫g(z)c k (z)dz and G k|k − 1[h] = ∫h X f k|k − 1(X|Z (k − 1))δX.

To derive a closed-form update equation for the multipath multitarget PHD filter, we assume that the predicted multitarget distribution is a Poisson process:

$$ {G}_{k\Big|k-1}\left[h\right]={e}^{\mu s\left[h\right]-\mu } $$
(16)

where s[h] = μ − 1h(x)D k|k − 1(x)dx and μ = ∫D k|k − 1(x)dx. Thus, Eq. (15) can be written as

$$ \begin{array}{c}F\left[g,h\right]= \exp \Big(\lambda c\left[g\right]-\lambda +\mu s\left[h{q}_{D,k}^4\right]+\mu s\left[h{q}_{D,k}^3{p}_{D,k}\left({\displaystyle \sum_{i=1}^4{p}_{g,i}}\right)\right]+\mu s\left[h{q}_{D,k}^2{p}_{D,k}^2\left({\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{p}_{g,i}{p}_{g,j}}}\right)\right]\\ {}+\mu s\left[h{q}_{D,k}{p}_{D,k}^3\left({\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{p}_{g,i}{p}_{g,j}{p}_{g,m}}}}\right)\right]+\mu s\left[h{p}_{D,k}^4\left({\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{p}_{g,i}{p}_{g,j}{p}_{g,m}{p}_{g,n}}}}}\right)\right]-\mu \Big)\end{array} $$
(17)

Next, we deduce the formula for \( {\scriptscriptstyle \frac{\delta F}{\delta Z}}\left[0,h\right] \). Note that we set g = 0 because the formula for \( {\scriptscriptstyle \frac{\delta F}{\delta Z}}\left[g,h\right] \) is very cumbersome without setting it.

Lemma 1. The derivatives of F[g, h] is given by the formula

$$ \frac{\delta F}{\delta Z}\left[0,h\right]=F\left[0,h\right]\cdot {\displaystyle {\prod}_Z\cdot {\displaystyle \sum_{\wp \angle Z}{\displaystyle \prod_{W\in \wp }{d}_W\left[0,h\right]}}} $$
(18)

where the notation ‘Z’ is shorthand for “ partitions Z into cells W”, Z  = ∏ zZ λc(z),

$$ {d}_W\left[0,h\right]=\left\{\begin{array}{c}\hfill 1+\mu s\left[h{q}_{D,k}^3{p}_{D,k}\left({\displaystyle \sum_{i=1}^4{\ell}_{z_1,i}}\right)\right]\kern16em ifW=\left\{{z}_1\right\}\hfill \\ {}\hfill \mu s\left[h{q}_{D,k}^2{p}_{D,k}^2\left({\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{\displaystyle \sum_{a=1}^2{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^2{\ell}_{z_a,i}{\ell}_{z_b,j}}}}}\right)\right]\kern11em ifW=\left\{{z}_1,{z}_2\right\}\hfill \\ {}\hfill \mu s\left[h{q}_{D,k}{p}_{D,k}^3\left({\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{\displaystyle \sum_{a=1}^3{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^3{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^3{\ell}_{z_a,i}{\ell}_{z_b,j}{\ell}_{z_c,m}}}}}}}\right)\right]\kern5em ifW=\left\{{z}_1,{z}_2,{z}_3\right\}\hfill \\ {}\hfill \mu s\left[h{p}_{D,k}^4\left({\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{\displaystyle \sum_{a=1}^4{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^4{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^4{\displaystyle \sum_{\begin{array}{l}d=1\\ {}a\ne d\\ {}b\ne d\\ {}c\ne d\end{array}}^4{\ell}_{z_a,i}{\ell}_{z_b,j}{\ell}_{z_c,m}{\ell}_{z_d,n}}}}}}}}}\right)\right]\kern2em ifW=\left\{{z}_1,{z}_2,{z}_3,{z}_4\right\}\hfill \end{array}\right. $$
(19)

and where

$$ {\ell}_{z,i}(x)=\frac{g_{k,i}\left(z\Big|x\right)}{\lambda c(z)}\begin{array}{cc}\hfill \hfill & \hfill i=1,\cdots, 4.\hfill \end{array} $$
(20)

See Appendix 1 for the proof of Lemma 1.

Lemma 2. The multipath posterior PGFL G k|k [h] is given as

$$ {G}_{k\Big|k}\left[h\right]={F}_0\left[h\right]\cdot \frac{{\displaystyle {\sum}_{\wp \angle {Z}_k}{\displaystyle {\prod}_{W\in \wp }{d}_W\left[0,h\right]}}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W\left[0,1\right]}}} $$
(21)

where

$$ {F}_0\left[h\right]= \exp \left(\mu s\left[\left(h-1\right){q}_{D,k}^4\right]\right) $$
(22)

and d W [0, h] is the Eq. ( 19 ).

See Appendix 2 for the proof of Lemma 2.

According to Lemmas 1 and 2, we can obtain the MP-PHD filter update equation in Proposition 1 as follows:

Proposition 1. The update equation for the MP-PHD filter is

$$ {D}_{k\Big|k}\left(x\Big|{Z}^{(k)}\right)\cong {L}_{Z_k}\left(x\Big|{Z}^{\left(k-1\right)}\right)\cdot {D}_{k\Big|k-1}\left(x\Big|{Z}^{\left(k-1\right)}\right) $$
(23)

where

$$ {L}_{Z_k}\left(x\Big|{Z}^{\left(k-1\right)}\right)={q}_{D,k}^4(x)+{\displaystyle \sum_{\wp \angle {Z}_k}{\omega}_{\wp }}\cdot {\displaystyle \sum_{W\in \wp}\frac{b_W}{d_W}} $$
(24)

here

$$ {\omega}_{\wp }=\frac{{\displaystyle {\prod}_{W\in \wp }{d}_W}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W}}} $$
(25)

where

$$ {d}_W=\left\{\begin{array}{c}\hfill 1+{D}_{k\Big|k-1}\left[{q}_{D,k}^3{p}_{D,k}\left({\displaystyle \sum_{i=1}^4{\ell}_{z_1,i}}\right)\right]ifW=\left\{{z}_1\right\}\hfill \\ {}\hfill {D}_{k\Big|k-1}\left[{q}_{D,k}^2{p}_{D,k}^2\left({\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{\displaystyle \sum_{a=1}^2{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^2{\ell}_{z_a,i}{\ell}_{z_b,j}}}}}\right)\right]ifW=\left\{{z}_1,{z}_2\right\}\hfill \\ {}\hfill {D}_{k\Big|k-1}\left[{q}_{D,k}{p}_{D,k}^3\left({\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{\displaystyle \sum_{a=1}^3{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^3{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^3{\ell}_{z_a,i}{\ell}_{z_b,j}{\ell}_{z_c,m}}}}}}}\right)\right]ifW=\left\{{z}_1,{z}_2,{z}_3\right\}\hfill \\ {}\hfill {D}_{k\Big|k-1}\left[{p}_{D,k}^4\left({\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{\displaystyle \sum_{a=1}^4{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^4{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^4{\displaystyle \sum_{\begin{array}{l}d=1\\ {}a\ne d\\ {}b\ne d\\ {}c\ne d\end{array}}^4{\ell}_{z_a,i}{\ell}_{z_b,j}{\ell}_{z_c,m}{\ell}_{z_d,n}}}}}}}}}\right)\right]\kern08em ifW=\left\{{z}_1,{z}_2,{z}_3,{z}_4\right\}\hfill \end{array}\right. $$
(26)
$$ {b}_W=\left\{\begin{array}{c}\hfill {q}_{D,k}^3{p}_{D,k}\left({\displaystyle \sum_{i=1}^4{\ell}_{z_1,i}}\right)\kern27.5em ifW=\left\{{z}_1\right\}\hfill \\ {}\hfill {q}_{D,k}^2{p}_{D,k}^2\left({\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{\displaystyle \sum_{a=1}^2{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^2{\ell}_{z_a,i}{\ell}_{z_b,j}}}}}\right)\kern20.5em ifW=\left\{{z}_1,{z}_2\right\}\hfill \\ {}\hfill {q}_{D,k}{p}_{D,k}^3\left({\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{\displaystyle \sum_{a=1}^3{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^3{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^3{\ell}_{z_a,i}{\ell}_{z_b,j}{\ell}_{z_c,m}}}}}}}\right)\kern14.9em ifW=\left\{{z}_1,{z}_2,{z}_3\right\}\hfill \\ {}\hfill {p}_{D,k}^4\left({\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{\displaystyle \sum_{a=1}^4{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^4{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^4{\displaystyle \sum_{\begin{array}{l}d=1\\ {}a\ne d\\ {}b\ne d\\ {}c\ne d\end{array}}^4{\ell}_{z_a,i}{\ell}_{z_b,j}{\ell}_{z_c,m}{\ell}_{z_d,n}}}}}}}}}\right)ifW=\left\{{z}_1,{z}_2,{z}_3,{z}_4\right\}\hfill \end{array}\right. $$
(27)

See Appendix 3 for the proof of Proposition 1.

4 Gaussian mixture implementation of the MP-PHD filter

According to the GM implemented in the standard PHD filter in [25], a GM-PHD filter recursion can be derived for the multipath target tracking. Due to the MP-PHD filter prediction equation for the OTHR are identical to the standard PHD filter, the GM prediction equation of MP-PHD filter are the same as the standard GM-PHD filter, where the details can be found in [25]. Therefore, in this section, we only consider the GM update equation of the MP-PHD filter.

4.1 Gaussian mixture MP-PHD

To derive a closed-form solution to the update equation of the MP-PHD filter, we assume that each target follows a linear-Gaussian transition and measurement model, i.e.,

$$ {f}_{k\Big|k-1}\left(x\Big|\zeta \right)=N\left(x;{F}_{k\Big|k-1}{x}_{k\Big|k-1},{Q}_{k\Big|k-1}\right) $$
(28)
$$ {g}_{k,i}\left(z\Big|x\right)=N\left(z;{H}_{k,i}x,{R}_{k,i}\right)\begin{array}{cc}\hfill, \hfill & \hfill i\hfill \end{array}=1,\cdots, 4 $$
(29)

where N(; m, P) denotes a Gaussian density with mean m and covariance P. F is the state transition matrix, Q is the process noise covariance, H k,i is the observation matrix, and R k,i is the observation noise covariance. The survival and detection probability are assumed as a constant, i.e., p S,k (x) = p S,k , p D,k (x) = p D,k . To facilitate the derivation of the closed-form solution, we define an intermediate operator \( {G}_{k,z}^i \) by

$$ \left({G}_{k,{z}_a}^i\phi \right)(x)={g}_{k,i}\left({z}_a\Big|x\right)\phi (x) $$
(30)

If g k,i (z a |x) = N(z a ; H k,i x, R k,i ) and ϕ(x) = N(z; m ϕ , P ϕ ), then \( \left({G}_{k,{z}_a}^i\phi \right)(x) \) is a Gaussian density

$$ \left({G}_{k,{z}_a}^i\phi \right)(x)={\rho}_{k,i}\left({z}_a\right)N\left(x;{m}_{k,i},{P}_{k,i}\right) $$
(31)

where

$$ {\xi}_{k\Big|k-1,i}={H}_{k,i}{m}_{\phi } $$
(32)
$$ {S}_{k\Big|k-1,i}={R}_{k,i}+{H}_{k,i}{P}_{\phi }{H}_{k,i}^T $$
(33)
$$ {\rho}_{k,i}\left({z}_a\right)=N\left({z}_a;{\xi}_{k\Big|k-1,i},{S}_{k\Big|k-1,i}\right) $$
(34)
$$ {m}_{k,i}={m}_{\phi }+{K}_{k,i}\left({z}_a-{\xi}_{k\Big|k-1,i}\right) $$
(35)
$$ {P}_{k,i}=\left(I-{K}_{k,i}{H}_{k,i}\right){P}_{\phi } $$
(36)
$$ {K}_{k,i}={P}_{\phi }{H}_{k,i}^T{\left({S}_{k\Big|k-1,i}\right)}^{-1} $$
(37)

We suppose that the predicted PHD filter at time k − 1 has the following GM representation:

$$ \begin{array}{c}{D}_{k\Big|k-1}(x)={\displaystyle \sum_{l=1}^{J_{k\Big|k-1}}{\omega}_{k\Big|k-1}^{(l)}}\cdot N\left(x;{m}_{k\Big|k-1}^{(l)},{P}_{k\Big|k-1}^{(l)}\right)\\ {}={\displaystyle \sum_{l=1}^{J_{k\Big|k-1}}{\omega}_{k\Big|k-1}^{(l)}}\cdot {\phi}_{k\Big|k-1}^{(l)}(x)\end{array} $$
(38)

where \( {\omega}_{k\Big|k-1}^{(l)} \) is the weight of the lth component and J k|k − 1 is the predicted number of the components. Then, the update equation at time k can be expressed as a GM of

$$ {D}_{k\Big|k}(x)={D}_{k\Big|k}^{ND}(x)+{\displaystyle \sum_{\wp \angle {Z}_k}{\displaystyle \sum_{W\in \wp }{D}_{k\Big|k}^D\left(x,W\right)}} $$
(39)

The GM components \( {D}_{k\Big|k}^{ND}(x) \), handling the no detections cases, are given by

$$ {D}_{k\Big|k}^{ND}(x)={\displaystyle \sum_{l=1}^{J_{k\Big|k-1}}{\omega}_{k\Big|k}^{(l)}}\cdot N\left(x;{m}_{k\Big|k}^{(l)},{P}_{k\Big|k}^{(l)}\right) $$
(40)
$$ {\omega}_{k\Big|k}^{(l)}={q_{D,k}}^4{\omega}_{k\Big|k-1}^{(l)} $$
(41)
$$ {m}_{k\Big|k}^{(l)}={m}_{k\Big|k-1}^{(l)} $$
(42)
$$ {P}_{k\Big|k}^{(l)}={P}_{k\Big|k-1}^{(l)} $$
(43)

The GM components \( {D}_{k\Big|k}^D\left(x,W\right) \), handling detected target cases, are given by

$$ {D}_{k\Big|k}^D(x)={\displaystyle \sum_{l=1}^{J_{k\Big|k-1}}{\omega}_{k\Big|k}^{(l)}}\cdot N\left(x;{m}_{k\Big|k}^{(l)},{P}_{k\Big|k}^{(l)}\right) $$
(44)
$$ {\omega}_{k\Big|k}^{(l)}={\omega}_{\wp}\cdot {\left(\frac{b_W}{d_W}\right)}^{(l)}\cdot {\omega}_{k\Big|k-1}^{(l)} $$
(45)
$$ {\omega}_{\wp }=\frac{{\displaystyle {\prod}_{W\in \wp }{d}_W}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W}}} $$
(46)

when W = {z 1},

$$ {d}_W=1+{\displaystyle \sum_{l=1}^{J_{k\Big|k-1}}\frac{{q_{D,k}}^3{p}_{D,k}}{\lambda c\left({z}_1\right)}{\omega}_{k\Big|k-1}^{(l)}}{\displaystyle \sum_{i=1}^4\left({G}_{k,{z}_1}^i{\phi}_{k\Big|k-1}^{(l)}\right)(x)} $$
(47)
$$ {\left(\frac{b_W}{d_W}\right)}^{(l)}=\frac{{q_{D,k}}^3{p}_{D,k}{\displaystyle \sum_{i=1}^4\left({G}_{k,{z}_1}^i{\phi}_{k\Big|k-1}^{(l)}\right)(x)}}{\lambda c\left({z}_1\right)+{q_{D,k}}^3{p}_{D,k}{\displaystyle \sum_{r=1}^{J_{k\Big|k-1}}{\omega}_{k\Big|k-1}^{(r)}}{\displaystyle \sum_{i=1}^4\left({G}_{k,{z}_1}^i{\phi}_{k\Big|k-1}^{(r)}\right)(x)}} $$
(48)

when W = {z 1, z 2}

$$ {d}_W=\frac{{q_{D,k}}^2{p_{D,k}}^2}{\lambda c\left({z}_1\right)\lambda c\left({z}_2\right)}{\displaystyle \sum_{l=1}^{J_{k\Big|k-1}}{\omega}_{k\Big|k-1}^{(l)}{\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{\displaystyle \sum_{a=1}^2{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^2\left({G}_{k,{z}_a}^i\circ {G}_{k,{z}_b}^j{\phi}_{k\Big|k-1}^{(l)}\right)(x)}}}}} $$
(49)
$$ {\left(\frac{b_W}{d_W}\right)}^{(l)}=\frac{{\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{\displaystyle \sum_{a=1}^2{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^2\left({G}_{k,{z}_a}^i\circ {G}_{k,{z}_b}^j{\phi}_{k\Big|k-1}^{(l)}\right)(x)}}}}}{{\displaystyle \sum_{r=1}^{J_{k\Big|k-1}}{\omega}_{k\Big|k-1}^{(r)}{\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{\displaystyle \sum_{a=1}^2{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^2\left({G}_{k,{z}_a}^i\circ {G}_{k,{z}_b}^j{\phi}_{k\Big|k-1}^{(r)}\right)(x)}}}}}} $$
(50)

when W = {z 1, z 2, z 3}

$$ {d}_W=\frac{q_{D,k}{p_{D,k}}^3}{\lambda c\left({z}_1\right)\lambda c\left({z}_2\right)\lambda c\left({z}_3\right)}{\displaystyle \sum_{l=1}^{J_{k\Big|k-1}}{\omega}_{k\Big|k-1}^{(l)}{\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{\displaystyle \sum_{a=1}^3{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^3{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^3\left({G}_{k,{z}_a}^i\circ {G}_{k,{z}_b}^j\circ {G}_{k,{z}_c}^m{\phi}_{k\Big|k-1}^{(l)}\right)(x)}}}}}}} $$
(51)
$$ {\left(\frac{b_W}{d_W}\right)}^{(l)}=\frac{{\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{\displaystyle \sum_{a=1}^3{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^3{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^3\left({G}_{k,{z}_a}^i\circ {G}_{k,{z}_b}^j\circ {G}_{k,{z}_c}^m{\phi}_{k\Big|k-1}^{(l)}\right)(x)}}}}}}}{{\displaystyle \sum_{r=1}^{J_{k\Big|k-1}}{\omega}_{k\Big|k-1}^{(r)}{\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{\displaystyle \sum_{a=1}^3{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^3{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^3\left({G}_{k,{z}_a}^i\circ {G}_{k,{z}_b}^j\circ {G}_{k,{z}_c}^m{\phi}_{k\Big|k-1}^{(r)}\right)(x)}}}}}}}} $$
(52)

when W = {z 1, z 2, z 3, z 4}

$$ {d}_W=\frac{{p_{D,k}}^4}{\lambda c\left({z}_1\right)\lambda c\left({z}_2\right)\lambda c\left({z}_3\right)\lambda c\left({z}_4\right)}{\displaystyle \sum_{l=1}^{J_{k\Big|k-1}}{\omega}_{k\Big|k-1}^{(l)}{\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{\displaystyle \sum_{a=1}^4{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^4{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^4{\displaystyle \sum_{\begin{array}{l}d=1\\ {}a\ne d\\ {}b\ne d\\ {}c\ne d\end{array}}^4\left({G}_{k,{z}_a}^i\circ {G}_{k,{z}_b}^j\circ {G}_{k,{z}_c}^m\circ {G}_{k,{z}_d}^n{\phi}_{k\Big|k-1}^{(l)}\right)(x)}}}}}}}}} $$
(53)
$$ {\left(\frac{b_W}{d_W}\right)}^{(l)}=\frac{{\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{\displaystyle \sum_{a=1}^4{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^4{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^4{\displaystyle \sum_{\begin{array}{l}d=1\\ {}a\ne d\\ {}b\ne d\\ {}c\ne d\end{array}}^4\left({G}_{k,{z}_a}^i\circ {G}_{k,{z}_b}^j\circ {G}_{k,{z}_c}^m\circ {G}_{k,{z}_d}^n{\phi}_{k\Big|k-1}^{(l)}\right)(x)}}}}}}}}}{{\displaystyle \sum_{r=1}^{J_{k\Big|k-1}}{\omega}_{k\Big|k-1}^{(r)}{\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{\displaystyle \sum_{a=1}^4{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^4{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^4{\displaystyle \sum_{\begin{array}{l}d=1\\ {}a\ne d\\ {}b\ne d\\ {}c\ne d\end{array}}^4\left({G}_{k,{z}_a}^i\circ {G}_{k,{z}_b}^j\circ {G}_{k,{z}_c}^m\circ {G}_{k,{z}_d}^n{\phi}_{k\Big|k-1}^{(r)}\right)(x)}}}}}}}}}} $$
(54)

where denotes a composition.

4.2 Implementation issues

  1. 1)

    Nonlinear measurement model: The above Gaussian mixture MP-PHD filter can be extended to a nonlinear measurement model using the EKF. Due to that the dynamic model in the OTHR is linear, we only apply EKF to deal with the nonlinear measurement model in the update step of MP-PHD filter. Use the approximations (55) and (56) in place of the originals (32) and (33), and use the linearizations in (57) to calculate (36) and (37)

$$ {\xi}_{k\Big|k-1,i}={h}_{k,i}\left({m}_{\phi },0\right) $$
(55)
$$ {S}_{k\Big|k-1,i}={U}_{k,i}{R}_{k,i}{U}_{k,i}^T+{H}_{k,i}{P}_{\phi }{H}_{k,i}^T $$
(56)

where

$$ {H}_{k,i}=\frac{\partial {h}_{k,i}\left(x,0\right)}{\partial x}\left|{}_{x={m}_{\phi }}\right. $$
(57)
  1. 2)

    Managing mixture components: Similar to the standard GM-PHD filter, the techniques of merging and pruning must be used to reduce the exponential growth of the number of Gaussian components, which the merging and pruning methods specific to Gaussian components can be found in [25].

  2. 3)

    Computational complexity: As is shown in Eq. (24), the MP-PHD filter requires all partitions of the current measurements for updating, Z k denotes that is one partition of the measurement set Z k and W denotes that W is one cell of . Note that W does not include null set. For instance, Z k  = {z 1, z 2, z 3} can be partitioned as follows [1921]

$$ {\wp}_1=\left\{\left\{{z}_1\right\},\left\{{z}_2\right\},\left\{{z}_3\right\}\right\},\ {\wp}_2=\left\{\left\{{z}_1,{z}_2\right\},\left\{{z}_3\right\}\right\},\ {\wp}_3=\left\{\left\{{z}_1\right\},\left\{{z}_2,{z}_3\right\}\right\}, $$
$$ {\wp}_4=\left\{\left\{{z}_1,{z}_3\right\},\left\{{z}_2\right\}\right\},\ {\wp}_5=\left\{\left\{{z}_1,{z}_2,{z}_3\right\}\right\}, $$

whereas with the number of measurements growing, the number of partitions grows very large, which is computationally infeasible.

From the above analysis, we can see that the computational load required by the proposed MP-PHD filter is similar to that of the extended targets PHD filter [1822]. Note that some methods, such as K-means++ method [19] and spectral clustering [20], have been suggested for the implementation of the extended targets PHD filter to reduce the number of partitions. However, these methods that applied in the extended targets are measurement dependent. Since the measurement function of each measurement mode can be very different in OTHR. These methods may fail with the multipath tracking problem in OTHR. To limit the number of partitions in the MP-PHD filter, some approximations are necessary in this paper. First, we can use the gating technology with multiple validation gates to reduce the number of measurements for each propagation model before the update step at each time [11], and then we can use the similar method of “effective binary partition” that was implemented in [26] to reduce the number of partitions. For space considerations which are not the focus of this paper, we omit the details. Note that we used only a simple partition technique in our implementation, and the development of more efficient implementations will be a future task.

5 Numerical simulations

The performance of the proposed MP-PHD filter will be evaluated through two numerical examples, including the single target and multitarget tracking in this section. And we use the optimal subpattern assignment (OSPA) metric to evaluate the performance of the proposed algorithm [27].

5.1 Experiment 1

In this subsection, we have designed an experiment to compare the MP-PHD filter with the standard single-path PHD filter in a single target tracking scenario. Note here that the standard single-path PHD filter uses one specific mode (mode EE used in this experiment) to track the targets, and we use the gating technology to reduce the number of measurements for propagation model EE before the update step at each time. The two filters are simulated with the same environment in OTHR. It is assumed that a single non-maneuvering target is in the presence of clutter. Clutter is generally modeled as a Poisson RFS with the intensity function κ k (y) = λVu(y), where u() indicates the uniform density over the region [1000, 1400]km × [0.069813, 0.17453]rad, V = 41.8868 km  rad is the “volume” of this surveillance region, and λ = 4.7743(km  rad)− 1 is the average number of the clutter returns per unit volume. This translates to 200 clutter measurements per scan. And each propagation model has the same detection probability and target survive probability that are set to p D,k  = 0.6 and p S,k  = 0.95. To keep computations tractable, the Gaussian component pruning is applied in this numerical example with a pruning threshold τ = 10− 5, a merging threshold U = 20, and a maximum number of Gaussian components J max = 100. It is assumed that a single target moving in the surveillance region with the initial target state x 0 = (1100 km, 0.15 km/s, 0.10472 rad, 8.72665e−05 rad/s). Other simulation parameters are defined as the same as [11]. The true target trajectory and OTHR multipath detections without clutter under p D,k  = 0.6 are shown in Fig. 2. It is shown that the receiver can obtain more than one measurement from the same target at one scan.

Fig. 2
figure 2

The scenario of true target trajectory and OTHR multipath detections. a True target trajectory in ground coordinates. b OTHR multipath detections without clutter in slant coordinates

The target tracking results of the MP-PHD filter and the standard single-path PHD filter for one simulation are shown in Fig. 3. As the position estimates are shown in (a) and (b), we can see that the MP-PHD filter provides more accurate position estimates for almost all the time. Moreover, for the simulation results, shown in (c) and (d), using the MP-PHD, there is only one false estimated number at time 8.

Fig. 3
figure 3

Target tracking produced by a MP-PHD filter and a single-path PHD filter. a Standard single-path PHD filter. b MP-PHD filter. c Standard single-path PHD filter. d MP-PHD filter

5.2 Experiment 2

In this scenario, we assumed that the number of targets is constant. The performance of the GM-MP-PHD filter is compared with the MD-JPDA filter. Two filters are simulated with the same environment in OTHR tracking system, and we also use the gating technology with multiple validation gates to reduce the number of measurements in the MD-JPDA filter. The experiment goes on for 600 s (sampling period T = 20 s), and we assumed that two targets appear with initial state target 1 x 1 = (1130 km, 0.12 km/s, 0.10472 rad, 8.72665e−05 rad/s) and target 2 x 2 =(1125 km, −0.1 km/s, 0.11472 rad, 7.92665e−05 rad/s), respectively. And each propagation model has the same target survive probability and detection probability that are set to p S,k  = 0.95 and p D,k  = 0.6. Other parameters are adopted as the same as the experiment 1.

The position root mean squared error (RMSE) of the two filters are shown in Fig. 4. It is shown that the proposed MP-PHD filter has a notable performance gain over the MD-JPDA filter. This is due to the fact that the proposed MP-PHD filter encapsulates all the information about the targets by enumerating all of the association between propagation paths and measurements. However, only a subset of the associations (not all the permutations of the propagation paths) is used in the MD-JPDA filter, this is the reason degrades the performance of MD-JPDA filter.

Fig. 4
figure 4

Position RMSE for a MP-PHD filter and a MD-JPDA filter with OTHR data

To evaluate the computational time of the proposed algorithm, the averaged time is computed in MATLAB7.1 on an Intel CORE i5 CPU computer with 2 GB of RAM. The proposed GM-MP-PHD consumed approximately 8.6 s per MC run and the MD-JPDA consumed approximately 1.1 s. Although the proposed MP-PHD filter has high computational time, it has better performance than the MD-JPDA filter.

5.3 Experiment 3

To validate the tracking performance of MP-PHD filter in multipath multitarget tracking, in this experiment, we have designed a multiple track scenario including three non-maneuvering targets where the target number is changing in OTHR. The experiment goes on for 800 s with the sampling period T = 20 s, and it is assumed that target 1 and target 2 appear in the course of the whole experiment with initial state x 1 = (1100 km, 0.15 km/s, 0.10472 rad, 8.72665e−05 rad/s) and x 2 = (1170 km, −0.14 km/s, 0.11472 rad, 7.72665e−05 rad/s), respectively, target 3 appears at t = 180 s and disappears at t = 480 s with the initial states x 3 = (1170 km, −0.05 km/s, 0.15701 rad, −8.72665e−05 rad/s). In this experiment, there are no spawning targets for simplification. Other simulation conditions are the same as the experiment 1.

  1. 1)

    One simulation: The true target trajectory is shown in Fig. 5, and the position estimates of the MP-PHD filter for one simulation compared with the standard single-path PHD filter are shown in Fig. 6, the estimation of target number is shown in Fig. 7. As shown in Figs. 6 and 7, the MP-PHD filter can accurately estimate the target state and the target number of the multipath multitarget tracking in OTHR. In comparison, the simulation result shows that the MP-PHD filter provides more accurate position estimates for almost all the times.

    Fig. 5
    figure 5

    The true target trajectory

    Fig. 6
    figure 6

    The true and estimation target trajectory (p D,k  = 0.6)

    Fig. 7
    figure 7

    The true and estimated target number (p D,k  = 0.6)

  2. 2)

    100 Monte Carlo simulations: We compare the MP-PHD filter with the standard single-path PHD filter over 100 Monte Carlo simulations with different detection probability p D,k  = 0.6 and p D,k  = 0.99, respectively. We use the OSPA metric to evaluate the performance, the parameters of OSPA distance are set to the cutoff parameter c = 15 km and the order parameter p = 2 (see [27] for more details).

Figure 8 shows the average of OSPA distance comparison between the MP-PHD filter and the standard single-path PHD filter, and the average of estimated number is shown in Fig. 9. The simulation results demonstrate that the tracking performance of the MP-PHD filter surpasses the performance of the standard single-path PHD filter. Moreover, it can accurately estimate the target state and the target number with low detection probability under the multitarget tracking circumstance. As shown in Figs. 8 and 9, the standard single-path PHD filter cannot deal with the multitarget tracking problem effectively under low detection probability. This is due to the fact that the MP-PHD filter use four path measurements to track the targets, and then, it can more effectively use the measurements from the targets under low detection probability.

Fig. 8
figure 8

The average of OSPA distance

Fig. 9
figure 9

Average of estimated target number

6 Conclusions

In this paper, the MP-PHD filter, which tracks the multipath multitarget in OTHR, has been proposed. The MP-PHD filter is based on the theory of FISST, the FISST is employed to derive the update equation, and then a closed-form solution of MP-PHD filter in the form of GM has been proposed. The EKF is used to deal with the nonlinear problem of the measurement model. Simulation results show that the proposed MP-PHD filter has a notable performance gain over the MD-JPDA filter.

Similar to the standard PHD filter, the MP-PHD filter should be improved for some problems. For example, as the number of targets in the surveillance region is large, the estimated target number is unreliable. To solve this problem, a modified PHD filter name, cardinalized PHD filter, is proposed by Mahlar [28, 29]. Our future work could derive a CPHD filter which can be applied to the multitarget tracking in OTHR.

References

  1. DJ Percival and KAB White, Multipath track fusion for over-the-horizon radar, in: Proceedings of SPIE Signal Data Process of Small Targets, San Diego, CA, United States, July 27, 1997, pp. 363-374.

  2. DJ Percival and KAB White, Multihypothesis fusion of multipath over-the-horizon radar tracks, in: Proceedings of SPIE Signal Data Process of Small Targets, Orlando, FL, April 13, 1998, pp. 440-451.

  3. MG Rutten and D J Percival, Joint Ionospheric and Target State Estimation for Multipath OTHR Track Fusion, in: Proceedings of SPIE Signal Data Process of Small Targets, vol. 4473, San Diego, CA, USA , July 29, 2001, pp. 118-129.

  4. GW Pulford, OTHR multipath tracking with uncertain coordinate registration. IEEE Transactions on Aerospace and Electronic Systems 40(41), 38–56 (2004)

    Article  Google Scholar 

  5. H Liu, Q Pan, Y Liang, Y Cheng, M Cheng, Comments on “multipath data association tracker for over-the-horizon radar”. IEEE Transactions on Aerospace and Electronic Systems 41(3), 1147–1148 (2005)

    Article  Google Scholar 

  6. MG Rutten, S Maskell, M Briers and NJ Gordon, Multipath track association for over-the-horizon radar using Lagrangian relaxation, in: Proceedings of SPIE Signal Data Process of Small Targets, Orlando, FL, April 12, 2004, pp. 452-463.

  7. Y Bar-Shalom, E Tse, Tracking in a cluttered environment with probabilistic data association. Automatica 11(5), 451–460 (1975)

    Article  MATH  Google Scholar 

  8. SB Colegrove, AW Davis, JK Ayliffe, Track initiation and nearest neighbours incorporated into probabilistic data association. Journal Elect. Electronics Eng 6(3), 191–198 (1986)

    Google Scholar 

  9. T Kirubarajan, Y Bar-Shalom, Probabilistic data association techniques for target tracking in clutter. Proceedings of the IEEE 92(3), 536–557 (2004)

    Article  Google Scholar 

  10. H Liu, Y Liang, Q Pan and Y Cheng, A multipath viterbi data association algorithm for OTHR, in: Proceedings of International Conference on Radar, Shanghai, China, 2006, pp. 1–4

  11. GW Pulford, RJ Evans, A multipath data association tracker for over-the-horizon radar. IEEE Transactions on Aerospace and Electronic Systems 34(4), 1165–1183 (1998)

    Article  Google Scholar 

  12. SS Blackman, Multiple hypothesis tracking for multiple target tracking. IEEE Aerospace and Electronic Systems Magazine 19(1), 5–18 (2004)

    Article  Google Scholar 

  13. HAP Blom, EA Bloem, Interacting multiplemodel joint probabilistic data association, avoiding track coalescence, in Proceedings of the 41th IEEE Conference on Decision Control, vol. 3, 2002, pp. 3408–3415

    Google Scholar 

  14. D Musicki, R Evans, Joint integrated probabilistic data association: JIPDA. IEEE Transactions on Aerospace and Electronic Systems 40(3), 1093–1099 (2004)

    Article  Google Scholar 

  15. R Mahler, Multi-target Bayes filtering via first-order multi-target moments. IEEE Transactions on Aerospace and Electronic Systems 39(4), 1152–1178 (2003)

    Article  Google Scholar 

  16. B Habtemariam, R Tharmarasa, T Thayaparan et al., A multiple-detection joint probabilistic data association filter. IEEE Journal Selected Topics Signal Processing 7(3), 461–471 (2013)

    Article  Google Scholar 

  17. R Mahler, “Second-generation” PHD/CPHD filters and multitarget calculus, in: Proceedings of SPIE Signal Data Process of Small Targets, San Diego, CA, August 02 , 2009.

  18. R Mahler, PHD filters for nonstandard target, I: extended targets, in: Proceedings of the 12th International Conference on Information Fusion, Seattle, WA, July 2009, pp. 915–921

  19. K Granstrom, C Lundquist, U Orguner, A Gaussian mixture PHD filter for extended target tracking, in: Proceedings of the 13th International Conference on Information Fusion, Edinburgh, UKL, July 2010, pp. 1–8

  20. Yang J, Liu F, Ge H, et al., Multiple extended target tracking algorithm based on GM-PHD filter and spectral clustering, EURASIP Journal on Advances in Signal Processing 2014:117

  21. Y Zhang, H Ji, Gaussian mixture reduction based on fuzzy ART for extended target tracking. Signal Process 97, 232–241 (2014)

    Article  Google Scholar 

  22. R Mahler, PHD filters for nonstandard targets, II: unresolved targets, in: Proceedings of the 12th International Conference on Information Fusion, Seattle, WA, USA, July 2009 pp. 922–929

  23. R Mahler, A El-Fallah, CPHD and PHD filters for unknown backgrounds, III: tractable multitarget filtering in dynamic clutter, in Proceedings of SPIE, vol. 7698, 2010

    Google Scholar 

  24. R Mahler, CPHD filters for superpositional sensors, in Proceedings of SPIE, vol. 7445, 2009

    Google Scholar 

  25. B-N Vo, W-K Ma, The Gaussian mixture probability hypothesis density filter. IEEE Transactions on Signal Processing 54(11), 4091–4104 (2006)

    Article  Google Scholar 

  26. J Xu, F Huang, Z Huang, The multi-sensor PHD filter Analytic implementation via Gaussian mixture and effective binary partition, in: Proceedings of the 16th International Conference on Information Fusion, Istanbul, TUR, July 2013 pp. 945–952

  27. D Schuhmacher, B-T Vo, B-N Vo, A consistent metric for performance evaluation of multi-object filters. IEEE Transactions on Signal Processing 56(8), 3447–3457 (2008)

    Article  MathSciNet  Google Scholar 

  28. B-T Vo, B-N Vo, A Cantoni, Analytic implementations of the cardinalized probability hypothesis density filter. IEEE Transactions on Signal Processing 55(7), 3553–3567 (2007)

    Article  MathSciNet  Google Scholar 

  29. R Mahler, PHD filters of higher order in target number. IEEE Transactions on Aerospace and Electronic Systems 43(4), 1523–1543 (2007)

    Article  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yong Qin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Appendices

Appendices

1.1 Appendix 1: Proof of Lemma 1

This lemma is proved by mathematical induction. For the initial induction step, assume Z = {z 1}. In this case there is only one partition,  = {{z 1}}, and we get

$$ \begin{array}{c}\frac{\delta F}{\delta {z}_1}\left[g,h\right]\left|{}_{g=0}\right.=F\left[g,h\right]\lambda c\left({z}_1\right)\Big(1+\mu s\left[h{q}_{D,k}^3{p}_{D,k}\left({\displaystyle \sum_{i=1}^4{\ell}_{z_1,i}}\right)\right]+\mu s\left[h{q}_{D,k}^2{p}_{D,k}^2\left({\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4\left({\ell}_{z_1,i}{p}_{g,j}+{p}_{g,i}{\ell}_{z_1,j}\right)}}\right)\right]\\ {}+\mu s\left[h{q}_{D,k}{p}_{D,k}^3\left({\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4\left({\ell}_{z_1,i}{p}_{g,j}{p}_{g,m}+{p}_{g,i}{\ell}_{z_1,j}{p}_{g,m}+{p}_{g,i}{p}_{g,j}{\ell}_{z_1,m}\right)}}}\right)\right]\left)+\mu s\right[h{p}_{D,k}^4\\ {}\left({\ell}_{z_1,1}{p}_{g,2}{p}_{g,3}{p}_{g,4}+{p}_{g,1}{\ell}_{z_1,2}{p}_{g,3}{p}_{g,4}+{p}_{g,1}{p}_{g,2}{\ell}_{z_1,3}{p}_{g,4}+{p}_{g,1}{p}_{g,2}{p}_{g,3}{\ell}_{z_1,4}\right)\left]\right)\left|{}_{g=0}\right.\\ {}=F\left[0,h\right]\lambda c\left({z}_1\right)\left(1+\mu s\left[h{q}_{D,k}^3{p}_{D,k}\left({\displaystyle \sum_{i=1}^4{\ell}_{z_1,i}}\right)\right]\right)\\ {}=F\left[0,h\right]\lambda c\left({z}_1\right){d}_{\left\{{z}_1\right\}}\left[0,h\right]\end{array} $$
(58)

Now, assume that we have established Eq. (18) for Z = {z 1, , z m } with |Z| = m > 1. We are to establish Eq. (18) for Z = {z 1, , z m , z m + 1} and |Z| = m + 1. Using the product rule for functional derivatives, we get

$$ \begin{array}{c}\frac{\delta^2F}{\delta Z\delta {z}_{m+1}}\left[g,h\right]\left|{}_{g=0}\right.=\frac{\delta F}{\delta {z}_{m+1}}\left[g,h\right]\left|{}_{g=0}\right.\cdot {\displaystyle {\prod}_Z\cdot {\displaystyle \sum_{\wp \angle Z}\left({\displaystyle \prod_{W\in \wp }{d}_W\left[0,h\right]}\right)}}\\ {}+F\left[0,h\right]\cdot {\displaystyle {\prod}_Z\cdot {\displaystyle \sum_{\wp \angle Z}\frac{\delta }{\delta {z}_{m+1}}\left({\displaystyle \prod_{W\in \wp }{d}_W\left[g,h\right]}\right)}}\left|{}_{g=0}\right.\end{array} $$
(59)

On the one hand,

$$ \frac{\delta F}{\delta {z}_{m+1}}\left[g,h\right]\left|{}_{g=0}\right.=F\left[0,h\right]\lambda c\left({z}_{m+1}\right){d}_{\left\{{z}_{m+1}\right\}}\left[0,h\right] $$
(60)

On the other hand,

$$ \frac{\delta }{\delta {z}_{m+1}}{\displaystyle \prod_{W\in \wp }{d}_W\left[g,h\right]}\left|{}_{g=0}\right.=\left({\displaystyle \prod_{W\in \wp }{d}_W\left[0,h\right]}\right)\cdot {\displaystyle \sum_{W\in \wp}\frac{1}{d_W\left[0,h\right]}}\frac{\delta {d}_W}{\delta {z}_{m+1}}\left[g,h\right]\left|{}_{g=0}\right. $$
(61)

where

$$ \frac{\frac{\delta {d}_W}{\delta {z}_{m+1}}\left[g,h\right]\left|{}_{g=0}\right.}{\lambda c\left({z}_{m+1}\right)}=\left\{\begin{array}{c}\hfill \mu s\left[h{q}_{D,k}^2{p}_{D,k}^2\left({\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{\displaystyle \sum_{a=1}^2{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^2{\ell}_{z_a,i}{\ell}_{z_b,j}}}}}\right)\right]\kern11em ifW=\left\{{z}_2\right\}\hfill \\ {}\hfill \mu s\left[h{q}_{D,k}{p}_{D,k}^3\left({\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{\displaystyle \sum_{a=1}^3{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^3{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^3{\ell}_{z_a,i}{\ell}_{z_b,j}{\ell}_{z_c,m}}}}}}}\right)\right]\kern5em ifW=\left\{{z}_2,{z}_3\right\}\hfill \\ {}\hfill \mu s\left[h{p}_{D,k}^4\left({\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{\displaystyle \sum_{a=1}^4{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^4{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^4{\displaystyle \sum_{\begin{array}{l}d=1\\ {}a\ne d\\ {}b\ne d\\ {}c\ne d\end{array}}^4{\ell}_{z_a,i}{\ell}_{z_b,j}{\ell}_{z_c,m}{\ell}_{z_d,n}}}}}}}}}\right)\right]\kern2em ifW=\left\{{z}_2,{z}_3,{z}_4\right\}\hfill \end{array}\right. $$
(62)

In the above equation, we use z 1 to replace z m + 1 for simplicity. And it is easy to find that \( \frac{\delta F}{\delta {z}_{m+1}}\left[g,h\right]\left|{}_{g=0}\right.=\lambda c\left({z}_{m+1}\right){d}_{W\cup \left\{{z}_{m+1}\right\}}\left[0,h\right] \). Putting this all together, Eq. (59) can be written as

$$ \begin{array}{c}\frac{\delta^2F}{\delta Z\delta {z}_{m+1}}\left[g,h\right]\left|{}_{g=0}\right.=F\left[0,h\right]\cdot {\displaystyle {\prod}_{Z\cup \left\{{z}_{m+1}\right\}}\cdot {\displaystyle \sum_{\wp \angle Z}\left({\displaystyle \prod_{W\in \wp \cup \left\{\left\{{z}_{m+1}\right\}\right\}}{d}_W\left[0,h\right]}\right)}}\\ {}+F\left[0,h\right]\cdot {\displaystyle {\prod}_{Z\cup \left\{{z}_{m+1}\right\}}\cdot {\displaystyle \sum_{\wp \angle Z}\left({\displaystyle \prod_{W\in \wp }{d}_W\left[0,h\right]}\right)}}\cdot \left({\displaystyle \sum_{W\in \wp}\frac{d_{W\cup \left\{{z}_{m+1}\right\}}\left[0,h\right]}{d_W\left[0,h\right]}}\right)\end{array} $$
(63)

Note that all partitions of the Z {z m + 1} have the following forms. First, take a partition of the set Z and add the cell {z m + 1} to obtain a new partition 1(z m + 1) =  {{z m + 1}} of Z {z m + 1}. This action is what is denoted mathematically in the first line of Eq. (63). Secondly, remove a cell W from the partition and replace it with W {z m + 1} of a new partition 2(z m + 1, W) = {W {z m + 1}}  V − {W}{V} of Z {z m + 1}. This action is what is denoted by the product

$$ \left({\displaystyle \prod_{W\in \wp }{d}_W\left[0,h\right]}\right)\cdot \left({\displaystyle \sum_{W\in \wp}\frac{d_{W\cup \left\{{z}_{m+1}\right\}}\left[0,h\right]}{d_W\left[0,h\right]}}\right) $$
(64)

in the second line of Eq. (63). Consequently, Eq. (63) becomes

$$ \frac{\delta^2F}{\delta Z\delta {z}_{m+1}}\left[g,h\right]\left|{}_{g=0}\right.=F\left[0,h\right]\cdot {\displaystyle {\prod}_{Z\cup \left\{{z}_{m+1}\right\}}\cdot {\displaystyle \sum_{\wp \angle Z\cup \left\{{z}_{m+1}\right\}}\left({\displaystyle \prod_{W\in \wp }{d}_W\left[0,h\right]}\right)}} $$
(65)

and this ends the inductive step.

1.2 Appendix 2: Proof of Lemma 2

From Eq. (10), we get

$$ \begin{array}{c}{G}_{k\Big|k}\left[h\right]=\frac{\frac{\delta F}{\delta {Z}_k}\left[0,h\right]}{\frac{\delta F}{\delta {Z}_k}\left[0,1\right]}\\ {}=\frac{F\left[0,h\right]}{F\left[0,1\right]}\cdot \frac{{\displaystyle {\sum}_{\wp \angle {Z}_k}{\displaystyle {\prod}_{W\in \wp }{d}_W\left[0,h\right]}}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W\left[0,1\right]}}}\\ {}= \exp \left(\mu s\left[\left(h-1\right){q}_{D,k}^4\right]\right)\cdot \frac{{\displaystyle {\sum}_{\wp \angle {Z}_k}{\displaystyle {\prod}_{W\in \wp }{d}_W\left[0,h\right]}}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W\left[0,1\right]}}}\end{array} $$
(66)

1.3 Appendix 3: Proof of proposition 1

The posterior multipath PHD filter can be calculated as

$$ {D}_{k\Big|k}(x)=\frac{\delta {G}_{k\Big|k}}{\delta x}\left[h\right]\left|{}_{h=1}\right. $$
(67)

Thus, from Eq. (21), we must first determine

$$ \frac{\delta {G}_{k\Big|k}}{\delta x}\left[h\right]=\frac{\delta {F}_0}{\delta x}\left[h\right]\cdot \frac{{\displaystyle {\sum}_{\wp \angle {Z}_k}{\displaystyle {\prod}_{W\in \wp }{d}_W\left[0,h\right]}}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W\left[0,1\right]}}}+{F}_0\left[h\right]\cdot \frac{{\displaystyle {\sum}_{\wp \angle {Z}_k}\frac{\delta }{\delta x}{\displaystyle {\prod}_{W\in \wp }{d}_W\left[0,h\right]}}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W\left[0,1\right]}}} $$
(68)

On the one hand, use Eq. (22) and get

$$ \frac{\delta {F}_0}{\delta x}\left[h\right]={F}_0\left[h\right]\cdot {D}_{k\Big|k-1}(x)\cdot {q}_{D,k}^4(x) $$
(69)

On the other hand,

$$ \frac{\delta }{\delta x}{\displaystyle \prod_{W\in \wp }{d}_W\left[0,h\right]}={\displaystyle \prod_{W\in \wp }{d}_W\left[0,h\right]}{\displaystyle \sum_{W\in \wp}\frac{\frac{\delta {d}_W}{\delta x}\left[0,h\right]}{d_W\left[0,h\right]}} $$
(70)

where from Eq. (19), we have

$$ \frac{\frac{\delta {d}_W}{\delta x}\left[0,h\right]}{D_{k\Big|k-1}(x)}={b}_W=\left\{\begin{array}{c}\hfill {q}_{D,k}^3(x){p}_{D,k}(x)\left({\displaystyle \sum_{i=1}^4{\ell}_{z_1,i}}(x)\right)i\kern21.1em fW=\left\{{z}_1\right\}\hfill \\ {}\hfill {q}_{D,k}^2(x){p}_{D,k}^2(x)\left({\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{\displaystyle \sum_{a=1}^2{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^2{\ell}_{z_a,i}(x){\ell}_{z_b,j}(x)}}}}\right)\kern13.1em ifW=\left\{{z}_1,{z}_2\right\}\hfill \\ {}\hfill \begin{array}{c}\hfill {q}_{D,k}(x){p}_{D,k}^3(x)\Big({\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{\displaystyle \sum_{a=1}^3{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^3{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^3{\ell}_{z_a,i}(x)}}}}}}\hfill \\ {}\hfill {\ell}_{z_b,j}(x){\ell}_{z_c,m}(x)\Big)\hfill \end{array}\kern6em ifW=\left\{{z}_1,{z}_2,{z}_3\right\}\hfill \\ {}\hfill \begin{array}{c}\hfill {p}_{D,k}^4(x)\Big({\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{\displaystyle \sum_{a=1}^4{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^4{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^4{\displaystyle \sum_{\begin{array}{l}d=1\\ {}a\ne d\\ {}b\ne d\\ {}c\ne d\end{array}}^4{\ell}_{z_a,i}(x)}}}}}}}}\hfill \\ {}\hfill {\ell}_{z_b,j}(x){\ell}_{z_c,m}(x){\ell}_{z_d,n}(x)\Big)\hfill \end{array}\kern3em ifW=\left\{{z}_1,{z}_2,{z}_3,{z}_4\right\}\hfill \end{array}\right. $$
(71)

Substituting Eqs. (69) and (71) into Eq. (68) yields

$$ \frac{\frac{\delta {G}_{k\Big|k}}{\delta x}\left[h\right]}{F_0\left[h\right]\cdot {D}_{k\Big|k-1}(x)}={q}_{D,k}^4(x)\cdot \frac{{\displaystyle {\sum}_{\wp \angle {Z}_k}{\displaystyle {\prod}_{W\in \wp }{d}_W\left[0,h\right]}}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W\left[0,1\right]}}}+\frac{{\displaystyle {\sum}_{\wp \angle {Z}_k}{\displaystyle {\prod}_{W\in \wp }{d}_W\left[0,h\right]}{\displaystyle {\sum}_{W\in \wp}\frac{b_W}{d_W\left[0,h\right]}}}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W\left[0,1\right]}}} $$
(72)

Setting h = 1 in the above equation then yields

$$ \begin{array}{c}\frac{D_{k\Big|k}(x)}{D_{k\Big|k-1}(x)}={q}_{D,k}^4(x)+\frac{{\displaystyle {\sum}_{\wp \angle {Z}_k}{\displaystyle {\prod}_{W\in \wp }{d}_W\left[0,1\right]}{\displaystyle {\sum}_{W\in \wp}\frac{b_W}{d_W\left[0,1\right]}}}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W\left[0,1\right]}}}\\ {}={q}_{D,k}^4(x)+{\displaystyle \sum_{\wp \angle {Z}_k}\frac{{\displaystyle {\prod}_{W\in \wp }{d}_W\left[0,1\right]}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W\left[0,1\right]}}}}\cdot {\displaystyle \sum_{W\in \wp}\frac{b_W}{d_W\left[0,1\right]}}\\ {}={q}_{D,k}^4(x)+{\displaystyle \sum_{\wp \angle {Z}_k}{\omega}_{\wp }}\cdot {\displaystyle \sum_{W\in \wp}\frac{b_W}{d_W}}\end{array} $$
(73)

where

$$ {\omega}_{\wp }=\frac{{\displaystyle {\prod}_{W\in \wp }{d}_W}}{{\displaystyle {\sum}_{\wp \hbox{'}\angle {Z}_k}{\displaystyle {\prod}_{W\in \wp \hbox{'}}{d}_W}}} $$
(74)
$$ {d}_W={d}_W\left[0,1\right]=\left\{\begin{array}{c}\hfill 1+{D}_{k\Big|k-1}\left[{q}_{D,k}^3{p}_{D,k}\left({\displaystyle \sum_{i=1}^4{\ell}_{z_1,i}}\right)\right]\kern16em ifW=\left\{{z}_1\right\}\hfill \\ {}\hfill {D}_{k\Big|k-1}\left[{q}_{D,k}^2{p}_{D,k}^2\left({\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=i+1}^4{\displaystyle \sum_{a=1}^2{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^2{\ell}_{z_a,i}{\ell}_{z_b,j}}}}}\right)\right]\kern11em ifW=\left\{{z}_1,{z}_2\right\}\hfill \\ {}\hfill {D}_{k\Big|k-1}\left[{q}_{D,k}{p}_{D,k}^3\left({\displaystyle \sum_{i=1}^2{\displaystyle \sum_{j=i+1}^3{\displaystyle \sum_{m=j+1}^4{\displaystyle \sum_{a=1}^3{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^3{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^3{\ell}_{z_a,i}{\ell}_{z_b,j}{\ell}_{z_c,m}}}}}}}\right)\right]\kern5em ifW=\left\{{z}_1,{z}_2,{z}_3\right\}\hfill \\ {}\hfill {D}_{k\Big|k-1}\left[{p}_{D,k}^4\left({\displaystyle \sum_{i=1}{\displaystyle \sum_{j=2}{\displaystyle \sum_{m=3}{\displaystyle \sum_{n=4}{\displaystyle \sum_{a=1}^4{\displaystyle \sum_{\begin{array}{l}b=1\\ {}a\ne b\end{array}}^4{\displaystyle \sum_{\begin{array}{l}c=1\\ {}a\ne c\\ {}b\ne c\end{array}}^4{\displaystyle \sum_{\begin{array}{l}d=1\\ {}a\ne d\\ {}b\ne d\\ {}c\ne d\end{array}}^4{\ell}_{z_a,i}{\ell}_{z_b,j}{\ell}_{z_c,m}{\ell}_{z_d,n}}}}}}}}}\right)\right]\kern1.6em ifW=\left\{{z}_1,{z}_2,{z}_3,{z}_4\right\}\hfill \end{array}\right. $$
(75)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Qin, Y., Ma, H., Chen, J. et al. Gaussian mixture probability hypothesis density filter for multipath multitarget tracking in over-the-horizon radar. EURASIP J. Adv. Signal Process. 2015, 108 (2015). https://doi.org/10.1186/s13634-015-0294-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13634-015-0294-y

Keywords