International Journal of Antennas and Propagation

Volume 2015 (2015), Article ID 802471, 8 pages

http://dx.doi.org/10.1155/2015/802471

## A New Jammer Suppression Method in MIMO Radar Using Matrix Pencil Method and Generalized Likelihood Ratio Test

^{1}Department of Electronic Convergence Engineering, Kwangwoon University, Seoul 139-701, Republic of Korea^{2}Agency for Defense Development, Daejeon 305-600, Republic of Korea

Received 13 October 2014; Accepted 5 February 2015

Academic Editor: Luciano Tarricone

Copyright © 2015 Kang-in Lee et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We propose a new algorithm to suppress the jammer signals and estimate the direction of arrival (DOA) of the signal of interest (SOI) for collocated MIMO radar by using the matrix pencil method (MPM) and the generalized likelihood ratio test (GLRT). The conventional GLRT divides the visible region into small angle samples, suppresses the jammer signals at each angle sample, and then estimates the DOA of the SOI. In the proposed algorithm, we extract the eigenvalues of received signals regardless of the SOI and jammer by using the MPM, which contain the information of the DOA of SOIs or jammers. Then, in order to suppress the jammers, we apply the GLRT to the extracted DOAs instead of to the entire visible region. By applying the MPM again to the received signals in which the jammer signals are suppressed, we can estimate the DOAs of the SOI. Since the proposed algorithm does not depend on the number of angle samples, it shows fast and accurate results regardless of the angle resolution. In order to verify the proposed algorithm, we compared the results with the results of the conventional GLRT and show the computing time.

#### 1. Introduction

In recent times, multi-input multioutput (MIMO) radar technologies have been of interest in many areas of military and civilians, and it has been verified that MIMO radar has a lot of potential advantages over conventional phased-array radar [1–8]. In military area, the MIMO radar has an important role of detecting and tracking the low-observable targets and performing high resolution imaging. In also commercial areas such as the ITS (intelligent transport system), the collision avoidance system of car, and the weather broadcasting radar, the MIMO radar has attracted considerable attention.

A MIMO radar system is composed of transmitters and receivers, which can be collocated or distributed spatially [3]. Typically, it is assumed that at each receiver a technique exists for unambiguously separating the reflected signals of interest (SOI) from each transmitter, by utilizing orthogonal waveforms and a matched filter bank.

In [5, 9], many researchers have shown that, by utilizing orthogonal waveforms, a MIMO radar system with spatially diverse transmitters and receivers, can provide advantages in target detection and parameter estimation compared to a traditional phased array system [10]. For collocated transmit and receive antennas, the MIMO radar has been shown to get higher resolution than that of a phased array radar using the same number of physical antenna elements. This is because the MIMO radar has more degrees of freedom (DOF) than a phased array system with a single transmitting element [3, 4]. In general, the received signal of the radar is composed of the signal of interest (SOI), interference or jammer, and noise. In order to extract the SOI from the received signal, the jammer signal should be suppressed efficiently and various suppression algorithms have been studied [11–15].

In this paper, we propose a new jammer suppression algorithm for the MIMO radar by combining the MPM and the GLRT. The conventional GLRT divides the entire angle region or the visible region into an angle sample , where is the number of angle samples [1, 16–20]. It then performs suppression of the jammer signals or estimation of the DOAs of the SOI, by comparison with the given threshold value at each angle sample. If the angle of the DOA is an integer multiple of the angle sample (which in this paper is called the matched angle), the performance of the conventional GLRT shows good results; otherwise, the results may be of relatively poor quality. Although we can solve this problem by using a smaller angle sample, it increases the computing time.

In this paper, we first extract the number and DOA of signals above a threshold by using MPM. The extracted signals can be SOIs or jammers. Then, the jammer signals can be suppressed, by applying the GLRT to the extracted DOAs. In general, since the number of the extracted signals is much smaller than , the proposed algorithm can reduce the number of operations of the GLRT compared with the conventional algorithm. The DOAs of the SOI can then be estimated by applying the MPM again to the output signals of the virtual array of the MIMO radar.

In this paper, we verify the performance of the proposed algorithm under various environments. In the first case, we test it according to the level of the jammer-to-noise ratio (JNR). In the other case, when the DOAs of the signals match the angle sample or do not, we compare the results of the proposed algorithm with those of the conventional one. Since the proposed algorithm estimates the DOAs of all the signals regardless of the angle sample, it shows good performance even when the DOA is not matched to the angle sample.

Also, we derive the computational complexity of the proposed algorithm and compare it with that of the conventional one.

#### 2. Signal Model

In this paper, we assume the following:(i)the transmitted probing signal is assumed to be a narrowband signal [1, 21];(ii)there is nondispersive propagation;(iii)transmitters and receivers are collocated;(iv)the targets are represented by nonmoving points.

In general, MIMO radar is composed of transmitters with orthogonal waveforms and receivers with matched filters. At the target location, the transmitted signals are given by where, is the discrete-time baseband signal sequence from the th transmitter, is the azimuth angle of the target, is the carrier angular frequency of the radar, is the propagation time from the th transmitter to the target, and is the number of data samples of the transmitted signal.

