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

Authors: Yu Yao; Lenan Wu; Chenmei Zhang
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

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 andK 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.


---

## Body

## 1. 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–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–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 K 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. Section2 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 K 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.


## 2. 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 Figure1, where the transmitted signal st is EBPSK modulated signal with random phase encoding, and the echo signal is interfered by additive white Gaussian noise (AWGN) nt and Weibull (K) distribution clutter hCt.Figure 1
The diagram of JRC system model.The echo signal would be shown as(1)rt=st-t0+2vc0t-t01+hCt+nt,where t0=2r0/c0 is the time corresponding to target initial position, v denotes instantaneous velocity of maneuvering target, and c0 represents velocity of light. Motion expression of point target is(2)vk+1=vk+at,xk+1=xk+vkt+12at2,where v represents the target velocity, x describes the target location, and a denotes the target acceleration. For ease of description, assume that a⌢t is maneuvering acceleration predicted value of previous time. Random maneuvering acceleration is still a first-order time related process:(3)x¨t=a⌢t+at,a˙t=-αat+wt,where at represents zero mean colored acceleration noise, α describes the maneuver time constants, and wt denotes AWGN. And writing up the expression as an equation of motion state, the current statistics adaptive model of maneuvering target would be described as(4)x˙tx¨tx⃛t=I3TI3T22I30I3TI300I3xtx˙tx¨t+00αa⌢t+001wt.For maneuvering target detection, a new scheme with Kalman filter is proposed in the paper, and its flowchart is shown in Figure2, 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.Figure 2
The flowchart of signal processing in JRC system.


## 3. 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 presentH1 and absent H0, and then the echo signal is(5)H0:r=n,r~N0,R0,H1:r=S+n,r~N0,R1,where R1=R0+covS. Under the condition of H0 and H1, the distributions of echo signals are also Gaussian distribution, and probability distributions of echo signals would be written as(6)pr;H0=1πMdet⁡R0exp⁡-rHR0-1r,pr;H1=1πMdet⁡R1exp⁡-rHR1-1r.And the logarithmic likelihood ratio is(7)lr=ln⁡pr;H1pr;H0=ln⁡det⁡R0det⁡R1+rHR0-1r-rHR1-1r=ln⁡det⁡R0det⁡R1+rHR0-1-R0+σx2SSH-1r=ln⁡det⁡R0det⁡R1+rHR0-1S21+σx2SR0-1SH.Equation (7) can be simplified, and CFAR target detection algorithm is(8)Tr=rHR0-1S2≷H0H1γ.Set the false alarm probability Pfa=PTr>γ∣H0. From literature [6], according to Neyman-Pearson lemma, target detection probability would be expressed as(9)Pd=PTr>γ∣H1=Pfaσ02/σ12,where σ12=σ02+SHR0-1S2 and σ02=SHR0-1S. Hence, maximization CFAR target detection probability algorithm is equivalent to(10)d2=SHR0-1S2SHR0-1S=SHR0-1S.From (10), we can get diagonal matrix R0=diag⁡PCf-K/2+⋯+PCfK/2-1, where PCf 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(11)x^i∣i-1=e2aT2/λx^i-1∣i-1,whereλ denotes the wavelength of echo signal.(2)
According to the prediction results, the MSE matrix of echo signal is(12)Pi∣i-1=ΦPi-1∣i-1ΦH+Q,Q=EunuHn.(3)
Define the Kalman filtering gain matrix as(13)Ki=Pi∣i-1HHHPi∣i-1HH+RC-1,RC=Ex^nx^n,H=QS.(4)
Update echo signal:(14)x^i∣i=x^i∣i-1+Kix^i-Hx^i∣i-1.(5)
Update the MSE matrix:(15)Pi∣i=Pi∣i-1-KiHPi∣i-1.(6)
Run steps (2)–(5) iteratively.


