## FCFGS-CV-Based Channel Estimation for Wideband MmWave Massive MIMO Systems with Low-Resolution ADCs

In-soo Kim and Junil Choi

Abstract -In this paper, the fully corrective forward greedy selection-cross validation-based (FCFGS-CV-based) channel estimator is proposed for wideband millimeter wave (mmWave) massive multiple-input multiple-output (MIMO) systems with lowresolution analog-to-digital converters (ADCs). The sparse nature of the mmWave virtual channel in the angular and delay domains is exploited to convert the maximum a posteriori (MAP) channel estimation problem to an optimization problem with a concave objective function and sparsity constraint. The FCFGS algorithm, which is the generalized orthogonal matching pursuit (OMP) algorithm, is used to solve the sparsity-constrained optimization problem. Furthermore, the CV technique is adopted to determine the proper termination condition by detecting overfitting when the sparsity level is unknown.

with low-resolution ADCs. The maximum a posteriori (MAP) channel estimation problem is formulated based on the on-grid compressed sensing framework. To promote the sparsity of the mmWave virtual channel in the angular and delay domains, high grid resolution is employed to reduce the off-grid error between continuous and discretized angles and delays. Then, the problem is converted to an optimization problem with a concave objective function and sparsity constraint. To solve this sparsity-constrained optimization problem, the FCFGS algorithm [6] is adopted, which is the generalized orthogonal matching pursuit (OMP) algorithm [7]. In addition, the CV technique [8] is used to detect overfitting, thereby properly terminating FCFGS when the sparsity level is unknown.

## I. INTRODUCTION

The millimeter wave (mmWave) band in the range of 30 -300 GHz offers wide bandwidth, thereby enabling the use of high data rates [1]. Furthermore, a large number of antennas can be compactly deployed due to the small wavelength, which is known as massive multiple-input multiple-output (MIMO). One problem of mmWave massive MIMO, however, is the prohibitive analog-to-digital converter (ADC) power consumption incurred by the high sampling frequency and large number of radio frequency (RF) chains. To deal with this issue, one possible solution is to use low-resolution ADCs because the ADC power consumption increases exponentially with the ADC resolution [2].

In [3]-[5], compressed sensing-based algorithms for MIMO systems with low-resolution ADCs were proposed. The generalized approximate message passing (GAMP) and vector AMP (VAMP) algorithms [3] are low complexity belief propagationbased (BP-based) algorithms. However, GAMP and VAMP break down when the sensing matrix is ill-conditioned. The turbo principle-based generalized expectation consistent signal recovery (GEC-SR) [4] and hierarchical model-based sparse Bayesian learning (SBL) [5] algorithms are relatively robust to the condition number of the sensing matrix but require matrix inversion, which results in high complexity.

In this paper, we propose the fully corrective forward greedy selection-cross validation-based (FCFGS-CV-based) channel estimator for wideband mmWave massive MIMO systems

The authors are with the Department of Electrical Engineering, POSTECH, Pohang, Korea (e-mail: insookim@postech.ac.kr; junil@postech.ac.kr).

This work was supported by the Institute for Information &amp; communications Technology Promotion(IITP) grant funded by the Korea government(MSIT) (No.2016-0-00123, Development of Integer-Forcing MIMO Transceivers for 5G &amp; Beyond Mobile Communication Systems).

Notation: a , a , and A denote a scalar, vector, and matrix. ‖ a ‖ 0 and ‖ a ‖ represent the 0 -norm and 2 -norm of a . ‖ A ‖ F is the Frobenius norm of A . The Kronecker product of A and B is denoted as A ⊗ B . The support of a is represented as supp( a ) . The row restriction of A to the index set I is A I . [ n ] denotes [ n ] = { 1 , · · · , n } . C ([ n , k ] ) represents the set of all k -combinations of [ n ] .

## II. SYSTEM MODEL

Assume a single-cell uplink massive MIMO system with an M -antenna base station and K single-antenna users. 1 A pair of B -bit ADCs is equipped at each receive antenna to reduce the power consumed by the base station. The system operates in the mmWave wideband with D delay taps.

The received signal y [ n ] ∈ C M at time n is

