Next Article in Journal
AHSWFM: Automated and Hierarchical Surface Water Fraction Mapping for Small Water Bodies Using Sentinel-2 Images
Next Article in Special Issue
Data-Driven Joint Beam Selection and Power Allocation for Multiple Target Tracking
Previous Article in Journal
Unsupervised Infrared Small-Object-Detection Approach of Spatial–Temporal Patch Tensor and Object Selection
Previous Article in Special Issue
Smoothing Linear Multi-Target Tracking Using Integrated Track Splitting Filter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tracking of Maneuvering Extended Target Using Modified Variable Structure Multiple-Model Based on Adaptive Grid Best Model Augmentation

1
School of Information Engineering, Henan University of Science and Technology, Luoyang 471023, China
2
School of Information and Communication Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
3
School of Automation, Nanjing University of Science and Technology, Nanjing 210094, China
*
Author to whom correspondence should be addressed.
Submission received: 20 February 2022 / Revised: 19 March 2022 / Accepted: 23 March 2022 / Published: 28 March 2022
(This article belongs to the Special Issue Radar Signal Processing for Target Tracking)

Abstract

:
Maneuvering extended target tracking, an important but challenging research field, has attracted increasing attention in the field of radar signal processing. Variable structure multiple-model (VSMM) estimation is the current mainstream tracking algorithm, it possesses high tracking performance but depends largely on the model set used. However, the existing model set design methods still cannot offer higher tracking performance due to the complex maneuverability of the target. The current best model augmentation (BMA) algorithm is an efficient and universal model set design method, but it cannot adapt well to complex maneuvering situations because of its fixed basic and candidate model set. Hence, this paper proposes a modified BMA algorithm for VSMM based on adaptive grid, it can make full use of the previous information to realize an adaptation and time variation of the model set. Overall, two different scenarios were considered, wherein the digital simulation shows that the proposed method can provide a more reliable estimation of the kinematic state and shape of the extended target.

1. Introduction

With the emergence of the new generation of aviation weapons, such as the fourth-generation fighters and unmanned combat aircraft, modern warfare will involve targets with stronger mobility, increasingly sophisticated escape methods, and low detection capabilities [1]. To achieve accurate guidance and attack solutions for this type of strong maneuvering target, it is essential to estimate its state quickly and accurately. Meanwhile, modern radar sensors have shown the characteristics of high resolution and high detection accuracy, multiple measurements can be captured on the target surface at each sampling moment, thus the shape cannot be ignored anymore [2]. That is to say, the state that needs to be estimated includes not only its kinematic state but also its extension state. In this paper, a new model set design method for the variable structure multiple-model estimation has been proposed to improve the tracking performance of the maneuvering extended target in air combats.
Any target tracking algorithm is built on the basis of a dynamic model. Selecting the appropriate dynamic models is an important entry point for target tracking, which can effectively improve its estimation performance [3]. The development of a dynamic model is a process from simple to complex. In order to accommodate the complex maneuverability of the target, scholars have proposed models, such as the constant velocity (CV) model [4], the constant acceleration (CA) model [4], the constant turn (CT) model [5], the semi-Markov model [6], and the “current” statistics (CS) model [7]. However, in the field of modern target tracking, a single mathematical model cannot be used for describing the kinematic state due to its complexity and variability.
Blom and Bar-Shalom first proposed the interactive multiple-model (IMM) algorithm [8]. IMM is a fixed structure multiple-model (FSMM) algorithm. At any time, it uses several fixed models to perform parallel filtering, and the state estimation result is the weighted value of each output. Clearly, it is difficult to make use of the real-time information of the system. Only when the model set matches the true mode of the current system, it can be considered as the optimal estimation. Once the true system mode deviates from the current model structure, it might cause tracking failure. In order to obtain high accuracy, a large number of fixed models are often required to implement coverage of the true system mode. However, this inevitably will increase the amount of computation and lead to competition between models, thus resulting in decreased tracking accuracy [9].
In order to overcome this dilemma, Professor Li and colleagues proposed the variable structure multiple-model (VSMM) algorithm in the 1990s for the first time [10]. Compared to the pre-existing FSMM algorithm, the VSMM algorithm inherits the estimation and fusion steps of the multiple-model (MM) algorithm but involves an additional step of model set design. That is to say, its recursive process involves two tasks: model set adaptation (MSA) and MM sequence condition estimation, respectively. Among them, the purpose of MSA is to realize the dual self-adaptation of the model structure and the model number. The latter is responsible for estimating the target state by using multi-model estimation methods under the conditions of a given model set. It is obvious that the performance of the MM estimation algorithm depends largely on the model set used. Therefore, the research on model set design methods is of great significance.
The existing model set design methods are mainly divided into two aspects: (i) to design an over-complete model set and subsequently perform adaptive switching and selection among the model sets [11], and (ii) to generate a new model set online in real-time based on the target state estimation [12,13,14].
The adaptive strategies of the model group switching (MGS) algorithm and the likely-model set (LMS) algorithm adopt the first design idea. MGS adopts a model activation and termination approach, but its tracking performance is affected by the division and topology of the model set. Once the target maneuvers too fast, its performance will become poor. LMS can certainly reduce some computation while retaining the estimation accuracy, but it cannot activate models outside the total model set. Therefore, when the target maneuvers, there will be a high peak error. The expected-mode augmentation (EMA) algorithm is a model activation method, which adopts the second design idea. It activates an expected model, which is a probability-weighted sum of the basic models. The tracking performance of EMA is obviously better than that of the IMM algorithm, but it is only suitable for the model set with parameter change and cannot be directly used in the model set with structural change. The adaptive grid (AG) algorithm is a model self-learning algorithm in essence. It does not need to set the corresponding models in advance. The parameters of each model at k 1 are weighted to generate new model parameters at k to realize the time variation of the model set. Such algorithms and designs have mainly been included in a few previous studies [15,16,17].
Lan et al. [18] first proposed the best model augmentation (BMA) algorithm, and its adaptive strategy was realized based on Kullback–Leibler (KL) information [19]. In this algorithm, the set of basic models, as well as candidate models, have been assumed to be known. The algorithm activates a “best” model from the candidate model set and adds it to the basic model set. The “best” here refers to the minimum difference between the model and the true system mode, which is measured by KL divergence. The smaller the value of KL, the smaller is the difference. BMA overcomes the defects of EMA and is widely used in situations with different model structures and parameters. However, there are still several problems that can be improved: (i) The KL criterion need to be frequently calculated each time the best model is activated from the candidate model set, resulting in a large computational burden. (ii) The models in the basic model set and the candidate model set are fixed, which means that the model parameters cannot be adjusted in real time according to the maneuvering state. (iii) The prior known information of the model set is introduced, that is, all dynamic models are known, thus the correction capability of the best model to the basic model set is limited.
To deal with the problem, this paper proposes a new model set design method by modifying the BMA algorithm. The main idea of the proposed method is as follows: Firstly, the basic model set is represented by a coarse grid without any prior setting, and the corresponding adaptive mechanism is applied to realize the time variation of the basic model set. Secondly, a uniformly quantized mode space is generated, centered on the best model at the above moment, and this is taken as the candidate model set at this moment. Due to the time-varying model set, the uncertainty caused by the assumption of prior knowledge in the traditional method can be effectively avoided, and the correction ability of the best model to the basic model set can also be further improved. Finally, the complex KL criterion in the traditional BMA algorithm is discarded, and a simpler mechanism based on similarity distance is adopted to select the best model. In the simulation experiment part, the single maneuvering extended target has been taken as the tracking target, and a digital simulation has been carried out in two scenarios. The simulation results show that this method can effectively avoid tracking failure and improve the estimation accuracy under strong maneuvering conditions.
The remainder of this paper is organized as follows. In the next section, we mainly complete the problem formulation and propose the modified BMA algorithm for VSMM based on adaptive grid. Section 3 does the numerical simulation about the proposed method, wherein its tracking performances are illustrated in comparison to some existing approaches. Section 4 discusses the experimental results. Section 5 concludes the paper.

