Next Article in Journal
Cryptographically Secure Pseudo-Random Number Generator IP-Core Based on SHA2 Algorithm
Next Article in Special Issue
Efficient Two-Dimensional Direction Finding Algorithm for Rectilinear Sources Under Unknown Mutual Coupling
Previous Article in Journal
A Reliable Prognosis Approach for Degradation Evaluation of Rolling Bearing Using MCLSTM
Previous Article in Special Issue
An Array Switching Strategy for Direction of Arrival Estimation with Coprime Linear Array in the Presence of Mutual Coupling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Null Broadening Beamforming Based on Covariance Matrix Reconstruction via Virtual Interference Sources

1
School of Electronic Engineering, Xidian University, Xi’an 710071, China
2
School of Engineering, Rocket Force University of Engineering, Xi’an 710025, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(7), 1865; https://doi.org/10.3390/s20071865
Submission received: 9 February 2020 / Revised: 23 March 2020 / Accepted: 25 March 2020 / Published: 27 March 2020

Abstract

:
When jammers move rapidly or an antenna platform travels at high speed, interference signals may move out of the null width in the array beampattern. Consequently, the interference suppression performance can be significantly degraded. To solve this problem, both the null broadening technique and robust adaptive beamforming are considered in this paper. A novel null broadening beamforming method based on reconstruction of the interference-plus-noise covariance (INC) matrix is proposed, in order to broaden the null width and offset the motion of the interfering signals. In the moving case, a single interference signal can have multiple directions of arrival, which is equivalent to the existence of multiple interference sources. In the reconstruction of the INC matrix, several virtual interference sources are set up around each of the actual jammers, such that the nulls can be broadened. Based on the reconstructed INC and signal-plus-noise covariance (SNC) matrices, the steering vector of the desired signal can be obtained by solving a new convex optimization problem. Simulation results show that the proposed beamformer can effectively broaden the null width and deepen the null depth, and its performance in interference cancellation is robust against fast-moving jammers or array platform motion. Furthermore, the null depth can be controlled by adjusting the power parameters in the reconstruction process and, if the direction of interference motion is known, the virtual interference sources can be set to achieve better performance.

1. Introduction

As antenna arrays have been applied to radar, sonar, radio astronomy, medical imaging, and in other areas, beamforming is widely used for interference suppression, target detection, and direction of arrival (DOA) estimation [1,2]. The purpose of adaptive beamforming is to enhance the desired signal while simultaneously suppressing interference in the spatial domain [3]. Under mismatches, such as wavefront distortion, array calibration errors, look direction errors, and incoherent local scattering, traditional beamformers suffer from severe performance degradation [4]. To address this problem, various robust adaptive beamformers have been proposed. In recent years, due to its good robustness against mismatches, a novel robust adaptive beamforming method based on covariance matrix reconstruction and steering vector estimation has attracted much attention [5,6,7,8,9]. However, in the presence of fast-moving interference sources or antenna platform movement, such as that which occurs in highly dynamic conditions, the directions of arrival (DOAs) of the interference signals may range over multiple resolution units and the beamforming weights cannot be updated in real time, which will cause the interference to move out of the nulls and decrease the anti-interference capacity of adaptive beamforming techniques. An alternative method is to keep the interference in the nulls by broadening the null widths; that is, a wider attenuation is formed in the region where the interference may appear, in order to effectively suppress the interference direction perturbation. Some research has been carried out in the area of null broadening techniques, including virtual interference sources, derivative constraints, and projection transformations.
Therefore, null broadening techniques can provide effective and robust approaches to overcoming this problem. In References [10,11,12,13,14], these approaches were based on the covariance matrix tapering and effectively broadened the null width. However, this technique decreases the array gain and the depth may become shallower [15]. Semidefinite programming was introduced to improve the performance of adaptive null broadening beamforming in References [16,17]. Moreover, Qian et al. [18] first imposed nulls toward the regions of nonstationary interference via the reconstruction of the interference-plus-noise covariance matrix, and a similarity constraint was enforced at the design stage to restrict the shape of the beam pattern. Then, the adaptive weight vector could be computed via maximizing a new signal-to-interference-plus-noise ratio (SINR) criterion subject to the similarity constraint. In mathematics, the optimisation problem was a nonconvex fractional quadratically constrained quadratic programming (QCQP) problem, which was converted into a convex optimisation problem using semidefinite programming techniques. In Reference [19], by combining with white-noise constraint adaptive processing, the null widening technique could provide effective nulling under snapshot-deficient conditions. Nevertheless, compared to the previous methods, these approaches only result in a small improvement in performance. A new null broadening beamformer based on projection and diagonal loading has been proposed in Reference [20]. Although this approach successfully broadens the null width, its shortcoming is shallower nulls. In Reference [21], by means of the covariance matrix of the auxiliary elements, a null broadening method was realized based on the sidelobe canceller. This approach demonstrated good performance in practical applications, but has increased hardware complexity. In References [22,23], an adaptive null broadening technique based on reconstruction of the covariance matrix was proposed. However, these approaches have poor robustness against various mismatches in the antenna array. Ors et al. [24] proposed first- and second-order iterative null broadening beamforming, and illustrated that the second order required fewer iterations compared to the first order. Nevertheless, the robustness of the beamforming was not considered either. Xiao et al. [25] proposed a method based on optimization algorithm and neural networks to the complicated wide null problem, and the performance of the data-based general regression neural network wide nulling model was demonstrated by numerical experiment. In Reference [26], Z. Liu et al. proposed a computationally efficient null widening method for sidelobe canceller, which is a covariance matrix taper-based method and puts fictitious interferers into snapshots to broaden the sharp null. Based on Mailloux’s methodology, the covariance matrix and cross-correlation vector were tapered via random disturbance. Compared with the existing methods, the method required much less computation, but its performance is similar. However, the performance of the proposed method in the presence of array errors needs to be further studied. In order to solve the problem that the performance of existing null broadening beamforming will decrease when array geometry errors are present, Zhao et al. [27] proposed a null broadening beamforming only against array geometry errors, which was not applied against other errors, such as unknown array steering vector errors, and incoherent local scattering errors. Yang et al. [28] proposed robust wideband adaptive beamforming with null broadening and constant beamwidth via adding virtual interferers around original interferers. For that algorithm, the amount of virtual interference was an important parameter for the taper matrix, but it was difficult to determine and the method of adding virtual interferers was relatively complex. The interference-plus-noise covariance matrix was reconstructed with the simplified power spectral density function and broadened notches corresponding to the time-varying DOAs of the moving interference signals [29]. This method could impose nulls towards the regions of the moving interference signals and obtain lower sidelobes in the beampattern. However, except for DOA mismatch, its robustness against other array errors needs to be further verified and improved. In addition, Wang et al. [30] and Cong et al. [31] proposed Null broadening beamforming methods for anti-jamming in GNSS receivers. However, for their methods, the interferers were assumed to be in a certain distribution, for example, Gauss or trigonometric distribution, which was often not satisfied.
In this paper, we combine the null broadening technique with robust adaptive beamforming and propose a novel adaptive null broadening beamformer based on covariance matrix reconstruction. First, a new interference-plus-noise covariance (INC) matrix is reconstructed by placing several virtual interference sources around each of the actual jammers. Then, the signal-plus-noise covariance (SNC) matrix is reconstructed using the Capon spatial spectrum [32]. Next, the steering vector of the desired signal can be estimated by solving a new convex optimization problem. Finally, based on the reconstructed INC matrix and the estimated steering vector of the desired signal, the weight vector can be obtained according to the minimum variance distortionless response (MVDR) beamformer. Simulation results demonstrate that the robust adaptive beamformer can effectively broaden the null width and has good performance.
The main contributions of this paper can be summarized as follows:
1) Compared with other existing null broadening beamforming techniques, we propose an algorithm based on INC matrix reconstruction by setting up several virtual interference sources around the actual jammers directly in the Capon spatial spectrum, which can simultaneously broaden the nulls and eliminate the self-null phenomenon of the desired signal. Furthermore, the null depth and width can be controlled by setting the parameters of the virtual interference sources.
2) Based on the reconstructed INC and SNC matrices, we design a convex optimization problem to obtain the steering vector of the desired signal, in order to further improve the robustness of algorithm.
The rest of this paper is organized as follows. In Section 2, the signal model, including the antenna array and the received signal, is given. In Section 3, the proposed algorithm is introduced, which consists of INC matrix reconstruction, SNC matrix reconstruction, and steering vector estimation. Section 4 presents the simulation setup and results. Finally, Section 5 gives some concluding remarks.

