International Journal of Antennas and Propagation

Volume 2017 (2017), Article ID 3265236, 9 pages

https://doi.org/10.1155/2017/3265236

## Adaptive Antenna Null Broadening Beamforming against Array Calibration Error Based on Adaptive Variable Diagonal Loading

^{1}College of Information and Communications Engineering, Harbin Engineering University, Harbin, China^{2}Department of Electronics, Carleton University, Ottawa, ON, Canada

Correspondence should be addressed to Yu Zhao; nc.ude.uebrh@6180uyoahz

Received 7 June 2017; Accepted 29 August 2017; Published 8 October 2017

Academic Editor: Ikmo Park

Copyright © 2017 Wenxing Li 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

An approach for null broadening beamforming is proposed based on adaptive variable diagonal loading (VDL) and combined with the covariance matrix taper (CMT) approach, aiming at improving the robustness of adaptive antenna null broadening beamforming when array calibration error exists. Hence, it is named VDL-CMT. In this novel approach, the signal-to-noise ratio in the tapered sample covariance matrix is estimated and the VDL factor can be obtained adaptively. Then, the covariance matrix of the CMT approach is loaded with the obtained VDL factor. According to simulation results, in the case of array calibration error, robustness of the VDL-CMT is significantly improved and its performance is better than that of the existing adaptive antenna null broadening beamforming approaches.

#### 1. Introduction

Adaptive antenna beamforming is a technique used to receive desired signals and suppress interferences by adjusting the weight vectors of arrays adaptively. It has been widely used in radar, wireless communication, radio astronomy, and many other areas [1–3]. The performance of the adaptive antenna beamforming will be severely degraded if the weight vector of an array cannot be able to adapt sufficiently fast to the jammer motion. Forming broad nulls around the directions of interferences is an effective means to settle this problem [4, 5].

There exist many adaptive antenna null broadening beamforming approaches in literature [6–14]. Covariance matrix taper (CMT) [6–9] is a classical and widespread approach of null broadening. The computational complexity of the CMT approach is extremely low, while the depth of broad nulls is shallow. Quadratic constraint sector suppressed (QCSS) [10, 11] is also a null broadening beamforming approach. The null depth of the QCSS is much deeper in comparison with CMT. However, the solving process of the QCSS approach is complicated. Linear constraint sector suppressed (LCSS) approach [12] is proposed on the basis of the QCSS approach. The nonlinear quadratic constraint is transformed into a set of linear constraints; hence the solving process is greatly simplified. A null broadening beamforming approach based on projection and diagonal loading (PDL) is proposed in [13], where the CMT and projection are combined. The null depth of the PDL approach is deeper than that of the CMT approach, and the computational complexity is low. The null broadening beamforming approach of covariance matrix reconstruction and similarity constraint (CMRSC) is proposed in [14]. The interference-plus-noise covariance matrix of array is reconstructed by the CMRSC. The depth of the broad nulls of CMRSC is much deeper than that of the CMT approach, but the computational complexity is high.

The approaches of LCSS [12], PDL [13], and CMRSC [14] are all proposed in recent years. Without array calibration error, deeper broad nulls and higher output signal-to-interference-plus-noise ratio (SINR) can be obtained by virtue of these approaches in comparison with the CMT approach. Nevertheless, all three approaches construct the correlation matrix of the steering vector based on the presumed steering vector, so they are greatly sensitive to array calibration error. In practice, the array calibration error may exist, such as the antenna location error and mutual coupling error [15]. In the case of the array calibration error, the performance of the above three adaptive antenna null broadening beamforming approaches will be severely degraded. On the contrary, the CMT approach [7] does not construct the correlation matrix of the steering vector. Therefore, its performance is better than that of the three approaches above in overcoming the array calibration error.

