Next Article in Journal
An Advanced Soil Classification Method Employing the Random Forest Technique in Machine Learning
Previous Article in Journal
Multi-Objective Optimization Study of Annular Fluid Flow Structure in Cordless Core Drilling Tools
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tracking Extended Targets: Novel TPMB Filter Driven by Model and Data Collaboration

School of Electronics and Information Engineering, Liaoning University of Technology, Jinzhou 121001, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(16), 7201; https://doi.org/10.3390/app14167201
Submission received: 2 July 2024 / Revised: 10 August 2024 / Accepted: 14 August 2024 / Published: 16 August 2024

Abstract

:
In most filtering algorithms involving measurement data association, handling the complex computations due to multiple hypotheses is necessary. This paper introduces a novel Trajectory Poisson Multi-Bernoulli (TPMB) filter for tracking extended targets, facilitated by a synergy between the model and the data. This filter can track extended targets under unknown process and measurement noise. Initially, on the model-driven side, we compute multi-model transition probabilities using the posterior probabilities from models at two consecutive time points with the targets in high maneuverability state. The accuracy of the tracking algorithm is improved by calculating the improved Interacting Multiple Model (IMM) transition probability at each time step. For the data-driven aspect, the Gate-control Belief Propagation (GBP) is set in the message- passing algorithm to reduce the running time of false hypothesis associations. Thus, it is unnecessary to consider all message information when computing the likelihood matrix for target-measurement associations. Subsequently, the posterior density function of the Adaptive Square Root Cubature Kalman Filter (ASCKF) is constructed to adaptively estimate unknown process and measurement noises, while importance sampling in the current particle filter further mitigates particle degradation. Experiments demonstrate that our algorithm reduces the running time of data associations, alleviates particle degradation, and more accurately tracks maneuvering targets under nonlinear conditions and estimates their states.

1. Introduction

In the evolution of Multi-Target Tracking (MTT) technology, foundational methods centered with data association, including Probabilistic Data Association (PDA), Joint Probabilistic Data Association (JPDA), and Multiple Hypothesis Tracking (MHT), have been crucial [1]. These methods decompose the multi-target tracking problem into simpler single-target tracking problems through data association. The Bayesian multi-target filter, which propagates the multi-target posterior density over time, offers an optimal Bayesian solution for multi-target tracking challenges. The theory of Random Finite Set (RFS) has been widely applied in fields such as radar tracking [2], air traffic management [3], and autonomous driving [4], constructing optimal Bayesian filters. Representative algorithms include the Probability Hypothesis Density (PHD) filter, the Cardinalized PHD (CPHD) filter, and the Cardinalized Balanced Multi-Target Multi-Bernoulli (CBMeMBer) filter. These foundational theoretical filters, such as the PHD and CPHD filters, directly propagate the posterior first moment [5]. However, they encounter challenges in tracking trajectories.
For modern high-resolution millimeter-wave radars and lidars, the traditional point target assumption is no longer applicable, as multiple measurements can be generated at each time step. Therefore, Extended Object Tracking (EOT) has become a mainstream challenge in applications [6]. For elliptical extended targets, commonly used measurement models include the random matrix model [7], set cluster process [8], physics-based model [9], and Poisson space models [10]. Among these, the Poisson space model is the most widely used. In this model, target detection is modeled by a non-homogeneous Poisson Point Process (PPP). The PPP measurement likelihood has a simple factorization that avoids the explicit representation of assumptions about the association between measurements and targets [11]. At each time step, a Poisson distribution of random measurements is generated, spatially distributed around the measurements of the target. To address the issue of tracking multiple targets, the RFS framework is used for modeling. The RFS framework offers a systematic approach to Multi-Target Tracking (MTT) by accommodating changes in the number and state of targets over time. Both targets and measurements are modeled as random finite sets, and the PPP extended target model has been utilized in several computable complex models. The ensemble approach within this framework allows for an optimal Bayesian representation of the birth and death of targets. According to the literature, the state-of-the-art algorithm for tracking an unknown number of elliptical targets is the Poisson Multi-Bernoulli Mixture (PMBM) filter. This method acknowledges that the number, size, shape, and kinematic states of targets are interrelated and influence each other. However, due to the numerous prior conditions and associations that need to be considered, such associations present computational challenges.
According to [12,13], the most advanced algorithm for tracking an unknown number of extended targets is the PMBM filter. For estimating target trajectories, MTT methods based on random vectors link the current state estimates of an object with its past states, or represent the emergence of new targets. In the RFS-based MTT methods, one effective approach is to assign a unique label to each trajectory, allowing each object to be identified as it changes over time [14]. Another effective method is to calculate the multi-target posterior probability over the set of trajectories, which contains all information about multiple trajectories and does not require the use of labels [15]. In multi-target tracking scenarios, the Trajectory Probability Poisson Multi-Bernoulli Mixture (TPMBM) filter and its approximate filter [16], the Trajectory Probability Poisson Multi-Bernoulli (TPMB) filter, have emerged [17].
For methods based on gating and sampling, the assumption that negligible weights can be truncated avoids the exhaustive data association often required for large datasets [18]. Therefore, in scenarios with high uncertainty in data association, performance degradation can occur due to information loss. To prevent the loss of information associated with explicit enumeration and local association assumptions, the method of computing the multi-target marginal posterior for data association is used, which marginalizes the unknown data association uncertainties [10]. Currently, the best-performing multi-EOT algorithm, as proposed in [19], is the Sum-Product Algorithm (SPA) based on particle Belief Propagation (BP), which is implemented through a factor graph of the joint posterior of multi-target states and measurement association variables. Simulation results in [20] demonstrate that the TPMB filter based on factor graphs outperforms the TPMBM filter that uses gating clustering and truncated sampling in terms of estimation errors and runtime. The BP algorithm recursively iterates to obtain global hypothesis data association information, transmitting this association information between factor nodes and variable nodes via a directed graph. In the process of approximating PMB with PMBM, Kullback–Leibler Divergence (KLD) minimization is used to marginalize the uncertainty of data hypotheses, thereby reducing the amount of data processing. This process has led to the development of a closed-form filter that initiates the recursive process from the joint posterior factor graph model of object sets and association variables, along with the message-passing equations.
This paper presents a model-data-driven message-passing algorithm that adaptively adjusts the motion models of targets in different dynamics through an Interactive Multiple Model (IMM) to adapt to the switching motion state at different time steps [21]. In the message-passing process using factor graphs, considering the measurement association factor’s message as an approximation of the posterior probability requires computation for each message-passing variable, which is complex. The proposed algorithm reduces the complexity of data association and minimizes message passing by gating hypothetical surviving targets, thereby enhancing algorithm efficiency. In the filter, a new particle filtering method is introduced that constructs the importance density function for each particle’s state via the Adaptive Square Root Cubature Kalman Filter (ASCKF), optimizing the covariance matrices of process noise and measurement noise to improve the tracking performance of high-maneuver targets under noise and clutter influences [22,23]. This paper is an important extension of [19,22], including the following contributions:
  • In addressing the complexity of model transition probability calculations in the IMM algorithm, we utilize changes in the posterior probability between two consecutive moments to reflect the transition probabilities between models. On the other hand, the standard IMM algorithm considers process noise as a constant during model switching, affecting the accuracy of the filtering algorithm. The improved IMM algorithm adaptively adjusts the magnitude of process noise, reducing errors in model switching.
  • In addressing the complex message passing between variable nodes and factor nodes in factor graphs, we have improved the posterior association step of the factor graph. The message-passing process is streamlined through adaptive gating, eliminating the need to consider all inter-node relationships when computing association nodes, thereby reducing running time.
  • To address the issue of particle degradation observed in the methods using particles, we employ the SCKF to perform particle importance sampling and use the Gate-control Belief Propagation (GBP) for the improved IMM-TPMB-GBP filter, updating the particle filter based on these sampling results. Due to the decrease in tracking accuracy of adaptively modeled maneuvering targets under noise and clutter, we have developed the ASCKF. This involves adaptively estimating the covariance of process and measurement noise, ensuring that the algorithm’s performance does not depend on noise statistics and is more robust when these statistics are unknown.
The rest of the paper is organized as follows. Section 2 introduces the background on state and trajectory parameters, the trajectory set, and TPMB. Section 3 provides the BP cycle for TPMB marginal posterior probabilities and the principles of factor graphs, further detailing the adaptive gating improvement for factor graphs and the ASCKF mixed particle filtering method. The implementation of ASCKF for unknown noise is presented in Section 4. Simulation results are provided in Section 5, and the conclusion is given in Section 6.

2. Background

In this section, we first introduce the state and trajectory parameters of targets, then model the multi-target extended model state parameters and measurements, and finally present the TPMB approximation of the PMBM posterior density.
A finite spatial set is D, which is denoted as F ( D ) , and the cardinality of its subset A F ( D ) is | A | . We use ⨄ to denote the union of disjoint sets, f , g to denote the inner product f ( x ) g ( x ) d x , and the multi-target index is represented as f A . For certain real-valued functions, the product x A f ( x )   with   f = 1 is conventionally represented. Additionally, δ x [ · ] and δ x ( · ) , respectively, denote the Kronecker delta centered at x , and the Kronecker delta function. Here, x in δ x [ · ] represents a subspace V within the Cartesian product space, x X V .

2.1. State and Trajectory Variables