## 4. 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 Figure3. 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.Figure 3
Spatial whitening filtering.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 theN pulse returns from the moving target be represented by ar,ejϕar,ej2ϕar,…,ejN-1ϕar, where ar is the spatially whitened output of the first pulse. Also after this transformation, let cnr be the spatially whitened clutter on the nth pulse. After the spatial whitening transformation, we desire to devise a detector to distinguish between the two hypotheses:(16)H1:zr=arsϕ+cr,H0:zr=cr.Let znr be the output of the spatial whitening transformation for the nth pulse. Set zr=z1r,z2r,…,zNrT, cr=c1r,c2r,…,cNrT, and sϕ=(1ejϕej2ϕ⋯ejN-1ϕ)T.Assuming thatar is unknown. The joint PDF of the elements of zr under each hypothesis is given by(17)pzr∣H0=∏ΩrcτNrdetR0exp-1τrzHrR0-1zr,pzr∣H1=∏ΩrcτNrdetR0×exp-1τrzHr-arsϕHR0-1zr-arsϕ,where det⁡ denotes determinant and c is the PDF normalization constant. In this case though, τNr(n=1,2,…,N) are assumed to be independent Gaussian RVs (random variables) from range cell to range cell. The RVsτNr 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 K 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 τr under each hypothesis is given by(18)H0:τMLr=1NzHrR0-1zr,H1:τMLr=1Nzr-arsϕHrR0-1zr-arsϕ.The ML estimate of ar under H1 is given by(19)aMLr=sϕHR0-1zrsϕHR0-1sϕ.It can be shown that an equivalent form for the GLRT is(20)λ=-2N-1∑Ωrln⁡1-sϕHR0-1zr2zHrR0-1zrsϕHR0-1sϕ.And the GLRT takes the form(21)λ≷H0H1T.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(22)z^i∣i-1=e2aT2/λz^i-1∣i-1,whereλ denotes the wavelength of echo signal.(2)
According to the prediction results, the MSE matrix of echo signal is(23)Pi∣i-1=ΦPi-1∣i-1ΦH+Q,Q=EunuHn.(3)
Define the Kalman filter gain matrix as(24)Ki=Pi∣i-1HHHPi∣i-1HH+RC-1,RC=Ez^nz^n,H=QS.(4)
Update echo signal:(25)z^i∣i=z^i∣i-1+Kiz^i-Hz^i∣i-1.(5)
Update the MSE matrix(26)Pi∣i=Pi∣i-1-KiHPi∣i-1.Run steps (2)–(5) iteratively.


## 5. CFAR Detection Probability in Clutter

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


### 5.1. OSGO/OSSO-CFAR Detection inK Distribution

The cumulative distribution function (CDF) ofK distribution would be described as(27)Fx=∫0xfxdx=∫0x2cΓvcx2vKv-1cxdx=∫0cx2ΓvX2vKv-1XdX=1-2Γvcx2vKvcx.The Ordered Statistic with Greatest Option- (OSGO-) CFAR detector is shown in Figure4 (the Ordered Statistic with smallest Option- (OSSO-) CFAR has a similar structure). D denotes detecting units, both sides of the D are reference units (also known as reference sliding windows), and both of the lengths of leading and lagging edge reference sliding windows are M/2. Estimation of clutter power level Z can be derived from reference sliding windows. The function sorts M/2 leading and lagging edge reference units based on OSGO-CFAR and selects the maximum value from X,Y as estimation of clutter power level Z, respectively (it selects the minimum value based on OSSO-CFAR).Figure 4
OSGO-CFAR detector.Threshold valueT 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.Figure4 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 Z. The second step is to multiply this estimation Z by a scaling factor T depending first on the estimation method applied and secondly on the false alarm rate required. The resulting product TZ 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 followK distribution; that is to say, D,X1,…,XM follow K distribution, and the PDF of the kth sample from the overall in uniform background is expressed as(28)
f
k
x
=
k
M
2
k
F
k
-
1
x
f
x
1
-
F
x
M
/
2
-
k
.The CDF of the kth sample is(29)
F
k
x
=
∫
-
∞
x
f
k
x
d
x
=
∫
0
x
k
M
2
k
F
k
-
1
x
f
x
1
-
F
x
M
/
2
-
k
d
x
,where fx and Fx denote PDF and CDF of K distribution, respectively. And for OSGO-CFAR, the PDF of Z=max⁡(X,Y) is expressed as(30)fzz=fXzFYz+fzFz=2fxxFkxx=z.From (28) and (30), false alarm probability of OSGO-CFAR detector would be described as(31)Pfa=PD≥TZ=∫0∞∫T∞fxfzzdxdz=∫0∞1-FTzfzzdz=∫0∞2ΓvTcz2vKvTczckM2k2cz/2vΓvM/2-k+1·Kv-1czKvM/2-kcz1-2Γvcz2vKvczk-1·∫0tckM2k2cx/2vΓvM/2-k+1Kv-1cxKvM/2-kcx1-2Γvcx2vKvcxk-1dxdz,where T denotes threshold factor. From (31), when shape factor v is known, false alarm probability Pfa has nothing to do with scaling factor c. Therefore, OSGO-CFAR detector has the constant false alarm performance in K distribution clutter, and so does OSSO-CFAR detector. The relationship between Pfa and threshold coefficients T is taken by numerical integral method.The optimal values ofk need to be calculated in order to optimize performance of detector in clutter. Based on CFAR detection, the aim of selecting value of k is to make detection probability Pd maximization. We use adaptive detection threshold (ADT) as measurement criteria in this paper. The sum of sinusoidal signal and K distribution calls K-Rice distribution and its PDF is(32)fxH1=∫0∞xrexp⁡-x2+A22rIoAxrfrdr,where r follows Gamma distribution and A is the magnitude of echo signal. The equation applies to unfluctuation target and Swerling 1 model target. For Swerling 2 model target, A follow Rayleigh distribution. A unified formula of detection probability of CFAR detector is(33)Pd=∫0∞∫T∞fxH1fzdxdz.


