Next Article in Journal
FIgLib & SmokeyNet: Dataset and Deep Learning Model for Real-Time Wildland Fire Smoke Detection
Previous Article in Journal
Assessing Storm Response of Multiple Intertidal Bars Using an Open-Source Automatic Processing Toolbox
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Parameter Estimation for Uniformly Accelerating Moving Target in the Forward Scatter Radar Network

State Key Laboratory of Complex Electromagnetic Environment Effects on Electronics and Information System, National University of Defense Technology, Changsha 564211, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(4), 1006; https://doi.org/10.3390/rs14041006
Submission received: 15 December 2021 / Revised: 14 February 2022 / Accepted: 16 February 2022 / Published: 18 February 2022

Abstract

:
Passive radar based on the global navigation satellite positioning system (GNSS) has become the focus of attention in the field of radar. A parameter estimation method is proposed in the forward scatter radar (FSR) network based on GNSS to extend the application scenarios. For uniformly accelerating moving targets, only the instant times when the target crosses the individual baselines are used to retrieve the target motion parameters. The target position, velocity, and acceleration information can be obtained. Firstly, the minimum network configuration is derived theoretically. Then, the effects of crossing time error, station location error, transmitting/receiving station deployment, and target height on the accuracy are analyzed through Monte Carlo simulations. Finally, the simulation results indicate that the target position estimation error is in the order of 100 m. This paper provides the fundamental theory of aerial target positioning with a GNSS-based FSR network.

1. Introduction

Passive radar based on GNSS has been researched for more than twenty years. GNSS-based passive radar has the advantages of wide coverage [1,2] and good concealment. However, the great challenge is that the scattering signal is extremely weak. Long-term accumulation and super large aperture antenna are needed for long-distance target detection. This leads to high complexity and huge cost.
As for the detection capability of the GNSS-based passive radar system in practical scenarios, some experimental verification results of various countries are summarized in Table 1, from when global position system (GPS) signals were used to detect civil aircraft, military aircraft, and Russian space stations in 1995 [3]. With the support of the European navigation system authority, a project named “Galileo-based passive radar system for maritime surveillance-spyglass” was established in 2015. In the actual experiments, two different maritime targets were detected by up to 12 GNSS satellites belonging to different satellite constellations (GPS, Global Navigation Satellite System (GLONASS), and Galileo) simultaneously. Among these satellite constellations, the time difference of arrival (TDOA) location method is used to obtain the target position. The long-term moving target indicator (MTI) processing of merchant ships is realized, and the detection distance can achieve 3 km [4,5,6]. Table 1 indicates that the detection distance is very close, and the application scenarios are very limited. There is an urgent need to find a new target location method and expand the application scenarios.
In bistatic radar, when the bistatic angle reaches 180 degrees, the radar cross section (RCS) of the target increases sharply [12,13,14,15,16,17,18,19]. This characteristic can be utilized in FSR. Regarding the parameter estimation in FSR, Professors Cherniakov [20] and Hu [21] have studied some algorithms. However, these algorithms require the arrival azimuth or pitch angle of the forward scattering echo. Therefore, the large aperture multi-beam antenna must be adopted, which causes a high cost-effectiveness ratio in the practical application.
If only the flight altitude of the aircraft target (usually less than 10 km) has to be determined, the problem of the low power limitation for GNSS-based passive radar can be reduced. Studies have shown that there was an apparent disturbance of echo when the target crossed the baselines [22,23,24,25]. By observing disturbances of echoes at multiple receivers, we can detect moving targets.
In this paper, the GNSS-based FSR architecture is proposed. Available target position estimation can be obtained by using the instant times that the target crosses the baselines. It is possible to find an effective method to estimate the crossing time by analyzing the characteristics of echo in the time domain, frequency domain, and time-frequency domain. In addition, the receivers (such as GPS receivers) have a simple structure when only measuring the time. A parameter estimation method based on crossing time was studied in [26], but the transmitters and receivers were required to be deployed at equal intervals and in parallel. The method in this paper removes this limitation. In addition, the significance of the crossing time is clarified, and the minimum configuration of the system is derived theoretically. Various error sources affecting the performance of parameter estimation are analyzed through Monte Carlo simulations, including the crossing time error, station location error, transmitting/receiving station deployment, and target height. According to the simulation results, the application scenario of the system is given, which provides a reference for the deployment of the actual FSR network.
The rest of this paper is organized as follows. Section 2 defines the uniformly accelerated motion observation model of the FSR network. Section 3 describes the parameter estimation method. Section 4 shows some simulation results with different errors which shows the potential application scenario of the system. Our conclusions are given in Section 5.

2. Observation Model in FSR Network

As shown in Figure 1, the FSR network consists of the radiation source of widely existing satellite navigation signals and sparse receiving arrays on the ground. FSR is a particular case of a multistatic radar configuration. It takes full advantage of the forward scattering area of the bistatic radar cross section for target detection, and has a strong anti-stealth ability. Moreover, the network has the following prominent advantages:
(1)
Satellite navigation signals exist widely, and any point on Earth is illuminated typically by 20–30 satellites simultaneously from different angles;
(2)
Satellite navigation signals have high time-frequency synchronization accuracy;
(3)
There is no need to make a single passive antenna very large, so it is simple and easy to deploy;
(4)
Only the height information needs to be detected, which alleviates the problem of long-distance detection.
Figure 1. FSR network structure.
Figure 1. FSR network structure.
Remotesensing 14 01006 g001
When the target crosses the connection between the transmitter station and the receiver, the ideal forward scattering state is formed. The existence of the target can be confirmed by the disturbance of the receiving signal, and the movement of the target can be perceived continuously by the detection of multiple receivers. As for the processing of the receiving signal, the conventional GNSS signal processing method could be used. We can extract the C/A code correlation value for every signal integration interval and confirm the satellite serial number. Furthermore, there is no range resolution, but rather Doppler resolution in the forward scattering region. The signal’s Doppler resolution has been proven to be a linear function of time [25] and approaches zero when the target crosses the baseline, and this feature can be used to estimate the crossing time.

