Maneuvering Target Detection Based on JRC System in Gaussian and Non-Gaussian Clutter

Aimed at the problem of detecting maneuvering targets in the Gaussian and sea clutter environments and based on the established motion state model, this paper proposed a new scheme that uses a joint radar-communication (JRC) system with Kalman filter to accurately detect the target with the generalized likelihood ratio test (GLRT) theory and a constant false alarm rate (CFAR) based threshold. Also, the theoretical threshold and probability function of GLRT target detection based on CFAR were given. Moreover, target detection probability of the new JRC system in Weibull and K distribution clutter is deduced. In addition to theoretical considerations, simulations and measurement results of the new JRC systems demonstrate excellent detection performance for maneuvering targets in the Weibull and K distribution channel.


Introduction
In recent years, maneuvering target detection problem has attracted substantial interest, which is very meaningful for modern radars to improve their detection performance [1][2][3][4].The complexity of target detection is not only related to the target itself, but also influenced by the electromagnetic environment and clutter background.There are various kinds of weak moving targets: (1) small size target with weak radar returns; (2) low-observable stealth target with low radar cross section (RCS); (3) remote target or target in strong clutter with low signal-to clutter/noise ratio (SCR/SNR); and (4) high speed or highly mobile target with limited available energy.The above four kinds of targets can be summarized as "far-range, low-observable, and highly maneuvering" targets.They share a lot in common that, both in time and in frequency domain, the SCR/SNR is too low to compete with clutter and electronic interference, resulting in poor detection performance [5,6].
The joint radar-communications (JRC) system in the literature [7][8][9] showed excellent constant false alarm rate (CFAR) target detection performance when performing radar measurements for uniform motion targets.A simple amplitude threshold detector in previous JRC system can demodulate the EBPSK modulated pulse echoes, resulting in a simple receiver structure.The impacting filter [10,11] can be digitally implemented, which is beneficial for chip integration.However, target detection performance of the previous JRC system becomes unsatisfactory for maneuvering target detection in non-Gaussian channel.
As the radar range resolution increases, the background clutter may no longer be modeled accurately as a Gaussian random variable (RV).At the higher range resolution, the radar system receives target-like spikes that give rise to non-Gaussian observations.In many cases, it has been found that the spiky clutter returns can be suitably modeled by a multidimensional Gaussian mixture distribution or more specifically a spherically invariant random vector (SIRV) [12].Such well-known clutter distributions are known as the  and Weibull [13,14].In [15] it was shown that from a reasonable and simple phenomenological model of clutter, the clutter returns from multiple pulses of a given range can be modeled as a SIRV distribution.
In this paper, the generalized likelihood ratio test (GLRT) target detection algorithm based on CFAR is used to detect the presence of target.Aimed at the problem of detecting maneuvering targets with high speed in sea clutter and based on the established motion state model, the JRC system with Kalman filter accurately detects the target with GLRT theory and sets the detection threshold.If the power of echo signal after Kalman filtering is greater than the threshold, then the radar system detects the presence of target, otherwise not.The rest of this paper is organized as follows.Section 2 introduces the model of EBPSK modulated JRC system.Section 3 describes target detection based on GLRT in Gaussian clutter.Section 4 illustrates target detection based on GLRT in Weibull and  distribution clutter via Kalman filtering.In Section 5, the CFAR target detection probability of the JRC system in sea clutter is deduced.Some indicative simulation results and performance analysis are presented in Section 6.And finally, Section 7 gives the conclusion of the paper.

The Model of EBPSK Modulated JRC System
The EBPSK modulated JRC system model for detection and estimation of maneuvering target in non-Gaussian clutter is shown in Figure 1, where the transmitted signal () is EBPSK modulated signal with random phase encoding, and the echo signal is interfered by additive white Gaussian noise (AWGN) () and Weibull () distribution clutter ℎ  ().
The echo signal would be shown as where  0 = 2 0 / 0 is the time corresponding to target initial position, V denotes instantaneous velocity of maneuvering target, and  0 represents velocity of light.Motion expression of point target is where V represents the target velocity,  describes the target location, and  denotes the target acceleration.For ease of description, assume that  () is maneuvering acceleration predicted value of previous time.Random maneuvering acceleration is still a first-order time related process: where () represents zero mean colored acceleration noise,  describes the maneuver time constants, and () denotes AWGN.And writing up the expression as an equation