At time step k , the single-target state x k X d x includes its kinematic state and any potentially identified target extension states. At time step k , the set of target states x k is a random finite set, and x k = { x k 1 , , x k n k } F ( X d x ) . The measurements z k = { z k 1 , , z k m k } collected by sensors at time step k consist of the jth single-target measurement z k j X d z produced at time step k and j { 1 , , m k } . The sequence from the initial moment to the current time k is denoted as z k .
A target’s motion trajectory is its finite sequence of states over continuous time steps, represented as X = ( t , x 1 : ν ) . Here, t is the initial time step of the trajectory, v is the duration of the trajectory steps, and x 1 : ν = ( x 1 , , x ν ) is the finite sequence representing the target state trajectory. The parameter variable t , v belongs to the set I ( k ) = ( t , ν ) : 0 t k , 1 v k t + 1 . Thus, the single-target trajectory reaching a finite time step at moment k belongs to space T ( k ) = ( t , ν ) I ( k ) { t } × X v d x , with ⨄ representing the union of disjoint sets. The trajectory set up to the current moment k can be represented as X k F ( T ( k ) ) . The single trajectory function can be decomposed into p ( X ) = p ( x 1 : ν | t , ν ) p ( t , ν ) , where p ( t , ν ) is defined on I ( k ) and X T ( k ) .

2.2. Multi-Target Models for Sets of Trajectories

  • Multi-target dynamic models: At time step k , existing targets x k 1 survive to the next time step with probability p S ( x k 1 ) , or perish with probability 1 p S ( x k 1 ) . If a target survives, its Markov state transition function is g k ( · | x k 1 ) . At each time step, newly detected targets x k are independent of other surviving targets and follow a PPP with a birth intensity of λ k B ( x k ) .
  • Multi-target measurement models: The measurement set z k for targets is divided into two parts: one consists of actual measurement sequences generated by extended single-target detection by sensors, modeled as PPPs with measurement likelihoods of l k ( z k | x k ) produced by each x k X d x and with a Poisson intensity of γ k ( x k ) . The other consists of clutter measurements induced by measurement noise, modeled as a PPP with a Poisson intensity of λ k C ( z k ) = γ k C μ k C ( z k ) , a Poisson rate of γ k C and clutter intensity of μ k C ( z k ) .
  • Dynamic model for sets of all trajectories: For the set of surviving trajectories X k 1 at time step k 1 , each trajectory X = ( t , x 1 : ν ) X k 1 has a survival probability p S ( X ) = p S ( x v ) . If the trajectory continues to survive in the next time step, its transition probability density is
    g k ( t ¯ , y 1 : ν ¯ | X ) = δ t [ t ¯ ] δ ν + 1 [ ν ¯ ] δ x 1 : ν ( y 1 : ν ¯ 1 ) g k ( y ν ¯ | x ν )
    Newly initiated trajectory sets X k independently follow a PPP with an intensity of
    λ k B t , x 1 : ν = δ k [ t ] δ 1 [ ν ] λ k B x ν
    Given the complete set X k of trajectories, both surviving and newly initiated, at time step k for each X = t , x 1 : v X k , t + v 1 k , the survival probability is 100%, i.e., p S ( X ) = 1 , and the trajectory’s transition density is
    g k t ¯ , y 1 : ν ¯ | X = δ t [ t ¯ ] × δ ν [ ν ¯ ] δ x 1 : ν y 1 : ν ¯ ω ¯ < k 1 δ ν [ ν ¯ ] δ x 1 : ν y 1 : ν ¯ 1 p S x ν ω ¯ = k 1 δ ν + 1 [ ν ¯ ] δ x 1 : ν y 1 : ν ¯ 1 p S x ν g k y ν ¯ | x ν ω ¯ = k
    Thus, the current moment’s survival probability of a trajectory depends solely on the final state of the trajectory in the previous moment and is independent of other moments. The probability density of the measurement sequence w k produced by a single trajectory ( t , x 1 : ν ) X k is l k ( w k | t , x 1 : ν ) = l k ( w k | x ν ) .

2.3. Standard Bayesian Model for Sets of Trajectories

By extending an auxiliary variable u U k | k = { 0 , 1 , , n k | k } over the single-target trajectory space, ( u , X ) U k | k × T ( k ) is used. When u = 0 represents that the target trajectory was not detected, it corresponds to a PPP; when u = { 1 , , n k | k } , it indicates that the target trajectory corresponds to the ith Bernoulli component. After adding the auxiliary variable, we use X ˜ F ( U k | k × T ( k ) ) to represent the set of target trajectories with the auxiliary variable. w k | k a denotes the weight, ⨄ denotes the union of disjoint sets; f k | k p ( Y ) is the Poisson part probability density function, representing the trajectories of objects assumed to exist but undetected; λ k | k ( · ) is the Poisson intensity for trajectories; f k | k m b ( V ) is the multi-Bernoulli part probability density function, representing the probability of latent trajectories being detected at least once by some time k ; in the multi-Bernoulli probability density Equation (5), n k | k indicates the number of Bernoulli components that exist up to some time k , with each Bernoulli component having h k | k i independent hypotheses; a represents the global hypothesis index. For the ith Bernoulli component density, there are a i local hypotheses and a i 1 , , h k | k i . For a global hypothesis a = ( a 1 , , a n k | k ) , a local hypothesis a i is assigned to each Bernoulli component, and A k | k is the set of global hypothesis indices; that is, a A k | k .
w k | k a = i = 1 n k | k w k | k i , a i b A k | k i = 1 n k | k w k | k i , b i
where w k | k i , a i represents the weight of the local hypothesis Bernoulli component, and w k | k a represents the weight of the global hypothesis Bernoulli component. Normalization ensures that the sum of the probabilities of the global Bernoulli components equals 1; that is, a A k | k w k | k a = 1 .

2.4. PMBM Conjugate Prior

TPMB is a special filter that approximates the TPMBM posterior probability density. In this approximation, trajectory-oriented MBMs are merged into a single MB by introducing the parameter variable, thus handling the mixture of densities. The preferred method for approximating PMBM is minimizing the Kullback–Leibler divergence (KLD).
Given the multi-target probability density defined in Equation (5), after the inclusion of parameter variables, the multi-target probability density f ˜ k | k ( · ) can be represented as follows:
f ˜ k | k X ˜ k = l = 1 n k | k X ˜ l Y ˜ = X ˜ k f ˜ k | k p Y ˜ a A k | k w k | k a i = 1 n k | k f ˜ k | k i , a i X ˜ i = f ˜ k | k p Y ˜ k a A k | k w k | k a i = 1 n k | k f ˜ k | k i , a i X ˜ i
For the given surviving targets X ˜ k i = { ( u , X ) X ˜ k : u = i } , the undetected targets Y ˜ k = { ( u , X ) X ˜ k : u = 0 } are included, and their predicted multi-target probability density after incorporating the parameter variables is the following:
f ˜ k | k p X ˜ = e λ k | k ( x ) d x λ ˜ k | k · X ˜
λ ˜ k | k ( u , X ) = δ 0 u λ k | k ( X )
f k | k i , a i X ˜ = 1 r k | k i , a i X = r k | k i , a i p k | k i , a i ( X ) δ i u X = { ( u , X ) } 0 otherwise
where the auxiliary variable u = 0 may represent multiple undetected trajectories, and it does not necessarily imply that f k | k · equals zero.

2.5. Trajectory PMB Approximation

According to reference [16], the TPMBM approximates the TPMB method, and this paper also utilizes special multi-trajectory states and measurement models. According to Equation (10), we derive the TPMB prediction equation.
Given the multi-trajectory set prediction density f ˜ k | k X ˜ k and measurement set z k = { z k 1 , , z k m k } at time step k , from Equation (5) to Equation (8), it is known that the number of updated Bernoulli components n k | k = n k | k 1 + m k , and the updated Poisson intensity is as follows:
λ k | k ( X ) = l k ( X ) λ k | k 1 ( X )
For each Bernoulli component predicted at the previous measurements i 1 , , n k | k , the current moment’s new measurements generate m k Bernoulli components. At the same time, each Bernoulli component predicted from the previous moment generates m k new local hypotheses at the current moment. These newly generated local hypotheses may be false detection or updates for non-empty subsets of measurements z k . After augmenting the model using auxiliary variables, it is necessary to define the global hypothesis set of generalized measurement model parameter variables. We give the augmented representation method of measurement M ( k , j ) = j for z k j . For the Bernoulli hypothesis component i, a set of measurement collections represented by a local hypothesis a i is M k i , a i M k :
M k = M k 1 { ( k , j ) j { 1 , , m k } }
A k = a 1 , , a n k | k : a i h i ,   i = 1 n k | k M i , a i = m k , M i , a i M j , a j = ,   i j
where m k = { 1 , , m k } .
For each Bernoulli component i { 1 , , n k | k 1 } in the predicted probability density f k | k 1 ( · ) and the new measurements detected at the current moment generating m k Bernoulli components, there are h i = m k + 1 local hypotheses corresponding to either a false detection or an update at the current moment. For the false detection hypothesis of the Bernoulli component i 1 , , n k | k 1 , it can be represented as M ( i , 1 ) = .
l k | k i , 0 = f k | k 1 i , l k ( w k j · )
w k | k i , 1 = w k | k 1 i , 1 1 r k | k 1 i , 1 + r k | k 1 i , 1 l k | k i , 1
r k | k i , 1 = r k | k 1 i f k | k 1 i , l k ( w k j · ) w k | k i , 1
f k | k i , 1 ( X ) = l k ( w k j X ) f k | k 1 i ( X ) f k | k 1 i , l k ( w k j · )
Consider that the measurement set is empty in Equations (12)–(15). We use w k 1 , , w k 2 m k 1 to represent the non-null measurement set z k in Equations (16)–(18). For the predicted Bernoulli component i 1 , , n k | k 1 at time step k 1 , it survives to the next moment and is associated with measurement z k j , which yields, j 1 , , 2 m k 1 , M ( i , j ) = 1 , r k | k i , j + 1 = 1 .
l k | k i , j + 1 = f k | k 1 i , l k w k j ·
w k | k i , j + 1 = w k | k 1 i r k | k 1 i l k | k i , j + 1
f k | k i , j + 1 ( X ) = l k v k j X f k | k 1 i ( X ) l k | k i , j + 1
Then, we construct the set S i to associate with the new Bernoulli component i 1 , , m k , with S i = z k i w j = 1 i 1 S i z k i w , and use w k i , ι to represent the Bernoulli component of the measurement set ι , where i = n k | k 1 + j , j { 1 , , m k } , ι { 1 , , 2 j 1 } . For the new Bernoulli components i = n k | k 1 + j generated by new measurements z k j at time step k , there are two single-target hypotheses for each Bernoulli component: one corresponds to a previously undetected and non-existent Bernoulli component, represented as M ( i , 1 ) = , r k | k i , 1 = 0 , w k | k i , 1 = 1 ; the other corresponds to a Bernoulli component that existed in the previous moment.
l k | k i , 2 = λ k | k 1 , l k w k i , ι ·
w k | k i , 2 = λ C w k j + λ k | k 1 , l w k i , ι ·
r k | k i , 2 = λ k | k 1 , l w k i , ι · λ C w k i , ι + λ k | k 1 , l w k i , ι ·
f k | k i , a i ( X ) = l k w k i , ι X λ k | k 1 ( X ) l k | k i
The set M k i , j represents the local hypotheses a i related to the Bernoulli component i measurement indicators. After the Bayesian update, the updated density of PMBM is obtained through the aforementioned process. The established generalized measurement of the model’s trajectory posterior set, utilizing the posterior conjugacy of PMBM, provides convenience for using general nonlinear target measurement models.

3. Model-Driven IMM-TPMB Filter

The TPMBM obtains a computationally less intensive TPMB posterior density through KLD approximation. However, during the approximation process, the uncertainty in the marginal density of the PMB trajectory set increases, and the uncertainty of the global hypotheses has been marginalized. For methods based on gating and sampling, negligible hypotheses can be ignored by truncating weights to avoid exhaustive data association. Nevertheless, in scenarios with high data association uncertainty, performance can degrade due to information loss. Therefore, to prevent the loss of information from local association hypotheses caused by explicit enumeration, a data association method that computes the multi-target marginal posterior is used to marginalize the unknown data association uncertainty. The model established in this paper is based on the factor graph probability model in [17,18], which enables SPA to have good computational complexity and utilizes the accuracy and flexibility of an SCKF-based implementation.
We define the measurement-vector association factor for multi-target trajectories. Due to the uncertainty in the posterior marginal association probability for each measurement j 1 , , m k , the indicator factor for the Bernoulli component generated by each measurement is defined as follows:
β k = β k j β k j 1 , , n k | k 1 0 otherswise
where β k j = i > 0 only when the measurement z k j is associated with the ith potential target; when β k j = 0 , it indicates that the jth measurement is a false alarm.
The joint posterior probability of the trajectory set and the measurement vector association factor can be factored into a factor graph. According to [10], the joint posterior of the trajectory set at time step k and the measurement vector association model during the TPMB model update process is
f ˜ k | k X ˜ k , β k f ˜ k | k p Y k I i = 1 n k | k 1 f _ k | k 1 i X ˜ k i j = 1 m k p _ s , k X ˜ k i , β k j ; z k j II i = n k | k 1 + 1 n k | k f ¯ k | k 1 i X ˜ k i p ¯ b , k X ˜ k i , β k i n k | k 1 ; z k i n k | k 1 j = 1 i n k | k 1 1 p _ s , k X ˜ k i , β k j ; z k j ] III
The joint posterior expression consists of three parts: the first part represents the updated posterior probability density of the undetected trajectory set; the second part represents the joint posterior of the survival probability density of the Bernoulli component i 1 , , n k | k 1 from the previous moment and the measurement association vector at this moment; the third part represents the joint posterior of the probability density function of the newly detected trajectories and the measurement association vector after receiving m k measurements.
According to Equation (24), the factor graph is depicted with circles representing factor nodes and squares representing variable nodes. For undetected trajectories, only an additional factor node and a variable node are used to represent them, independent of the factor graph’s message-passing process. For the surviving trajectories at time step k 1 , it is predicted that the next moment will have n k | k 1 Bernoulli components, represented here by N factor nodes. Among these N factor nodes, an equal number of m k variable nodes are generated, which connect to the measurement vector association nodes as the posterior marginal probability density of this Bernoulli component. For the m k Bernoulli components generated by the measurements at time step k , it is indicated that M factor nodes are generated. In these M factor nodes, each factor node produces a different number of posterior marginal probability density nodes.
In the factor graph message-passing model, as the starting variable nodes in the message-passing process, f ¯ k | k 1 i ( x ) and f _ k | k 1 i ( x ) correspond to the prior probability densities of undetected and surviving targets, respectively, while f ˜ k | k p ( x k u ) represents the posterior probability density of undetected targets and newly detected targets.
f _ k | k 1 i X ˜ = r k | k 1 i f k | k 1 i ( X ) e γ k ( X ) δ i [ u ] X ˜ = { ( u , X ) } 1 r k | k 1 i X ˜ = 0 otherwise
f ¯ k | k 1 i X ˜ = λ k | k 1 ( X ) e γ k ( X ) δ i [ u ] X ˜ = { ( u , X ) } 1 X ˜ = 0 otherwise
In Figure 1, the factor nodes p _ S , k X ˜ k i , β k j ; z k j and p ¯ S , k X ˜ k i , β k j ; z k j in the information propagation are represented as follows:
p _ s , k X ˜ k i , β k j ; z k j = l k ( z k j x ) γ k ( X ) / λ k c ( z k j ) X ˜ k i = { ( u , X ) } , β k j = i 1 β k j i 0 otherwise
p ¯ s , k X ˜ k i , β k j ; z k j = l k ( z k j x ) γ k ( X ) / λ k c ( z k j ) X ˜ k i = { ( u , X ) } , β k j = i 1 X k i = , β k j i 0 otherwise

3.1. Implementation of an IMM-TPMB Filter Based on Two Consecutive Time Points

In multi-maneuvering target tracking scenarios, the motion states at different times may be described by multiple motion models, and the filtering results must consider the state outputs from various models. The IMM method is an effective solution for such problems. The IMM method uses multiple different models to model the target’s state, and each model could be the optimal one at the current moment.
Assuming the number of tracking models is ψ , the motion matching model at time step k 1 is M k 1 ξ , and at time step k it is M k ξ , with a model mixing probability of μ k 1 ξ | ξ , it can be represented as follows:
μ k 1 ξ | ξ = P M k 1 ξ M k ξ , Z k 1 = 1 / c ¯ ξ P M k 1 ξ M k ξ , Z k 1 P M k 1 ξ Z k 1 = 1 / c ¯ ξ π k 1 ξ | ξ μ k 1 ξ
where the Markov transition probability from ξ to ξ at time step k 1 for model π k 1 ξ | ξ , μ k 1 ξ represents the probability of model ξ being the matching model at time step k 1 , and c ¯ ξ represents the normalization constant, which can be expressed as follows:
c ¯ ξ = ξ = 1 ψ π k 1 ξ | ξ μ k 1 ξ , ξ = 1 , , ψ
Using the initial state mean x ^ k 1 i of model i , the initial conditions for model ξ are calculated.
x ^ k 1 0 ξ = E X k 1 M k ξ , Z k 1 = ξ = 1 N x ^ k 1 ξ μ k 1 ξ | ξ
P k 1 0 ξ = ξ = 1 N μ k 1 ξ | ξ P k 1 ξ + x ^ k 1 ξ x ^ k 1 0 | ξ x ^ k 1 ξ x ^ k 1 0 | ξ T
where x ^ k 1 0 ξ and P k 1 0 ξ represent the mixed state matrix and mixed covariance input to the filter, respectively. The filtering algorithm is used for prediction and updating. After each time update, it is necessary to recompute the Markov matrix. However, the computation of the transition matrix is complex and susceptible to noise. To address this issue, this section adopts the posterior probability density of each model to calculate the transition probabilities π k 1 ξ | ξ . As the transition probabilities are closely related to state changes, and the state changes are mainly determined by the posterior probability p k ξ of model M k ξ at time step k , the transition probabilities between models can be reflected using changes in the posterior probabilities between two consecutive moments.
Given the posterior probability density μ k 1 ξ of model M k 1 ξ at time step k 1 , the relative rate of change of probability between two adjacent moments Δ k ξ | ξ can be given by the following proportional form:
Δ k ξ | ξ = μ k ξ μ k 1 ξ
As indicated by Equation (33), when the probability density μ k ξ of model M k ξ at time step k increases, Δ k ξ | ξ > 1 . This means that the posterior probability density from other models at k 1 moments transitioning to the current model gradually increases; the converse is also true. Since the posterior probability density μ k ξ changes over time, π k ξ | ξ can be represented by the model transition probability from the previous moment π k 1 ξ | ξ and the relative rate of change of probability Δ k ξ | ξ . Assuming a normalization constant ρ k ξ | ξ = Δ k ξ | ξ ξ = 1 M Δ k ξ | ξ π k 1 ξ | ξ , the update of the transition probability π k ξ | ξ can be represented as follows:
π k ξ | ξ = Δ k ξ | ξ π k 1 ξ | ξ ξ = 1 M Δ k ξ | ξ π k 1 ξ | ξ = ρ k ξ | ξ π k 1 ξ | ξ
According to Equation (34), π k ξ | ξ is represented as the product of ρ k ξ | ξ and π k 1 ξ | ξ , and there exists ξ = 1 M Δ k ξ | ξ = 1 , thus satisfying the Markov transition matrix conditions. Obviously, the above parameters can inherently match the maneuvering dynamics, including the CT motion and the motion switching. As a result, the maneuvering targets with different dynamics can be detected effectively.
On the other hand, the intensity of target maneuvers directly determines the magnitude of process noise. The standard IMM algorithm considers process noise as a constant, which affects the accuracy of the filtering algorithm. Moreover, the process noise changes with unknown noise influences. Therefore, it is necessary to adapt the IMM algorithm to adjust the process noise to accommodate scenarios where the target is maneuvering or non-maneuvering and the noise is unknown. When the target is maneuvering, increase the process noise V k to compensate for modeling errors, resulting in smaller filtering estimation errors; when the target is non-maneuvering, reduce V k to maintain tracking accuracy. Assuming μ max , k 1 ξ is the maximum value of the posterior probability density of each model at time step k 1 , then the covariance matrix Q k ξ of the process noise V k at time step k can be represented as follows:
Q k ξ = 1 μ max , k 1 ξ 2 Q 0 ξ , ξ = ξ 1 + μ max , k 1 ξ 2 Q 0 ξ , ξ ξ
where Q 0 ξ is the initial value of the process noise covariance matrix for model M k 1 ξ . According to Equation (35), when the models described by M k 1 ξ and M k ξ are the same, M k 1 ξ matches the current target motion model, and Q k ξ takes the minimum value; when M k 1 ξ and M k ξ describe different models, M k 1 ξ does not match the current target motion model, and Q k ξ takes the maximum value. Considering the impact from the process noise, the adaptive adjustment of Q k ξ ensures consistency with the varying intensity of target maneuvers.

