Robust mainlobe interference suppression for coherent interference environment

Abstract

A mainlobe interference suppression method is proposed in this paper, which can still work when the signal of interest (SOI) is present in the training data. In this method, the iterative adaptive approach (IAA) is applied to spatial spectrum estimation at first. Then, IAA spatial spectrum is used to reconstruct the interference-plus-noise covariance matrix (INCM). Next, the eigenvector associated with mainlobe interference in INCM is determined, and the eigen-projection matrix can be calculated to suppress the mainlobe interference. Meanwhile, the sidelobe-interference-plus-noise covariance matrix (SINCM) can be reconstructed. Finally, the adaptive weight vector is obtained. One main advantage is that the proposed method can deal with coherent mainlobe interference and sidelobe interferences simultaneously. Simulation results demonstrate the effectiveness and robustness of the proposed method.

Wang et al. EURASIP Journal on Advances in Signal Processing (2016) 2016:135 DOI 10.1186/s13634-016-0434-z

EURASIP Journal on Advances in Signal Processing

http://crossmark.crossref.org/dialog/?doi=10.1186/s13634-016-0434-z&domain=pdf

Web End = http://crossmark.crossref.org/dialog/?doi=10.1186/s13634-016-0434-z&domain=pdf

Web End = Robust mainlobe interference suppression for coherent interference environment

Yasen Wang*, Qinglong Bao and Zengping Chen

Abstract

A mainlobe interference suppression method is proposed in this paper, which can still work when the signal of interest (SOI) is present in the training data. In this method, the iterative adaptive approach (IAA) is applied to spatial spectrum estimation at first. Then, IAA spatial spectrum is used to reconstruct the interference-plus-noise covariance matrix (INCM). Next, the eigenvector associated with mainlobe interference in INCM is determined, and the eigen-projection matrix can be calculated to suppress the mainlobe interference. Meanwhile, the sidelobe-interference-plus-noise covariance matrix (SINCM) can be reconstructed. Finally, the adaptive weight vector is obtained. One main advantage is that the proposed method can deal with coherent mainlobe interference and sidelobe interferences simultaneously. Simulation results demonstrate the effectiveness and robustness of the proposed method.

Keywords: Mainlobe interference suppression, Iterative adaptive approach, Covariance matrix reconstruction, Eigen-projection

1 Introduction

Adaptive beamformers have wide applications in radar, navigation, sonar, radio astronomy, microphone array speech processing, and so on. Conventionally, adaptive beamformers are applied to suppress interferences that exist in sidelobe area [1, 2]. This kind of interference is called sidelobe interference. However, when the interference from mainlobe area exists, it will induce peak offset and distortion of the mainlobe and heightening of the sidelobe level. This kind of interference is mainlobe interference. Nowadays, as the electromagnetic environment becomes more complex, the mainlobe interference appears more frequently. Alternative solutions such as the blocking matrix preprocessing (BMP) method [3] and the eigen-projection matrix preprocessing (EMP) method [4] for suppressing the mainlobe interference are proposed. However, the former could work well provided that the direction of mainlobe interference and giving rise to the reduction of the effective degrees of freedom, and the latter relies on the absence of the signal of interest (SOI) in the training data. Furthermore, the problem of mainlobe peak offset cannot be overcome

by both methods. Later, a modified EMP method [5] is proposed for solving the mainlobe peak offset problem. Another new method based on EMP method and similarity constraints [6] can suppress multiple mainlobe interferences. Both of them make the assumption that there is no SOI appearing in the training data, which is not always available in many applications such as passive location, mobile communications, and radio astronomy[7]. For example, in the case of passive radar, the transmit signal may already exist when the receiver starts working. Besides, coherent interference is another problem in adaptive beamforming. In applications like passive radar [8] and meter-wave radar [9], it is likely that target echoes are masked by multipath echoes, which would lead to coherent interference. When this coherent interference appears in mainlobe area, the methods mentioned above are no longer effective.

Recently, some new beamforming techniques have been proposed [1012]. Among them, an interference-plus-noise covariance matrix (INCM) reconstruction processing is put forward in [8] to eliminate the SOI component when calculating the adaptive weight vector. However, the performance of this procedure degrades when coherent interference exists. In [13] and [14], an iterative adaptive approach (IAA) is proposed, which has

