Next Article in Journal
Determination of Excitation Amplitude and Phase for Wide-Band Phased Array Antenna Based on Spherical Wave Expansion and Mode Filtering
Next Article in Special Issue
k-Level Extended Sparse Array Design for Direction-of-Arrival Estimation
Previous Article in Journal
A 0.00426 mm2 77.6-dB Dynamic Range VCO-Based CTDSM for Multi-Channel Neural Recording
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Oblique Projection-Based Covariance Matrix Reconstruction and Steering Vector Estimation for Robust Adaptive Beamforming

Guangxi Key Laboratory of Wireless Wideband Communication & Signal Processing, Guilin University of Electronic Technology at Guilin, Guilin 541004, China
*
Author to whom correspondence should be addressed.
Electronics 2022, 11(21), 3478; https://doi.org/10.3390/electronics11213478
Submission received: 13 September 2022 / Revised: 16 October 2022 / Accepted: 24 October 2022 / Published: 26 October 2022
(This article belongs to the Special Issue Recent Advances and Applications of Array Signal Processing)

Abstract

:
Adaptive beamforming can efficiently contract interference and noise. Due to high sensitivity of the beamformer to model mismatch, the capability of interference reduction will critically degrade when the signal model mismatch occurs, particularly when the sampling sequence contains the desired signal. For the purpose of enhancing the robustness of beamformers to signal model mismatch, we propose a new robust adaptive beamforming (RAB) method. Firstly, the precise steering vector (SV) associating with the desired signal is estimated by employing the minimum norm of subspace projection (MNSP) approach. Secondly, the nominal interference SVs are estimated via the maximum entropy power spectrum. Subsequently, the corrected interference SVs and powers are obtained by oblique projection. Finally, the interference-plus-noise covariance matrix (INCM) is reconstructed, and the proposed RAB is obtained. Multiple simulations are carried out and demonstrate the robustness of the proposed RAB method.

1. Introduction

Adaptive beamforming is a spatial filter which is designed to form directional beams to transmit or receive a desired signal from a designated direction, meanwhile attenuating interference signals from other directions. Adaptive beamforming has been widely used in wireless communication, microphone array speech processing, radar, sonar, medical imaging, radio astronomy and other areas [1]. The Standard Capon Beamformer (SCB) is an optimal spatial filter that maximizes the array output signal to the interference-plus-noise ratio (SINR), provided that the true covariance matrix and the signal steering vector are accurately known [2]. However, the existence of systematic model mismatch, such as array calibration error, finite snapshots, and others, is inevitable [1,3]. Adaptive beamformers are sensitive to model mismatch, especially when the desired signal is present in the sampling sequences [4]. Therefore, robust adaptive beamforming has been an intensive research topic, and various robust adaptive beamforming techniques have been proposed in the past decades [5].
Traditional RAB techniques include the diagonal loading (DL) technique and the eigenspace-based (ESB) technique [3,6]. The DL-based beamformers are derived by adding a scaled identity matrix on the sample covariance matrix (SCM) [3,7]. Although it has been proved to be effective in mitigating the effects of the finite sample size and SV mismatch [2,8], it is difficult to choose the optimal DL level in different scenarios [3]. Multiple methods were proposed in [2,9,10,11] to compute the DL levels automatically without specifying any additional user parameters [3]. However, these methods suffer from performance degradation at high signal-to-noise ratio (SNR) [12], which is a common problem with this category of method. The ESB beamformers are performed by projecting the nominal SV onto an estimated subspace projector to mitigate the adverse impact caused by noisy disturbance [13,14,15]. The performance of the ESB beamformer is restricted at low SNR [8]. The author in [13] presents a modified ESB method which enhances the SCM for better performance, and a method in [15] adopts the knowledge of the signal angular region as a criterion to construct signal subspace. Although these methods lead to better performance at low SNRs, the problem of performance degradation aroused by large SV mismatches and high input SNRs appears. In [14], the proposed ESB method estimates the desired signal SV by adding an inequality constraint to an optimization problem, which exploits the orthogonality between the precise SV and noise subspace to ensure the validity of the estimation. This method is claimed to be effective in restraining arbitrary SV mismatch. However, there is still a common problem with ESB methods of low SNR performance degradation.
The uncertain-set-based technique utilizes a presumed spherical or ellipsoidal uncertainty set to constrain the signal SV mismatches and correct the nominal SV by solving an optimization problem [8]. The uncertain-set-based beamformers roughly include the following methods [6]: worst-case-based method [1,16], doubly constrained method [9,17], probabilistically constrained method [10,11] and linear programming method [18]. Compared with traditional DL methods, uncertainty-set-based methods are proposed based on clear theoretical analysis and thus have better robustness. Similar to the DL method, there is no reference criterion for the selection of the boundary of the uncertain-set. In fact, the uncertain-set technique transforms the DL level selection problem into an uncertain-set boundary selection problem. Therefore, the uncertainty-set-based method is essentially equivalent to a DL-based method or belongs to its extended class [19], and the uncertainty-set-based method is no longer valid when the SV mismatch is large or the precise SV mismatch does not locate in the presumed set [8,20].
The above methods mainly focus on correcting signal SV and enhancing the covariance matrix [21]. Since the desired signal component in the SCM is not eliminated [22], the above methods always suffer from performance degradation at high input SNRs [6]. To avoid this issue, the authors in [23] pioneer a novel RAB method called INCM reconstruction based method. This method has been proved to effectively eliminate the desired signal components. Gu modifies the above method in [5] and improves the performance, where the INCM is reconstructed by integrating over the complement of the desired signal angular region [6]. Unfortunately, this method performs poorly when the array calibration error exists in the systematic model, and the complexity of this method is very high. Authors in [4,24,25] represent multiple methods to reduce complexity. However, good performance can be obtained only in certain conditions. To resist multiple types of mismatches, authors construct a new constraint in [19] and add it to a QCQP optimization problem to estimate and correct the desired signal SV mismatch. For reducing INCM reconstruction errors, the residual noise is considered and eliminated from the Capon power spectrum estimator to reconstruct INCM and the desired signal SV more precisely in [7]. It is true that the incident signal may be covered by noise in the Capon power spectrum when the SV mismatch is large enough [6]. To overcome this drawback, the methods proposed in [26,27] substitute the Capon power spectrum with the maximum entropy power spectrum. Compared with [19], these methods can effectively reduce the complexity while maintaining good performance. Different from interference angular region integral-based or interference SV estimation-based INCM reconstruction RAB methods [6], the methods in [6,20,21] reconstruct INCM by using a subspace projection operation to filter out the desired signal component from the SCM. These methods achieve good performance with low complexity.
In this study, we propose a new INCM reconstruction-based RAB method [3]. Firstly, we estimate the desired signal SV via employing the MNSP approach. Secondly, we estimate nominal interference SVs via the maximum entropy power spectrum. Subsequently, we establish an oblique projection matrix based on the corrected desired signal SV and the signal subspace of the SCM to correct each interference nominal SV and power. Theoretical analysis and multiple simulations are carried out to illustrate that the proposed RAB method performs well and can effectively combat various mismatches.
This paper is arranged into the following five sections. The necessary signal model and background associated with adaptive beamformer are described in detail in Section 2. In Section 3, the implementation step of this method is recounted in detail, including estimating the desired signal SV by MNSP and correcting the interference SV by oblique projection. In Section 4, multiple simulations and analysis are presented [28]. Finally, conclusions of the study are presented in Section 5. Acronyms used in this paper are las follow: robust adaptive beamforming (RAB), standard Capon beamformer (SCB), signal to interference-plus-noise ratio (SINR), direction of arrival (DOA), signal-to-noise ratio (SNR), interference-to-noise-ratio (INR), steering vector (SV), the minimum norm of subspace projection (MNSP), interference-plus-noise covariance matrix (INCM), sample covariance matrix (SCM), diagonal loading (DL), eigenspace-based (ESB), quadratically constrained quadratic programs (QCQP), uniform linear array (ULA), minimum variance distortionless response (MVDR), independently and identically distributed (i.i.d).