### 5.2. CFAR Detection in Weibull Distribution

The PDF of random variableZ=Xk is(34)fkz=kNkcBzBc-1exp⁡-zBcN-k+11-exp⁡-zBck-1,z≥0,where k∈1,2,…,M and B and c denote scale and shape parameter, respectively. Shape parameters 1.160≤C≤1.783,(35)Pfa=∫0∞1-F0Tzfkzdz=kMk∫0∞cBzBc-1exp⁡zBcN-k+1+Tc1-exp⁡-zBck-1dz.From (35), the value of factor T has nothing to do with average interference power B of Weibull distribution sample. So it is a CFAR method, and false alarm probability would be simplified as(36)Pfa=kMkΓN-k+1+TCΓkΓN+1+TC.Through the constant false alarm processing, the decision threshold S (factor T) would be calculated for a given false alarm probability Pfa, the length of reference units M and k.Set false alarm probabilityPfa=1×10-6, K distribution clutter average power μ=2, and the length of reference units M=32. Shape parameters v=1.5 and 5. From [10], when k=3/4M~7/8M, target detection performance of OSSO-CFAR detector is optimal. When k=3/8M~7/16M, target detection performance of OSGO-CFAR detector is optimal. For multiple target detection, the detection capability of OSSO/OSGO-CFAR detector would get stronger as the value of k gets smaller.


## 5.1. OSGO/OSSO-CFAR Detection inK Distribution

The cumulative distribution function (CDF) ofK distribution would be described as(27)Fx=∫0xfxdx=∫0x2cΓvcx2vKv-1cxdx=∫0cx2ΓvX2vKv-1XdX=1-2Γvcx2vKvcx.The Ordered Statistic with Greatest Option- (OSGO-) CFAR detector is shown in Figure4 (the Ordered Statistic with smallest Option- (OSSO-) CFAR has a similar structure). D denotes detecting units, both sides of the D are reference units (also known as reference sliding windows), and both of the lengths of leading and lagging edge reference sliding windows are M/2. Estimation of clutter power level Z can be derived from reference sliding windows. The function sorts M/2 leading and lagging edge reference units based on OSGO-CFAR and selects the maximum value from X,Y as estimation of clutter power level Z, respectively (it selects the minimum value based on OSSO-CFAR).Figure 4
OSGO-CFAR detector.Threshold valueT 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.Figure4 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 Z. The second step is to multiply this estimation Z by a scaling factor T depending first on the estimation method applied and secondly on the false alarm rate required. The resulting product TZ 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 followK distribution; that is to say, D,X1,…,XM follow K distribution, and the PDF of the kth sample from the overall in uniform background is expressed as(28)
f
k
x
=
k
M
2
k
F
k
-
1
x
f
x
1
-
F
x
M
/
2
-
k
.The CDF of the kth sample is(29)
F
k
x
=
∫
-
∞
x
f
k
x
d
x
=
∫
0
x
k
M
2
k
F
k
-
1
x
f
x
1
-
F
x
M
/
2
-
k
d
x
,where fx and Fx denote PDF and CDF of K distribution, respectively. And for OSGO-CFAR, the PDF of Z=max⁡(X,Y) is expressed as(30)fzz=fXzFYz+fzFz=2fxxFkxx=z.From (28) and (30), false alarm probability of OSGO-CFAR detector would be described as(31)Pfa=PD≥TZ=∫0∞∫T∞fxfzzdxdz=∫0∞1-FTzfzzdz=∫0∞2ΓvTcz2vKvTczckM2k2cz/2vΓvM/2-k+1·Kv-1czKvM/2-kcz1-2Γvcz2vKvczk-1·∫0tckM2k2cx/2vΓvM/2-k+1Kv-1cxKvM/2-kcx1-2Γvcx2vKvcxk-1dxdz,where T denotes threshold factor. From (31), when shape factor v is known, false alarm probability Pfa has nothing to do with scaling factor c. Therefore, OSGO-CFAR detector has the constant false alarm performance in K distribution clutter, and so does OSSO-CFAR detector. The relationship between Pfa and threshold coefficients T is taken by numerical integral method.The optimal values ofk need to be calculated in order to optimize performance of detector in clutter. Based on CFAR detection, the aim of selecting value of k is to make detection probability Pd maximization. We use adaptive detection threshold (ADT) as measurement criteria in this paper. The sum of sinusoidal signal and K distribution calls K-Rice distribution and its PDF is(32)fxH1=∫0∞xrexp⁡-x2+A22rIoAxrfrdr,where r follows Gamma distribution and A is the magnitude of echo signal. The equation applies to unfluctuation target and Swerling 1 model target. For Swerling 2 model target, A follow Rayleigh distribution. A unified formula of detection probability of CFAR detector is(33)Pd=∫0∞∫T∞fxH1fzdxdz.


