Low-Complexity Robust Adaptive Beamforming Based on INCM Reconstruction via Subspace Projection

Adaptive beamforming is sensitive to steering vector (SV) and covariance matrix mismatches, especially when the signal of interest (SOI) component exists in the training sequence. In this paper, we present a low-complexity robust adaptive beamforming (RAB) method based on an interference–noise covariance matrix (INCM) reconstruction and SOI SV estimation. First, the proposed method employs the minimum mean square error criterion to construct the blocking matrix. Then, the projection matrix is obtained by projecting the blocking matrix onto the signal subspace of the sample covariance matrix (SCM). The INCM is reconstructed by replacing part of the eigenvector columns of the SCM with the corresponding eigenvectors of the projection matrix. On the other hand, the SOI SV is estimated via the iterative mismatch approximation method. The proposed method only needs to know the priori-knowledge of the array geometry and angular region where the SOI is located. The simulation results showed that the proposed method can deal with multiple types of mismatches, while taking into account both low complexity and high robustness.


Introduction
Adaptive beamforming adjusts the weight vector according to the application environment to enhance the SOI by suppressing interference and noise, and it has been widely used in radar, sonar, microphone array speech processing, wireless communication, radio astronomy, and other areas [1][2][3]. Generally, the standard Capon beamformer (SCB) obtains the maximum-array-output-signal-to-interference-plus-noise ratio (SINR) if the covariance matrix and SOI SV are accurately known [4]. However, severe performance degradation may occur in the presence of SV and INCM mismatches due to the fact of array calibration errors, finite snapshots, and other factors, especially when the SOI component is presented in the INCM [5,6]. Therefore, various RAB algorithms have been proposed to ensure the robustness of beamformers over the past years. In general, these RAB methods can be classified into the following types [6,7]: diagonal loading (DL) technique, eigenspace-based (ESB) technique, uncertain-set based technique, and covariance matrix reconstruction-based technique.
DL is one of the most classical RAB methods for improving the robustness of a beamformer, which is derived by imposing a quadratic constraint either on the norm of the weight vector or on its SOI SV [8,9]. However, its major challenge is that it is difficult to choose the optimal DL level in different scenarios. To overcome this drawback, parameterfree methods in [4,[10][11][12][13] can automatically compute the DL level without specifying In this paper, a low-complexity RAB method based on INCM reconstruction and SOI SV estimation is proposed. Unlike previous methods, the INCM in the proposed method was reconstructed utilizing the orthogonality of subspace. First, based on the idea of the matrix filter in [36,37], the minimum mean square error criterion was employed to construct the blocking matrix. Then, we performed eigen-decomposition on the SCM and obtained the orthogonal projection matrix of the signal subspace. By projecting the blocking matrix onto the orthogonal projection matrix, a projection matrix was obtained. Subsequently, the INCM was reconstructed by replacing the eigenvector columns of the SCM, which can span to the signal subspace with the corresponding eigenvectors of the projection matrix. Finally, the SOI SV was estimated by employing the iterative mismatch approximation method presented in [29]. The theoretical analysis and simulation results demonstrated that the proposed method can efficiently deal with multiple types of mismatches.
The rest of this paper is organized as follows. The signal model and necessary background regarding the adaptive beamforming method is introduced in Section 2. In Section 3, the proposed RAB methods are described in detail, and the feasibility analysis of the blocking matrix is performed. The simulation results are provided in Section 4. Finally, conclusions are drawn in Section 5.

Signal Model and Background
Consider a uniform linear array (ULA) composed of M omnidirectional sensors that are illuminated by L + 1 far-field uncorrelated narrowband signals, which consist of one SOI and L interferences. The array complex sample vector at time k can be presented as: where x s (k) = s 0 (k)a 0 and x i (k) = ∑ L l=1 s l (k)a l , respectively, represent the M × 1 vector of the SOI and interference signal component in the received data. s l (k) and a l (l = 0, . . . , L) are the lth incident signal waveform and corresponding SV. x n (k) is additive complex Gaussian noise with a zero mean and a variance of σ 2 n , which is uncorrelated with all the other signals. In this paper, the sensors were spaced at half of the wavelength. The nominal SV from θ can be written as: where (·) T denotes the transpose, and λ and d, respectively, denote signal wavelength and distance between two adjacent sensors. The output of the beamformer is written as: where w = (w 1 , . . . , w M ) T is the complex beamformer weight vector, and (·) H is the Hermitian transpose. The optimal beamformer weight vector, w, can be calculated by maximizing the output SINR, which is defined as follow: where σ 2 0 = E |s 0 (k)| 2 and a 0 denote the SOI power and SV, E[·] denotes the expectation operator of the stochastic variable. R i+n ∈ C M×M denotes the precise INCM which can be written as: where σ 2 l = E |s l (k)| 2 and I denote the lth interference power and identity matrix, respectively, and σ 2 n is the noise power. The main purpose of the optimal beamformer is to maximize the output SINR while keeping the SOI undistorted at the same time, which is the so-called the minimum variance distortionless response (MVDR) problem [2]: The optimal beamformer weight vector is given by: In practical applications, the precise INCM R i+n is always unavailable, and it is usually replaced by SCMR:R where K denotes the number of snapshots. As K increases,R converge to the actual one. It has been proved that replacing R i+n by SCMR does not change the optimal output SINR [2]. Substituting the actual SOI SV a 0 by the nominal SV a 0 based on the known array structure, the optimal weight vector becomes the sample covariance inversion (SMI) beamformer: With the optimal weight vector, the Capon spatial power spectrum, is employed as a power estimator over all directions [28]: where a(θ) is the nominal SV associated with θ ∈ (−90 • , 90 • ).

Proposed Method
The main idea of the proposed method is to utilize the reconstructed INCM and the corrected SOI SV to derive the optimal weight vector. Depending on the minimum mean square error criterion, a blocking matrix is obtained. The orthogonality of the subspace is employed to derive the projection matrix, and the INCM is reconstructed by replacing the eigenvector columns of the SCM such that they can span to the signal subspace with the corresponding eigenvectors of the projection matrix. The optimal weight vector is obtained along with the SOI SV estimated by the iterative mismatch approximation method.

INCM Reconstruction
Different from Capon power spectrum integration and interference estimation based INCM reconstruction methods, we present a blocking matrix based on the matrix filter principle in [36,37], which is suitable for suppressing signals illuminating within a spe-Sensors 2021, 21, 7783 5 of 13 cific angular region. Consider a blocking matrix G ∈ C M×M , the property of G can be described as: where Θ s denotes the stopband angular region, and the direction of arrival (DOA) of the SOI lies in it. Θ p denotes the passband angular region which contains the locations of interference. Sampling Θ s and Θ p uniformly with the N s and N p sampling points, the corresponding angular sequences can be presented as θ p p = 1, . . . , N p and θ s (s = 1, . . . , N s ). The blocking matrix design problem based on the minimum mean square error criterion can be described as: min where · F denotes Frobenius norm. P = a Θ p , a(Θ s ) ∈ C M×(N p +N s ) denotes the nominal manifold matrix, andP = a Θ p , 0 M×N s ∈ C M×(N p +N s ) denotes the desired manifold matrix. The solution to (12) can be found by taking the gradient of F(G) = G H P −P F and making it equal to zero: Performing eigen decomposition onR yields: where Γ 1 ≥ Γ 2 ≥ · · · Γ M−1 ≥ Γ M denotes the eigenvalues ofR arranged in descending order. The minimum eigenvalue can be approximately considered as the estimation of the noise power σ 2 n [6]. u i ∈ C M×1 is an eigenvector associated with Γ i . U s = (u 1 , . . . , u L+1 ) ∈ C M×(L+1) and U n = (u L+2 , . . . , u M ) ∈ C M×(M−L−1) , respectively, denote the signal subspace eigenvectors and the noise subspace eigenvectors. Γ s = diag(Γ 1 , . . . , Γ L+1 ) ∈ C M×(L+1) and Γ n = diag(Γ L+2 , . . . , Γ M ) ∈ C M×(M−L−1) are diagonal matrices. According to the properties of the eigen subspace, we have: where span{u 1 , . . . , u L+1 } denotes the spanned subspace generated by the vector group {u 1 , . . . , u L+1 }. Then, any accurate SV a l can be expressed as a linear combination of columns of U s [27], which means that the projection matrix Φ = U s U H s G is orthogonal to a 0 . Q = Φ H a(θ) 2 is used to measure the orthogonality between projection matrix Φ and a(θ). Assume that the SOI and interference impinge on the half-wavelength spacing ULA with M = 10 from θ 0 = 3 • , θ 1 = −35 • , and θ 2 = 42 • , the stopband and passband region are set as . It can be observed from Figure 1 that Q will be much smaller when θ = θ 0 than θ = θ 1,2 . This means Φ and a(θ 0 ) are orthogonal or approximately orthogonal. In addition, the value of Q corresponding to θ 1,2 is equal to a(θ 1,2 ) 2 . Hence, projection matrix Φ collects the spatial information of interference and removes the spatial information of the SOI.
Consider that the interference SV lies in the signal subspace spanned by the dominant eigenvectors of Φ. Then, employing eigen decomposition on Φ to obtain the signal subspace: Consider that the interference SV lies in the signal subspace spanned by the dominant eigenvectors of . Then, employing eigen decomposition on to obtain the signal subspace: where = [ 1 , … , ] = [ , ] and = { 1 , … , ), respectively, denote unitary and diagonal matrices. , = 1, … , denotes the eigenvector corresponding to , which are arranged in descending order. In addition, contains + 1 eigenvectors columns corresponding to the + 1 largest eigenvalues, and { 1 , … , +1 } can be considered as the signal subspace.
Substituting back into (14) and replacing , then we can obtain the reconstructed INCM: plotting the power spectrum of ̃+ by replacing ̃+ with ̂ in (10), which is written as: The spatial power spectrum distribution based on (10) and (18) is drawn in Figure 2. The SOI was assumed to be impinging from 0 = 3 ∘ with a fixed = 30 dB and that two interferences were impinging from 1 = −35 ∘ and 2 = 42 ∘ with a fixed interference-to-noise rate (INR) = 20 dB, respectively. This shows that the blocking matrix can effectively filter the SOI components and interference components are retained. Therefore, ̃+ can be used as an INCM to derive the beamformer. Substituting B s back into (14) and replacing U s , then we can obtain the reconstructed INCM: plotting the power spectrum of R i+n by replacing R i+n withR in (10), which is written as: The spatial power spectrum distribution based on (10) and (18) is drawn in Figure 2. The SOI was assumed to be impinging from θ 0 = 3 • with a fixed SNR = 30 dB and that two interferences were impinging from θ 1 = −35 • and θ 2 = 42 • with a fixed interferenceto-noise rate (INR) I NR = 20 dB, respectively. This shows that the blocking matrix can effectively filter the SOI components and interference components are retained. Therefore, R i+n can be used as an INCM to derive the beamformer.

SOI SV Estimation and Beamformer Weight Vector Calculation
For estimating the SOI SV, we employed the iterative mismatch approximation method proposed in [29] to correct the presumed SOI SV. The iterative mismatch approx-

SOI SV Estimation and Beamformer Weight Vector Calculation
For estimating the SOI SV, we employed the iterative mismatch approximation method proposed in [29] to correct the presumed SOI SV. The iterative mismatch approximation method depends on searching for the SV mismatch in the margin of the amplitude and phase error. In the presence of SV mismatches, the actual SOI SV can be written as a 0 = a 0 + e = α • a 0 • e jβ . Employing the principle of estimating the signal steering vectors via maximizing the beamformer output power [16], the estimated SOI SV a 0 can be obtained by solving the pseudo-optimization problem: where • denotes the Hadamard product, and the power spectrum associated with θ q , α, β is and Φ denote the predefined boundary values of the amplitude and phase mismatch, respectively. α = (α 1 , . . . , α M ) T and β = (β 1 . . . , β M ) T denote the amplitude and phase mismatch vectors, respectively. The iterative mismatch approximation method is described in Algorithm 1. In addition, the initial nominal SOI SV is associated with the middle value of Θ s . Algorithm 1. Iterative mismatch approximation method.

9:
Zoom out of the amplitude/phase error area built byΘ α andΘ β to E ap (idx) 10: UpdateΘ α andΘ β 11: end 12: Then, the estimated SOI SV a 0 can be corrected as a q . Substituting (17) together with a 0 back into the Capon beamformer (7), the robust adaptive beamforming based on INCM reconstruction via the projecting matrix and SV estimation can be written as: In a general case, the estimated SOI SV may be imprecise. Hence, we took a 0 as an input parameter and performed multiple iterations to improve accuracy. The iteration was terminated when the following conditions were satisfied: where ϕ and σ 2 q denote a predefined saturation coefficient and power corresponding to a q , respectively. The method we propose is summarized in Algorithm 2.

Algorithm 2. Proposed RAB method.
1: Calculate the SCMR using (8) and eigen decomposeR to obtain U S and σ 2 n ; 2: Obtain the blocking matrix G using (13)  Assuming that depth = D, the computational complexity of estimating the SOI SV would be approximately O M 3 + 3DM flops. In practice, calculating G is independent ofR and can be seen as a pre-processing operation of the proposed method. With K M ∼ = D, the overall complexity of proposed method is approximately O M 2 K flops. In [6], the major computations were conducted to solve a QCQP problem and estimate interference SVs.
Suppose that S denotes the number of search points in -Θ, the computational complexity is O max M 2 S, M 3.5 . In [30], the sample points, N, for the line search significantly affected the adjustment of the SVs. Hence, when N > M, the computational complexity is O LN M 2 . The computational complexity of the MEPS-IPNC in [31] was O M 2 L , where L is a small multiple of M. In [32], the computational complexity was O max M 2 J, M 2 S , where S and J are the number of the sampling points in Θ s and Θ s ∪ Θ i . In [33], the computational complexity was declared to be O M 3 . In [34], the computational complexity was O max M 2 J, M 2 S , where S and J are the number of the sampling points in Θ and -Θ. In [35], the computational complexity was declared to be O M 2 I(L + 1) , where I is the number of sample points for Θ and Θ i−l . Clearly, since the power spectrum calculation, power spectrum integration, and QCQP problem solving are avoided in our proposed method, it has the lowest computational complexity. Furthermore, prior information regarding the SOI angular region and array geometry are needed.

Simulation
In this section, a half-wavelength spacing ULA with M = 10 was considered. We assumed that there existed an SOI impinging from the direction of θ 0 = 3 • and two interferences impinging from θ 1 = −35 • and θ 2 = 42 • . The additive noise was presumed to be a complex circularly symmetric Gaussian zero-mean unit-variance spatially and temporally white process. All these sources were narrowband and assumed to be independent to the noise. To obtain each output SINR point, 200 Monte Carlo trials were used in each simulation. The proposed method was compared with the RAB method based on INCM reconstruction and steering vector estimation (INCM-SVE) in [6], subspace-decomposition and SV adjustment (SDA) in [30], MEPS-IPNC in [31], INCM reconstruction via orthogonality of subspace (INCM-OS) in [32], desired signal eigenvalue replacement (DSEB) in [33], SPCMR in [34] and INCM reconstruction via the intersection of subspaces (INCM-IS) in [35]. For all methods involved in the comparisons, the angular region was presumed to be Θ = Θ s = (θ 0 − 6 • , θ 0 + 6 • ) and -Θ = Θ p = (−90 • , θ 0 − 6 • ) (θ 0 + 6 • , 90 • ). and the interference angular region to be Θ i = (θ i − 6 • , θ i + 6 • ). The number of non-dominant eigenvectors of the matrix C was set as L = 7, and the RCB boundary was = √ 0.1 in [6]. The N = 7 dominant eigenvectors of matrix B were employed for B 1 in [32]. ξ = 0.95 as in [33]. The constant satisfying µ = 0.9 and τ = √ 0.1 were as in [34]. Sampling points L = 5M and S = 20 were as in [31]. The scale factor µ = 0.1 and the sampling points N = 2 µ/0.01 were as in [30]. For the proposed method, the amplitude and phase mismatch boundary were set as = 0.3 and Φ = 6 • , respectively, the depth of the iteration was set as depth = 10, and the saturation value was ϕ = 0.05. Furthermore, the MATLAB CVX toolbox was used to solve the QCQP optimization problem in [6]. In our simulations, the optimal output SINR can be calculated by:

Example 1: Mismatch Due to the Amplitude and Phase Error of the SV
In the first example, the influence of the SVs with arbitrary amplitude and phase errors on the beamformer output SINR was considered. The relationship between the mth element of the nominal SV and the actual SV was modeled as a m = α m a m e jβ m , where the arbitrary amplitude error, α m , and phase error, β m , on each array sensor, respectively, followed the Gaussian distribution N 1, 0.05 2 and N 0, (5 • ) 2 [6]. Figure 3a depicts the output SINR of the tested methods versus the input SNR for the fixed number of snapshots K = 100. It was observed that the proposed method had a similar performance among the tested methods except in [32,33] at high SNRs. In addition, the performance of the proposed method was only lower than in [6] when the SNR was low. However, the computational complexity of our method was obviously lower than that in [6]. In Figure 3b, the output SINRs are shown versus the number of snapshots for the fixed SNR = 30 dB and I NR = 20 dB. The proposed method had a similar performance to the tested methods in [6,30,34,35], and the number of snapshots did not affect the performance of our proposed method.

Example 1: Mismatch Due to the Amplitude and Phase Error of the SV
In the first example, the influence of the SVs with arbitrary amplitude and phase errors on the beamformer output SINR was considered. The relationship between the th element of the nominal SV and the actual SV was modeled as = ̅ , where the arbitrary amplitude error, , and phase error, , on each array sensor, respectively, followed the Gaussian distribution (1, 0.05 2 ) and (0, (5 ∘ ) 2 ) [6]. Figure 3a depicts the output SINR of the tested methods versus the input SNR for the fixed number of snapshots = 100. It was observed that the proposed method had a similar performance among the tested methods except in [32,33] at high SNRs. In addition, the performance of the proposed method was only lower than in [6] when the SNR was low. However, the computational complexity of our method was obviously lower than that in [6]. In Figure  3b, the output SINRs are shown versus the number of snapshots for the fixed = 30 dB and = 20 dB. The proposed method had a similar performance to the tested methods in [6,30,34,35], and the number of snapshots did not affect the performance of our proposed method.

Example 2: Mismatch Due to the Random Look Direction Error
In the second example, the influence of the random look direction errors on the beamformer output SINR was considered. Assuming that the look direction mismatch of both the SOI and interferences were uniformly distributed in (−5 • , 5 • ). That is means that the DOA of the SOI was uniformly distributed in (−2 • , 8 • ), and the DOAs of the two interference were uniformly distributed in (−40 • , 30 • ) and (37 • , 47 • ). Note that the random DOAs of the SOI and interferences changed in each trial while remaining constant over snapshots. Figure 4a shows the output SINRs of the tested methods versus the input SNRs with the fixed snapshots K = 100. It was observed that our proposed method was only inferior to that in [6] in the performance at low SNR and inferior to that in [6,35] at high SNRs. Figure 4b depicts the output SINRs of the tested methods against the snapshot number at SNR = 30 dB and I NR = 20 dB. It was observed that the performance of our proposed method was similar to that in [30] when K > 40. In addition, the methods in [33,34] were significantly affected by mismatches due to the look direction error.
snapshots. Figure 4a shows the output SINRs of the tested methods versus the input SNRs with the fixed snapshots = 100. It was observed that our proposed method was only inferior to that in [6] in the performance at low SNR and inferior to that in [6,35] at high SNRs. Figure 4b depicts the output SINRs of the tested methods against the snapshot number at = 30 dB and = 20 dB. It was observed that the performance of our proposed method was similar to that in [30] when > 40. In addition, the methods in [33,34] were significantly affected by mismatches due to the look direction error.

Example 3: Mismatch Due to the Incoherent Local Scattering Error
In the third example, the influence of the incoherent local scattering error on the beamformer output SINRs was considered. The SOI was assumed to have a time-varying signature, which was modeled as: where 0 denotes the SOI SV. ̅( )( = 1,2,3,4) denotes the incoherent scattering signal SV, and the DOAs, , are independently distributed in a Gaussian distribution drawn from a random generator ( 0 , 4 ∘ ) in each trial. ( ) are independently and identically distributed zero-mean complex Gaussian random variables independently drawn from a random generator, (0,1). In addition, changes from trial to trial, while it remains fixed over the samples. At the same time, ( ) changes both from trial to trial and from sample to sample. In this case, the SOI covariance matrix is no longer a rank-one matrix and the output SINR should be expressed as [5]: The optimal weight vector can be obtained by maximizing the SINR [5]: where {⋅} denotes the principal eigenvector of a matrix. Figure 5a shows the output SINRs of the tested methods versus the input SNRs with the fixed snapshots = 100. It was observed that the performance of our proposed method was similar to that in [6,30] at high SNRs and only lower than in [33] at low SNRs. However, the method in [33] had

Example 3: Mismatch Due to the Incoherent Local Scattering Error
In the third example, the influence of the incoherent local scattering error on the beamformer output SINRs was considered. The SOI was assumed to have a time-varying signature, which was modeled as: where a 0 denotes the SOI SV. a θ p (p = 1, 2, 3, 4) denotes the incoherent scattering signal SV, and the DOAs, θ p , are independently distributed in a Gaussian distribution drawn from a random generator N(θ 0 , 4 • ) in each trial. s p (k) are independently and identically distributed zero-mean complex Gaussian random variables independently drawn from a random generator, N(0, 1). In addition, θ p changes from trial to trial, while it remains fixed over the samples. At the same time, s p (k) changes both from trial to trial and from sample to sample. In this case, the SOI covariance matrix is no longer a rank-one matrix and the output SINR should be expressed as [5]: The optimal weight vector can be obtained by maximizing the SINR [5]: where P{·} denotes the principal eigenvector of a matrix. Figure 5a shows the output SINRs of the tested methods versus the input SNRs with the fixed snapshots K = 100. It was observed that the performance of our proposed method was similar to that in [6,30] at high SNRs and only lower than in [33] at low SNRs. However, the method in [33] had severe performance degradation at high SNRs. Figure 5b depicts the output SINRs of the tested methods against the snapshot number at SNR = 30 dB and I NR = 20 dB. It was observed that the proposed method had a similar performance with the tested methods in [6,30,35], and the number of snapshots did not affect the performance of our proposed method.
severe performance degradation at high SNRs. Figure 5b depicts the output SINRs of the tested methods against the snapshot number at = 30 dB and = 20 dB. It was observed that the proposed method had a similar performance with the tested methods in [6,30,35], and the number of snapshots did not affect the performance of our proposed method.

Example 4: Mismatch Due to the Coherent Local Scattering Error
In the fourth example, the influence of the coherent local scattering mismatch on the beamformer output SINRs was considered. The coherent local scattering mismatch usually occurs in multipath propagation scenarios. Assume that the SOI is distorted by local scattering and consists of four coherent paths; the actual SV is formed as: where 0 denotes the SOI SV. ̅( )( = 1,2,3,4) denotes the coherent signal path from . are independently distributed in a Gaussian distribution drawn from a random generator, ( 0 , 4 ∘ ), in each trial.
denotes the path phase and uniformly distributed in (0,2 ) from trial to trial. and change from trial to trial, while it remains fixed over the samples. Figure 6a shows the output SINRs of the tested methods versus the input SNRs with the fixed snapshots = 100. It was observed that the performance of the optimal beamformer had an approximately 6 dB increment in output SINRs due to the extra paths. The performance of our proposed method was similar to that in [6,30,35] at high SNRs. The method in [6] achieved the best performance at the cost of the highest complexity compared with the others. Figure 6b depicts the output SINRs of the tested methods against the snapshot number at = 30 dB and = 20 dB. It was observed that the proposed method had a small impact on the number of snapshots.

Example 4: Mismatch Due to the Coherent Local Scattering Error
In the fourth example, the influence of the coherent local scattering mismatch on the beamformer output SINRs was considered. The coherent local scattering mismatch usually occurs in multipath propagation scenarios. Assume that the SOI is distorted by local scattering and consists of four coherent paths; the actual SV is formed as: where a 0 denotes the SOI SV. a θ p (p = 1, 2, 3, 4) denotes the coherent signal path from θ p . θ p are independently distributed in a Gaussian distribution drawn from a random generator, N(θ 0 , 4 • ), in each trial. Φ p denotes the path phase and uniformly distributed in (0, 2π) from trial to trial. θ p and Φ p change from trial to trial, while it remains fixed over the samples. Figure 6a shows the output SINRs of the tested methods versus the input SNRs with the fixed snapshots K = 100. It was observed that the performance of the optimal beamformer had an approximately 6 dB increment in output SINRs due to the extra paths. The performance of our proposed method was similar to that in [6,30,35] at high SNRs. The method in [6] achieved the best performance at the cost of the highest complexity compared with the others. Figure 6b depicts the output SINRs of the tested methods against the snapshot number at SNR = 30 dB and I NR = 20 dB. It was observed that the proposed method had a small impact on the number of snapshots.
severe performance degradation at high SNRs. Figure 5b depicts the output SINRs of the tested methods against the snapshot number at = 30 dB and = 20 dB. It was observed that the proposed method had a similar performance with the tested methods in [6,30,35], and the number of snapshots did not affect the performance of our proposed method.

Example 4: Mismatch Due to the Coherent Local Scattering Error
In the fourth example, the influence of the coherent local scattering mismatch on the beamformer output SINRs was considered. The coherent local scattering mismatch usually occurs in multipath propagation scenarios. Assume that the SOI is distorted by local scattering and consists of four coherent paths; the actual SV is formed as: where 0 denotes the SOI SV. ̅( )( = 1,2,3,4) denotes the coherent signal path from . are independently distributed in a Gaussian distribution drawn from a random generator, ( 0 , 4 ∘ ), in each trial.
denotes the path phase and uniformly distributed in (0,2 ) from trial to trial. and change from trial to trial, while it remains fixed over the samples. Figure 6a shows the output SINRs of the tested methods versus the input SNRs with the fixed snapshots = 100. It was observed that the performance of the optimal beamformer had an approximately 6 dB increment in output SINRs due to the extra paths. The performance of our proposed method was similar to that in [6,30,35] at high SNRs. The method in [6] achieved the best performance at the cost of the highest complexity compared with the others. Figure 6b depicts the output SINRs of the tested methods against the snapshot number at = 30 dB and = 20 dB. It was observed that the proposed method had a small impact on the number of snapshots.

Conclusions
In this paper, we proposed a low-complexity RAB method based on INCM reconstruction via subspace projection. In this method, the component of the SOI in the SCM was eliminated by replacing the eigenvector columns in the SCM such that they could span to the signal subspace with the corresponding eigenvectors in the projection matrix.
Meanwhile, the SOI SV was estimated by employing the iterative mismatch approximation method. Since the calculation of the blocking matrix can be seen as pre-processing, the complexity of INCM reconstruction only depends on a countable number of matrices' eigen decomposition and multiplication. Both the analysis and simulation illustrate that the proposed method is robust to various types of mismatches while maintaining low complexity.