 Open Access
 Total Downloads : 451
 Authors : A. Sai Raghu, T. Srikanth, Dr. Ch. Sumanth Kumar
 Paper ID : IJERTV3IS042325
 Volume & Issue : Volume 03, Issue 04 (April 2014)
 Published (First Online): 03052014
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Modulation Classification of LPI Radar using Higher Order Statistics(HOS)

Sai Raghu
Department of ECE GIT, GITAM University Visakhapatnam, India

Srikanth
SCD, Signal Processing Wing DLRL – DRDO
Hyderabad, India
Dr. CH. Sumanth Kumar Associate Professor, Department of ECE
GIT, GITAM University Visakhapatnam, India
Abstract The main objective this paper is to detect the LPI radar signal using the Higher Order Statistics (Bispectrum) technique and automatically identify the type of modulation used. Normal spectral analysis techniques are failed to identify LPI signals in the noisy environment. QMFB, Wignerville techniques can detect the LPI signals but they unable to identify the modulation used in the signal and the phase of the signal is not preserved in these techniques. This paper reports the results of HOS technique (Bispectrum) applied to LPI radar signals. Various intra pulse modulations such as Barker signals, P1, P2, P3 and P4 poly phase coded signals are analyzed using HOSA technique to produce 2D signatures of these signals. The results obtained clearly indicate the capability to automatically identify the type of LPI signal even with SNR as low as 3dB.
KeywordsLPI Radar;HOS ; identify

INTRODUCTION
In the modern battlefield, radars face increasingly serious threats from Electronic attack and ARMs. An advanced feature of modern radar systems is, it has the capability to see and not to be seen. Modern radars like Low Probability of Intercept radar has a powerful detection capability while simultaneously itself being not easily detected by electronic reconnaissance equipment.
LPI(Low Probability of Intercept) is the property of radar that because of its very low power, frequency variability,lower sidelobes, wide bandwidth and other design attributes makes its detection difficult by means of passive intercept receivers and ESM (Electronic Support Measure) systems [1]. These combined features prevent the LPI radar detection by modern intercept receivers. Many users of radar today are specifying LPI and LPID as important tactical requirements. Intercept receivers use different strategies to detect the radar based on carrier frequency, antenna scan rate, bandwidth, modulation period etc. Randomly varying of these parameters creates confusion to the intercept receivers [2]. Modern ESM (Electronic Support Measure) system is used for the detection of these LPI radar signals. ESM has to perform the task of detection, parameter identification and classification of these signals in a complex environment of high noise interference and multiple signals.

LPI RADAR WAVEFORM
LPI uses the continuous wave wide bandwidth signals. There are several LPI radar techniques available to the modern radar designer that may be used singly or in various
combinations, depending on the application. Reducing the radars peak effective radiated power (ERP) by using some form of pulse compression technique is the most common LPI radar technique. This spread the radars signal over a wide bandwidth and a period of time. The wideband CW techniques include [1]:

Barker Codes

Phase modulation (polyphase codes)
1. Barker Codes
The binary phasecodes sequence of 0, values that results in equal sidelobe after passage through the matched filter is called Barker code. The longest barker code is of length 13. When longer pulse compression is desired some of the pseudo random code is used. The peak side lobe levels of the linear recursive sequence and of barker codes greater than 5 are lower than 13.2 dB. However the sideobes of the Barker codes can be further lowered by employing a mismatched filter and accepting a slight loss in the peak signaltonoise ratio. Table 1 shows the Barker codes for different lengths.
Table 1 Barker codes
Known Barker Codes of Length (L)
L
Codes
2
+1 1
+1 1
3
+1 +1 1
4
+1 1 +1 +1
+1 1 1 1
5
+1 +1 +1 1 +1
7
+1 +1 +1 1 1 +1 1
11
+1+1 +1 1 1 1 +1 1 1 +1 1
13
+1 +1 +1 +1 +1 1 1 +1 +1 1 +1 1 +1
Phase Shift Keying (PSK) continuous wave (CW) have recently been a topic of active investigation in LPI radar due to their wide bandwidth and low side lobes. Ployphase codes having lengths going into tens or hundreds of subcodes are used to reduce the probability of interception.
A. P1 Code
Different phase codes are generated with different phases but with equal amplitude by varying the synchronous oscillator frequency. By placing the synchronous oscillator at
Assuming that the first sample of I and Q are taken at the leading edge of the waveform, the phases of successive samples taking = (1/B) apart are,
(1)2
the center frequency of the step chirp IF waveform and by sampling the base band waveform at the Nyquist rate, the
=
, where = 1,2, ,
polyphase code called P1 may be obtained. The number of elements in p1 code is2.
If i is the number of the samples in a given frequencyandj is the number of the frequency, the phase of the sample of the frequency is given by the equation
.(3)
D. P4 Code
The P4 code is conceptually derived from the same waveform as the P3 code. In this case the local oscillator is frequency is set equal to 0 + /2 in the IQ detectors. With this frequency, the phases of successive samples taking
,
= 2 1 [ 1 + 1 ] (1)
apart are
Where
= 1,2, ,
and
= 1,2, ,
= 1 2 1 , where = 1,2, ,
An 8phase P1 code of length 64 is considered in this paper.
B. P2 Code
.(4)