2. Problem Background

Assume that a uniform linear array (ULA) comprised of N omni-directional sensors receives M + 1 far-field narrowband source signals. The first signal is the desired signal, impinging upon the array from the direction θ s , while the remaining M interfering signals come from directions θ 1 , θ 2 , , θ M , respectively. The array sample data at the k th snapshot can be modeled as
x ( k ) = x s ( k ) + x int ( k ) + n ( k ) = x s ( k ) a ( θ s ) + m = 1 M x int ( k ) a ( θ m ) + n ( k ) ,
where a ( θ s ) is the steering vector of the desired signal, a ( θ m ) is the steering vector of the m th interference signal, and x s ( k ) = x s ( k ) a ( θ s ) , x int ( k ) = m = 1 M x int ( k ) a ( θ m ) , and n ( k ) denote the desired signal, interference signals, and additive white Gaussian noise vectors, respectively. In addition, we assume that they are statistically independent from each other. The steering vector of the ULA can be written as
a ( θ ) = [ 1 , e j 2 π sin θ d λ , , e j 2 π sin θ ( N 1 ) d λ ] T ,
where ( · ) T is the transpose product, d is the sensor interval, λ is the wavelength, and θ is the DOA of the impinging signal. The theoretical covariance matrix of the array data can be obtained by
R x = E { x ( k ) x H ( k ) } = m = 0 M δ m 2 a ( θ m ) a H ( θ m ) + δ n 2 I ,
where ( · ) H is the Hermitian transpose, E ( · ) is the mathematical expectation, { σ m 2 } m = 0 M are the powers of the impinging signals, σ n 2 is the noise power, and I is the N × N identity matrix. In practice, R x is usually unavailable, and is often replaced by the sample covariance matrix
R ˜ x = 1 K k = 1 K x ( k ) x H ( k ) ,
where K is the number of snapshots. The output of the adaptive beamformer is expressed as
y ( k ) = w H x ( k ) ,
where w is the N × 1 weight vector. With the MVDR criterion, the optimal beamformer can be obtained by solving the following problem
min w w H R int + n w subject to w H a ( θ s ) = 1 ,
where R int + n denotes the INC matrix. Therefore, the optimal weight vector can be expressed as
w o p t = R int + n 1 a ( θ s ) a H ( θ s ) R int + n 1 a ( θ s ) .
However, in practice, the INC matrix cannot be obtained directly, and the sample covariance matrix R ˜ x is usually used to replace the INC matrix, which may cause the self-nulling phenomenon of the desired signal and result in significant performance deterioration, especially with high SNR [4].

3. The Proposed Algorithm

In this section, a new adaptive null broadening beamforming method is proposed. The main idea is to reconstruct a new INC matrix by utilizing virtual interference sources and estimate the steering vector of the desired signal by solving a new optimization problem. As the sample covariance matrix is usually inevitably mixed into the desired signal, which can cause the self-null phenomenon with a high signal-to-noise ratio (SNR) and result in significant performance deterioration, it is necessary to remove the desired signal component from the sample covariance matrix and reconstruct the INC matrix. On the other hand, to avoid convergence to the DOA of the interfering signal, the SNC matrix is also reconstructed to effectively improve the pointing accuracy of the estimated steering vector. Moreover, the method used to reconstruct the INC matrix is also used for SNC matrix reconstruction. Therefore, this section is divided into three parts: The first part discusses the INC matrix reconstruction, the second reconstructs the SNC matrix, and the last estimates the steering vector of the desired signal based on the INC and SNC matrices. Furthermore, the complete algorithmic flow is detailed in the third part.

3.1. Inc Matrix Reconstruction

The Capon spatial spectrum [32] estimator is
P ( θ ) = 1 a H ( θ ) R ˜ x 1 a ( θ ) .
Using the Capon spatial spectrum (8), an INC matrix reconstruction approach has been proposed [4,5], which uses the following
R ˜ i + n = Θ ¯ P ( θ ) a ( θ ) a H ( θ ) d θ = Θ ¯ a ( θ ) a H ( θ ) a H ( θ ) R ˜ x 1 a ( θ ) d θ ,
where Θ ¯ denotes the angular region excluding the desired signal region Θ , Θ Θ ¯ covers the whole spatial domain, and Θ Θ ¯ is empty. Considering the structured feature of the integral, the reconstructed INC matrix can be approximately calculated by a discrete sum in the discrete time domain as
R ˜ i + n k = 1 Θ ¯ Θ ¯ Δ θ Δ θ P ( θ k ) a ( θ k ) a H ( θ k ) Δ θ ,
where Δ θ is the angle interval, which should be small.
In the INC matrix reconstruction, the range of virtual interference sources can be determined according to the motion of the antenna platform or jammers. Assume that M regions of virtual interference sources are located in [ ϑ 1 β 11 , ϑ 1 + β 12 ] , …, [ ϑ M β M 1 , ϑ M + β M 2 ] , respectively, and the spatial spectra are correspondingly set as κ 1 P ( ϑ 1 ) , , κ M P ( ϑ M ) . In the array beampattern, the numerical pairs ( β 11 , β 12 ) , …, ( β M 1 , β M 2 ) determine the regions of null broadening (which may be different in practice) and the respective null widths are β 12 + β 11 , …, β M 2 + β M 1 . In addition, β i j ( 1 i M , 1 j 2 ) is usually greater than zero. It should be noted that the relative motion between the interference sources and the antenna array is known, which usually could be obtained in advance by DOA estimated techniques and the inertial measurement system, the virtual interference sources can be set up according to the incoming direction of the interference signals to achieve better performance in interference suppression. In particular, it can be realized by controlling these parameters β i j . The parameters κ 1 , , κ M determine the depth of the nulls, which is closely related to the performance of interference suppression. Meanwhile, κ m will cause the mainlobe to widen and the sidelobe to be higher in the beampattern, which are considered to be negative effects. The parameters κ m are introduced to better suppress interference, where κ m can be determined by the value of P ( ϑ m ) . Specifically, the vector composed of P ( ϑ m ) , ( 1 m M ) is arranged in ascending order; that is, its first and last elements, respectively, correspond to κ 1 = 1 and κ M = 10 , and the remaining values of κ m can be obtained by linear fitting. Therefore, by (10), the INC matrix based on virtual interference sources in the Capon spatial spectrum can be reconstructed as
R ^ i + n = i = 1 Θ ¯ 1 Θ ¯ 1 Δ θ Δ θ κ i p ( θ i ) a ( θ i ) a H ( θ i ) Δ θ + k = 1 Θ ¯ 2 Θ ¯ 2 Δ θ Δ θ p ( θ k ) a ( θ k ) a H ( θ k ) Δ θ ,
where Θ ¯ 1 = ϑ 1 β 11 , ϑ 1 + β 12 ϑ 2 β 21 , ϑ 2 + β 22 ϑ M β M 1 , ϑ M + β M 2 denotes the angular sector including all interference and virtual interference signals, and Θ ¯ 2 denotes the angular sector excluding the desired, interference and virtual interference signals. Θ ¯ = Θ ¯ 1 Θ ¯ 2 , Θ ¯ 1 Θ ¯ 2 is empty.

3.2. Snc Matrix Reconstruction

In order to accurately estimate the steering vector of the desired signal, the SNC matrix can be reconstructed by the Capon spectrum. In the same manner as (9), we can reconstruct the SNC matrix as
R ˜ s + n = Θ ¯ a ( θ ) a H ( θ ) a H ( θ ) R ˜ x 1 a ( θ ) d θ ,
where Θ ¯ denotes the angular region excluding the interference signal regions Θ , Θ Θ ¯ covers the whole spatial domain, and Θ Θ ¯ is empty.

3.3. Desired Signal Steering Vector Estimation