Diagonal loading (DL) [16] is a widespread robust beamforming technique. The robustness of adaptive antenna beamforming can be effectively improved by DL with a proper loading factor that is difficult to achieve [17, 18]. The CMT approach always uses a small fixed loading factor [9, 13, 14], so its robustness against the array calibration error is very limited. In order to improve the robustness of adaptive antenna null broadening beamforming against array calibration error, an adaptive variable diagonal loading based covariance matrix taper (VDL-CMT) null broadening beamforming approach is put forward. The signal-to-noise ratio (SNR) in the tapered sample covariance matrix is estimated and the variable diagonal loading (VDL) factor can be obtained adaptively. Then, the covariance matrix of the CMT approach is loaded with the obtained VDL factor, and the optimal weight vector of the adaptive antenna is calculated by the use of the loaded covariance matrix. The advancement and effectiveness of the VDL-CMT approach are verified by theoretical analysis and simulation results.

#### 2. The Existing Null Broadening Beamforming Approaches

##### 2.1. Signal Model

Consider a uniform linear array (ULA) with elements and the element spacing is equal to one-half wavelength. is the number of far-field narrowband excitation sources which are incoherent. The received data of the antenna array can be expressed aswhere is the matrix of array steering vectors, is the complex sources envelope, and is the noise of the antenna array.

The covariance matrix of the array can be written aswhere denotes the expectation operation, denotes the conjugate transposition, represents the matrix of sources, is the power of noise, and is an identity matrix of dimension.

The standard minimum variance distortionless response (MVDR) beamformer can be expressed as follows:where is the weight vector of the array, is the steering vector of the desired signal, and is the covariance matrix of interference-plus-noise. Due to the difficulty of directly calculating , this matrix is commonly replaced by the sample covariance matrix with limited number of snapshots aswhere is the number of snapshots. The optimal weight vector of the MVDR beamformer can be expressed aswhere denotes the inverse operation. The output SINR is used to measure the output performance of the adaptive beamforming, which can be expressed aswhere is the input power of the desired signal.

##### 2.2. The CMT Approach

CMT [7] is a classical and widespread null broadening beamforming approach. The sample covariance matrix is tapered by the matrix aswhere represents Hadamard product and denotes the tapered sample covariance matrix. The element of th row and th column of can be expressed aswhere , and determines the width of the broad nulls. Then, the optimal weight vector of the array can be calculated using .

The broad nulls can be formed around the interferences adaptively without the prior information about the directions of the interferences, and the computational complexity of the CMT approach is very low. However, a single interference is replaced by a set of virtual interferences with equal power to form the broad nulls, so the null depth of the CMT approach is shallow. This also leads to rising of sidelobes and deteriorates output performance.

##### 2.3. The LCSS Approach

The average output power of the predefined sector is constrained by the LCSS approach [12], which can be expressed aswhere is the predefined output power and is the correlation matrix of the steering vector that can be expressed as equidistant sampling points in sector are used to construct . The ideal goal to suppress the output power in sector is to obtain a zero power response; that is, . can be decomposed as based on its eigenvalues, where is the matrix of eigenvectors and is a diagonal matrix containing the eigenvalues of in a decreasing order. Assume that has a rank equal to , and is the matrix of eigenvectors corresponding to the large eigenvalues. The quadratic constraint can be achieved if , so the broad nulls can be formed through a set of linear constraints as follows:where , , and denotes transposition. The solution to (11) iswhere is the number of linear constraints that minimize the average output power of sector . The depth of the broad nulls can be set according to the number of linear constraints. In the case without the array calibration error, the performance of the LCSS approach is better than that of the CMT approach. However, the correlation matrix of the steering vector in the LCSS approach is constructed based on the presumed steering vector and is very sensitive to the array calibration error. Therefore, the performance of the LCSS approach will be severely degraded when the array calibration error exists.

##### 2.4. The PDL Approach

The broad nulls can be formed by the PDL approach [13] through projection. The correlation matrix of the steering vector is constructed in the predefined null broadening sectors aswhere is the number of the predefined null broadening sectors. Combined with the CMT approach, is tapered as

Based on its eigenvalues, is decomposed aswhere is the eigenvalues of , and is the corresponding eigenvectors. The projection matrix is constructed by the eigenvector corresponding to the large eigenvalues as

Then, we project onto the sample covariance matrix as

The optimal weight vector of the PDL approach can be calculated using . The correlation matrix of the steering vector is constructed according to the presumed steering vector as formula (13) and very sensitive to the array calibration error. Thus the performance of the PDL approach will be severely degraded in the case of array calibration error.

##### 2.5. The CMRSC Approach