3.2. Implementation of an IMM-TPMB Filter Based on Gate-Control

3.2.1. Predict

  • Undetected and newly detected targets: at time step k , the PPP prediction intensity for all undetected and newly detected targets within the sensor-monitored area is as follows:
    λ k | k 1 ( t , x 1 : v ) = l = 1 L k | k 1 u w k | k 1 0 , ( l ) δ t k | k 1 0 , l [ t ] δ k t k | k 1 0 , ( l ) + 1 [ v ] δ x k | k 1 0 , ( l ) [ x 1 : v ] + l = 1 L k | k 1 b w k b , ( l ) δ k [ t ] δ 1 [ v ] δ X k b , ( l ) [ x k ]
    X ˜ k | k 1 0 , ξ , ( l ) = F ( · X ˜ k 1 | k 1 0 , ξ , ( l ) )
    L k | k 1 u = L k | k 1 u + L b
    w k | k 1 0 , ξ , ( l ) = P S ( X ˜ k 1 | k 1 0 , ξ , ( l ) ) w k | k 1 0 , ξ , ( l )
    w k | k 1 b , ( l ) = λ k B , 1 L b
    where the probability density of the single-target particle belief during propagation is f ( x ) = l = 1 L w ( l ) δ x ( l ) [ x ] , the particle weight indexed by l is w ( l ) 0 , and δ i u = 1 only when u = i ; otherwise, it is 0. L k b indicates the number of newly detected target particles, and the target component is in the form of particles w 0 , ξ , ( l ) , x 0 , ξ , ( l ) l = 1 L .
  • Surviving targets: for the detected surviving targets, the target state of the ith MBM component at time step k is represented as follows:
    f k | k 1 i ( t , x 1 : v ) = l = 1 L k | k 1 d w k | k 1 i , ( l ) δ t k | k 1 0 , ξ , ( l ) [ t ] δ v k | k 1 0 , ξ , ( l ) [ v ] δ x k | k 1 0 , ξ , ( l ) [ x 1 : v ] + l = 1 L k | k 1 d 1 P S ( X ¯ k 1 | k 1 i , ( l ) ) w k 1 | k 1 i , ( l ) δ t k 1 | k 1 0 , ξ , ( l ) [ t ] × δ k t k 1 | k 1 0 , ξ , ( l ) [ v k 1 | k 1 i , ( l ) ] δ X k 1 | k 1 0 , ξ , ( l ) [ x 1 : v ]
    X ¯ k | k 1 ξ , ( l ) = m = 1 M μ k ξ | ξ F k ξ | ξ X ¯ k 1 | k 1 0 , ξ , ( l )
    r k | k 1 i = r k 1 | k 1 i l = 1 L k | k 1 d P S X ¯ k 1 | k 1 i , ξ , ( l ) w k 1 | k 1 i , ξ , ( l )
    L k | k 1 d = L k 1 | k 1 d
    w k | k 1 i , ξ , ( l ) = P S ( X ¯ k 1 | k 1 i , ξ , ( l ) ) w k 1 | k 1 i , ξ , ( l ) l = 1 L k | k 1 d P S ( X ¯ k 1 | k 1 i , ξ , ( l ) ) w k 1 | k 1 i , ξ , ( l )
    In the case of undetected and newly detected targets, there is no need to consider the target model transition probability. F · represents the motion model transition function, which performs maneuvering target particle sampling through multi-model mixing prediction.

3.2.2. Update

In the BP iterative propagation algorithm, the optimization of the propagated information is achieved through P iterations, leading to the convergence of belief for all nodes. At this point, the label of each node aligns with the label of the node with the highest posterior probability density among the related nodes, representing the optimal label. In the update step, we improve the way of information transmission and reduce the amount of data association calculation by using the GBP. The specific improvement method is in step 3. The update step in model-driven scenarios includes six stages: particle initialization, importance sampling of particles, iterative information transfer, belief calculation, pruning and resampling, and target state estimation.
  • Particle initialization: The Poisson intensity is represented as a set of particles denoted by t k | k 1 0 , ξ , ( l ) , X k | k 1 0 , ξ , ( l ) w k | k 1 0 , ξ , ( l ) e γ k X ¯ k | k 1 0 , ξ , ( l ) l = 1 L k | k 1 u . Updating for new, missed, and surviving targets requires considering whether targets are alive at the current moment.
  • Particle importance sampling: Utilize the ASCKF for importance sampling of particles. After iterating the prediction and updating process of prior and posterior information at time step k 1 , the particles at time step k follow the distribution function
    x k i , ξ , ( l ) q ( X ¯ k i , ξ , ( l ) , Z k ) = N ( x ^ k i , ξ , ( l ) , P k i , ξ , ( l ) )
    Here, N ( x ^ k i , ξ , ( l ) , P k i , ξ , ( l ) ) denotes a Gaussian distribution with mean x ^ k i , ξ , ( l ) and variance P k i , ξ , ( l ) , as determined by the ASCKF.
  • Iterative information transfer: Measurement assessment requires calculating the information m k j transferred from each factor node p _ k to variable node β k j in the vector measurement. Factor node p _ k represents the process of information transfer between surviving trajectories and both sides.
    m i , j ( p , l ) ( i ) = 1 λ k C ( w k j ) u i , j ( p ) l = 1 L i w i , j ( p , l ) δ k t k | k i , ξ , ( l ) [ v k | k i , ξ , ( l ) ] × γ k ( X ¯ k | k i , ξ , ( l ) ) l k ( w k ξ X ¯ k | k i , ξ , ( l ) )
    γ j , i ( p , l ) ( X ¯ k i ) = m k j ˜ p ¯ k X ¯ k i , ξ , ( l ) , β i , j p · ; w k j v j , i ( p ) β i , j p · = l k w k j | X ¯ k | k ξ , ( l ) γ k X ¯ k | k ξ , ( l ) λ k c ( z k j ) + i = 1 n k | k m i , j ( p ) + 1 X k ξ , ( l ) = { x } i = 1 n k | k m i , j ( p ) ( i ) + 1 X k ξ , ( l ) =