3. Proposed Parameter Estimation Method

The coordinates of the satellites and receiving stations are three-dimensional space coordinates in the real application scenario. Describing the target position in a three-dimensional space can reflect the practical value of the algorithm much better, as shown in Figure 1, assuming that the baselines that the target has crossed correspond to N transmitters and K receivers. According to the order in which the target crosses, the transmitters are numbered T n ( x t n , y t n , z t n ) , (n = 1, 2, …, N), and the receivers are numbered R k ( x r k , y r k , z r k ) , (k = 1, 2, …, K). The coordinate system is taken to be centered in R1. The coordinate of the intersection of the target and the first baseline is ( x 0 , y 0 , z 0 ) , and ( x 0 , y 0 , z 0 ) used in Section 4. To reduce the impact of satellite motion, the observation time should be limited, and assume that the crossed baselines are approximately in the same plane. The error caused by this approximation can be converted into the crossing time error and the station location error.
Then, we convert the three-dimensional (3D) coordinate ( x , y , z ) to a two-dimensional (2D) coordinate (x,y) as shown in Figure 2. Assuming the equation of the plane is a x + b y + c z + d = 0 , the receiver R1 is set as the origin O, with the receivers R1, R2, …, RK placed along the x-axis, while the line passing through the origin and being perpendicular to the x-axis is taken as the y-axis. The direction vectors of the x-axis and y-axis are x = (mx,nx,px) and y = (my,ny,py), respectively. Then, the 2D coordinate can be obtained by calculating the distance from each point to the x-axis and y-axis.
Taking the solution of the ordinate of T n ( x t n , y t n , z t n ) as an example, the orthogonal projection of Tn on the x-axis is
{ x d = m x t y d = n x t z d = p x t
where t = m x x t n + n x y t n + p x z t n m x 2 + n x 2 + p x 2 .
Therefore, the vertical ordinate of T n ( x t n , y t n , z t n ) in the 2D coordinate system is y t n = ( x t n x d ) 2 + ( y t n y d ) 2 + ( z t n z d ) 2 .
For the consideration of simplification, the coordinate of the nth transmitter is Tn = [xtn,ytn]T, (n = 1, 2, …, N), and the coordinate of the kth receiver is Rk= [xrk,yrk]T, (k = 1, 2, …, K). Assuming that the target moves in a straight line with uniform acceleration a in the xOy plane, the initial velocity is v, vx is the velocity in the x-direction, vy is the velocity in the y-direction, and the initial position is (x0,y0). The target coordinates change with time as follows:
x ( t ) = x 0 + v x t + 1 2 a x t 2 , y ( t ) = y 0 + v y t + 1 2 a y t 2
The pair (n,k) corresponds to the nth transmitter and the kth receiver. The start time (tnk = 0) is the moment when the target crosses the first baseline. The target crosses the individual baseline (n,k) at time tnk, and thus the following equation can be derived:
( y t n y r k ) ( x 0 + v x t n k + 1 2 a x t n k 2 ) ( y t n y r k ) x t n = ( x r k x t n ) y t n ( x r k x t n ) ( y 0 + v y t n k + 1 2 a y t n k 2 )
where tnk represents the time required for the target to move from the initial position to the current position. For the first baseline, tnk = 0 and the initial position (x0,y0) satisfies the following equation:
( y t n y r k ) x 0 + ( x r k x t n ) y 0 ( y t n y r k ) x t n = x r k y t n x t n y r k
In a matrix form, Equation (3) can be expressed as
[ y t n y r k x r k x t n ( y t n y r k ) t n k ( x r k x t n ) t n k 1 2 ( y t n y r k ) t n k 2 1 2 ( x r k x t n ) t n k 2 ] [ x 0 y 0 v x v y a x a y ] = x r k y t n x t n y r k
where the six unknowns (target position, velocities, and accelerations) are arranged in a vector. Assuming that P intersections are available, and PNU, with NU being the number of unknowns, the unknown vector will be the solution of a linear equation set, as shown in Equation (6)
[ y t n 1 y r k 1 x r k 1 x t n 1 0 0 0 0 y t n 2 y r k 2 x r k 2 x t n 2 ( y t n 2 y r k 2 ) t n 2 k 2 ( x r k 2 x t n 2 ) t n 2 k 2 1 2 ( y t n 2 y r k 2 ) t n 2 k 2 2 1 2 ( x r k 2 x t n 2 ) t n 2 k 2 2 y t n 3 y r k 3 x r k 3 x t n 3 ( y t n 3 y r k 3 ) t n 3 k 3 ( x r k 3 x t n 3 ) t n 3 k 3 1 2 ( y t n 3 y r k 3 ) t n 3 k 3 2 1 2 ( x r k 3 x t n 3 ) t n 3 k 3 2 y t n P y r k P x r k P x t n P ( y t n P y r k P ) t n P k P ( x r k P x t n P ) t n P k P 1 2 ( y t n P y r k P ) t n P k P 2 1 2 ( x r k P x t n P ) t n P k P 2 ] [ x 0 y 0 v x v y a x a y ] = [ x r k 1 y t n 1 x t n 1 y r k 1 x r k 2 y t n 2 x t n 2 y r k 2 x r k 3 y t n 3 x t n 3 y r k 3 x r k P y t n P x t n P y r k P ]
For the sake of simplicity, Equation (6) can be expressed in a matrix form as follows.
C X = F
where
C = [ y t n 1 y r k 1 x r k 1 x t n 1 0 0 0 0 y t n 2 y r k 2 x r k 2 x t n 2 ( y t n 2 y r k 2 ) t n 2 k 2 ( x r k 2 x t n 2 ) t n 2 k 2 1 2 ( y t n 2 y r k 2 ) t n 2 k 2 2 1 2 ( x r k 2 x t n 2 ) t n 2 k 2 2 y t n 3 y r k 3 x r k 3 x t n 3 ( y t n 3 y r k 3 ) t n 3 k 3 ( x r k 3 x t n 3 ) t n 3 k 3 1 2 ( y t n 3 y r k 3 ) t n 3 k 3 2 1 2 ( x r k 3 x t n 3 ) t n 3 k 3 2 y t n P y r k P x r k P x t n P ( y t n P y r k P ) t n P k P ( x r k P x t n P ) t n P k P 1 2 ( y t n P y r k P ) t n P k P 2 1 2 ( x r k P x t n P ) t n P k P 2 ] X = [ x 0 y 0 v x v y a x a y ] ,   and   F = [ x r k 1 y t n 1 x t n 1 y r k 1 x r k 2 y t n 2 x t n 2 y r k 2 x r k 3 y t n 3 x t n 3 y r k 3 x r k P y t n P x t n P y r k P ]
F is constructed by the locations of the navigation satellites and the receivers, which can be obtained by the mature methods in the navigation field. Since the measurement accuracy of the navigation satellites and the receiver locations in F is high, this equation can be solved by the least squares method. When C T C is full rank, the estimated value of target motion parameters is given by
X = ( C T C ) 1 C T F
Equation (8) can be solved by the square root method, the successive over-relaxation (SOR) method, the Newton iterative method, and so on. When the Newton iterative method is used, the initial iterative value can be obtained by randomly selecting enough equations.
According to the definition of matrix C, when the abscissa of the receiving station and the transmitting station are equal ( x r k i = y t n i ), the element in the second column of the i-th row of the matrix is zero. If that is satisfied for I = 1~P, it means that the system is composed of multiple baselines perpendicular to the horizontal plane and parallel to each other. At this time, parameter estimation cannot be carried out.
Therefore, the estimated result of the target initial position is given by
{ x 0 = [ ( C T C ) 1 C T F ] 1 y 0 = [ ( C T C ) 1 C T F ] 2
For a “one-transmitter and six-receiver” system, the rank of C is five (rank deficit of 1). Thus, it is not possible to use T-RN or TN-R mode for parameter estimation in the FSR network. The minimum configuration of the system is “two-transmitter and three-receiver” or “three-transmitter and two-receiver”. The details are given in the Appendix A.
Then, the initial position estimation result in 2D can be converted into 3D space as follows:
{ x 0 = x 0 cos α x + y 0 cos α y y 0 = x 0 cos β x + y 0 cos β y z 0 = x 0 cos γ x + y 0 cos γ y
where ( α x , β x , γ x ) and ( α y , β y , γ y ) correspond to the direction angles of the x-axis and y-axis direction vectors in the 3D coordinate system, respectively.
From (7), it can be inferred that the factors that affect the estimation accuracy are composed of the crossing time error and station location error. Considering the technical performance of the actual system, the error sources and approximate error magnitude are given in Table 2 (assuming each error is independent of each other and obeys the zero-mean Gaussian distribution).
Where the position of the transmitting/receiving station can be obtained by theodolite or navigation satellite positioning in advance. The errors ex and ey measured by theodolite are about 5~8 m, and eh is about 1 m. The positioning error of GPS is generally about 3 m.
The measurement accuracy of crossing time is analyzed in 23 through incoherent processing. The extracted multiple crossing times are regarded as independent Gaussian random variables. The Cramer Rao lower bound is expressed as
σ n , k = 1 β n , k SNR n , k
where βn,k is the effective bandwidth of the target signal on the baseline, and SNRn,k is the signal-to-noise ratio on the baseline (n, k) after coherent integration within the processing interval. n is the transmitting station number and k is the receiving station number. Considering the FSR network based on Beidou, βn,k = 20.46 MHz, SNRn,k = 13 dB, so σn,k ≈ 13.6 ns. Taking other factors in Table 2 into account, the error standard deviation is enlarged to 1300 µs.

4. Simulation Verification and Error Analysis

In this section, the impacts of the crossing time error and the location errors of the receiver and the transmitter on the parameter estimation accuracy are analyzed through Monte Carlo simulation, and the target positioning capability of the actual system is determined. The system simulation parameters are set as shown in Table 3.
Here, the transmitting/receiving error is a three-dimensional space error. For example, the standard deviation of satellite position error is set to 1000 m, and the random error with a standard deviation of 1000 m is added to the coordinates of the satellite in three directions. Moreover, it is assumed that each error follows a Gaussian distribution with a mean of zero. Monte Carlo simulation is used for statistical error analysis.
The target position estimate refers to the coordinates of the intersection of the target trajectory and the first baseline. In addition, the error of each dimension of the three-dimensional coordinates is defined as follows:
R M S E x 0 = 1 N i = 1 N ( x i x 0 ) 2
R M S E y 0 = 1 N i = 1 N ( y i y 0 ) 2
R M S E z 0 = 1 N i = 1 N ( z i z 0 ) 2
where ( x 0 , y 0 , z 0 ) is the 3D coordinate of the real initial baseline intersection, ( x i , y i , z i ) is the estimated value of three-dimensional coordinates obtained from the i-th Monte Carlo simulation, and N is the number of Monte Carlo simulations.
A.
Influence of measurement error on estimation accuracy
Figure 3 shows the root mean square errors (RMSEs) of the target location, velocity, and acceleration with different crossing time errors, receiver and transmitter station site errors. It can be obtained that all the RMSEs increase with the crossing time error, receiver and transmitter station site errors. The RMSEs of the initial location are shown in Figure 3a,c,e. In accordance with the systematic error in Table 2, it can be seen that the error is less than 100 m when the crossing time error is less than 1.3 ms, the error is less than 280 m when the receiver station site error is less than 5 m, and the error is less than 100 m when the transmitter station site error is less than 1000 m. In the real scene, the current position measurement accuracy of the navigation satellites and receivers can achieve the meter level. Through some special processing, higher accuracy can be obtained. In Figure 3b,d,f, the RMSEs of the velocity and acceleration estimates exceed 200 m/s and 200 m/s2, respectively, when the crossing time error reaches 1.3 ms and the receiver station site error reaches 5 m. Those estimates are larger than the true values.
In general, the RMSE of the target position estimate is in the order of 100 m. Additionally, the estimation error of the coordinate component in the z-direction is much larger than that in the x- and y-directions, which is the main component of the position error. Moreover, the measurement error greatly affects the estimation accuracy of velocity and acceleration. In the case of a large measurement error, the RMSE of acceleration estimate reaches an unacceptable level in the z-direction. In this condition, there is less reference significance to use the velocity and acceleration estimates in this model.
B.
Influence of system structure on estimation accuracy
The system sensitivity analysis consists of the different geometric structures of the FSR network and different heights. The following simulation shows the influence of the spacing between receivers, the spacing between transmitters, and target height on the positioning. The error setting is shown in Table 3, and the minimum system structure of T2–R3 is adopted. In the simulation, the variation range is 100–1000 m for receiver spacing, 10–1000 km for transmitter spacing, and 1000–10,000 m for target height. The results are shown in Figure 4 and Figure 5. As shown in a, the RMSE of positioning decreases gradually with the receiver spacing. When the spacing exceeds 500 m, the error change is no longer obvious. Moreover, a larger receiver spacing is unfavorable for data transmission and synchronization. In Figure 4b, as the transmitter spacing increases, the RMSE of positioning decreases. In addition, the error reduction is no longer obvious when the spacing exceeds 70 km. To obtain a high-precision measurement of the crossing time and small RMSE of positioning, satellites far from each other with better signal quality should be selected as transmitters. In Figure 4c, the positioning error in the z-direction increases gradually with the target height, and we approximate the height as the z coordinate of the initial position of the target; however, the positioning RMSE remains within the acceptable range, and the error is less than 5% when converted into relative error, R M S E z 0 / z 0 = 1 N i = 1 N ( z i z 0 ) 2 / z 0 420.5 / 10 , 000 = 4.205 % , where N is the number of Monte Carlo simulations, and z is the z-axis coordinate of the initial position of the target obtained by each simulation. Therefore, high-precision positioning can be achieved for general target heights.
Figure 5 shows the RMSEs of the target velocity and acceleration. As shown in Figure 5a,b, the RMSEs of velocity and acceleration estimates decline with the receiver spacing and transmitter spacing. The RMSE of acceleration estimate in each direction reduces to approximately 10 m/s2, while the RMSE of velocity estimate in the z-direction exceeds 100 m/s. It is difficult to meet the practical requirements by changing the observation configuration. In Figure 5c, the RMSEs of velocity and acceleration estimates in the x-direction and ydirection hardly change with the increase in target height, while that in the z-direction increases gradually. The RMSE of the acceleration estimate is more than 100 m/s2 in the z-direction, so the estimated results cannot be utilized for the target track.
C.
Influence of target motion on estimation accuracy
The following simulation shows the influence of the target initial velocity and acceleration on the positioning. The error setting is shown in Table 3, and the minimum system structure of T2–R3 is adopted. In the simulation, the variation range is 100–1000 m/s for initial velocity and 10–100 m/s2 for acceleration. Figure 6a shows the RMSE of the coordinate component estimate in each direction with different initial velocities. The RMSE of positioning decreases gradually with the initial velocity increasing. When the initial velocity exceeds 500 m/s, the error change is no longer obvious. Additionally, the RMSE in the z-direction is basically around 260 m. Figure 6b shows the RMSE of positioning with different accelerations. The initial velocity is set to 350 m/s. As we can see, the RMSE of positioning gradually increases with increasing acceleration; however, the increase is limited and the impact can be ignored. Therefore, the target acceleration has little effect on the initial position estimation accuracy.
Summarizing the above simulation analysis, we can obtain the following conclusions:
(1)
The location of a moving target with uniform acceleration can be estimated only using the multiple crossing times in a FSR network. Under the condition of existing systematic error, the RMSE of positioning is less than 100 m.
(2)
The system structure has an impact on target positioning and motion parameter estimation. It is better to set the receiver spacing to about 500 m and choose the transmitters with good signal quality and large spacings.
(3)
Target velocity and acceleration have an impact on target positioning and motion parameter estimation. A high velocity leads to a shorter time for the target to continuously cross multiple baselines and a higher accuracy of target position estimation, while the impact of acceleration can be ignored.
(4)
The velocity and acceleration estimate can be referenced when the systematic error is small. However, it is necessary to use position tracking filtering to obtain higher accuracy when the error is large.

5. Conclusions

In this paper, a motion parameter estimation method only using the crossing times is proposed in the FSR network with GNSS satellite as the radiation source. Additionally, the minimum network configuration is derived theoretically. The effect of error sources and system structure on the parameter estimation accuracy is analyzed. The simulation results indicate that the initial position estimation error is in the order of 100 m for moving target with uniform acceleration under the condition of actual measurement error. This magnitude of error meets the accuracy requirements of early warning detection. The proposed method provides a potential application scenario for passive radar based on GNSS. However, there is a problem that the error values of velocity and acceleration estimates are large. We could obtain higher accuracy by continuously positioning and tracking filtering. In future research, we will further quantitatively analyze the impact of the inaccurate crossing (the target trajectory has a certain distance from the baseline) on the crossing time measurement accuracy. Based on that, the acceptable inaccurate crossing could be determined, which can provide theoretical support for the optimal station layout. A Multi-baseline FSR experimental system is currently being designed and it will be used for the experimental verification of the proposed method later.

Author Contributions

Conceptualization, X.A. and Y.Z.; methodology, X.A. and Y.Z.; validation, Y.Z. writing—original draft preparation, X.A. and Y.Z.; writing—review and editing, Z.X. and F.Z.; supervision, F.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [National Natural Science Foundation of China] grant number [62071475].

Data Availability Statement

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Assuming that a FSR network is composed of only one transmitter and multiple receivers, multiple equations for the motion parameters of the uniformly accelerated motion target can be obtained. Since there are six unknown parameters, at least six equations are required. The simplest system is the “one-transmitter and six-receiver” system, as shown in Figure A1.
Figure A1. FSR system structure diagram of one-transmitter and six-receiver.
Figure A1. FSR system structure diagram of one-transmitter and six-receiver.
Remotesensing 14 01006 g0a1
According to the target motion parameter equation derived in this paper, the positioning matrix C is
C = [ y t y r 1 x r 1 x t 0 0 0 0 y t y r 2 x r 2 x t ( y t y r 2 ) t 12 ( x r 2 x t ) t 12 0.5 ( y t y r 2 ) t 12 2 0.5 ( x r 2 y t ) t 12 2 y t y r 3 x r 3 x t ( y t y r 3 ) t 13 ( x r 3 x t ) t 13 0.5 ( y t y r 3 ) t 13 2 0.5 ( x r 3 y t ) t 13 2 y t y r 4 x r 4 x t ( y t y r 4 ) t 14 ( x r 4 x t ) t 14 0.5 ( y t y r 4 ) t 14 2 0.5 ( x r 4 y t ) t 14 2 y t y r 5 x r 5 x t ( y t y r 5 ) t 15 ( x r 5 x t ) t 15 0.5 ( y t y r 5 ) t 15 2 0.5 ( x r 5 y t ) t 15 2 y t y r 6 x r 6 x t ( y t y r 6 ) t 16 ( x r 6 x t ) t 16 0.5 ( y t y r 6 ) t 16 2 0.5 ( x r 6 y t ) t 16 2 ]
The first four columns of the matrix can be transformed as follows
[ y t y r 1 x r 1 x t 0 0 y t y r 2 x r 2 x t ( y t y r 2 ) t 12 ( x r 2 x t ) t 12 y t y r 3 x r 3 x t ( y t y r 3 ) t 13 ( x r 3 x t ) t 13 y t y r 4 x r 4 x t ( y t y r 4 ) t 14 ( x r 4 x t ) t 14 y t y r 5 x r 5 x t ( y t y r 5 ) t 15 ( x r 5 x t ) t 15 y t y r 6 x r 6 x t ( y t y r 6 ) t 16 ( x r 6 x t ) t 16 ] c 2 × ( y t y 0 ) c 4 × v y [ y t y r 1 ( x r 1 x t ) ( y t y 0 ) 0 0 y t y r 2 ( x r 2 x t ) ( y t y 0 v y t 12 ) ( y t y r 2 ) t 12 ( x r 2 x t ) t 12 y t y r 3 ( x r 3 x t ) ( y t y 0 v y t 13 ) ( y t y r 3 ) t 13 ( x r 3 x t ) t 13 y t y r 4 ( x r 4 x t ) ( y t y 0 v y t 14 ) ( y t y r 4 ) t 14 ( x r 4 x t ) t 14 y t y r 5 ( x r 5 x t ) ( y t y 0 v y t 15 ) ( y t y r 5 ) t 15 ( x r 5 x t ) t 15 y t y r 6 ( x r 6 x t ) ( y t y 0 v y t 16 ) ( y t y r 6 ) t 16 ( x r 6 x t ) t 16 ] c 2 c 3 × v x [ y t y r 1 ( x r 1 x t ) ( y t y 0 ) 0 0 y t y r 2 ( x r 2 x t ) ( y t y 0 v y t 12 ) ( y t y r 2 ) v x t 12 ( y t y r 2 ) t 12 ( x r 2 x t ) t 12 y t y r 3 ( x r 3 x t ) ( y t y 0 v y t 13 ) ( y t y r 3 ) v x t 13 ( y t y r 3 ) t 13 ( x r 3 x t ) t 13 y t y r 4 ( x r 4 x t ) ( y t y 0 v y t 14 ) ( y t y r 4 ) v x t 14 ( y t y r 4 ) t 14 ( x r 4 x t ) t 14 y t y r 5 ( x r 5 x t ) ( y t y 0 v y t 15 ) ( y t y r 5 ) v x t 15 ( y t y r 5 ) t 13 ( x r 5 x t ) t 13 y t y r 6 ( x r 6 x t ) ( y t y 0 v y t 16 ) ( y t y r 6 ) v x t 16 ( y t y r 6 ) t 14 ( x r 6 x t ) t 14 ] c 2 c 1 × ( x 0 x t ) [ y t y r 1 ( x r 1 x t ) ( y t y 0 ) ( y t y r 1 ) ( x 0 x t ) 0 0 y t y r 2 ( x r 2 x t ) ( y t y 0 v y t 12 ) ( y t y r 2 ) ( x 0 x t + v x t 12 ) ( y t y r 2 ) t 12 ( x r 2 x t ) t 12 y t y r 3 ( x r 3 x t ) ( y t y 0 v y t 13 ) ( y t y r 3 ) ( x 0 x t + v x t 13 ) ( y t y r 3 ) t 13 ( x r 3 x t ) t 13 y t y r 4 ( x r 4 x t ) ( y t y 0 v y t 14 ) ( y t y r 4 ) ( x 0 x t + v x t 14 ) ( y t y r 4 ) t 14 ( x r 4 x t ) t 14 y t y r 5 ( x r 5 x t ) ( y t y 0 v y t 15 ) ( y t y r 5 ) ( x 0 x t + v x t 15 ) ( y t y r 5 ) t 15 ( x r 5 x t ) t 15 y t y r 6 ( x r 6 x t ) ( y t y 0 v y t 16 ) ( y t y r 6 ) ( x 0 x t + v x t 16 ) ( y t y r 6 ) t 16 ( x r 6 x t ) t 16 ] = [ y t y r 1 0 0 0 y t y r 2 0 ( y t y r 2 ) t 12 ( x r 2 x t ) t 12 y t y r 3 0 ( y t y r 3 ) t 13 ( x r 3 x t ) t 13 y t y r 4 0 ( y t y r 4 ) t 14 ( x r 4 x t ) t 14 y t y r 5 0 ( y t y r 5 ) t 15 ( x r 5 x t ) t 15 y t y r 6 0 ( y t y r 6 ) t 16 ( x r 6 x t ) t 16 ]
The rank of the positioning matrix is 5, which is less than the number of unknown parameters. Moreover, for the number of receiver stations N > 6, the rank of the positioning matrix is also 5, so neither T-RN nor TN-R mode can be used for motion parameter estimation. The minimum configuration of the FSR network can be determined as a “two-transmitter and three-receiver” or a “three-transmitter and two-receiver” structure.

References

  1. Hamdollahzadeh, M.; Amiri, R.; Behnia, F. Moving Target Localization in Bistatic Forward Scatter Radars: Performance Study and Efficient Estimators. IEEE Trans. Aerosp. Electron. Syst. 2019, 56, 1582–1594. [Google Scholar] [CrossRef]
  2. Gao, C.; Yang, D.; Wang, B.; Zhu, Y. Target moving trajectory estimation by multiple receivers based on GPS forward scattering radar. J. Eng. 2019, 2019, 7755–7760. [Google Scholar] [CrossRef]
  3. Koch, V.; Westphal, R. New approach to a multistatic passive radar sensor for air/space defense. IEEE Aerosp. Electron. Syst. Mag. 1995, 10, 24–32. [Google Scholar] [CrossRef]
  4. Clemente, C.; Soraghan, J.J. GNSS-Based Passive Bistatic Radar for Micro-Doppler Analysis of Helicopter Rotor Blades. IEEE Trans. Aerosp. Electron. Syst. 2014, 50, 491–500. [Google Scholar] [CrossRef]
  5. Ma, H.; Antoniou, M.; Stove, A.G.; Winkel, J.; Cherniakov, M. Maritime Moving Target Localization Using Passive GNSS-Based Multistatic Radar. IEEE Trans. Geosci. Remote Sens. 2018, 56, 4808–4819. [Google Scholar] [CrossRef]
  6. Antoniou, M.; Stove, A.; Tzagkas, D.; Cherniakov, M.; Ma, H. Marine Target Localization with Passive GNSS-Based Multistatic Radar: Experimental Results. In Proceedings of the 2018 International Conference on Radar (RADAR), Brisbane, Australia, 27–31 August 2018; pp. 1–5. [Google Scholar]
  7. Ma, H.; Antoniou, M.; Pastina, D.; Santi, F.; Pieralice, F.; Bucciarelli, M.; Cherniakov, M. Maritime Moving Target Indication Using Passive GNSS-Based Bistatic Radar. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 115–130. [Google Scholar] [CrossRef]
  8. Santi, F.; Pieralice, F.; Pastina, D. Joint Detection and Localization of Vessels at Sea With a GNSS-Based Multistatic Radar. IEEE Trans. Geosci. Remote Sens. 2019, 57, 5894–5913. [Google Scholar] [CrossRef]
  9. Santi, F.; Pastina, D.; Bucciarelli, M. Experimental Demonstration of Ship Target Detection in GNSS-Based Passive Radar Combining Target Motion Compensation and Track-before-Detect Strategies. Sensors 2020, 20, 599. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Gomez-Del-Hoyo, P.; Jarabo-Amores, M.-P.; Mata-Moya, D.; Del-Rey-Maestre, N.; Rosado-Sanz, J. First Approach on Ground Target Detection with GPS based Passive Radar: Experimental Results. In Proceedings of the 2019 Signal Processing Symposium (SPSympo), Krakow, Poland, 17–19 September 2019; pp. 71–75. [Google Scholar]
  11. Hu, C.; Liu, C.; Wang, R.; Chen, L.; Wang, L. Detection and SISAR Imaging of Aircrafts Using GNSS Forward Scatter Radar: Signal Modeling and Experimental Validation. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 2077–2093. [Google Scholar] [CrossRef]
  12. Hu, C.; Liu, C.; Zeng, T. Bistatic forward scattering radar detection and imaging. J. Radars 2016, 5, 229–243. (In Chinese) [Google Scholar]
  13. Zeng, T.; Hu, C.; Cherniakov, M.; Zhuo, C.; Mao, C. Joint parameter estimation and Cramer-Rao bound analysis in ground-based forward scatter radar. EURASIP J. Adv. Signal Process. 2012, 2012, 80. [Google Scholar] [CrossRef] [Green Version]
  14. Gashinova, M.; Daniel, L.; Sizov, V.; Hoare, E.; Cherniakov, M. Phenomenology of Doppler forward scatter radar for surface targets observation. IET Radar Sonar Navig. 2013, 7, 422–432. [Google Scholar] [CrossRef] [Green Version]
  15. Contu, M.; De Luca, A.; Hristov, S.; Daniel, L.; Stove, A.; Gashinova, M.; Cherniakov, M.; Pastina, D.; Lombardo, P.; Baruzzi, A.; et al. Passive Multifrequency Forward-Scatter Radar Measurements of Airborne Targets Using Broadcasting Signals. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 1067–1087. [Google Scholar] [CrossRef]
  16. Daud, N.A.M.; Rashid, N.E.A.; Othman, K.A.; Ahmad, N. Analysis on Radar Cross Section of different target specifications for Forward Scatter Radar (FSR). In Proceedings of the 2014 Fourth International Conference on Digital Information and Communication Technology and its Applications (DICTAP), Bangkok, Thailand, 6–8 May 2014; pp. 353–356. [Google Scholar]
  17. Blyakhman, A.; Burov, V.; Myakinkov, A.; Ryndyk, A. Detection of unmanned aerial vehicles via multi-static forward scattering radar with airborne transmit positions. In Proceedings of the 2014 International Radar Conference, Lille, France, 13–17 October 2014; pp. 1–5. [Google Scholar]
  18. Musa, S.A.; Rsa, R.A.; Sali, A.; Isma’Il, A.; Rashid, N.E.A. DVBS based Forward Scattering Radar for Drone Detection. In Proceedings of the 2019 20th International Radar Symposium (IRS), Ulm, Germany, 26–28 June 2019; pp. 1–8. [Google Scholar]
  19. Radmard, M.; Bayat, S.; Farina, A.; Hajsadeghian, S.; Nayebi, M.M. Satellite-based forward scatter passive radar. In Proceedings of the 2016 17th International Radar Symposium (IRS), Krakow, Poland, 10–12 May 2016; pp. 1–4. [Google Scholar]
  20. Hu, C.; Sizov, V.; Antoniou, M.; Gashinova, M.; Cherniakov, M. Optimal Signal Processing in Ground-Based Forward Scatter Micro Radars. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 3006–3026. [Google Scholar] [CrossRef]
  21. Hu, C.; Zeng, T.; Zhou, C. Accurate three-dimensional tracking method in bistatic forward scatter radar. EURASIP J. Adv. Signal Process. 2013, 2013, 66. [Google Scholar] [CrossRef] [Green Version]
  22. Hu, C.; Wang, L.; Liu, C. SISAR imaging method based on GNSS signal: Theory and experimental results. In Proceedings of the 2016 CIE International Conference on Radar (RADAR), Guangzhou, China, 10–13 October 2016; pp. 1–5. [Google Scholar]
  23. Pastina, D.; Contu, M.; Lombardo, P.; Gashinova, M.; De Luca, A.; Daniel, L.; Cherniakov, M. Target motion estimation via multi-node forward scatter radar system. IET Radar Sonar Navig. 2016, 10, 3–14. [Google Scholar] [CrossRef]
  24. De Luca, A.; Contu, M.; Hristov, S.; Daniel, L.; Gashinova, M.; Cherniakov, M. FSR velocity estimation using spectrogram. In Proceedings of the 2016 17th International Radar Symposium (IRS), Krakow, Poland, 10–12 May 2016; pp. 1–5. [Google Scholar]
  25. Liu, C.; Hu, C.; Wang, R.; Nie, X.; Liu, F. GNSS forward scatter radar detection: Signal processing and experiment. In Proceedings of the 2017 18th International Radar Symposium (IRS), Prague, Czech Republic, 28–30 June 2017; pp. 1–9. [Google Scholar]
  26. Gashinova, M.; Daniel, L.; Cherniakov, M.; Lombardo, P.; Pastina, D.; De Luca, A. Multistatic Forward Scatter Radar for accurate motion parameters estimation of low-observable targets. In Proceedings of the 2014 International Radar Conference, Lille, France, 13–17 October 2014; pp. 1–4. [Google Scholar]
Figure 2. Positioning plane geometry structure.
Figure 2. Positioning plane geometry structure.
Remotesensing 14 01006 g002
Figure 3. Influence of error sources on parameter estimation accuracy. (a) Initial position estimation. (b) Velocity and acceleration estimation. (c) Initial position estimation. (d) Velocity and acceleration estimation. (e) Initial position estimation. (f) Velocity and acceleration estimation.
Figure 3. Influence of error sources on parameter estimation accuracy. (a) Initial position estimation. (b) Velocity and acceleration estimation. (c) Initial position estimation. (d) Velocity and acceleration estimation. (e) Initial position estimation. (f) Velocity and acceleration estimation.
Remotesensing 14 01006 g003
Figure 4. Influence of system structure on the accuracy of the initial position. (a) Influence of receiver spacing. (b) Influence of transmitter spacing. (c) Influence of target height.
Figure 4. Influence of system structure on the accuracy of the initial position. (a) Influence of receiver spacing. (b) Influence of transmitter spacing. (c) Influence of target height.
Remotesensing 14 01006 g004
Figure 5. Influence of system structure on parameter estimation accuracy. (a) Influence of receiver spacing. (b) Influence of transmitter spacing. (c) Influence of target height.
Figure 5. Influence of system structure on parameter estimation accuracy. (a) Influence of receiver spacing. (b) Influence of transmitter spacing. (c) Influence of target height.
Remotesensing 14 01006 g005
Figure 6. Influence of target motion on parameter estimation accuracy. (a) Influence of initial velocity (a = 50 m/s2). (b) Influence of acceleration (v = 350 m/s).
Figure 6. Influence of target motion on parameter estimation accuracy. (a) Influence of initial velocity (a = 50 m/s2). (b) Influence of acceleration (v = 350 m/s).
Remotesensing 14 01006 g006
Table 1. Experimental results of GNSS-based passive radar. (LNA: low noise amplifier; MTI: moving target indicator; CAF: cross ambiguity function; SISAR: shadow inverse synthetic aperture radar).
Table 1. Experimental results of GNSS-based passive radar. (LNA: low noise amplifier; MTI: moving target indicator; CAF: cross ambiguity function; SISAR: shadow inverse synthetic aperture radar).
NumberCountryObjectDetection MethodReceive Gain (dB)Accumulation TimeDetection Range
1Germany [3]Space StationSignal disturbance
Analysis
−5~−320 ms400 km
2Britain [4]Merchant ship (180 × 25 m)Long integration time MTI processing1550 s3 km
3Britain [7]HelicopterMicro-Doppler signature Analysis25 (with LNA)21 ms2.4 km
4Italy [8,9]Merchant ship (150 × 23 m)Long integration time;
Track-before-detect
1510 s1.8 km
5Spain [10]CarCoherent integration of consecutive CAFs35 (with LNA)1s60 m
6China [11]A321/A330SISAR Imaging420 ms1 km
Table 2. Systematic errors.
Table 2. Systematic errors.
Error SourceMain FactorStandard Deviation
Crossing time measurement errorSpace noise environment,
target crossing inaccurately,
RCS fluctuation,
system synchronization error, etc.
<1300 µs
Receiver station site errorStation position measurement error3~5 m
Transmitter station site errorSatellite motion, satellite position measurement error<1000 m
Table 3. Simulation parameter configuration.
Table 3. Simulation parameter configuration.
TargetFlight trajectory (2D) (m)Y = 200x + 7000
Velocity (m/s)350
Acceleration (m/s2)50
ReceiverR1 (m)(0,0,0)
R2 (m)(400,400,10)
R3 (m)(800,800,−10)
TransmitterT1 (km)(100,200,20,200)
T2 (km)(200,300,20,200)
ErrorsThe standard deviation of crossing time error (us)300
The standard deviation of receiver station site error (m)5
The standard deviation of transmitter station site error (m)1000
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ai, X.; Zheng, Y.; Xu, Z.; Zhao, F. Parameter Estimation for Uniformly Accelerating Moving Target in the Forward Scatter Radar Network. Remote Sens. 2022, 14, 1006. https://doi.org/10.3390/rs14041006

AMA Style

Ai X, Zheng Y, Xu Z, Zhao F. Parameter Estimation for Uniformly Accelerating Moving Target in the Forward Scatter Radar Network. Remote Sensing. 2022; 14(4):1006. https://doi.org/10.3390/rs14041006

Chicago/Turabian Style

Ai, Xiaofeng, Yuqing Zheng, Zhiming Xu, and Feng Zhao. 2022. "Parameter Estimation for Uniformly Accelerating Moving Target in the Forward Scatter Radar Network" Remote Sensing 14, no. 4: 1006. https://doi.org/10.3390/rs14041006

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