$$r \, \text{AMP} \quad \text{$line receive} \, \text{signal} \, y [ n ] \in \mathbb { C } \dashv \text{ at time $n$ is} \\ \text{gation} _ { \quad } \\ \text{d. The} \quad \text{$y[n]=\underbrace { [H[0] \ \cdots \ \mathbb { H } [ D - 1 ] ] } _ { = \mathbb { H } } \underbrace { \left [ \substack { s[n] \\ \vdots \\ s [ n - D + 1 ] \right ] } _ { = s _ { n } } } _ { = s _ { n } } + \mathbf v [ n ] \ ( 1 ) \\ \text{matrix} \quad \text{where} \, \mathbf H [ d ] \in \mathbb { C } ^ { M \times K } \text{ is the $d$-th delay tap channel, $s[n] \in \mathbb { C } ^ { K } \\ \quad \text{$is the training} \, signal at time $n$ satisfying $\mathbb{E} \{s[n]s[n]^{H}]$} \right ] =$$

where H [ d ] ∈ C M K × is the d -th delay tap channel, s [ n ] ∈ C K is the training signal at time n satisfying E { s [ n ] s [ n ] H } = ρ I K with ρ representing the signal-to-noise ratio (SNR), and v [ n ] ∼ CN ( 0 M , I M ) is the additive white Gaussian noise (AWGN) at time n .

1 To compensate the severe path loss in the mmWave band, the users would deploy multiple antennas as well. In practice, however, the limited hardware cost would allow the users to deploy only a single RF chain. Assuming that the users perform proper analog beamforming, the base station would effectively see single-antenna users.

The channels in the mmWave band are characterized by a small number of paths. The k -th column of H [ d ] , which is the channel of the k -th user, is [9]

$$\mathbf h _ { k } [ d ] = \sum _ { \ell = 1 } ^ { L _ { k } } \alpha _ { k, \ell } p _ { k } ( d T - \tau _ { k, \ell } ) \mathbf a ( \theta _ { k, \ell } ) \quad ( 2 ) \quad \mathbf x = \mathbf v \quad \text{estim} \\ \text{ere $L$} \text{$L$} \text{ is the number of $maths_{\ } \alpha_{\nu,\varepsilon} \in \mathbb{C}$ is the $\ell$-th}$ math} \quad \text{Re}$$

where L k is the number of paths, α k,/lscript ∈ C is the /lscript -th path gain, θ k,/lscript ∈ [ -π/ , π/ 2 2] is the /lscript -th angle-of-arrival (AoA), τ k,/lscript ∈ [0 , ( D -1) T ] is the /lscript -th delay, T is the sampling period, a ( θ ) ∈ C M is the array response vector, and p k ( ) t is the pulse shaping filter normalized to satisfy E {‖ h k [ d ] ‖ } 2 = M .

In the channel estimation phase of length N , the received signal Y is

$$\text{signal $Y$ is} \\ = & H \underbrace { [ s _ { 0 } \cdots s _ { N - 1 } ] } _ { = S } + \underbrace { [ v [ 0 ] \cdots v [ N - 1 ] ] } _ { = V }. \quad ( 3 ) } _ { = V }. \\ \text{At each receive antenna, the real and imaginary parts of $Y$ are} \\ \cdot \cdot \hat { \hat { \ } } \hat { \ } \hat { \ } \cdot.$$

At each receive antenna, the real and imaginary parts of Y are quantized by a pair of B -bit ADCs. The quantized received signal ˆ Y is

$$\hat { Y } = Q ( \text{HS} + V ) \text{ \quad \ \ } ( 4 ) \text{ \ \ word}$$

where Q( ) · is the B -bit quantization function defined as

$$\hat { y } = Q ( y ) \iff \begin{cases} \Re e ( \hat { y } ^ { \text{lo} } ) \leq \Re e ( y ) < \Re e ( \hat { y } ^ { \text{up} } ) \\ \text{Im} ( \hat { y } ^ { \text{lo} } ) \leq \text{Im} ( y ) < \text{Im} ( \hat { y } ^ { \text{up} } ) \end{cases}$$

with ˆ y lo and ˆ y up the lower and upper thresholds associated with ˆ y . In other words, the real and imaginary parts of ˆ y lo , ˆ y up , and ˆ y correspond to one of the 2 B quantization intervals.

To promote the sparsity in the angular and delay domains, the virtual channel representation of H is employed. To transform h k [ d ] , define the AoA dictionary B ∈ C M R × AoA and delay dictionary p k [ d ] ∈ R R delay as

$$\mathbf B & = \left [ \mathbf a ( \hat { \theta } _ { 1 } ) \ \cdots \ \mathbf a ( \hat { \theta } _ { R _ { \Lambda \circ \Lambda } } ) \right ], \quad & ( 6 ) \\ \mathbf p _ { k } [ d ] & = \left [ p _ { k } ( d T - \hat { \tau } _ { 1 } ) \ \cdots \ \ p _ { k } ( d T - \hat { \tau } _ { R _ { \text{delay} } } ) \right ] ^ { \mathbf T } \quad & ( 7 ) \quad & \text{as in } (.$$

with R AoA ≥ M discretized AoAs and R delay ≥ D discretized delays. Then, the relationship between h k [ d ] and its virtual channel X k ∈ C R AoA × R delay is

$$\Phi _ { k } [ d ] & = \left [ p _ { k } ( d T - \hat { \tau } _ { 1 } ) \quad \cdots \quad p _ { k } ( d T - \hat { \tau } _ { R _ { \text{delay} } } ) \right ] ^ { \mathbf T } \quad ( 7 ) \quad \text{ as in } ( 5 \\ \text{with } R _ { \text{AoA} } \geq \, M \text{ discretized } \, \text{AoAs and } R _ { \text{delay} } \, \geq \, D \text{ dis- } \quad \text{used in} \,.$$

$$\mathbf h _ { k } [ d ] = \mathbf B \mathbf X _ { k } \mathbf p _ { k } [ d ], \quad \ \ ( 8 ) _ { \quad \ h e c a }$$

which means that and

$$\text{which means that} \\ \text{H} [ d ] = B \underbrace { [ X _ { 1 } \ \cdots \ X _ { K } ] } _ { = X } \underbrace { \left [ \begin{matrix} \mathbf P _ { 1 } [ d ] \\ & \ddots \\ & \mathbf P _ { K } [ d ] \end{matrix} \right ] } _ { = \mathbf P _ { [ d ] } } \begin{matrix} \text{in Secti} \\ & \text{from (1} \\ & \ell ( x ) = } \\ & = \\ \text{and} \end{matrix} } \text{$H = BX \left [ \mathbf P [ 0 ] \ \cdots \ \mathbf P [ D - 1 ] \right ]. \quad ( 1 0 ) \quad \text{where} \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d    } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { \ } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak {d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { D } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { d } \ \mathfrak { | } \ \mathfrak { d } \ \mathfrak { | } \ \mathfrak { d } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { \ } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak {| } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { { \ } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathfrak { | } \ \mathf ]$$

In the sequel, the shorthand notations L = ∑ K k =1 L k and R = R AoA R delay K will be used.

$$H = B X \underbrace { \left [ \mathbf P [ 0 ] \ \cdots \ \mathbf P [ D - 1 ] \right ] } _ { = \mathbf P }. \quad \text{(10)} \quad \text{where} \quad \text{of} \ \mathcal { N } ( \\ e \, \text{sequel, the shorthand notations } L = \sum _ { k = 1 } ^ { K } L _ { k } \text{ and } R = \quad \text{$^{2}_{\text{In}g}$} \\. \, R \,. \,. \, K \text{ will be needed}$$

To facilitate the analysis, (4) is vectorized in conjunction with (10) as

$$\hat { y } = Q ( \mathbf A x + \mathbf v )$$

where ˆ = vec( ˆ ) y Y ∈ C MN , A = S P T T ⊗ B ∈ C MN R × , x = vec( X ) ∈ C R , and v = vec( V ) ∈ C MN . The goal is to estimate x from ˆ y .

Remark 1: The fact which should be emphasized from (10) is that an L -sparse X exists 2 as R AoA →∞ and R delay →∞ with α k,/lscript as its elements [10]. In other words, the off-grid error is negligible when the grid resolution is high.

## III. PROPOSED APPROXIMATE MAP CHANNEL ESTIMATION

## A. FCFGS-Based Channel Estimation

In this paper, the parameters of H are assumed to be independent and identically distributed (i.i.d.) as

$$\alpha _ { k, \ell } \sim \mathcal { C N } ( 0, 1 ),$$

$$\theta _ { k, \ell } \sim \text{Uniform} ( [ - \pi / 2, \pi / 2 ] ), \text{ \quad \ \ } ( 1 3 )$$

$$\tau _ { k, \ell } \sim \text{Uniform} ( [ 0, ( D - 1 ) T ] ) \quad \quad ( 1 4 )$$

for all ( k, /lscript ) . To build the MAP channel estimation framework, the likelihood function is formulated based on the real counterparts of ˆ y , A , and x with the subscript R , which are

$$\hat { y } _ { \text{R} } & = \left [ \text{Re} ( \hat { y } ) ^ { \text{T} } \ \text{Im} ( \hat { y } ) ^ { \text{T} } \right ] ^ { \text{T} } \\ & = \begin{bmatrix} \hat { y } _ { \text{R}, 1 } & \cdots & \hat { y } _ { \text{R}, 2 M N } \end{bmatrix} ^ { \text{T} },$$

with the lower and upper thresholds associated with ˆ y R defined as

$$\hat { y } _ { R } & = \left [ \L Re ( \hat { y } ) ^ { T } \, \text{ Im} ( \hat { y } ) ^ { T } \right ] ^ { T } \\ & = \left [ \hat { y } _ { R, 1 } \, \cdots \, \hat { y } _ { R, 2 M N } \right ] ^ { T }, \quad \text{ (15)} \\ \text{ A} _ { R } & = \begin{bmatrix} \L Re ( \text{A} ) & - \text{Im} ( \text{A} ) \\ \text{Im} ( \text{A} ) & \text{Re} ( \text{A} ) \end{bmatrix} \\ & = \left [ \text{a} _ { R, 1 } & \cdots \, \text{ a} _ { R, 2 M N } \right ] ^ { T }, \quad \text{ (16)} \\ \text{ x} _ { R } & = \left [ \text{Re} ( \text{x} ) ^ { T } \, \text{ Im} ( \text{x} ) ^ { T } \right ] ^ { T } \quad \text{ (17)} \\ \text{the lower and upper thresholds associated with } \hat { y } _ { R } \\ \text{d as} \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$\hat { y } _ { R } ^ { \text{lo} } & = \left [ \hat { y } _ { R, 1 } ^ { \text{lo} } \quad \cdots \quad \hat { y } _ { R, 2 M N } ^ { \text{lo} } \right ] ^ { \text{T} }, \quad \quad ( 1 8 ) \\ \hat { y } _ { R } ^ { \text{up} } & = \left [ \hat { y } _ { R, 1 } ^ { \text{up} } \quad \cdots \quad \hat { y } _ { R, 2 M N } ^ { \text{up} } \right ] ^ { \text{T} } \quad \quad ( 1 9 )$$

as in (5). From now on, the complex and real counterparts are used interchangeably. To proceed with the likelihood function, note that conditioned on x ,

$$\lambda \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha,
 \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi,
 \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda. \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \lambda, \eta, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha,\alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha. \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \alpha, \beta, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi,\xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi. \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \x, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \Xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi, \xi,$$

$$\text{Ax} + \mathbf v \sim \mathcal { C N } ( \text{Ax}, \mathbf I _ { M N } ) \text{ \quad \ \ } ( 2 0 )$$

because v ∼ CN ( 0 MN , I MN ) is independent of x as assumed in Section II. Then, the likelihood function /lscript ( x ) is formulated from (11) and (20) as [11]

$$\Big | _ { \quad ( 9 ) } \quad & \ell ( x ) = \Pr \left [ \hat { y } | x \right ] \\ = \prod _ { i = 1 } ^ { 2 M N } \left ( \Phi \left ( \frac { \hat { y } _ { R, i } ^ { u p } - a _ { R, i } ^ { T } x _ { R } } { \sqrt { 1 / 2 } } \right ) - \Phi \left ( \frac { \hat { y } _ { R, i } ^ { l o } - a _ { R, i } ^ { T } x _ { R } } { \sqrt { 1 / 2 } } \right ) \right ) \\ ( 1 0 ) \quad & \text{where } \Phi ( \cdot ) \text{ denotes the cumulative distribution function (CDF)}$$

where Φ( ) · denotes the cumulative distribution function (CDF) of N (0 , 1) .

2 In general, X is not unique as long as R AoA &gt; M and R delay &gt; D . One possible X is X ∗ = B H ( BB H ) -1 H P ( T P ) -1 P T , which is the minimum 2 -norm solution.

To approximate the distribution of x based on Remark 1 , consider the discrete analogue of (13) and (14). In other words, consider an L -sparse vector with CN ( 0 L , I L ) as its elements and uniformly distributed discretized AoAs and delays, which is distributed as

$$x _ { \text{supp} ( \mathbf x ) } \sim \mathcal { C N } ( 0 _ { L }, \mathbf I _ { L } ), \quad \ \ ( 2 2 ) \quad \ \ _ { 4 } ^ { 3 \colon }.$$

$$\text{supp} ( \mathbf x ) \sim \text{Uniform} ( C ( [ R ], L ) ). \quad \text{(23)}$$

According to Remark 1 , the deviation of (22) and (23) from the distribution of x tightens as R AoA and R delay increase with the decreasing off-grid error. Therefore, the MAP channel estimation framework is proposed based on (22) and (23). The grid resolution is configured as R AoA /greatermuch M and R delay /greatermuch D to reduce the off-grid error, which is a widely adopted criterion [12]. However, the off-grid error cannot be eliminated unless R AoA →∞ and R delay →∞ . In short, the mismatch between (22), (23), and the distribution of x can be reduced as the grid resolution is increased but not eliminated.

Now, combining (21), (22), and (23), the MAP estimator ˆ x of x is formulated as 3

$$\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \\ \hat { x } & = \underset { x \in \mathbb { C } ^ { R } } { \arg \max } ( \ell ( x ) \cdot e ^ { - \| x \| ^ { 2 } } ) \ s. t. \ \| x \| _ { 0 } \leq L \quad \text{ (24)} \quad \text{verfitt}$$

$$\overset { ( a ) } { = } \underset { x \in \mathbb { C } ^ { R } } { \arg \max } ( \ell ( { \mathbf x } ) \cdot e ^ { - \| { \mathbf x } \| ^ { 2 } } ) \text{ s.t. } \| { \mathbf x } \| _ { 0 } \leq L$$

where (a) follows from ‖ x supp( x ) ‖ = ‖ x ‖ . From (24), observe that the objective function is log-concave because /lscript ( x ) has the form of Φ( b -x ) -Φ( a -x ) with b &gt; a , which is logconcave, and e -‖ ‖ x 2 has the minus of 2 -norm as its exponent, which is concave [11]. At this point, define the logarithm of the objective function as f ( x ) with its gradient ∇ f ( x ) , which are

$$f ( { \mathbf x } ) = \log \ell ( { \mathbf x } ) - \| { \mathbf x } \| ^ { 2 }, \text{ \quad \ \ } ( 2 5 )$$

$$\nabla f ( { \mathbf x } ) = \nabla \log \ell ( { \mathbf x } ) - 2 { \mathbf x }. \quad \quad ( 2 6 ) \quad \text{In ad}$$

Then, the goal is to solve

$$\hat { x } = \underset { x \in \mathbb { C } ^ { R } } { \arg \max } f ( x ) \ s. t. \ \| x \| _ { 0 } \leq L, \quad \quad ( 2 7 ) \ \text{ as}$$

an optimization problem with a concave objective function and sparsity constraint.

In general, (27) is NP-hard to solve because of its sparsity constraint. To approximately solve (27) iteratively, the FCFGS algorithm [6] is adopted, which generalizes the OMP algorithm [7] to convex objective functions as presented in Algorithm 1. Further explanation of the algorithm is provided in the next subsection but in short, Line 4 selects the largest element of the gradient of the objective function to update the support of the estimate in Line 5. Then, Line 6 performs convex optimization to update the estimate, whose convergence is guaranteed because the objective function and support constraint are convex.

Before moving on, we mention that the off-grid error incurs the leakage effect [9], which leads to an approximately (instead of exactly) sparse x . As a result, the recovery guarantees of FCFGS established from [6] break down, which may degrade the accuracy of FCFGS.

3 The constants independent of x are neglected for simplicity, which refer to the constants in the Gaussian probability density function (PDF) and probability mass function (PMF) of Uniform( C ([ R , L ] )) .

## Algorithm 1 FCFGS-CV algorithm to solve (27)

$$\varprojlim _ { \alpha \alpha } 6 \colon \quad b = \arg \max _ { - \mu } f _ { \mathcal { E } } ( \mathbf x ) \ s. t. \, \text{supp} ( \mathbf x ) \subseteq \mathcal { I }$$

$$1, \overline { \text{Algorithm} } 1 \, \text{FCFG-CV algorithm} \\ \text{is}, \, \text{$f_{\mathcal{E}(\ell),\,f_{\mathcal{C}\nu(\ell)$} }$} \\ \text{Output} \, \text{$\hat{x}$} \\ \, 1 \, b = 0 _ { R }, \, \epsilon = - \infty \\ \, 2 \, \text{while} \, \text{$f_{\mathcal{C}\nu(\mathcal{B})\,>\epsilon$} \, \mathrm do} \\ \, 3 \, \text{$\hat{x}=b$} \\ \, 4 \, \text{$\hat{i}=argmax|\nabla f_{\mathcal{E}(\hat{x})|} \\ \, \quad \, \, \,$$

## B. CV-Based Termination Condition

In practice, the knowledge of ‖ x ‖ 0 is critical to determining the proper termination condition of FCFGS. However, ‖ x ‖ 0 is difficult to acquire because L is often unknown. Furthermore, ‖ x ‖ 0 may deviate from L when the grid resolution is low. To determine the proper termination condition of FCFGS, the CV technique [8] is adopted, which is a model validation technique designed to assess the estimation quality, thereby preventing overfitting.

FCFGS-CV runs estimation using estimation data, whereas validation is performed based on CV data. The disjoint nature of estimation and CV data enables validation to properly assess the quality of the estimation data-based estimate. To proceed, ˆ y is divided to the estimation measurement ˆ y E ∈ C |E| and CV measurement ˆ y CV ∈ C |CV| where the disjoint sets E and CV partition [ MN ] . The estimation and CV sensing matrices are

$$A _ { \mathcal { E } } \in \mathbb { C } ^ { | \mathcal { E } | \times R } \ \ \colon \text{ estimation sensing matrix,} \quad ( 2 8 )$$

$$\mathsf \rangle \quad \mathsf A _ { \mathcal { C V } } \in \mathbb { C } ^ { | \mathcal { C V } | \times R } \colon \, \mathsf C \mathsf V \, \text{sensing matrix}. \quad \quad \ \ ( 2 9 )$$

In addition, define the estimation objective function f E ( x ) and CV objective function f CV ( x ) using ˆ y E , A E , ˆ y CV , and A CV as

$$f _ { \mathcal { E } } ( \mathbf x ) & = \log \Pr \left [ \hat { y } _ { \mathcal { E } } | \mathbf x \right ] - \| \mathbf x \| ^ { 2 }, \quad \quad ( 3 0 ) \\ f _ { \mathcal { C } \nu } ( \mathbf x ) & = \log \Pr \left [ \hat { y } _ { \mathcal { C } \nu } | \mathbf x \right ] - \| \mathbf x \| ^ { 2 }, \quad \quad ( 3 1 ) \\ \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots$$

whose likelihood functions are formulated as in (21). Then, ˆ x is updated based on f E ( ) · , whereas the quality of ˆ x is assessed using f CV ( ) · . The proposed FCFGS-CV algorithm is presented in Algorithm 1. In short, Line 2 tests whether f CV ( ) · increases or not to detect overfitting, which is widely used as an indicator of termination. To demonstrate how CV indicates the proper termination timing, define the normalized mean squared error (NMSE) as

$$J \varepsilon \vee \varepsilon ^ { + } \lfloor \varepsilon \lfloor \varepsilon \rfloor \right | \left \| \varepsilon \right \| \right | \,, \quad \vee \vee \rfloor \\ f _ { \mathcal { C } \nu } ( x ) & = \log \Pr \left [ \hat { y } _ { \mathcal { C } \nu } | x \right ] - \left \| x \right \| ^ { 2 }, \quad \left ( 3 1 \right ) \\ e \, \text{ likelihood functions are formulated as in } ( 2 1 ). \, \text{Then,} \, \hat { x } _ { \mathcal { C } \nu } \right \| \, \text{the, and, } \, \hat { e } \, \hat { e } \, \text{ is accessed}$$

$$\text{as} \quad \text{NMSE} = \mathbb { E } \left \{ \frac { \| \hat { H } - \mathbf H \| _ { F } ^ { 2 } } { \| \mathbf H \| _ { F } ^ { 2 } } \right \} \quad \text{(3)}$$

where ˆ H = BXP ˆ , which measures the accuracy of ˆ x . Then, an illustration of how the NMSE, f CV (ˆ) x , and f E (ˆ) x evolve with the iteration of FCFGS-CV for one problem instance from the simulation results in Section IV is shown in Fig. 1. According to Fig. 1, the decreasing point of f CV ( ) · clearly indicates when the minimum NMSE occurs, whereas f E ( ) · provides no clear indication of overfitting.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-84183337-e995-4d9b-b03f-d147fb389d36.jpg)

Fig. 1. NMSE versus iteration of FCFGS-CV with normalized f CV (ˆ) x and x Fig. 3. NMSE versus N for B = 1 2 3 4 , , , with SNR = -10 dB.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-cf9e42f6-75ba-48d6-be87-04517039566b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-58c2a570-b109-4ac6-8a7e-d7752622f3c0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-0691e8ef-7645-4e20-a657-4eb5f3a7eec1.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-b5bb6212-0d0e-4268-9553-f4f15258eb4a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-df9d643b-434e-4395-945d-c65ac4223eac.jpg)

Fig. 2. NMSE versus SNR for B = 1 2 3 4 , , , with N = 160 . In each subfigure, B = ∞ is shown as a reference. Fig. 4. NMSE versus N for B = 1 2 3 4 , , , with SNR = 0 dB.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-e9430ac5-c159-4a23-9af6-401f0321b493.jpg)

f E (ˆ) for B = 2 and SNR = 0 dB.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-40f51791-3b08-4a88-b59c-522d3b2d45b9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-770131ff-6533-4883-8234-c8c0b2a2d38c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-116524f1-6262-4952-aba7-05d352f8d4e3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-8c2ce13e-12f0-41ec-acd0-fd0956144e86.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-1d3336fe-c517-471a-b2fe-df9300a51e3d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00532v1-figure-31be5141-384c-4e17-b985-4941540141af.jpg)

## IV. SIMULATION RESULTS

In this section, we evaluate the performance of the proposed FCFGS-CV-based channel estimator based on its accuracy and complexity. The NMSE is the measure of accuracy. The system parameters are selected as M = 64 and K = 4 , while N varies from simulation to simulation. The training signals are configured as circularly shifted Zadoff-Chu (ZC) sequences of length N with cyclic prefixes of length D -1 . The base station employs a uniform linear array (ULA) structure with half-wavelength inter-element spacings, while the users adopt the raised-cosine (RC) pulse shaping filter with a roll-off factor of 0 35 . . A pair of B -bit uniform quantizers is deployed for the real and imaginary parts at each receive antenna. The channel parameters are D = 8 with L k = 2 for all k . The dictionaries are configured as R AoA = 2 M and R delay = 2 D to promote the sparsity of x . To implement CV, among N training signals, ( N -KD ) signals are used for estimation, while the remaining KD signals are used for CV.

In the simulation results, the GAMP [3], generalized VAMP (GVAMP) [3], GEC-SR [4], and generalized SBL (Gr-SBL) [5] algorithms for generalized linear models (GLMs) [13], [14] are adopted as benchmarks, which are state-of-the-art compressed sensing-based algorithms. Gr-SBL is configured as R AoA = 2 M and R delay = 2 D . In contrast, GAMP, GVAMP, and GEC-SR are selected as R AoA = M and R delay = D because these algorithms diverge when R AoA /greatermuch M and R delay /greatermuch D with ill-conditioned sensing matrices.

In Fig. 2, the NMSEs of FCFGS-CV and other benchmarks are provided for various SNRs with B = 1 2 3 4 , , , and fixed N = 160 . FCFGS-CV outperforms other benchmarks for all SNRs and B because its high grid resolution guarantees the validity of the assumption established from (22) and (23). In contrast, other benchmarks suffer from the mismatch between the postulated distributions (Gaussian mixture and Student-t) and distribution of x . Before moving on, we point out that the NMSE increases in the high SNR regime when B = 1 for all algorithms because the coarse quantization of y eliminates the

magnitude information of x .

In Figs. 3 and 4, the NMSEs are analyzed for various N with B = 1 2 3 4 , , , and fixed SNR = -10 0 , dB. Similar to Fig. 2, FCFGS-CV shows the best performance among all algorithms for all N and B , which reconfirms the superior performance of FCFGS-CV over other benchmarks. Therefore, FCFGS-CV is accurate for not only large N (= 160) but moderate N (= 80) as well.

Now, the asymptotic complexity of FCFGS-CV is evaluated. Assuming that R = O ( MDK ) with t -iteration gradient descent method used in Line 6 of Algorithm 1, the complexities of the k -th iteration of FCFGS-CV, GAMP, GVAMP, GEC-SR, and Gr-SBL are O ( M NDK 2 + MNkt ) , O ( M NDK 2 ) [9], O ( M NDK 2 ) [9], O ( M ND K 3 2 2 ) [4], and O ( M ND K 3 2 2 ) [15]. The complexities of FCFGS-CV, GAMP, and GVAMP have the same order, whereas the complexities of GEC-SR and Gr-SBL are increased by a factor of O ( MDK ) compared to those of FCFGS-CV, GAMP, and GVAMP. Therefore, FCFGSCV is as efficient as GAMP and GVAMP, which is appealing because AMP-based algorithms are considered as practical channel estimators for massive MIMO systems in the research perspective [9]. In practice, the fact that the complexity of FCFGS-CV increases with M as O ( M 2 ) may be problematic, so reducing its complexity is left as an interesting future work.

## V. CONCLUSION

In this paper, we proposed the FCFGS-CV-based channel estimator for wideband mmWave massive MIMO systems with low-resolution ADCs. The MAP channel estimation framework was formulated as a sparsity-constrained optimization problem with a concave objective function. Then, the FCFGS algorithm was adopted in conjunction with the CV technique, which was revealed to be accurate and efficient in the simulation results.

## REFERENCES

- [1] A. L. Swindlehurst, E. Ayanoglu, P. Heydari, and F. Capolino, 'Millimeter-wave massive MIMO: the next wireless revolution?' IEEE Commun. Mag. , vol. 52, no. 9, pp. 56-62, Sep. 2014.
- [2] R. H. Walden, 'Analog-to-digital converter survey and analysis,' IEEE J. Sel. Areas Commun. , vol. 17, no. 4, pp. 539-550, Apr. 1999.
- [3] X. Meng, S. Wu, and J. Zhu, 'A unified bayesian inference framework for generalized linear models,' IEEE Signal Process. Lett. , vol. 25, no. 3, pp. 398-402, Mar. 2018.
- [4] H. He, C. Wen, and S. Jin, 'Bayesian optimal data detector for hybrid mmWave MIMO-OFDM systems with low-resolution ADCs,' IEEE J. Sel. Topics Signal Process. , vol. 12, no. 3, pp. 469-483, June 2018.
- [5] X. Meng and J. Zhu, 'A generalized sparse bayesian learning algorithm for 1-bit DOA estimation,' IEEE Commun. Lett. , vol. 22, no. 7, pp. 1414-1417, July 2018.
- [6] S. Shalev-Shwartz, N. Srebro, and T. Zhang, 'Trading accuracy for sparsity in optimization problems with sparsity constraints,' SIAM Journal on Optimization , vol. 20, no. 6, pp. 2807-2832, 2010.
- [7] J. A. Tropp and A. C. Gilbert, 'Signal recovery from random measurements via orthogonal matching pursuit,' IEEE Trans. Inf. Theory , vol. 53, no. 12, pp. 4655-4666, Dec. 2007.
- [8] P. Boufounos, M. F. Duarte, and R. G. Baraniuk, 'Sparse signal reconstruction from noisy compressive measurements using cross validation,' in 2007 IEEE/SP 14th Workshop on Statistical Signal Processing , Aug. 2007, pp. 299-303.
- [9] J. Mo, P. Schniter, and R. W. Heath, 'Channel estimation in broadband millimeter wave MIMO systems with few-bit ADCs,' IEEE Trans. Signal Process. , vol. 66, no. 5, pp. 1141-1154, Mar. 2018.
- [10] A. M. Sayeed, 'Deconstructing multiantenna fading channels,' IEEE Trans. Signal Process. , vol. 50, no. 10, pp. 2563-2579, Oct. 2002.
- [11] A. Mezghani, F. Antreich, and J. A. Nossek, 'Multiple parameter estimation with quantized channel output,' in 2010 International ITG Workshop on Smart Antennas (WSA) , Feb. 2010, pp. 143-150.
- [12] J. Lee, G. Gil, and Y. H. Lee, 'Channel estimation via orthogonal matching pursuit for hybrid MIMO systems in millimeter wave communications,' IEEE Trans. Commun. , vol. 64, no. 6, pp. 2370-2386, June 2016.
- [13] S. Rangan, 'Generalized approximate message passing for estimation with random linear mixing,' in 2011 IEEE International Symposium on Information Theory Proceedings , July 2011, pp. 2168-2172.
- [14] P. Schniter, S. Rangan, and A. K. Fletcher, 'Vector approximate message passing for the generalized linear model,' in 2016 50th Asilomar Conference on Signals, Systems and Computers , Nov. 2016, pp. 15251529.
- [15] Y. Ding, S. Chiu, and B. D. Rao, 'Bayesian channel estimation algorithms for massive MIMO systems with hybrid analog-digital processing and low-resolution ADCs,' IEEE J. Sel. Topics Signal Process. , vol. 12, no. 3, pp. 499-513, June 2018.
In-soo Kim Junil Choi
FCFGS-CV-Based Channel Estimation for Wideband MmWave Massive MIMO Systems with Low-Resolution ADCs
In this paper, the fully corrective forward greedy selection-cross
validation-based (FCFGS-CV-based) channel estimator is proposed for wideband
millimeter wave (mmWave) massive multiple-input multiple-output (MIMO) systems
with low-resolution analog-to-digital converters (ADCs). The sparse nature of
the mmWave virtual channel in the angular and delay domains is exploited to
convert the maximum a posteriori (MAP) channel estimation problem to an
optimization problem with a concave objective function and sparsity constraint.
The FCFGS algorithm, which is the generalized orthogonal matching pursuit (OMP)
algorithm, is used to solve the sparsity-constrained optimization problem.
Furthermore, the CV technique is adopted to determine the proper termination
condition by detecting overfitting when the sparsity level is unknown.