𝑥 (𝑡)
For maneuvering target detection, a new scheme with Kalman filter is proposed in the paper, and its flowchart is shown in Figure 2, where echo signal preprocess and target detection in the new JRC system are as follows: (1) A state model of the target is established, which is in accord with the actual maneuvering situation and convenient for digital processing.(2) Firstly, the JRC system sends EBPSK modulated pulse with random phase coded.(3) The JRC system preliminarily detects the target with GLRT theory and sets the detection threshold based on CFAR.If the power of echo signal is greater than the threshold, then the radar system detects the presence of target, otherwise not.(4) When the GLRT detection shows that target is present, based on the location (range), velocity, and acceleration correlation of maneuvering target, a method to process the echo via Kalman filtering based on their phase information correlation is proposed.(5) Based on the established motion state model, the JRC system with Kalman filter accurately detects the target with GLRT theory and sets the detection threshold based on CFAR.If the power of echo signal after Kalman filtering is greater than the threshold, then the radar system detects the presence of target, otherwise not.

Target Detection Based on GLRT in Gaussian
CFAR target detection algorithm based on Neyman-Pearson criterion is used as the rule of EBPSK modulated JRC system in this paper.Assume that the target is present  1 and absent  0 , and then the echo signal is where . Under the condition of  0 and  1 , the distributions of echo signals are also Gaussian distribution, and probability distributions of echo signals would be written as And the logarithmic likelihood ratio is Equation ( 7) can be simplified, and CFAR target detection algorithm is Set the false alarm probability  fa = ((r) >  |  0 ).From literature [6], according to Neyman-Pearson lemma, target detection probability would be expressed as where Hence, maximization CFAR target detection probability algorithm is equivalent to From (10), we can get diagonal matrix R 0 = diag(  ( −/2 ) + ⋅ ⋅ ⋅+  ( /2−1 )), where   () denotes power spectral density function of clutter.
With the time and velocity correlation of the maneuvering target, the EBPSK modulated JRC system uses Kalman filtering to process echo signals when the GLRT detection shows the presence of target.And the estimation performance can be improved by taking the advantage of prediction and estimation at the same time.The Kalman filtering process can be described as follows: (1) According to the phase correlation of echo signal, the prediction of echo signal is where  denotes the wavelength of echo signal.

Target Detection in Non-Gaussian Clutter
We assume that the individual pulses are passed through a spatial whitening filter and thereafter detection takes places, depicted in Figure 3. Hence, the clutter component of the samples in range is statistically decorrelated.The spatial whitening filter essentially removes the correlation from range sample to range sample due to overlapping of the waveform returns from the individual scatterers.After spatial whitening the texture component on each range sample is an independent RV with respect to the other range samples.The reason behind spatially whitening the input clutter is that we can conveniently write down its space dependent multidimensional probability density function (PDF).Because all of the RVs are conditionally Gaussian, spatially whitening the clutter allows us to write down the joint space PDF conditioned on the texture RV as a product of the PDFs conditioned on the texture RV at each range.
It is pointed out that matched filtering in range on receiver is not used because the target and clutter have the same power spectrum.
Let the  pulse returns from the moving target be represented by ((),   (),  2 (), . . .,  (−1) ()), where () is the spatially whitened output of the first pulse.Also after this transformation, let   () be the spatially whitened clutter on the th pulse.After the spatial whitening transformation, we desire to devise a detector to distinguish between the two hypotheses: Let z  () be the output of the spatial whitening transformation for the th pulse.Set z() = ( 1 (),  2 (), . . .,   ())  , c() = ( 1 (),  2 (), . . .,   ())  , and Assuming that () is unknown.The joint PDF of the elements of z() under each hypothesis is given by where det() denotes determinant and  is the PDF normalization constant.In this case though,   () ( = 1, 2, . . ., ) are assumed to be independent Gaussian RVs (random variables) from range cell to range cell.The RVs   () are used to model the large power fluctuations associated with the various clutter levels found in different range cells.Standard clutter models of amplitude distribution, such as  and Weibull, are included in the class of Rayleigh mixture distributions which are related to the class of SIRV (spherically invariant random vector) distribution.The maximum likelihood (ML) estimate of () under each hypothesis is given by The ML estimate of () under  1 is given by It can be shown that an equivalent form for the GLRT is And the GLRT takes the form Thus, for a chosen false alarm probability the threshold is easily computable.
The JRC system with Kalman filter accurately detects the target with GLRT theory and sets the detection threshold based on CFAR.The processing of Kalman filtering would be described as follows: (1) According to the prediction results, pulse echo is where  denotes the wavelength of echo signal.
(2) According to the prediction results, the MSE matrix of echo signal is (3) Define the Kalman filter gain matrix as (4) Update echo signal: (5) Update the MSE matrix Run steps (2)-( 5) iteratively.

CFAR Detection Probability in Clutter
Compared to Rayleigh distribution, the envelopes of the PDFs of sea clutter have a longer tail, such as Weibull and  distribution.From literature [16], the PDFs of Weibull and  distribution have two unknown parameters.We discuss presentative OSGO/OSSO-CFAR detector in Weibull and  distribution clutter in this section.

OSGO/OSSO-CFAR Detection in 𝐾 Distribution.
The cumulative distribution function (CDF) of  distribution would be described as The Ordered Statistic with Greatest Option-(OSGO-) CFAR detector is shown in Figure 4 (the Ordered Statistic with smallest Option-(OSSO-) CFAR has a similar structure). denotes detecting units, both sides of the  are reference units (also known as reference sliding windows), and both of the lengths of leading and lagging edge reference sliding windows are /2.Estimation of clutter power level  can be derived from reference sliding windows.The function sorts /2 leading and lagging edge reference units based on OSGO-CFAR and selects the maximum value from ,  as estimation of clutter power level , respectively (it selects the minimum value based on OSSO-CFAR).
Threshold value  can be calculated according to the corresponding constant false alarm algorithm and set value of false alarm probability.Estimation of clutter power level has self-adaption when the clutter intensity changes.
Figure 4 shows a certain reference window and the signal processing structure for target detection as a flowchart.The first step is to measure the mean clutter power level .The second step is to multiply this estimation  by a scaling factor  depending first on the estimation method applied and secondly on the false alarm rate required.The resulting product  is directly used as the threshold value.An estimation procedure is proposed which derives the clutter power estimation from the so-called ordered statistic.The two clutter power estimators are then combined into one single value either by further averaging or by maximum selection.
When maneuvering target is absent, background clutters follow  distribution; that is to say, , 1, . . .,  follow  distribution, and the PDF of the th sample from the overall in uniform background is expressed as The CDF of the th sample is where () and () denote PDF and CDF of  distribution, respectively.And for OSGO-CFAR, the PDF of  = max(, ) is expressed as = .(30) From ( 28) and (30), false alarm probability of OSGO-CFAR detector would be described as where  denotes threshold factor.From (31), when shape factor V is known, false alarm probability  fa has nothing to do with scaling factor . Therefore, OSGO-CFAR detector has the constant false alarm performance in  distribution clutter, and so does OSSO-CFAR detector.The relationship between  fa and threshold coefficients  is taken by numerical integral method.The optimal values of  need to be calculated in order to optimize performance of detector in clutter.Based on CFAR detection, the aim of selecting value of  is to make detection probability   maximization.We use adaptive detection threshold (ADT) as measurement criteria in this paper.The sum of sinusoidal signal and  distribution calls -Rice distribution and its PDF is where  follows Gamma distribution and  is the magnitude of echo signal.The equation applies to unfluctuation target and Swerling 1 model target.For Swerling 2 model target,  follow Rayleigh distribution.A unified formula of detection probability of CFAR detector is )  ()  . (33)

CFAR Detection in Weibull
Distribution.The PDF of random variable  =   is where  ∈ {1, 2, . . ., } and  and  denote scale and shape parameter, respectively.Shape parameters 1.160 ≤  ≤ 1.783, From (35), the value of factor  has nothing to do with average interference power  of Weibull distribution sample.So it is a CFAR method, and false alarm probability would be simplified as Through the constant false alarm processing, the decision threshold  (factor ) would be calculated for a given false alarm probability  fa , the length of reference units  and .