Also, means the Hermitian transpose, the transmitted sequence vector, and the steering vector of the transmitting antenna, which are represented by where means the transpose of a vector. At the receiver, the received signal sequence vector is given by where is the propagation time from the target to the th receiver, is the th target location, is the complex amplitude proportional to RCS of the th target, is the spatially and temporally white circularly symmetric complex Gaussian noise sequence vector with mean zero [1], is the number of targets of interest, and represents the jamming signal vector. Also, is the steering vector of the receiving antenna and denotes the complex conjugate of a vector.

In (3), the received jamming signal can be written as where is the th jammer location, the complex amplitude of the th jammer, and the discrete-time baseband random signal of the th jammer. denotes the number of jammer signals.

#### 3. Jammer Suppression and DOA Estimation Using MPM and GLRT

The matrix pencil method (MPM) approximates a given data set or function with the sum of the complex exponential functions. Applying the MPM to a phased array, it can extract the DOAs and the number of SOIs [22–24].

Consider a uniformly spaced linear array (ULA) antenna. Then, the received signal at the th antenna element of the ULA can be represented as where is the number of SOIs, the wavelength, the distance between the adjacent antenna elements, the complex amplitude, and the DOA of the th SOI. In this paper, is called a pole. In a phased array, the pole vector corresponds to the steering vector that includes the DOAs of received signals. The MPM estimates the pole vector by applying singular value decomposition (SVD) to the Hankel matrix of (7).

In this paper, by using the MPM, we extract the DOAs of signals greater than the threshold from the received signal vector in the MIMO radar. Note that the extracted signals may not include a few jammers or SOIs according to the jammer-to-noise ratio (JNR) level.

When (5) is substituted into (3), (3) can be rewritten as where and mean the distance between adjacent antennas of the receiver and the transmitter, respectively. Note that in the ULA, and . From (9), we can derive the following Hankel matrix at each data sample as [22]: In (10), is the received signal at the th receiver with the th data sample and means the pencil parameter and its size is about 1/3 the size of (9). Using the block matrix of (10), we construct a new matrix as [23]

Applying the MPM to (11), we extract the pole of the received signals above the threshold value as where is the th DOA of the estimated signals and is the number of the extracted poles. Using extracted pole , the received signals can be reconstructed as

In general, ; and in the case of high JNR, converges to . In (13), the complex amplitude vector can be estimated by using the least square method (LSM) as where In (15), means the pseudo inverse of a matrix.

Applying the GLRT at the estimated DOA , we determine whether the signal from is a jammer or an SOI. In this paper, in order to suppress the jammer signal, we use the following equation as the GLRT [5]:

Equation (17) represents the probability of the existence of the SOI at the testing angle . While in the conventional GLRT method, is the integer multiplication of the angle sample , the proposed method makes it possible to be an arbitrary angle estimated by applying the MPM to (11). In (17), denotes the sample covariance matrix of the received signal and the residual covariance matrix as follows [1, 5, 25]:

If the extracted angle satisfies , it is considered as the DOA of a jammer. Note that when the JNR is very high, the first step MPM may fail to extract the DOA of the SOI due to its eigenvalue being lower than the threshold . In this case, the proposed method can reconstruct the DOA of the SOI, by applying the MPM again to the received signal vector in which the jammer signals are suppressed.

Since the DOAs of the jammer signal are classified, we can find the corresponding amplitude to jammer’s DOA in (15). Using these poles and amplitudes, we can obtain the received signal vector without the jammer signals as where,

In (19), since is composed of only SOIs without jamming signals, it can be rewritten as

By applying the same number of matched filters as transmitters to each receiver, we obtain the receiving sequence matrix:

If we choose an orthonormal sequence for , the matrix becomes an identity matrix . Also, we assume the noise to be uncorrelated to .

In the case of a uniformly spaced linear array (ULA), (22) can be rewritten as the following: where and . Assuming , (23) can be simplified as where .

Applying the second step MPM to (24), we can obtain the DOAs of the SOI. The reason why the proposed method needs two applications of the MPM is because the first step MPM may fail to extract the DOAs of the SOI under very high JNR. Suppressing the jammer signals discriminated by the GLRT method, the second step of MPM can successfully extract the DOAs of the SOI from (24).

#### 4. Computational Complexity

In this section, we compare the computational complexity between the proposed method and the conventional GLRT method. Assuming that is the number of angle samples, the computational complexity of the conventional method is approximately given by

The computational complexity of the proposed method depends on the computation of the singular value decomposition (SVD) for the MPM. Assuming that and are the pencil parameters of the first and second MPM, respectively, the computational complexity of the proposed method is approximately given by [26]:

In this paper, we use and , and therefore the effect of is dominant. Figures 2 and 3 show a comparison of the computational complexity between the proposed method and the conventional GLRT method according to the number of virtual arrays in the MIMO radar. For the conventional GLRT method, we use . In Figure 2, when and , the computational complexity of the proposed method is smaller than that of the conventional method. Figure 3 shows that when and , the complexity of the proposed method is smaller.

Therefore, the proposed algorithm is effective when the number of the virtual arrays in the MIMO radar is limited according to the number of data samples .

#### 5. Numerical Examples

Figure 1 is simulation geometry for collocated MIMO radar. In this figure, means units of SOI’s DOA, and means the DOA of the th jammer signal. The array and the array have the same phase center and it is assumed that the array elements locate on the same straight line.