2. VSMM Algorithm Based on AG-BMA

2.1. Problem Formulation

The task of maneuvering extended target tracking involves estimating the target state based on the constructed mathematical model and the measured value obtained by the sensor. In this paper, the establishment of the mathematical model and the generation of the measurement points have been carried out in the two-dimensional Cartesian coordinate system. In addition, this paper mainly focuses on designing the model set for extended target tracking and proposes a modified best model augmentation algorithm based on adaptive grid (AG-BMA), in which the clutter and missed detection problems are not involved.
Assuming that the state vector, x k = x k m , x k e T , of the target includes the motion variable, x k m , and the shape variable, x k e . The kinematic state of the target is represented by a random variable x k m = x , v x , y , v y T , where x , y T and v x , v y T represent the position and velocity of the target centroid, respectively.
For a linear jump Markov system, the transfer process between states can be formulated as follows:
x k + 1 m x k + 1 e = F k m 0 0 F k e x k m x k e + Γ k m i 0 0 Γ k e i w k m w k e , k N , i = 1 , 2 , r
where F k m and F k e represent the state transition matrix of the motion variable and the morphological variable, Γ k m i and Γ k e i represent the covariance matrices of the process noise, w k m N w ¯ k m , Q k m and w k e N w ¯ k e , Q k e are the independent process noise.

2.1.1. Star-Convex Extended Target

In many real applications, there is a pressing need for incorporating the target shape information into target tracking to improve performance. In addition, knowledge of the target shape is also beneficial for target classification and identification. For extended target tracking, many approaches mainly differ in target shape models for the unknown true shape. However, they usually consider using some basic shapes, such as ellipse, rectangle, and so on. In recent years, there is a novel extended target model which describes the target approximately using star-convex shape [20]. It has an advantage over the other shape modeling methods, because more detailed shape information can be described and utilized, e.g., the real shape of the aircraft can be approximated by the star-convex more accurately than other shapes.
In this paper, the random hypersurface model (RHM) has been used to represent the shape parameters of a star-convex extended target [21]. By changing these parameters, the shape can be fitted arbitrarily [22]. The contour, S x k , of the star-convex extended target can be described by a radial function, r ϕ , the size of which represents the distance between each contour point and the centroid x k c [23]. The details of the model are shown in Figure 1.
Then, the target contour can be described as follows:
S ¯ x k = s ^ k r ϕ e ϕ + x k c ϕ 0 , 2 π
where the scaling factor, s ^ k , is uniformly distributed in the interval 0 , 1 , ϕ represents the angle formed by the ray from the centroid, x k c , to the contour point and the X-axis, e ϕ = cos ϕ sin ϕ T .
The radial function is expanded as a Fourier series, resulting in the following expressions:
r B k , ϕ = a k 0 + j = 1 , , N F a k j cos j ϕ + b k j sin j ϕ = R ϕ B k
B k = a k 0 , a k 1 , b k 1 , , a k N F , b k N F T
R ϕ = 1 , cos ϕ , sin ϕ , , cos N F ϕ , sin N F ϕ
Therefore, the shape parameters of the target can be obtained after parameterization by the relevant Fourier coefficients.
x k e = a k 0 , a k 1 , b k 1 , , a k N F , b k N F T

2.1.2. Measurement Model

The measurement sources drawn from the surface of the extended target are uniformly distributed. Additionally, it is the location of the measurement source with noise. In this paper, the measurement source has been modeled with RHM. It has been assumed that each scan will result in random n k independent position measurements, which can be defined as follows:
z k = H k x k + υ k , k N , i = 1 , 2 , r
where H k is the observation matrix, z k = z x , z y T represents the position of measurement points, and υ k N υ ¯ k , R k is the independent measurement noise. Both υ k and w k belong to the unrelated zero-mean Gaussian white noise sequences.

2.1.3. Model Jump Sequences

The jump between models follows the first-order Markov process, and the Markov chain is homogeneous. Assuming that the model matching mode, s i , at k 1 is m i and the model matching mode, s j , at k is m j , then the transition probability from the model m i to the model m j is denoted as π i j .
P m k = m j | m k 1 = m i = π i j , k N , i , j = 1 , 2 , , r
The implementation steps of the VSMM algorithm, based on AG-BMA, proposed in this work are shown in Figure 2. In this chapter, the tasks to be completed can be divided into the following two parts:
(1)
MSA M k = M k b m k b e s t . where M k b is the basic model set, m k b e s t is the activated best model, and M k represents the valid model set used for MM estimation;
(2)
MM estimation M k , M k 1 . It is responsible for a recursive calculation of the model set M k 1 and M k at k 1 and k, respectively.
Figure 2. Framework of the VSMM algorithm based on AG-BMA.
Figure 2. Framework of the VSMM algorithm based on AG-BMA.
Remotesensing 14 01613 g002
The former is responsible for determining a set of models for MM estimation, whereas the latter is responsible for estimating the target state under the conditions of a given model set. Next, the detailed process of MSA and MM estimation will be given in Section 2.2 and Section 2.3, respectively.

2.2. MSA

