{'fraction_non_alphanumeric': 0.03013045878490281, 'fraction_numerical': 0.00993694976319842, 'mean_word_length': 4.575571035153086, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 9, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Multisensor Management Method for Ground Moving Target Tracking Based on Doppler Blind Zone Information

Tracking ground moving target with sensors proves to be a challenge due to the uncertainty of target motion area, the existence of Doppler blind zone (DBZ), and the complex terrain. In this paper, a multisensor management method based on DBZ information is presented, in which the available sensors are selected to obtain the best operational revenues for ground moving target tracking. First, the ground target motion model is established considering the off-road/on-road state based on road topology information. Second, the sensor measurement model is given combined with the DBZ information, and a decorrelation method of measurement noise is proposed. Third, a target state estimation algorithm is derived using particle filter, in which the DBZ information is regarded as prior information. Then, combined with the variable structure interacting multiple model method, an estimation algorithm for tracking maneuvering target is proposed. Furthermore, an optimization model of nonmyopic sensor management is constructed to obtain the best sensor management scheme. Finally, the advancement and effectiveness of the proposed management method are verified in the simulations.


Introduction
With the development of sensing technology and information fusion technology, multisensor systems have been widely used in the military field [1][2][3][4][5]. How to determine a reasonable and effective sensor management method to obtain the best operational revenues has become a research hotspot for many scholars.
In terms of decision-making methods in sensor management, there are two kinds of methods, including the myopic sensor management method [6] and the nonmyopic sensor management method [7]. The myopic method decides the management scheme based on predicting the one-step revenue in the future. On the contrast, the nonmyopic method decides the management scheme based on predicting the multistep revenues, which can obtain a better management performance than the myopic method, but with a large computing time.
At present, scholars usually focus their research on the sensor management method based on the optimization indicator, whose connotation is to set an objective function closely related to the optimization indicator to maximize the desired operational revenues [8,9]. According to the selected optimization indicator, the sensor management methods can be divided into three main categories: information-indicator-based management methods, riskindicator-based management methods, and task-indicatorbased management methods. The information-indicatorbased methods focus on managing sensors to maximize the information gain, thus reducing the uncertainty in the observation. Commonly used information indicators are the Shannon entropy [10,11], the Kullback-Leibler divergence [12], and the Rényi information divergence [13,14]. However, the disadvantage of the methods is that the meaning of the used information indicators is too abstract to describe their concrete physical meaning, which may make it difficult for commanders to understand their connotation accurately. The risk-indicator-based management methods are mainly used in scenarios where sensor resources are scarce. The methods consider that the sensor management is performed to reduce the losses due to measurement uncertainty, rather than to maximize measurement performance. Commonly used risk indicators include the threat assessment risk [15], the target loss risk [16,17], and the sensor radiation risk [18][19][20]. The task-indicator-based management methods focus on the combat task, and the relevant indicators are directly related to the tasks. Unlike the risk-indicator-based method, the methods are applicable when the sensor resources are redundant. Typical indicators are the covariance matrix of target state [21,22], the target detection probability [23], and the posterior Cramér-Rao lower bound (PCRLB) [24,25]. Among the three kinds of sensor management methods mentioned above, methods for target tracking have been most studied, and the representative work is shown in [19,25,26]. In [19], a nonmyopic sensor management method is applied to track aerial targets for the tradeoff between the tracking accuracy and the radiation risk. Mohammadi and Asif [25] present a dynamic sensor scheduling for tracking problems in distributed sensor networks, in which the PCRLB is used to quantify the tracking accuracy. Gostar et al. [26] present a constrained sensor management method for the multitarget tracking based on information divergence, and the labeled multi-Bernoulli filter was proposed for estimating the target state.
However, most of the existing studies about sensor management methods mainly focus on aerial targets and neglect ground targets. In the actual battlefield, there is inevitably a need for tracking ground targets, whose motions are uncertain and often subject to the terrain information [27][28][29]. Inground target tracking, the most widely used terrain information, is the road information, and the representative work is shown in [30][31][32][33][34][35]. In [30], an on-road target tracking is proposed, in which the linear roads are mapped to the ground coordinate system and the road network is established as a one-dimensional linear mode. Koch et al. [31] propose a road map extraction method to reduce the error of the road map. In [32], the road network is modeled as constant curvature segments, and the nonlinear roadconstrained Kalman filtering is used in target state estimation. In [33], the widths of road segments are considered in the state model of road target which influences the transition of the target state. Zheng and Gao [34] introduce the Gaussian mixture probability hypothesis density filter in the ground target tracking and proposed a multitarget tracking method in clutter combined with the road constraints. In [35], a ground moving target indication (GMTI) tracking problem is presented, and a variable structure interacting multiple model (VSIMM) particle filter was applied in the switching problem of road segments.
Meanwhile, the ground target detection sensors, such as the GMTI radar, mostly used the pulsed Doppler techniques to reduce the effect of clutters [36]. However, the Doppler blind zone (DBZ) problem has become a great challenge in target tracking. The low-Doppler targets, whose Doppler magnitude falls (i.e., the radial velocity of the target) below the minimum detectable velocity (MDV), cannot be detected by the sensor [37]. This results in a significant decrease in tracking accuracy. For the target tracking in the presence of the DBZ, most previous works focus on improving the filter and data fusion algorithms [38][39][40], but few works utilize sensor management to solve it.
To solve the problem mentioned above, we propose a multisensor management method for ground moving target tracking based on the DBZ information. The study focuses on the following points.
First, a more realistic ground moving target tracking scenario is considered in which the target can move in different motion areas and its movement is constrained by road topology information. That is to say, the target motion state is divided into two categories: off-road and on-road.
Second, a sensor measurement model under measurement uncertainty is presented considering the DBZ information. Meanwhile, Lei and Han in [41] indicate that the measurement noise of the range and the radial velocity are statistically correlated. To solve this problem, we propose a decorrelation method of measurement noise.
Third, combined with the particle filter (PF) algorithm, a target state estimation algorithm based on the DBZ information is proposed, in which the DBZ is regarded as prior information. Then, a VSIMM-PF-DBZ algorithm is proposed for tracking the maneuvering target combining the VSIMM method.
Finally, we introduce the PCRLB to quantify the tracking accuracy in the future and establish a nonmyopic sensor management optimization model to obtain the optimal management scheme.
The research framework of this paper is as follows. Section 2 describes the sensor management problem. The ground target motion model and the sensor measurement model are established in Section 3 and Section 4, respectively. In Section 5, the target state estimation algorithm based on the DBZ information is proposed. In Section 6, the VSIMM-PF-DBZ algorithm is presented. The sensor management optimization model is given in Section 7. Section 8 makes some simulations to illustrate the effectiveness and advancement of the proposed management method. Finally, conclusions are given in Section 9.

Problem Description
Sensors are Doppler radars in this paper. The target tracking scene is shown in Figure 1. Assume that there areNsensors to track a ground moving target, which can move on the road or off-road. The data processing center unifies the measurements obtained by the sensors, searches the sensor management scheme with the best revenues, and sends the corresponding commands to control the working of the sensor system. Note that we use the nonmyopic sensor management method in which the management scheme is decided based on the cumulative revenues over a time horizon in the future [7,42].
At time k, the sensor management scheme is denoted as A k = ða n k Þ N×1 , where a n k indicates the working state of sensor n. If sensor n is activated to track the target at time k, then a n k = 1; otherwise, a n k = 0. Then, the management scheme sequence over the H − step time horizon is denoted as For the convenience of presentation, we set the following two constraints related to A k:k+H−1 : The above constraints can be written as where tðA k+h−1 Þ represents the continuous working time of the switched sensor at time k + h − 1.

Ground Target Motion Model
The ground target state at time k is X k = ½x k , _ x k , y k , _ y k T , where x k and y k are the position coordinates in the realworld coordinate system and _ x k and _ y k are the corresponding velocities. When a target moves on the road, its motion state is mainly affected by the road topology [43]. For the off-road target, its motion is relatively free. Therefore, different state models are required to describe the characteristics of the target motion in different areas.
3.1. Off-Road State Model. The off-road state model can be described as [19] where F off -road i is the state transition matrix of motion model i, G off -road i is the corresponding process noise gain matrix, ε off -road k is the zero means Gaussian process noise with covariance matrix Q off -road k = diag ðσ x , σ y Þ, and τ is the sampling interval. In this paper, two typical motion models are considered, including the constant velocity (CV) model and the constant turn (CT) model. Thus, the corresponding state transition matrixes and process noise gain matrices are as follows [24]. where φ is the turn rate of the target.

On-Road State
Model. The information of road segments can be collected from geographic information systems (GIS) [43]. Then, the mathematical model of the road network is established by using relevant information to represent the road network as a connection of many road segments, as shown in Figure 2. When the target is moving continuously on the road, it can be considered as moving along the road centerline without deviating largely normal to it. When the target is moving along the centerline of road segment l whose start point is ðx start l , y start l Þ and end point is ðx end l , y end l Þ, its state constraint can be described as with where h½ _ x k , _ y k T , ς ! i represents the angle between two vectors, ς ! is the direction vector of the road segment direction, a l , b l , and c l are the coefficients of the centerline equation. Equation (6) considers the state constraint when the target moves along the centerline, but the target may slightly deviate to the centerline because of the presence of the process noise. Combined with equation (6), the target state where d½ðx k , y k Þja l x + b l y + c l = 0 represents the distance from point ðx k , y k Þ to line a l x + b l y + c l = 0 and Δd and Δv are the deviation threshold of the distance and velocity angle, respectively. While satisfying the above state constraint, the on-road state model can be described as In this paper, we consider that the on-road target mentioned only moves in CV model. Therefore, F on-road We denote the variances of the process noise along the road and orthogonal to the road as σ 2 ∥ and σ 2 ⊥ (σ 2 ⊥ ≪ σ 2 ∥ ), respectively. The covariance matrix Q on-road k after projecting it to the X-Y coordinate system can be written as [33] Q on-road where θ l represents the orientation angle of road segment l. Note that the state transition of the on-road target is also determined by the road segment on which they are located, while the state transition of the off-road target is mainly determined by its motion model (CV or CT).

Motion Area Transitions.
Obviously, the motion areas of the target can be switched between on-road and off-road. In this paper, parameterm k is used to indicate the area in which the target is moving at timek, wherem k = 0represents the target moving on the off-road area andm k = 1represents the target moving on the road. The transition process of motion area can be approximated as a Markov process, and the cor-responding transition matrix can be stated as with Here, d entry k and d exit k represent the distance to the nearest entry point and nearest exit point of the road network, respectively. ϑ is a probability parameter (in this paper, ϑ = 0:04).

Sensor Measurement Equation.
The Doppler sensor can obtain the distance, azimuth, and radial velocity information of the target, and the corresponding measurement equation can be expressed as [19]  with where the measurement information ½r n k φ n k _ r n k T represent the range, azimuth, and radial velocity of the target at time k and the noise υ n = ½υ n r υ n φ υ n _ r T represent the Gaussian measurement noise with zero means. x n s and y n s are the position coordinates of sensor n. If the target is not detected successfully, the sensor cannot obtain the measurement information. When the DBZ is not considered, the detection probability of sensor n can  Journal of Sensors be written as [44] with where SNR is the signal to noise ratio and p n f , SNR n min , and R n max represent the false probability, the minimum detectable SNR, and the maximum detection distance of sensor n, respectively.
When the DBZ is considered, the target with radial velocity below than MDV will not be detected, which is an important priori information for sensor management. The corresponding detection probability can be calculated by [37] where V n MDV is the MDV of sensor n and n c ðX k Þ is called the clutter notch function (see [37] for details).
Combined with the detection probability, the measurement equation can be written as where η represents a random number, which is taken as 0 or 1 according to the detection probability p n d ðX k Þ, that is

Decorrelation Method of the Measurement Noise.
According to [41], the range measurement noise υ n r and the radial velocity measurement noise υ n _ r are statistically correlated, which cannot be ignored in the process of target tracking. Hence, we propose a decorrelation method of measurement noise to improve the accuracy of target tracking in this paper.
We define σ 2 ðυ n r Þ, σ 2 ðυ n φ Þ, and σ 2 ðυ n _ r Þ as the variances of υ n φ , υ n φ , and υ n _ r , respectively. The correlation coefficient of σ 2 ðυ n r Þ and σ 2 ðυ n _ r Þ is denoted as ϖ. Then, the covariance matrix of measurement noise can be written as Rewrite equation (20) as with To eliminate the correlation of the elements in R n , we use the Cholesky decomposition method and construct the decomposition matrix with where I 2×2 is an identity matrix. Multiplying both sides of equation (18) simultaneously left by B n , the decorrelated measurement equation is obtained as followsZ whereψ n k represents the pseudomeasurement of the radial velocity and ψ k and υ n ψ are the measurement value and noise ofψ n k , respectively, which can be calculated by Obviously, υ n ψ is the Gaussian noise with zero means, whose variance can be expressed as Therefore, the new measurement equation can be The corresponding covariance matrix of the measurement noise is

Target State Estimation Algorithm Based on the DBZ Information
In the process of sensor management, estimating the target state at each moment is a prerequisite for system decisionmaking. When the DBZ is considered, the sensors will lose the measurements in some cases, and the target state only can be obtained by prediction based on the recurrence of the target motion models. At this time, the tracking accuracy will be greatly reduced, which is not conducive to the stable tracking of the target. It is known that the target state in the DBZ satisfies a certain constraint relationship, that is, the radial velocity of the target is below than the MDV of the sensor, which can be used to improve the tracking performance as prior information when the sensor cannot obtain the measurements. Therefore, we propose a particle filter algorithm based on the DBZ information (PF-DBZ) for target state estimation in this paper.
For the convenience of expression, the multimodel and motion area transitions of the target are not considered in this section.

5.1.
Process of the PF-DBZ. The generated particles are divided into two categories: one is the particles outside the blind zone called unconstrained particles with a number of M UN , denoted as fX i UN,k : i = 1, 2, ⋯, M UN g; another is the particles in the DBZ with a number of M DBZ , called DBZ particles, denoted as fX i DBZ,k−1 : i = 1, 2, ⋯, M DBZ g, which are constrained by the DBZ information. At time k, the occurrence probability of the two kinds of particles is p UN,k and p DBZ,k , respectively. Then, the processes of PF-DBZ are as follows.
For the unconstrained particles, the prediction statẽ X i UN,k ði = 1, 2, ⋯, M UN Þ can be calculated based on the state model mentioned in Section 3: For the DBZ particles, if p DBZ,k−1 > 0, the prediction statẽ X i DBZ,k ði = 1, 2, ⋯, M DBZ Þ can also be obtained by equation (31), but the state must meet the following DBZ information constraint.
At time k, whether sensor n can obtain the measurement is an uncertain event. Therefore, different cases will be discussed as follows.
(2) Normalize the weights: Journal of Sensors Case 1. Measurement Z n k is existing.
Case 2. Z n k is not existing and p DBZ,k−1 = 0.
Case 3. Z n k is not existing and p DBZ,k−1 > 0 (continuous missing detection).
The state estimation and covariance matrix after weighted average are given by whereφ n k is the prediction value of the azimuth. When _ r n k < −V n MDV , the correction method is expressed as 5.3. Judgment Method of Measurement Loss Causes. In ground moving target tracking, the occlusion of obstacles may also cause continuous missed detection. Therefore, it is necessary to distinguish whether the target is occluded by obstacles or enters the DBZ. According to the motion characteristics of the ground target, the radial velocity has a general tendency to decrease gradually before the target enters the DBZ. Therefore, a radial velocity sliding window is used to judge measurement loss causes in this paper. The length of the sliding window is set to 5, and the historical radial velocity f_ r n k−1 , _ r n k−2 , _ r n k−3 , _ r n k−4 , _ r n k−5 g is stored in the sliding window at time k. Then, we define ∇_ r n k−i as the radial velocity change rate, which is expressed as where τ is the sampling interval.
Using the three-fourths judgment rule, when three inequalities in equation (40) satisfy the condition, it is judged that the reason of losing measurement is entering the DBZ. In this case, the PF-DBZ and correction method of target state proposed in this paper can be used. Otherwise, it shows that the target is occluded by obstacles, and the target state only can be obtained by prediction based on the recurrence of the target motion models.

VSIMM-PF-DBZ Algorithm
The target state estimation algorithm mentioned in Section 5 is suitable for the target with a single motion model. However, according to Section 2, the motion of the ground target involves the switching of the motion model, the located road segment, and the motion area, which is difficult to describe its motion characteristics by only a single model. Hence, we combine the PF-DBZ algorithm with the VSIMM (see [45] for details) and propose a VSIMM-PF-DBZ algorithm to estimate the state of maneuvering target. In the VSIMM method, the motion model set of the target is variable at When the target moves in the road network, its motion is mainly constrained by the road topology, so the road information can be used to assist the target estimation. Since the road information can be obtained from GIS in advance and generally does not change in target tracking, the motion model set update strategy can be designed based on the road information.

Motion Model Set Update
Strategy. The motion model of the target moving on a certain road segment will generally not change, and only at intersections will the target switch the road segment and its motion model will be changed. Therefore, to update the model set of the next time, it is necessary to determine which road segment the target will be located at the next time and whether the target is close to the intersection.
As seen from Figure 3, the possible location of the target can be obtained according to the target stateX n k−1 . Then, it can be judged whether the target is close to the intersection, such as the ellipse region in Figure 3. The elliptic region can be expressed as an inequality where P pos k−1 is the position component of P k−1 and β represents the gate threshold (in this paper, β = 6).
If intersection i with coordinates ðx i intersection , y i intersection Þ satisfies inequality (41), it means that the target is close to the intersection, and the model set needs to be updated. The processes of updating J on-road k−1 are as follows.
Step 1. Get the target state ðX k−1 , P k−1 Þ and model set J on-road Step 2. Judge whether the target has closed to the intersection according to equation (41). If the target is close to the intersection, the corresponding motion models of all the links connected to the intersection are added to the model set J on-road k ; otherwise, go to Step 3.
Step 3. Determine the road segment to which the target belongs at time k according toX    Step 1. Initialization.
Obtain the target state ðX k−1 , P k−1 Þ, the motion model set J k−1 , and the parameter m k−1 . Then, according to Section 3.3, calculate the motion area transition matrix T k , and obtain pðm k = 0jm k−1 Þ and pðm k = 1jm k−1 Þ.
Step 2. Update the motion model set.
According to the motion model set update strategy proposed in Section 6.1, obtain the model set J on-road k of on-road state.
Use the PF-DBZ algorithm in the framework of the VSIMM [45] to predict and update the target states and corresponding covariance matrices of all motion models. Then, the target state estimation of the on-road state and off-road state can be obtained, which are stated as ðX

Sensor Management Optimization Model
The core of sensor management is predicting the revenue corresponding to different decision schemes in the future and selecting the schemes with the best revenue [19]. Therefore, it is important to select an optimization index for quantifying the revenue. The trace of PCRLB can predict the theoretical lower bound of the target state estimation error, which is often used to reflect the sensor tracking performance [24]. Hence, in this paper, the trace of PCRLB is used as the optimization index in the sensor management process.
According to the relevant theory of PCRLB [24], the following inequality exists.
where ΨðX∧ k Þ −1 represents the PCRLB ofX k , which is the inverse of the Fisher information matrix ΨðX k Þ. ΨðX k Þ satisfies the following recurrence relation: where the symbol ∇ represents the second-order derivative and D Z k is the Fisher information gain. When the sensor cannot obtain the measurement, there is D Z k = 0. Furthermore, for the Gaussian system discussed in this paper, there is where h _n is the Jacobian matrix of measurement function H _n ðX k Þ and F k and Q k represent the state transition matrix and covariance matrix of the process noise corresponding to the real target motion model at time k, respectively. Obviously, it is impossible to obtain the real target motion model at time k in the process of decision-making, so the motion model corresponding to the maximum distribution probability at the current time is selected to do prediction [24].
In this paper, the sensor management method is nonmyopic, in which the management scheme is decided based on the cumulative revenues over a time horizon in the future. When the decision step is H, combining the sensor management scheme A k:k+H−1 in the time domain ½k, k + H − 1 and the constraints equation (1), the objectivê where tr½ΨðX∧ k , A k Þ −1 represents the trace of PCRLB after executing sensor management scheme A k and γð0 ≤ γ ≤ 1Þ is a discount coefficient that indicates. The optimal solution A opt k:k+H−1 of the objective function is the best sensor management scheme in the time domain ½k, k + H − 1.

Numerical Simulations
As can be seen from Figure 4, in our simulations, six Doppler sensors are used to track a ground moving target that can move on the road or off-road. Besides, there are some obstacles in the battlefield which may obstruct the detection of the All sensors have the same MDV V MDV = 25m/s, the minimum dwell time t min = 3s, the false alarm probability p f = 10 −6 , the minimum detectable SNR SNR min = 15db, and the maximum detection distance R max = 10km. The other parameters of the sensors are shown in Table 1.
The other parameters used in simulations are displayed as follows: the discount coefficient γ = 0:9, the correlation coefficient ϖ of the measurement noise is 0.5, and the number of the unconstrained and DBZ particles is 100 in the PF-DBZ algorithm.
All the results were obtained by 100 trials of Monte Carlo simulations. Figure 5 shows the radial velocities of the target relative to all sensors. Obviously, for each sensor, there is a situation where the radial velocity of the target is below than MDV, which means that the target will enter the DBZ of the sensor at some times and a single sensor cannot achieve continuous tracking of the target. Therefore, it is necessary for tracking the target to use an effective sensor management method.

Determination of the Decision
Step. The decision step H is a very important parameter in sensor management, which can directly affect the effectiveness of target tracking [24]. Hence, the performances of target tracking under different H are studied in our simulations, and the corresponding result is used as a basis to determine the value of H. Figures 6 and 7 show the comparisons of the root mean square error (RMSE) and the root time average square error (RTAMSE) of the target under different H, respectively.

Journal of Sensors
Especially when H = 4, the corresponding errors are larger than the errors under H = 3. The reason for these is that with the increase of H, the prediction error of the target state in the future will also increase, resulting in the increase of the target tracking error. At the same time, the occurrence of measurement uncertainty events such as the target enters the DBZ which will also lead to the increase of the tracking error in the multistep prediction. Therefore, it is not the case that a larger decision step is better for tracking. Based on the above analysis, we choose H = 3 in next simulations.    14 Journal of Sensors (1) The myopic sensor management method (MSMM), in which the management scheme is decided based on one-step revenues in the future (2) The closest sensor management method (CSMM), in which the sensor closest to the target is selected for tracking [19] (3) The fixed sensor management method (FSMM), in which the fixed sensor is selected for tracking in the simulation duration (in our simulations, sensor 1 is selected) The occurrence probabilities of blind zones under different sensor management methods in 100 repeated simulations are shown in Figure 8. For the convenience of presentation, we call the target to enter the visual blind zone (VBZ) when it is blocked by obstacles. Compared with the other methods, we can see that the DBZ and VBZ appear the least times under the PSMM. This is due to the PSMM considering not only the one-step revenue of the sensor system but also the future expected revenue in decision-making, which can better predict the target state and thus select the nonblind sensor for tracking. The MSM only considers the one-step revenue in decision-making which is myopic for tracking, so the corresponding blind zone occurrence probability is higher than those of PSMM. The CSMM and FSMM do not select the sensor based on the tracking revenue in the future, which makes the blind zones occurred many times in the tracking. Figures 9 and 10 show the comparisons of the RMSE and RTAMSE under different sensor management methods, respectively. Obviously, the PSMM can obtain the best tracking performance compared with the other methods. Besides, the position RMSE curves and velocity RMSE curves under PSMM are relatively stable in the simulation duration, and the corresponding error is basically not sharply increased, which means that it can achieve continuous and stable tracking of the target by the PSMM. Figure 11 shows the sensor management schemes of the PSMM, which is the optimal solution of the objective function (47). It can be seen that the minimum continuous working time of the selected sensor satisfies the constraint that the minimum dwell time t min = 3s in tracking, which verified that the management scheme is reasonable and effective.
Furthermore, in order to show that the proposed target state estimation algorithm based on the DBZ information is effective and advanced, the VSIMM-PF-DBZ algorithm is compared with the VSIMM-PF algorithm [43] in which the DBZ information is not be used. It can be seen from Figure 8(d) that the DBZ will appear continuously during 4-8 s under the FSMM. Therefore, we compare the two algorithms during 0-20 s under the FSMM in Figure 12, so as to study the tracking performance of them. Obviously, the tracking errors of the VSIMM-PF-DBZ are less than VSIMM-PF when the DBZ appears (4-8 s), which indicates that the estimation algorithm proposed in this paper can reasonably use the DBZ information to track the target.
8.3. Analysis of the Proposed Motion Model. In this paper, a more realistic target motion model is established based on road topology information. To analyze the effect of utilizing road topology information in terms of tracking performance, the proposed motion model is compared with the model without utilizing any road topology information. Note that if no information is used, equation (2) is always used to describe the target motion in proposed tracking algorithm. Figure 13 shows the comparisons of the RMSE under different motion models. Obviously, the position and velocity RMSE under the proposed motion model are less than those under the motion model without road topology information. Especially when the target moves on the road (0~45 s and 86~100 s), the difference between the results with the two kinds of models is obvious. From the comparisons, we can see that the tracking accuracy is improved by utilizing road topology information to describe the target motion.

Analysis of the Proposed Decorrelation Method of the Measurement Noise.
To analyze the performance of the proposed decorrelation method of the measurement noise, the results of tracking errors corresponding to using and not using decorrelation methods are compared. Note that the correlation coefficient ϖ is set as two values, 0.1 and 0.8, and the other parameters remain unchanged. Figures 14 and 15 show the comparison of the RMSE when ϖ = 0:1 and ϖ = 0:8, respectively. As can be seen, there is little difference between the results of the two methods when ϖ = 0:1. On the contrary, when ϖ = 0:8, the results with decorrelation method are much better than those without decorrelation method. The comparison results show that the proposed decorrelation method is effective, especially when the correlation coefficient is large.

Conclusions
In this paper, a sensor management method for ground moving target tracking is proposed, in which the DBZ information is used to improve the tracking accuracy in the presence of the DBZ. Firstly, a more realistic motion model of the ground target is established, in which the target motion state is divided into two categories: off-road and on-road. Secondly, a sensor measurement model under measurement uncertainty is presented, and a decorrelation method is proposed to solve the problem that the measurement noises are statistically correlated. Third, a target state estimation algorithm based on the DBZ information is proposed, in which the DBZ information is fully utilized. Furthermore, combined with the VSIMM, an estimation algorithm is given to track the maneuvering target with multiple models. Finally, the PCRLB is used to quantify the tracking accuracy in the future, and a nonmyopic sensor management optimization model is established. Simulation results indicate that the proposed sensor management method can track the target accurately by selecting the suitable sensor management scheme at each time. Meanwhile, the proposed motion model and decorrelation method of the measurement noise are verified to be effective. As future work, we will study the sensor management method for multitarget tracking in the presence of the DBZ.

Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.