Due to many mismatches, such as gain and phase perturbations, incoherent local scattering, signal look direction error, array geometry errors, the steering vector of the desired signal is often inaccurate. To obtain the optimal weight vector, we introduce the convex theory and obtain the accurate steering vector by solving the optimization problem. In the presence of steering vector mismatches, the mismatch vector a a ( θ s ) with a denoting the actual desired signal steering vector can be further decomposed into two components, that is, a a ( θ s ) = e + e | | . One denoted by e is orthogonal to a ( θ s ) , and the other denoted by e | | is parallel to a ( θ s ) , and e | | does not affect the beamforming quality and does not impact the output SINR [5]. Therefore, a new constraint can be established as
( a ^ a ( θ s ) ) H a ( θ s ) = 0 ,
where a ^ = a ( θ s ) + e . Furthermore, the spatial power estimate must satisfy
1 a ^ H R ˜ s + n 1 a ^ 1 a H ( θ s ) R ˜ s + n 1 a ( θ s ) .
Therefore, the optimization problem can be written as
min a ^ a ^ H R ˜ x 1 a ^ , subject to ( a ^ a ( θ s ) ) H a ( θ s ) = 0 , a ^ H R ˜ s + n 1 a ^ a H ( θ s ) R ˜ s + n 1 a ( θ s ) .
By solving the optimization problem, the desired signal steering vector can be obtained. Set a ^ = x + j y , a ( θ s ) = c + j b , and a H ( θ s ) R ˜ s + n 1 a ( θ s ) = Φ 1 , then the objective function can be further written as
( x T j y T ) R ˜ x 1 ( x + j y ) = x T R ˜ x 1 x + y T R ˜ x 1 y .
The first constraint can be expressed as
c T x + b T y c T c b T b = 0 b T x c T y = 0 ,
and the second constraint is equivalent to the following inequality,
x T R ˜ s + n 1 x + y T R ˜ s + n 1 y Φ 1 0 .
According to (16)–(18), the optimization problem can be further written as
min x , y x T R ^ x 1 x + y T R ^ x 1 y subject to c T x + b T y c T c b T b = 0 b T x c T y = 0 x T R ˜ s + n 1 x + y T R ˜ s + n 1 y Φ 1 0 .
With Karush–Kuhn–Tucker (KKT) conditions [33], the Lagrange function can be expressed as
L ( x , y , λ 1 , λ 2 , λ 3 ) = x T R ˜ x 1 x + y T R ˜ x 1 y + λ 1 ( c T x + b T y c T c b T b ) + λ 2 ( b T x c T y ) + λ 3 ( x T R ˜ s + n 1 x + y T R ˜ s + n 1 y Φ 1 ) .
Hence, we can obtain
x L ( x , y , λ 1 , λ 2 , λ 3 ) = 2 R ˜ x 1 x + λ 1 c + λ 2 b + 2 λ 3 R ˜ s + n 1 x = 0 y L ( x , y , λ 1 , λ 2 , λ 3 ) = 2 R ˜ x 1 y + λ 1 b λ 2 c + 2 λ 3 R ˜ s + n 1 y = 0 λ 1 L ( x , y , λ 1 , λ 2 , λ 3 ) = c T x + b T y c T c b T b = 0 λ 2 L ( x , y , λ 1 , λ 2 , λ 3 ) = b T x c T y = 0 λ 3 L ( x , y , λ 1 , λ 2 , λ 3 ) = x T R ˜ s + n 1 x + y T R ˜ s + n 1 y Φ 1 = 0 .
By solving (21), we can obtain x and y , following which the steering vector of the desired signal can be obtained as a ^ = x + j y . On the other hand, the optimization problem (15) can be easily solved using convex optimization software, such as CVX [34]. To date, both the INC matrix and the steering vector of the desired signal have been obtained. From (7), (11) and (15), the weight vector of the adaptive null broadening beamformer can be written as
As shown in Figure 1, the proposed approach can be illustrated step-by-step. Therefore, the proposed adaptive null broadening beamforming technique based on covariance matrix reconstruction via virtual interference sources can be summarized as follows:
Step 1. Calculate the Capon spatial spectrum (8) and search for the peaks ( ϑ 1 , , ϑ M ) from the interfering signal region;
Step 2. Reconstruct the new INC matrix (11) based on the virtual interference sources, where β i j controls the null width and κ i controls the null depth;
Step 3. Estimate the steering vector a ^ by solving the optimization problem (15); and
Step 4. Calculate the weight vector w o p t n u l l (22).
In the proposed algorithm, the computational complexity is dominated by solving the optimization problem, which can be efficiently solved using an interior point method with a computational complexity cost of O ( N 3.5 ) . The robust adaptive beamforming algorithm in Reference [5] has a complexity of O ( N 3.5 ) , due to the use of a quadratically constrained quadratic programming problem. Therefore, from a computational complexity point of view, the cost of the proposed algorithm is similar to that of the robust adaptive beamforming algorithm in Reference [5], and is comparable with that of other robust beamforming algorithms [5,9].
w o p t n u l l = R ^ i + n 1 a ^ a ^ H R ^ i + n 1 a ^ .

4. Simulation Results

In our simulations, a ULA consisting of N = 10 omni-directional sensors with an inter-sensor spacing of half a wavelength was considered. We assumed that the additive noise obeys a Gaussian distribution with unit covariance. In the simulations, the proposed algorithm was compared with the following methods: (i) the robust adaptive beamformer (Rec-ISVPE) of Reference [5], (ii) the Laplace-based null broadening algorithm of Reference [10], (iii) the Mailloux-based null broadening algorithm of Reference [11], (iv) the improved null broadening beamforming based on covariance matrix reconstruction (Rec-MVDR) of Reference [23], and (v) the analytical approach to null broadening (AANB) of Reference [15]. Among these methods, methods (i) and (iv) use INC matrix reconstruction, while the others are null broadening beamforming algorithms. For all tested beamformers, to effectively verify the performance of null broadening beamformers, the width of nulls is set to 4 , and the control parameter b of Rec-MVDR is set to 1.2. The optimal weight vector w o p t is obtained by (7), which is calculated from the exact INCM and the actual desired signal steering vector. Then, the output optimal SINR can be obtained by SIN R opt = σ s 2 w o p t H a ( θ s ) σ s 2 w o p t H a ( θ s ) w o p t H R int + n w o p t w o p t H R int + n w o p t , where σ s 2 denotes the desired signal power [5,35]. For each scenario (except for that in Section 4.1), 200 Monte Carlo runs were performed.