2. Signal Model and Background

Consider a half-wavelength-spaced ULA made up of M-omnidirectional-sensors [28]. L far-field uncorrelated narrowband signals, including one desired signal and L 1 interferences, illuminate the ULA [6,29]. The complex sampling sequences received by the ULA at k time-slot are expressed as [30]:
x k = x s k + x i k + x n k ,
where x s k = s 0 k a 0 denotes an M × 1 vector of the desired signal, x i k = l = 1 L 1 s l k a l denotes interference signal vector, and s l k ,   l = 0 , , L 1 denotes the l th signal vector, and its corresponding SV is a l . x n k represents additive complex Gaussian noise with a zero mean and a variance of σ n 2 [6]. The nominal SV associated with DOA θ can be written as:
a θ = 1 , e j 2 π d λ s i n θ , , e j M 1 2 π d λ s i n θ T = 1 , e j π s i n θ , , e j M 1 π s i n θ T ,
where T represents the transpose operator, λ is the carrier wavelength, and d is the interelement spacing [31]. The output sequence of beamformer can be written as:
y k = w H x k ,
where w = w 1 , , w M T is the complex weight vector of the beamformer [30], H represents the Hermitian transpose [32], and w can be calculated by maximizing the output SINR [6]. When precise σ 0 2 and a 0 is previously known, the definition of output SINR is formulated as follow:
S I N R σ 0 2 w H a 0 2 w H R i + n w ,
where E is the expectation operator of the stochastic variable [6]. The power of the desired signal is expressed as its mathematical expectation; that is; σ 0 2 = E s 0 k 2 , and the corresponding SV is a 0 . R i + n M × M is precise INCM, and its definition is formulated as [31]:
R i + n = E x i k + x n k x i k + x n k H = l = 1 L 1 σ l 2 a l a l H + E x n k x n H k = R i + σ n 2 I ,
where σ l 2 = E s l k 2 is the l th interference power. I M × M , and σ n 2 represents identity matrix and noise power. The maximization of (4) is equivalent to MVDR beamformer or called the Capon beamformer [33,34,35], and it is written as:
min w w H R i + n w     s . t . w H a 0 = 1 .
The solution of the Capon beamformer is written as:
w o p t = R i + n 1 a 0 a 0 H R i + n 1 a 0 .
Unfortunately, R i + n and a 0 cannot be obtained in practical scenarios. Therefore, they are usually replaced by the SCM R ^ = 1 K k = 1 K x k x H k with K sample snapshots and nominal SV a ¯ 0 [7].
It has been proved that replacing R i + n by the SCM R ^ does not change the optimal output SINR [6]. At this point, the optimal weight vector becomes the sample covariance inversion (SMI) beamformer [7]:
w S M I = R ^ 1 a ¯ 0 a ¯ 0 H R ^ 1 a ¯ 0 .
The SCM R ^ converge to the precise one only when K . Due to the mismatches of a ¯ 0 and R ^ , the performance of the SMI beamformer will severely deteriorate, especially in the case of high SNR [33,36].