HIGHER ORDER STATISTICS
P2 code is derived in the similar way as the P1 code is derived. Starting phase of P2 code is different but the phase increments in each group are same as that of P1 code. The P2 code also has a length or compression ratio of Nc=M*M. P2 code is valid only for even number of N, and each group of code is symmetric about 0 phases. The P2 code is given by
, = 2 1 [2 1 ]
.(2)
Where = 1,2, , and = 1,2, ,
P2 code has frequency symmetry as the P1 code and it also contain property of being palindromic code since the phases are symmetric in the center of the code.
In this paper 8phase P2 code of length 64 is considered. Figure 2 shows the P2 signals.
The LPI waveform data is processed by the Higher Order Spectral Analysis algorithm (Bispectrum) to modulationinformation from the input signal. This includes the use of Bispectrum technique [3]. All these techniques use 256 point FFT. The images (2D plots) produced by the above algorithms are unique for each LPI signal and serves a signature[6]. An experienced operator can also quickly identify the type of modulation code by looking at these plots
A. Bispectrum Etimtion Usisng Direct class of Conventional Methods
Let 1 , 2 , () be the available set of observationsfor Bispectrum estimation. Let us assume that is th sampling frequency and 0 = /0 is the required spacingbetween frequency samples in the Bispectrum domainalong horizontal or vertical directions. Thus 0is the totalnumber of frequency samples.

Segment the data into K segments of M samples each, i.e., N = KM, and subtract the average value of each segment. If necessary, add zeros at each segment to obtain a convenient length M for the FFT.

Assuming that {() (k),k=0,1,2,M1} are the data of segment {i}, generate the DFT coefficients
=
1 exp 2 , (5)
=0
= 0.1 . , 1 and = 1,2, ,
C. P3 Code
Figure 1 P2 signals

In general, = Ã— , where is a positive integer (assumed odd number), i.e., = 2 + 1.In other words size of determines the size of higher order spectrum smoothing over neighboring frequencies. Since is evenand
This code is derived by converting a linear frequency modulation waveform to base band using a local oscillator on one end of the frequency sweep and sampling the I and Q
is odd, we compromise on the value of (closest integer). Estimate the thorder moment spectrum,
, , by frequency domain averaging
1 1
video at the Nyquist rate. If it is assumed that the waveform has a pulse width T in frequency f=f0+kt where k is a constant, the bandwidth B of the signal will be approximately
B=kT.
,
1
1
=
1
+
1
1=
1 1
1=
1
()
+ +
1 + 1
+ + 1
1
+
1
= 1,2, , K
..(6)
coefficient. An experienced operator can detect the signals by looking at the plots, because the plots serves as a signature and unique for each modulation.
,
= 1 .
Figure 6 and 10 shows the Bispectra of P1and b7 signal.
1
1
1 1
1
It is observed that the uniqueness is maintained for all the
1
+ +
1
, = 1,2, .
.(7)
signals even with SNRs of 15 dB and 3 dB to 6 dB and hence facilitates identification [5]. As an example, figures 7,8 shows the calculated Bispectra of p2 signal with SNRs of 0 dB and 6 dB respectively. Figure 9 shows the Bispectra of P4
For third order momentum n=3 i.e., Bispectrum
, = 1 . + . . (8)
modulation.
0
1 2
2
1
2 1 2
B. Modulation Detection
LPI signals use different poly phase modulations. In order

The nth order moment spectrum of the given data is the average over the K pieces
to detect the type of modulation correlation technique is used. A set of reference signals are used for this technique which
, ,
= 1 ( , , )
. . (9)
contains poly phase codes P1, P2, P3 and P4 and Barker
codes of b7, b11 and b13.
1
1
1
=1
1
Figure 4 and 5 shows the b13 and P1 reference signals.
Where 2 ,and = 1,2, , 1
, = 1 , . . (10)
START
1 2
3 1 2
=1


IMPLEMENTATION RESULTS

GUI Implementation
A GUI (Graphical User Interface) has been developed to analyze the LPI signals (at IF level of 160 MHz, Sampled at 500 MHz) of 10uS duration (5000 samples) generated from the Signal Source / samples generated through software.
The following signals are considered

Barker 5, 7, 13 2. P1, P2, P3, P4

For testing the algorithms of HOSA a large database of LPI waveforms is required. The waveform data Polyphase codes P1, P2, P3, and P4 of length 64, carrier frequency 3000, sampling frequency 8000 and length of phase codes considered in this paper is 8.
Figure 2 shows the GUI.
Fig 2 Graphical User Interface
From the output of the Bispectrum the received signal modulation is detected automatically by using the correlation
SET DEFAULT PARAMETERS
fft length,nfft=256, length of Rao gabor window=5,
samples for segment,overlap parameters
CREATE 2D WINDOW
(Rao Gabor)
COMPUTE BISPECTRUM BY ACCUMULATING TRIPLE PRODUCTS H(f1).H(f2).H*(f1+f2)
PERFORM FREQUENCY DOMAIN SMOOTHING
DRAW CONTOUR PLOT OF MAGNITUDE BISPECTRUM
STOP
Fig 3 flow chart forcross correlation computation
Correlation coefficient is computed by cross correlating the Bispectrum signal with the different reference signals i.e., poly phase codes P1, P2, P3 and P4, Barker codes b5, b7, b13 and we get different correlation coefficient values for each modulation. The maximum value of correlation coefficient gives the particular modulation. By considering the maximum value of correlation coefficient it can able to detect the type of modulation used in the LPI signal automatically. Figure 11 shows the automatic detection of the modulation.
Fig 4 reference signal of b13 signal
`
Fig 5 Reference signal of P2 signal
Fig 6 Bispectrum signal of P1 signal
Fig 7 Bispectrum signal of P2 signal
Fig 8Bispectrum signal of P2 signal SNR 6 dB
Fig 9 Bispectrum signal of P4 signal
Fig 10 Bispectrum signal of b7 signal
Fig 11 automatic detection of P7 modulation


CONCLUSION


Bispectrum techniques of LPI signal clearly demarcate one signal from the other even on addition of Gaussian noise. All the poly phase signals have distinct features in bispectrum contour plot. These features are quite visible even in the low SNR of 0dB. This system would be very efficient in detecting
and classifying LPI radar signals automatically and it could be integrated in an ordinary Radar Warning Receivers (RWRs) to make them capable of detecting LPI signals.
Bispectrum gives the non stationary information of the signal. With this feature, it can be applied to find unintentional modulations on pulse to characterize the radar hardware.
ACKNOWLEDGMENT
The authors wish to acknowledge with thanks, the support given by GITAM University Visakhapatnam and DLRL DRDO, Hyderabad, A.P, India, in carrying out this work.
REFERENCES

Philip Pace, Detecting and Classifying LPI Radar, Artech House, Norwood, 2003.

F. Taboada, P.E.Pace et al., Intercept Receiver signal processing techniques to detect low probability of intercept radar signals, Center for Joint Services EW, Naval Postgraduate School, CAK. Elissa, Title of paper if known, unpublished.

Ananthram Swami, Jerry M. Mendel and Chrysostomos L. (Max) Nikias, HigherOrder Spectral Analysis Toolbox for use with Matlab, Mathworks.

C.L. Nikias and A.P.Petropulu., Higher Order Spectral Analysis, PTR Printice Hall, New Jersey, 1993.

Georgios B. Giannakis, Amod V. DandawateHigherorder statistics based input/output system identification and application to noise cancellation, University of Virginia,1989.

L.Anjaneyulu,N.S.Murthy,N.V.S.N.Sarma Identification of LPI Radar Signals by Higher Order Spectra and Neural Network Techniques IEEE, December2008.