## 5.2. CFAR Detection in Weibull Distribution

The PDF of random variableZ=Xk is(34)fkz=kNkcBzBc-1exp⁡-zBcN-k+11-exp⁡-zBck-1,z≥0,where k∈1,2,…,M and B and c denote scale and shape parameter, respectively. Shape parameters 1.160≤C≤1.783,(35)Pfa=∫0∞1-F0Tzfkzdz=kMk∫0∞cBzBc-1exp⁡zBcN-k+1+Tc1-exp⁡-zBck-1dz.From (35), the value of factor T has nothing to do with average interference power B of Weibull distribution sample. So it is a CFAR method, and false alarm probability would be simplified as(36)Pfa=kMkΓN-k+1+TCΓkΓN+1+TC.Through the constant false alarm processing, the decision threshold S (factor T) would be calculated for a given false alarm probability Pfa, the length of reference units M and k.Set false alarm probabilityPfa=1×10-6, K distribution clutter average power μ=2, and the length of reference units M=32. Shape parameters v=1.5 and 5. From [10], when k=3/4M~7/8M, target detection performance of OSSO-CFAR detector is optimal. When k=3/8M~7/16M, target detection performance of OSGO-CFAR detector is optimal. For multiple target detection, the detection capability of OSSO/OSGO-CFAR detector would get stronger as the value of k gets smaller.


## 6. Simulation Results

In this section, the performance of the JRC system with Kalman filtering is simulated. In Table1, a summary of the most important parameters of the simulation model is provided.Table 1
Simulation parameters.
Parameter value
Carrier frequency:f1=10 GHz
Intermediate frequency:f=100 MHz
Sample frequency:fs=10 GHz
EBPSK modulation parameters:N=7K=2
False alarm probability:Pfa=10-6
Target detection probability:Pd=0.9
The number of pulses:N=10
Pulse width:T=200 μs
Simulated Orbit 1:
Initial velocity: 450 m/s; accelerated velocity: 40 m/s2
Simulated Orbit 2:
Initial velocity: 450 m/s; centripetal acceleration: 20 m/s2Simulated 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/s2. 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 x-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/s2. 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 Figures5–8. The results are obtained from 4000 Monte Carlo simulations under different SNRs.Figure 5
Simulated Orbit 1 target detection.Figure 6
Simulated Orbit 2 target detection.Figure 7
Target detection in Weibull distribution clutter.Figure 8
Target detection inK distribution clutter.Simulation 1. Simulated Orbit 1 target detection is shown in Figure 5.Figure5 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 probabilityPd 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 probabilityPd 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 Figure6.Figure6 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 probabilityPd 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 M=32, and shape parameter c=1.8. Take k=3M/4 in the following simulation.Figure7 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, whenPd=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, whenPd=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 (
K
Distribution Clutter). A K 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 M=32, and shape parameter c=1.5. Take k=3M/8 in the following simulation.Figure8 shows CFAR target detection probability comparisons of the JRC systems with or without Kalman filtering in K distribution clutter channel.For Simulated Orbit 1 maneuvering target, whenPd=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, whenPd=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 inK distribution clutter is slightly worse than that in Weibull distribution clutter.


## 7. 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 modulated echo pulse. Furthermore, CFAR target detection probability of the new JRC system in Weibull andK distribution clutter is deduced. The theoretical and simulation results presented confirm the improvement of maneuvering target detection performance in the Gaussian and sea clutter channel via Kalman filtering.



---

*Source: 471305-2015-12-28.xml*