Simulation Results
In this section, the performance of the JRC system with Kalman filtering is simulated.In Table 1, a summary of the most important parameters of the simulation model is provided.
Simulated Orbit 1.At the beginning, the target bears a constant velocity, and its initial value is 450 m/s.The target is moving in uniformly accelerated motion between 51st and 100th scanning intervals, and accelerated velocity is 40 m/s 2 .And the target is moving in uniform motion again between 101st and 150th scanning intervals.Simulated Orbit 2. At first, the target is moving in constant velocity along the -axis, the initial velocity is 450 m/s, the target is moving in circular accelerated motion from the 51st scan cycle, and centripetal acceleration is 20 m/s 2 .And the target is moving in uniform motion again after the end of a period.In the following simulations, the relations between CFAR target detection probabilities and the SNR are shown in Figures 5-8.The results are obtained from 4000 Monte Carlo simulations under different SNRs.
Simulation 1. Simulated Orbit 1 target detection is shown in Figure 5.
Figure 5 shows Simulated Orbit 1 target detection probability comparisons for EBPSK modulated pulse echo signal with or without Kalman filtering at the receiver.
For Simulated Orbit 1 maneuvering target, when the detection probability   is 90%, the SNR performance of the JRC system with Kalman filtering would be improved by approximately 2.4 dB as compared with the proposed system in AWGN channel.For uniform motion targets, when the detection probability   is 90%, the SNR performance of JRC system did not improve after Kalman filtering as compared with the proposed system in AWGN channel.
Simulated Orbit 2 target detection simulation is shown in Figure 6.
Figure 6 shows Simulated Orbit 2 target detection probability comparisons for EBPSK modulated pulse echo signal with or without Kalman filtering at the receiver.
For Simulated Orbit 2 maneuvering target, when the detection probability   is 90%, the SNR performance of the JRC system with Kalman filtering would be improved by approximately 2 dB as compared with the proposed system in AWGN channel.And detection performance of maneuvering target 2 is better than maneuvering target 1.

Simulation 2 (Weibull Distribution Clutter).
A Weibull distribution clutter model with 32 ranging units is established.Low-power is located between range units 1 and 16.Clutter edger is located in range unit 17.High-power clutter is located between range units 17 and 32.The length of reference units  = 32, and shape parameter  = 1.8.Take  = 3/4 in the following simulation.
Figure 7 shows CFAR target detection probability comparisons of the JRC systems with or without Kalman filtering in Weibull distribution clutter channel.
For Simulated Orbit 1 maneuvering target, when   = 90%, the SNR performance of the JRC system with Kalman filtering would be improved by approximately 2.5 dB as compared with the proposed system.
For uniform motion targets, when   = 90%, the SNR performance of the JRC system with Kalman filtering would be improved by approximately 0.5 dB as compared with the proposed system.Target detection performance of the JRC system in Weibull distribution clutter is much worse than in AWGN.

Simulation 3 (𝐾 Distribution Clutter).
A  distribution clutter model with 32 ranging units is established.Low-power is located between range units 1 and 16.Clutter edger is located in range unit 17.High-power clutter is located between range units 17 and 32.The length of reference units  = 32, and shape parameter  = 1.5.Take  = 3/8 in the following simulation.
Figure 8 shows CFAR target detection probability comparisons of the JRC systems with or without Kalman filtering in  distribution clutter channel.
For Simulated Orbit 1 maneuvering target, when   = 90%, the SNR performance of the JRC system with Kalman filtering would be improved by approximately 3 dB as compared with the proposed system.
For uniform motion targets, when   = 90%, the SNR performance of the JRC system with Kalman filtering would be improved by approximately 0.6 dB as compared with the proposed system.
Therefore, target detection performance of the JRC system in  distribution clutter is slightly worse than that in Weibull distribution clutter.

Conclusions
The detection and estimation problems of maneuvering target are considered, and a new JRC working scheme with Kalman filter is proposed in this paper.The established motion state model of the maneuvering target is given.The generalized likelihood ratio test (GLRT) target detection algorithm based on constant false alarm rate (CFAR) is used to detect the presence of target.When target appears, the Kalman filtering approach is used to estimate EBPSK

Figure 2 :
Figure 2: The flowchart of signal processing in JRC system.