Open AccessArticle

DOA Estimation for Massive MIMO Systems with Unknown Mutual Coupling Based on Block Sparse Bayesian Learning

Yang Liu

Na Dong

¹,

Xiaohui Zhang

Xin Zhao

^1,*

Yinghui Zhang

and

Tianshuang Qiu

College of Electronic Information Engineering, Inner Mongolia University, Hohhot 010021, China

Faculty of Electronic Information and Electrical Engineering, Dalian University of Technology, Dalian 116024, China

Author to whom correspondence should be addressed.

Sensors 2022, 22(22), 8634; https://doi.org/10.3390/s22228634

Submission received: 9 October 2022 / Revised: 6 November 2022 / Accepted: 7 November 2022 / Published: 9 November 2022

(This article belongs to the Section Physical Sensors)

Download

Browse Figures

Figure 1
The impact of array imperfections on the spatial spectra. "> Figure 2
The spatial spectra of the MUSIC algorithm. "> Figure 3
The spatial spectra of the BSBL algorithm and the proposed ImBSBL algorithm in an ideal environment. "> Figure 4
The spatial spectra of the BSBL algorithm and the proposed ImBSBL algorithm in the presence of mutual coupling. "> Figure 5
The comparison results of the BSBL algorithm and the proposed ImBSBL algorithm. "> Figure 6
The RMSE results of DOA estimations obtained by different algorithms versus the SNR under different numbers of snapshots. "> Figure 7
The RMSE results of DOA estimations obtained by different algorithms versus the snapshot number. "> Figure 8
The RMSE results of DOA estimations obtained by different algorithms versus the antenna number under different SNR values. "> Figure 9
The RMSE results of DOA estimations obtained by different algorithms versus the angle separation. "> Figure 10
The RMSE results of the proposed algorithm. (a) RMSE versus SNR under different snapshots. (b) RMSE versus the SNR under different array element configurations. (c) RMSE versus the snapshots number under different SNR. ">

Versions Notes

Abstract

Obtaining accurate angle parameters using direction-of-arrival (DOA) estimation algorithms is crucial for acquiring channel state information (CSI) in massive multiple-input multiple-output (MIMO) systems. However, the performance of the existing algorithms deteriorates severely due to mutual coupling between antenna elements in practical engineering. Therefore, for solving the array mutual coupling, the array output signal vector is modeled by mutual coupling coefficients and the DOA estimation problem is transformed into block sparse signal reconstruction and parameter optimization in this paper. Then, a novel sparse Bayesian learning (SBL)-based algorithm is proposed, in which the expectation-maximum (EM) algorithm is used to estimate the unknown parameters iteratively, and the convergence speed of the algorithm is enhanced by utilizing the approximate approximation. Moreover, considering the off-grid error caused by discretization processes, the grid refinement is carried out using the polynomial roots to realize the dynamic update of the grid points, so as to improve the DOA estimation accuracy. Simulation results show that compared with the existing algorithms, the proposed algorithm is more robust to mutual coupling and off-grid error and can obtain better estimation performance.

Keywords:

massive MIMO; DOA estimation; array mutual coupling; sparse bayesian learning (SBL)

1. Introduction

Massive multiple-input multiple-output (MIMO) technology has been considered one of the most promising technologies for the fifth-generation (5G) and beyond wireless communication networks, because it can obtain a high spatial multiplexing gain and can significantly improve the spatial resolution of array elements [1,2,3]. Accurate angle information obtained through the direction-of-arrival (DOA) estimation plays a crucial role in distinguishing target users from interfering users, reducing pilot pollution, and realizing accurate transmission of effective information in massive MIMO systems [4,5]. Therefore, studies on DOA estimation for massive MIMO systems have important research value and application significance in the field of wireless communications [6,7].

Traditional subspace-based DOA estimation algorithms, such as multiple signal classification (MUSIC) and estimation of signal parameters via rotational invariance techniques (ESPRIT), have been proven to provide high-precision angle estimation, but their computational complexity is relatively high in massive MIMO systems [8,9]. In view of that, numerous improved algorithms have been developed based on these algorithms in recent decades [10,11]. However, as the electromagnetic environment becomes increasingly complex and the direction-finding scenario worsens, these algorithms cause wastage of resources and the severe deterioration of direction-finding performance under low signal-to-noise ratio (SNR), small snapshots, and even large array errors. In addition, these algorithms cannot completely deal with the above-mentioned problems due to the limitation of the subspace framework.

Compressed sensing technology has good applicability and can accurately recover the original signal from low-dimension measurement data using non-zero coefficients [12]. It should be noted that, in compressed sensing technology, the sampling rate is not related to the highest frequency, and a sampling method depends on the sparsity and structure information of a signal. Recently, it has been proposed to use Bayesian compressed sensing [13] to address problems in DOA estimation. Compared with greedy and convex relaxation algorithms, sparse Bayesian learning (SBL) does not require strict information about the number of sources and can accurately reconstruct sparse signals. In view of the two major problems (i.e., small array aperture and large modeling error) existing in the off-grid SBL DOA estimation methods for nested arrays, [14] provided a novel grid refinement method based on a new data formulation to eliminate the modeling error caused by off-grid gaps. In [15,16,17], the authors developed Bayesian learning-based strategies for solving the DOA estimation problem. In [15], it has been proven that SBL can perform remarkably better than traditional beamforming and MUSIC. Since the most frequently used real-valued transformation is suitable for uniform linear arrays only, an innovative real-valued transformation with a lower computational burden was developed in [16], using the virtual steering of linear arrays. In addition, two efficient SBL-based methods for DOA estimation with arbitrary linear arrays were proposed. To diagnose the blocked elements timely under a low SNR, [17] developed an adaptive iterative diagnostic approach based on clustering block SBL, which has a high estimation accuracy in solving the coupling of DOA and antenna blockage estimation.

In practical massive MIMO systems, the performance of conventional DOA estimation algorithms degrades severely or even fails completely due to the mutual coupling between antenna elements [18,19]. The DOA estimation problem for MIMO radar systems with unknown mutual coupling was transformed into a sparse reconstruction problem, in which the off-grid effect and mutual coupling effect are considered simultaneously [20]. Then, a novel sparse Bayesian learning with the mutual coupling (SBLMC) method was proposed. To this end, a

ε