4.1. Performance of Null Broadening

In the first example, the beampatterns of null broadening were examined. Assumed that the desired signal with fixed 5 dB SNR comes from the direction θ s = 20 , and two interfering signals impinge upon the array antenna from the directions θ 1 = 30 and θ 2 = + 80 , respectively, with 40 dB interference-to-noise ratios (INRs). For all tested null broadening algorithms, we set up the regions of null broadening to be [ 32 , 28 ] and [ 38 , 42 ] . In the reconstruction of the INC matrix, the region of the desired signal was set to be located in Θ = [ 2 , + 8 ] , and the power parameters were set as κ 1 , 2 = 1 . In the simulations, the interval of the virtual interference sources was set as Δ η = 0.5 , which meant that there were eight virtual interferers in the vicinity of each interfering source.
Figure 2a demonstrates, from the array beampatterns, that null broadening was achieved. Compared to the Rec-ISVPE beamformer, which is a robust and adaptive beamforming technique without wide nulling, the proposed null broadening method could widen the main lobe and raise the sidelobes. As shown in Figure 2b,c, the beampattern of the proposed method was very close to that of Rec-MVDR, and their null depths significantly outperformed the other tested algorithms. In order to better compare the effect of null broadened, the numeric widths of nulls obtained by different algorithms in Figure 2 were measured, with the results reported in Table 1. These numeric widths were the averages of the two nulls at 30 and 80 . The results of Table 1 demonstrated that the proposed and Rec-MVDR methods could effectively broaden the nulls low to the beampattern gain at −150 dB, with the null widths similar to the target width 4 , while other tested beamformers could not. Moreover, the null at angle 80 was close to the end of the array beampattern, which caused a small degree of distortion of the beampattern, and the lowest beampattern gain of the nulls was below −150 dB, where the actual width was a little smaller than 4.2 and more similar to 4 . Furthermore, the beampattern of Rec-MVDR also confirmed these analyses. Therefore, the proposed method demonstrated good performance of null broadening.

4.2. Interference Suppression of Fast-Moving Jammers

In the second example, we examined the performance of the considered methods in suppressing fast-moving interfering signals. The directions of the two interfering signals were assumed to be random and uniformly distributed in [ 32 , 28 ] and [ 38 , 42 ] , respectively, while the desired signal impinged upon the array from the angular sector [ 8 , 12 ] . Moreover, the DOAs of the interferers varied in each trial, remaining fixed from snapshot to snapshot. In the performance comparison of output SINR versus input SNR, the input INRs and the number of snapshots were fixed as 40 dB and 60, respectively. In the performance comparison of output SINR versus the number of snapshots, the input SNR and INR were set as 10 dB and 40 dB, respectively.
Figure 3a shows that the performance of the proposed method was very close to the optimal SINR, while the input SNRs varied in a large range (i.e., from −30 dB to 50 dB). As shown from Figure 3b, the number of snapshots had no effect on the output SINR, and the performance of the proposed method was obviously superior to the others. Therefore, when interference motion leads to mismatches between the DOAs and the adaptive nulls, the proposed null broadening beamformer can effectively suppress interference.

4.3. Mismatch Due to Signal Look-Direction Error

In the third example, the impact of the signal look-direction error on the output SINR was examined. We assumed that the random DOA mismatches of the desired signal were uniformly distributed in [ 4 , 4 ] , and the motion of the strong interferers was the same as in the previous simulation. Note that the random DOAs of the desired and interference signals changed from trial to trial, but remained fixed from snapshot to snapshot. When examining the effect of the input SINR, the input INRs and the number of snapshots were fixed as 40 dB and 80, respectively. The input SNR was set to 10 dB in the performance comparison of output SINR versus the number of snapshots.
Figure 4a compares the output SINR versus the input SNR. It can been seen that the performance of the proposed algorithm was better than that of the other beamformers and was close to the optimal SINR. In Figure 4b, the output SINR versus the number of snapshots is illustrated. The number of snapshots had no obvious effect on the performance of all tested beamformers. Thus, it has been demonstrated that the performance of the proposed beamformer against DOA errors is obviously better than that of Rec-MVDR. Compared to the remaining beamformers, the proposed algorithm displayed significantly better performance.

4.4. Mismatch Due to Array Geometry Errors

In the fourth example, the effects of array geometry mismatches on the tested beamformers was investigated. Mismatches in antenna geometry are usually modelled by sensor position errors, which can be assumed to be subject to a Gaussian distribution. In the simulations, the mismatches were modelled as a Gaussian process with a mean of 0 and a variance of 0.002. At the same time, the DOAs of the interferers changed from run to run, but were kept fixed from snapshot to snapshot. In addition, the other parameters were the same as those in the third example.
From Figure 5a, it can be seen that the performance of the Laplace- and Mailloux-based methods was superior to that of the tested null broadening algorithms at low SNRs, but dramatically degraded with an increase of input SNR. For the proposed beamformer, array geometry mismatches hardly degraded the output SINR. From Figure 5a, at high SNRs, the proposed beamformer still outperformed the other tested methods. Figure 5b illustrates that the number of snapshots had no or very little effect on the performance of the tested beamformers. Therefore, due to array geometry errors, the performance of the proposed beamformer declined, to a certain degree, but its overall performance was still acceptable.

4.5. Mismatch Due to Incoherent Local Scattering