3. Proposed RAB Method

In this section, the implementation steps of our proposed method are described in detail. For deriving optimal weight vector, reconstructed INCM and the estimated desired signal SV need to be obtained first [37]. The desired signal SV is estimated via employing the MNSP approach [37]. The INCM is reconstructed by two steps: interference DOAs estimation and interference SVs and powers correction.

3.1. MNSP for Desired Signal SV Estimation

In this subsection, we suppose that the DOA of the desired signal θ 0 lies in the angular region Θ s , and all DOAs of interference θ l lie in Θ l , which is the complement region of Θ s [3]. For estimating the desired signal SV, we employ MNSP to correct the nominal desired signal SV. Utilizing the orthogonality of noise subspace with incident signal SV [38], MNSP approach minimizing the norm of the projection of SV on the noise subspace estimates the desired signal SV [20]. Eigen-decompose the matrix R ^ :
R ^ = i = 1 M γ i u i u i H = U Γ U H = U s Γ s U s H + U n Γ n U n H ,
where γ 1 γ 2 γ M 1 γ M are the eigenvalues of R ^ arranged in descending order [20], and u i i = 1 , , M are eigenvectors of R ^ associating with γ i . γ M usually considered as a rough estimation of the noise power σ ˜ n 2 [19]. U s = u 1 , , u L M × L and U n = u L + 1 , , u M M × M L consist of L dominant eigenvectors and M L residual eigenvectors, respectively. Utilizing the eigenvalue-decomposition, we can estimate the column number of U s quickly by min L i = 1 L γ i i = 1 M γ i ζ . Γ s = d i a g γ 1 , , γ L M × L and Γ n = d i a g γ L + 1 , , γ M M × M L are eigen-diagonal matrices. According to the properties of the eigen subspace [6], we have:
s p a n u 1 , , u L = s p a n a 0 , , a L 1 ,
where s p a n u 1 , , u L and s p a n a 0 , , a L 1 are signal subspaces spanned by the L dominant eigenvector group u 1 , , u L and precise signal SV group a 0 , , a L 1 , respectively. The signal subspace is orthogonal to the noise subspace spanned by U n [39]. For any precise SV, a l is a linear combination of eigenvector columns U s and orthogonal to noise subspace. It means that U n H a l = 0 . However, the SV mismatch destroys its orthogonality and U n H a l = 0 becomes U n H a ¯ l 0 . Assume that e is the mismatch vector between nominal SV and the precise one; precise SV is expressed as a 0 = a ¯ 0 + e [8]. Therefore, the MNSP approach corrects the a ¯ l by estimating a mismatch vector e to reestablish the orthogonality between U n H and a ˜ l = a ¯ l + e . The desired signal SV can be estimated by solving the following optimization problem [14,40]:
m i n e a ¯ 0 + e H U n U n H a ¯ 0 + e s . t . a ¯ 0 H e = 0 a ¯ 0 + e H U n 2 a ¯ 0 H U n 2 e 2 ξ ,
where e is the orthogonal component of the mismatch vector e [8]. The parallel component is referred to as e . Since e is actually a scaling version of a ¯ 0 , it can be easily proved that e has no impact on output SINR by substituting (7) into (4). Ignoring e , the corrected desired signal SV is formulated as a ˜ 0 = a ¯ 0 + e . In (11), the equality constraint is introduced to maintain the orthogonality between a ¯ 0 and e [5,41]. The first inequality constraint is introduced to ensure that a ˜ 0 = a ¯ 0 + e converges to the desired signal angular region Θ s . The second inequality constraint is used to limit the e 2 to be a small range which guarantees the validity of desired signal SV estimation at low SNR [41]. This optimization problem (11) is a QCQP problem and can be easily solved using a convex optimization toolbox, such as CVX [5,15,19]. Thanks to the idea in [42], we replace R ^ by R ^ X = 1 2 R ^ + I V R * ^ I V against the coherent scattering mismatch in this paper. I V is a M × M anti-identity matrix. The Hermite’s Toeplitz matrix R ^ X is an unbiased estimation of R ^ [42].
Assume that a signal arrives at a half-wavelength-spaced ULA made up of M = 10 sensors from θ 0 = 3 . The values of U n H a θ 2 and U n H a θ 0 2 are drawn in Figure 1. It can be observed that the existence of model mismatch destroys the orthogonality of U n and a ¯ 0 , and (11) can effectively correct the nominal signal SV a ¯ 0 .

3.2. INCM Reconstruction and Beamformer Weight Vector Calculation