-negative meta-surface superstrate and I-shaped resonators were introduced to reduce mutual coupling [21,22]. Considering that the spatial covariance matrix (SCM), which plays an important role in DOA estimation, cannot be obtained in hybrid massive MIMO systems because the received signals cannot reach the receiver directly, the SCM reconstruction was performed using different algorithms [23,24]. Since random sampling based two-dimensional (2D) DOA estimation in massive MIMO systems can lead to matrix completion, an efficient matrix completion method based on the alternating direction method of multipliers was proposed in [25] to reduce power consumption. In [26], a new higher-order propagator method (HOPM) was designed to estimate the tensor-based signal subspace for 2D DOA estimation, which can achieve high estimation accuracy at low computational complexity. Moreover, a higher-order unitary PM (HOUPM) was proposed to improve the performance further. In [27], a deep-learning-based DOA estimation algorithm was proposed for hybrid massive MIMO systems with a uniform circular array (UCA). This system can realize comparable and, in certain cases, even superior performance at lower computational complexity compared with the maximum likelihood (ML) method. For the 2D DOA and polarization estimation of non-circular signals in polarized massive MIMO systems with uniform rectangular arrays (URAs), an efficient quaternion non-circular MUSIC (QNC-MUSIC) algorithm, which has a high estimation accuracy but low computational complexity, especially for a large source number was proposed in [28]. In [29], a spatial spectrum fusion estimation and localization (SSFEAL) algorithm for massive MIMO systems with UCAs, which can obtain decimeter-level localization accuracy and meet the requirement of future networks, was presented. To solve the two main challenges in Terahertz dynamic array-of-subarray (DAoAS) systems, [30] developed a DAoAS-MUSIC algorithm and a deep convolutional neural network (DCNN)-based algorithm for DOA estimation. Both proposed algorithms can realize super-resolution DOA estimation, and the DCNN-based algorithm can perform better than the DAoAS-MUSIC at a high SNR. In [31], an efficient method for performing the 2D DOA estimation of incoherently distributed sources with gain-phase perturbations was proposed. The proposed method can improve the estimation accuracy efficiently and also is unaffected by the array gain phase disturbance. To overcome the challenges brought by unknown mutual coupling and small samples in accurate DOA estimation for massive MIMO systems, an improved DOA estimation method based on eigenvalue comparison and real-valued MUSIC was proposed in [32]. This method can improve both the performance and the computational efficiency of DOA estimation. It is worth noting that although there have been many studies on DOA estimation algorithms for massive MIMO systems, most of them have ignored the impact of mutual coupling error on DOA estimation performance and have only slightly considered the conditions of low SNR and a small number of snapshots.

To address the shortcomings of the existing work, this paper develops an innovative algorithm based on the block SBL (BSBL) to solve the DOA estimation for massive MIMO systems with unknown array mutual coupling under the conditions of low SNR and small snapshots. Considering the existence of array errors in the actual environment, the array defects are addressed, and the DOA estimation problem is transformed into a block sparse signal reconstruction and parameter optimization problem. Then, structural information of a signal is used, and the cost function of the proposed algorithm is optimized, thus significantly improving the performance of the sparse reconstruction algorithm. The proposed algorithm is verified by simulations, and simulation results demonstrate the excellent performance of the proposed algorithm. The main contributions of this paper can be summarized as follows:

Considering the array mutual coupling error, which can severely deteriorate the DOA estimation performance, an array output vector is constructed based on the block sparsity of the original signal, where the array flow matrix of the signal is improved by the mutual coupling coefficient. Therefore, the DOA estimation problem is regarded as a block sparse signal reconstruction and parameter optimization problem.
The expectation-maximization (EM) algorithm is used to estimate hyperparameters and noise parameters to reconstruct the sparse signal. To enhance algorithm convergence speed, the boundary of the cost function is obtained by adopting the approximate approximation, and the hyperparameter representing the signal correlation is optimized, which improves the convergence speed to a large extent.
Considering the high computational complexity and low estimation accuracy of the equal-spacing angle space division, this study refines the angle space grids by finding the roots of the polynomial to realize the dynamic update of grid points in a discrete space; then a spatial screening of the updated grid is performed to achieve DOA estimation. In addition, the threshold is set to determine the grid points to be updated, thus improving the estimation efficiency. After the updating process, the grid points are closer to the real DOAs, which efficiently eliminates the off-grid error and improves the estimation accuracy. Moreover, an improved SBL-based algorithm (ImBSBL) is proposed for DOA estimation. Simulation results illustrate that compared with the existing algorithms, the proposed algorithm is more robust to mutual coupling and the off-grid error and can obtain better estimation performance.

The rest of this paper is organized as follows. In Section 2, the signal model is introduced. In Section 3, a DOA estimation algorithm for massive MIMO systems based on addressing unknown mutual coupling under a low SNR and small snapshots is introduced. The simulation results are given in Section 4. Finally, the conclusions are drawn in Section 5.

2. System Model

Consider a massive MIMO system with a uniform linear array (ULA) consisting of M antenna elements, receiving K far-field narrowband signals from distinct directions

θ_{k}

(k = 1, 2, \dots, K)

. From the perspective of spatial sparsity, the whole space can be treated as an angle set

{\{{\hat{θ}}_{i}\}}_{i = 1}^{N}

that evenly divides the DOA range

[- \frac{π}{2}, \frac{π}{2}]

into N segments,

N (N ≫ K)

, where

{\hat{θ}}_{i}

denotes the possible angle of the incident signal. When the incident signal approaches the array from angle

θ_{k} = {\hat{θ}}_{i}

, the spatial energy at that angle accumulates, forming a significant peak point, while the spatial energy at all other angles is zero. Thus, the array output vector can be expressed by:

y (t) = \bar{A} s (t) + n (t),

(1)

where

y (t)

and

n (t)

are the

M \times 1

observation matrix and observation noise matrix, respectively;

s (t)

is an

N \times 1

sparse target signal.

This makes it the best represented in the span of the base matrix,

\bar{A} \in C^{M \times N}

is the estimation problem to solve. To reconstruct the target signal accurately, the base matrix

\bar{A}

is required to be an overcomplete matrix (i.e.,

N > M

), and it is defined as follows:

\bar{A} = [a ({\hat{θ}}_{1}), \dots, a ({\hat{θ}}_{N})],

(2)

where

a ({\hat{θ}}_{i}) = {[1, e^{- j 2 π \frac{d}{λ} sin {\hat{θ}}_{i}}, \dots, e^{- j 2 π \frac{d}{λ} (M - 1) sin {\hat{θ}}_{i}}]}^{T}

and

{[\cdot]}^{T}

represents transposition.

Since a unique solution cannot be obtained when

s (t)

is solved according to Equation (1), it is necessary to constrain the signal part of measurement data. In the presence of noise

n (t)

with zero mean and variance

σ^{2}

, the DOA estimation problem can be regarded as a

ℓ_{0}

-norm optimization problem, which is defined by:

\hat{s} (t) = arg min_{s} {∥ s ∥}_{0} s . t . ∥ y - \bar{A} s ∥ \leq σ^{2} .

(3)

3. Improved BSBL-Based DOA Estimation Algorithm for Massive MIMO System with Unknown Mutual Coupling

3.1. Mutual Coupling Modeling

Mutual coupling error is common in applications using antenna arrays, and it typically occurs between adjacent elements, particularly when a system operates at a high frequency [16]. The mutual coupling matrix of a ULA can be modeled as a symmetric Toeplitz matrix, which is given by:

C = [\begin{matrix} c_{1} & c_{2} & \dots & c_{L} & \dots & c_{M} \\ c_{2} & c_{1} & c_{2} & ⋱ & ⋮ \\ ⋮ & c_{2} & c_{1} & c_{2} & c_{L} \\ c_{L} & c_{2} & ⋱ & ⋱ & ⋮ \\ ⋮ & ⋱ & ⋱ & c_{1} & c_{2} \\ c_{M} & \dots & c_{L} & \dots & c_{2} & c_{1} \end{matrix}],

(4)

where

c_{i} (i = 1, 2, \dots, M)

denotes the mutual coupling coefficient that satisfies the condition of

1 = c_{1} > \dots > |c_{L}| > |c_{L + 1}| = \dots = |c_{M}| = 0

; L is the degree of freedom of the mutual coupling matrix, which is unknown in practice. Considering the influence of mutual coupling, Equation (1) can be modified as follows:

y (t) = C \bar{A} (t) + n (t) .

(5)

Furthermore, to decompose the angle information and mutual coupling coefficient of signal Equation (5), according to the theorem [33], the following equation can be derived:

Ca = Q (a) c,

(6)

where

C

is a mutual coupling matrix, and

C = Toeplitz {c} \in C^{M \times M}

, where

c

is the vector consisting of the first-row elements of

C

;

a

is a complex vector,

a \in C^{M \times 1}

, which represents the steering vector.

To solve an

M \times L

transformation matrix

Q (a)

that can realize the decoupling of position factor and mutual coupling factor,

C

is decomposeed into

C = WZ

, where

W = [W_{1}, . ., W_{L}]

and

Z = c \otimes I

denote

M \times M L

and

M L \times M

matrices, respectively;

I

represents the identity matrix; ⊗ stands for the Kronecker product.

The

(i, j)

-th element in the l-th subblock of

W

can be calculated by:

{[W_{l}]}_{i j} = \{\begin{matrix} 1, if {[C]}_{i j} = c_{l} \\ 0, otherwise \end{matrix},

(7)

Therefore, Equation (6) can be rewritten as follows:

Ca = W c \otimes a = W [\begin{matrix} c_{1} a \\ c_{2} a \\ ⋮ \\ c_{L} a \end{matrix}] = [W_{1} a, W_{2} a, \dots, W_{L} a] [\begin{matrix} c_{1} \\ c_{2} \\ ⋮ \\ c_{L} \end{matrix}] .

(8)

According to Equation (6) and Equation (8),

Q (a)

is a function of the incident signal with

\hat{θ}

, so it can be expressed as

Q (\hat{θ})

, and then it holds that

Q (\hat{θ}) = [W_{1} a, W_{2} a, \dots, W_{L} a]

. Thus, Equation (5) can be rewritten as follows:

\begin{matrix} y (t) & = C [a ({\hat{θ}}_{1}), \dots, a ({\hat{θ}}_{N})] s (t) + n (t) \\ = [Q ({\hat{θ}}_{1}) c, \dots, Q ({\hat{θ}}_{N}) c] s (t) + n (t) \\ = D (\hat{θ}) x (t) + n (t) \end{matrix},

(9)

where

D (\hat{θ}) = [Q ({\hat{θ}}_{1}), \dots, Q ({\hat{θ}}_{N})]

is an array flow pattern matrix considering the mutual coupling factor, and the received signal of the array is given by:

x (t) = s (t) \otimes c = [\begin{matrix} s_{1} (t) c \\ s_{2} (t) c \\ ⋮ \\ s_{N} (t) c \end{matrix}] .

(10)

According to Equation (10),

x (t) = {[x_{1}^{T} (t), \dots, x_{N}^{T} (t)]}^{T}

is the signal with a block structure. In addition, since each row in

s (t)

corresponds to a signal in a potential incident direction and the element in a row is non-zero when this direction is the real DOA direction,

s (t)

indicates the row is sparse. Therefore, this paper regards the DOA estimation problem as a problem of block sparse signal reconstruction and parameter optimization.

3.2. Sparse Bayesian Solution

Since block sparsity can reduce the search degree of freedom in the signal space and can improve the reconstruction accuracy and robustness of an algorithm [34], sparse reconstruction algorithms can accurately reconstruct sparse signals. In view of this, this work employs the EM algorithm to reconstruct the sparse signal and realize the DOA estimation for the array output signal model, as shown in Equation (9). For convenience, Equation (9) can be briefly rewritten as follows:

y = D x + n,

(11)

where

x

denotes the block sparse signal.

To avoid over-matching, an SBL method can fully mine and use the prior information on data and constrain the estimated signals, which not only makes the signal more easily meet the sparse condition but also contributes to the high-precision recovery of a sparse reconstruction algorithm [35]. Assuming that the blocks are not correlated with each other; then the block signal

x = {[x_{1}^{T}, \dots x_{g}^{T}]}^{T}

can be modeled as follows:

p (x ∣ γ_{i}, B_{i}) = C N (x; 0, Γ),

(12)

where CN denotes the complex Gaussian distribution,

γ_{i}

controls the sparsity of a block,

B_{i}

captures the intra-block correlation structure [36], and

Γ = {diag}^{- 1} (γ_{1} B_{1}, \dots, γ_{g} B_{g})

. The observation noise vector

n

is assumed to have the distribution of

p (n) \sim C N (0, β^{- 1} I)

, where

β

represents the noise variance. Using the block model Equation (11), the Gaussian likelihood function with mean

D x

and variance

β^{- 1}

can be obtained as follows:

p (y ∣ x; β) = C N (y; D x, β^{- 1} I),

(13)

where

y

is the observation vector,

y \in C^{M \times 1}

, and it obeys the Gaussian distribution.

The posterior probability can be obtained by Equations (12) and (13) using Gauss identity property as follows:

\begin{matrix} p (x ∣ y; \{γ_{i}, B_{i}\}, β) \cdot p (y ∣ \{γ_{i}, B_{i}\}, β) \\ = & p (y ∣ x, β) p (x ∣ \{γ_{i}, B_{i}\}) \\ = & C N (y ∣ D x, β^{- 1} I) C N (x ∣ 0, Γ) \\ = & C N (x ∣ μ, Σ) C N (y ∣ 0, T) \end{matrix},

(14)

where

T = β^{- 1} I + D Γ D^{H},

(15)

Σ^{- 1} = Γ^{- 1} + D^{H} β D,

(16)

μ = Σ D^{H} β y

(17)

where

{[\cdot]}^{H}

represents conjugate transpose.

However, there will occur a large numerical error in calculating

Γ^{- 1}

when

γ_{i}

is too small or close to zero; thus, Equation (16) is modified as follows:

Σ = Γ - Γ D^{H} T^{- 1} D Γ .

(18)

Furthermore, to make the calculation of

μ

not related to the update of

Σ

μ

is simplified as follows:

\begin{matrix} μ & = (β^{- 1} Γ^{- 1} + D^{H} D) D^{H} y \\ = D^{H} {(β^{- 1} Γ^{- 1} + D D^{H})}^{- 1} y \\ = D^{H} Γ T^{- 1} y \end{matrix},

(19)

where

μ

and

Σ

denote the mean and covariance, respectively, and they are estimated by the EM method.

To find a sparse solution, the global minimum of the cost function [37] is calculated by:

\begin{matrix} L (\{γ_{i}, B_{i}\}, β) & = - 2 log p (y ∣ (\{γ_{i}, B_{i}\}, β)) \\ = log | T | + y^{H} T^{- 1} y \end{matrix} .

(20)

Next, using the Sherman–Morrison matrix identity formula,

{(A + XB X^{T})}^{- 1} = A^{- 1} - A^{- 1} X {(B^{- 1} + X^{T} A^{- 1} X)}^{- 1} X^{T} A^{- 1},

(21)

|A + X B X^{T}| = | B ∥ A | |B^{- 1} + X^{T} A^{- 1} X|,

(22)

where

A

and

B

are invertible matrices, Equation (20) can be written as follows:

L = (\sum_{i = 1}^{g} log |γ_{i} B_{i}| - N log | β | + log |Σ^{- 1}|) + ({β ∥ y - D μ ∥}_{2}^{2} + μ^{H} Γ^{- 1} μ) .

(23)

Then, the solution to paraments

γ_{i}

B_{i}

, and

β

is obtained. First, the partial derivative of Equation (23) is given by:

\frac{\partial L}{\partial γ_{i}} = \frac{d_{i}}{γ_{i}} - Tr [\frac{B_{i}^{- 1} (μ_{i} μ_{i}^{H} + Σ_{i})}{γ_{i}^{2}}],

(24)

\frac{\partial L}{\partial B_{i}} = B_{i}^{- 1} - γ_{i}^{- 1} B_{i}^{- 1} [Σ_{i} + μ_{i} μ_{i}^{H}] B_{i}^{- 1},

(25)

\frac{\partial L}{\partial β} = - \frac{N}{β} + Tr [Σ D^{H} D] + {∥ y - D μ ∥}_{2}^{2},

(26)

where

Tr [•]

represents trace,

μ_{i}

denotes the i-th row of

μ

, and

Σ_{i}

denotes the i-th main diagonal element of

Σ

Assume that values of Equations (24)–(26) are zeros; then, it can be written as follows:

γ_{i} = \frac{1}{d_{i}} Tr [B_{i}^{- 1} (Σ_{i} + μ_{i} μ_{i}^{H})],

(27)

B_{i} = \frac{[Σ_{i} + μ_{i} μ_{i}^{H}]}{γ_{i}},

(28)

β = \frac{N}{{∥ y - D μ ∥}^{2} + Tr [Σ D^{H} D]},

(29)

where

γ_{i}

and

β

the parameters that affect the convergence speed of the SBL algorithm.

Compared with the influence of parameter

β

on algorithm performance, parameter

γ_{i}

has a greater influence on algorithm performance; namely, different update approaches of

γ_{i}

affect algorithm performance, altering its convergence speed performance [18]. Moreover, the updated

β

and

B_{i}

both affect the reconstruction performance of the sparse recovery algorithm, but the effect of

β

is more significant. If optimal

β

cannot be obtained, the reconstruction and recovery performances of an algorithm cannot be improved even if other parameters are optimized. A number of studies have shown that the parameter

β

in an EM algorithm has poor robustness under a certain SNR because the non-diagonal elements in block elements interfere with the elements on the main diagonal [35]. Parameter

β

is obtained by

β = \frac{N}{{∥ y - D μ ∥}_{2}^{2} + \sum_{i = 1}^{g} Tr [Σ_{i} D_{i}^{H} D_{i}]} .

(30)

It should be noted that

B_{i}

affects only the local convergence performance, and the EM algorithm specifies different constraints for different blocks, which can lead to the overfitting problem. Therefore, according to the parameter average idea, so

B_{i} = B

, and it holds that:

B = \frac{1}{g} \sum_{i = 1}^{g} \frac{Σ_{i} + μ_{i} {(μ_{i})}^{H}}{γ_{i}} .

(31)

Since

B_{i} = c c^{H}

contains the mutual coupling information under unknown mutual coupling, the mutual coupling coefficient can be estimated using the first column of

B

, and the mutual coupling matrix can be obtained. Finally, when the parameters

β

and

{\{γ_{i}, B_{i}\}}_{i}

are estimated, the Maximum-A-Posteriori (MAP) estimate of

x

, which is denoted by

\hat{x}

, can be obtained by:

\hat{x} \leftarrow Γ D^{H} {(β^{- 1} I + D Γ D^{H})}^{- 1} y .

(32)

3.3. Optimized Hyperparameter Updating

In the previous subsection, a sparse Bayesian solution is described in detail, and it is explained that the EM algorithm is applied to obtain the updated parameters and sparse recovery result of a signal. However, convergence speed affects algorithm efficiency significantly. Particularly, in massive MIMO, enhancing the convergence speed can effectively reduce algorithm complexity. Numerous studies have shown that parameter

γ_{i}

affects the convergence speed of an EM algorithm [38,39]. Therefore, this section optimizes parameter

γ_{i}

and then refines the grid by finding roots of polynomials. In view of this, this study proposes an improved ImBSBL algorithm for DOA estimation, whose main aim is to find an alternative method through approximate approximation and obtain the boundary of the cost function.

The following analysis will be conducted based on Equation (20), whose first term is concave while the second term is convex. First, the upper bound of the first term is obtained by:

\begin{matrix} log |β^{- 1} I + D Γ D^{H}| & \leq log |β^{- 1} I + D Γ^{*} D^{H}| \\ + \sum_{i = 1}^{g} Tr (T^{* - 1} D_{i} B_{i} D_{i}^{H}) (γ_{i} - γ_{i}^{*}) \\ = \sum_{i = 1}^{g} Tr (T^{* - 1} D_{i} B_{i} D_{i}^{H}) γ_{i} \\ + log |T^{*}| - \sum_{i = 1}^{g} Tr (T^{* - 1} D_{i} B_{i} D_{i}^{H}) γ_{i}^{*} \end{matrix} .

(33)

where

{[\cdot]}^{*}

represents conjugation.

Then, by substituting Equation (33) into Equation (20), the following inequality is obtained:

\begin{matrix} L (γ) & \leq \sum_{i = 1}^{g} Tr (T^{* - 1} D_{i} B_{i} D_{i}^{H}) γ_{i} + log |T^{*}| \\ - \sum_{i = 1}^{g} Tr (T^{* - 1} D_{i} B_{i} D_{i}^{H}) γ_{i}^{*} \\ + y^{H} {(β I + D Γ D^{H})}^{- 1} y \\ = \tilde{L} (γ) \end{matrix} .

(34)

The abovementioned inequality indicates that

\tilde{L} (γ)

is a convex function, so its minimum can be obtained at the minimum of the variable

γ

and it should satisfy the condition of

L (γ_{min}) \leq \tilde{L} (γ_{min}) \leq \tilde{L} (γ^{*}) = L (γ^{*})

. Thus, minimizing the right side of Equation (34) will reduce the maximum value of the cost function.

Furthermore, the second term of Equation (20) is given by:

y^{H} {(β^{- 1} I + D Γ D^{T})}^{- 1} y = min_{x} [\frac{1}{β} {∥ y - D x ∥}_{2}^{2} + x^{H} Γ^{- 1} x] .

(35)

Thus, the substitution function is defined as follows:

\begin{matrix} \tilde{L} (γ) & = \sum_{i = 1}^{g} Tr (T^{* - 1} D_{i} B_{i} D_{i}^{H}) γ_{i} + log |T^{*}| \\ - \sum_{i = 1}^{g} Tr (T^{* - 1} D_{i} B_{i} D_{i}^{H}) γ_{i}^{*} \\ + min_{x} \frac{1}{β} {∥ y - Φ x ∥}_{2}^{2} + x^{H} Γ^{- 1} x \end{matrix} .

(36)

Next, let

\frac{\partial \tilde{L} (γ)}{\partial γ_{i}} = 0

; then, the update of

γ_{i}

is given by:

γ_{i} \leftarrow \sqrt{\frac{μ_{i}^{H} B_{i}^{- 1} μ_{i}}{Tr [D_{i}^{H} T^{- 1} D_{i} B_{i}]}} .

(37)

Since the sparse model is divided based on the grid, when the incident angle of a signal is not in the divided grid, there will be an off-grid error in the DOA estimation result. To reduce the influence of the off-grid error on the estimation result, this paper refines the gird using the approach of finding roots of polynomials. To obtain the grid parameter

{\hat{θ}}_{i} (i = 1, \dots, N)

, logarithmic terms irrelevant to formula derivation are ignored, and the following equation is maximized:

\begin{matrix} E_{p (x ∣ y; \{γ_{i}, B_{i}\}, β)} [log p (y ∣ x; β)] \\ = - β E_{p (x ∣ y; \{γ_{i}, B_{i}\}, β)} [{∥ y - D x ∥}_{2}^{2}] \\ = {- β ∥ y - D μ ∥}_{2}^{2} - β Tr [D Σ D^{H}] \end{matrix} .

(38)

Then, to refine

{\hat{θ}}_{i}

, let the partial derivative of Equation (38) with respect to the steering vector

V_{{\hat{θ}}_{i}} = e^{- j 2 π d / λ sin {\hat{θ}}_{i}}

be zero, which can be expressed by:

\frac{{\partial ∥ y - D μ ∥}_{2}^{2}}{\partial V_{{\hat{θ}}_{i}}} + \frac{\partial Tr [D Σ D^{H}]}{\partial V_{{\hat{θ}}_{i}}} = 0 .

(39)

According to the idea of a polynomial root, Equation (39) can be expressed in the form of a polynomial with respect to

{\hat{θ}}_{i}

using algebraic manipulation. Since Equation (39) contains the first derivative of the direction vector with respect to

V_{{\hat{θ}}_{i}}

, its polynomial form has

M - 1

roots in the complex field, and roots of the refined grid parameters should have unit amplitude. However, the roots cannot fall in the unit circle exactly due to the noise in a real environment. Assuming that

V_{{\hat{θ}}_{i}^{*}}

denotes the root closest to the unit circle, then, it can be written as:

{\hat{θ}}_{i^{*}}^{new} = arcsin (- \frac{λ}{2 π d} \cdot angle (V_{{\hat{θ}}_{i}^{*}})) .

(40)

Thus, the parameters for refining the grid are obtained. Considering that grid refinement for real DOA, which is close to or equal to the initial grid point, will increase the estimation error of the algorithm, a threshold should be defined to judge whether to perform grid refinement. The initial grid needs to be refined by

{\hat{θ}}_{i^{*}}^{n e w}

when

{\hat{θ}}_{i^{*}}^{new} \in [\frac{{\hat{θ}}_{i^{*} - 1} + {\hat{θ}}_{i^{*}}}{2}, \frac{{\hat{θ}}_{i^{*}} + {\hat{θ}}_{i^{*} + 1}}{2}]

; otherwise, the refinement is not necessary. Since it is not effective to refine all grid points in each iteration, a threshold

η

is set, and the Frobenius norm is used to select the first

η (1 \leq η \leq M)

largest elements of

μ

and obtain their indexes, which are used to determine the grid to be updated [40]. In this study, the number of sources is set to

η = M

For clarity, the proposed ImBSBL algorithm for DOA estimation with unknown mutual coupling is summarized in Algorithm 1. It is worth noting that this research does not involve the synchronization of complex dynamic networks. Different from [14], our work not only takes the off-grid factor into account, but also deals with the influence of array mutual coupling on DOA estimation. In addition, compared with [20], the proposed algorithm improves the convergence speed of the EM algorithm when estimating parameters by approximate approximation.

Algorithm 1: The Proposed ImBSBL Algorithm for DOA Estimation with Unknown Mutual Coupling and Off-grid Error.

1:: Input: Received signal $x$ ;
2:: Initialization: Initialize the number of iterations $κ$ and parameters $γ^{(κ)}$ , $B^{(κ)}$ , and $β^{(κ)}$ ;
3:: while not converged do
4:: Calculate $μ^{(κ + 1)}$ and $Σ^{(κ + 1)}$ using Equation (19) and Equation (18), respectively.
5:: Update $γ^{(κ + 1)}$ , $B^{(κ + 1)}$ , $β^{(κ + 1)}$ according to Equation (37), Equation (31) and Equation (30), respectively;
6:: Update ${\hat{θ}}_{i^{*}}^{(κ + 1)}$ by Equation (40);
7:: if $∥γ^{(κ + 1)} - γ^{(k)}∥ / ∥γ^{(k)}∥ < 10^{- 4}$ or $κ = 1000$ then
8:: break
9:: else
10:: Set $κ = κ + 1$ ;
11:: Return to Step 4.
12:: end if
13:: end while
14:: Output: $μ^{(κ + 1)}$ and ${\hat{θ}}_{i^{*}}^{(κ + 1)}$ .
15:: Perform one dimensional spectrum search on new ${\hat{θ}}_{i^{*}}^{(κ + 1)}$ to find the K maximum values to perform the DOA estimation.

4. Simulation Results

Several experiments were performed to illustrate the superiority of the proposed algorithm over the existing algorithms, including the root SBL (ROSBL) algorithm based on root finding [40], the classical high-resolution high-precision MUSIC algorithm [41], the multi-resolution focal underdetermined system solver (MFOCUSS) algorithm [42] in the convex optimization algorithm, and the sparse Bayesian array calibration (SBAC) algorithm based on the signal observation error model [43].

4.1. Effect of Mutual Coupling on DOA Estimation

The mutual coupling relationship between the elements was described by [44]

c_{n} = \{\begin{matrix} (1 + ξ) e^{j ϕ} 10 \frac{α_{c} (1 + 0.5 n)}{20}, n < 5 \\ 0, otherwise \end{matrix},

(41)

where

ξ \sim u ([- 0.05, 0.05])

φ \sim u ([0, 2 π])

and

α_{c}

are used to measure the mutual coupling effect between adjacent elements. A ULA with

M = 16

elements is used to receive

K = 2

far-field signals arriving at the incident angles of

θ = 0^{\circ}

and

θ = 25^{\circ}

. The SNR is set to 5 dB, the number of snapshots is set to 100, and the grid interval is

1^{\circ}

. The spatial spectrum of the MFOCUSS, MUSIC, OMP, and ImBSBL algorithms are presented in Figure 1, where it can be observed that the proposed ImBSBL algorithm has two sharp spectral peaks, which are close to the actual incident positions. Compared with the OMP algorithm, the MFOCUSS algorithm is more accurate but produces many smaller spectral peaks. The MUSIC algorithm performs the worst among all algorithms, having no sharp spectral peak, and the peak is far away from the real position. The reason for these results is that, among the compared algorithms, only the proposed ImBSBL algorithm considers the mutual coupling factor and off−grid error, which could severely deteriorate the DOA estimation performance. Therefore, its estimation performance is optimal.

The incident angle are set to different values, including

θ = - 15^{\circ}

θ = - 10^{\circ}

θ = 0^{\circ}

θ = 10^{\circ}

and

θ = 15^{\circ}

, while the other parameters are set as defined above. The DOA estimation spatial spectrum of the MUSIC algorithm in the presence of mutual coupling and in an ideal environment obtained by 100 Monte−Carlo runs is presented in Figure 2. The performances of the BSBL algorithm and the proposed ImBSBL algorithm in an ideal environment and in the case of mutual coupling are presented in Figure 3 and Figure 4, respectively. As shown in Figure 2, the conventional MUSIC algorithm could accurately estimate the real incident angles of signals from multiple sources in an ideal environment. However, the mutual coupling has a destructive effect and severely affects the resolution and estimation accuracy of the MUSIC algorithm. Furthermore, as presented in Figure 3, although the BSBL algorithm could accurately estimate the real incident angles of signals in the ideal environment, its performance significantly degrades in the presence of mutual coupling. Compared with the BSBL algorithm, the proposed ImBSBL algorithm is more robust to mutual coupling, and could accurately estimate the true incident angles of signals from multiple sources with a lower side lobe and smaller power leakage without pseudo peaks. This is because the mutual coupling factors are decoupled by utilizing the symmetric Toeplitz property and the banded structure of the mutual coupling matrix. Furthermore, the dynamic update of grid points by exploiting the polynomial roots can also improve the DOA estimation precision.

The comparison results regarding the running time and root mean squared error (RMSE) performance of the BSBL and the proposed ImBSBL algorithms are presented in Figure 5. As shown in Figure 5, the running times of both algorithms increase substantially with the number of antennas. However, since the parameters related to the convergence speed are optimized, the proposed ImBSBL algorithm has a faster convergence speed, especially when the antenna number is large. Moreover, as shown in Figure 5b, the RMSE of the proposed ImBSBL algorithm is slightly higher than that of the BSBL algorithm when the SNR is less than 5 dB. However, the two curves overlap when the SNR is higher than 5 dB. Thus, the proposed algorithm is more suitable than the BSBL algorithm for practice applications with a large number of antennas.

4.2. Rmse Performance Analysis

The RMSE results of different DOA estimation algorithms versus the input SNR under different snapshot are presented in Figure 6. The results are obtained under the incident angles of the signals of

- 11 . 3^{\circ}

and

3 . 7^{\circ}

and the number of snapshots of

T = 30

T = 100

, and

T = 800

. The SNR varies from 0 to 25 dB, while the other settings are the same as in Section 4.1 As shown in Figure 6, the RMSE performance is significantly improved with the increase in snapshot number for all algorithms; in addition, the proposed ImBSBL algorithm achieves higher estimation accuracy performance than other algorithms. In particular, the proposed algorithm could maintain excellent estimation performance even under low SNR and fewer snapshots. This can be explained by the fact that the block sparsity is used to model the signals, which can improve the reconstruction accuracy and robustness to noise. Moreover, the grid points are dynamically updated by finding the polynomial roots that can reduce the impact of the off-grid error on the estimation accuracy, coupling on estimation results. Since the ROSBL and MFOCUSS algorithms are not robust to mutual coupling, their performances are relatively poor. Although the SBAC algorithm could suppress array defects, its performance is inferior to that of the proposed algorithm because it is based on the point sparse model.

The RMSE results of the algorithms versus the snapshots number are presented in Figure 7. In this experiment, the snapshot number changes from 100 to 1100, the number of antennas is 64 and the other settings are the same as in Section 4.1. The results indicate that the average estimation accuracy of all algorithms increases monotonically with the snapshot number. When the number of snapshots is smaller than 220 and 250, the SBAC algorithm and the ROSBL algorithm fail to separate two signals, respectively, and their DOA estimation precisions are still rougher than that of the proposed ImBSBL algorithm.

The RMSE performance of the algorithms versus the antenna number under different SNR values are presented in Figure 8. In this experiment, the number of snapshots is 1000, the grid interval is

0 . 5^{\circ}

, the SNR is set to 0 dB and 10 dB, and the other settings are the same as in Section 4.1. When the SNR is fixed, the estimation performance of all methods improves significantly with the increase in antenna number. However, the SBAC algorithm and the proposed ImBSBL algorithm have smaller estimation errors than the ROSBL and MFOCUSS methods. This is because the SBAC and ImBSBL algorithms are robust to the effect of mutual coupling on the estimation results, so their performances are relatively good.

However, since the SBAC algorithm uses the point sparsity of signals for modeling, its performance is slightly worse than that of the ImBSBL algorithm, exploiting the block sparsity of signals for modeling. Thus, the proposed algorithm is suitable for large-scale MIMO systems. In addition, the estimation performances of all methods are improved when the SNR increases, but the proposed algorithm achieves better estimation performance than the other algorithms.

The RMSE performances of the algorithms versus the angle separation are presented in Figure 9. In this experiment, DOAs of two signals are set to

θ_{1} = - 11 . 3^{\circ}

and

θ_{2} = θ_{1} + Δ θ^{\circ}

Δ θ^{\circ}

varies from

4^{\circ}

14^{\circ}

, and the number of antennas is 64; the other settings are the same as those in Section 4.1. As shown in Figure 9, there is an RMSE performance gap between the SBAC and ROSBL algorithms, and HEIR estimation accuracies are much lower than those of the other two algorithms; however, the other two algorithms’ performances are significantly reduced when the angle interval is larger than

8^{\circ}

. In addition, regardless of the angle interval value, the proposed ImBSBL algorithm performs better than the MFOCUSS algorithm. This could be explained by the fact that the proposed algorithm is robust to mutual coupling and off-grid errors and explicitly use the sparsity structure of incident signals.

In Figure 10, the RMSE results of the proposed ImBSBL algorithm obtained by 500 independent Monte−Carlo experiments under different conditions are presented. The RMSE results versus the SNR under different snapshots numbers and array element configurations are presented in Figure 10a,b, respectively. In Figure 10c, the RMSE results versus snapshots number under different SNR are displayed. As shown in Figure 10a,c, when the SNR or snapshot number is fixed, the RMSE value of the proposed algorithm gradually decreases with the snapshots number or SNR, respectively. However, the proposed ImBSBL algorithm could maintain good estimation performance even under a small snapshot number and a low SNR. Moreover, as shown in Figure 10b, the estimation performance of the proposed algorithm could be improved by increasing the number of antennas. Consequently, the proposed algorithm achieves higher estimation accuracy under low SNR and fewer snapshots than the other algorithms; thus, it is more suitable for massive MIMO systems compared with the other algorithms.

5. Conclusions

This paper studies the DOA estimation problem in massive MIMO systems with unknown mutual coupling and off-grid error. Considering the mutual coupling between antennas, the studied problem can be regarded as a problem of block sparse signal reconstruction and parameter optimization. To solve this problem, this paper proposes an improved block sparse Bayesian learning algorithm to estimation the target DOAs. In the proposed algorithm, the cost function boundary is determined by exploiting the approximate approximation. Then, the parameters affecting the convergence speed are optimized to improve the convergence speed of the EM algorithm. Moreover, the effect of the off-grid error on the estimation results of the EM algorithm is considered, and dynamic updating of grid points is realized by finding the polynomial roots. Simulation results confirmed that the proposed algorithm is more robust to mutual coupling and off-grid error compared with the existing algorithms and can significantly outperform the existing DOA estimation algorithms.

Author Contributions

Conceptualization, Y.L. and X.Z. (Xin Zhao); methodology, Y.L., N.D. and X.Z. (Xin Zhao); validation, Y.L., X.Z. (Xin Zhao) and Y.Z.; investigation, N.D., X.Z. (Xin Zhao) and Y.Z.; data curation, X.Z. (Xiaohui Zhang); writing—original draft preparation, Y.L., X.Z. (Xin Zhao) and N.D.; writing—review and editing, Y.L., X.Z. (Xin Zhao) and T.Q.; visualization, Y.L., N.D. and X.Z. (Xiaohui Zhang); supervision, Y.L., Y.Z. and T.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China under grants 62071257 and 62161037, and is supported in part by the Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region under Grant NJYT-20-A11, and in part by the Natural Science Foundation of Inner Mongolia Autonomous Region under Grants 2021JQ07.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the support of National Natural Science Foundation of China, and thank all the reviewers for their valuable comments and suggestions during the revision of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

Zhang, J.; Zhang, J.; Ng, D.W.K.; Jin, S.; Ai, B. Improving Sum-Rate of Cell-Free Massive MIMO With Expanded Compute-and-Forward. IEEE Trans. Signal Process. 2022, 70, 202–215. [Google Scholar] [CrossRef]
Xu, K.; Zheng, F.C.; Cao, P.; Xu, H.; Zhu, X.; Xiong, X. DNN-Aided Codebook Based Beamforming for FDD Millimeter-Wave Massive MIMO Systems Under Multipath. IEEE Trans. Veh. Technol. 2022, 71, 437–452. [Google Scholar] [CrossRef]
Dicandia, F.A.; Fonseca, N.J.G.; Bacco, M.; Mugnaini, S.; Genovesi, S. Space-Air-Ground Integrated 6G Wireless Communication Networks: A Review of Antenna Technologies and Application Scenarios. Sensors 2022, 22, 3136. [Google Scholar] [CrossRef] [PubMed]
Abdallah, A.; Mansour, M.M. Efficient Angle-Domain Processing for FDD-Based Cell-Free Massive MIMO Systems. IEEE Trans. Commun. 2020, 68, 2188–2203. [Google Scholar] [CrossRef] [Green Version]
Ma, J.; Zhang, J.; Yang, Z.; Qiu, T. Off-Grid DOA Estimation Using Sparse Bayesian Learning for MIMO Radar under Impulsive Noise. Sensors 2022, 22, 6268. [Google Scholar] [CrossRef] [PubMed]
Gong, P.; Chen, X. Computationally Efficient Direction-of-Arrival Estimation Algorithms for a Cubic Coprime Array. Sensors 2022, 22, 136. [Google Scholar] [CrossRef] [PubMed]
Fang, Y.; Zhu, S.; Gao, Y.; Lan, L.; Zeng, C.; Liu, Z. Direction-of-Arrival Estimation of Coherent Signals for Uniform Linear Antenna Arrays With Mutual Coupling in Unknown Nonuniform Noise. IEEE Trans. Veh. Technol. 2022, 71, 1656–1668. [Google Scholar] [CrossRef]
Xu, W.; Chen, B.; Li, Y.; Hu, Y.; Li, J.; Zeng, Z. Dir-MUSIC Algorithm for DOA Estimation of Partial Discharge Based on Signal Strength Represented by Antenna Gain Array Manifold. Sensors 2022, 22, 5406. [Google Scholar] [CrossRef] [PubMed]
Liu, Z.; Wu, J.; Yang, S.; Lu, W. DOA Estimation Method Based on EMD and MUSIC for Mutual Interference in FMCW Automotive Radars. IEEE Geosci. Remote Sens. Lett. 2022, 19, 1–5. [Google Scholar] [CrossRef]
Jiang, X.; Qian, S. DOA estimation of coherent signals based on modified music algorithm. In Proceedings of the 2021 IEEE 3rd International Conference on Civil Aviation Safety and Information Technology (ICCASIT), Changsha, China, 20–22 October 2021; pp. 918–921. [Google Scholar]
Zhang, Q.; Liu, Y.; Long, X.; Song, K.; He, X.; Ren, X.; Qiu, T. A Cyclostationarity Based Esprit Algorithm for DOA Estimation of Uniform Circular Array. In Proceedings of the 2021 IEEE Statistical Signal Processing Workshop (SSP), Rio de Janeiro, Brazil, 11–14 July 2021; pp. 216–220. [Google Scholar]
Daponte, P.; De Vito, L.; Picariello, F.; Rapuano, S.; Tudosa, I. Compressed Sensing Technologies and Challenges for Aerospace and Defense RF Source Localization. In Proceedings of the 2018 5th IEEE International Workshop on Metrology for AeroSpace (MetroAeroSpace), Rome, Italy, 20–22 June 2018; pp. 634–639. [Google Scholar]
Lin, Z.; Chen, Y.; Liu, X.; Jiang, R.; Shen, B. Adaptive Beamforming Design of Planar Arrays Based on Bayesian Compressive Sensing. IEEE Sens. J. 2021, 21, 5185–5194. [Google Scholar] [CrossRef]
Chen, F.; Dai, J.; Hu, N.; Ye, Z. Sparse Bayesian learning for off-grid DOA estimation with nested arrays. Digit. Signal Process. 2018, 82, 187–193. [Google Scholar] [CrossRef]
Pandey, R.; Nannuru, S.; Siripuram, A. Sparse Bayesian Learning for Acoustic Source Localization. In Proceedings of the ICASSP 2021–2021 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Toronto, ON, Canada, 6–11 June 2021; pp. 4670–4674. [Google Scholar]
Dai, J.; So, H.C. Real-Valued Sparse Bayesian Learning for DOA Estimation With Arbitrary Linear Arrays. IEEE Trans. Signal Process. 2021, 69, 4977–4990. [Google Scholar] [CrossRef]
Liu, J.; Li, X.; Fan, T.; Lv, S.; Shi, M. Adaptive Antenna Diagnosis Based on Clustering Block Sparse Bayesian Learning. IEEE Commun. Lett. 2022, 26, 434–438. [Google Scholar] [CrossRef]
Xiong, X.; Zhang, M.; Shi, H.; Zhang, A.; Xu, Z. SBL-Based 2-D DOA Estimation for L-Shaped Array With Unknown Mutual Coupling. IEEE Access 2021, 9, 70071–70079. [Google Scholar] [CrossRef]
Ufiteyezu, E.; Yun, L. Dynamic Power Allocation for Cooperative Multi- Antenna Networks in Presence of Mutual Coupling. In Proceedings of the 2021 IEEE 6th International Conference on Computer and Communication Systems (ICCCS), Chengdu, China, 23–26 April 2021; pp. 846–855. [Google Scholar]
Chen, P.; Cao, Z.; Wang, X. Off-Grid DOA Estimation Using Sparse Bayesian Learning in MIMO Radar With Unknown Mutual Coupling. IEEE Trans. Signal Process. 2019, 67, 208–220. [Google Scholar] [CrossRef] [Green Version]
Tang, J.; Faraz, F.; Chen, X.; Zhang, Q.; Li, Q.; Li, Y.; Zhang, S. A Metasurface Superstrate for Mutual Coupling Reduction of Large Antenna Arrays. IEEE Access 2020, 8, 126859–126867. [Google Scholar] [CrossRef]
Bhattacharjee, S.; Ghosh, C.K. Reduction of Mutual Coupling between two adjacent Microstrip antennas using I-shaped Resonators. In Proceedings of the 2020 IEEE International Conference on Advent Trends in Multidisciplinary Research and Innovation (ICATMRI), Buldhana, India, 30 December 2020; pp. 1–5. [Google Scholar]
Fu, Z.; Liu, Y.; Yan, Y. Fast Reconstruction and Iterative Updating of Spatial Covariance Matrix for DOA Estimation in Hybrid Massive MIMO. IEEE Access 2020, 8, 213206–213214. [Google Scholar] [CrossRef]
Yan, Y.; Liu, Y.; Duan, H. Adaptive Reconstruction for Spatial Covariance Matrix in Hybrid Massive MIMO Systems. IEEE Access 2021, 9, 145060–145065. [Google Scholar] [CrossRef]
Liu, Q.; Li, X.; Cao, H. Two-Dimensional Localization: Low-Rank Matrix Completion With Random Sampling in Massive MIMO System. IEEE Syst. J. 2021, 15, 3628–3631. [Google Scholar] [CrossRef]
Ahmed, T.; Zhang, X.; Hassan, W.U. A Higher-Order Propagator Method for 2D-DOA Estimation in Massive MIMO Systems. IEEE Commun. Lett. 2020, 24, 543–547. [Google Scholar] [CrossRef]
Hu, D.; Zhang, Y.; He, L.; Wu, J. Low-Complexity Deep-Learning-Based DOA Estimation for Hybrid Massive MIMO Systems With Uniform Circular Arrays. IEEE Wireless Commun. Lett. 2020, 9, 83–86. [Google Scholar] [CrossRef]
Wan, L.; Liu, K.; Liang, Y.C.; Zhu, T. DOA and Polarization Estimation for Non-Circular Signals in 3-D Millimeter Wave Polarized Massive MIMO Systems. IEEE Trans. Wirel. Commun. 2021, 20, 3152–3167. [Google Scholar] [CrossRef]
He, D.; Chen, X.; Pei, L.; Zhu, F.; Jiang, L.; Yu, W. Multi-BS Spatial Spectrum Fusion for 2-D DOA Estimation and Localization Using UCA in Massive MIMO System. IEEE Trans. Instrum. Meas. 2021, 70, 1–13. [Google Scholar] [CrossRef]
Chen, Y.; Yan, L.; Han, C.; Tao, M. Millidegree-Level Direction-of-Arrival Estimation and Tracking for Terahertz Ultra-Massive MIMO Systems. IEEE Trans. Wirel. Commun. 2022, 21, 869–883. [Google Scholar] [CrossRef]
Tian, Y.; Liu, W.; Xu, H.; Liu, S.; Dong, Z. 2-D DOA Estimation of Incoherently Distributed Sources Considering Gain-Phase Perturbations in Massive MIMO Systems. IEEE Trans. Wirel. Commun. 2022, 21, 1143–1155. [Google Scholar] [CrossRef]
Tian, Y.; Wang, R.; Chen, H.; Qin, Y.; Jin, M. Real-Valued DOA Estimation Utilizing Enhanced Covariance Matrix With Unknown Mutual Coupling. IEEE Commun. Lett. 2022, 26, 912–916. [Google Scholar] [CrossRef]
Friedlander, B.; Weiss, A.J. Direction finding in the presence of mutual coupling. IEEE Trans. Antennas Propag. 1991, 39, 273–284. [Google Scholar] [CrossRef]
Wang, M.; Gao, F.; Shlezinger, N.; Flanagan, M.F.; Eldar, Y.C. A Block Sparsity Based Estimator for mmWave Massive MIMO Channels With Beam Squint. IEEE Trans. Signal Process. 2020, 68, 49–64. [Google Scholar] [CrossRef] [Green Version]
Zhang, Z.; Rao, B.D. Extension of SBL Algorithms for the Recovery of Block Sparse Signals With Intra-Block Correlation. IEEE Trans. Signal Process. 2013, 61, 2009–2015. [Google Scholar] [CrossRef] [Green Version]
Pan, T.; Wu, C.; Chen, Q. Sparse Reconstruction Using Block Sparse Bayesian Learning With Fast Marginalized Likelihood Maximization for Near-Infrared Spectroscopy. IEEE Trans. Instrum. Meas. 2022, 71, 1–10. [Google Scholar] [CrossRef]
Tipping, M.E. Sparse Bayesian learning and the relevance vector machine. J. Mach. Learn. Res. 2001, 1, 211–244. [Google Scholar]
Pan, Y.; Tai, N.; Cheng, S.; Yuan, N. Joint estimation of DOA and mutual coupling via block sparse Bayesian learning. In Proceedings of the 2015 IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC), Ningbo, China, 19–22 September 2015; pp. 1–6. [Google Scholar]
Wang, H.; Wang, X.; Huang, M.; Lan, X.; Wan, L. Off-grid DOA Estimation for Temporally Correlated Source via Robust Block-SBL in Mutual Coupling. In Proceedings of the 2022 Photonics and Electromagnetics Research Symposium (PIERS), Hangzhou, China, 25–27 April 2022; pp. 633–639. [Google Scholar]
Dai, J.; Bao, X.; Xu, W.; Chang, C. Root Sparse Bayesian Learning for Off-Grid DOA Estimation. IEEE Signal Process. Lett. 2017, 24, 46–50. [Google Scholar] [CrossRef] [Green Version]
Schmidt, R. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef] [Green Version]
Dou, H.; Liang, X.; Zhang, W. Research on DOA Estimation Method Based on RM-FOCUSS Improved Algorithms. In Proceedings of the 2019 IEEE 2nd International Conference on Information Communication and Signal Processing (ICICSP), Weihai, China, 28–30 September 2019; pp. 92–96. [Google Scholar]
Liu, Z.M.; Zhou, Y.Y. A Unified Framework and Sparse Bayesian Perspective for Direction-of-Arrival Estimation in the Presence of Array Imperfections. IEEE Trans. Signal Process. 2013, 61, 3786–3798. [Google Scholar] [CrossRef]
Wang, Q.; Dou, T.; Chen, H.; Yan, W.; Liu, W. Effective Block Sparse Representation Algorithm for DOA Estimation With Unknown Mutual Coupling. IEEE Commun. Lett. 2017, 21, 2622–2625. [Google Scholar] [CrossRef]

Figure 1. The impact of array imperfections on the spatial spectra.

Figure 2. The spatial spectra of the MUSIC algorithm.

Figure 3. The spatial spectra of the BSBL algorithm and the proposed ImBSBL algorithm in an ideal environment.

Figure 4. The spatial spectra of the BSBL algorithm and the proposed ImBSBL algorithm in the presence of mutual coupling.

Figure 5. The comparison results of the BSBL algorithm and the proposed ImBSBL algorithm.

Figure 6. The RMSE results of DOA estimations obtained by different algorithms versus the SNR under different numbers of snapshots.

Figure 7. The RMSE results of DOA estimations obtained by different algorithms versus the snapshot number.

Figure 8. The RMSE results of DOA estimations obtained by different algorithms versus the antenna number under different SNR values.

Figure 9. The RMSE results of DOA estimations obtained by different algorithms versus the angle separation.

Figure 10. The RMSE results of the proposed algorithm. (a) RMSE versus SNR under different snapshots. (b) RMSE versus the SNR under different array element configurations. (c) RMSE versus the snapshots number under different SNR.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Liu, Y.; Dong, N.; Zhang, X.; Zhao, X.; Zhang, Y.; Qiu, T. DOA Estimation for Massive MIMO Systems with Unknown Mutual Coupling Based on Block Sparse Bayesian Learning. Sensors 2022, 22, 8634. https://doi.org/10.3390/s22228634

AMA Style

Liu Y, Dong N, Zhang X, Zhao X, Zhang Y, Qiu T. DOA Estimation for Massive MIMO Systems with Unknown Mutual Coupling Based on Block Sparse Bayesian Learning. Sensors. 2022; 22(22):8634. https://doi.org/10.3390/s22228634

Chicago/Turabian Style

Liu, Yang, Na Dong, Xiaohui Zhang, Xin Zhao, Yinghui Zhang, and Tianshuang Qiu. 2022. "DOA Estimation for Massive MIMO Systems with Unknown Mutual Coupling Based on Block Sparse Bayesian Learning" Sensors 22, no. 22: 8634. https://doi.org/10.3390/s22228634

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

Article Menu

DOA Estimation for Massive MIMO Systems with Unknown Mutual Coupling Based on Block Sparse Bayesian Learning

Abstract

1. Introduction

2. System Model

3. Improved BSBL-Based DOA Estimation Algorithm for Massive MIMO System with Unknown Mutual Coupling

3.1. Mutual Coupling Modeling

3.2. Sparse Bayesian Solution

3.3. Optimized Hyperparameter Updating

4. Simulation Results

4.1. Effect of Mutual Coupling on DOA Estimation

4.2. Rmse Performance Analysis

5. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI