Next Article in Journal
A Method for Determination and Compensation of a Cant Influence in a Track Centerline Identification Using GNSS Methods and Inertial Measurement
Previous Article in Journal
Effect of Nanodisks at Different Positions on the Fano Resonance of Graphene Heptamers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of CS-Based Channel Estimation for Millimeter Wave Massive MIMO Systems

College of Communications Engineering, Army Engineering University of PLA, Nanjing 210007, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(20), 4346; https://doi.org/10.3390/app9204346
Submission received: 25 August 2019 / Revised: 28 September 2019 / Accepted: 11 October 2019 / Published: 15 October 2019

Abstract

:
Compressed sensing (CS) has great potential in channel estimation for millimeter wave (mmWave) massive multiple-input multiple-output (MIMO) systems. To solve such a CS-based channel estimation problem, three categories of algorithms, namely convex relaxation algorithms, greedy iteration algorithms and Bayesian inference algorithms, are widely used. In this paper, with a unified massive MIMO framework, comprehensive comparisons among three categories of algorithms are presented in the perspective of the estimated accuracy, which is affected by the received signal-to-noise ratio (SNR), the number of resolvable paths, angular quantization error, the number of pilot symbols and hardware impairments. Specifically, it shows that convex relation algorithms achieve the best estimation accuracy at the high SNR range and it is mainly affected by the received SNR and transmitter’s hardware impairments. At the low SNR range, greedy iteration algorithms outperform others and the estimated accuracy is then limited by the angle quantization error. Furthermore, a tradeoff between the estimated error and complexity is achieved in Bayesian inference algorithms, a;though its estimated error is sensitive to the number of available pilot symbols.

1. Introduction

Millimeter wave (mmWave) communication has been widely considered as a key technology for B5G due to its potential to provide gigabits-per-second data rates by exploiting the large unused bandwidth [1]. Unfortunately, it is challenging for long range wireless communication because of the huge paths loss in mmWave band. To compensate the paths loss and extend the coverage, massive MIMO antennas are supposed to be deployed in mmWave communication systems. However, in fell-digital MIMO structure, the radio frequency (RF) chains will increase as the antennas become large scale, which results in unaffordable hardware cost and power consumption [2]. Thus, massive MIMO with phase shifter network based hybrid precoding and electromagnetic lens based beamspace has been considered [1,2,3]. Generally, it is difficult for such an efficient structure to obtain the full channel state information (CSI), since the limited number of RF chains only combines the measurement signal with smaller dimension compared to the entire channel matrix. Meanwhile, many pilot symbols need to be assigned during the channel estimation phase, leading to overwhelming pilot overhead. Fortunately, some inherent properties including low rank and structural sparsity of mmWave massive MIMO channels have been demonstrated in practical measurement environments [4,5]. By exploiting low rank and structural sparsity, some channel estimation algorithms have been proposed to efficiently recover the sparse channel matrix [6,7,8,9,10,11,12,13,14], including blind, semi-blind and pilot-aided estimation techniques.
For blind and semi-blind channel estimation techniques, the basic idea is to exploit some inherent features in a modulated signal or use a decision feedback method rather than many pilots to acquire the CSI [6,7,8]. Specifically, the blind 1-norm regularized channel estimation algorithm was proposed to solve the narrowband channel estimation problem in [6]. In [7], a novel blind signal detection scheme that simultaneously estimates the channel and data was proposed. In [8], based on certain available prior data obtained from the output of a soft-input soft-output channel decoder, a data-aided channel estimation scheme was proposed. In such a channel estimation scheme, the CSI can be coarsely estimated via a few pilot symbols, and then its precision can be further improved by data symbols. However, for such blind or semi-blind channel estimation schemes, only after a long time receiving the data symbols can the CSI be estimated accurately. In addition, the great computational complexity will lead to a long time delay, thus such algorithms cannot be applied to estimate the channel with relatively short coherence time or meet the demand of low time delay.
To efficiently acquire the CSI, pilot-aided channel estimation algorithms are widely studied [9,10,11,12,13,14], which can be divided into three categories, i.e., convex relaxation algorithms, greedy iterative algorithms and Bayesian inference algorithms. Specifically, based on the frequency-flat geometric channel model, the CSI of complex gains, angles of arrival (AoAs) and angles of departure (AoDs) were estimated through gradient descent algorithm in [9]. In [10], by exploiting sparse and low-rank structures of the mmWave massive MIMO channel with angle spread, the two-stage compressed sensing algorithm was proposed to reduce the sample complexity. In [11], the angle spread at receiver was estimated by greedy iterative algorithms and the normalized mean squared error (NMSE) was compared among different greedy iterative algorithms. In [12], the frequency selective channel was recovered with orthogonal matching pursuit (OMP) algorithm in the single carrier (SC) system and orthogonal frequency division multiplexing (OFDM) system, respectively. In [13], the off-grid sparse Bayesian learning algorithm was proposed to estimate CSI in an uplink multiuser mmWave massive MIMO system. In [14], the Bayesian compressed sensing (BCS) algorithm was adopted to recover the channel matrix based on a geometric channel model, where the transmitter’s hardware impairments were taken into consideration. Note that the prior works have investigated these channel estimation algorithms based on different assumptions and drawn some conclusions. However, these conclusions are only for specific parameter ranges under different scenarios. To the best of our knowledge, there is no comprehensive comparison among CS-based channel estimation algorithms with a unified framework for mmWave massive MIMO systems.
Inspired by the above analysis, in this paper, we formulate the sparse channel estimation as a normal CS problem and provide a general overview of the current mmWave massive MIMO channel estimation approaches. Different from the previous works, the purpose of our work is a comprehensive comparison among existing algorithms with a unified framework, including their basic assumptions, application conditions and key results. The main contributions and key results of this work are listed as follows.
  • In this paper, with a unified massive MIMO framework, we compare the NMSE performance among three categories of algorithms, which is affected by the received SNR, the number of resolvable paths and pilot symbols, angular quantization error, hardware impairments and computational complexity. Through comprehensive comparison, the characteristics and application conditions of each algorithm are revealed and the factors that affect the estimated error and computational complexity are also presented.
  • Through theoretical analysis and simulation results, we show that convex relation algorithms achieve the best estimation accuracy at the high SNR range and it is mainly affected by the received SNR and transmitter’s hardware impairments. At the low SNR range, greedy iteration algorithms outperform others and the estimated accuracy is then limited by the angle quantization error. Furthermore, a tradeoff between the estimated error and complexity is achieved in Bayesian inference algorithms, although its estimated error is sensitive to the number of available pilot symbols.
  • We also analyze the overall computational complexities of three categories of algorithms and visually represent them by the running time. Through illustrating the runtime of different algorithms versus the sparseness, we show that the computational complexity in the convex relaxation algorithm is the highest, and it even squarely increases with the sparseness in the gradient descent-based convex algorithm, while that in greedy iteration algorithms is minimum and grows linearly with the sparseness. In contrast to them, the computational complexities of Bayesian inference algorithm decreases as the sparseness increases.