To perform INCM reconstruction, the interference SVs and powers need to be estimated first. Based on the advantages of high resolution and low complexity, the maximum entropy power spectrum is adapted to searching each interference DOA in the angular region Θ l . The maximum entropy power spectrum is written as:
P ^ m e p s θ = 1 ϵ p a ¯ H θ R ^ 1 u 1 2 , θ Θ l ,
where u 1 = 1   0   0 T M × 1 , and ϵ p = 1 u 1 T R ^ 1 u 1 . Due to the impact of noise, peaks located in Θ l may contain some fake peaks that do not correspond to any real interference [19]; σ ˜ n 2 is set as a threshold for eliminating fake peaks. So far, we obtain P peaks corresponding to interference DOAs by searching J points in Θ l , and P nominal interference SVs are a ¯ p ,   p = 1 , 2 , P . It has been proved that the Capon, as well as maximum entropy, power spectrum would cause a peak power under-estimation problem [20,43] in the case of model mismatch. We adopt the oblique projection matrix to correct interference SVs for overcoming this drawback.
Given a matrix A M × N N M , A represents a linear subspace of M spanned by the columns of A [29], P A = A A + represents a projection matrix associated with A , and P A = I P A represents an orthogonal projection matrix associated with A [44]. The oblique projection matrix with range space A and null space B is denoted by [45]:
E A B = A A H P B A 1 A H P B ,
where · + represents a pseudo-inverse operator. E A B M × M is an idempotent and non-symmetric matrix, and E A B has the following resolution and nulling properties [45]:
E A B A = A ;     E A B B = 0 ;     P A B = E A B + E B A .
According to the definition of oblique projection, let A = U S , B = a ˜ 0 and P B = I B B + = I B B H B 1 B H . The oblique projection matrix associated with interferences is E A B = U S U S H P B U S 1 U S H P B . Then, the p th interference SV is corrected as a ˜ p = E A B a ¯ p , and the power corresponding to the p th interference is obtained by σ ˜ p 2 = 1 ϵ p a ˜ p H R ^ 1 u 1 2 [46]. The oblique projection is mainly used to correct the nominal interference SV mitigating the impact of desired signal components on the p th interference SV [20,47]. Therefore, the oblique projection can correct interference SVs more accurately compared with orthogonal projection. Based on each corrected interference SV and power, the more precise INCM is reconstructed as (5) [21]:
R ˜ i + n = i = 1 P σ ˜ i 2 a ˜ i a ˜ i H + σ ˜ n 2 I .
Substituting (15) together with a ˜ 0 calculated by (11) back into the Capon beamformer (7), the proposed method base on INCM reconstruction via oblique projection and desired signal SV estimation via MNSP is written as [6,19]:
w ˜ = R ˜ i + n 1 a ˜ 0 a ˜ 0 H R ˜ i + n 1 a ˜ 0 .
In our proposed method, computation complexity is principally dominated by calculating and eigen-decomposing the SCM R ^ , desired signal SV estimation via MNSP, interference peaks searching and interference SVs correction. Calculating the SCM R ^ costs about O M 2 + M / 2 flops. Matrix inversion and eigen-decomposition cost about O M 3 flops. Estimating the desired signal SV via MNSP requires solving a QCQP problem, and the complexity is about O M 3.5 flops. The interference peaks searching and interference SV correction cost about O m a x 2 M J ,   2 M 3 flops, where J is the searching points in the angular region of Θ l [20]. In general, J < M 2 . Therefore, the overall complexity of this proposed method is roughly O m a x M 3.5 ,   2 M 3 . By reason of the foregoing, the proposed method is summarized as follows (Method 1):
Method 1. The proposed INCM reconstruction-based RAB method.
1: Calculate the SCM R ^ and R ^ X , eigen-decompose R ^ X , and obtain U s , U n , σ ˜ n 2 ;
2: Correct the desired signal SV by solving the QCQP problem (11);
3: Obtain the interference SV estimates a ¯ 1 ,   , a ¯ P using the maximum entropy power spectrum (12);
4: Obtain each corrected interference SV and corresponding power by employing the oblique projection matrix, and reconstruct the INCM via (15);
5: Substitute R ˜ i + n and a ˜ 0 back into (7) to obtain the weight vector.

4. Numerical Simulation and Analysis

In this section, we perform several numerical simulations under the condition of a half-wavelength ULA with M = 10 omnidirectional sensors [6,20]. We assume one desired signal impinging from θ 0 = 3 , two interferences impinging from θ 1 = 36 and θ 2 = 43 . The additive noise is modeled as a complex circularly symmetric Gaussian zero-mean spatially and temporally white process [5]. All signals, including noise, are assumed to be independent of each other. The input INR is fixed at I N R = 20 d B . In the performance comparison of mean output SINR versus the input SNR [5], the snapshots are set as K = 100 , and S N R = 30 d B in the performance comparison of output SINR versus different snapshots. To obtain each output SINR point, 200 Monte Carlo trials are used [25].
The methods involved in the comparison are as follows: subspace projection and covariance matrix reconstruction (SPCMR) beamformer [20], INCM reconstruction via subspace projection (INCM-SP) beamformer [6], INCM reconstruction and steering vector estimation (INCM-SVE) beamformer [19], INCM reconstruction via orthogonality of subspace (INCM-OS) beamformer [21], INCM reconstruction via the intersection of subspaces (INCM-IS) beamformer [33], INCM reconstruction with iterative mismatch approximation (INCM-IMA) beamformer [3] and INCM reconstruction via power method processing and spatial spectrum matching (PMPSSM) beamformer [48]. Among these methods involved in the comparison, all the angular regions are set as Θ s = Θ = θ 0 5 , θ 0 + 5 , Θ l = Θ ¯ = 90 , θ 0 5 θ 0 + 5 , 90 and Θ i = θ i 5 , θ i + 5 . The L = 7 residual eigenvectors of the matrix C are used, and the boundary of RCB is ϵ = 0.1 in INCM-SVE [6]. The N = 7 principal eigenvectors of B are utilized for B 1 in INCM-OS. The constant parameter is μ = 0.9 and τ = 0.1 in SPCMR. The boundaries of amplitude and phase mismatch are set as ϵ = 0.3 and Φ = 6 , the search depth is set as d e p t h = 10 , and the saturation value is φ = 0.05 in INCM-SP and INCM-IMA [3]. The number of iterations is set as k = 4 in PMPSSM. The determinate threshold ρ = ρ l = 0.9 in INCM-IS. In the proposed method, ζ = 0.95 and ξ = 0.1 . In the following simulations, the optimal output SINR is calculated by [19]:
S I N R o p t = σ 0 2 a 0 H R i + n 1 a 0 .