* Correspondence: mailto:[email protected]

Web End [email protected] College of Electronic Science and Engineering, National University of Defense Technology, Changsha 410073, China

The Author(s). 2016 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Wang et al. EURASIP Journal on Advances in Signal Processing (2016) 2016:135 Page 2 of 7

high accuracy of direction of arrival (DOA) and power estimation even when coherent sources and finite-sampling effects exist. In this paper, the IAA spatial spectrum is used to reconstruct the INCM. Then, inspired by the method proposed in [5], mainlobe and sidelobe interferences can be suppressed using the reconstructed INCM. Compared to the existing methods, the proposal in this paper can handle the coherent main-lobe interference even when the SOI is present in the training data.

2 Signal model

Consider an array composed of M omnidirectional antenna elements, the output of a narrowband adaptive beamformer at the kth snapshot y(k) is given by y(k) = wHx(k), where w CM 1 is the complex beamforming weight vector, x(k) stands for the signal received by the array antenna, and ()H denotes the Hermitian transpose. Assume there are Q sidelobe interferences and one mainlobe interference, where Q + 1 < M; then, x(k) is given by

x k

xs0 k

xint k

n k

X

Q1

q0

^

PCapon

1

4

where d() is the SV associated with the array structure and a hypothetical direction . Whenever the training data contains the SOI, SMI beamformer is in essence a minimum power distortionless response (MPDR) beam-former. However, the MPDR beamformer is more sensitive than MVDR beamformer when SV error exists [15].

3 Proposed method

In this paper, a mainlobe coherent interference suppression algorithm is proposed based on IAA spatial spectrum, INCM reconstruction, eigen-projection processing, and sidelobe-interference-plus-noise covariance matrix (SINCM) reconstruction.

3.1 IAA algorithm

IAA is a nonparametric adaptive algorithm. This algorithm assumes the presence of a large number of potential targets uniformly located on a fine grid with P M points in it. Let P be a P P diagonal matrix whose diagonal contains the power at each angular grid point. Besides, define A

a1 a2 aP

with ap denoting the SV

of the pth source and s k

s1 k

s2 k

sP k

T is

the source waveform vector. Define the interference and

noise covariance matrix Gp to be Gp R sp

2a

paHp where

dH

^

R1d

a q

s

q k n k

1

where sq(k) and a(q) are the qth signal waveform and the corresponding steering vector (SV), respectively. q stands for the qth signal DOA, and n(k) denotes the

additive white Gaussian noise. xs0(k) = a(0)s0(k) and xint

k

XQ1

q1

R = APAH; then, the cost function is given by

a q

s

XKk1y k sp k ap

2

Gp1 5

q k denote the SOI and interference vectors, respectively. sQ + 1(k) denotes the mainlobe interference waveform. The INCM can be written as

Rin E xint k

n k

xint k

n k

H

n o

where x

k k2Gp1 xHGp

2

The solution to the minimum variance distortionless response (MVDR) beamformer is given by

wopt

R1ina 0

a 0

HR1ina 0

1 x . The minimization of (5) leads to the IAA algorithm. The initialization of IAA is done by the data-independent delay-and-sum (DAS) beamformer. The steps of IAA can be summarized as follows.

Initialization:Set ^sp k aHpy k

=M, p = 1, , P, k = 1, , K

^

Pp 1=KXK

k1

^

CM M

sp k

2;

p 1; ; P

Iteration:

R A

3

Since Ri + n is unavailable when the SOI is present in the training data, it is usually replaced by the sample co

variance matrix ^

R 1K X

K

k1

^

PAH

for p = 1, , P

wp

R1ap

aH

p R1ap; ^

xH k

, where K is the number of snapshots. Then, the adaptive beamformer

wSMI

^

Pp wHp

x k

^

Rwp

is the sample matrix inversion (SMI)

beamformer. The corresponding Capon spatial spectrum distribution can be calculated as

R1a 0

a 0

H

^

R1a 0

end for Termination:

^

pi represents the vector containing the power estimates at the ith iteration in IAA algorithm. Thus, the iteration is

Wang et al. EURASIP Journal on Advances in Signal Processing (2016) 2016:135 Page 3 of 7

terminated when the relative change ^

p i1

p i ^

F

sidelobe interferences Q can be determined. Besides, when more than one mainlobe interference exists, the number of mainlobe interference can also be determined through the same method.

The correlation coefficient of two vectors x1 and x2 is

defined as cor x1; x2

xH1x2 x1

k k x2

k k . Then, the correlation coefficients of SOI SV a ^

0

and interference eigenvec

tors ei(i = 1, , Q + 1) are calculated. These correlation coefficients can be maximized when the interference is a mainlobe interference. Therefore, the mainlobe interference eigenvector em can be determined as

cor em; a ^

0

max

is

less than a specified tolerance, where F denotes the Frobenius norm.

An accurate spatial spectrum can be obtained from the results of IAA, and the SOI DOA ^

0 is also determined by searching for the peak of IAA spatial spectrum in the possible angle sector. Besides, the performance of IAA is not affected much by the presence of coherent sources [14]. For each IAA iteration, the complexity is O(M2P). Thus, a bigger P results in a higher accuracy as well as a larger calculation.

3.2 INCM reconstruction

According to [7], INCM can be reconstructed as

^

Rin Z

ei cor ei; a ^

0

11

d 6

where is the angular sector in which only the desired signal is located and denotes the complement sector of . When mainlobe interference exists, it also locates in . is easily determined by the results of IAA spatial spectrum. Then, the INCM can be reconstructed as

^

Rin X

p

^

Ppd p

Then, the corresponding eigenvalue m can be obtained. When multiple mainlobe interferences like J mainlobe interferences exist, em can be determined by searching for the biggest J values of the correlation coefficients.

The eigen-projection matrix can be calculated as

^

B Iem eHmem

1eH

m

^

PCapon d

dH

12

dH

p

7

where I is the identity matrix. ^

B is adopted to suppress the mainlobe interference in echo data X as follows

Y

^

BX 13

3.4 SINCM reconstruction

d represents the average value of the noise eigenvalues. In order to remove the effect of background noise perturbation, d is used to replace the noise eigenvalues. d can be obtained by

d ^

2n

where ^

Pp is obtained through IAA.

3.3 Eigen-projection processing

The reconstructed INCM is eigen-decomposed as

^

Rin X

M

i1

ieieHi EssEHs EnnEHn 8

where {i, i = 1, , M} is the eigenvalues of ^

Rin in the descending order and ei is the ith eigenvector. Es = [e1, e2, , eQ + 1] spans the interference subspace and En = [eQ + 2, , eM] spans the noise subspace. s = diag[1, , Q + 1] and n = diag[Q + 2, , M].

The estimation of Q is a public problem and is very critical. For an array with known array structure and SOI DOA, the mainlobe width is fixed. For example, for a uniform linear array (ULA), the mainlobe width BW0 can be expressed as

BW0 2 arcsin

Md sin0

9

Therefore, the mainlobe angle sector can be defined as

0BW0

2 ; 0 BW0

Q2 M

MQ1 14

m is the eigenvalue corresponding to em. According to [4] and [3], a noticeable part of mainlobe interference will be eliminated if its eigenvalue m is replaced by d.

Thus, the SINCM can be reconstructed as

^

Rre E

^

E 15 where ^

diag 1; ; m1; ^

2n; m1; ; Q1; ^

2n; ; ^

2n

is the reconstructed eigenvalue matrix and the eigenvector matrix E Es En

remains unchanged.

Based on the reconstructed SINCM ^

Rre , the weight

vector of the proposed algorithm is calculated as

wre re

^

2

10

R1rea ^

0

16

denotes the complement sector of . By using and the results of IAA spatial spectrum, the number of

Then, the output of the adaptive beamformer is calculated as

Wang et al. EURASIP Journal on Advances in Signal Processing (2016) 2016:135 Page 4 of 7

BX 17 Processing block diagram of the proposed algorithm is described in Fig. 1. The known parameters of this algorithm contain array geometry and a presumed SOI direction sector. Based on these parameters, the proposed method can suppress coherent mainlobe and sidelobe interferences when SOI appears in the training data. Moreover, the proposed method improves robustness of the beamformer from the two aspects. First, since the INCM reconstruction procedure converts a MPDR beamformer to a MVDR beamformer, the proposed method is less sensitive to SV error when SOI appears in the training data. Second, due to the high accuracy of IAA spatial spectrum, the presumed SOI direction error can be corrected.