In the fifth example, the impact of the steering vector mismatch due to incoherent scattering was considered. Assuming that the steering vector of the desired signal changes over time, the desired signal can be modeled as [3]
x ˜ s ( k ) = x s ( k ) a ( θ s ) + p = 1 4 s p ( k ) a ( θ p ) ,
where a denotes the steering vector of the desired signal corresponding to the direct angle θ s , a ( θ p ) , p = 1 , 2 , 3 , 4 denote the steering vectors of the incoherently scattered signals, the directions θ p , p = 1 , 2 , 3 , 4 are independently distributed by a Gaussian distribution with mean θ s = 10 and standard deviation 2 in each run of the simulation, and s p ( k ) , p = 1 , 2 , 3 , 4 are independently and identically distributed complex random variables drawn from a random generator r a n d ( 0 , 1 ) . In these scenarios, the signal covariance matrix R s is no longer a rank-one matrix, and the output SINR should be rewritten as [35]
S I N R o p t = w H R s w w H R int + n w ,
which can be maximized by [35]
w o p t = P { R int + n 1 R s } ,
where P { · } stands for the principal eigenvector of a matrix, which is the eigenvector corresponding to the maximal eigenvalue. The DOAs ϑ 1 and ϑ 2 of the two interferers are uniformly distributed in [ 32 , 28 ] and [ 38 , 42 ] , respectively. In addition, it must be noted that ϑ 1 , ϑ 2 , and θ p , p = 1 , 2 , 3 , 4 change from run to run, but are kept fixed from snapshot to snapshot.
It can be seen, from Figure 6a, that the output SINR of the proposed beamformer was very close to the optimal SINR and its robustness against incoherent local scattering was better than the other tested methods. Figure 6b indicates the output SINR versus number of snapshots, which did not affect the performance of any of the algorithms. In general, the proposed algorithm has good robustness against mismatching due to incoherent local scattering.

4.6. Performance Comparison against Multiple Interferers

In the last example, the performance against multiple interferers is examined. Assumed that five interfering sources with random waveforms were uniformly distributed in [ 62 , 58 ] , [ 37 , 33 ] , [ 17 , 23 ] , [ 37 , 43 ] , [ 52 , 58 ] , respectively, and the INRs of five interferers were set as 40 dB, 20 dB, 25 dB, 30 dB, 35 dB, respectively. The desired signal with input SNR = 10 dB came from the direction 15 . The weight vectors of all tested beamformers could be obtained by broadening the nulls at the central angles.
The output SINR of the tested beamformers versus input SNR for the number of snapshots K = 60 was shown in Figure 7a, while Figure 7b compared the output SINR versus the number of snapshots at a SNR of 10 dB. The simulation results showed that the output SINRs of all tested algorithms decreased in the case of multiple interferers. However, the proposed beamformer effectively suppressed multiple interferers and obviously outperformed the other tested algorithms. It has been demonstrated that the number of interferers has little impact on the anti-jamming performance of the tested beamformers, and the proposed null broadening beamforming technique is robust against multiple interferers with different INRs.

5. Conclusions

The paper proposes a novel robust beamforming technique using null broadening for mitigating intereference caused by fast-moving jammers or antenna platform motion. The weight vector can be obtained by introducing the reconstructed INC matrix and the estimated the steering vector of the desired signal into the traditional MVDR beamformer. The DOAs of interfering signals can be estimated with the Capon spectrum estimator, and a new INC matrix can be reconstructed by placing several virtual interference sources around each of the actual jammers. The null depth can be controlled by the spatial power parameter in the process of reconstruction. Using convex optimization software, the steering vector of the desired signal can be estimated by solving a new optimization problem based on the reconstructed INC and SNC matrices. The simulation results demonstrated that the performance of the proposed null broadening beamformer dramatically outperforms other methods including classic robust adaptive beamforming and null broadening techniques, and can achieve good robustness against fast-moving jammers or antenna platform motion. Moreover, if the DOAs of the interferers or the direction of antenna platform motion is known, more virtual interference sources can be set along the actual interference DOAs, in order to achieve better performance.

Author Contributions

Conceptualization, J.Y. and J.L.; Methodology, J.Y. and J.L.; software, J.Y. and J.L.; validation, J.Y. and X.L.; formal analysis, J.Y.; investigation, J.Y. and J.L.; data curation, J.Y. and J.L.; writing–original draft preparation, J.Y. and J.L.; writing–review and editing, J.Y., J.L., and G.L.; visualization, X.L.; supervision, G.L.; project administration, J.Y.; funding acquisition, J.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 61501471.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, C.; Gu, Y.; Fan, Y.; Shi, Z.; Mao, Z.; Zhang, Y.D. Direction-of-Arrival estimation for coprime array via virtual array interpolation. IEEE Trans. Signal Process. 2018, 66, 5956–5971. [Google Scholar] [CrossRef]
  2. Shi, Z.; Zhou, C.; Gu, Y.; Goodman, N.A.; Qu, F. Source estimation using coprime array: A sparse reconstruction perspective. IEEE Sens. J. 2017, 17, 755–765. [Google Scholar] [CrossRef]
  3. Singh, R.-K.; Basu, A.; Koul, S.-K. Efficient Null Broadening and Steering Using Slot Antenna Array for Radar Applications. In Proceedings of the 2016 Asia-Pacific Microwave Conference (APMC), New Delhi, India, 5–9 December 2016; pp. 1–4. [Google Scholar]
  4. Zheng, Z.; Zheng, Z.; Wang, Z.; Zhang, Z. Covariance Matrix Reconstruction With Interference Steering Vector and Power Estimation for Robust Adaptive Beamforming. IEEE Trans. Veh. Technol. 2018, 67, 8495–8502. [Google Scholar] [CrossRef]
  5. Gu, Y.; Leshem, A. Robust adaptive beamforming based on interference covariance matrix reconstruction and steering vector estimation. IEEE Trans. Signal Process. 2012, 60, 3881–3885. [Google Scholar]
  6. Gu, A.; Goodman, N.A.; Hong, S.; Li, Y. Robust adaptive beamforming based on interference covariance matrix sparse reconstruction. Signal Process. 2014, 96, 375–381. [Google Scholar] [CrossRef]
  7. Zhou, C.; Gu, Y.; He, S.; Shi, S. A robust and efficient algorithm for coprime array adaptive beamforming. IEEE Trans. Veh. Technol. 2018, 67, 1099–1112. [Google Scholar] [CrossRef]
  8. Zhang, Z.; Liu, W.; Leng, W.; Wang, A.; Shi, H. Interference-plus-noise covariance matrix reconstruction via spatial power spectrum sampling for robust adaptive beamforming. IEEE Signal Process. Lett. 2016, 23, 121–125. [Google Scholar] [CrossRef] [Green Version]
  9. Shen, F.; Chen, F.; Song, J. Robust adaptive beamforming based on steering vector estimation and covariance matrix reconstruction. IEEE Commun. Lett. 2015, 19, 1636–1639. [Google Scholar] [CrossRef]
  10. Ma, Y.-X.; Lu, D.; Wang, W.-Y.; Wang, L.; Wu, R.-B. A high-dynamic null-widen GPS anti-jamming algorithm based on statistical model of the changing interference DOA. In China Satellite Navigation Conference (CSNC) 2014 Proceedings: Volume I; Springer: Berlin/Heidelberg, Germany, 2014; pp. 694–702. [Google Scholar]
  11. Mailloux, R.-J. Covariance matrix augmentation to produce adaptive array pattern troughs. Electron. Lett. 1995, 31, 771–772. [Google Scholar] [CrossRef]
  12. Guerci, J.-R. Theory and application of covariance matrix tapers for robust adaptive beamforming. IEEE Trans. Signal Process. 1999, 47, 977–985. [Google Scholar] [CrossRef]
  13. Zatman, M.; Guerci, J.-R. Comment on theory and application of covariance matrix tapers for robust adaptive beam forming. IEEE Trans. Signal Process. 2000, 48, 1796–1800. [Google Scholar] [CrossRef]
  14. Rubsamen, M.; Gerlach, C.; Gershman, A.-B. Low-rank covariance matrix tapering for robust adaptive beamforming. In Proceedings of the 2008 IEEE International Conference on Acoustics, Speech and Signal Processing, Las Vegas, NV, USA, 31 March–4 April 2008; pp. 2333–2336. [Google Scholar]
  15. Zhou, Z.; Zhang, Z.; Lai, S.-M. An Analytical Approach to Null Broadening. In Proceedings of the 2018 IEEE International Conference on Computer and Communication Engineering Technology (CCET), Beijing, China, 18–20 August 2018; pp. 181–185. [Google Scholar]
  16. Turner, J.-S. New directions in communications. IEEE J. Sel. Areas Commun. 1995, 13, 11–23. [Google Scholar]
  17. Qian, J.; He, Z.; Zhang, Y.; Wu, B. Null Broadening Adaptive Beamforming Based on Semidefinite Programming. In Proceedings of the 2016 IEEE 13th International Conference on Signal Processing (ICSP), Chengdu, China, 6–10 November 2016; pp. 1601–1604. [Google Scholar]
  18. Qian, J.; He, Z.; Xie, J.; Zhang, Y. Null Broadening Adaptive Beamforming Based on covariance matrix reconstruction and similarity constraint. EURASIP J. Adv. Signal Process. 2017, 2017, 1–10. [Google Scholar] [CrossRef] [Green Version]
  19. Song, H.; Kuperman, W.-A.; Hodgkiss, W.-S. Null broadening with snapshot-deficient covariance matrices in passive sonar. IEEE J. Ocean. Eng. 2003, 28, 250–263. [Google Scholar] [CrossRef]
  20. Li, W.-X.; Mao, X.-J.; Sun, Y.-X. A new algorithm for null broadening beamforming. J. Electron. Inf. Technol. 2014, 36, 2882–2888. [Google Scholar]
  21. Liu, Z.; Su, Z.; Hu, Q. Robust Sidelobes Cancellation Algorithm with Null Broadening. J. Electron. Inf. Technol. 2016, 38, 565–570. [Google Scholar]
  22. Yang, X.; Li, S.; Long, T.; Sarkar, T.-K. Adaptive null broadening method in wideband beamforming for rapidly moving interference suppression. Electron. Lett. 2018, 54, 1003–1005. [Google Scholar] [CrossRef]
  23. Li, W.; Yang, B. An Improved Null Broadening Beamforming Method based on Covariance Matrix Reconstruction. In Proceedings of the 2017 International Applied Computational Electromagnetics Society Symposium—Italy (ACES), Firenze, Italy, 26–30 March 2017; pp. 1–2. [Google Scholar]
  24. Ors, B.; Suleesathira, R. First and Second Order Iterative Null Broadening Beamforming. In Proceedings of the 2019 3rd International Conference on Imaging, Signal Processing and Communication (ICISPC), Singapore, 27–29 July 2019; pp. 47–51. [Google Scholar]
  25. Xiao, Z.; Lu, Y. Data-Based Model for Wide Nulling Problem in Adaptive Digial Beamforming Antenna Array. IEEE Antennas Wirel. Propag. Lett. 2019, 18, 2249–2253. [Google Scholar] [CrossRef]
  26. Liu, Z.; Zhao, S.; Zhang, G.; Jiao, B. Robust Adaptive Beamforming for Sidelobe Canceller with Null Widening. IEEE Sens. J. 2019, 19, 11213–11220. [Google Scholar] [CrossRef]
  27. Zhao, Y.; Li, W.; Mao, X.; Zhang, N. A null broadening beamforming against array calibration errors. J. Harbin Eng. Univ. 2018, 39, 163–168. [Google Scholar]
  28. Yang, X.; Li, S.; Sun, Y.; Long, T.; Sarkar, T.-K. Robust Wideband Adaptive Beamforming with Null Broadening and Constant Beamwidth. IEEE Trans. Antennas Propag. 2019, 67, 5380–5389. [Google Scholar] [CrossRef]
  29. Mohammadzadeh, S.; Kukrer, O. Robust adaptive beamforming for fast moving interference based on the covariance matrix reconstruction. IET Signal Proc. 2019, 13, 486–493. [Google Scholar] [CrossRef]
  30. Wang, H.; Liu, G.; Fan, Z.; Yao, Z. A null widening anti-jamming algorithm for GNSS receivers. J. Harbin Ins. Technol. 2019, 51, 94–98. [Google Scholar]
  31. Cong, Y.; Feng, D.; Li, H. High Dynamic GNSS Anti-jamming Algorithms Based on Nulling Widening and Deepening. J. Jilin Univ. 2020, 38, 1–8. [Google Scholar]
  32. Capon, J. High-resolution frequency-wavenumber spectrum analysis. Proc. IEEE 1969, 57, 1408–1418. [Google Scholar] [CrossRef] [Green Version]
  33. Gass, S.I.; Harris, C.M. Karush-Kuhn-Tucker (KKT) conditions. In Encyclopedia of Operations Research and Management Science; Gass, S.I., Harris, C.M., Eds.; Springer: New York, NY, USA, 2001. [Google Scholar] [CrossRef]
  34. Grant, M.; Boyd, S.; Ye, S. CVX: Matlab Software for Disciplined Convex Programming, Version 2.1. Mar. 2017. Available online: http://cvxr.com/cvx (accessed on 16 May 2019).
  35. Vorobyov, S.A.; Gershman, A.B.; Luo, Z.-Q. Robust adaptive beamforming using worst-case performance optimization: A solution to the signal mismatch problem. IEEE Trans. Signal Process. 2003, 51, 313–324. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flow of the proposed beamformer algorithm.