The interference-plus-noise covariance matrix is reconstructed by the CMRSC approach [14] in the predefined sector , and the broad nulls can be formed as demonstrated below. Firstly, calculate the maximum power in sector as the interference power, which can be expressed as

Reconstruct the interference covariance matrix according to and the presumed steering vector as

Then, is decomposed based on its eigenvalues and the power of noise can be calculated according to the average value of the small eigenvalues corresponding to noise. The power of noise can be expressed aswhere is the number of sources. The interference-plus-noise covariance matrix is reconstructed as

The optimal weight vector can be calculated using . The interference-plus-noise covariance matrix of the CMRSC approach is reconstructed according to the presumed steering vector as formula (19), so its performance will be severely degraded in the case of the calibration error.

#### 3. The Proposed Approach

##### 3.1. The Method to Obtain the VDL Factor

It is generally known that the improvement in robustness of the CMT approach by DL with a fixed loading factor is limited. The main reason is displayed in this paper by analyzing the effect of DL on the components of the tapered sample covariance matrix , and then the VDL-CMT approach is proposed.

The covariance matrix of the CMT approach with DL can be expressed as [7, 9, 16]where is the sample covariance matrix, is an identity matrix of dimension, and is the loading factor that is expressed as the multiple of noise power . is decomposed based on its eigenvalues aswhere is the matrix of eigenvectors, is a diagonal matrix containing the eigenvalues, denotes the signal-plus-interference subspace, is a diagonal matrix containing the large eigenvalues corresponding to , denotes the noise subspace, and is a diagonal matrix containing the small eigenvalues corresponding to . can be expressed aswhere is the power of the desired signal in , is the power of interference, is the power of noise, denotes the desired signal component, denotes the interference component, and denotes the noise component.

Ideally, can be obtained by tapering the interference-plus-noise covariance matrix , and the performance of the beamforming will be extremely good. However, is commonly replaced by the sample covariance matrix in practical applications [19]. The output performance and robustness of the CMT approach will be degraded, because has two drawbacks compared with . Firstly, is obtained through limited number of snapshots, which causes the disturbance of small eigenvalues corresponding to the noise subspace. The sidelobes will rise and the performance of the beamforming will be degraded [16, 17]. Secondly, includes the desired signal component, which will degrade the performance of beamforming, especially in nonideal situations [19–21]. Therefore, if the disturbance of the small eigenvalues corresponding to the noise subspace is suppressed, the sidelobes will be depressed, and the output performance of beamforming will be improved. If the desired signal component or its proportion is decreased, the output performance and robustness of beamforming will be improved. The covariance matrix of the CMT approach with DL can be expressed as

It can be known from formula (25) that white noise is added to the covariance matrix by DL. The noise component is strengthened, exerting different impacts on each component of the covariance matrix. The disturbance of the small eigenvalues corresponding to the noise subspace is suppressed, so should be bigger to further suppress the disturbance of the small eigenvalues. The SNR in the covariance matrix is decreased. That is to say, the proportion of the desired signal component is decreased. Therefore, should be bigger to further reduce the SNR in the covariance matrix. The proportion of the interference component is decreased, which reduces the null depth, so should be smaller to prevent the null depth from being reduced. Therefore, the components of have various requirements for , and they are contradictory. The performance improvement for the beamforming is limited if is too small. On the contrary, the performance of the beamforming may be degraded if is too large and so a proper loading factor is difficult to achieve. The CMT approach always uses a small fixed loading factor [9, 13, 14]. In the following, we make analysis on the main reason why the robustness improvement of the CMT approach by DL with a fixed loading factor is limited.

The SNR in is

The SNR in is

It can be obtained from formula (27) that if is fixed, will always increase with the increase of , so the SNR in the covariance matrix cannot be effectively decreased. That is the main reason why the robustness improvement of the CMT approach by DL with a fixed loading factor is limited. If can be adjusted adaptively with the change of , the performance of the CMT approach may be further improved. Thus the VDL-CMT approach is proposed in this paper, where the VDL factor is adjusted adaptively according to , and the SNR in the covariance matrix can be effectively decreased.