4.1. Simulation Example 1

In this subsection, we consider the detrimental impact on output SINR of the tested methods when amplitude and phase mismatch existed in the signal model. From the above description, the incident signal SV mismatch is modelled as a i = a ¯ i + e ,   i = 0 , , L 1 . To facilitate the analysis, it is indispensable to formally transform the mismatch model of the signal SV. Suppose that α m and β m , respectively, represent the amplitude and phase mismatch on m th array-sensor which are extracted from the distribution N 1 , 0.05 2 and N 0 , 5 2 [19], and the m th element of the precise SV is formulated as a m = α m a ¯ m e j β m [6,19]. The output SINR versus input SNR is delineated in Figure 2a. Our proposed method performs similarly to SPCMR, INCM-SP, INCM-IMA and INCM-SVE at high SNRs. Figure 2b depicts the output SINRs versus snapshot points. The proposed method usually performs better than the other tested methods, and the proposed method can availably resist the finite sampling mismatch [8,36].

4.2. Simulation Example 2

In this subsection, we consider the detrimental impact on output SINR of the tested methods when random look direction mismatch existed in the signal model. Suppose that the random look direction mismatch of incident signal extracts from uniform distribution U 5 , 5 [3]. It is important to note that the amplitude and phase mismatch persist in subsequent simulations. In addition, the random DOAs change in each independent trial while remaining constant over samples [19,31]. The output SINR versus input SNR is delineated in Figure 3a [49]. The performance of our proposed method outperforms other tested methods at high SNRs and is slightly lower than INCM-SVE and INCM-IMA at low SNRs. Figure 3b illustrates that the proposed method outperforms other methods at different snapshots, even with random look direction mismatch.

4.3. Simulation Example 3

In this subsection, we consider the detrimental impact on output SINR of the tested methods when incoherent local scattering mismatch existed in the signal model. Suppose that the desired signal has a time-varying signature [7], which is modelled as:
x ^ s k = s 0 k a 0 + p = 1 4 s p k a ¯ θ p ,
where a 0 and a ¯ θ p represent the desired signal SVs corresponding to the direct-path and incoherent scattering paths [21]. In each trial, the DOAs θ p are independently distributed in a Gaussian distribution drawn from a random generator N θ 0 , 4 [6]; s p k are i.i.d zero-mean complex Gaussian random variables drawn from a random generator N 0 , 1 [31], and s p k changes both from trial to trial and from sample to sample [19]. Therefore, the rank of desired signal covariance matrix R s is larger than one, and the output SINR should be expressed as [6]:
S I N R o p t = w H R s w w H R i + n w .
By maximizing the above SINR in (19), we can obtain the weight vector w o p t through the following equation:
w o p t = P R i + n 1 R s ,
where P denotes the prime eigenvector of a matrix [4,8]. In Figure 4a, we can observe that our proposed method has similar performance with INCM-SP, INCM-SVE and INCM-IS at high SNRs and has the best performance at low SNRs. Figure 4b shows the deviations between the tested methods under the conditions of different snapshots. The proposed method obtains the similar output SINR with INCM-SP, INCM-SVE and INCM-IS. Moreover, the IMCMPM-SSM method suffers severe performance degradation in this case.

4.4. Simulation Example 4

In this subsection, we consider the detrimental impact on output SINR of the tested methods when coherent local scattering mismatch existed in the signal model. The coherent local scattering mismatch is usually encountered in multipath propagation scenarios [20]. Suppose that the desired signal is distorted by local scattering mismatch and consists of four coherent paths [6]. The precise desired signal SV becomes:
a ^ 0 = a 0 + p = 1 4 e j ϕ p a ¯ θ p ,
where a 0 and a ¯ θ p represent desired signal SV and coherent scattering path from θ p , respectively. In each trial, θ p ,   p = 1 , 2 , 3 , 4 are independently distributed in a Gaussian distribution drawn from a random generator N θ 0 , 4 [6]. The path phases ϕ p extract from uniform distribution N 0 , 2 π in each trial.   θ p and ϕ p remain fixed over the samples [3,6]. Moreover, the norm of a ^ 0 is no longer equal to M , and optimal output SINR is increased. It can be observed from Figure 5a that the optimal curve has an approximately 6 d B increment due to the extra coherent scattering paths [3,6]. The performance of our proposed method outperforms other tested methods at high SNRs and is merely lower than INCM-SVE, INCM-IMA and INCMPM-SSM at low SNRs. Figure 5b illustrates that our proposed method generally outperforms others under the condition of finite sampling [8]. Moreover, the IMCMPM-SSM method has good performance compared to the condition of incoherent local scattering mismatch.

5. Conclusions

This paper utilizes the MNSP approach and established oblique projection matrix to estimate desired signal SV and reconstruct INCM respectively. An MNSP approach minimizes the norm of the projection of SV on the noise subspace to estimate the desired signal SV. An oblique projection matrix is used to alleviate the leakage problem of spatial power spectrum estimation and correct interference SVs. The proposed method can treat various types of SV mismatches while merely requiring a small amount of knowledge about the array configuration and angular region in which the desired signal is located. Multiple simulations are carried out and illustrate that this proposed RAB method performs well and can effectively combat various mismatches.

Author Contributions

All authors have significantly contributed to the research presented in this manuscript; Y.D. presented the main idea and wrote the manuscript; Y.G. edited the manuscript; X.Y. and W.C. reviewed and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Guangxi Major Projects of Science and Technology under Grant AA21077006; in part by Guangxi Key Laboratory of Wireless Wideband Communication and Signal Processing under Grants GXKL06200117, GXKL06200128.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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]
  2. Du, L.; Li, J.; Stoica, P. Fully Automatic Computation of Diagonal Loading Levels for Robust Adaptive Beamforming. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 449–458. [Google Scholar] [CrossRef]
  3. Duan, Y.; Zhang, S.; Cao, W. Covariance matrix reconstruction with iterative mismatch approximation for robust adaptive beamforming. J. Electromagn. Waves 2021, 35, 2468–2479. [Google Scholar] [CrossRef]
  4. Gu, Y.J.; Goodman, N.A.; Hong, S.H.; Li, Y. Robust adaptive beamforming based on interference covariance matrix sparse reconstruction. Signal Process. 2014, 96, 375–381. [Google Scholar] [CrossRef]
  5. Gu, Y.J.; 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] [CrossRef]
  6. Duan, Y.L.; Yu, X.H.; Mei, L.R.; Cao, W.P. Low-Complexity Robust Adaptive Beamforming Based on INCM Reconstruction via Subspace Projection. Sensors 2021, 21, 7783. [Google Scholar] [CrossRef] [PubMed]
  7. Zhu, X.Y.; Ye, Z.F.; Xu, X.; Zheng, R. Covariance Matrix Reconstruction via Residual Noise Elimination and Interference Powers Estimation for Robust Adaptive Beamforming. IEEE Access 2019, 7, 53262–53272. [Google Scholar] [CrossRef]
  8. Huang, L.; Zhang, J.; Xu, X.; Ye, Z. Robust Adaptive Beamforming With a Novel Interference-Plus-Noise Covariance Matrix Reconstruction Method. IEEE Trans. Signal Process. 2015, 63, 1643–1650. [Google Scholar] [CrossRef]
  9. Li, J.; Stoica, P.; Wang, Z.S. Doubly constrained robust Capon beamformer. IEEE Trans. Signal Process. 2004, 52, 2407–2423. [Google Scholar] [CrossRef]
  10. Vorobyov, S.A.; Rong, Y.; Gershman, A.B. Robust adaptive beamforming using probability-constrained optimization. In Proceedings of the 2005 IEEE/SP 13th Workshop on Statistical Signal Processing (SSP), Bordeaux, France, 17–20 July 2005; Volumes 1–2, pp. 869–874. [Google Scholar]
  11. Vorobyov, S.A.; Chen, H.H.; Gershman, A.B. On the Relationship Between Robust Minimum Variance Beamformers With Probabilistic and Worst-Case Distortionless Response Constraints. IEEE Trans. Signal Process. 2008, 56, 5719–5724. [Google Scholar] [CrossRef]
  12. Liu, J.B.; Gui, G.; Xie, W.; Cong, X.C.; Wan, Q.; Adachi, F. Robust Widely Linear Beamforming via an IAA Method for the Augmented IPNCM Reconstruction. IEICE Trans. Fundam. Electron. 2017, E100a, 1562–1566. [Google Scholar] [CrossRef]
  13. Huang, F.; Sheng, W.X.; Ma, X.F. Modified projection approach for robust adaptive array beamforming. Signal Process. 2012, 92, 1758–1763. [Google Scholar] [CrossRef]
  14. Jia, W.M.; Jin, W.; Zhou, S.H.; Yao, M.L. Robust adaptive beamforming based on a new steering vector estimation algorithm. Signal Process. 2013, 93, 2539–2542. [Google Scholar] [CrossRef]
  15. Huang, L.; Zhang, B.; Ye, Z.F. Robust Adaptive Beamforming Using a New Projection Approach. In Proceedings of the 2015 IEEE International Conference on Digital Signal Processing (DSP), Singapore, 21–24 July 2015; pp. 1181–1185. [Google Scholar]
  16. Yu, Z.L.; Gu, Z.H.; Zhou, J.J.; Li, Y.Q.; Ser, W.; Er, M.H. A Robust Adaptive Beamformer Based on Worst-Case Semi-Definite Programming. IEEE Trans. Signal Process. 2010, 58, 5914–5919. [Google Scholar] [CrossRef]
  17. Beck, A.; Eldar, Y.C. Doubly constrained robust capon beamformer with ellipsoidal uncertainty sets. IEEE Trans. Signal Process. 2007, 55, 753–758. [Google Scholar] [CrossRef] [Green Version]
  18. Jiang, X.; Zeng, W.J.; Yasotharan, A.; So, H.C.; Kirubarajan, T. Robust Beamforming by Linear Programming. IEEE Trans. Signal Process. 2014, 62, 1834–1849. [Google Scholar] [CrossRef]
  19. Zheng, Z.; Zheng, Y.; Wang, W.Q.; Zhang, H.B. Covariance Matrix Reconstruction With Interference Steering Vector and Power Estimation for Robust Adaptive Beamforming. IEEE Trans. Veh. Technol. 2018, 67, 8495–8503. [Google Scholar] [CrossRef]
  20. Ai, X.Y.; Gan, L. Robust Adaptive Beamforming With Subspace Projection and Covariance Matrix Reconstruction. IEEE Access 2019, 7, 102149–102159. [Google Scholar] [CrossRef]
  21. Zhu, X.; Xu, X.; Ye, Z. Robust adaptive beamforming via subspace for interference covariance matrix reconstruction. Signal Process. 2020, 167, 107289. [Google Scholar] [CrossRef]
  22. Chen, P.; Yang, Y.X. An uncertainty-set-shrinkage-based covariance matrix reconstruction algorithm for robust adaptive beamforming. Multidimens. Syst. Signal Process. 2021, 32, 263–279. [Google Scholar] [CrossRef]
  23. Mallipeddi, R.; Lie, J.P.; Suganthan, P.N.; Razul, S.G.; See, C.M.S. Robust adaptive beamforming based on covariance matrix reconstruction for look direction mismatch. PIER Lett. 2011, 25, 34–46. [Google Scholar] [CrossRef] [Green Version]
  24. Ruan, H.; de Lamare, R.C. Robust Adaptive Beamforming Using a Low-Complexity Shrinkage-Based Mismatch Estimation Algorithm. IEEE Signal Process. Lett. 2014, 21, 60–64. [Google Scholar] [CrossRef] [Green Version]
  25. Shen, F.; Chen, F.F.; Song, J.Y. Robust Adaptive Beamforming Based on Steering Vector Estimation and Covariance Matrix Reconstruction. IEEE Commun. Lett. 2015, 19, 1636–1639. [Google Scholar] [CrossRef]
  26. Mohammadzadeh, S.; Nascimento, V.H.; Lamare, R.C.d.; Kukrer, O. Maximum Entropy-Based Interference-Plus-Noise Covariance Matrix Reconstruction for Robust Adaptive Beamforming. IEEE Signal Process. Lett. 2020, 27, 845–849. [Google Scholar] [CrossRef]
  27. Mohammadzadeh, S.; Nascimento, V.H.; Lamare, R.C.d.; Kukrer, O. Low-Cost Adaptive Maximum Entropy Covariance Matrix Reconstruction for Robust Beamforming. In Proceedings of the 2020 54th Asilomar Conference on Signals, Systems, and Computers, California, CA, USA, 1–4 November 2020; pp. 1462–1466. [Google Scholar]
  28. Zhang, Y.P.; Li, Y.J.; Gao, M.G. Robust adaptive beamforming based on the effectiveness of reconstruction. Signal Process. 2016, 120, 572–579. [Google Scholar] [CrossRef] [Green Version]
  29. Zhuang, J.; Manikas, A. Interference cancellation beamforming robust to pointing errors. IET Signal Process. 2013, 7, 120–127. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, X.-J.; Feng, D.-Z.; Nie, W.-K.; Lin, S.-Y.; Hu, H.-S. Robust adaptive beamforming algorithm based on damped singular value decomposition regularization. Digit. Signal Process. 2022, 122, 103356. [Google Scholar] [CrossRef]
  31. Zheng, Z.; Yang, T.; Wang, W.Q.; So, H.C. Robust Adaptive Beamforming via Simplified Interference Power Estimation. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 3139–3152. [Google Scholar] [CrossRef]
  32. Shao, X.; Hu, T.; Zhang, J.; Pan, Y.; Zhang, Z.; Wei, Y. A low-sidelobe adaptive beamforming algorithm with single snapshot. In Proceedings of the IET International Radar Conference (IET IRC 2020), Chongqing, China, 4–6 November 2020; pp. 682–686. [Google Scholar]
  33. Wang, R.; Wang, Y.; Han, C.; Gong, Y.; Wang, L. Robust Adaptive Beamforming Based on Interference Covariance Matrix Reconstruction and Steering Vector Estimation. In Proceedings of the 2021 IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC), Xi’an, China, 17–19 August 2021; pp. 1–5. [Google Scholar]
  34. Liao, B.; Chan, S.C.; Tsui, K.M. Recursive Steering Vector Estimation and Adaptive Beamforming under Uncertainties. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 489–501. [Google Scholar] [CrossRef] [Green Version]
  35. Gu, Y.J.; Shi, Z.G.; Chen, K.S.; Li, Y. Robust adaptive beamforming for steering vector uncertainties based on equivalent DOAs method. Prog. Electromagn. Res. 2008, 79, 277–290. [Google Scholar] [CrossRef]
  36. Huang, J.S.; Su, H.T.; Yang, Y. Low-complexity robust adaptive beamforming method for MIMO radar based on covariance matrix estimation and steering vector mismatch correction. IET Radar Sonar Navig. 2019, 13, 712–720. [Google Scholar] [CrossRef]
  37. Xie, Z.; Zhu, J.H.; Fan, C.Y.; Huang, X.T. An Improved Sub-Array Adaptive Beamforming Technique Based on Multiple Sources of Errors. J. Mar. Sci. Eng. 2020, 8, 757. [Google Scholar] [CrossRef]
  38. Lee, J.H.; Cheng, K.P.; Wang, C.C. Robust adaptive array beamforming under steering angle mismatch. Signal Process. 2006, 86, 296–309. [Google Scholar] [CrossRef]
  39. Xu, X.; Ye, Z.F.; Zhang, Y.F. DOA Estimation for Mixed Signals in the Presence of Mutual Coupling. IEEE Trans. Signal Process. 2009, 57, 3523–3532. [Google Scholar] [CrossRef]
  40. Yang, B.; Li, W.X.; Li, Y.Y.; Zhang, Q. Robust adaptive null broadening beamforming based on subspace projection. Int. J. Electron. 2022, 1–15. [Google Scholar] [CrossRef]
  41. Gong, C.; Huang, L.; Xu, D.; Ye, Z. Knowledge-aided robust adaptive beamforming with small snapshots. Electron. Lett. 2013, 49, 1258–1259. [Google Scholar] [CrossRef]
  42. Du, Y.; Cui, W.; Mei, F.; Xu, H.; Wang, Y. Coherent signals adaptive beamforming algorithm based on eigenvalue. In Proceedings of the 2021 2nd International Conference on Electronics, Communications and Information Technology (CECIT), Sanya, China, 27–29 December 2021; pp. 388–393. [Google Scholar]
  43. Stoica, P.; Li, J.; Tan, X. On spatial power spectrum and signal estimation using the Pisarenko framework. IEEE Trans. Signal Process. 2008, 56, 5109–5119. [Google Scholar] [CrossRef]
  44. Chen, P.; Zuo, L.; Wang, W. Adaptive Beamforming for Passive Synthetic Aperture with Uncertain Curvilinear Trajectory. Remote Sens. 2021, 13, 2562. [Google Scholar] [CrossRef]
  45. McCloud, M.L.; Scharf, L.L. A new subspace identification algorithm for high-resolution DOA estimation. IEEE Trans. Antennas Propag. 2002, 50, 1382–1390. [Google Scholar] [CrossRef]
  46. Yang, B.; Li, W.X.; Li, Y.Y.; Yue, C.Y. Novel robust adaptive beamforming against unknown mutual coupling. J. Electromagn. Waves 2021, 35, 2447–2467. [Google Scholar] [CrossRef]
  47. Yang, X.P.; Li, Y.Q.; Liu, F.F.; Lan, T.; Long, T.; Sarkar, T.K. Robust Adaptive Beamforming Based on Covariance Matrix Reconstruction With Annular Uncertainty Set and Vector Space Projection. IEEE Antennas Wirel. Propag. Lett. 2021, 20, 130–134. [Google Scholar] [CrossRef]
  48. Mohammadzadeh, S.; Nascimento, V.H.; Lamare, R.C.D.; Kukrer, O. Robust Adaptive Beamforming Based on Power Method Processing and Spatial Spectrum Matching. In Proceedings of the ICASSP 2022—2022 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Singapore, 23–27 May 2022; pp. 4903–4907. [Google Scholar]
  49. Lu, J.; Yang, J.; Yao, Z.C.; Zhu, F.C.; Yu, Z.Y.; Liu, G.B. Robust wideband adaptive beamforming based on covariance matrix reconstruction in the spatial-frequency domain. IET Microw. Antennas Propag. 2021, 15, 309–322. [Google Scholar] [CrossRef]
Figure 1. Comparison of (a) U n H a θ 2 versus θ with fixed incident signal S N R = 5 d B ; (b) U n H a θ 0 2 versus incident signal SNR.
Figure 1. Comparison of (a) U n H a θ 2 versus θ with fixed incident signal S N R = 5 d B ; (b) U n H a θ 0 2 versus incident signal SNR.
Electronics 11 03478 g001
Figure 2. Example 1: output SINR vs.: (a) input SNR; (b) number of snapshots.
Figure 2. Example 1: output SINR vs.: (a) input SNR; (b) number of snapshots.
Electronics 11 03478 g002
Figure 3. Example 2: output SINR vs.: (a) input SNR; (b) number of snapshots.
Figure 3. Example 2: output SINR vs.: (a) input SNR; (b) number of snapshots.
Electronics 11 03478 g003
Figure 4. Example 3: output SINR vs.: (a) input SNR; (b) number of snapshots.
Figure 4. Example 3: output SINR vs.: (a) input SNR; (b) number of snapshots.
Electronics 11 03478 g004
Figure 5. Example 4: output SINR vs.: (a) input SNR; (b) number of snapshots.
Figure 5. Example 4: output SINR vs.: (a) input SNR; (b) number of snapshots.
Electronics 11 03478 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Duan, Y.; Gong, Y.; Yang, X.; Cao, W. Oblique Projection-Based Covariance Matrix Reconstruction and Steering Vector Estimation for Robust Adaptive Beamforming. Electronics 2022, 11, 3478. https://doi.org/10.3390/electronics11213478

AMA Style

Duan Y, Gong Y, Yang X, Cao W. Oblique Projection-Based Covariance Matrix Reconstruction and Steering Vector Estimation for Robust Adaptive Beamforming. Electronics. 2022; 11(21):3478. https://doi.org/10.3390/electronics11213478

Chicago/Turabian Style

Duan, Yanliang, Yanping Gong, Xiaohui Yang, and Weiping Cao. 2022. "Oblique Projection-Based Covariance Matrix Reconstruction and Steering Vector Estimation for Robust Adaptive Beamforming" Electronics 11, no. 21: 3478. https://doi.org/10.3390/electronics11213478

APA Style

Duan, Y., Gong, Y., Yang, X., & Cao, W. (2022). Oblique Projection-Based Covariance Matrix Reconstruction and Steering Vector Estimation for Robust Adaptive Beamforming. Electronics, 11(21), 3478. https://doi.org/10.3390/electronics11213478

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