In p { 1 , , P } iterations, based on the factor graph formula, the running time of calculating particle weights during the information propagation between each variable node and factor node increases. As the number of targets increases, we also consider the propagation time between factor nodes and variable nodes. The proposed GBP is set to reduce the running time of false hypothesis associations. Therefore, the main computational cost is effectively saved by only achieving the related weights. As a result, the tracking accuracy is well ensured.
There is no need to consider the information propagated across all nodes; instead, only the information within the gating needs to be calculated, thus reducing the weight calculation. The likelihood function for target generation is represented as follows:
l k Z k x k | k i , ξ , ( l ) = e γ k ( x k | k i , ξ , ( l ) ) z k I w k γ k x k | k i , ξ , ( l ) l k w k j x k | k i , ξ , ( l ) d ( x , z ) δ 0 o t h e r w i s e
where d ( x , z ) represents the Euclidean distance between the surviving target particles and the current measurement location, with δ being the gating threshold for particle-measurement association. After one iteration, particle weights converge towards values close to the true state of the target. In subsequent iterations, the information u i , j ( p + 1 , l ) propagated from the variable node x ˜ k l to the factor node p _ k represents the unnormalized Bernoulli component density.
When i { 1 , , n k | k 1 } , j { 1 , , m k } , predict the next moment’s surviving target generation of particles and update with current measurements:
u i , j ( p + 1 ) = r k | k 1 i n = 1 M l = 1 L i w i , j ( p + 1 , n , l ) + 1 r k | k 1 i
w i , j ( p + 1 , l ) = n = 1 M w i , j ( 1 , n , l ) j = 1 m k l k ( z k j x k | k i , ξ , ( l ) ) γ k ( x k | k i , ξ , ( l ) ) λ k C ( z k j ) i = 1 n k | k m i , j ( p ) ( i ) + 1 + 1
When i { n k | k 1 + 1 , , n k | k } , j { 1 , , i n k | k 1 } , predict the next moment’s missed and new target generation of particles and update with current measurements. The information u i , j ( p + 1 ) and particle weights are expressed as follows:
u i , j ( p + 1 ) = l = 1 L i w i , j ( p + 1 , l ) + 1
w i , j ( p + 1 , l ) = w i , j ( 1 , l ) j = 1 j j i n k | k 1 l k ( z k j | x k | k i , ξ , ( l ) ) γ k ( x k | k i , ξ , ( l ) ) λ k C ( z k j ) i = 1 n k | k m i , j ( p ) ( i ) + 1 + 1
4.
Belief calculation: After multiple iterations of the above steps, the information and example weights converge to the true posterior probability density of the target. The belief f ˜ k | k p ( · ) for the target state in each MBM component is represented as a set r k | k i , w k | k i , ξ , ( l ) , t k | k i , ξ , ( l ) , X k | k i , ξ , ( l ) l = 1 L i of particles, and its calculation method is similar to information u i , j ( p + 1 ) . Additionally, subsequent iterations still require calculating particle weights and the existence probability of Bernoulli to ensure the validity of trajectory belief, followed by weight normalization.
For new targets i in the detection area at time step k , when i { 1 , , n k | k 1 } , then the updated existence probability of the target is as follows:
r k | k i , l = r k | k i l = 1 L I w i ( P + 1 , l ) r k | k 1 i l = 1 L I w i ( P , l ) + 1 r k | k 1 i
For undetected potential targets i , when i { n k | k 1 + 1 , , n k | k } , then the updated existence probability of the target is the following:
r k | k i , l = l = 1 L w i ( P + 1 , l ) r k | k 1 i l = 1 L w i ( P , l ) + 1
The updated particle weight is normalized as follows:
w k | k ( i , ξ , l ) = w i ( P + 1 , l ) l = 1 L I w i ( P + 1 , l )
5.
Pruning and resampling: To reduce the number of global hypotheses, pruning and resampling of the MBM components are necessary. By approximating MBMs with low existence probabilities as PPPs, the algorithm increases computational speed without significantly impacting accuracy. The calculation of distortion caused by minimizing the KLD under the PPP approximation is as follows:
D K L f k | k 1 s , f | | f k | k 1 u , b = r k | k 1 + 1 r k | k 1 log 1 r k | k 1
6.
Target state estimation: Select the MB components with existence probabilities above a set threshold to estimate the current target state at time step k , from a set of weighted particle collections, represented as the outputted target trajectory X k i , ξ , ( l ) = 1 l L k i w k | k i , ξ , ( l ) l L k i w k | k i , ξ , ( l ) X ¯ k | k i , ξ , ( l ) . New target estimates are achieved using a Poisson point process, outputting a set of new target trajectories, denoted as trajectory X b = t b , x 1 : v , where t b = k is the birth time, and v = 1 is the trajectory length.

4. IMM-TPMB Filter Driven by Model and Data Collaboration

Addressing the particle degeneration issue identified in [20], this paper proposes a SCKF-based particle importance sampling method for IMM-TPMB-GBP. The improved IMM-TPMB-GBP particle filter employs the ASCKF to obtain state estimates and their covariance matrices at each moment, using these as the mean and variance in the original filter to establish a proposal distribution function. Innovatively, in scenarios where noise statistics are unknown, the ASCKF provides an adaptive noise ASCKF, enhancing the accuracy of particle weights. While particle filters typically face degeneration, the ASCKF performs moderately in systems with non-Gaussian noise estimation. Utilizing the advantages of ASCKF in estimating unknown noise and nonlinear estimation, this paper combines the improved ASCKF with IMM-TPMB-GBP, proposing a new state estimation method. Further, the ASCKF in this paper adaptively estimates the covariance of process noise and measurement noise in interacting multiple models, not relying on noise system metrics, thus offering robust performance.
We omit the auxiliary variables l of the target state and use the node i to represent the highest posterior probability density. In a nonlinear model, the target state equation and the measurement equation for a motion model ξ are formulated as follows:
X k = F X k 1 ξ X k 1 ξ + G ω k 1 Z k = x k 2 + y k 2 arctan y k / x k + v k
where G = T 2 / 2 T 0 0 0 0 T 2 / 2 T denotes the input matrix, assuming the position and velocity components are x k , y k and x ˙ k , y ˙ k , respectively, T represents the sampling interval, X k can be expressed as X k = x k , x ˙ k , y k , y ˙ k T , ω k is process noise with zero mean and the covariance matrix Q k . Z k represents the measurement set at time step k , V k is the measurement noise with zero mean and the covariance matrix R k .

4.1. Predict

Initialize the state mean x ^ 0 and the error covariance P ^ 0 and decompose the Gaussian terms’ covariance in the target trajectory composed by Equation (44) through Cholesky decomposition c h o l · .
S 0 = c h o l P 0 T
Use the third-order spherical-radial volume criteria to generate volume points χ i , k 1 , where χ i , k 1 = S k 1 ξ i + x ^ k 1 , i = 1 , 2 , , 2 n . The volume points after propagation are yielded through the nonlinear state equation:
χ i , k | k 1 = f χ k 1 | k 1 ξ
Calculate the predicted state mean and square root covariance from the volume points after propagation:
x ^ k | k 1 = 1 2 n i = 1 2 n χ i , k | k 1 S k | k 1 = T r i a X k | k 1 , Q k X k | k 1 = 1 2 n χ 1 , k | k 1 x ^ k | k 1 , χ 2 , k | k 1 x ^ k | k 1 , , χ 2 n , k | k 1 x ^ k | k 1
where T r i a · represents the triangulation operation. The volume points after propagation are yielded through the measurement equation:
χ i , k | k 1 = h χ i , k | k 1
Calculate the measurement mean and the square root covariance matrix:
z ^ k | k 1 = 1 2 n i = 1 2 n χ i , k | k 1 S z z , k | k 1 = T r i a Z k | k 1 , S R Z k | k 1 = 1 2 n χ 1 , k | k 1 z ^ k | k 1 , χ 2 , k | k 1 z ^ k | k 1 , , χ 2 n , k | k 1 z ^ k | k 1
Calculate the cross-covariance between the state and measurement:
S x z , k | k 1 = X k | k 1 Z k | k 1 T
Calculate the Kalman gain:
K k = S x z , k | k 1 S z z , k | k 1 T S z z , k | k 1

4.2. Update

Calculate the updated state mean and the square root of the state covariance:
x ^ k | k = x ^ k | k 1 + K k z k z ^ k | k 1 S k | k = T r i a X k | k 1 K k Z k | k 1 , K k R k
During the model update, calculate the probability of the ξ t h model:
Λ k ξ = P Z k M k ξ , Z k 1 = N γ k j : 0 , S z z , k | k 1 = 1 2 π S z z , k | k 1 exp 1 2 γ k ξ T S z z , k | k 1 1 γ k ξ
Calculate the conditional probability:
μ k ξ = P M k ξ Z k = 1 c Λ k ξ c ¯ ξ
As known from Equations (61) and (66), Q k and R k significantly influence the values of S k | k 1 and S k | k . However, in the target tracking process, the noise data change over time. Any mismatch between the actual noise affecting the system and the noise assumed in the SCKF can degrade the performance of the SCKF, as it may also diverge. Therefore, it is necessary to precisely estimate the matrices Q k and R k . The posterior density function can be calculated as follows:
Π = p X k , Q k , R k Z k = P Z k X k , Q k , R k P X k Q k , R k p Q k , R k P Z k
where p Q k , R k is considered a constant dependent on prior information, and P Z k does not involve an optimization problem; hence, Equation (69) is equivalent to the following:
Π = P Z k X k , Q k , R k P X k Q k , R k p Q k , R k
where P X k | Q k , R k can be calculated using the conditional probability multiplication theorem as follows:
P X k Q k , R k = P x 0 m = 1 k P x m x m 1 , Q k = 1 2 π n / 2 P 0 | 0 1 / 2 exp 1 2 x 0 x ^ 0 P 0 | 0 1 2 × j = 1 k 1 2 π n / 2 Q k 1 / 2 exp 1 2 x m f x m 1 Q k 1 2 = M 1 Q k k / 2 exp 1 2 x 0 x ^ 0 P 0 | 0 1 2 + j = 1 k x m f x m 1 Q k 1 2
where M 1 = 1 2 π n / 2 P 0 | 0 1 / 2 is a constant, n represents the dimension of process noise, and P Z k | X k , Q k , R k can be expressed as follows:
P Z k X k , Q k , R k = j = 1 k P z m x m , R k = j = 1 k 1 2 π q / 2 R k 1 / 2 exp 1 2 z m h z m R k 1 2 = M 2 R k k / 2 exp 1 2 j = 1 k z m h z m R k 1 2
where q represents the dimension of measurement noise and M 2 = 1 2 π q k / 2 is a constant. By substituting Equations (71) and (72) into Equation (70), Π can be calculated:
Π = M 1 M 2 P 0 1 / 2 Q k k / 2 p Q k , R k × exp 1 2 x 0 x ^ 0 P 0 | 0 1 2 + m = 1 k x m f x m 1 Q k 1 2 + m = 1 k z m h z m R k 1 2 = C Q k k / 2 R k k / 2 exp 1 2 m = 1 k x m f x m 1 Q k 1 2 + m = 1 k z m h z m R k 1 2
where C = M 1 M 2 P 0 | 0 1 / 2 p Q k , R k exp 1 2 x 0 x ^ 0 P 0 | 0 1 2 . According to Equation (73), calculate the partial derivatives of the function Π with respect to Q k and R k , where the extremum points are as follows:
Q k = 1 k m = 1 k x ^ m f x ^ m 1 x ^ m f x ^ m 1 T R k = 1 k m = 1 k z m h x ^ m z m h x ^ m T
According to the calculation process of the ASCKF, f x ^ m 1 and h x ^ m can be represented as follows:
Q k = 1 k m = 1 k x ^ m 1 2 n i = 1 2 n f χ i , m 1 | m 1 x ^ m 1 2 n i = 1 2 n f χ i , m 1 | m 1 T R k = 1 k m = 1 k z m 1 2 n i = 1 2 n h χ i , m | m 1 z m 1 2 n i = 1 2 n h χ i , m | m 1 T
By solving the first-order partial derivatives of the function Π with respect to Q k and R k , a suboptimal estimate of the unknown covariance matrices of process and measurement noise can be obtained, which is then incorporated into the nonlinear filter ASCKF for recursive computation. Table 1 gives the algorithm process. Combined with the algorithm proposed in Section 3, this method does not depend on noise statistics and exhibits good robustness. Usually, the detection probability and clutter rate are two important factors that affect tracking performance during the process of target estimation. Therefore, we analyze measurement noise and process noise individually, and then obtain the modification of Q k and R k in the ASCKF. Both the over-estimation and the under-estimation derived from the factors can be corrected. On the other hand, the IMM filtering framework ensures the motion switching between two continuing dynamics, which can distinguish the close targets. Especially, these targets can be further distinguished by using the proposed filter based on the novel filtering mechanism.
This algorithm replaces the covariance matrix recursion of the CKF with the square root form of the covariance matrix, ensuring the efficiency of the filtering process while avoiding numerical instabilities caused by rounding errors, thus enhancing the precision and stability of the filtering algorithm. Moreover, ensuring the positive definiteness of the covariance matrix enhances its numerical stability. The performance of the proposed method, compared to CKF and UKF, features smaller estimation errors and faster convergence rates, and it effectively tracks the motion state of maneuvering targets, enhancing the system’s accuracy.

5. Simulation Results and Discussions

In the numerical study, a typical 2-dimensional scenario is performed to evaluate both the reliability and efficiency of the proposed filter, which is driven by model and data collaboration for extended target tracking. The experimental environment was the following: IntelTM CoreTM i5, RAM 16 GB, WindowsTM 10, and MATLABTM R2022b.

5.1. Simulation Scenario

This algorithm primarily addresses the tracking of extended targets under the conditions of multi-model switching and unknown noise influences. In this section, we present the results of 100 Monte Carlo simulations conducted in the first scenarios. We tested the proposed IMM-TPMB-GBP-ASCKF filter in two target motion scenarios and compared it with the GGIW-TPMB filter from [24], which used clustering and assignment, as well as with the TPMB-BP [20], TPMB-GBP, and IMM-TPMB-GBP-CKF filters.
The linear programming (LP) metric was used to compare the performance of these filters [25,26]. Thus, we set different numbers of maneuvering targets in two scenarios, each under different motion models. The target motion models include Constant Velocity (CV) and Constant Turning (CT), with their corresponding motion equations as follows:
F CV = 1 0 T 0 0 1 0 T 0 0 1 0 0 0 0 1 ,   F CT = 1 0 sin ( θ T ) / θ [ 1 cos ( θ T ) ] / θ 0 1 [ 1 cos ( θ T ) ] / θ sin ( θ T ) / θ 0 0 cos ( θ T ) / θ sin ( θ T ) / θ 0 0 sin ( θ T ) / θ cos ( θ T ) / θ ,   Q CV = Q CT = T 2 / 3 0 T 2 / 2 0 0 T 2 / 3 0 T 2 / 2 T 2 / 2 0 T 0 0 T 2 / 2 0 T σ 2
Each object consists of a two-dimensional position and velocity, which we represent using the target matrix x k = [ x k , x ˙ k , y k , y ˙ k ] T . ( x k , y k ) denotes the plane position coordinates, and ( x ˙ k , y ˙ k ) is used to represent the planar position. We consider a monitoring area of 200   m × 200   m , with a total motion duration of 100   s . For the following two scenarios, we set the sampling period T = 1   s and the mean clutter rate γ = 10 . The Poisson rate of PPP birth is 0.01, the clutter is distributed in the region. For the birth density, the velocity follows a Gaussian distribution with zero mean and covariance 100 I 2 , where I 2 represents an identity matrix of size 2. For the GGIW implementations, the position is Gaussian distributed with zero mean and covariance 150 I 2 . In the two experimental scenarios, different numbers of target trajectories are set to test the algorithm’s performance. Assume a detection probability of P D = 0.9 , a survival probability of P S = 0.95 , and a forgetting factor for the measurement rate state set to η = 1.25 . The experimental results indicate that using 1000 particles in our filter achieves a good balance between the runtime and estimation performance.

5.2. Simulation 1

For Simulation 1, we set up a challenging experiment to demonstrate the tracking error of target trajectories and the precision of extended state tracking under a fixed clutter rate, 10, by four filters. The parameters of target trajectories are shown in Table 2. Note that the targets in three different motion states appear at different locations within the detection area at the 10th s, 20th s, and 30th s. With a total of 100 steps run, each trajectory initialized from the starting point to the midway step, with the remaining steps generated forwards and backwards.
In this scenario, we concentrate on estimating trajectories and compare the estimated extended states against the actual values. Among the various performance metrics for trajectory estimation, we employ the Linear Programming (LP) and the Gaussian Wasserstein Distance (GWD) metric d · , · [27]. We set the cut-off distance to 20, establish the order p = 1 , and set the track switch cost to 2 in order to apply these metrics at each time step.
The target trajectories are shown in Figure 2. We know that the available measurement information contains both real targets and random clutter. From this figure, we observe three trajectories. Clearly, the three targets have nonlinear motion states, performing CV and CT motions in stages within the detection area and the model has different angular velocities during CT motion. During three stages, the respective motion states are designated as CV, CT, and CV in turn, where the turning velocity is given before the scenario. The clutter obeys Poisson distribution, and the measurement noise and covariance are unknown to the filter.
In Figure 3a, it is evident that the filters proposed in this paper offer excellent extended state estimation and stable tracking performance for elliptical extended targets. In Figure 3b, we note that the cardinality statistics change as time progresses. Obviously, all algorithms perform well in the potential estimation. The TPMB-BP, TPMB-GBP, IMM-TPMB-GBP-CKF algorithms, and the algorithm discussed in this paper all closely approximate the true number of targets, while the classical GGIW-TPMB algorithm underestimates the number of highly maneuverable targets; significant deviations in the potential estimation by the GGIW-TPMB algorithm from the true values are evident at the 30th s, 40th s, and 65th s. The TPMB-BP and TPMB-GBP algorithms have a faster computation speed compared to IMM-TPMB-BP-CKF, and both algorithms essentially provide similar estimates of target potential. Notably, the standard filters are affected by multiple measurements from random clutter, resulting in some positional estimation biases. The filters proposed in this paper effectively eliminate the impact of unknown noise on tracking performance, providing accurate estimates.
The track error of the trajectory in Scenario 1 under the LP metric is presented in Table 3. For the extended target, the cost of state error dominates the overall LP estimation error. When the targets are in the same clutter, the measurement errors increase. For the most advanced TPMB-BP filter, we improved its correlation process and obtained the Gate-control BP association method that reduces the LP measurement error with adaptive gating. Among the five filters, the IMM-TPMB-GBP-ASCKF filter, with its superior noise adaptability, achieves the lowest LP metric. All five filters perform well when tracking targets at adequate intervals. However, as the targets move closer to each other, they become indistinguishable. Thanks to its adaptive noise estimation, the improved ASCKF enables the algorithm proposed in this study to achieve the lowest target state cost and false detection cost among the four filters. In such scenarios, the cost of losing targets predominates in the LP estimation error. Moreover, in situations where highly maneuvering targets frequently change their motion states, the improved algorithm allows the trajectory to achieve reduced switch cost during transitions.
Undoubtedly, the TPMB-BP filter is the current optimal filter, which avoids explicit assumption enumeration in data association compared to traditional filters, improves data association accuracy, and reduces the running time of the data association process. Our work is to further reduce the running time of the data association process and improve the anti-interference ability against noise based on the TPMB-BP filter, and its performance is demonstrated in Simulation 2.

5.3. Simulation 2

For Scenario 2, to further validate the performance of the proposed algorithm, we designed a set of comparative experiments under the low clutter rate of 10 and the high clutter rate of 30, respectively. The parameters of target trajectories are shown in Table 4.
Figure 4 shows that all five filters can converge to the actual number of targets. It was observed that as the clutter density increases, the state estimation error of each algorithm also increases. The reason for the increase in state estimation error in standard filters is due to the generation of random clutter by the filter and environment while tracking maneuvering targets. In contrast, the other four filters exhibit better robustness to the state transitions of targets, resulting in more accurate state estimations. More importantly, this paper employs a method using interactive gating with GBP data association and an ASCKF that adapts the covariance matrix for unknown noise, addressing inaccuracies in state estimation under high clutter and unknown noise conditions. However, the cardinality estimated by standard filters is the least stable. During low clutter phases, false alarms occur due to the filter’s inability to adaptively model estimates for non-linear targets. Since the targets maintain non-maneuverability, standard filters mix random clutter and unknown noise, leading to inaccurate state estimation. Conversely, at the 55th s mark, there is a missed detection because two of the targets switch from CT to CV motion models, and the filters fail to immediately adapt to the motion changes. The IMM-TPMB-GBP-CKF, TPMB-GBP, and TPMB-BP algorithms exhibit lower state estimation errors and show stable performance throughout the detection period, though they also face the issues of increased estimation error with higher clutter rates. In contrast, the algorithm proposed in this work maintains low state estimation errors under both low and high clutter rates. It is less affected by increases in clutter and motion model transitions and controls the overestimation of target states under high clutter conditions.
The GOSPA distances for the five filters are shown in Figure 5. As expected, throughout the detection period, the algorithm proposed in this paper consistently shows the lowest GOSPA distance. At the 10th s, 20th s, 45th s, and 65th s, the peaks in clutter density from the standard filter are caused by false alarms, mode transitions, and missed targets, involving both cardinality and location errors. The experiments reveal that due to the inability of the filters to switch motion models during the target’s motion cycle, significant tracking mismatches occur, leading to increased state estimation errors. Since the targets include CV and CT models, which appear with equal probability at certain times, the proposed filter automatically adjusts the gating thresholds based on the target state model, reducing the number of incorrect associations and association duration, thereby enhancing the accuracy of target measurement associations and correcting unstable cardinality estimates. Considering the actual tracking reliability of the filter proposed in this paper under different clutter rates, four types of filters were operated under two different clutter rates. This figure displays a comparison chart of the GOSPA errors of the five filters under these two clutter rates. Notably, as the clutter rate increases, the filter proposed in this paper provides stable cardinality estimates. In situations with noise and unknown covariance, the algorithm proposed demonstrates stable estimation performance. Moreover, in the ASCKF, it also performs well with noise having smaller covariance. This means that the filter proposed in this paper does not exhibit significant changes under various clutter environments.
Finally, due to the specificity of the GOSPA evaluation algorithm, it can be differentiated into two parts: the number of missed targets and the number of false targets. The missed and false target state estimates are compared under different clutter rates. Figure 6 shows the missed target state estimation under different clutter rates. It can be seen that all four filters experience a sharp increase in the number of missed targets at the 10th s, as targets begin to appear, due to the classical filter’s lack of a multimodal assumption and the increasing number of targets under the influence of clutter. However, the remaining three filters, having a multimodal assumption, adjust their matching models with changes in the motion state of the models, thus significantly reducing missed targets compared to the classical filter. Compared with TPMB-GBP and TPMB-BP, the improved TPMB-GBP reduced the GOSPA error. Due to the lack of a mixture of adaptive multi-model algorithm for high maneuvering targets, the performance difference between the two algorithms is relatively small.
Table 5 shows that TBP-GBP has shorter running time under different clutter and particles. The IMM-TPMB-GBP-CKF, TPMB-GBP, and TPMB-BP algorithms show an increase in missed counts under different clutter levels, but the IMM-TPMB-GBP-CKF algorithm reduces the number of misses in high clutter rates due to adaptive gating adjustments and also reduces operation time. The algorithm proposed in this paper, with the inclusion of ASCKF, reduces target location errors and clutter interference, and the number of missed targets does not significantly change with increasing clutter.
The false target state estimation, shown in Figure 7a, shows a significant increase at the 10th s as the number of emerging targets grows. It can be observed that at the 20th s, as the target motion models switch, the number of false targets reaches its peak, with the algorithm proposed in this paper having the lowest false target state estimation among the five filters. Studies show that, compared to the four filters analyzed alongside the proposed algorithm, this algorithm exhibits the lowest GOSPA error values and adapts to non-linear target estimation under the influence of different clutter rates.
In order to verify the average running time of the proposed filter in this article under different particle and clutter conditions, we set up a set of comparative experiments to verify its effectiveness. Due to the fact that the GGIW-TPMB filter is a Gaussian implementation method, we only compare the average running time of the other four filters under different particles. Note that TPMB-BP is a particle implementation based on the confidence propagation data association, which results in a large computational workload due to multiple iterations of calculating example weights. On the basis of the same tracking effect, the proposed algorithm reduces the computation of particle likelihood in the target measurement matching process by integrating the TPMB-GBP algorithm for comparison, thereby reducing the overall time. With the increase in trajectory, iteration count, and clutter rate, the time consumption of the IMM-TPMB-GBP-ASCKF filter is only slightly higher than that of the IMM-TPMB-GBP-CKF filter, which is reasonable due to the additional noise estimation process of the proposed filter. As the noise increases, the running time of the four filters significantly increases. Due to the improved ASCKF filter’s good adaptability to noise, it can effectively filter out clutter, and the time consumption does not change significantly compared to CKF. Additionally, under high clutter density conditions, the IMM-TPMB-GBP-ASCKF filter demonstrates higher accuracy. The results indicate that the proposed IMM-TPMB-GBP-ASCKF filter maintains robust tracking performance for highly maneuverable targets and under the impact of unknown noise, while incurring a moderate computational cost.

6. Conclusions

This paper introduces a new robust TPMB filter designed for conditions with uncertain process and measurement noise. By enhancing the adaptive multi-model switching method and performing process noise estimation, we provide robust recursive formulations for TPMB. In terms of model-driven aspects, we have improved the IMM by increasing the matching degree of target motion mode switching through model probability transformation at two adjacent moments. From a data-driven perspective, the GBP in the factor graph for the association of measurement vectors allows the model to adaptively adjust the gating threshold based on the association vector, effectively reducing the computational complexity of data association. Additionally, we present an improved method of importance sampling for IMM-TPMB’s ASCKF particles, optimizing the particle degradation issue by estimating the process and measurement noise in the context of target positioning. Simulation studies show that, compared to the GGIW-TPMB, TPMB-BP, TPMB-GBP, and IMM-TPMB-GBP-CKF, the IMM-TPMB-GBP-ASCKF achieves the best overall performance in terms of filtering efficiency and computation time across different target numbers, motion states, and clutter rates.
Future work includes conducting more extensive tests on different clutter models. Although the algorithm proposed in this paper mitigates particle degradation and running time to some extent, the running time continues to increase with the number of targets, and measurement estimation errors arise due to closely spaced measurements. For scenarios with a high number of measurements and closely spaced measurements, an improved measurement tracking gating method can be employed for processing.

Author Contributions