The minimum loading factor that can effectively decrease the SNR in the covariance matrix and suppress the disturbance of the small eigenvalues corresponding to the noise subspace is selected as the VDL factor for the VDL-CMT approach. When there is Gaussian white noise, the small eigenvalues disturb around the noise power . The disturbance of the small eigenvalues can be partly suppressed when . As the lower bound, 2 is a widely used loading factor, which can be used in both the MVDR beamforming and the CMT approach. In order to effectively decrease the SNR in the covariance matrix to a certain degree even if the input SNR is very high, can be adjusted adaptively according to as

Formula (28) is obtained based on qualitative analysis. As the lower bound, when is equal to , the SNR in can be expressed as

It can be observed from formula (29) that, no matter how changes, there is always and . That is to say, the power of the desired signal in the covariance matrix is close to that of the noise. The SNR in the covariance matrix can be decreased effectively.

From the above analysis, we learn that the SNR in the covariance matrix can be effectively decreased when , and the disturbance of the small eigenvalues corresponding to the noise subspace can be suppressed when . The minimum VDL factor under the two constraints is a compromising for the VDL-CMT approach, which can be expressed as

The VDL factor for the VDL-CMT approach is adjusted adaptively according to , by which the SNR in the covariance matrix can be effectively decreased. Thus the robustness of the VDL-CMT approach is improved compared with the CMT approach using a fixed loading factor. It is important to emphasize that is the SNR in the tapered sample covariance matrix , but not the input SNR of the array. Then, should be estimated to get .

##### 3.2. The Method to Estimate

There is no existing method to estimate . In this paper, is obtained based on the relationship between and the SNR in the sample covariance matrix , which will be derived below.

Firstly, the SNR in is estimated according to literature [22]. Similar to formula (23), is decomposed based on its eigenvalues as

The presumed steering vector of the desired signal is , which is projected to as is sorted in a decreasing order as . The eigenvectors and the eigenvalues are all sorted corresponding to as and . The maximum of the projections is obtained when is the eigenvector corresponding to the desired signal. The SNR in can be expressed as [22]

Then, we derive the relationship between and . According to formulas (18) and (19) in literature [9], can be expressed aswherewhere . When and , no broad nulls will be formed. When , the broad nulls will be formed with width , and . Phase “dithered” is introduced into the steering vectors of through by the CMT approach [9]. According to formula (35), the original single phase becomes phases, which realizes the “dithered” phase. Therefore, the original single interference is replaced by virtual interferences with equal power, and broad nulls can be formed with the width . The desired signal is affected the same as interferences, so the relationship between and is can be obtained according to formulas (33) and (36), so the VDL factor for the adaptive VDL-CMT approach can be obtained according to and formula (30).

##### 3.3. The VDL-CMT Approach

Solution procedures of the VDL-CMT approach are summarized as below.

*Step 1. *Obtain the tapered sample covariance matrix based on formulas (7) and (8).

*Step 2. *Estimate according to formulas (31), (32), and (33).

*Step 3. *Estimate according to formula (36).

*Step 4. *Obtain according to formula (30).

*Step 5. *Obtain the loaded covariance matrix according to and formula (25).

*Step 6. *Calculate the optimal weight vector using .

The VDL-CMT approach can be done by the above steps. The performance of the VDL-CMT approach is examined by simulations next.

#### 4. Simulation Results

Assume that the incoherent narrowband excitation sources are in the far field from the antenna array. We take account of a ULA with 12 array elements, and the element space is equal to one-half wavelength. The input SNR of array is 5 dB. The input interference to noise ratio is 30 dB. Additive noise is modeled as Gaussian white noise. The desired signal illuminates on the antenna array in the direction of . The interferences are from the directions of and . Assume that the width of the broad nulls is set to . The number of snapshots is 100. The VDL-CMT approach is compared with the LCSS approach [12], the PDL approach [13], the CMRSC approach [14], and the CMT approach [7] using the fixed loading factor 2. Four examples are used to verify the performance of the VDL-CMT approach. In all simulations, 100 Monte Carlo runs are used to obtain the results.

##### 4.1. In the Case without Array Calibration Error

The presumed array steering vector is the same as the actual situation. Figure 1 shows the beam patterns obtained by the five approaches. Figure 2 shows the output SINR versus the input SNR for the five approaches and the theoretical optimal result when the input SNR changes.