When the AG-BMA method is used for the model set design, its recursive flow can be described, as shown in Figure 3. where M k 2 b , M k 1 b , and M k b are the basic model set at k 2 , k 1 , and k, respectively. M k 2 c , M k 1 c , and M k c represent the candidate model set at k 2 , k 1 , and k, respectively. m k 2 b e s t , m k 1 b e s t , and m k b e s t denote the best model activated from M k 2 c , M k 1 c , and M k c , respectively. M k 2 , M k 1 , and M k represent the valid model set.
In this section, the adaptive methods of basic model set, candidate model set, and valid model set are introduced in detail, respectively. Given space constraints, the following only takes the CA model framework as an example to introduce the corresponding steps.

2.2.1. Basic MSA

(1)
Initialize a rough grid.
Assuming that the acceleration (m/s 2 ) is a x , a y T , and it ranges from a x max , a y max T to a x max , a y max T . In order to cover the mode space to the maximum extent, the maximum values of acceleration are taken as the initial value of the grid. In addition, the model interval is set as λ k , the threshold parameter of the impossible model is t 1 , and t 2 is the threshold parameter of the model jump. In other words, if the model probability is less than t 1 , it is considered that the model does not match the current maneuvering state. If the model probability is greater than t 2 , it can be considered that the model has a high degree of matching with the current maneuvering state.
D 0 = a 0 L = a x max , 0 T , a 0 R = a x max , 0 T a 0 C = 0 , 0 T a 0 T = 0 , a y max T , a 0 D = 0 , a y max T
where a 0 C is the initial center model, a 0 L , a 0 R , a 0 T , and a 0 D represent the initial left, right, upper, and lower edge models, respectively.
(2)
The adaptive rules.
From k 1 to k, the following will introduce the adaptive rules for the grid.
Part 1: Adaptation of the grid center. Suppose that the model set parameters at k 1 are a k 1 L , a k 1 R , a k 1 C , a k 1 T , a k 1 D and their corresponding model probabilities are μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T and μ k 1 D , respectively. Then, the central value parameters of the grid are
a k C = μ k 1 L a k 1 L + μ k 1 R a k 1 R + μ k 1 C a k 1 C + μ k 1 T a k 1 T + μ k 1 D a k 1 D
Part 2: Adaptation of the grid edge. If, μ k 1 C = max μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T , μ k 1 D , it indicates that the central model is closest to the true system mode. At this point, if μ k 1 L , μ k 1 R , μ k 1 T , and μ k 1 D are all smaller than t 1 , it indicates that the left, right, upper, and lower edge models deviate from the current state of the system. Then, the model interval is reduced to λ k 1 / 2 . Otherwise, the model interval remains unchanged. The jump rule is formulated as follows:
a k L = a x ( k ) C λ k 1 L / 2 , 0 T , μ k 1 L < t 1 a x ( k ) C λ k 1 L , 0 T , o t h e r s a k R = a x ( k ) C + λ k 1 R / 2 , 0 T , μ k 1 R < t 1 a x ( k ) C + λ k 1 R , 0 T , o t h e r s a k T = 0 , a x ( k ) C + λ k 1 T / 2 T , μ k 1 T < t 1 0 , a y ( k ) C + λ k 1 T T , o t h e r s a k D = 0 , a x ( k ) C λ k 1 D / 2 T , μ k 1 D < t 1 0 , a y ( k ) C λ k 1 D T , o t h e r s
If μ k 1 L = max μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T , μ k 1 D , it indicates that the left edge model is closest to the true system mode, all the other three edge models keep the original model interval. At this time, if μ k 1 L is greater than t 2 , it can be considered that the difference between the left edge model and the system true mode is the smallest. In other words, the velocity in the X direction is constantly decreasing and the left edge model jumps 2 λ k 1 L to the left. Otherwise, let the left edge model jump λ k 1 L to the left. The jump rule is formulated as follows:
a k L = a x ( k ) C 2 λ k 1 L , 0 T , μ k 1 L > t 2 a x ( k ) C λ k 1 L , 0 T , o t h e r s a k R = a x ( k ) C + λ k 1 R , 0 T a k T = 0 , a y ( k ) C + λ k 1 T T a k D = 0 , a y ( k ) C λ k 1 D T
If μ k 1 R = max μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T , μ k 1 D , it indicates that the right edge model is closest to the true system mode, the other three edge models keep the original model interval. At this point, the jump rule of the right edge model is similar to that of the left edge model, which can be formulated as follows:
a k L = a x ( k ) C λ k 1 L , 0 T a k R = a x ( k ) C + 2 λ k 1 R , 0 T , μ k 1 R > t 2 a x ( k ) C + λ k 1 R , 0 T , o t h e r s a k T = 0 , a y ( k ) C + λ k 1 T T a k D = 0 , a y ( k ) C λ k 1 D T
If μ k 1 T = max μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T , μ k 1 D , it indicates that the upper edge model is closest to the true system mode, the other three edge models keep the original model interval. At this time, if μ k 1 T is greater than t 2 , it can be considered that the upper edge model has the highest matching degree with the true system mode. In other words, the velocity in the Y direction is constantly increasing, and the upper edge model jumps up 2 λ k 1 T . Otherwise, let the upper edge model jump up λ k 1 T . The jump rule can be formulated as follows:
a k L = a x ( k ) C λ k 1 L , 0 T a k R = a x ( k ) C + λ k 1 R , 0 T a k T = 0 , a y ( k ) C + 2 λ k 1 T T , μ k 1 T > t 2 0 , a y ( k ) C + λ k 1 T T , o t h e r s a k D = 0 , a x ( k ) C λ k 1 D T
If μ k 1 D = max μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T , μ k 1 D , it indicates that the lower edge model is closest to the true system mode. Its adaptive rules are similar to those of the upper edge model, which can be formulated as follows:
a k L = a x ( k ) C λ k 1 L , 0 T a k R = a x ( k ) C + λ k 1 R , 0 T a k T = 0 , a x ( k ) C + λ k 1 T T a k D = 0 , a y ( k ) C 2 λ k 1 D T , μ k 1 D > t 2 0 , a y ( k ) C λ k 1 D T , o t h e r s
where λ k 1 L = min a x k 1 C a x k 1 L , ϑ w , λ k 1 R = min a x k 1 R a x k 1 C , ϑ w , λ k 1 T = min a y k 1 T a y k 1 C , ϑ w , λ k 1 D = min a y k 1 C a y k 1 D , ϑ w , ϑ w refers to the minimum distance between each model, which not only ensures the comprehensiveness of the mode space but also maintains the independence between models.
The above describes how to use the grid to realize the adaptation of the basic model set, and the inspiration of its adaptive strategy mainly comes from literature [24]. Based on the adaptive mechanism, a basic model set M k b = a k L , a k R , a k C , a k T , a k D that adapt well to the maneuvering state can be generated in real-time. On the one hand, the distribution of these models is relatively concentrated. On the other hand, the number of models is relatively small.

2.2.2. Candidate MSA

The determination of the candidate model set can also utilize the AG technology, where the parameters of each model slide adaptively over a continuous region. Therefore, only a small number of models related to the system mode are required at any one time.
In this paper, five CA models are selected to form the candidate model set, and each model corresponds to different maneuvering levels. Assuming that the best model activated at k 1 is m k 1 b e s t , and the candidate model set at time k is M k c , then horizontal, as well as vertical adaptation are required. The candidate model set at k is
M k c = a x k 1 b e s t θ w , 0 T a x k 1 b e s t + θ w , 0 T a x k 1 b e s t , a y k 1 b e s t T 0 , a y k 1 b e s t + θ w T 0 , a y k 1 b e s t θ w T
where θ w is the distance between the models in the candidate model set, θ w 0.5 ϑ w . In this way, the candidate model set can be fine-tuned and closer to the true system mode, resulting in the powerful modification of the best model to the basic model set.

2.2.3. Valid MSA

The valid model set used for MM estimation at k can be obtained by M k = M k b m k b e s t . where, the best model, m k b e s t , is activated from the candidate model set, M k c . Additionally, the “best” here refers to the minimum difference between the model and the true system mode. In this paper, the complex KL criterion in the traditional BMA algorithm is discarded, and a simpler mechanism based on similarity distance is adopted to quantify this difference. The smaller the value of the Euclidean distance, the smaller is the difference [25]. The process can be formulated as follows:
m k b e s t = arg min m k j M k c s k m k j
Among them, since the true mode, s k , of the system is often unavailable, it can only be approximated with some prior information. In this paper, the expected model has been used to approximate the true system mode [24]. Based on the idea of finding the expected model in the EMA algorithm, it can be thought of as a weighted sum of probabilities, μ k 1 j .
s k = m k 1 j M k 1 m k 1 j μ k 1 j , k N
where μ k 1 j is the probability of the model m k 1 j , and M k 1 is the valid model set at k 1 .
In all, the complete pseudo codes for the above MSA algorithm can be summarized in Algorithm 1.
Algorithm 1: Pseudo codes for the proposed MSA algorithm.
➀ Basic model set adaption
D 0 = a 0 L = a x max , 0 T , a 0 R = a x max , 0 T a 0 C = 0 , 0 T a 0 T = 0 , a y max T , a 0 D = 0 , a y max T , where D 0 represents the initial model set.
Part 1: Adaptation of the grid center.
Suppose that the model set parameters at k 1 are a k 1 L , a k 1 R , a k 1 C , a k 1 T , a k 1 D and their
corresponding model probabilities are μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T and μ k 1 D , respectively. Then,
a k C = μ k 1 L a k 1 L + μ k 1 R a k 1 R + μ k 1 C a k 1 C + μ k 1 T a k 1 T + μ k 1 D a k 1 D .
Part 2: Adaptation of the grid edge.
When μ k 1 C = max μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T , μ k 1 D , then execute Equation (11);
When μ k 1 L = max μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T , μ k 1 D , then execute Equation (12);
When μ k 1 R = max μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T , μ k 1 D , then execute Equation (13);
When μ k 1 T = max μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T , μ k 1 D , then execute Equation (14);
When μ k 1 D = max μ k 1 L , μ k 1 R , μ k 1 C , μ k 1 T , μ k 1 D , then execute Equation (15).
➁ Candidate model set adaption
The candidate mode set can be formulated as Equation (16), where a x k 1 b e s t , a y k 1 b e s t T represents
the best model at the above moment.
➂ Determination of the valid model set
The approximate true system mode is s k = m k 1 j M k 1 m k 1 j μ k 1 j , k N ;
The activated best model can be obtained by m k b e s t = arg min m k j M k c s k m k j ;
The valid model set is M k = M k b m k b e s t .
Then, k = k + 1 , return to step ➀ and proceed to the next iteration.
To visualize the whole model set design process, the topology of the AG-BMA algorithm is shown in Figure 4. In the simulated well-known examples of [12,26,27], the maximum acceleration in any coordinate direction is assumed to be 4 g (g = 10 m/s 2 ). Thus, the mode space can be defined as follows:
A c = a x , a y : a x 2 + a y 2 40
Moreover, a more intuitive explanation of the position is given. Here, we use the circular areas (the center of the circle is 0 , 0 and the radius is 40 (m)) to represent the entire mode space A c . The coordinates a x , a y of any point in the region can represent a kind of motion mode, where a x and a y represent the acceleration (m/s 2 ) of the target in the X-axis and Y-axis direction, respectively. For example, the coordinate point 40 , 0 indicates a motion that the constant deceleration with an acceleration of 40 (m/s 2 ) in the X-axis direction and the constant velocity in the Y-axis direction. Figure 4a,b represent the positions of the candidate model sets at k 1 and k in the entire mode space, in which the symbol ‘∗’ denotes the true system mode and these black circles form the candidate model set. Figure 4c,d represent the positions of the basic model set at k 1 and k, where these white circles represent four edge models and one center model, respectively. Figure 4e,f represent the valid model set, in which the black circle denotes the best model augmented from the candidate model set. The models in the valid model set will be used to perform the MM estimation in Section 2.3. It is not difficult to see that the candidate model set and the basic model set have formed a double-layer coverage near the true system mode. The valid model set at the current moment is only composed of the best model and the basic model set, thus the true mode ‘∗’ of the system can be maximally approximated with fewer models.

2.3. MM Estimation

Based on the valid model set at k, MM M k , M k 1 is run to obtain the filtering results, x ^ k j , p k j , of each model and calculate the model probability, μ k j . The multi-model estimator uses the simplest reinitialization method. First of all, the last state estimation, x ^ k 1 | k 1 , and its covariance matrix, P k 1 | k 1 , are taken as the common initial inputs. Then, each model performs its own state estimation using the basic filtering algorithm and calculates the probability of each model. Finally, the state estimation result and its error covariance matrix are obtained by weighted summation. Its algorithm structure is shown in Figure 5.
For j = 1 , 2 , , r , the loop is as follows:

2.3.1. Reinitialization of the Model Conditions

Since each model filter can be a valid system model filter, the initial conditions of each filter are based on the synthesis of the filtering results of the previous model. Then, the input of the j t h model filter after reinitialization [28] is
x ^ ^ k 1 | k 1 ( j ) = d e f E x k 1 | m k ( j ) , Z k 1 = i = 1 r x ^ k 1 | k 1 ( j ) μ k 1 | k 1 ( i , j ) , j = 1 , 2 , r
P ^ k 1 | k 1 ( j ) = i = 1 r P ^ k 1 | k 1 ( i ) + x ^ k 1 | k 1 ( i ) x ^ ^ k 1 | k 1 ( j ) x ^ k 1 | k 1 ( i ) x ^ ^ k 1 | k 1 ( j ) T μ k 1 | k 1 ( i , j ) j = 1 , 2 , r
μ k 1 | k 1 ( i , j ) = P m k 1 i | m k j , Z k 1 = 1 c ¯ j π i j μ k 1 ( i ) , i , j = 1 , 2 , , r
Assuming that the matching model at k 1 is m k 1 i and the matching model at k is m k j , Equation (22) represents the mixed probability under the condition of Z k 1 . Where c ¯ j = i = 1 r π i j μ k 1 i , μ k 1 ( i ) is the matching probability of the model i and the true system mode at k 1 .

2.3.2. Model Conditional Filtering

Given the initial state and the covariance matrix, the next step is to update the state estimation according to the new measurements z k obtained at k. Assume that the model m k 1 i is adopted at k 1 and the model m k j is adopted at the moment k. Further, taking x ^ ^ k 1 | k 1 ( j ) from Equation (20) and P ^ k 1 | k 1 ( j ) from Equation (21) as the initial values, the state estimation, x ^ k | k ( i , j ) , and its covariance matrix, P k | k ( i , j ) , are calculated by using the models matching with the true system mode and their corresponding filtering equation.
(1) State prediction:
For j = 1 , 2 , , r , calculate separately
x ^ k | k 1 j = F k 1 j x ^ ^ k 1 | k 1 ( j ) + Γ k 1 j w ¯ k 1 j
P k | k 1 j = F k 1 j P ^ k 1 | k 1 ( j ) F k 1 j T + Γ k 1 j Q k 1 j Γ k 1 j T
(2) Measurement prediction residual, z ˜ k j , and its covariance calculation:
For j = 1 , 2 , , r ,
z ˜ k j = z k H k j x ^ k | k 1 j υ ¯ k j
S k j = H k j P k | k 1 j H k j T + R k j
The likelihood function matching m k j is calculated as follows:
k j = p z k | m k j , Z k 1 p z k | m k j , x ^ ^ k 1 | k 1 j , S k j P ^ k 1 | k 1 j
Under the Gauss hypothesis, the likelihood function can be calculated as follows:
k j = p z ˜ k j | m k j , Z k 1 = G a u s s 2 π S k j 1 2 exp 1 2 z ˜ k j T S k j 1 z ˜ k j
(3) Update of the filtering results:
For j = 1 , 2 , r , the filter gain array, the updated status estimation result, and its error covariance are calculated as follows:
K k j = P k | k 1 j H k j T ( S k j ) 1
x ^ k | k j = x ^ k | k 1 j + K k j z ˜ k j
P k | k j = P k | k 1 j K k j S k j K k j T
where K k j is the Kalman gain.
(4) Update of the model probabilities:
For j = 1 , 2 , r , the model probability is calculated as follows:
μ k j = P m k j | Z k = 1 c k j c ¯ j , j = 1 , 2 , r
where c ¯ j = i = 1 r π i j μ k 1 i , c = i = 1 r k i c ¯ i , k j represents the likelihood function matching model m k j .
(5) Estimation of fusion:
While synthesizing the estimation results of each model, the updated model probability is used as the weight coefficient. Then, the overall estimation and its error covariance matrix at is formulated as follows:
x ^ k | k = j = 1 n x ^ k | k j μ k j
P k | k = j = 1 r P k | k j + x ^ k | k x ^ k | k j x ^ k | k x ^ k | k j T μ k j
where x ^ k | k j is the updated state estimation result of the j t h filter.
As for the moment at k, the VSMM algorithm based on AG-BMA method completes an iteration process. Additionally, then, k = k + 1 , return to step (1) and proceed to the next iteration.

3. Simulation Results

In order to prove that the proposed method is effective in improving the estimation performance of the kinematic state and the extension state, an extended target with star-convex shape, as shown in Figure 6, was used in this paper. Assuming that the measurement sources are all distributed on the target surface, and the radar sensor with high resolution capability and detection accuracy obtains the information. The tracking object of this paper is a single maneuvering extended target, which does not involve clutter and missed detection.
The performance of the proposed VSMM algorithm based on AG-BMA was investigated over a deterministic maneuver scenario and a random scenario. The two scenarios enable the evaluation of the kinematic state and the shape estimation performance of the algorithm, which can be judged based on the root-mean-square error and Hausdorff distance, respectively. In addition, the proposed algorithm has been compared with the IMM algorithm and two other VSMM algorithms to illustrate its effectiveness.
Here, we need to emphasize the formula of Hausdorff distance d H S ¯ C x k , S ¯ C ^ x ^ k between the estimated shape S ¯ C x k and the real shape S ¯ C ^ x ^ k , where C and C ^ represent the estimated target and the real target, respectively. The radial function that used to describe star-convex is characterized by its extension parameter vector x k e , i.e.,
r x k e , ϕ = R ϕ x k e
where ϕ is a continuous set and can be replaced by a discrete set with the idea of angular uniform discretization.
r ˙ x k e , ϕ = R ϕ x k e , ϕ i 2 π N s , i = 1 , 2 , , N s
where N s is the number of samples. Accordingly, S ¯ C x k and S ¯ C ^ x ^ k can be expressed as
S ˙ C x k = s ^ k r ˙ x k e , ϕ e ϕ + x k c , ϕ i 2 π N s , i = 1 , 2 , , N s
S ˙ C ^ x ^ k = s ^ k r ˙ x ^ k e , ϕ e ϕ + x ^ k c , ϕ i 2 π N s , i = 1 , 2 , , N s
Then, the formula of Hausdorff distance is [29]
d H S ˙ C x k , S ˙ C ^ x ^ k = max d S ˙ C x k , S ˙ C ^ x ^ k , d S ˙ C ^ x ^ k , S ˙ C x k
d S ˙ C x k , S ˙ C ^ x ^ k = max a S ˙ C x k d E a , S ˙ C ^ x ^ k
d S ˙ C ^ x ^ k , S ˙ C x k = max a ^ S ˙ C ^ x ^ k d E a , S ˙ C x k
where a and a ^ represent a point of the estimated shape and the real shape, respectively. d E a , S ˙ C ^ x ^ k denotes the Euclidean distance between point a and each point in S ˙ C ^ x ^ k , d E a , S ˙ C x k denotes the Euclidean distance between point a ^ and each point in S ˙ C x k .

3.1. Deterministic Scenario

In order to demonstrate the general applicability and universality of the proposed algorithm, a deterministic scenario (DS) was introduced. The simulation was carried out under the framework of the CT and CV models, in which the structure and parameters were all changed.
Table 1 lists the kinematic state of the target, namely the serpentine turning process. The parameter settings of the CT models with different turning rates were as follows:
X = x , v x , y , v y , w T X 0 = 1000 m , 150 m / s , 1000 m , 100 m / s , 0 rad / s T H k j = 1 0 0 0 0 1 0 0 F k j = 1 sin w j T w j 0 cos w j T w 0 0 cos w j T 0 sin w j T 0 0 1 cos w j T w j 1 sin w j T w j 0 0 0 sin w j T 0 cos w j T 0 0 0 0 0 1 w k N 0 , 0.01 2 , υ k N 0 , 5
When using the proposed AG-BMA-based VSMM algorithm for state estimation, the adaptive process of the model set was different from that of the CA model framework introduced in Section 2.2. Suppose that the turning rate of the target is w and its range of variation is w max , w max , then the boundary values of the turning rate are taken as the initial value of the grid, and only horizontal self-adaptation is required.
When the basic model set was used, the number of models was set to 3. The model set at the initial time was set to D 0 = w 0 L = w max , w 0 C = 0 , w 0 R = w max , where w max = 15 rad / s . The adaptive steps of the grid center and the grid distance were similar to those of the constant acceleration scenarios, and, thus, a detailed description has not been given here. Eventually, an adaptive model set, M k b = w k L , w k C , w k R , was obtained, which was generated in real-time according to the maneuvering state of the target.
Similarly, three CT models were selected to form the candidate model set, and the turning rate of each model corresponded to different maneuvering levels. If the best model activated at k 1 is m k 1 b e s t , then the set of candidate models at k can be defined as M k c = w k 1 b e s t θ w , w k 1 b e s t , w k 1 b e s t + θ w . In this paper, θ w was taken to be 5 rad/s in order to achieve fine tuning.
Figure 7 depicts the comparative results between the three algorithms, namely, IMM, BMA, and AG-BMA, over 100 Monte Carlo runs. Hausdorff distance and root-mean-square error are used to evaluate the performance of the algorithm. The smaller their values, the better is the tracking performance of the corresponding algorithm.
The root-mean-square error values (unit: m, m/s) and the computational load of the three algorithms are given in Table 2 and Table 3, respectively.
Figure 7a,b show the estimation performance of the extension state. From the results, it can be inferred that BMA provides the poor shape estimation performance. Furthermore, the tracking performance of the three algorithms is largely consistent in the non-maneuvering stage. However, in the maneuvering stage, the AG-BMA algorithm exhibits superior performance in shape estimation. As shown in Figure 7a, the shape estimated by the AG-BMA algorithm is closer to the true shape of the target.
Figure 7c,d show the root-mean-square error comparison of the centroid position and velocity resulting from the application of the three algorithms. From the results, it can be seen that the proposed AG-BMA algorithm can better improve the estimation performance of the kinematic state. This is mainly because the AG-BMA algorithm can automatically adjust the parameters of each model in the basic and candidate model set depending on the maneuvering condition in order to efficiently match the true system mode. However, the results show that even without activating a new model, the IMM algorithm still has a smaller estimation error than BMA sometimes. The reasons for this result can be summarized as follows: On the one hand, the BMA algorithm can achieve optimal estimation only when the true system mode is in the whole mode space. On the other hand, this could be because the model contained in the IMM algorithm happens to be able to describe part of the true system mode.

3.2. Random Scenario

To provide a performance comparison over an ensemble of maneuver trajectories that is as fair as possible, the algorithms were tested on the random scenario (RS) that has been developed in previous studies [26,30]. In the entire simulation process, the target performed random acceleration motion, and thus it can be regarded as a continuous strong maneuvering process. Here, the acceleration vector is a two-dimensional semi-Markov process and it jumps from one state of acceleration of magnitude α and phase θ to another after a random time interval. Obviously, such strong maneuvering scenarios put higher demands on the performance of the MM estimators.
The specific model assumptions in this scenario were as follows: The program executed 90 simulation steps each time and the acceleration changed 10 times in the entire process. Then, the residence time, τ k , of the state a = a k was a random number up to 90 and satisfied the condition N = 1 10 τ k = 90 . The phase angle, θ k , obeyed a Gaussian distribution with mean, θ ¯ , and variance, σ θ 2 . In this simulation scenario, the parameters were set as follows:
X = x , v x , a x , y , v y , a y T X 0 = 1000 m , 100 m / s , 10 m / s 2 , 1000 m , 100 m / s , 10 m / s 2 T H k j = 1 0 0 0 0 0 0 0 0 1 0 0 F k j = 1 T T 2 2 0 0 0 0 1 T 0 0 0 0 0 1 0 0 0 0 0 0 1 T T 2 2 0 0 0 0 0 0 0 0 1 0 T 1 w k N 0 , 0.01 2 , υ k N 0 , 5
Figure 8 depicts the comparative results of the four algorithms, IMM7, IMM13, BMA, and AG-BMA, over 100 Monte Carlo runs for an random scenario. The Hausdorff distance and the root-mean-square error were used as the evaluation metrics for accessing the estimation results of the extension state and kinematic state, respectively. The smaller their values, the better is the tracking performance of the corresponding algorithm.
The root-mean-square error and the computational load corresponding to the four algorithms are given in Table 4 and Table 5, respectively.
Figure 8a,b show the performance of four algorithms for shape estimation. The smaller the Hausdorff distance is, the better the shape estimation is. Obviously, the extension state estimated by AG-BMA algorithm is closer to the true shape of the target.
From the results presented in Table 4 and Table 5, it can be seen that the original BMA algorithm exhibits the poor estimation accuracy and also has a large computational burden. The reasons for this result can be summed up as follows: On one hand, the KL information needs to be frequently calculated each time the best model is activated. On the other hand, it is due to the uncertainty caused by its fixed and prior known model set. Although the IMM7 algorithm requires the shortest computational time, its estimation error is bigger than IMM13 and AG-BMA. IMM13 can cover the entire model space to a large extent due to a large number of models, but it is still less effective than AG-BMA, and it also requires the largest amount of calculation. Undoubtedly, the key to improving the tracking accuracy is to ensure that the dynamic model used is as close as possible to the true system mode of the target. Fortunately, the proposed AG-BMA algorithm can construct an adaptive model set online based on the information of the previous moment. Thus, the AG-BMA algorithm can obtain the best estimation performance of kinematic state among the four algorithms.

4. Discussion

From the performance comparison in Section 4, we can draw a conclusion that the proposed AG-BMA algorithm can obtain a better estimation of kinematic state while requiring a slightly longer calculation time. This is mainly because the proposed algorithm can adjust the model parameters online according to the information of the previous moment, which helps the dynamic model approximates the true system mode to a greater extent.
In addition, although we only try to obtain a more matched dynamic model, the proposed method still shows better shape estimation performance. In reality, this can be attributed to the inherent coupling relationship between shape and motion. For example, the target might always be aligned with its motion direction. It is obvious that a more detailed shape model improves the estimation quality for the kinematic state as it accurately captures the target. Meanwhile, an accurate estimation of the target position will also help to describe the shape of the target. In other words, the estimates for the extension parameters, as well as the kinematic parameters can benefit from this coupling. This well explains why a more suitable dynamic model can also describe the shape more accurately instead of just improving the estimation accuracy of kinematic state.

5. Conclusions

To improve the tracking accuracy of the maneuvering extended target, a new model set design method, namely, AG-BMA, has been proposed in this work. This method utilizes an adaptive grid to achieve a dual adaptation of the basic and candidate model sets so that it can approximate the true system mode to a large extent with a small number of models. Furthermore, it can effectively overcome the model mismatch caused by strong maneuvering. In this paper, the proposed algorithm was simulated in a deterministic scenario and a random scenario to prove the effectiveness of the method. From the simulation results, it can be concluded that the proposed AG-BMA algorithm can effectively improve the estimation accuracy of the extended shape and centroid motion.
How to designing an effective model set for multiple maneuvering extended targets tracking is still challenging. The proposed model set designing method (i.e., AG-BMA) has room for further development and improvement. As our research directions to next step, we will consider modifying AG-BMA and apply it to tracking of multiple maneuvering extended targets. In this context, more work remains to be completed.

Author Contributions

L.S. and J.Z. conceived of the idea and developed the proposed approaches. H.Y. advised the research and helped edit the paper. Z.F. and Z.H. improved the quality of the manuscript and the completed revision. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by the National Thirteen-Five Equipment Pre-Research Foundation of China (No. 61403120207), the Aeronautical Science Foundation of China (No. 20185142003), Natural Science Foundation of Henan Province, China (No. 222300420433), the Science and Technology Innovative Talents in Universities of Henan Province (No. 21HASTIT030), Young Backbone Teachers in Universities of Henan Province (No. 2020GGJS073).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, L.; Mao, D.; Niu, J.; Wu, Q.; Ji, Y. Continuous tracking of targets for stereoscopic HFSWR based on IMM filtering combined with ELM. Remote Sens. 2020, 12, 272. [Google Scholar] [CrossRef] [Green Version]
  2. Sun, L.; Yu, H.; Lan, J.; Fu, Z.; He, Z.; Pu, J. Tracking of Multiple Maneuvering Random Hypersurface Extended Objects Using High Resolution Sensors. Remote Sens. 2021, 13, 2963. [Google Scholar] [CrossRef]
  3. Wei, H.; Huang, Y.; Hu, F.; Zhao, B.; Guo, Z.; Zhang, R. Motion Estimation Using Region-Level Segmentation and Extended Kalman Filter for Autonomous Driving. Remote Sens. 2021, 13, 1828. [Google Scholar] [CrossRef]
  4. Bensheng, Q.; Haiyi, L.; Youquan, C.; Diancheng, Z. An adaptive parameter model for maneuvering target tracking. In Proceedings of the ICSP’98—1998 Fourth International Conference on Signal Processing (Cat. No. 98TH8344), Beijing, China, 12–16 October 1998; pp. 465–469. [Google Scholar]
  5. Yadav, A.; Awasthi, P.; Naik, N.; Ananthasayanam, M. A constant gain Kalman filter approach to track maneuvering targets. In Proceedings of the 2013 IEEE International Conference on Control Applications (CCA), Hyderabad, India, 28–30 August 2013; pp. 562–567. [Google Scholar]
  6. Moose, R.L.; Vanlandingham, H.F.; McCabe, D. Modeling and estimation for tracking maneuvering targets. IEEE Trans. Aerosp. Electron. Syst. 1979, AES-15, 448–456. [Google Scholar] [CrossRef]
  7. Yang, L.N.; Yuan, Z.; Yan, B.D.; Xu, L. Application of IMM algorithm based on model CS in target tracking. Jisuanji Gongcheng Yingyong 2010, 46, 230–232. [Google Scholar]
  8. Blom, H.A.P. An efficient filter for abruptly changing systems. In Proceedings of the 23rd IEEE Conference on Decision and Control, Las Vegas, NV, USA, 12–14 December 1984; pp. 656–658. [Google Scholar]
  9. Li, X.R. Multiple-model estimation with variable structure. II. Model-set adaptation. IEEE Trans. Autom. Control 2000, 45, 2047–2060. [Google Scholar]
  10. Li, X.R.; Bar-Shalom, Y. Multiple-model estimation with variable structure. IEEE Trans. Autom. Control 1996, 41, 478–493. [Google Scholar]
  11. Dong, P.; Jing, Z.; Li, M.; Pan, H. The variable structure multiple model GM-PHD filter based on likely-model set algorithm. In Proceedings of the 2016 19th International Conference on Information Fusion (FUSION), Heidelberg, Germany, 5–8 July 2016; pp. 2289–2295. [Google Scholar]
  12. Li, X.R.; Jilkov, V.P.; Ru, J. Multiple-model estimation with variable structure-part VI: Expected-mode augmentation. IEEE Trans. Aerosp. Electron. Syst. 2005, 41, 853–867. [Google Scholar]
  13. Lan, J.; Li, X.R. Equivalent-model augmentation for variable-structure multiple-model estimation. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 2615–2630. [Google Scholar] [CrossRef]
  14. Mao, Y.; Zhou, X.; Zhang, J. Passive radar tracking of a maneuvering target using variable structure multiple-model algorithm. In Proceedings of the Fifth International Conference on Machine Vision (ICMV): Algorithms, Pattern Recognition, and Basic Technologies, Wuhan, China, 20–21 October 2012. [Google Scholar]
  15. Xu, L.; Li, X.R.; Duan, Z. Hybrid grid multiple-model estimation with application to maneuvering target tracking. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 122–136. [Google Scholar] [CrossRef]
  16. Zhang, Y.; Guo, C.; Hu, H.; Liu, S.; Chu, J. An algorithm of the adaptive grid and fuzzy interacting multiple model. J. Mar. Sci. Appl. 2014, 13, 340–345. [Google Scholar] [CrossRef]
  17. Wu, P.; Li, X.; Zhang, L.; Bo, Y. Tracking algorithm with radar and infrared sensors using a novel adaptive grid interacting multiple model. IET Sci. Meas. Technol. 2014, 8, 270–276. [Google Scholar] [CrossRef]
  18. Lan, J.; Li, X.R.; Mu, C. Best model augmentation for variable-structure multiple-model estimation. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 2008–2025. [Google Scholar] [CrossRef]
  19. Kullback, S.; Leibler, R.A. On information and sufficiency. Ann. Math. Stat. 1951, 22, 79–86. [Google Scholar] [CrossRef]
  20. Baum, M.; Hanebeck, U.D. Shape tracking of extended objects and group targets with star-convex RHMs. In Proceedings of the 14th International Conference on Information Fusion, Chicago, IL, USA, 5–8 July 2011; pp. 1–8. [Google Scholar]
  21. Baum, M.; Hanebeck, U.D. Extended object tracking with random hypersurface models. IEEE Trans. Aerosp. Electron. Syst. 2014, 50, 149–159. [Google Scholar] [CrossRef] [Green Version]
  22. Chen, H.; Du, J. An adaptive contour algorithm for irregular shape extended target. In Proceedings of the 2017 IEEE International Conference on Mechatronics and Automation (ICMA), Takamatsu, Japan, 6–9 August 2017; pp. 734–739. [Google Scholar]
  23. Sun, L.; Yu, H.; Fu, Z.; He, Z.; Zou, J. Modeling and Tracking of Maneuvering Extended Object With Random Hypersurface. IEEE Sens. J. 2021, 21, 20552–20562. [Google Scholar] [CrossRef]
  24. Yimao, L. Research on Variable Structure Multi-Model Maneuvering Target Tracking Algorithm. Ph.D. Thesis, Xi Dian University, Xi’an, China, 2014. [Google Scholar]
  25. Liu, M.; Chen, H.; Qiu, Z.; Ren, X. Moving Target Location Method Based on Euclidean Distance and Camshift Algorithm. In Proceedings of the 2018 Eighth International Conference on Instrumentation & Measurement, Computer, Communication and Control (IMCCC), Harbin, China, 19–21 July 2018; pp. 558–563. [Google Scholar]
  26. Li, X.R.; Zhang, Y. Multiple-model estimation with variable structure. V. Likely-model set algorithm. IEEE Trans. Aerosp. Electron. Syst. 2000, 36, 448–466. [Google Scholar]
  27. Li, X.R.; Jilkov, V.P. Expected-mode augmentation for multiple-model estimation. In Proceedings of the 2001 International Conference on Information Fusion, Montreal, QC, Canada, 7–10 August 2001; pp. 7–10. [Google Scholar]
  28. Seah, C.E.; Hwang, I. State estimation for stochastic linear hybrid systems with continuous-state-dependent transitions: An IMM approach. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 376–392. [Google Scholar] [CrossRef]
  29. Sun, L.; Zhang, S.; Ji, B.; Pu, J. Shape Estimation Evaluation of Extended Objects Based on Modified Hausdorff Distance. Acta Opt. Sin. 2017, 37, 0728003. [Google Scholar]
  30. Li, X.R.; Jilkov, V.P.; Ru, J.; Bashi, A. Expected-mode augmentation algorithms for variable-structure multiple-model estimation. IFAC Proc. Vol. 2002, 35, 175–180. [Google Scholar] [CrossRef]
Figure 1. Representation of a star-convex shape using a radial function. (a) Star-convex extended target. (b) Radial function.
Figure 1. Representation of a star-convex shape using a radial function. (a) Star-convex extended target. (b) Radial function.
Remotesensing 14 01613 g001
Figure 3. Recursive flow of model set changes in AG-BMA algorithm.
Figure 3. Recursive flow of model set changes in AG-BMA algorithm.
Remotesensing 14 01613 g003
Figure 4. Topology of the model set based on the AG-BMA algorithm.
Figure 4. Topology of the model set based on the AG-BMA algorithm.
Remotesensing 14 01613 g004
Figure 5. Structure of the MM algorithm.
Figure 5. Structure of the MM algorithm.
Remotesensing 14 01613 g005
Figure 6. Star-convex extended target and its measurements.
Figure 6. Star-convex extended target and its measurements.
Remotesensing 14 01613 g006
Figure 7. Performance comparison of the three algorithms. (a) Trajectory of the target. (b) Hausdorff distance. (c) RMSE of position. (d) RMSE of velocity.
Figure 7. Performance comparison of the three algorithms. (a) Trajectory of the target. (b) Hausdorff distance. (c) RMSE of position. (d) RMSE of velocity.
Remotesensing 14 01613 g007
Figure 8. Performance comparison of the four algorithms. (a) Trajectory of the target (b) Hausdorff distance. (c) RMSE of acceleration. (d) RMSE of position. (e) RMSE of velocity.
Figure 8. Performance comparison of the four algorithms. (a) Trajectory of the target (b) Hausdorff distance. (c) RMSE of acceleration. (d) RMSE of position. (e) RMSE of velocity.
Remotesensing 14 01613 g008
Table 1. Maneuverable process of the target in DS.
Table 1. Maneuverable process of the target in DS.
Step[0,10)[10,30)[30,50)[50,70)[70,80)[80,100)[100,110]
w ( rad / s ) 0 8 080 8 0
Table 2. Average root-mean-square error corresponding to the three algorithms.
Table 2. Average root-mean-square error corresponding to the three algorithms.
RMSEIMMBMAAG-BMA
Position (m)6.1455.0324.294
Velocity (m/s)0.4020.5360.133
Table 3. Computational time taken by the three algorithms in DS.
Table 3. Computational time taken by the three algorithms in DS.
AlgorithmsIMMBMAAG-BMA
Time (s)1.9433.4143.400
Table 4. Average root-mean-square error corresponding to the four algorithms.
Table 4. Average root-mean-square error corresponding to the four algorithms.
RMSEIMM7IMM13BMAAG-BMA
Acceleration (m/s 2 )14.50613.44214.6657.404
Position (m)40.18135.94448.94020.226
Velocity (m/s)23.71023.69636.36716.576
Table 5. Computational time taken by the four algorithms in RS.
Table 5. Computational time taken by the four algorithms in RS.
AlgorithmsIMM7IMM13BMAAG-BMA
Time (s)3.0175.5444.1593.998
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, L.; Zhang, J.; Yu, H.; Fu, Z.; He, Z. Tracking of Maneuvering Extended Target Using Modified Variable Structure Multiple-Model Based on Adaptive Grid Best Model Augmentation. Remote Sens. 2022, 14, 1613. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14071613

AMA Style

Sun L, Zhang J, Yu H, Fu Z, He Z. Tracking of Maneuvering Extended Target Using Modified Variable Structure Multiple-Model Based on Adaptive Grid Best Model Augmentation. Remote Sensing. 2022; 14(7):1613. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14071613

Chicago/Turabian Style

Sun, Lifan, Jinjin Zhang, Haofang Yu, Zhumu Fu, and Zishu He. 2022. "Tracking of Maneuvering Extended Target Using Modified Variable Structure Multiple-Model Based on Adaptive Grid Best Model Augmentation" Remote Sensing 14, no. 7: 1613. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14071613

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