Two-Dimensional DOA Estimation Using Arbitrary Arrays for Massive MIMO Systems

With the quick advancement of wireless communication networks, the need for massive multiple-input-multiple-output (MIMO) to offer adequate network capacity has turned out to be apparent. As a portion of array signal processing, direction-of-arrival (DOA) estimation is of indispensable significance to acquire directional data of sources and to empower the 3D beamforming. In this paper, the performance of DOA estimation for massive MIMO systems is analyzed and compared using a low-complexity algorithm. To be exact, the 2D Fourier domain line search (FDLS) MUSIC algorithm is studied to mutually estimate elevation and azimuth angle, and arbitrary array geometry is utilized to represent massive MIMO systems. To avoid the computational burden in estimating the data covariance matrix and its eigenvalue decomposition (EVD) due to the large-scale sensors involved in massive MIMO systems, the reduced-dimension data matrix is applied on the signals received by the array. The performance is examined and contrasted with the 2D MUSIC algorithm for different types of antenna configuration. Finally, the array resolution is selected to investigate the performance of elevation and azimuth estimation. The effectiveness and advantage of the proposed technique have been proven by detailed simulations for different types of MIMO array configuration.


Introduction
The use of multiple antennas at the two finishes of wireless links is the consequence of the natural progression of more than four decades of advancement of adaptive antenna technology.Key administrations, for example, e-banking, e-learning, and e-well-being, will proceed to prosper and end up being more mobile.On-demand data and entertainment will steadily be conveyed over mobile and wireless communication systems.These expansions will prompt an avalanche of mobile and wireless traffic volume, anticipated to increase a thousandfold throughout the following decade [1].Late advances have shown that multiple-inputmultiple-output (MIMO) wireless systems can accomplish great increments in the overall system performance.Such frameworks are expected for the development of new generations of mobile radio systems for future wireless communication standards and applications.Subsequently, it has gotten the consideration not just of the worldwide research and development community but also of the wireless communications industry [2].There are several signal processing functions performed in MIMO systems, among which 3D beamforming for link reliability enhancement has gotten impressive consideration.To empower the 3D beamforming, precise estimation of both elevation and azimuth angles of signals is of key significance [3].
Over the last four decades, numerous high-resolution methods for estimating the DOA parameters of multiple narrow-band far-field signal sources, such as multiple signal classification (MUSIC), estimation of signal parameters via rotational invariance techniques (ESPRIT), weighted subspace fitting (WSF), and maximum likelihood (ML), have been proposed.A large number of these methods are applicable only to specific array geometries such as uniform linear arrays (ULAs), uniform circular arrays (UCAs), or uniform rectangular arrays (URAs).However, in practical cases, the limitations of an array platform usually lead to an unrealistic choice of array geometry from some specific classes.Furthermore, nonuniform array geometries are empowered to accomplish considerably enhanced resolution performance as compared to uniform array geometries with a similar number of sensors [4].
The automatic weighted subspace fitting (AWSF) algorithm for DOA estimation is demonstrated in [5].The estimation accuracy is enhanced, but no implementation has already been done in MIMO systems to assess the performance.Algorithms based on sparse representative (SR) procedures [6,7] deal with the computational complexity reduction of the 3D DOA estimation and handle the issues of the assistant calibration for smart transportation frameworks, yet it is substantial for a predefined (URA, ULA, planar array, etc.) array configuration as a vast majority of DOA algorithms [8,9].The DOA estimation task has been recast to a probabilistic system in [8,9], identifying the smallest angular regions where the approaching signal is in all likelihood.Although effective for some applications, the method is insufficient for high-resolution analysis since spatially close signals cannot be properly sensed.Therefore, a multiresolution technique has been applied to improve the system [10].Notwithstanding the positive and appealing highlights of previously mentioned approaches, every one of them shares a similar bottleneck.Undoubtedly, they need the assessment of the covariance matrix estimated from the measurements of every sensor at different snapshots.This includes a conspicuous increment in the receiver complexity and a delay in the DOA recovery [11].In [12], a two-stage full-dimension DOA estimation scheme based on the MUSIC algorithm is given.It requires an exhaustive multidimensional peak search, resulting in a relatively high computational complexity as most 2D DOA estimation algorithms.The performance of DOA estimation for massive MIMO systems is assessed in [13] using a lowcomplexity algorithm; good performance was achieved, but this method is unsubstantial for arbitrary array structure.Therefore, development of low-complexity DOA estimation algorithms for massive MIMO systems in arbitrary array configuration becomes extremely necessary for a practical implementation perspective.
In this paper, 2D Fourier domain line search MUSIC (FDLSM) is introduced to jointly estimate elevation and azimuth information for massive MIMO systems.The motivation for applying this algorithm is threefold [14,15].
(1) It reduces the computational complexity contrasted with the conventional 2D space search MUSIC or polynomial rooting strategies, by using 2D fast Fourier transform (2D FFT) that avoids the polynomial rooting step by replacing it with a computationally simple line search procedure.
(2) It provides an improved DOA estimation performance as compared to the manifold separation technique (MST).
(3) It is formulated in terms of a reduced-dimension technique throughout the whole estimation process of the data covariance matrix, thus decreasing the computational burden of estimating the covariance matrix and its eigenvalue decomposition (EVD).
This paper fundamentally concentrates on a lowcomplexity 2D DOA estimation using the 2D FDLSM algorithm with a reduced-dimension transformation in massive MIMO systems.Signal-to-noise ratio (SNR), source number, and antenna configuration are chosen to assist the evaluation of algorithm performance.Moreover, to the best of our knowledge, 2D FDLSM is used for the first time to examine the elevation and azimuth estimation in massive MIMO systems.
The rest of this paper is sorted out as follows.The system model and 2D FDLSM are described in Section 2 and Section 3, respectively.The computer simulation and performance evaluation are presented in Section 4. Section 5 concludes the paper.

System Model
We consider a MIMO system equipped with M transmit array sensors and N receive array sensors at the base station.Both transmit and receive arrays are assumed to be closely located in the arbitrary space so that any target located in the far field can be seen at the same direction by both arrays.All the array sensors are omnidirectional.We assume that the number of incident sources is known and there exist P distinct uncorrelated signals coming from directions of (θ p , ϕ p ), p = 1, 2, … , P, where θ p and ϕ p are the elevation and azimuth angle of the pth signal, respectively.The system model is depicted in Figure 1.International Journal of Antennas and Propagation At the transmit side, M transmit sensors from an M sensor (arbitrary) array and its steering vector a t θ, ϕ is given by where λ is the signal wavelength, j = −1 and (x ti , y ti , z ti ) are coordinates of the ith transmit sensor, and ⋅ T denotes the transpose.
For simplicity of notation and without loss of generality, we assume that the steering vector of the receive array is also a function of (θ, ϕ), given by where (x ri , y ri , z ri ) are coordinates of the ith receive sensor.
The output of the matched filters at the receiver is given by [16] x τ = A θ, ϕ s τ + n τ , 3 where is the vector of signal waveforms with τ denoting the time interval, n τ is the B × 1 (B = MN) complex Gaussian white noise vector of the zero mean and covariance matrix σ 2 I B , and A θ, ϕ = a 1 , a 2 , … , a P is the B × P steering matrix composed of P steering vectors with being the B × 1 steering vector of the pth signal.⊗ represents the Kronecker product.Assume that all impinging signals and noises are uncorrelated with each other.Then, the B × B data covariance matrix can be expressed as [16] where L is the snapshot number and ⋅ H stands for the Hermitian.

2D DOA Estimation Using Arbitrary Arrays for Massive MIMO Systems
In this section, a low-DOA estimation algorithm based on the FDLS MUSIC algorithm to jointly estimate the elevation and azimuth is presented.

Reduced-Dimension Beamspace (RDBS) Cramer-Rao
Bound (CRB).The length of the steering vector a p is B, which is too long and will add high computational burden.Furthermore, it will make the latter estimation of the covariance matrix in ( 5) and its EVD costly to implement because of the matrix size.So, the reduced-dimension transformation is necessary.It has been shown in [15] that there exist a transformation matrix T that makes CRB RDBS T = CRB ESP (see (5) for ESP data), which is based on the following proven theorems [15]: where Here, θ, φ , θ, ϕ , S, and ⊙ denote the estimated DOAs, the true DOAs, the signal covariance matrix, and the Schur product, respectively.
The expression of CRB RDBS T is obtained by substitution of A, W and R xx with T H A, T H W, and T H R xx T, respectively.Here, Ψ ⋅ denotes the range space of the matrix and ⊆ is read as "is contained in." The proofs of the above theorems can be found in [15] and the references therein.After various simulations, k = 4P was found to be the minimum value meeting the requirement in terms of accuracy.
According to [15], the matrix T ∈ ℂ B×k can be obtained by applying the singular value decomposition (SVD) of A W .Let A W = UΣV H  10 be the SVD of A W for the k largest singular values, where For the signal model in (11), the covariance matrix R yy can be estimated with L snapshots by The eigen decomposition of the sample covariance matrix (12) yields where the sample eigenvalues are again sorted in descending order (λ 1 ≥ λ 2 ≥ ⋯ ≥ λ k ); the matrices E s ≅ e 1 , e 2 , … , e P and E n ≅ e P+1 , λ P+2 , … , e k contain in their columns the signal and noise subspace eigenvectors of (13), respectively.Correspondingly, the diagonal matrices Λ s ≅ diag λ 1 , λ 2 , … , λ P and Λ n ≅ diag λ P+1 , λ P+2 , … , λ k are built from the signal and noise subspace eigenvalues of (13), respectively.
3.2.2D FDLS MUSIC-Based DOA Estimation.This section extends the idea of [14] into the 2D DFT to solve efficiently (14).The procedure is as follows: By combining (4) and ( 13) and using the transformation matrix T, the MUSIC null spectrum in [17] can be rewritten as

14
where ∥⋅ ∥ is the vector 2-norm.From ( 14), it can be seen that if θ, ϕ is the desired DOA, f θ, ϕ will be equivalent to zero.Given that f θ, ϕ is a periodic function in θ and ϕ with the period 2π, it can be expressed using the finite 2D Fourier series as [18] G mn e jmθ e jnϕ , 15 where the Fourier coefficients are given by with k 1 = q e Δθ, k 2 = q a Δϕ, Δθ = 2π/ 2Q e − 1 , and Δϕ = 2π/ 2Q a − 1 , where 2Q e − 1 and 2Q a − 1 represent the sample number in elevation and azimuth, respectively.We notice that 2D FFT can be used to obtain the Fourier coefficient matrix G ∈ C 2Q e −1 × 2Q a −1 rapidly.Once the Fourier coefficients are acquired, to ameliorate the resolution, zero padding is applied efficiently as follows: Generally, we pick J 0 ≫ 2Q e − 1 and J 1 ≫ 2Q a − 1.The new cost function f θ, ϕ of every grid can be computed effectively by applying 2D IFFT to the zero padded Fourier coefficient matrix G ∈ C J 0 ×J 1 in (17).Accordingly, no polynomial rooting is required any longer as the polynomial rooting step is replaced by a simple line search with θ and ϕ.
Recall that 1/f θ, ϕ will assume a very large value when θ, ϕ is equal to the DOA of one of the signals.Thus, the P signals can be easily obtained by taking the P greatest peaks of 1/f θ, ϕ .
We also recall that the computational complexities of the conventional 2D space search MUSIC algorithm and the 2D FFT-based technique are entirely unequal.For the conventional 2D space search MUSIC algorithm, to obtain suitable accuracy, the intrigued area has to be divided into exceptionally dense grids.For instance, if the intrigued area is θ ∈ 0 ∘ , 45 ∘ and ϕ ∈ −50 ∘ , 50 ∘ with 0 1 ∘ × 0 1 ∘ as the unit grid, we have to compute 451 × 1001 = 451451 times the matrix-vector product ∥E H n T H a θ, ϕ ∥ 2 in ( 14).However, in the 2D FFT technique, only 2Q e − 1 × 2Q a − 1 matrixvector products in (14) need to be computed (Q e = Q a = 77 was used for our test).For other grids, zero padding is applied straightforwardly; hence, the computational complexity is decreased significantly.Furthermore, this technique has a lower complexity than the 2D polynomial rooting strategy.For instance, to estimate the 2D DOAs in [19], additional time is required to solve the high-order bivariate nonlinear polynomials.Since rooting a bivariate nonlinear polynomial is costly and the 2D FFT/IFFT method for real-time implementation is promptly accessible, the 2D FFT-based approach comes out to be competitive.

4
International Journal of Antennas and Propagation 3.3.Summary.The major steps of the proposed algorithm for 2D DOA estimation in arbitrary geometry for the MIMO system are as follows: (1) Obtain the RDBS matrix T from the SVD of (10) (for k = 4P); then, apply the mapping x↦y = T H x to acquire the new representation of the array observations (11).
(2) From (11), construct the covariance matrix R yy in (12) and perform the eigenvalue decomposition to obtain the noise subspace E n .
(5) Employ the zero padding technique described in (17) to shape the extensive matrix G ∈ C J 0 ×J 1 .
(6) Apply the 2D IFFT to G to obtain the new cost function f θ, ϕ .
(7) Estimate the desired DOAs from the null spectrum 1/f θ, ϕ by taking the P largest peaks.
3.4.Computational Complexity.The order of computational complexity of the proposed FD line search MUSIC for MIMO systems with a reduced-dimension matrix is compared with that of the case where the reduced-dimension matrix is not used.The results are shown in Table 1.
The last row of Table 1 shows the total number of operations (addition and multiplication) needed for both cases (i.e., B = 256, P = 5, L = 100, and J 0 = J 1 = 2000).We can easily see that, despite the fact that our approach uses the SVD, the computational complexity is reduced by 53% as compared with the case where the mapping matrix is not used.Table 2 gives the elapsed time by the two methods, from which it can be easily seen that the enhancement of the speed of our method is considerable.

Simulation Results
To get a more quantitative understanding of how the 2D Fourier domain (FD) line search MUSIC algorithm performs in massive MIMO systems, Monte Carlo simulations are conducted in this section.In the simulation, the total number of antenna elements is set to 256, which is considered one possible dimension of massive MIMO [1].The case where there are 5 sources is used with ϕ = −70 ∘ , −40 ∘ , −13 ∘ , 60 ∘ , 80 ∘ and θ = 40 ∘ , 31 ∘ , 60 ∘ , 70 ∘ , 50 ∘ .The simulation parameters are listed in Table 3.

Performance Evaluation.
The root-mean-squared error (RMSE) has been used to measure the effectiveness of our proposed 2D DOA estimation algorithm, which is defined as [13] where θp,n is the estimated elevation angle θ p at the nth Monte Carlo trial, which is similar for φp,n .P is the estimated source number and N c represents the Monte Carlo simulation times.The results are shown in Figure 3.
From the RMSE results for massive MIMO DOA estimation that are illustrated in Figure 3, we can see that the RMSE values of our proposed 2D FDLS MUSIC are

Algorithm
Using reduced-dimension matrix Without using reduced-dimension matrix

Number of Monte Carlo 100
Directivity of the base station antennas Omnidirectional 5 International Journal of Antennas and Propagation below 0 065 with SNR ranging from 0 dB to 30 dB.Furthermore, the DOA can be clearly observed; therefore, we can conclude that the proposed algorithm works well in massive MIMO systems.Different curves correspond to different types of antenna configuration as indicated in the legend.As the SNR increases, the performance of each array configuration becomes much more better.It is interesting to see that the 64 × 4 and 4 × 64 arrays have produced excellent quality compared with other array configurations.
Figure 2 also compares our proposed algorithm with the 2D MUSIC (using uniform circular array (UCA) structure).It shows that for SNR ranging from 10 to 30 dB, both algorithms provide almost the same accuracy (for 64 × 4, 4 × 64, 8 × 32, and 32 × 8 configurations) while our approach uses less complexity as explained in Section 3.2.
From Figure 4, we can clearly observe that for the number of sensors varying from 50 to 256 elements, our approach still offers better performance than the traditional 2D MUSIC (UCA and SNR = 20 dB were used).

Performance Comparison between Elevation and
Azimuth Estimation. Figure 3 shows that two identical antenna configurations can have different performance.To better understand it, we have plotted the RMSE comparison between elevation and azimuth estimation of five sources with four types of array configuration (notice that 16 × 16 array configuration is not included due to its higher RMSE values); the results are shown in Figure 5.
It is easy to see from this figure that the overall performance of angle estimations weighs more in some cases while    8.It shows that the DOA estimation performance is getting better with the increase in the number of snapshots.That is because the error in estimating the signal covariance matrix in ( 12) is reduced with the increase in the snapshot number.

Conclusion
In this paper, a 2D DOA estimator in arbitrary array geometry for massive MIMO has been analyzed and modeled based on low-complexity 2D Fourier domain line

(Figure 1 :
Figure 1: System model of 2D massive MIMO with an MxN arbitrary array configuration.

2
spans the range space of A W , one may take T = U. Include the B × k matrix T, P < k ≤ B, and the mapping x↦y = T H x from ESP to RDBS, a new set of observation result, consisting of k-dimensional vectors.Then, in RDBS, the representation of the signals received by the array will be m u t h ( d e g r e e ) E le va ti o n (d eg re e) (a) 64 × 4 configuration 2D FDLS MUSIC for MIMO 2D Spectrum (dB) A z i m u t h ( d e g r e e ) E le v at io n (d eg re e) m u t h ( d e g r e e ) E le va ti on (d eg re e) (c) 32 × 8 configuration 2D FDLS MUSIC for MIMO 2D Spectrum (dB) A z i m u t h ( d e g r e e ) E le v a ti o n (d e g re e ) 8 × 32 configuration 2D FDLS MUSIC for MIMO 2D Spectrum (dB) A z im u t h ( d e g r e e ) E le v at io n (d e g re e )

Figure 2 :
Figure 2: 2D FDLS MUSIC spectrum for MIMO in different antenna configurations.

Figure 5 :
Figure 5: RMSE comparison between elevation and azimuth of five sources with eight types of array configuration.
Estimation results of elevation angle

Figure 8 :
Figure 8: RMSE for massive MIMO DOA estimation versus snapshot number for SNR = 10 dB with different types of array configuration.

6
International Journal of Antennas and Propagation the accuracy of azimuth/elevation angles is more valued in other cases.Thus, the antenna array configuration can be adjusted as an efficient way to better satisfy the requirements of DOA estimation.Based on the simulations above, the Figure 3: RMSE comparison between 2D MUSIC and 2D FDLS MUSIC algorithms for massive MIMO DOA estimation versus SNR for five sources with different types of array configuration.
Figure 7: RMSE for massive MIMO DOA estimation versus source number for SNR = 10 dB with different types of array configuration.