The rest of the this paper is organized as follows. After describing the system model and channel model in Section 2, we formulate the channel estimation as a compressed sensing (CS)-based sparse signal recovery problem in Section 3. A comparison of sparse signal recovery algorithms, namely convex relaxation algorithms, greedy iteration algorithms and Bayesian inference algorithms, are presented in Section 4. Simulation results are provided in Section 5. Finally, we draw conclusions in Section 6.
Throughout our discussions, bold uppercase A denotes a matrix and lowercase a denotes a vector. Superscripts A * , A T , A H , and A 1 denote the conjugate, the transpose, the conjugate transpose, and the inverse of a matrix A , respectively. v e c A denotes a vector obtained through the vectorization of a matrix A . For M × N matrices A and B , A B denotes the M 2 × N 2 matrix of Kronecker product between A and B . a 0 , a 1 , and a 2 are the l 0 , l 1 , and l 2 norms, respectively.

2. System Model

2.1. System Model

Without loss of generality, we consider a typical mmWave point-to-point hybrid massive MIMO system, as shown as Figure 1. Let N T and N R denote the number of antennas at the transmitter and receiver, respectively. For simplicity, we assume the number of RF chains deployed at both sides is equal, i.e., N R F T = N R F R = N R F . Since the RF chains are costly and due to power consumption, the number of RF chains is much smaller than that of antennas, i.e., min N T , N R N R F .
For such an mmWave massive MIMO system with hybrid structure, the transmitter and receiver communicate via N s streams where N s N R F . Specifically, the transmitter employs a baseband precoder F B B C N R F × N s followed by an analog phase shifter precoder F R F C N T × N R F . After the combined transmitter precoding matrix F = F R F F B B C N T × N s , the transmitted signals can be written as x = Fs C N T × 1 , where the ith element of x is the transmitted signal at the ith antenna, s is the transmitted symbol vector which satisfies E s · s H = P N s I N s , and P is the total power of transmitted symbols [12]. During the channel estimation phase, it is worth noting that the pilot power and pilot pattern also affect the performance of sparse channel estimation, which need to be carefully designed. However, to simplify expression, we assume the pilot symbols are evenly placed and transmitted with equal power.
At the receiver, the received signals on all antennas are combined with a combiner W C N R × N s composed of the RF combiner W R F C N R × N R F and baseband combiner W B B C N R F × N s , which can be expressed as
y = W H HFs + W H n
where H C N T × N R denotes the channel matrix discussed in the next subsection and n C N 0 , σ 2 I N s denotes the symmetric complex Gaussian distributed additive noise vector at the receiver.
In a practical communication system, the non-ideal characteristics of components, including nonlinear power amplifier [14], inphase/quadrature-phase imbalance [15], phase noise [16] and carrier frequency offset [17], will greatly affect the estimation performance. Especially in mmWave band, the use of cheap equipment is encouraged, which leads to unignorable hardware impairments. To analyze the effect of hardware impairments on the performance of channel estimation algorithms, we take the residual transmitter hardware impairments into consideration and model them as the additive distortion noise [18]. Specifically, the additive distortion noise denoted by e I C N T × 1 can be modeled as a Gaussian distributed random variable and its power is proportional to the pilot symbols power. That is, e I C N ( 0 , Λ ) and Λ = Δ k e q 11 , , q N T N T T , where q i i denotes the ith element of covariance matrix Q = Δ E x x H and k e denotes the proportionality coefficient [19]. In addition, we assume that the hardware impairments at the transmitter are uncorrelated with the transmitted signals x and ignore the hardware impairments at the receiver as they can be simply treated as additive noise during the channel estimation phase. Then, the combined signals at the receiver can be modeled as
y = W H H Fs + e I + W H n

2.2. Channel Model

In a practical massive MIMO system, the transmitters, such as base stations, are typically located at high elevation with fewer surrounding scatterers, while the receivers, such as user equipment, are surrounded by rich scatterers at low elevation. Thus, each path departing from the transmitter may arrive at several reception directions [11]. Moreover, the channel is time-varying as the multipath delay and frequency offset are common and it is often useful to represent the channel response in frequency domain, which can be expressed as [20,21]
H t , f = l = 1 L q = 1 L A O A a l , q e j 2 π υ l , q t τ l , q f a R θ l , q R a T H θ l T
where L denotes the number of paths that depart from the transmitter with different AoDs; L A o A denotes that each of departure path will contribute L A o A reception directions; a l , q , υ l , q , τ l , q and θ l , q R , respectively, denote the complex gain, Doppler shift, delay and the AoA of the qth arrival path of the lth departure path; θ l T 0 , 2 π denotes the AoD of the lth path; and a T θ l T and a R θ l , q R denote the antenna array response vector at the transmitter and receiver, respectively. For a uniform linear array (ULA) antenna, the array response vectors at transmitter and receiver can be presented as
a ( θ l , q T ) = 1 , e j 2 π d sin θ l , q T / λ , , e j 2 π ( N 1 ) d sin θ l , q T / λ T a ( θ l R ) = 1 , e j 2 π d sin θ l R / λ , , e j 2 π ( N 1 ) d sin θ l R / λ T
where λ is the wavelength of signal and d is the antenna spacing, which is typically assumed as d = 0.5 λ .
It is difficult to estimate such a channel matrix as five parameters need to be estimated in each path: the complex gain a l , the angle of arrival θ l , q R , the angle of departure θ l T , the Doppler shift υ l and the delay τ l . Besides, although there are only 2–8 paths on average in the mmWave band [22], the actual number of paths also needs to be estimated, which makes it more difficult to estimate the CSI. Therefore, it is necessary to simplify the channel model according to the practical scenarios.
Suppose that the channel is slowly varying over the signal duration of interest T, which means the channel is quasi-static and the Doppler frequency shift can be ignored [21,23]. Based on such assumption, we can rewrite Equation (3) as
H = l = 1 L q = 1 L A O A a l , q e j 2 π τ l , q f a R θ l , q R a T H θ l T
In addition, if we further assume that the bandwidth B of channel is sufficiently small (i.e., B τ l 1 l , l = 1 , , L ), the channels of mmWave massive MIMO system can be simplified as a narrowband frequency-flat physical channel model [11]. Then, Equation (5) can be simplified as
H = l = 1 L q = 1 L A O A a l , q a R θ l , q R a T H θ l T
Although there are abundant scatterers surrounding the ends, the penetration and reflection loss are also significant at mmWave band. Thus, a geometric channel model is widely used in the prior work [9,10,24], which can be expressed as
H = l = 1 L a l a R θ l R a T H θ l T
To simplify the expression, we adopt the multipath narrowband frequency-flat channel model in Equation (7) to compare three categories channel estimation algorithms. Note that Equation (7) can be further expressed by a more compact way
H = A R Δ d A T H
where Δ d C L × L is a diagonal matrix that contains non-zero elements of channel gains, and A R C N R × L and A T C N T × L contain the columns a R θ l R and a T θ l T , respectively.
With the multipath sparsity of physical channel, a discrete virtual angular domain (VAD) channel model is widely used [9,10,11,12,13]. Assume that the AoAs and AoDs are taken from a uniform grid of size G, i.e., θ l T , θ l R 0 , 2 π G , , 2 π G 1 G with G L . Then, Equation (8) can be rewritten as
H = A ˜ R Δ d v A ˜ T H + E Q
where Δ d v C G × G is a L-sparse VAD channel matrix with L non-zero elements in the positions corresponding to the AoAs and AoDs, A ˜ T C N T × G and A ˜ R C N R × G are the transformation dictionaries that conation the array response vectors corresponding to the angles in the grid, and E Q denotes the quantization error that is bounded, as analyzed in [25]. In addition, a suitable choice of G is crucial as a larger G can reduce the quantization error but increases the computational complexity, vice versa. Typically, it is assumed that the number of grids G is larger than the number of antennas [25].

3. Formulation of the Channel Estimation Problem via Compressed Sensing

In this section, we formulate a technique via CS theory to estimate the complex gain a l , the angle of arrival θ l R and the angle of departure θ l T with a unified framework, leveraging the sparsity of mmWave channel.
To formulate a sparse signal recovery problem of channel estimation, we substitute Equation (9) into Equation (2) and use the result of Kronecker product (i.e., v e c ( ABC ) = ( C T A ) v e c ( B ) ) to transform it as
v e c y = Fs + e I T W H ( A ˜ T * A ˜ R ) v e c ( Δ d v ) + ( Fs + e I ) T ( W ) H v e c ( E Q ) + v e c W H n
According to CS theory, sufficient measurement vectors are needed to guarantee the sparse vector recovered with high probability. Assume that M training frames are occupied and both sides use pseudorandomly precoders and combiners during the channel estimation phase. After collecting the combined signals at the receiver, Equation (10) can be extended to a standard expression based on CS theory.
v e c y 1 v e c y M Y = F 1 s 1 + e I 1 T W 1 H F M s M + e I M T W M H Φ A ˜ T * 1 A ˜ R 1 A ˜ T * M A ˜ R M Ψ v e c ( Δ d v ) 1 v e c ( Δ d v ) M h + v e c e N 1 v e c e N M E N
where Φ m = Δ ( F m s m + ( e I ) m ) T ( W ) m H , m 1 , M is defined as the measurement matrix, which is assumed to be known at the receiver for standard CS-based model. However, the transmitter hardware impairments will prevent the receiver from fully knowing the real measurement matrix. Ψ m = Δ A ˜ T * m A ˜ R m is defined as the orthogonal dictionary matrix. h m = Δ v e c ( Δ d v ) m C G 2 × 1 is a sparse vector that contains all complex channel gains. e N m = Δ ( F m s m + e I m ) T ( W m ) H v e c ( E Q m ) + v e c W m H n m is defined as the additive noise vector that contains the quantization error, hardware impairments and additive white Gaussian noise (AWGN) at the receiver.
For such a sparse signal recovery problem, it has been widely accepted [9,10,11,12,13,14] that the parameters, including the complex gain a l , the angle of arrival θ l R and the angle of departure θ l T , can be estimated by the following optimization problem
m i n h h 0 s . t . y Φ Ψ h 2 2 < ε
where ε is a parameter that denotes the tolerance of estimation error and the variance of the noise is a suitable choice for it.

4. A Comparison of Sparse Signal Recovery Algorithms

In this section, the properties of different algorithms are revealed and the factors that affect the estimated accuracy and computational complexity are also presented. Further, the estimated error and complexity are compared for three kinds of algorithms.

4.1. Convex Relaxation Algorithms

(1) CVX-based Algorithm: It is hard to solve an optimization problem when the objective function to minimize the l 0 norm is nonconvex. Fortunately, l 1 norm has been proved to approximate it well and widely used in prior works [9,10,11]. By transforming the objective function from l 0 norm to l 1 norm, the optimization problem in Equation (12) can be formulated as a convex optimization problem, which can be expressed as
h ^ = arg min h h 1 + λ y Φ Ψ h 2
where λ controls the tradeoff between the sparseness and estimation error and it is a key parameter that affects the rate of convergence, which needs to be updated in each iteration.
To work out such a convex optimization problem, there are numerous toolboxes such as CVX and SeDumi [26,27]. For the existing toolboxes, we can directly use them since the algorithms have been integrated into software such as MATLAB. However, the toolboxes are restricted by the hardware because the size of channel matrix may exceed the hardware memory when the transmitter and receiver deploy massive MIMO antennas. In addition, transforming the sparse channel recovery as a convex optimization problem will cause some performance loss and the quantization error of angles will limit the estimated performance. Besides, some extra procedures are included in the common toolboxes, which will increase the computational complexity of algorithm.
(2) IR-based Algorithm: Except the existing toolboxes, authors in [9] adopted the gradient descent and iterative reweighting (IR) method to solve the minimum l 0 norm problem. By replacing the objective function with a log-sum function, the optimization problem in Equation (12) can be rewritten as
min h , θ T , θ R S i h , θ T , θ R = Δ λ 1 h H D i h + y Φ Ψ h 2 2
where D i is a diagonal matrix and its elements are the estimation of path gains at the ith iteration. Then, by defining the maximal tolerance of channel gain error and the corresponding optimal value of objective function S at the ith iteration, we can update the θ ^ R ( i + 1 ) and θ ^ T ( i + 1 ) and compute the residue in each iteration via gradient descent method until the error is less than the threshold or iteration reaches the maximum.
By using the gradient descent and IR method, the quantization error of angles and limited hardware memory can be effectively avoided. However, the complexity increases as complex matrix operations such as derivation are involved, which cannot meet the requirement of low delay. To reduce the computational complexity, some preconditioning procedures can be adopted such as the suitable choice of initial candidates and the iteration step size [9].

4.2. Greedy Iterative Algorithms

The typical representation of the greedy iterative is OMP algorithm. Its main idea is to find the most relevant column contained in the measurement matrix to the residue in the ith iteration [25].
j i = arg max j Ω ^ i 1 r i H Φ Ψ j 2
where r i is the residual vector, Φ Ψ j denotes the jth column of the measurement matrix and Ω ^ denotes the estimation of the common support. Ω ^ can be updated by
Ω ^ i = Ω ^ i 1 j i
After the sparse common support Ω is estimated, the least square estimation algorithm can be used to estimate the full channel, which is given as follows:
h ^ = Ω H Ω 1 Ω y
Except OMP algorithm, there are some modified greedy iteration algorithms such as subspace pursuit (SP) algorithm to solve such a problem [11]. For SP algorithm, the relevant L columns are selected in each iteration to reduce iterations. For greedy iterative algorithms, a suitable criterion to end the procedure is crucial. If the sparseness of channel is known in advance, the algorithms will be stopped when the iterations reach the sparseness; otherwise, the maximum iteration will be set larger than the real channel sparseness to ensure the recovery accuracy. However, with regard to SP algorithm, there is an extra threshold to exceed since numerous columns are selected in each iteration.
The defect of greedy iteration algorithms is that they cannot avoid the quantification error of angle even when G is large enough. Especially at the high SNR range, the quantization error is a key factor that plays a decisive role. This indicates that the performance based on greedy iterative algorithms is similar because the quantification error of angle is a bottleneck for such a sparsity channel estimation problem.
Although the performance loses a lot by replacing the real angle of arrival and departure by discrete grid angle, the computational complexity is relatively lower than convex relaxation algorithms. This is due to the fact that the greedy iterative algorithms do not involve the complicated matrix operations such as derivation. Hence, the greedy iterative algorithms are candidates to solve such a sparse channel estimation problem as the low delay requirement for next wireless communication can be satisfied well.

4.3. Bayesian Inference Algorithms

The Bayesian inference algorithms are based on the Bayesian theorem and their main idea can be summarized as follows [28]. Firstly, assume the sparsity prior distribution of parameters that need to be estimated. Secondly, transform the prior to posterior distribution based on Bayesian theorem. Thirdly, find the optimal value of parameters to maximize the posterior probability. The prior assumption is crucial as it controls the computational complexity and recovery accuracy of algorithms. To enlarge the sparseness and formulate a closed-form expression after Bayesian inference, the hierarchical sparsity prior is widely used in the sparsity Bayesian learning, particularly in the relevance vector machine (RVM) [28]. Specifically, define a zero-mean Gaussian prior of each element in h , and then the probability density function (PDF) can be expressed as
p h | α = i = 1 L N h i | 0 , α i 1
where α i is the precision of a Gaussian density function. A Gamma prior is considered over α
p α | a , b = i = 1 L Γ α i | a , b
Based on Bayesian inference, we can express the posterior for h as a Gaussian distribution with the mean and variance given as follows
μ = α 0 Σ Φ T h Σ = α 0 Φ T Φ + A 1
where A = d i a g α 1 , α 2 , , α L . By marginalizing over the h , the logarithm marginal likelihood for a and α 0 can be expressed as
L α , α 0 = log p y | h , α 0 p h | α d h = 1 2 K log 2 π + log C + y T C 1 y
where C = σ 2 I + Φ A 1 Φ T . By using the type-II ML [29] to maximize Equation (21), the parameters α and α 0 can be updated as follows
α i n e w = γ i μ i 2 1 α 0 n e w = y Φ μ 2 2 K i γ i
where μ i is the ith element of μ i and γ i = 1 α i Σ i i .
It is worth noting that the channel estimation in mmWave massive MIMO involves complex signal recovery so that the traditional Bayesian inference algorithm designed to recover the real signal in [30,31] needs to be improved to fit the complex signal. A feasible method is to slip the real and imaginary components as follows
R y I y = R Φ Ψ I Φ Ψ I Φ Ψ R Φ Ψ R h I h + R e N I e N
where R · denotes the real component and I · denotes the imaginary component. As a result, the dimensions will increase so a simplified model for Bayesian inference is essential and the main motivation of adopting the RVM is because of its highly efficient computation.
Based on above analysis, the Bayesian compressed sensing (BCS) algorithm can be used to recover the sparse channel through proper conversion. In addition, as the BCS algorithm highly depends on the statistical distribution of parameters and more available data at receiver can fit the prior distribution better, the sparseness and the number of pilot symbols are key factors that affect the recovery accuracy. Moreover, the BCS algorithm may break out early when the recovery accuracy reaches the threshold.
Remark 1: The overall computational complexities of three categories of algorithms are compared in Table 1. For convex relaxation such as IR-based algorithm, its computational complexity will grow squarely with the sparseness increase, since there are gradients of θ T and θ R that need to be calculated in each iteration [9]. For greedy iteration algorithms, its computational complexity will grow linearly with the sparseness due to the number of iterations depends on the sparseness [25]. With respect to Bayesian inference algorithms, its computational complexity will mainly depend on the amount of available pilot symbols at receiver where more available pilot symbols will let the algorithm break out early [30]. In Table 1, we can see that the computational complexity of all algorithms is quite high when deployed on massive MIMO antennas. The CS-based algorithms proposed in the previous works mainly aim at improving the estimation accuracy rather than reducing the computational complexity. However, reducing the computational complexity without significantly decreasing the accuracy is a direction that is worth studying.

5. Simulation Results

In this section, we evaluate the complexity and NMSE performance of three categories of algorithms through Monte Carlo simulation with the following parameters: (1) Both transmitter and receiver are equipped with ULAs where N T = N R = 16 and N R F R = N R F T = N R F = 4 . (2) The sin θ appearing in the array response vector is uniformly distributed in [−1,1) [25]. (3) The received SNR is defined as S N R n = 10 log 10 P / σ 2 d B . In addition, the hardware impairment at transmitter is defined as S N R e = 10 log 10 1 / k e d B , where the ratio coefficient k e reflects the degree of hardware impairment. Without loss of generality, we set k e = 0.05 , which is larger than that in the long-term evolution (LTE) standard as the signals in mmWave band are more vulnerable to hardware impairments [14]. In addition, we take the NMSE as the performance metric to evaluate the estimated quality, which is defined as
N M S E = Δ E H ^ H 2 2 / H 2 2
where H ^ denotes the estimated channel.

5.1. Comprehensive Comparison of Estimation Quality

In Figure 2, we compare the NMSE performance of different methods with respect to received SNR. We can observe that the performance of greedy iterative algorithms are best at the low SNR range. This is because the noise will be treated as the non-zero elements, which badly destroy the sparsity of the channel and reduce the amount of available pilot symbols. However, the gradient descent and iterative reweighting method are best at the high SNR range. It can also be observed that the performance of greedy iterative algorithms will remain unchanged with the SNR improving at the high SNR range because the estimated error caused by noise becomes smaller and smaller, but the quantization error decided by the resolution of angle quantization remains unchanged. Meanwhile, we can see the NMSE of Bayesian inference algorithm linearly decreases with the SNR improving, approximatively.
In Figure 3, we illustrate the NMSE performance versus the sparseness L. We can see that the NMSE of greedy iteration algorithms remains unchanged with the sparseness increasing, approximately, since the most relevant columns are selected in each iteration. The performance of IR algorithm also remains unchanged. Moreover, it is better than others at the low sparseness range. In addition, the performance of BCS and CVX algorithms improve as the sparseness increases. This is because the assumptions of prior distribution and the condition that let l 1 norm replace l 0 norm are fitted better when the sparseness increases. Furthermore, we can observe that the performance of BCS algorithm is better than greedy iterative algorithms when the sparseness is larger than 3 and the gap even reaches 8 dB when the sparseness is 7.
Figure 4 depicts pilot symbol frames versus the NMSE performance where each frame is assumed to be length of 16 pilot symbols. It is obvious that the NMSE becomes smaller with the increasing of pilot symbols for the three categories of algorithms. However, for the OMP and BCS algorithms, the NMSE curves become flat as the number of pilot symbol frames becomes large enough. The NMSE of channel estimated by convex relaxation algorithm such as basis pursuit (BP) algorithm improves with the pilot symbol frames increasing. The figure shows that the performance of BP and BCS algorithms even make the channel estimation results unusable at the low SNR range, but they are obviously better than OMP algorithm at the high SNR range. The figure also shows that the OMP algorithm is optimal in the range of low SNR with a few pilot symbols. After all, the estimation quality of OMP algorithm is least affected by the pilot symbols, while the BP algorithm is most affected.
In Figure 5, we illustrate the NMSE performance versus the angle quantization G. We can see that the performance of BP algorithm is hardly affected by angle quantization, which means it can avoid the quantization error as the AoAs and AoDs are accurately estimated by gradient descent method in each iteration. However, the estimation quality of BP algorithm is worse than OMP and BCS algorithms when the angle quantization is larger than 20. For the OMP algorithm, the MNSE decreases as angle quantization increases, which means the angle quantization is a key factor that affects the estimation performance of OMP algorithm. As to the BCS algorithm, the NMSE curve goes down rapidly and then stays flat, which indicates there is an optimal angle quantization G for it.
In Figure 6, we show the effect of transmitter’s hardware impairments on channel estimation quality, where the noise contained in the SNR is the AWGN without considering the angle quantization error and the hardware impairments. We can see that the estimation quality of all algorithms is affected by the transmitter’s hardware impairments, especially at the high SNR range. This is because the transmitter’s hardware impairment can be regarded as the extra noise and it is the main component of noise at the high SNR range. It is worth mentioning that the hardware impairments will increase with the increase of transmission power, which means simply increasing the transmission power cannot avoid the hardware impairments. Therefore, estimation and compensation are necessary for a practical mmWave massive MIMO system to combat them.
In Figure 7, we illustrate the success probability of channel estimation versus the number of antennas at the transmitter. Without loss of generality, we set 2 N R F T = N T , G T = 2 N T and γ t h = 10 dB, which indicates the number of RF chains N R F and the uniform grid of angle quantization G linearly increase when the antennas as transmitter sincrease, and the channel estimation is successful when the estimated error is less than −10 dB. We can observe that the success probability increases and then remains unchanged as the number of antennas increases for OMP algorithm because the angle quantization error at transmitter limits the success probability when the number of antennas is small, but the angle quantization error at receiver is the decisive factor when the antennas of transmitter is large enough. The same phenomenon can be observed in the BCS algorithm. However, different from the OMP algorithm, the number of available pilot symbols is the main factor that limits the success probability of BCS algorithm. As for BP algorithm, the number of antennas has little effect on the success probability, since the NMSE of BP algorithm is mainly affected by SNR at the receiver.
Based on above simulation results and analysis, a comprehensive comparison of different algorithms can be formed. Specifically, the convex relaxation algorithms are mainly affected by received SNR and transmitter’s hardware impairments, while the angle quantization error is the key factor that limits the performance of greedy iteration algorithms and all referred factors have an impact on Bayesian inference algorithms.

5.2. Computation Complexity versus Sparseness

Regarding to the computation complexity, a visual method is to compare the runtime of each estimator [11]. We performed simulations in MATLAB R2015b using a 3.0 GHz Inter Core i5-7400 CPU with 8 GB of memory and used 5000 Monte Carlo simulation trials to average the results. To ensure the fairness, we only treated the direct computational procedure as the runtime [11].
In Figure 8, we can observe that the runtime of IR algorithm increases squarely with the sparseness, approximately. At the low sparseness range, calculating the gradient has no significant contribution to runtime but it is the main source when the sparseness increases. The runtime of CVX algorithm remains unchanged but much longer than others as the CVX is a common toolbox and its runtime is not affected by the sparseness. It is worth mentioning that the runtime of BCS algorithm decreases with the sparseness increasing because the algorithm will break out ahead when the accuracy is satisfied. On the other hand, it can been observed that the runtime of greedy iterative algorithms is shortest and grows linearly with the sparseness increasing, since there are no complicated operations of the matrix. Among the greedy iterative algorithms, the runtime of OMP algorithm is longest as only one column is selected in each iteration.

6. Conclusions

In this paper, with a unified massive MIMO framework, we overview the CS-based mmWave channel estimation algorithms, namely convex relaxation algorithms, greedy iteration algorithms and Bayesian inference algorithms. The properties of three categories of algorithms are presented and how the estimation quality is affected by the sparseness, pilot symbols, angle quantization and transmitter’s hardware impairment is also compared. Acquiring a tradeoff between the performance and computation complexity, the greedy iterative algorithm SP outperforms others at the low SNR region but the BCS algorithm is the best at the high SNR range. Furthermore, it is possible to extend our work to other sparse signal recovery problems as the three categories of algorithms are common methods.

Author Contributions

X.L. and W.Y. conceived the model; X.L. performed the performance analysis and simulation results; W.Y. and Y.C. analyzed the numerical and simulation; X.L. wrote the paper; and W.Y. and X.G. provided some suggestions and revised the paper.

Funding

This work was supported by the National Natural Science Foundation of China (Nos. 61771487 and 61471393).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xiao, M.; Mumtaz, S.; Huang, Y.; Dai, L.; Li, Y. Millimeter Wave Communications for Future Mobile Networks. IEEE J. Sel. Areas Commun. 2017, 35, 1909–1935. [Google Scholar] [CrossRef] [Green Version]
  2. Yang, L.; Zeng, Y.; Zhang, R. Channel Estimation for Millimeter Wave MIMO Communications with Lens Antenna Arrays. IEEE Trans. Veh. Technol. 2018, 67, 3239–3251. [Google Scholar] [CrossRef]
  3. Gao, Z.; Dai, L.; Mi, D.; Wang, Z. Mmwave massive MIMO-based Wireless Backhaul for the 5G Ultra-Dense Network. IEEE Wirel. Commun. 2015, 22, 13–21. [Google Scholar] [CrossRef]
  4. Chen, J.; Yin, X.; Cai, X.; Wang, S. Measurement-based Massive MIMO Channel Modeling for Outdoor LoS and NLoS Environments. IEEE Access 2017, 5, 2126–2140. [Google Scholar] [CrossRef]
  5. Lin, X.; Wu, S.; Jiang, C.; Kuang, L.; Yan, J.; Hanzo, L. Estimation of Broadband Multiuser Millimeter Wave Massive MIMO-OFDM Channels by Exploiting Their Sparse Structure. IEEE Trans. Wirel. Commun. 2018, 17, 3959–3973. [Google Scholar] [CrossRef]
  6. Mezghani, A.; Swindlehurst, A. Blind Estimation of Sparse Broadband Massive MIOM Channels with Ideal and One-Bit ADCs. IEEE Trans. Signal Process. 2017, 66, 2972–2983. [Google Scholar] [CrossRef]
  7. Zhang, J.; Yuan, X.; Zhang, Y. Blind Signal Detection in Massive MIMO: Exploiting the Channel Sparsity. IEEE Trans. Commun. 2018, 66, 700–712. [Google Scholar] [CrossRef]
  8. Ma, J.; Ping, L. Data-aided Channel Estimation in Large Antenna Systems. IEEE Trans. Signal Process. 2014, 62, 3111–3124. [Google Scholar]
  9. Hu, C.; Dai, L.; Mir, T.; Gao, Z.; Fang, J. Super-resolution Channel Estimation for mmWave Massive MIMO with Hybrid Precoding. IEEE Trans. Veh. Technol. 2017, 67, 8954–8958. [Google Scholar] [CrossRef]
  10. Li, X.; Fang, J.; Li, H.; Wang, P. Millimeter Wave Channel Estimation via Exploiting Joint Sparse and Low-Rank Structures. IEEE Trans. Wirel. Commun. 2018, 17, 1123–1133. [Google Scholar] [CrossRef]
  11. Tsai, C.-R.; Liu, Y.-H.; Wu, A.-Y. Efficient Compressive Channel Estimation for Millimeter Wave Large-Scale Antenna Systems. IEEE Trans. Signal Process. 2018, 66, 2414–2428. [Google Scholar]
  12. Kiran, V.; Ahmed, A.; Nuria, G.P.; Robert, J.; Heath, W. Chnanel Estimation for Hybrid Architecture-based Wideband Millimeter Wave System. IEEE J. Sel. Areas Commun. 2017, 35, 1996–2009. [Google Scholar]
  13. Wen, C.K.; Jin, S.; Wong, K.K.; Chen, J.C.; Ting, P. Channel Estimation for Massive MIMO Using Gaussian Mixture Bayesian Learning. IEEE Trans. Wirel. Commun. 2015, 14, 1356–1368. [Google Scholar] [CrossRef]
  14. Wu, Y.; Gu, Y.; Wang, Z. Channel Estimation for mmWave MIMO with Transmitter Hardware Impairments. IEEE Commun. Lett. 2017, 22, 320–323. [Google Scholar] [CrossRef]
  15. Xiong, Y.; Wei, N.; Zhang, Z.; Li, B.; Chen, Y. Channel Estimation and IQ Imbalance Compensation for Uplink Massive MIMO Systems with Low-Resolution ADCs. IEEE Access 2017, 5, 6372–6388. [Google Scholar] [CrossRef]
  16. Wang, Y.F.; Lee, J.H. A Simple Phase Noise Suppression Scheme for Massive MIMO Uplink Systems. IEEE Trans. Veh. Technol. 2017, 66, 4769–4780. [Google Scholar] [CrossRef]
  17. Hojatian, H.; Omidi, M.J.; Saeedi-Sourck, H.; Farhang, A. Joint CFO and Channel Estimation in OFDM-based Massive MIMO Systems. In Proceedings of the International Symposium on Telecommunications, Chania, Greece, 4–6 October 2016; pp. 343–348. [Google Scholar]
  18. Bjrnson, E.; Hoydis, J.; Kountouris, M.; Debbah, M. Massive MIMO Systems with Non-Ideal Hardware: Energy Efficiency, Estimation, and Capacity Limits. IEEE Trans. Inf. Theory 2013, 60, 7112–7139. [Google Scholar] [CrossRef]
  19. Studer, C.; Wenk, M.; Burg, A. MIMO Transmission with Residual Transmit-RF Impairments. In Proceedings of the International ITG Workshop on Smart Antennas, Berlin, Germany, 16–18 February 2009; pp. 189–196. [Google Scholar]
  20. Bajwa, W.U.; Haupt, J.; Sayeed, A.M.; Nowak, R. Compressed Channel Sensing: A New Approach to Estimating Sparse Multipath Channels. Proc. IEEE 2010, 98, 1058–1076. [Google Scholar] [CrossRef]
  21. Heath, R.W.; Gonz lez-Prelcic, N.; Rangan, S.; Roh, W.; Sayeed, A.M. An Overview of Signal Processing Techniques for Millimeter Wave MIMO Systems. IEEE J. Sel. Top. Signal Process. 2017, 10, 436–453. [Google Scholar] [CrossRef]
  22. Zhang, J.; Dai, L.; Li, X.; Liu, Y.; Hanzo, L. On Low-Resolution ADCs in Practical 5G Millimeter Wave Massive MIMO Systems. IEEE Commun. Mag. 2018, 56, 2–8. [Google Scholar] [CrossRef]
  23. Zhang, Y.; Wang, D.; Wang, J.; You, X. Channel Estimation for Massive MIMO-OFDM Systems by Tracking the Joint Angle-Delay Subspace. IEEE Access 2016, 4, 10166–10179. [Google Scholar] [CrossRef]
  24. Alkhateeb, A.; Ayach, O.E.; Leus, G.; Heath, R.W. Channel Estimation and Hybrid Precoding for Millimeter Wave Cellular Systems. IEEE J. Sel. Top. Signal Process. 2014, 8, 831–846. [Google Scholar] [CrossRef] [Green Version]
  25. Junho, L.; Gye-Tae, G.; Yong, L. Chnanel Estimation via Orthogonal Matching Pursuit for Hybrid MIMO Systems in Millimeter Wave Communications. IEEE Trans. Commun. 2016, 64, 2370–2386. [Google Scholar]
  26. Ren, X.; Chen, W.; Tao, M. Position-based Compressed Channel Estimation and Pilot Design for High-Mobility OFDM Systems. IEEE Trans. Veh. Technol. 2015, 64, 1918–1929. [Google Scholar] [CrossRef]
  27. Nguyen, S.L.H.; Ghrayeb, A. Compressive Sensing-based Channel Estimation for Massive Multiuser MIMO Systems. In Proceedings of the IEEE Wireless Communications and Networking Conference, Shanghai, China, 7–10 April 2013; pp. 2890–2895. [Google Scholar]
  28. Ji, S.; Xue, Y.; Carin, L. Bayesian Compressive Sensing. IEEE Trans. Signal Process. 2008, 56, 2346–2356. [Google Scholar] [CrossRef]
  29. Chen, C.Y.; Wu, W.R. Joint AoD, AoA, and Channel Estimation for MIMO-OFDM Systems. IEEE Trans. Veh. Technol. 2018, 67, 5806–5820. [Google Scholar] [CrossRef]
  30. Carlin, M.; Rocca, P.; Oliveri, G.; Viani, F. Directions-of-Arrival Estimation Through Bayesian Compressive Sensing Strategies. IEEE Trans. Antennas Propag. 2013, 61, 3828–3838. [Google Scholar] [CrossRef]
  31. Rocca, P.; Hannan, M.A.; Salucci, M.; Massa, A. Single-Snapshot DoA Estimation in Array Antennas with Mutual Coupling through a Multi-Scaling Bayesian Compressive Sensing Strategy. IEEE Trans. Antennas Propag. 2017, 65, 3203–3213. [Google Scholar] [CrossRef]
Figure 1. Illustration of a hybrid MIMO architecture with the fully connected phase shifting network.
Figure 1. Illustration of a hybrid MIMO architecture with the fully connected phase shifting network.
Applsci 09 04346 g001
Figure 2. NMSE of the estimated CSI versus received SNR with sparseness L = 5 , pilot symbol frames M = 16 , angle quantization G = 32 and without hardware impairment.
Figure 2. NMSE of the estimated CSI versus received SNR with sparseness L = 5 , pilot symbol frames M = 16 , angle quantization G = 32 and without hardware impairment.
Applsci 09 04346 g002
Figure 3. NMSE of the estimated CSI versus sparseness with received S N R = 10 dB, pilot symbol frames M = 16 , angle quantization G = 32 and without hardware impairment.
Figure 3. NMSE of the estimated CSI versus sparseness with received S N R = 10 dB, pilot symbol frames M = 16 , angle quantization G = 32 and without hardware impairment.
Applsci 09 04346 g003
Figure 4. NMSE of the estimated CSI versus received pilot symbols with sparseness L = 5 , received S N R = 10 dB, angle quantization G = 32 and without hardware impairment.
Figure 4. NMSE of the estimated CSI versus received pilot symbols with sparseness L = 5 , received S N R = 10 dB, angle quantization G = 32 and without hardware impairment.
Applsci 09 04346 g004
Figure 5. NMSE of the estimated CSI versus received angle quantization (G) with received S N R = 10 dB, sparseness L = 5 , pilot symbol frames M = 40 and without hardware impairment.
Figure 5. NMSE of the estimated CSI versus received angle quantization (G) with received S N R = 10 dB, sparseness L = 5 , pilot symbol frames M = 40 and without hardware impairment.
Applsci 09 04346 g005
Figure 6. NMSE of the estimated CSI versus received SNR with hardware impairment = 13 dB, sparseness L = 5 , angle quantization G = 32 and pilot symbol frames M = 40 .
Figure 6. NMSE of the estimated CSI versus received SNR with hardware impairment = 13 dB, sparseness L = 5 , angle quantization G = 32 and pilot symbol frames M = 40 .
Applsci 09 04346 g006
Figure 7. Success probability of channel estimation versus the number of antennas ( N T ) with received S N R = 10 dB, N R = 4 , N R F R = 4 , γ t h = 10 dB, sparseness L = 5 and pilot symbol frames M = 16 and without hardware impairment.
Figure 7. Success probability of channel estimation versus the number of antennas ( N T ) with received S N R = 10 dB, N R = 4 , N R F R = 4 , γ t h = 10 dB, sparseness L = 5 and pilot symbol frames M = 16 and without hardware impairment.
Applsci 09 04346 g007
Figure 8. Runtime of the estimation algorithms versus the sparseness with received S N R = 10 dB, angle quantization G = 32 , pilot symbol frames M = 16 and without hardware impairment.
Figure 8. Runtime of the estimation algorithms versus the sparseness with received S N R = 10 dB, angle quantization G = 32 , pilot symbol frames M = 16 and without hardware impairment.
Applsci 09 04346 g008
Table 1. Computational complexity comparison of the three categories of algorithms.
Table 1. Computational complexity comparison of the three categories of algorithms.
AlgorithmsComputational Complexity
Convex Relaxation Algorithms [9] O ( N R F R N R F T ( N T + N R ) L 2 )
Greedy Iterative Algorithms [25] O ( N T N R L )
Bayesian Inference Algorithms [30] O ( N R l o g 2 N R )

Share and Cite

MDPI and ACS Style

Lu, X.; Yang, W.; Cai, Y.; Guan, X. Comparison of CS-Based Channel Estimation for Millimeter Wave Massive MIMO Systems. Appl. Sci. 2019, 9, 4346. https://doi.org/10.3390/app9204346

AMA Style

Lu X, Yang W, Cai Y, Guan X. Comparison of CS-Based Channel Estimation for Millimeter Wave Massive MIMO Systems. Applied Sciences. 2019; 9(20):4346. https://doi.org/10.3390/app9204346

Chicago/Turabian Style

Lu, Xingbo, Weiwei Yang, Yueming Cai, and Xinrong Guan. 2019. "Comparison of CS-Based Channel Estimation for Millimeter Wave Massive MIMO Systems" Applied Sciences 9, no. 20: 4346. https://doi.org/10.3390/app9204346

APA Style

Lu, X., Yang, W., Cai, Y., & Guan, X. (2019). Comparison of CS-Based Channel Estimation for Millimeter Wave Massive MIMO Systems. Applied Sciences, 9(20), 4346. https://doi.org/10.3390/app9204346

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

Article Metrics

Back to TopTop