Advances in Antenna Array Processing for Radar 2016View this Special Issue
Robust Adaptive Beamforming Using a Low-Complexity Steering Vector Estimation and Covariance Matrix Reconstruction Algorithm
A novel low-complexity robust adaptive beamforming (RAB) technique is proposed in order to overcome the major drawbacks from which the recent reported RAB algorithms suffer, mainly the high computational cost and the requirement for optimization programs. The proposed algorithm estimates the array steering vector (ASV) using a closed-form formula obtained by a subspace-based method and reconstructs the interference-plus-noise (IPN) covariance matrix by utilizing a sampling progress and employing the covariance matrix taper (CMT) technique. Moreover, the proposed beamformer only requires knowledge of the antenna array geometry and prior information of the probable angular sector in which the actual ASV lies. Simulation results demonstrate the effectiveness and robustness of the proposed algorithm and prove that this algorithm can achieve superior performance over the existing RAB methods.
Aiming at receiving a signal from a certain direction and suppressing interferences and noise, adaptive beamforming has found widespread application in many fields ranging from radar, sonar, wireless communication, and radio astronomy to medical imaging, speech processing, and so forth [1–4]. Beamformers can be regarded as spatial filters, which have proper performance under ideal circumstances. The standard Capon beamformer (SCB) is an optimal spatial filter that maximizes the array output signal-to-interference-plus-noise ratio (SINR) . However, it is sensitive to array steering vector (ASV) uncertainty and direction of arrival (DOA) estimation error for the desired signal (DS). The performance also degrades due to the presence of the DS component in the training data and small sample size. In order to improve the robustness of SCB, various robust adaptive beamforming (RAB) techniques have been proposed in the past decades .
The existing RAB techniques mainly consist of diagonal loading (DL) technique , eigenspace-based (ESB) technique , uncertainty set-based technique [9, 10], and the reconstruction-based technique . They have effective performance in the presence of different mismatches in typical situations. However, there exist some drawbacks that limit their application range, which include the difficulty to choose the optimal DL factor, the high probability of subspace swap for ESB beamformers at low signal-to-noise ratio (SNR), the cost of reducing output SINR due to the expansion of the uncertainty set when the SV mismatch is large, the ad hoc nature, and the high computational cost.
Recently, a new approach to the design of RAB based on the interference-plus-noise (IPN) covariance matrix reconstruction has been introduced in . This method utilizes the Capon spectrum to integrate over an angular sector separated from the direction of DS. The reconstruction-based beamformer is effective at removing the DS from the sample covariance matrix but suffers from its high computational complexity. The main computational cost is due to the large number of samples involved in both spectrum estimation and summation. On the other hand, the recent beamformers in [12, 13] enhance the robustness against ASV errors. However, the utilization of complicated optimization software restricts their applications in practical situations. To improve the reconstruction-based RAB technique, a Low-Complexity Shrinkage-Based Mismatch Estimation (LOSCME) algorithm is presented in . But it can perform effectively only when the interfering sources are weak.
To reduce the computational complexity of the RAB method in  and eliminate its dependence on the optimization software, in this paper, a novel reconstruction-based RAB algorithm is proposed. This algorithm is characterized by lower complexity and a closed-form formula estimation of the actual ASV. Moreover, it requires very little prior information and has a superior performance to previously reported RAB algorithms. The only prior information required is the knowledge of the antenna array geometry and the coarse angular sector in which the actual ASV lies. Three steps are needed to achieve this algorithm. Firstly, a subspace-based method is introduced to obtain a closed-form formula for ASV estimation, avoiding using the optimization software. Then, the IPN covariance matrix is reconstructed based on a spatial power spectrum sampling (SPSS) method , which is realized by a proposed sample equation. Finally, the covariance matrix taper (CMT) technique  is utilized to adopt the relatively small size of array systems in practice. Simulation results will be provided to prove the effectiveness and robustness of the proposed beamformer.
The rest of this paper is organized as follows. The data model of array output and some necessary backgrounds about adaptive beamformer are described in Section 2. In Section 3, a novel RAB algorithm is proposed by employing ASV estimation and IPN covariance matrix reconstruction. Simulation results are presented in Section 4. Finally, conclusion is drawn in Section 5.
2. Signal Model
Consider a uniform linear array (ULA) of omnidirectional antenna elements impinged by narrowband uncorrelated far-field signals. The vector representing the received signal at the th snapshot can be modeled as where , , and denote the statistically independent vectors of the DS, interference, and noise, respectively. The beamformer output is given by , where is the complex beamformer weighting vector, and and denote the Hermitian transpose and transpose, respectively. The beamformer output SINR is defined aswhere is the assumed DOA of the DS and is the ASV, which has the general form . Meanwhile, is the power of the DS, and is the IPN covariance matrix.
By maximizing the output SINR of the beamformer, the optimal weight vector can be obtained byThe solution is known as the Capon beamformer :
In practice, theoretical covariance matrix is usually unavailable and sample covariance matrix (5) is used as an approximation:where is the number of data snapshots.
3. Proposed Approach
The performance of standard beamformers degrades dramatically when the ASV errors exist and the DS with a high SNR is present in the training snapshots. To remove the DS from the sample covariance matrix, recently, an IPN covariance matrix reconstruction method was proposed .
In , the IPN covariance matrix can be reconstructed by using the spatial spectrum of the array aswhere is the Capon spectrum , is the complement sector of and stands for the assumed direction range of the DS; that is, , and covers whole spatial domain.
This method has mainly two aspects of drawbacks. Firstly, concerning the mismatch between the actual ASV and the nominal ASV, the IPN matrix may not be reconstructed accurately. Secondly, its high computational complexity restricts its practical performance . In the following part, more precise estimation of the actual ASV can be achieved and the computational cost of IPN matrix reconstruction can be reduced.
3.1. Desired Signal Array Steering Vector Estimation
Similar to the reconstruction of IPN covariance matrix as (6), a new matrix can be constructed bywhere denotes a probability density function, but the choice of is more flexible, which means that can be chosen to be independent of , for example, , or adjust values depending on the prior information. Then, can be eigendecomposed aswhere are the eigenvalues of in descending order, are diagonal matrices, are the eigenvectors associated with , and . Following the approach in , let be the smallest integer such that , where is a predetermined threshold. Then, the eigenvectors associated with largest eigenvalues of can be chosen to form the column orthogonal matrix . As derivation in , the actual ASV lies in the subspace spanned by the columns of .
Similar to (8), in order to obtain the eigenvector, the sample covariance matrix can be decomposed aswhere are the eigenvalues of in descending order, and are diagonal matrices, are the eigenvectors associated with , and and denote the signal-plus-interference (SPI) subspace and noise subspace, respectively. It is well known that the actual ASV of DS belongs to the subspace spanned by the columns of .
As mentioned above, two constraints can be imposed on :where and are the coefficient vectors. Then, the actual ASV should be a vector lying within the intersection of . This intersection can be obtained by employing the theorem of sequential vector space projections , and the ASV of DS is proved to be estimated bywhere , , and denotes the eigenvector associated with the maximal eigenvalue of a matrix. This estimation method can be more efficient when the look direction error is large, and the method avoids the need of optimization software by using a closed-form formula.
3.2. Spatial Power Spectrum Sampling Based IPN Covariance Matrix Reconstruction
The main computational cost of the method in  is the integration approximation by summation, where (the number of sampled values) times spectrum estimation and vector multiplication operations have to be performed. To reduce the cost, the complex spectrum estimation process should be eliminated.
From the derivation above, can be regarded as an inverse discrete Fourier transform (IDFT) of an -point rectangular function in the frequency domain. When is large enough, can be approximated as a sinc function; that is, ; will approximate a Kronecker delta function; that is,The function is called the selecting property of the steering vector in .
Denote the zeros of by , and consider (14); is obtained when ; then there are such values in the set denoted by . Then can be easily obtained by .
Define a matrix using :where is a specified angular sector. Consider that can be formed by integrating the spatial spectrum though the whole region asThen, let ; (18) can be obtained asAs discussed above, when is large enough, is a Hermitian matrix and can be approximated as
In this way, the IPN covariance matrix is estimated by the sample matrix without calculating the spatial power spectrum. In practice, can be replaced by , which yields .
However, when , there will be a large estimation error, because the sampling spacing is not dense enough. For the purpose of improving the performance, the covariance matrix tapering technique introduced in  is employed. The tapered matrix is given by , where “” denotes the Hadamard product and is the taper matrix. Here, the Mailloux-Zatman (MZ) taper is used, which is defined aswhere corresponds to the width of the dithering area. Hence, should be adopted to taper the sample covariance matrix and estimated IPN covariance matrix, respectively; that is, the reconstruction of IPN covariance matrix can be obtain by
3.3. The Proposed Beamforming Algorithm
Based on the discussions above, the proposed beamforming algorithm can be summarized by the following steps:
Step 1. Construct the two subspaces and by eigendecomposing and , respectively. Then estimate the ASV of DS using (11).
Step 3. Substituting the reconstructed IPN covariance matrix and estimated ASV of DS into the Capon beamformer, the weight vector of the proposed approach is computed as
In the proposed algorithm, the main computational complexity lies in the eigendecomposition operation and the matrix inversion operation, both of which have complexity of . That means the overall computational complexity is . Considering that the number of sampling points required in (7) and (16) is reasonably much less than that required in the beamformer in , the proposed beamformer has a lower cost than that of . Additionally, compared to the existing RAB methods in [12, 13], which have complexity equal or higher than , the proposed algorithm also avoids the need of optimization software and, thus, is easier to apply in practice.
In the simulations, without loss of generality, a ULA with omnidirectional sensors is considered. It is assumed that there is one DS from the presumed direction . Two uncorrelated interferences with random waveforms arrive from DOA angles of and , respectively. The noise is modeled as zero-mean and unity variance spatially and temporally white Gaussian noise. The interference-to-noise ratio (INR) at each sensor is set to be fixed at 30 dB. For each scenario, 500 Monte Carlo trials are performed.
The proposed beamformer is compared to the diagonal loading sample matrix inversion (LSMI) beamformer , the ESB beamformer , the worst-case-based beamformer , the sequential quadratic programming (SQP) based beamformer , the reconstruction-based beamformer , and the beamformer in . The diagonal loading factor of the LSMI beamformer is chosen as twice the noise power. The dimension of the signal-plus-interference subspace is assumed to be always estimated correctly for the eigenspace-based beamformer. The value is selected for the worst-case-based beamformer as it has been recommended in , while the value and 20 dominant eigenvectors of the matrix are used in the SQP based beamformer. For the proposed beamformer, the reconstruction-based beamformer, and the beamformer in , the angular sector of the DS is presumed to be . In this paper, , , and are chosen in (11), (12), and (20), respectively.
Example 1 (mismatch due to signal direction error). The look direction mismatch of the DS is assumed to be random and uniformly distributed in for each simulation run, which means the mismatch changes from run to run but keeps fixed from snapshot to snapshot. The array normalized beam patterns for the beamformers in a single simulation run where and are displayed in Figure 1. It can be seen that, in the presence of pointing error, the proposed beamformer and ESB beamformer are able to point the main lobe to the actual direction, while the LSMI beamformer and worst-case-based beamformer cannot. Furthermore, the proposed beamformer has lower side lobes and deeper nulls at the directions of interferences, which make the proposed beamformer suppress the noise and interferences effectively.
Considering the influence of random look direction error on array output SINR, the performance curves versus the input SNR and the number of snapshots are drawn in Figures 2(a) and 2(b), respectively. The number of snapshots is fixed to be and the INR is kept at 30 dB in Figure 2(a), while a fixed SNR for the DS at 10 dB and a fixed INR at 30 dB are considered in Figure 2(b). From the results shown in Figure 2, it can be found that the output SINR of the proposed algorithm is closer to the optimal SINR in a large range of SNR from to dB and for all values of number of snapshots from 10 to 100, which illustrates its high dynamic range. That means the proposed algorithm outperforms the other tested beamformers in the scenario where only signal direction error exists.
The SINR performance versus the pointing error is also investigated, and the results are shown in Figure 3 with fixed SNR and INR at 10 dB and 30 dB, respectively. It can be see that the proposed beamformer effectively deals with a wide range of pointing errors and achieves the best performance among the tested beamformers even when the pointing error is as large as ±4°.
Example 2 (mismatch due to ASV random error). In this simulation example, the ASV of DS is assumed to be randomly distributed in an uncertainty set, which can be modeled aswhere represents the nominal SV corresponding to the direction and is the random error vector, which is drawn in each simulation run independently from an uncertainty set as follows :where and are uniformly distributed in the intervals and , respectively. Figures 4(a) and 4(b) correspond to the performance of the investigated methods versus the input SNR and the number of snapshots. It can be seen from the figures that the proposed beamformer and the reconstruction-based beamformer outperform the other tested beamformers. Table 1 shows the deviations from the optimal SINR for the beamformers at and , respectively. The performance improvement is a direct result of the ASV estimation and DS elimination, especially at high SNR values.
Furthermore, the proposed algorithm performs almost as well as the reconstruction-based beamformer for the output SINR but enjoys a faster convergence rate because of the lower computational cost without complex integral computation. Since the ASV mismatch is comprehensive and arbitrary-type, the proposed beamformer is proved to be effective against the random error of ASV.
Example 3 (mismatch due to incoherent local scattering). In this example, it is assumed that the desired signal has a time-varying ASV which can be modeled as where , without loss of generality, , denotes the incoherently scattered signal paths. Assume that the directions are randomly distributed in a Gaussian distribution with mean and standard deviation ; and are independently and identically distributed complex Gaussian random variables drawn from which change from snapshot to snapshot.
In this scenario, the signal covariance matrix is no longer a rank-one matrix; thus the output SINR should be rewritten in a more general form :which is maximized by where denotes the same operation as in (11). Figures 5(a) and 5(b) show the performance curves versus the input SNR and the number of snapshots. It can be seen from the figures that the proposed beamformer presents an effective performance when handling incoherent local scattering. Similar to Example 2, in detail, there is about 0.8 dB performance loss for the proposed algorithm comparing to the reconstruction-based beamformer. The main reason is that the DS may leak into the complement sector due to the incoherent local scattering, which affects the accuracy of the ASV estimation.
Example 4 (impacts of some factors on performance). The main purpose of this example is to study the impacts of some factors on performance. and are two main factors in the proposed algorithm. The former affects the accuracy of the construction of subspace in which the actual ASV lies, and the latter decides the choice of taper matrix which plays a key role in the reconstruction of IPN covariance matrix.
For the purpose of studying the impacts of the two factors, the model of the mismatch is set to be the same as the first example. The number of snapshots is fixed to be ; the INR and SNR are kept at 20 dB and 10 dB, respectively. The performances of the proposed algorithm versus and are displayed in Figures 6(a) and 6(b), respectively. From Figure 6(a), it can be seen that the output SINR increases as gets larger, but the growth range is quite small. Considering that the computational cost will not significantly increase as gets larger, setting in the proposed algorithm is appropriate. Similarly, the result of the test of is shown in Figure 6(b), in which the choice of can be found near the peak of the performance curve. To ensure a satisfying output performance, can be decided by a test for a typical array before the algorithm operates. The procedure is offline and will not dramatically increase the system cost.
In this paper, a novel low-complexity RAB method is proposed which is easier to realize in practical situations and more robust to the look direction mismatch than other existing algorithms. The ASV is estimated by a closed-form formula so as to avoid utilizing the optimization software, and the IPN covariance matrix is reconstructed by a sampling progress. The proposed beamformer only requires prior knowledge of the antenna geometry and the angular sector in which the ASV is located. Simulation results have demonstrated that the proposed beamformer can achieve superior performance over the existing state of the art RAB methods. To simplify the illustration, the influence of the element pattern, the polarization, and the mutual coupling is not considered in this paper. However, these elements will be investigated in the future study.
The authors declare that they have no competing interests.
H. L. Van Trees, Optimum Array Processing, John Wiley & Sons, Hoboken, NJ, USA, 2002.
W. Liu and S. Weiss, Wideband Beamforming: Concepts and Techniques, John Wiley & Sons, Chichester, UK, 2010.