Conceptualization, Y.Z. and B.L.; methodology, Y.Z. and B.L.; software, J.Z. and Z.L. (Zhengyuan Li); validation, Z.L. (Zhikang Li); formal analysis, Y.Z. and B.L.; investigation, Y.Z.; resources, Y.Z.; data curation, Z.L. (Zhikang Li); writing—original draft preparation Y.Z.; writing—review and editing, Y.Z. and Z.L. (Zhikang Li); visualization, J.Z. and Z.L. (Zhengyuan Li); supervision, J.Z. and Z.L. (Zhengyuan Li); project administration, J.Z.; funding acquisition, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research work was supported in part by the Fundamental Research Project of the Educational Department of Liaoning Province (No. JYTMS20230862), and the Natural Science Foundation of Liaoning Province (No. 2020-MS-292).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ren, X.Y.; Huang, Z.P.; Sun, S.Y.; Liu, D.Y.; Wu, J.L. An efficient MHT implementation using GRASP. IEEE Trans. Aerosp. Electron Syst. 2014, 50, 86–101. [Google Scholar] [CrossRef]
  2. Ishrat, M.; Liang, Q.L. Multistep information fusion for target detection using UWB radar sensor network. IEEE Sens. J. 2015, 15, 5927–5937. [Google Scholar]
  3. Li, J.; Zhan, W.; Hu, Y.; Tomizuka, M. Generic tracking and probabilistic prediction framework and its application in autonomous driving. IEEE Trans. Intell. Transp. Syst. 2020, 21, 3634–3649. [Google Scholar] [CrossRef]
  4. Ravindran, R.; Santora, M.J.; Jamali, M.M. Multi-object detection and tracking, based on DNN, for autonomous vehicles: A Review. IEEE Sens. J. 2021, 21, 5668–5677. [Google Scholar] [CrossRef]
  5. Chen, X.; Li, Y.; Li, Y.X.; Yu, J. PHD and CPHD algorithms based on a novel detection probability applied in an active sonar tracking system. Appl. Sci. 2017, 8, 36. [Google Scholar] [CrossRef]
  6. Zhang, L.; Lan, J. Extended object tracking using random matrix with skewness. IEEE Trans. Signal Process. 2020, 68, 5107–5121. [Google Scholar] [CrossRef]
  7. Koch, J.W. Bayesian approach to extended object and cluster tracking using random matrices. IEEE Trans. Aerosp. Elecrton Syst. 2008, 44, 1042–1059. [Google Scholar] [CrossRef]
  8. Honer, J.; Kaulbersch, H. Bayesian extended target tracking with automotive radar using learned spatial distribution models. In Proceedings of the 2020 IEEE International Conference on Multi-sensor Fusion and Integration for Intelligent Systems (MFI), Karlsruhe, Germany, 14–16 September 2020; pp. 316–322. [Google Scholar]
  9. Knill, C.; Scheel, A.; Dietmayer, K. A direct scattering model for tracking vehicles with high-resolution radars. In Proceedings of the 2016 IEEE Intelligent Vehicles Symposium (IV), Gothenburg, Sweden, 19–22 June 2016; pp. 298–303. [Google Scholar]
  10. Meyer, F.; Williams, J.L. Scalable detection and tracking of geometric extended objects. IEEE Trans. Signal Process. 2021, 69, 6283–6298. [Google Scholar] [CrossRef]
  11. Granström, K.; Fatemi, M.; Svensson, L. Poisson multi-Bernoulli mixture conjugate prior for multiple extended target filtering. IEEE Trans. Aerosp. Electron Syst. 2019, 56, 208–225. [Google Scholar] [CrossRef]
  12. Wang, Y.; Rao, P.; Chen, X. Robust PMBM filter with unknown detection probability based on feature estimation. Sensors 2022, 22, 3730. [Google Scholar] [CrossRef] [PubMed]
  13. Li, G.C.; Kong, L.J.; Yi, W.; Li, X.L. Multiple model Poisson multi-Bernoulli mixture filter for maneuvering targets. IEEE Sens. J. 2020, 21, 3143–3154. [Google Scholar] [CrossRef]
  14. Chen, Y.M.; Liu, W.F.; Wang, X.D. Multiple extended target tracking based on GLMB filter and Gibbs sampler. In Proceedings of the 2017 International Conference on Control, Automation and Information Sciences (ICCAIS), Chiang Mai, Thailand, 31 October–1 November 2017; pp. 26–31. [Google Scholar]
  15. Lu, X.C.; Jing, D.H.; Jiang, D.F.; Gao, Y.Y.; Yang, J.L.; Li, Y.; Li, W.D.; Tao, J.; Liu, M. Trajectory PHD filter for adaptive measurement noise covariance based on variational Bayesian approximation. Appl. Sci. 2022, 12, 6388. [Google Scholar] [CrossRef]
  16. Zhou, Y.Y.; Sun, C.; Xie, L. A trajectory-oriented Poisson Multi-Bernoulli mixture method for matched field tracking to achieve trajectory continuity. In Proceedings of the OCEANS 2022, Chennai, India, 21–24 February 2022; pp. 1–4. [Google Scholar]
  17. Xia, Y.X.; García-Fernández, Á.F.; Svensson, L. An efficient implementation of the extended object Trajectory PMB filter using blocked Gibbs sampling. In Proceedings of the 2023 26th International Conference on Information Fusion (FUSION), Charleston, SC, USA, 27–30 June 2023; pp. 1–8. [Google Scholar]
  18. Hu, X.L.; Zhang, Q.; Song, B.J.; Zhao, M.X.; Xia, Z.Q. Student-t mixture GLMB filter with heavy-tailed noises. In Proceedings of the 2022 IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC), Xi’an, China, 25–27 October 2022; pp. 1–6. [Google Scholar]
  19. Meyer, F.; Win, M.Z. Scalable data association for extended object tracking. IEEE Trans. Signal Inf. Process. Over Netw. 2020, 6, 491–507. [Google Scholar] [CrossRef]
  20. Xia, Y.X.; García-Fernández, Á.F.; Meyer, F.; Williams, J.L.; Granström, K.; Svensson, L. Trajectory PMB filters for extended object tracking using belief propagation. IEEE Trans. Aerosp. Electron Syst. 2023, 59, 9312–9331. [Google Scholar] [CrossRef]
  21. Liu, M.T.; Zeng, C.; Li, S.D.; Tao, H.H.; Liao, G.S.; Li, J. Robust TPMB filtering using sensors with limited sensing range under nonuniform clutter background. IEEE Sens. J. 2022, 23, 470–483. [Google Scholar] [CrossRef]
  22. Luo, Y.L.; Li, Z.M.; Liao, Y.R.; Wang, H.N.; Ni, S.Y. Adaptive Markov IMM based multiple fading factors strong tracking CKF for maneuvering hypersonic-target tracking. Appl. Sci. 2022, 12, 10395. [Google Scholar] [CrossRef]
  23. Song, R.; Chen, X.Y.; Fang, Y.C.; Huang, H.Q. Integrated navigation of GPS/INS based on fusion of recursive maximum likelihood IMM and square-root cubature Kalman filter. ISA Trans. 2020, 105, 387–395. [Google Scholar] [CrossRef] [PubMed]
  24. Yao, W.L.; Zhuang, G.T.; Li, B.Y.; Chen, A.L. Autonomous navigation based on information filter multi-model SCKF. In Proceedings of the 2017 Chinese Automation Congress (CAC), Jinan, China, 20–22 October 2017; pp. 4524–4527. [Google Scholar]
  25. Xia, Y.X.; Granström, K.; Svensson, L.; García-Fernández, Á.F.; Williams, J.L. Extended target Poisson multi-Bernoulli mixture trackers based on sets of trajectories. In Proceedings of the 22th International Conference on Information Fusion, Ottawa, ON, Canada, 2–5 July 2019; pp. 1–8. [Google Scholar]
  26. García-Fernández, Á.F.; Rahmathullah, A.S.; Svensson, L. A metric on the space of finite sets of trajectories for evaluation of multi-target tracking algorithms. IEEE Trans. Signal. Process. 2020, 68, 3917–3928. [Google Scholar] [CrossRef]
  27. Yang, S.; Baum, M.; Granström, K. Metrics for performance evaluation of elliptic extended object tracking methods. In Proceedings of the 2016 IEEE International Conference on Multisensor Fusion and Integration for Intelligent Systems (MFI), Baden-Baden, Germany, 19–21 September 2016; pp. 523–528. [Google Scholar]
Figure 1. Factor graph of the factorization of Equation (24), with circles and squares representing factor nodes and variable nodes.
Figure 1. Factor graph of the factorization of Equation (24), with circles and squares representing factor nodes and variable nodes.
Applsci 14 07201 g001
Figure 2. Tracks of 3-target estimates.
Figure 2. Tracks of 3-target estimates.
Applsci 14 07201 g002
Figure 3. Target extended state estimation at the 40th s and cardinality statistics for 5 filters.
Figure 3. Target extended state estimation at the 40th s and cardinality statistics for 5 filters.
Applsci 14 07201 g003
Figure 4. Cardinality statics for four filters under low and high clutter densities.
Figure 4. Cardinality statics for four filters under low and high clutter densities.
Applsci 14 07201 g004
Figure 5. GOSPA error for four filters under low and high clutter densities.
Figure 5. GOSPA error for four filters under low and high clutter densities.
Applsci 14 07201 g005
Figure 6. Number of missed targets for four filters under low and high clutter densities.
Figure 6. Number of missed targets for four filters under low and high clutter densities.
Applsci 14 07201 g006
Figure 7. Number of false targets for four filters under low and high clutter rates.
Figure 7. Number of false targets for four filters under low and high clutter rates.
Applsci 14 07201 g007
Table 1. Algorithm process.
Table 1. Algorithm process.
Pseudocode of IMM-TPMB-GBP-ASCKF Filter
Applsci 14 07201 i001
Table 2. Parameters of target trajectories in Simulation 1.
Table 2. Parameters of target trajectories in Simulation 1.
No.Prior Target StateFirst StageSecond StageThird Stage
1(15, 2.1, 17, 2.5)T20~40 s40~63 s--
2(150, −1.5, 200, −1)T10~40 s25~75 s--
3(150, 0, 18, 2)T35~55 s55~90 s90~100 s
Table 3. Mean tracking errors for each filter.
Table 3. Mean tracking errors for each filter.
FilterTotalStateMissFalseSwitch
GGIW-TPMB52.6543.4111.4117.390.158
TPMB-BP44.3538.756.2713.740.076
TPMB-GBP43.7936.926.1313.490.074
IMM-TPMB-GBP-CKF40.1234.145.8414.010.097
IMM-TPMB-GBP-ASCKF37.6033.743.738.680.079
Table 4. Averaged tracking errors for each filter.
Table 4. Averaged tracking errors for each filter.
NumberPrior Target StateMotionLasting TimeDeath Position
1(100, −1.125, 66, −0.55)TCV(3,83)(10,22)
2(138, −0.95, 152, 0.125)TCV(3,83)(62,162)
3(51, 0.75, 4, 1.25)TCV(6,86)(111,104)
4(57, −0.225, 15, 0.55)TCT(6,86)(23,135)
5(133, −0.25, 41, 0.2)TCT(9,89)(33,141)
6(3, 0.175, 12, 0.75)TCT(9,89)(33,112)
7(29, 1.5, 133, −0.75)TCT-CV(12,92)(170,115)
8(34, 0.6, 9, 0.55)TCT-CV(12,92)(94,101)
9(21, 0.25, 16, 0.6)TCT-CV(15,95)(112,116)
10(24, 0.125, 141, −0.75)TCT-CV(15,95)(35,53)
Table 5. Comparisons of averaged running time under different clutter levels and particle counts in Simulation 2.
Table 5. Comparisons of averaged running time under different clutter levels and particle counts in Simulation 2.
Particle NumberClutter RateTPMB-BPTPMB-GBPIMM-TPMB-GBP-CKFIMM-TPMB-GBP-ASCKF
500low20.63519.229171.960176.263
high34.87130.615179.284184.736
1000low29.25627.836174.642178.918
high50.69747.698193.856197.604
1500low39.19833.651181.284182.684
high59.98954.873205.207208.837
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhou, Y.; Li, B.; Zhang, J.; Li, Z.; Li, Z. Tracking Extended Targets: Novel TPMB Filter Driven by Model and Data Collaboration. Appl. Sci. 2024, 14, 7201. https://doi.org/10.3390/app14167201

AMA Style

Zhou Y, Li B, Zhang J, Li Z, Li Z. Tracking Extended Targets: Novel TPMB Filter Driven by Model and Data Collaboration. Applied Sciences. 2024; 14(16):7201. https://doi.org/10.3390/app14167201

Chicago/Turabian Style

Zhou, Yubin, Bo Li, Jinyu Zhang, Zhikang Li, and Zhengyuan Li. 2024. "Tracking Extended Targets: Novel TPMB Filter Driven by Model and Data Collaboration" Applied Sciences 14, no. 16: 7201. https://doi.org/10.3390/app14167201

APA Style

Zhou, Y., Li, B., Zhang, J., Li, Z., & Li, Z. (2024). Tracking Extended Targets: Novel TPMB Filter Driven by Model and Data Collaboration. Applied Sciences, 14(16), 7201. https://doi.org/10.3390/app14167201

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop