Multiple fluctuating targets track-before-detect using multi-Bernoulli filter in radar sensor

This paper addresses the detection and tracking of multiple fluctuating targets for a track-before-detect algorithm based on the Multi-Bernoulli (MB-TBD) filter in surveillance radar systems. MB-TBD usually considers target amplitude information and ignores the fact that radar measurements are complex-valued. In this paper, we first propose to utilize phase information to improve the discrimination of targets from noise. More precisely, complex likelihood ratios are used instead of squared modulus measurements likelihood ratios for fluctuations of types Swerling 0, 1, 3. Secondly, the traditional MB-TBD filter cannot solve the problem of coexistence between targets with stronger amplitude and weaker amplitude when multiple fluctuating targets are moving. To address this limitation, an adaptive birth distribution based on joint successive target cancellation and measurement likelihood ratio driven is proposed. Moreover, in order to reduce computational complexity, the Bernoulli components of the same targets are merged after the MB-TBD updating. Finally, the proposed algorithm is implemented using Sequential Monte Carlo technology. The simulation results show that in challenging scenarios, the performance of the improved algorithm is better than the traditional algorithm, and it has a good application prospect.

high computational complexity. A well-established real-time TBD approach is based on recursive implementations [8][9][10] of particle filtering (PF) algorithm. Nonetheless, the PF algorithm cannot solve the problem well that the number of targets in MTT is unknown and time-varying. Another class of real-time TBD tracking is one of the methods that use random finite set (RFS). The MTT algorithm based on RFS theory solves the problems of the timevarying number of targets, data associations uncertainty, and detection uncertainty. The multi-Bernoulli (MB) filter [11][12][13] was another RFS-based MTT method proposed after the probability hypothesis density (PHD) filter [14][15][16][17][18] and the cardinalized PHD (CPHD) [19][20][21] filter. Compared to the two algorithms mentioned above, the MB filter has an advantage in terms of accuracy and computational complexity. It is worth noting that the above filters cannot provide identity information for multi-object estimation, and the existence of this problem drives the development of label RFS [22][23][24].
The key to the realization of MB-TBD is the calculation of the measurement likelihood ratio function (MLRF) based on the hidden state condition. It is worth noting that most MB-TBD algorithms in radar system choose the MLRF that only considers the amplitude information, i.e., square modulus measurement likelihood ratio (SLR) [25,26]. For SLR that discard phase information, there are two ways to deal with the amplitude fluctuations. The first one is the overall likelihood ratio of the marginalized amplitude fluctuation density [27]. However, in practice this treatment is difficult to obtain an exact solution and can only be achieved using numerical approximation. The second way consists in marginalizing the likelihood in each cell independently [28]. The advantage of this heuristic second solution is that a closed form can be obtained. Unlike infrared and optical sensors, which receive real-valued measurements that contain only amplitude information, radar measurements are complex-valued, i.e., phase information in addition to amplitude information. The literature [29] shown that the phase loss leads to a filter sensitivity decrease.
In order to improve the recognition of targets from noise, Davey et al. [29] proposed a MLRF method containing phase information, called complex likelihood ratio (CLR). In addition, due to the change of target relative radar view and the influence of other unknown factors causing the target RCS to change, making the target amplitude information is fluctuating. Typical radar fluctuating targets include Swerling 0, 1, 3 [30]. The literatures [30,31] introduce CLR into PF-TBD and DP-TBD, respectively, and consider a single target with fluctuation types of Swerling 0, 1 and 3. Single-target tracking of Swerling 0, 1, 3 target types in the Rayleigh sea clutter model was realized by the use of finite difference with continuous-discrete filtering in [32]. It is worth noting that the measurements used in [25] are square measurements.
In the detection and tracking of fluctuating targets, the existing research mainly introduces CLR into the filter to solve the problem of the single fluctuating target. Less research has been done for multiple fluctuating targets with unknown and time-varying target numbers. The problem of the co-existence of stronger and weaker amplitude targets can occur during multiple fluctuating targets movements. The traditional MB-TBD filter ignores targets with weak amplitudes, which leads to errors in the estimation of the target number. What is more serious is that the birth prior information of the target of traditional MB-TBD is known, which is inconsistent with most practical situations. Although there are some measurement-based adaptive birth algorithms, such as the literature [25,33] and so on. However, during the fluctuation process of the target amplitude, the amplitude information may be completely obscured in the noise for several consecutive frames. The maximum echo information received by the receiver may not be the surveillance target. When the fluctuating target reappears, the traditional MB-TBD filter will fail because the target prior information is unknown at this time.
To address the above issues, this paper extends the radar CLR method to MB-TBD, and embeds the assumption that target amplitude fluctuation follows the Swerling 0, 1, 3 model. Swerling statistical models belong to the RCS fluctuation statistical models, which are universal and include more radar target types. For example, Swerling 1 typical target such as small jet and Swerling 3 typical target such as helicopter. Specific contributions of this paper are given as follows.
Firstly, the expression of SLR with amplitude fluctuations of type Swerling 0, 1, 3 is modified, the expression of CLR considering phase information is given, and CLR is used as the likelihood ratio function of MB-TBD instead of SLR. The simulation results show that the MB-TBD algorithm-based CLR under amplitude fluctuation has better detection tracking performance than SLR, and it can reduce the complexity of the algorithm and run faster.
Secondly, in order to solve the problem that prior knowledge for target births is unknown, joint measurement likelihood ratio driven and successive-target-cancellation (STC)-based adaptive birth distribution for MB-TBD are proposed (LABer-STC-TBD), which is able to pick out the measurements that are likely to arise from real targets, so that the filter iteration can be maintained even when the target echo strength is covered below the noise.
Finally, after the update process, to account for the non-overlapping assumption, the overlapping estimates are merged. The Bernoulli components of the same target are merged at the end of the MB-TBD update, i.e., the trajectories are managed to estimate the correct number of targets after the MB-TBD filter is updated. The proposed trajectory management algorithm is more accurate and has lesser target estimation error.
The rest of the paper is organized as follows. Section 2 presents the multi-target state model under RFS, TBD measurement model, and the multi-target tracking problem under RFS. Section 3 gives the formulas for the CLR and the SLR under the amplitude fluctuation type Swerling 0, 1, 3. Section 4 presents the MB-TBD filter, the adaptive birth distribution and the trajectory merging model under the amplitude fluctuation type, and gives the SMC implementation. The experimental results are followed in Sect. 5. Results and discussion are given in Sect. 6. Finally, we conclude in Sect. 7.

Multi-targets dynamic model
At time k, the kinematic state of the i-th individual target is represented by where x k,i , y k,i and ẋ k,i ,ẏ k,i , respectively, denote the ith target position and velocity, T denotes the matrix transpose. For the multi-target tracking problem, since the number of targets is an unknown and random variable, the ranking of target states has no physical significance, so the state of the multi-target X k can be expressed as a random finite set (RFS) [11], which is defined as follows where F(X ) denote the collection of all subsets of state space X , and N k denote the number of targets at time k.
The dynamics model of multi-target state can be expressed as: where X k is the multi-target state set at time k, S k|k−1 (ζ ) is surviving RFS of target at time k that evolves from a target with previous state ζ , B k|k−1 (ζ ) is the spawning RFS of target at time k, Ŵ k is the RFS of spontaneous births at time k. In addition, it is assumed that each RFS in (2) is independent of each other.

Measurement model
In this paper we considered radar sensor is an active monopulse radar, with a linear phased array antenna. In the TBD algorithm, the sensor provides a two-dimensional image of the surveillance area at intervals T, the surveillance region scenario considered in this paper is shown in Fig. 1. As shown in Fig. 1, every image consists of N c = N r × N θ resolution cells [30]. The measurement matrix can also be denote as a vector form k can be denoted as [30], where ω (j) k is the complex Gaussian noise with covariance Ŵ . ϕ k,i is the unknown phase assumed to be uniformly distributed over the interval [0, 2π) at time k. h k (j) (x k.i ) denotes the contribution intensity from the target cell j. H 1 denote the assumption that there are N k targets, and H 0 denote that there is no target. ρ k,i is the amplitude of the target, in this paper, which is modeled using the Swerling family of target amplitude fluctuation models, i.e., 1 Under the Swerling 0 model, the amplitude ρ k,i of each target is equal to the parameter ρ i . 2 For the Swerling 1 model, the amplitude ρ k,i of the target is assumed to submit to a Rayleigh distribution, E[ρ 2 k,i ] = 2σ 2 ρ i , the PDF of ρ k,i is as follows . 3 When the amplitude fluctuation type is Swerling 3 and ρ k,i follows a chi-square distribution with four degrees of freedom, E ρ 2 k,i = ν ρ i , the PDF of ρ k,i is as follows For equation 3, range and bearing radar covering the surveillance area in polar coordinates are considered. For the range, it is assumed that the transmitted pulse is a linear frequency modulated signal with bandwidth B and duration T e , and the range ambiguity function is given by the [34] where τ l = 2(r k − r l )/c , c is the speed of electromagnetic waves; r k = x 2 k,i + y 2 k,i , r l = r min + l + 1 2 r , l ∈ [0, N r − 1] , r = c 2B is range resolution. At the receiver side, the radar consists of a linear phased array, where the spacing of the N a antenna is 2 and is the wavelength of the carrier frequency. Then the bearing ambiguity function is given by the [31] as The ambiguity function in the range-bearing cell (l, m) is The set of complete measurements collected up to time is given by

The optimal Bayesian filter
The finite sets statistics (FISST) proposed by Mahler provides a mathematical framework for the problem of detection, tracking and classification of multi-target under unified Bayesian paradigm [11]. In FISST framework, given the measurements Z k at time k and the history of measurements Z 1:k−1 , then the optimal multi-target Bayesian filter propagates RFS based on the posterior probability density p k|k ( X k |Z 1:k ) , its prediction and updating steps are as follows where considering the function f: F(X ) → R , F(X ) denotes the set containing all subsets of X , the integral of the function f over a closed subset S of X is given by TBD uses the raw measurements generated by sensor. The measurement is an image composed of thousands or even millions of pixels, however the target does not affect all the pixels, but only occupies a small area of the image, which leads to an estimation of inefficient if we use p( Z k |X k ) directly. Define the proportion of the likelihood p( Z k |X k ) of the existing target to the likelihood p 0 (Z k ) of no target [29].
It can be seen from the (13) that the calculation of the likelihood ratio is limited to the area around the target, which improves the efficiency of the calculation.

Basic MB-TBD filter
A MB RFS X is a union of M independent Bernoulli RFS X that represents a fixed number M of targets. Therefore, the set r (i) , p (i) M i=1 denotes the MB RFS, where r (i) is existence probability and p (i) is the spatial distribution of the i-th Bernoulli distribution. The probability density π (X) of a MB RFS is given by [11] In [20], it is shown that the MB filtering method under RFS provides an accessible solution to the multi-target estimation problem for image measurements. If the posterior probability density π k−1 is in the form of MB, then the prediction π k|k−1 and the update π k|k in the case of non-overlapping targets will also to be in the form of MB. Hence, the prediction (10) and the update (11) can be approximated in the following way. MB-TBD prediction: Given the posterior MB -TBD parameters at time k − 1 , the predict MB-TBD parameters are where where given previous state ζ , f k|k−1 (· | ζ ) and p S,k (ζ ) denote the single target transition density and probability of target existence at time k, respectively. r is parameters of the MB RFS of births at time k. MB-TBD update: Given the predict MB-TBD parameters , the update MB-TBD parameters are where Each pixel cell received by the radar sensor is a complex random variable containing both amplitude and phase information. The existing MB-TBD algorithm is the problem of target tracking in the non-fluctuating target, where the implementation only considers the amplitude information and ignores the phase information. In this paper, we consider two implementations of MB-TBD in the case of target amplitude fluctuation, one is the SLR under the squared radar measurement, and the other is the CLR under the consideration of phase information.

Calculation of the SLR and CLR
The measurement (3) contains the unknown phase ϕ k,i and amplitude ρ k,i of the target. They are unknown, therefore, the likelihood p(Z k |X k ) or p |Z k | 2 X k cannot be directly calculated, and the phase ϕ k,i and amplitude ρ k,i need to be marginalized.

Calculation of the SLR
Denote the squared modulus vector of the complex measure is Assuming that pixel values are distributed independently and that the target's influence area does not overlap, then, When the target does not exist, the likelihood p 0 (Z k ) is [29] (1) For Swerling 0, the L z can be written as [29] where I 0 (·) is the modified Bessel function of the first kind.
(2) Modify the L z of the Swerling 1 model in the literature [30], then, then, where (3) Different from the expression of L z by Swerling 3 in [30], this paper re-derives Swerling 3, where the derivation of Eq. (25) is analogous to Eq. (23).

Calculation of the CLR
The SLR of the square measurement does not consider the spatial coherence of the complex measurement, and for the case of amplitude with no fluctuation, the SLR also needs to calculate a Bessel function with a significant point spread function contribution for each pixel. Experience shows that the calculation of the Bessel function is the most computationally expensive part of the TBD algorithm. In summary, SLR using squared measures is not optimal, in order to avoid these drawbacks, in [29][30][31][32], the joint likelihood of the entire sensor image is given, taking into account the spatial coherence of the phase and spatially correlated noise. The calculation of CLR not only improves the performance of target detection and tracking, but also effectively reduces the computational complexity. The likelihood p 0 (Z k ) of no target is [29] This subsection gives the expression for the CLR for the three Swerling types [30].

Adaptive birth distribution of MB-TBD filter with amplitude fluctuating targets
The traditional MB-TBD assumes that the birth RFS is known as a priori and initializes the multi-Bernoulli filter with a priori knowledge (KpMB-TBD). However, when the prior knowledge is completely unknown and the SNR is low, the continuous missdetection of MB components, may lead to the ineffective of KpMB-TBD. For the environments of low SNR, missed detection is a common phenomenon. If the birth target appears in a state area that is not covered by the true target birth intensity, it is difficult for KpMB-TBD to find the targets, even if a large number of birth targets cover the whole space. The general strategy is to let the birth distribution cover the total scenario. However, this requires a large number of particles to represent the birth model, this method is very inefficient although it is feasible. In [25,32], it is shown that for the MB filter, the adaptive birth distribution must follow the MB distribution, and the MB birth distribution π Ŵ,k at time k + 1 depends on the measurement Z k at the previous moment, i.e., Then the existence probability of each Bernoulli component is where Ŵ, k+1|k is the expected number of targets for newborn at time k + 1 and r Ŵ,max ∈ [0, 1] . Since the mean cardinality of MB RFS is given by the sum of the existence probabilities, the mean cardinality of the new birth targets is For target kinematics of amplitude fluctuation, the target intensity will change at any time. In the target tracking process, the weaker target is often ignored, which leads to an error in the estimation of the number of targets. Figures 2 and 3 show that in the case of SNR=8dB, considering the influence of noise on the target amplitude fluctuation type Swerling 3.  Fig. 3 top, and Fig. 3 bottom shows the strongest echo position in both backgrounds. Figures 2 and 3 intuitively reflects that for Swerling 3 type, the target echo is completely swamped by noise at part of the time, which may lead to the ineffective of filter when the a priori information is known. Also, for fluctuating targets, the maximum echo information detected by the receiver in a low SNR environment may be non-surveillance targets, which will affect the judgment of the filter. While for Swerling 0 type with no fluctuation of the target, the real target echo amplitude is stable. The above only draws the fluctuation of Swerling 3, for another fluctuation of Swerling 1, the target amplitude fluctuation range is still very wide. It is also important to pay attention to the fact that in amplitude fluctuation multitarget tracking, the amplitude among targets are different, and the target with stronger echoes will affect the target with weaker echoes, and the KpBer-TBD cannot solve the amplitude fluctuation problems of weak multi-target tracking.
In order to solve the above difficulties, this section draws on the idea of successivetarget-cancellation (STC) in [35], by making a balance between low false alarm rate and high detection probability, and proposes a MB-TBD based on measurement likelihood adaptive birth distribution (LABer-STC-TBD).
To facilitate the following calculations, transform vector Z k into matrix Z k , i.e., The idea of LABer-STC-TBD firstly select previous instants measurements to generate birth targets adaptively, and not all measurement information is used to drive the birth distribution, i.e., The threshold γ can be calculated by giving the false alarm P fa of the resolution cell, the threshold γ is used to avoid information redundancy and to discard measurements that may originate from noise or clutter. In this paper, a threshold γ more suitable for the case of target fluctuation is modified on the basis of the [36].
Since each Bernoulli component represents a target, and this paper considers target detection and tracking in the case of non-overlapping targets, the measurements are divided after selection so that the particles in each cluster are generated from measurements around the target diffusion location, which not only avoids particle clutter in the components but also reduces the overlapping components as much as possible. The adaptive birth distribution can be expressed as Using (37) to calculate the existence probability of each cluster, the Bernoulli components whose existence probability less than 0.5 are eliminated. The detected component is used to modify the existing measurement to eliminate the influence of the component on the detection of the remaining target. When the detected target acts on the resolution cell (l, m) , the measurement of the resolution cell (l, m) after eliminating the effect of this target At this point the sensor gets the repaired measurement set within the whole N r × N θ scenario as After returning to the above steps and increasing the false alarm rate until the presence probability of all Bernoulli components is below 0.5. It is important to note that (39) holds because (3) of the measurement equation is established and the assumption that the target contribution to the intensity is additive. A modification of (39) is needed for the measurements that are not additive models. Algorithm 1 gives the LABer-STC-TBD algorithm for adaptive birth multi-Bernoulli density r A ∩ B = ∅ , and for any element in A are in the circle whose radius is R (R is set by the target diffusion intensity) and the center is . The number of elements in Ẑ with 1 is eliminated, i.e., the positions with no energy scattering are eliminated. Figure 4 reflects the point spreading of the target. Figure 4a indicates that the target energy spreads to the neighboring pixels around the target location with little spillover to more distant locations, and the white circle in the image is centered on the target location. Figure 4b indicates the measurement from the receiver with SNR = 10 dB. As can be seen from Fig. 4, the real target affects the surrounding pixel locations, while the noise locations are disorderly, even if the SNR is so low that the target location intensity is treated as noise in the first cycle, there is a chance of detecting the real target after the repair of the measurement. In Step 5-6, the target distance and bearing are converted to positions in Cartesian coordinates, and the SMC implementation is considered to generate particles in the near of the real target. In Step 7, r (i) Ŵ,k = r start is to cope with the update in Step 8-13, not the existence probability of the birth component, which is calculated based on equation 37, and the birth particle in Step 7 is not in the position of selection measurement, but contains some disturbances, the reason is that the target will be affected by noise in the process of moving, which will add some disturbances to prevent the estimated target from leaving the real trajectory in the moving process. Remove the Bernoulli components with the existence probability lower than 0.5 in the Step 8-13 and Step 14-18, if it leads to the disappearance of all components, then the search is over, otherwise the target location information will be found and the surrounding pixels will be eliminated for the next iteration. It should be noted that in the first search process, if N i = 1, i = 1, 2 . . . , N , then it shows that there is no newborn target at this time, and the MB filter will continue to iterate to the next moment with surviving targets.

Track merging
After the update of traditional MB-TBD, the Bernoulli components whose existence probability lower than the threshold H merge will be discarded, so as to reduce the increase of Bernoulli components. To solve the problem of linear growth of Bernoulli components after updating, an algorithm for trajectory merging based on distance and particle weights is proposed.
The Bernoulli component within the threshold represents the same target, then the MB components can be expressed as: where M 1 + M 2 + · · · M N = M k|k , and r , j = 1, 2, . . . M k|k denotes the Bernoulli components of the same target, the M k|k components are divided into a total of M N classes, 1 ≤ M j ≤ M N . Subsequently, the components with the highest existence probability of the same target are reserved, then the Bernoulli components before selection are denoted as , and there will be only one component of the same target after selection, , and |·| denotes the cardinality of the set.
If the existing probabilities are the same between the components, i.e., r i j k = r i m k , and the distance between the Bernoulli components is within the threshold, then the two components are combined into a new one, the details are depending on p i j k and p i m k , the specific algorithm is given in Algorithm 2.
Algorithm 2 gives the implements of merging algorithm to merge trajectories. Suppose that the MB filter eliminate the components whose existing probability below the threshold H merge after the update step and the total number of Bernoulli components is within T max . Note: Step 2 truncates the particles in each component and discard particles with low weights to avoid errors in the target estimation process.
Step 3 is to resample the truncated particles and select the particles with high weights to represent the Bernoulli components.
Step 7-13 is to remerge the probability densities of the hypothetical components with the same probability within the threshold to obtain the new components, and the components outside the threshold are retained.

SMC implementation of MB-TBD filtering with target amplitude fluctuations
In this section, the SMC implementation of the amplitude fluctuation target MB-TBD is given.
SMC Prediction: Given a MB posterior density π k−1 = r at time k − 1 , each spatial probability density p (i) k−1 , i = 1, . . . , M k−1 can be represented by a set of particles with weights Then the predicted multi-Bernoulli density can be expressed as follows The birth Bernoulli component r is given by the LABer-STC-TBD algorithm. SMC Update: Given the predicted MB density π k|k−1 = r at time k, each spatial probability density p (i) k|k−1 , i = 1, . . . , M k|k−1 can be represented by a set of particles with weights w Then the updated multi-Bernoulli multi-target density π k = r can be expressed as in section 3.1, and for SLR L z |Z k | 2 x (i,j) k|k−1 in section 3.2. Select the likelihood ratio under different measurement type and fluctuation type as needed.
Resampling and Implementation Issue: Analogous to the standard multi-Bernoulli filter, each Bernoulli component resamples the particles follow the update step, in order to reduce the increasing number of trajectories, the components whose existence probability lower than the threshold H merge are discarded. However, this cannot accurately estimate the number of Bernoulli components, especially when the birth Bernoulli component can also accurately estimate the true position of the target, it will produce cardinality bias.

Experiments
This paper addresses the detection and tracking of fluctuating targets such as stealth aircraft for range-bearing surveillance radar. The Swerling model is more effective for missiles and aircraft. In this section, the effectiveness of the method is verified by Monte Carlo simulation experiments.

Measurement model parameters and multi-objective error estimation
Assume that the noise covariance is Ŵ = 2σ 2 n I N c , SNR = 10log 10 E(ρ 2 k ) 2σ 2 n , r min = 100 km , r max = 120 km , θ min = −75 • , θ max = 75 • , N r = 300 , N θ = 100 , σ 2 = 0.5 , B = 150 KHz , T e = 6.67 × 10 −5 s , N a = 55 , = 3 cm , c = 3 × 10 8 m/s . In this paper, E ρ 2 The Optimal Sub-pattern Assignment (OSPA) [30] is used to evaluate the performance of the algorithm, and the OSPA metric evaluates the estimation error of the x, y = min c, x − y , and n denote the set of permutations on {1, 2, . . . , n} . Then, for p ≥ 1 , c > 0 , if m ≤ n , OSPA is defined as follows [37] if m > n , d (c) In the simulation experiments of this paper, set p = 1 , c = 1000 . The smaller the OSPA value, the more accurately the number of targets and the state estimation are indicated.

Scenario 1: experiments comparing LABer-STC-TBD and KpMB-TBD algorithms
Suppose that targets move in a straight line at a constant speed, there are three targets in the whole scenario, the duration is K = 30 s , the targets are initialized at different positions, the target state includes plane position and velocity, the model is given by the following equation To verify the broadness and effectiveness of the LABer-STC-TBD algorithm as well as the track merging algorithm in the case of target amplitude fluctuations, the LABer-STC-TBD algorithm and the KpMB-TBD algorithm are used to detect and estimate multiple targets based on different amplitude fluctuations. The multi-Bernoulli density of the birth process is where p 1000, 500, 1000, 500] T 2 . The real trajectory on the two-dimensional plane is shown in Fig. 5, and the initial position of the target is shown in Table 1.
In SNR=9dB, the number of particles for each new birth Bernoulli component is 1000, and the algorithm simulations are compared under different amplitude fluctuations using 100 Monte Carlo experimental simulations. The LABer-STC-TBD algorithm is compared with the KpBer-TBD algorithm implementation under different amplitude fluctuations considering four filters as follows.  The results confirm that the LABer-STC-TBD algorithm can accurately estimate the positions and number of targets, although the estimation error of the KpBer-TBD algorithm is smaller than that of the LABer-STC-TBD algorithm at the initial moment, but the error increases rapidly and eventually show a divergence. The reason is that the KpBer-TBD algorithm is given the correct initialization and the prediction of the state transfer equation happens to be close to the true state for the initial few moments, after which the target state prediction decreases as noise as well as uncertainties such as clutter interfere, eventually leading to worse and worse tracking and higher and higher errors. The same is true for the estimation of the number of targets. In the process of target amplitude fluctuation, the target echoes will be completely annihilated in the noise part of the time, and the LABer-STC-TBD algorithm will use the previous moment's measurement information to search for the true target location as much as possible, making the estimated number accurate and the filter can keep iterating even if the noise drowns out all the target information at a certain moment. If the KpBer-TBD algorithm does not estimate the target at a certain moment effectively, and the target is moving faster, the target will not be estimated at the next moment, resulting in the failure of the filter. And Figs. 6, 7 and 8 also show that the filter with the CLR outperforms the filter with the SLR. Comparing these four filter algorithms from Figs.6, 7 and 8, the LABer-STC-TBD algorithm under the CLR has the best effect, proving that the phase information can improve the detection tracking performance of the MeMBer-TBD algorithm, and the selection of the measurement threshold of the birth algorithm reduces the computational complexity, and the selected measurements generated by the real targets are based on the target amplitude, to the largest degree. Using the likelihood ratio to select Bernoulli components to eliminate the wrong estimation at the previous moment, the STC idea is introduced to avoid the influence of strong and weak echoes of the target at the same moment and to maximize the exploitation of the target echo information. From Figs. 7b and 8b, it can be seen that the LABer-STC-TBD with SLR is not accurate for target number estimation in the target with amplitude fluctuation type Swerling 1,3, due to the fact that the square measure ignore the phase information of the target, and comparing with Fig. 6b, it can be seen that the loss of phase information is extremely obvious for fluctuating targets.

Scenario 2: experiments comparing the likelihood ratios of CLR and SLR in MB-TBD
Scenario 5.2 demonstrates the advantages of the LABer-STC-TBD algorithm in the MB-TBD with amplitude fluctuations, and the tracking effect of the CLR is better than the SLR in the target tracking process under the same conditions. To further compare the advantages and disadvantages of the CLR and SLR as well as the reasonableness and superiority of the LABer-STC-TBD algorithm in the MB-TBD, targets are assumed to make a turn at a constant speed. There are five targets in the whole scenario, the time duration is K = 100 s , the targets are initialized at different positions, and the target state variables is x k = x T k , ω k T include the plane position as well as the velocity is x T k = p x,k ,ṗ x,k , p y,k ,ṗ y,k T and the turn rate is ω k . The state transfer model is   Figure 9 shows the real trajectory on a two-dimensional plane. Different targets have different turn rates, and the initial position of the target is shown in Table 2, where wturn = 2π 180 . The targets are tracked in two scenarios with SNR = 7 dB and SNR = 5 dB, and the number of particles of each birth Bernoulli component is 1000, using 100 Monte Carlo experimental simulations to compare algorithmic simulations considering the LABer-STC-TBD algorithm under different amplitude fluctuations. The following two filters are considered. The results are shown in Figs. 10, 11 and 12.
The detection performance of different MB-TBD strategies for Swerling type 0, 1, 3 targets is shown from Figs. 10, 11 and 12, respectively. The advantage of the CLR algorithm at low SNR ratio is verified, and for all detection results, the LABer-STC-TBD using the CLR outperforms the LABer-STC-TBD using the SLR. The introduction of the phase information improves the system performance, which makes it possible to detect and track the target more accurately at low SNR. Moreover, under the same fluctuation, the CLR only needs to perform the calculation of Bessel function once, while the SLR needs to perform multiple Bessel functions, which effectively reduces the computational complexity. Figure 10 shows the performance of LABer-STC-TBD algorithm with no fluctuations in the target, and both SLR and CLR can estimate the target state and the number of targets more accurately. Figures 11 and 12 show performance of the LABer-STC-TBD algorithm with fluctuation type Swerling 1, 3, which reflects the advantage of the CLR over the SLR. However, in the case of low SNR, the estimated number of targets is often lower than the true number of targets, the reason is that the fluctuation of the target amplitude lead to the targets annihilated in the noise.