4 Simulation results

In simulations, we consider a uniform linear array with M = 16 omnidirectional antenna elements spaced half a wavelength. The SOI direction is 0, and the mainlobe interference direction is 5 with interference-to-noise ratio (INR) 5 dB. Two sidelobe interferences impinging

from 45 to 35 are considered with INR 25 and 35 dB, respectively. Noise is assumed to be the Gaussian white additive noise. The number of snapshots is fixed to be K = 100. The number of sidelobe interferences Q is assumed to be known in simulations. For each scenario, 200 Monte Carlo runs are performed.

The proposed method is compared to SMI, BMP, EMP, and the method proposed in [5] (EMP-CMR). The SOI is not present in the training data when using EMP and EMP-CMR methods. In BMP method, the exact mainlobe interference angle is known. For Capon and IAA, the scanning grid is uniform in the range from 90 to 90 with 1 increment between adjacent grid points. The presumed SOI angular sector is set as = [2, 2] in the proposed method.

First, the input signal signal-to-noise ratio (SNR) is fixed to be 0 dB. We consider the accuracy of IAA, MUSIC, Capon, and forward-backward spatial smoothing (MSS) Capon spatial spectrum [16] when the SOI and the interference are coherent. When using the forward-backward spatial smoothing method, the array is divided into 2 sub-arrays and each sub-array contains 15 sensors. In Fig. 2a, the SOI and all the interferences are uncorrelated. In Fig. 2b, only the SOI and the main-lobe interference are coherent. In Fig. 2c, the SOI, the mainlobe interference, and the sidelobe interference from 45 are coherent. Figure 2 shows that IAA spatial spectrum outperforms the MUSIC and Capon spatial spectrum when coherent interference exists. Besides, the noise floor of the MSS Capon spatial spectrum is significantly higher than the IAA spatial spectrum.

Next, we compare the output signal-to-interference-plus-noise (SINR) of different methods when input signal SNR varies. The presumed SOI direction is 0. Figure 3 shows that EMP-CMR method outperforms other methods when all the sources are uncorrelated. Besides, the output SINR of the proposed method is better than BMP, EMP, and SMI methods. Figure 4 presents the output SINR when the SOI and the mainlobe interference are coherent. For EMP and EMP-CMR methods, the SOI is absent in the training data. Thus, the adaptive weighted vectors of these two methods are not affected by the coherent signals in this scenario. From Figs. 3 and 4, we ob-serve that the performance of EMP and EMP-CMR remains unchanged. However, for the proposed method, its output SINR improves a lot due to better accuracy of IAA spatial spectrum. In Fig. 5, the adaptive array patterns of the quiescent (QUI), EMP-CMR, and the proposed method are depicted when SNR is 0 dB. We note that the mainlobe peak offset is eliminated in both EMP-CMR method and the proposed method. Figure 6 shows the output SINR when the SOI, the mainlobe interference, and the sidelobe interference impinging from 35 are all coherent. For EMP and EMP-CMR methods, the

Z wHreY wHre

^

Fig. 1 Processing block diagram of the proposed algorithm

Wang et al. EURASIP Journal on Advances in Signal Processing (2016) 2016:135 Page 5 of 7

herent sidelobe interferences. Therefore, the performance of these two methods degrades severely in Fig. 6. For the proposed method, its performance also degrades.

Fig. 4 Output SINR against input SNR for the case of two coherent sources

Wang et al. EURASIP Journal on Advances in Signal Processing (2016) 2016:135 Page 6 of 7

Fig. 5 Adaptive array patterns of QUI, EMP-CMR, and proposed method

where q denotes the norm of eq and is uniformly distributed in the interval [0,0.4] in each simulation run and qm m 1; 2; ; M

is phases of mth coordinate of

the random error vector eq, which are independently and uniformly taken from the interval [0, 2] in each simulation run. ^

q is uniformly distributed in the interval [2, 2] in each simulation run. This kind of SV mismatch model is universal because it can be caused by calibration error, local scattering error, look direction error, and so on. The other simulation conditions are the same to the simulation sets of Fig. 4. The output SINR versus the input SNR is shown in Fig. 7. From Fig. 7, we notice that the SV error has more side effects

on EMP-CMR method and the proposed method. However, the performance of EMP-CMR method and the proposed method is still better than other methods.

5 Conclusions

In this paper, we have developed a new mainlobe interference suppression method based on IAA spatial spectrum, INCM reconstruction, eigen-projection processing, and SINCM reconstruction, which annihilates the prior information restrictions existing in the previous methods for mainlobe interference suppression. Vigorous steps have been taken to reduce these restrictions, and simulation results demonstrated that the proposed method exhibits superior performance when coherent

Wang et al. EURASIP Journal on Advances in Signal Processing (2016) 2016:135 Page 7 of 7

mainlobe interference and sidelobe interferences exist. In addition, robustness of the proposed method against SV error is also verified in this paper.

Acknowledgements

This work was supported in part by the National Natural Science Foundation of China under Grant 61401489.

Competing interests

The authors declare that they have no competing interests.

Received: 5 March 2016 Accepted: 1 December 2016

References1. W Liu, J Liu, L Wang, K Duan, Z Chen, Y Wang, Adaptive array detection in noise and completely unknown jamming. Digital Signal Processing. 46, 4148 (2015)

2. W Liu, J Liu, X Hu, Z Tang, L Huang, Y Wang, Statistical performance analysis of the adaptive orthogonal rejection detector. IEEE Signal Processing Letters. 23(6), 873877 (2016)

3. RF Li, YL Wang, SH Wan, Robust adaptive beam forming under main lobe interference conditions. System Engineering and Electronics. 24(7), 47 (2002)

4. RF Li, YL Wang, SH Wan, Research of reshaping adapted pattern under mainlobe interference conditions. Modern Radar. 24, 5055 (2002)

5. XP Yang, Z Zhang, T Zeng, T Long, TK Sarkar, Mainlobe interference suppression based on eigen-projection processing and covariance matrix reconstruction. IEEE Antennas Wireless Propag. Lett. 13, 13691372 (2014)

6. J Qian, Z He, Mainlobe interference suppression with eigenprojection algorithm and similarity constraints. Electronics Lett. 52, 228230 (2016)

7. SA Vorobyov, Principles of minimum variance robust adaptive beamforming design. Signal Process. 93, 32643277 (2013)

8. Special issue on passive radar systemsIEE Proceedings on Radar, Sonar and Navigation. 152(3), 106-223 (2005)

9. Y Zhao, R Li, Y Wang, W Liu, Angles measurement of meter-wave radars by mainlobe multipath jamming suppression. J. Huazhong Univ. of Science and Technology (Natural Science Edition) 41(4), 5155 (2013)

10. Y Gu, A Leshem, Robust adaptive beamforming based on interference covariance matrix reconstruction and steering vector estimation. IEEE Signal Process. Mag. 60(7), 29742987 (2012)

11. J Liu, H Li, B Himed, Joint optimization of transmit and receive beamforming in active arrays. IEEE Signal Process Lett. 21(1), 3942 (2014)

12. J Liu, W Liu, H Liu, B Chen, XG Xia, F Dai, Average SINR calculation of a persymmetric sample matrix inversion beamformer. IEEE Signal Process. Mag. 64(8), 21352145 (2016)

13. T. Yardibi, J. Li, P. Stoica, M. Xue, A.B. Baggeroer. Source localization and sensing: a nonparametric iterative adaptive approach based on weighted least squares. IEEE Trans. Aerospace and Electronic Systems. (2008)

14. L Du, T Yardibi, J Li, P Stoica, Review of user parameter-free robust adaptive beamforming algorithms. Digital Signal Process. 19(4), 567582 (2009)15. L. Ehrenberg, S. Gannot A. Leshem, E. Zehavai. Sensitivity analysis of MVDR and MPDR beamformers. IEEE 26th Convention of Electr. Electron. Eng. in Israel, Eliant, Israel, Nov., 2010.

16. W Liu, L Zhang, Robust beamforming for coherent signals based on the spatial-smoothing technique. Signal Process. 92(11), 27472758 (2012)

17. L Huang, XX Zhang, Z Ye, Robust adaptive beamforming with a novel interference-plus-noise covariance matrix reconstruction method. IEEE Trans. Signal Process. 63, 16431650 (2015)