Figure 1. Flow of the proposed beamformer algorithm.
Sensors 20 01865 g001
Figure 2. Beampattern of null broadening: (a) Array pattern at θ 1 = 30 ; (b) Null broadening at θ 1 ; and (c) Null broadening at θ 2 = + 40 .
Figure 2. Beampattern of null broadening: (a) Array pattern at θ 1 = 30 ; (b) Null broadening at θ 1 ; and (c) Null broadening at θ 2 = + 40 .
Sensors 20 01865 g002
Figure 3. Interference Suppression: (a) Output signal-to-interference-plus-noise ratio (SINR) versus input signal-to-noise ratio (SNR); and (b) Output SINR versus snapshot number.
Figure 3. Interference Suppression: (a) Output signal-to-interference-plus-noise ratio (SINR) versus input signal-to-noise ratio (SNR); and (b) Output SINR versus snapshot number.
Sensors 20 01865 g003
Figure 4. Performance against direction of arrival (DOA) error: (a) Output SINR versus input SNR; and (b) Output SINR versus snapshot number.
Figure 4. Performance against direction of arrival (DOA) error: (a) Output SINR versus input SNR; and (b) Output SINR versus snapshot number.
Sensors 20 01865 g004
Figure 5. Performance against array geometry errors: (a) Output SINR versus input SNR; and (b) Output SINR versus snapshot number.
Figure 5. Performance against array geometry errors: (a) Output SINR versus input SNR; and (b) Output SINR versus snapshot number.
Sensors 20 01865 g005
Figure 6. Performance against Incoherent Local Scattering: (a) Output SINR versus input SNR; and (b) Output SINR versus snapshot number.
Figure 6. Performance against Incoherent Local Scattering: (a) Output SINR versus input SNR; and (b) Output SINR versus snapshot number.
Sensors 20 01865 g006
Figure 7. Performance comparsion against multiple interferers: (a) Output SINR versus input SNR; and (b) Output SINR versus snapshot number.
Figure 7. Performance comparsion against multiple interferers: (a) Output SINR versus input SNR; and (b) Output SINR versus snapshot number.
Sensors 20 01865 g007
Table 1. Null widths at different beampattern gains.
Table 1. Null widths at different beampattern gains.
Beampattern Gain
Beamformer−50 dB−100 dB−150 dB
Rec-ISVPE 2.4 <0.1 <0.1
Laplace 6.9 <0.1 can not reach
Mailloux 6.4 <0.1 <0.1
Rec-MVDR 15.2 6.9 4.2
AANB 5.8 can not reachcan not reach
Proposed method 15.1 6.9 4.2

Share and Cite

MDPI and ACS Style

Yang, J.; Lu, J.; Liu, X.; Liao, G. Robust Null Broadening Beamforming Based on Covariance Matrix Reconstruction via Virtual Interference Sources. Sensors 2020, 20, 1865. https://doi.org/10.3390/s20071865

AMA Style

Yang J, Lu J, Liu X, Liao G. Robust Null Broadening Beamforming Based on Covariance Matrix Reconstruction via Virtual Interference Sources. Sensors. 2020; 20(7):1865. https://doi.org/10.3390/s20071865

Chicago/Turabian Style

Yang, Jian, Jian Lu, Xinxin Liu, and Guisheng Liao. 2020. "Robust Null Broadening Beamforming Based on Covariance Matrix Reconstruction via Virtual Interference Sources" Sensors 20, no. 7: 1865. https://doi.org/10.3390/s20071865

APA Style

Yang, J., Lu, J., Liu, X., & Liao, G. (2020). Robust Null Broadening Beamforming Based on Covariance Matrix Reconstruction via Virtual Interference Sources. Sensors, 20(7), 1865. https://doi.org/10.3390/s20071865

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop