Open AccessArticle

An Off-Grid Turbo Channel Estimation Algorithm for Millimeter Wave Communications

Lingyi Han

Yuexing Peng

^1,*,

Peng Wang

² and

Yonghui Li

Wireless Signal Processing and Networks (WSPN) Lab, Key Lab of Universal Wireless Communications, Ministry of Education, Beijing University of Posts and Telecommunications, No. 10 Xitucheng Road, Beijing 100876, China

Huawei Technologies, Sweden AB, Kista 164 40, Sweden

Center of Excellence in Telecommunication, University of Sydney, Sydney NSW 2006, Australia

Author to whom correspondence should be addressed.

Sensors 2016, 16(10), 1562; https://doi.org/10.3390/s16101562

Submission received: 31 May 2016 / Revised: 3 September 2016 / Accepted: 7 September 2016 / Published: 22 September 2016

(This article belongs to the Special Issue Millimeter Wave Wireless Communications and Networks)

Download

Browse Figures

Figure 1
Block diagram of structure with beamforming at the base station (BS) and combining at the mobile station (MS). "> Figure 2
Flow chart of the iterative off-grid turbo channel estimation algorithm. "> Figure 3
Flow chart of the improved turbo compressed channel sensing algorithm. FFT and IFFT denote the fast Fourier transform processing and the inverse transform processing, respectively. The modules <math display="inline"> <semantics> <mi mathvariant="bold-italic">D</mi> </semantics> </math> and <math display="inline"> <semantics> <msup> <mi mathvariant="bold-italic">D</mi> <mi>H</mi> </msup> </semantics> </math> denote the transform from <math display="inline"> <semantics> <mi mathvariant="bold-italic">f</mi> </semantics> </math> to <math display="inline"> <semantics> <mi mathvariant="bold-italic">q</mi> </semantics> </math> and the reverse processing. "> Figure 4
The average angle estimation error (AAEE) versus signal-to-noise ratio (SNR) for the proposed algorithm for different values of N. "> Figure 5
The average probability of integral grid point estimation error (APIEE) versus signal-to-noise ratio (SNR) for the proposed algorithm for different values of N. ">

Versions Notes

Abstract

The bandwidth shortage has motivated the exploration of the millimeter wave (mmWave) frequency spectrum for future communication networks. To compensate for the severe propagation attenuation in the mmWave band, massive antenna arrays can be adopted at both the transmitter and receiver to provide large array gains via directional beamforming. To achieve such array gains, channel estimation (CE) with high resolution and low latency is of great importance for mmWave communications. However, classic super-resolution subspace CE methods such as multiple signal classification (MUSIC) and estimation of signal parameters via rotation invariant technique (ESPRIT) cannot be applied here due to RF chain constraints. In this paper, an enhanced CE algorithm is developed for the off-grid problem when quantizing the angles of mmWave channel in the spatial domain where off-grid problem refers to the scenario that angles do not lie on the quantization grids with high probability, and it results in power leakage and severe reduction of the CE performance. A new model is first proposed to formulate the off-grid problem. The new model divides the continuously-distributed angle into a quantized discrete grid part, referred to as the integral grid angle, and an offset part, termed fractional off-grid angle. Accordingly, an iterative off-grid turbo CE (IOTCE) algorithm is proposed to renew and upgrade the CE between the integral grid part and the fractional off-grid part under the Turbo principle. By fully exploiting the sparse structure of mmWave channels, the integral grid part is estimated by a soft-decoding based compressed sensing (CS) method called improved turbo compressed channel sensing (ITCCS). It iteratively updates the soft information between the linear minimum mean square error (LMMSE) estimator and the sparsity combiner. Monte Carlo simulations are presented to evaluate the performance of the proposed method, and the results show that it enhances the angle detection resolution greatly.

Keywords:

channel estimation; millimeter wave communication; off-grid algorithm

1. Introduction

Thanks to the large bandwidth available at millimeter wave (mmWave) frequencies, mmWave communication technology has become a promising technology to meet the experientially increasing demands of future wireless networks [1]. The small wavelength at the mmWave band enables the integration of a massive number of antennas at both the transmitter and receiver. The directional beamforming technique can be used to achieve a sufficient link margin to compensate for the severe propagation attenuation in mmWave bands [2,3]. Since the channel state information, especially the angle of departure (AoD) and angle of arrival (AoA), are essential for the beamforming and coherent detection, CE with high resolution and efficiency is one of the key requirements in mmWave systems. Due to the high cost and high power consumption of mmWave radio frequency (RF) devices, only a very limited number of RF chains can be integrated in the mmWave systems, which results in much less RF chains than antennas. Therefore, classic angle estimation methods with super-resolution, such as multiple signal classification (MUSIC) and estimation of signal parameters via rotation invariant technique (ESPRIT)-based algorithms [4,5], cannot be directly applied in mmWave systems because these methods require the channel observations at each receive antenna in the digital domain. However, in the mmWave systems, the received signals at each antenna are firstly combined in the analog domain and thus no direct output information from each antenna is available for these classic angle estimation methods. Therefore, new CE algorithms with high-resolution should be developed and tailored to mmWave communications.

Recently, many research efforts have been devoted to the CE for the mmWave systems [6,7,8,9,10]. As adopted in IEEE 802.15.3c standard [6], a polling mechanism is employed to select the best beam vector pair from the known codebooks with p beams at the transmitter and q beams at the receiver. This method consumes

p q

time slots to achieve angle resolution

O (1 / p)

and

O (1 / q)

. In [7], an adaptive compressive sensing (ACS) method has been developed, which iteratively bisections the beam space and tries to find the grid that the path falls in via space match filtering. It needs

2 L^{3} l o g_{2} (p / L) + 2 L^{3} l o g_{2} (q / L)

time slots to achieve the resolution

O (1 / p)

at the transmitter and

O (1 / q)

at the receiver for an mmWave channel with L paths. This large time slot payload is unacceptable in practical mobile communication systems with massive antenna arrays, especially in the outdoor wireless coverage scenarios. So far, most existing mmWave CE methods, such as [8,9,10], are based on the on-grid channel model, i.e., all path angles are assumed to be located at the quantization grids exactly, which is not the case in practical systems because the angles are usually continuously distributed. When path angles do not lie on the quantization grids, which we refer to as the “off-grid” channel model, a power leakage effect happens and results in the degradation of CE performance [11]. Since the product of the number of simultaneously supported multi-users and the number of the multi-streams depends on the RF chains, it is preferable that higher angle resolution is achieved by advanced channel estimation methods to separate densely distributed users such that spatial multiplexing methods can be employed to support more users [12]. In this ”off-grid” case, the gap to the theoretic lower bound is pretty large, and more antennas are required to enhance the resolution of angle estimation at the cost of more hardware and antennas [13]. Therefore, a new method is necessary to enhance the angle detection performance. So far, few works concentrate on the off-grid problem in the mmWave systems with RF constraints. In [14], the continuous angle estimation was formulated as a sparse regularization problem by taking into account the noises in both the observations and the dictionary by modeling the grid points as noisy, which is hard to be applied in CE because the spectral norm of the noise is much bigger than that of the dictionary. In [15], a continuous basis-based compressed sensing (CS) method was proposed to solve the reconstruction problem in sensing applications where the atoms of the dictionary are used to describe a continuous field, such as frequency and angle. However, the first order approximation of Taylor expansion for the array steering vector in the multi-antenna systems is hard to formulate.

In this paper, we study the off-grid problem and propose an enhanced method to improve the resolution of angle estimation in mmWave systems with massive antenna arrays and RF chain constraints. Firstly, we formulate the signal model to describe the off-grid problem in the angle quantization of the mmWave channel, which decomposes each continuously-distributed angle into the integral grid part and fractional off-grid part. Accordingly, an iterative off-grid turbo CE (IOTCE) algorithm is developed to renew and upgrade the CE between the integral grid and fractional off-grid angle information under the Turbo principle. By fully exploiting the sparse structure of mmWave channels, a soft-decoding based CS method, termed improved turbo compressed channel sensing (ITCCS), is developed to estimate the integral grid angle information, which iteratively updates the soft information between the linear minimum mean square error (LMMSE) estimator and the sparsity combiner. Simulation results show that the proposed CE method can achieve a higher resolution than the existing ones.

The rest of the paper is organized as follows. The mmWave system model is introduced in Section 2. Then, a new off-grid formulation is specified in Section 3 and the proposed IOTCE algorithm is detailed in Section 4. Section 5 presents the simulation results and Section 6 concludes the paper.

Notations: Fonts

A

a

, a and

A

denote a matrix, a column vector, a scalar and a set, respectively.

A^{T}

and

A^{H}

are the transpose and Hermitian of

A

1_{N} \in C^{N \times 1}

is an all-ones vector and

I

is an identity matrix.

e_{N} (i) \in C^{N \times 1}

is a unit vector with 1 at its i-th entry and 0 elsewhere. The operator ⊗ are used to denote the Kronecker product.

Vect (A)

is the matrix vectorization operator that transforms a matrix into a vector by stacking its column vectors, and

Diag (a)

stands for a diagonal matrix with the entries of

a

on its diagonal.

{[A]}_{R, P}

consists of row vectors and column vectors of the matrix

A

whose indices contained in the set

R

and

P

, respectively.

E (\cdot)

represents the expectation.

δ (\cdot)

is the Dirac delta function and

I (\cdot)

is the indicator function.

2. System Model

According to [16], typical mmWave channels feature much stronger line-of-sight (LoS) path than non-LoS (NLoS) paths, and the numerical results in [17] show that approximating the sparse mmWave channel with only one LoS path results in a little capacity loss when employing the steerable directional antennas. Thus, single-path mmWave channel model is used in this paper. When the mmWave channel cannot be approximated as a single-path channel, the proposed algorithm cannot be applied directly. New algorithms need to be developed to deal with independently but not identically distributed multi-paths. Consider a downlink mmWave system consisting of a base station (BS) and a mobile station (MS), where the BS is configured with

N_{B}

antennas and

K_{B}

RF chains with

K_{B} ≪ N_{B}

while the MS is configured with

N_{M}

antennas and

K_{M}

RF chains with

K_{M} ≪ N_{M}

, as shown in Figure 1. The BS applies the beamforming matrix

W_{B} \in C^{N_{B} \times K_{B}}

and the MS employs the matrix

W_{M} \in C^{N_{M} \times K_{M}}

to combine the received signal

y \in C^{K_{M} \times 1}

. The received signal at the MS can be written as:

\begin{matrix} y = W_{M}^{T} H W_{B} x + z \end{matrix}

(1)

where

H \in C^{N_{M} \times N_{B}}

is the channel matrix,

x \in C^{K_{B} \times 1}

represents the transmitted pilot signal, and

z \in C^{K_{N} \times 1}

is the Gaussian noise with

E ({zz}^{H}) = σ^{2} I

In this paper, we also assume:

uniform linear arrays (ULAs) with half-wavelength spacing are deployed at both the BS and MS;
far-field scattering and block-fading are held, which means the signal waves arrive at different antennas with the same fading amplitudes but distinct phases and the channel fadings are kept constant during the CE precedure.

Under these assumptions, the single-path channel between the BS and MS can be expressed as:

\begin{matrix} H = β a_{M} (θ_{M}) a_{B}^{T} (θ_{B}) \end{matrix}

(2)

where β is the fading weight,

θ_{B} \in [- \frac{π}{2}, \frac{π}{2})

, and

θ_{M} \in [- \frac{π}{2}, \frac{π}{2})

are, respectively, AoD and AoA.

a_{B}

and

a_{M}

are the array steering vector at the BS and the combining vector at the MS, respectively. Denoting

\frac{sin θ_{B}}{2}

ϕ_{B} \in [- \frac{1}{2}, \frac{1}{2})

and

\frac{sin θ_{M}}{2}

ϕ_{M} \in [- \frac{1}{2}, \frac{1}{2})

a_{B}

and

a_{M}

can be written as:

\begin{matrix} a_{B} (ϕ_{B}) = \frac{1}{\sqrt{N_{B}}} {[1 e^{j 2 π ϕ_{B}} \dots e^{j 2 π (N_{B} - 1) ϕ_{B}}]}^{T} \end{matrix}

(3)

\begin{matrix} a_{M} (ϕ_{M}) = \frac{1}{\sqrt{N_{M}}} {[1 e^{j 2 π ϕ_{M}} \dots e^{j 2 π (N_{M} - 1) ϕ_{M}}]}^{T} \end{matrix}

(4)

The

(m, n)

-th element of the channel can be rewritten as:

\begin{matrix} {[H]}_{m, n} & = β e^{j 2 π (m - 1) ϕ_{M}} e^{j 2 π (n - 1) ϕ_{B}} \\ = \int_{- \frac{1}{2}}^{\frac{1}{2}} \int_{- \frac{1}{2}}^{\frac{1}{2}} h ({\tilde{ϕ}}_{M}, {\tilde{ϕ}}_{B}) e^{j 2 π (m - 1) {\tilde{ϕ}}_{M}} e^{j 2 π (n - 1) {\tilde{ϕ}}_{B}} d {\tilde{ϕ}}_{M} d {\tilde{ϕ}}_{B} \end{matrix}

(5)

where

h ({\tilde{ϕ}}_{M}, {\tilde{ϕ}}_{B}) = β δ ({\tilde{ϕ}}_{M} - ϕ_{M}) δ ({\tilde{ϕ}}_{B} - ϕ_{B})

is the channel impulse response (CIR). By discretizing Equation (5) at the angle period

[- \frac{1}{2}, \frac{1}{2})

\begin{matrix} {[H]}_{m, n} (k, i) = \sum_{k = - N_{M} / 2}^{N_{M} / 2 - 1} \sum_{i = - N_{B} / 2}^{N_{B} / 2 - 1} h (k / N_{M}, i / N_{B}) e^{j 2 π (m - 1) (k + δ_{M}) / N_{M}} e^{j 2 π (n - 1) (i + δ_{B}) / N_{B}} \end{matrix}

(6)

where

δ_{M}

and

δ_{B}

are the discretization errors corresponding to the

ϕ_{M}

and

ϕ_{B}

, respectively. The on-grid algorithms, such as [6,7,8,9,10], can estimate the AoA/AoD effectively, by ignoring the discretization errors and assuming that the AoA/AoD are taken from a uniform grid of N points, i.e.,

ϕ \in {- \frac{1}{2}, - \frac{1}{2} + \frac{1}{N}, \dots, \frac{1}{2} - \frac{1}{N}},

where N stands for the number of antennas and ϕ stands for the angle. In practice, the number of antennas is limited and the AoA/AoD are actually continuous, which causes the leakage of power, and, in turn, leads to the degradation of CE performance.

3. Off-Grid Channel Formulation

In order to suppress the negative effect caused by the off-grid problem, a specific and operable virtual channel representation model is introduced. From this virtual model, an enhanced off-grid angle estimation model is developed.

The continuously-distributed AoA can be written as:

\begin{matrix} ϕ_{M} = - \frac{1}{2} + \frac{ϕ_{M, k}}{N_{M}} + \frac{δ_{M}}{N_{M}} \end{matrix}

(7)

where

ϕ_{M, k} = k - 1

is the nearest k-th (

k = 1, 2, \dots, N_{M}

) discrete quantization grid of

ϕ_{M}

, which is referred to as the integral grid angle, and

|δ_{M}| < \frac{1}{2}

stands for the deviation between

ϕ_{M}

and

ϕ_{M, k}

, termed fractional off-grid angle. Similarly, the discrete form of the AoD is formulated as:

\begin{matrix} ϕ_{B} = - \frac{1}{2} + \frac{ϕ_{B, i}}{N_{B}} + \frac{δ_{B}}{N_{B}} \end{matrix}

(8)

where

ϕ_{B, i} = i - 1

(

i = 1, 2, \dots, N_{B}

) and

|δ_{B}| < \frac{1}{2}

are, respectively, the integral grid part and fractional off-grid part of the angle. Then, the fractional off-grid angle vectors

d_{M} \in C^{N_{M} \times 1}

and

d_{B} \in C^{N_{B} \times 1}

can be rewritten as:

\begin{matrix} d_{M} (δ_{M}) = {[1 e^{j 2 π δ_{M} / N_{M}} \dots e^{j 2 π (N_{M} - 1) δ_{M} / N_{M}}]}^{T} \end{matrix}

(9)

\begin{matrix} d_{B} (δ_{B}) = {[1 e^{j 2 π δ_{B} / N_{B}} \dots e^{j 2 π (N_{B} - 1) δ_{B} / N_{B}}]}^{T} \end{matrix}

(10)

Accordingly, the channel matrix in Equation (2) can be reformulated as:

\begin{matrix} H & = β Diag {(d_{M} (δ_{M}))}^{T} a_{M} (- \frac{1}{2} + \frac{ϕ_{M, k}}{N_{M}}) a_{B}^{T} (- \frac{1}{2} + \frac{ϕ_{B, i}}{N_{B}}) Diag (d_{B} (δ_{B})) \\ = D_{M}^{T} H_{D} D_{B} \end{matrix}

(11)

where

H_{D} = β a_{M} (- \frac{1}{2} + \frac{ϕ_{M, k}}{N_{M}}) a_{B}^{T} (- \frac{1}{2} + \frac{ϕ_{B, i}}{N_{B}})

represents the channel matrix only containing the integral grid part, and the diagonal matrices

D_{M} = Diag (d_{M} (δ_{M})) \in C^{N_{M} \times N_{M}}

and

D_{B} = Diag (d_{B} (δ_{B})) \in C^{N_{B} \times N_{B}}

represent the fractional off-grid angle matrixes. From the virtual channel model in the [18],

H_{D}

can be rewritten as:

\begin{matrix} H_{D} = F_{M}^{T} {GF}_{B} \end{matrix}

(12)

where

F_{B} \in C^{N_{B} \times N_{B}}

is the discrete Fourier transform (DFT) matrix whose

(m_{B}, n_{B})

-th entry is

\frac{1}{\sqrt{N_{B}}} e^{j 2 π m_{B} \frac{n_{B}}{N_{B}}}

0 \leq m_{B}, n_{B} \leq N_{B} - 1

, and

F_{M} \in C^{N_{M} \times N_{M}}

is also a DFT matrix whose

(m_{M}, n_{M})

-th entry is

\frac{1}{\sqrt{N_{M}}} e^{j 2 π m_{M} \frac{n_{M}}{N_{M}}}

0 \leq m_{M}, n_{M} \leq N_{M} - 1

G \in C^{N_{M} \times N_{B}}

is referred to as the integral virtual channel matrix in angle domain with only one non-zero element, and

{[G]}_{k, i}

denotes the nearest integral grid point to the actual angle. The row and column labels of the non-zero element, k and i, are connected to the

ϕ_{M, k}

and

ϕ_{B, i}

via

\begin{matrix} ϕ_{M, k} = k - 1 \end{matrix}

(13)

\begin{matrix} ϕ_{B, i} = i - 1 \end{matrix}

(14)

Then, the virtual representation model of the mmWave channel with continuously-distributed angle can be written as:

\begin{matrix} H = D_{M}^{T} F_{M}^{T} G F_{B} D_{B} \end{matrix}

(15)

where

D_{M}

and

D_{B}

are unitary matrices. It is interesting that the channel matrix

H

can be regarded as a two-dimensional (2D) orthogonal transform from the spatial domain to the beamforming domain. More specifically, the 2D orthogonal transform becomes the 2D DFT when the AoA/AoD locate at the quantization grids exactly.

4. The Off-Grid Turbo Channel Estimation Algorithm

In order to facilitate the operation of the proposed CE algorithm, the codebook designed in [17] is employed to divide the AoA/AoD estimation into two steps to estimate AoA and AoD separately at the BS and MS in adjacent time slots. Moreover, it can also be adopted here to estimate the fractional off-grid parts of AoA and AoD. Without loss of generality, the codebook design for AoA estimation at the MS is introduced as an example, and the same process can be employed at the BS to estimate the AoD. The process of training mode design is summarized as follows:

Training mode at the transmitter. The training signal is designed to hold the form of $x = 1_{K_{B}}$ , and only one column of the beamforming matrix $W_{B}$ is chosen as ${[W_{B}]}_{:, n_{o}} = e_{N_{B}} (i_{o})$ , $n_{o} \in {1, 2, \dots, K_{B}}, i_{o} \in {1, 2, \dots, N_{B}}$ while all other columns are zero vectors, where ${[W_{B}]}_{:, n_{o}}$ denotes the $n_{o}$ -th column vector of matrix $W_{B}$ . The training signals $x$ are transmitted from $K_{B}$ RF chains and beamformed with the same $W_{B}$ in the successive R time slots.
Training mode at the receiver. The combining matrix at the r-th time slot, $W_{M}^{r}, r = 1, 2, \dots, R$ , is designed in a way that its n-th column vector is ${[W_{M}^{r}]}_{:, n} = e_{N_{M}} (i_{n}^{r}), n = 1, 2, \dots, K_{M},$ where ${i_{1}^{r}, i_{2}^{r}, \dots, i_{K_{M}}^{r}} ≜ P_{M}^{r} \subset R_{M}^{r},$ with $P_{M}^{r}$ being randomly chosen from $R_{M}^{r}$ and $R_{M}^{1} ≜ {1, 2, \dots, N_{M}}$ , $R_{M}^{2} = {1, 2, \dots, N_{M}} ∖ P_{M}^{1}, \dots, R_{M}^{r} = {1, 2, \dots, N_{M}} ∖ {P_{M}^{1} \cup \dots \cup P_{M}^{r - 1}}$ .

In this way, the received signal at the MS can be written as:

\begin{matrix} y_{r} = {(W_{M}^{r})}^{T} D_{M}^{T} F_{M}^{T} {GF}_{B} D_{B} W_{B} x + z \end{matrix}

(16)

Stacking

y_{r}

into a vector, we have

\begin{matrix} y & = [y_{1}^{T} y_{2}^{T} \dots y_{R}^{T}] \\ = W D F g + z \end{matrix}

(17)

where

W ≜ {[W_{M}^{1} W_{M}^{2} \dots W_{M K}^{R}]}^{T} \in C^{R K \times N},

with N and K being the number of antennas and RF chains, respectively.

F \in C^{N \times N}

is a DFT matrix, and

g \in C^{N \times 1}

is the integral virtual channel vector only containing the AoA of the LoS path seen by the MS.

D \in C^{N \times N}

is the fractional off-grid angle matrix corresponding to the AoA. Apparently,

W

is a permutation matrix created by choosing K rows from an identity matrix

I

From Equation (17), the CE is accomplished by estimating the integral virtual channel, which, in turn, defines integral grid angle, and the fractional off-grid angle matrix, which stands for the fractional off-grid angle. As shown in Figure 2, an iterative off-grid turbo CE (IOTCE) algorithm is developed, which consists of two components: the integral grid angle estimator and the fractional off-grid angle estimator.

4.1. The Integral Grid Angle Estimator

Equation (17) degenerates into a CS recovery algorithm with

W

being the sampling matrix and

D \times F

being the transform matrix when the fractional off-grid angle δ is given. As the application of the TCCS algorithm in [19], the ITCCS algorithm is adopted to iteratively estimate the integral virtual channel vector g from the Equation (17).

g \in C^{N \times 1}

is a sparse signal to be estimated with L non-zero elements, following the Bernoulli–Gaussian distribution with

λ = \frac{L}{N}

being the sparsity.

Let

f ≜ Fg

and

q ≜ D f

, and Equation (17) can be rewritten as:

\begin{matrix} y = W q + z \end{matrix}

(18)

As illustrated in Figure 3, the ITCCS algorithm consists of two modules: module A, the LMMSE estimator, and module B, the sparsity combiner. The LMMSE estimator provides an estimate of

q

by the linear minimum mean square estimation (LMMSE) method, from which an estimate of

g

can be obtained through constraints. Using the probability density estimation method, the sparsity combiner refines the estimate by exploiting the Bernoulli–Gaussian distribution. The specific realization of the algorithm can be seen in [19].

Therefore, by fixing δ and using the constraints Equations (9) and (10), we get the corresponding fractional off-grid angle vector

d (δ)

and matrix

D

, to estimate the integral virtual channel vector

g

at the sparsity

λ = \frac{L}{N}

(

L = 1

) using the ITCCS algorithm, where the value of δ is zero at the first iteration and is fed back from the fractional off-grid angle estimator at the iteration. For the purposes of this article, ITCCS

(W^{T}, Diag (d^{t}), \frac{1}{N}, σ^{2})

is used to make the description of the output, which is the estimate of

g

, and the parameters in brackets present the input of the algorithm.

4.2. The Fractional Off-Grid Angle Estimator

By fixing the integral virtual channel

g

, and as a result of the equivalence of an elementary row operation and premultiplication by a permutation matrix according to the theory of linear algebra, the Equation (17) can be rewritten as:

\begin{matrix} y & = W D f + z \\ = {[d (δ)]}_{R} \circ {[f]}_{R} + z \end{matrix}

(19)

where

R = {n_{1}, n_{2}, \dots, n_{K}}

is a set who is constituted by the row indices that are selected from

I

W

and the operator ∘ denotes the element-wise product. Equation (19) can be solved by the least square (LS) approach.

According to the proposed iteration, the IOTCE algorithm is illustrated in Algorithm 1, which presents a structure of a nested loop. In addition, the corresponding explanation of parameters in the IOTCE algorithm is listed in Table 1. The inner loop iterates the integral grid angle between sparsity and linearity, while the integral part and the fractional part are updated in turn during each iteration of the outer loop by following a fixed relationship. As the estimation of integral grid angle improves, the estimation of the fractional off-grid angle enhances via the outer loop.

Algorithm 1 The Iterative Off-grid Turbo Channel Estimation Algorithm

Input:

Received signal $y \in C^{R K \times 1}$ ;
Sampling matrix $W \in C^{R K \times N}$ , noise power $σ^{2}$ .

Output:

Estimate of $g$ , $g^{t}$ ;
Estimation of δ, $δ^{t}$ .

Iteration:

Step 1: Initialization

$\begin{matrix} t = 0, δ^{0} = 0 \end{matrix}$
Step 2: Estimation of g
- (1) Fixing δ,
  
  $\begin{matrix} d^{t} & = {[1 e^{j 2 π δ^{t} / N} e^{j 2 π 2 δ^{t} / N} \dots e^{j 2 π (N - 1) δ^{t} / N}]}^{T} \\ g^{t + 1} & = ITCCS (W^{T}, Diag (d^{t}), \frac{1}{N}, σ^{2}) \end{matrix}$
Step 3: Estimation of δ
- (1) Fixing $g$ ,
  
  $\begin{matrix} {[d^{t + 1}]}_{R} = y ⊘ {[g^{t + 1}]}_{R} \end{matrix}$
  
  where ⊘ denotes the element-wise division operation.
- (2) Calculating $δ^{t + 1}$ by linear fitting after phase extraction from the equation,
  
  $\begin{matrix} {[d^{t + 1}]}_{R} = {[e^{j 2 π (n_{1} - 1) δ^{t + 1} / N} e^{j 2 π (n_{2} - 1) δ^{t + 1} / N} \dots e^{j 2 π (n_{K} - 1) δ^{t + 1} / N}]}^{T} \end{matrix}$
Step 4: Iteration of loop
- IF iteration stopping condition is satisfied
  - output the estimated value.
- ELSE
  - $t = t + 1$ .
  - Then, and turn to Step 2.

5. Simulation Results

In this section, we evaluate the performance of the proposed algorithm in the mmWave communication systems. The simulation system consists of a BS and an MS, and the same ULA with

N_{B} = N_{M} = N

half-wavelength spaced antennas and

K_{B} = K_{M} = K

RF chains is configured at both the BS and MS. Typical mmWave channels featuring much stronger LoS paths than NLoS paths are simulated, and a one-path channel model is used to approximate the typical mmWave channel in the low signal-to-noise ratio (SNR) region. The fading weight of the path follows Rayleigh distribution with variance 1, and the continuously-valued AoA/AoD are uniformly-distributed in the range

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

The simulation results of the average angle estimation error (AAEE), defined as

E [| ϕ - \hat{ϕ} |],

with

\hat{ϕ}

being the estimate of the true angle ϕ, are presented in Figure 4, where the on-grid algorithm in [17] and the lower bound of the AAEE for on-grid algorithms, which is

\frac{1}{N}

, are also presented for reference in two cases of

N = 1024

and

N = 512

. From the curves in Figure 4, we can observe that the proposed method outperforms the on-grid method and achieves AAEE of 10⁻² at SNR γ = −4

dB

. As expected, the AAEE of the proposed algorithm exceeds the lower bound of the on-grid algorithms as SNR increases, which validates the AAEE performance of our algorithm. For example, the gaps between the simulated AAEE and the lower bound are 0.00139 and 0.0007 when γ = 20

dB

for

N = 512

and

N = 1024

, respectively. The gap shrinks from 0.00139 to 0.0007 by the increase in N, because the less sparsity (

λ = \frac{1}{N}

) can achieve a the better recovery performance at the given samplings according to the theory of CS [20].

The average probability of integral grid point estimation error (APIEE) performance, which is defined as

E [I (ϕ_{k} \neq {\hat{ϕ}}_{k})],

with

\hat{ϕ_{k}}

being the estimate of the real integral grid angle

ϕ_{k}

, is also simulated, and the simulation results are illustrated in Figure 5. As clearly shown, the bigger N leads to the better APIEE performance at a given K. For example, at the APIEE of 10⁻², the algorithm with

N = 1024

almost enhances by 4

dB

with respect to

N = 512

6. Conclusions

The shortage of spectrum at present will be alleviated by advances in the mmWave communication systems, and directional precoding/beamforming with large antenna arrays appears to be inevitable to support longer outdoor links and to provide sufficient received signal power than before. In this paper, an enhanced CE algorithm is developed for the off-grid problem in the mmWave systems with massive antenna arrays and RF chain constraints. First, we developed an off-grid formulation to catch the quantization problem. Then, an iterative method, which depends on the developed off-grid formulation is proposed to renew and upgrade the CE between the integral and fraction angle information under the Turbo principle. Simulation results show the efficiency and high resolution of the proposed method. For future work, it would be more practical to consider the multi-paths channel with independently but not identically distributed paths, and develop an efficient method to calculate the Log Likelihood Ratio (LLR) information.

Acknowledgments

The work of Yuexing Peng was supported by the NSFC under grant 61171106, the National Key Scientific Instrument and Equipment Development Project under Grant 2013YQ20060706 and the 863 Program with grant 2014AA01A705. The work of Yonghui Li was supported by the ARC grants DP150104019, LP150100994, and FT120100487 as well as by funding from the Faculty of Engineering & Information Technologies, The University of Sydney, under the Faculty Research Cluster Program.

Author Contributions

This article is the product of joint efforts. All the authors contributed to the conception, design and performance of the experiments, the analysis of the data and the writing of the paper. L. Han and Y. Peng conceived and designed the algorithm; P. Wang and Y. Li discussed and provided valuable comments on the algorithm design; L. Han performed the experiments and wrote the paper; Y. Peng, P. Wang and Y. Li revised the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

Smulders, P. Exploiting the 60 GHz band for local wireless multimedia access: Prospects and future directions. IEEE Commun. Mag. 2002, 40, 140–147. [Google Scholar] [CrossRef]
Lei, L.; Zhong, Z.; Zheng, K.; Chen, J.; Meng, H. Challenges on wireless heterogeneous networks for mobile cloud computing. IEEE Wirel. Commun. 2013, 20, 34–44. [Google Scholar] [CrossRef]
Zheng, K.; Zhao, L.; Mei, J.; Dohler, M.; Xiang, W.; Peng, Y. 10 Gb/s hetsnets with millimeter-wave communications: Access and networking—Challenges and protocols. IEEE Commun. Mag. 2015, 53, 222–231. [Google Scholar] [CrossRef]
Schmidt, R. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef]
Roy, R.; Kailath, T. ESPRIT estimation of signal parameters via rotational invariance techniques. IEEE Trans. Acoust. Speech Sign. Process. 1989, 37, 984–995. [Google Scholar] [CrossRef]
IEEE Std 802.15.3c. Available online: http://standards.ieee.org/getieee802/download/802.15.3c-2009.pdf (accessed on 12 October 2009).
Alkhateeb, A.; Ayach, O.E.; Leus, G.; Heath, R. Channel estimation and hybrid precoding for millimeter wave dellular systems. IEEE J. Sel. Top. Sign. Process. 2014, 8, 831–846. [Google Scholar] [CrossRef]
Ayach, O.E.; Rajagopal, S.; Abu-Surra, S.; Pi, Z.; Heath, R. Spatially sparse precoding in millimeter wave MIMO systems. IEEE Trans. Wirel. Commun. 2014, 13, 1499–1513. [Google Scholar] [CrossRef]
Peng, Y.; Li, Y.; Wang, P. An enhanced channel estimation method for millimeter wave systems with massive antenna arrays. IEEE Commun. Lett. 2015, 19, 1592–1595. [Google Scholar] [CrossRef]
Lee, J.; Gil, G.; Lee, Y. Exploiting spatial sparsity for estimating channels of hybrid MIMO systems in millimeter wave communications. In Proceedings of the 2014 IEEE Global Communications Conference (GLOBECOM), Austin, TX, USA, 8–12 December 2014; pp. 8–12.
Chi, Y.; Scharf, L.L.; Pezeshki, A.; Calderbank, A.R. Sensitivity to basis mismatch in compressed sensing. IEEE Trans. Sign. Process. 2011, 59, 2182–2195. [Google Scholar] [CrossRef]
Zheng, K.; Liu, F.; Lei, L.; Lin, C.; Jiang, Y. Stochastic performance analysis of a wireless finite-state Markov channel. IEEE Trans. Wirel. Commun. 2013, 12, 782–793. [Google Scholar] [CrossRef]
Zheng, K.; Fan, B.; Ma, Z.; Liu, G.; Shen, X.; Wang, W. Multihop cellular networks toward LTE-advanced. IEEE Veh. Technol. Mag. 2009, 4, 40–47. [Google Scholar] [CrossRef]
Zhu, H.; Leus, G.; Giannakis, G.B. Sparse regularized total least squares for sensing applications. In Proceedings of the 2010 IEEE 7th International Workshop on Signal Processing Advances in Wireless Communications (SPAWC), Marrakech, Morocco, 20–23 June 2010; pp. 1–5.
Ekanadham, C.; Tranchina, D.; Simoncelli, E.P. Recovery of sparse translation-invariant signals with continuous basis pursuit. IEEE Trans. Sign. Process. 2011, 59, 4735–4744. [Google Scholar] [CrossRef] [PubMed]
Rappaport, T.S.; Sun, S.; Mayzus, R.; Zhao, H.; Azar, Y.; Wang, K.; Wong, G.; Schulz, J.K.; Samimi, M.; Gutierrez, F. Millimeter wave mobile communications for 5G Cellular: It will work! IEEE Access 2013, 1, 335–349. [Google Scholar] [CrossRef]
Han, L.; Peng, Y.; Li, Y.; Zhao, H.; Zhao, J. Turbo compressed channel sensing for millimeter wave communications with massive antenna arrays and RF chain constraints. In Proceedings of the 2015 9th International Conference on Signal Processing and Communication Systems (ICSPCS), Cairns, Australia, 14–16 December 2015; pp. 1–6.
Sayeed, A.M. Deconstructing multiantenna fading channels. IEEE Trans. Sign. Process. 2002, 50, 2563–2579. [Google Scholar] [CrossRef]
Ma, J.; Yuan, X.; Ping, L. Turbo compressed sensing with partial DFT sensing matrix. IEEE Sign. Process. Lett. 2015, 22, 158–161. [Google Scholar] [CrossRef]
Davenport, M.; Laska, J.; Treichler, J.; Baraniuk, R. The pros and cons of compressive sensing for wideband signal acquisition: Noise folding versus dynamic range. IEEE Trans. Sign. Process. 2012, 60, 4628–4642. [Google Scholar] [CrossRef]

Figure 1. Block diagram of structure with beamforming at the base station (BS) and combining at the mobile station (MS).

Figure 2. Flow chart of the iterative off-grid turbo channel estimation algorithm.

Figure 3. Flow chart of the improved turbo compressed channel sensing algorithm. FFT and IFFT denote the fast Fourier transform processing and the inverse transform processing, respectively. The modules

D

and

D^{H}

denote the transform from

f

q

and the reverse processing.

D

and

D^{H}

denote the transform from

f

q

and the reverse processing.

Figure 4. The average angle estimation error (AAEE) versus signal-to-noise ratio (SNR) for the proposed algorithm for different values of N.

Figure 5. The average probability of integral grid point estimation error (APIEE) versus signal-to-noise ratio (SNR) for the proposed algorithm for different values of N.

Table 1. Explanation of parameters in the iterative off-grid turbo channel estimation algorithm.

**Table 1.** Explanation of parameters in the iterative off-grid turbo channel estimation algorithm.
Parameters	Explanation
N	the number of antennas at receiver
$g$	the integral virtual channel vector
$W$	the combining matrix at receiver
δ	the fractional off-grid angle
$d$	the corresponding fractional off-grid angle

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

Share and Cite

MDPI and ACS Style

Han, L.; Peng, Y.; Wang, P.; Li, Y. An Off-Grid Turbo Channel Estimation Algorithm for Millimeter Wave Communications. Sensors 2016, 16, 1562. https://doi.org/10.3390/s16101562

AMA Style

Han L, Peng Y, Wang P, Li Y. An Off-Grid Turbo Channel Estimation Algorithm for Millimeter Wave Communications. Sensors. 2016; 16(10):1562. https://doi.org/10.3390/s16101562

Chicago/Turabian Style

Han, Lingyi, Yuexing Peng, Peng Wang, and Yonghui Li. 2016. "An Off-Grid Turbo Channel Estimation Algorithm for Millimeter Wave Communications" Sensors 16, no. 10: 1562. https://doi.org/10.3390/s16101562

APA Style

Han, L., Peng, Y., Wang, P., & Li, Y. (2016). An Off-Grid Turbo Channel Estimation Algorithm for Millimeter Wave Communications. Sensors, 16(10), 1562. https://doi.org/10.3390/s16101562

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

Article Menu

An Off-Grid Turbo Channel Estimation Algorithm for Millimeter Wave Communications

Abstract

1. Introduction

2. System Model

3. Off-Grid Channel Formulation

4. The Off-Grid Turbo Channel Estimation Algorithm

4.1. The Integral Grid Angle Estimator

4.2. The Fractional Off-Grid Angle Estimator

5. Simulation Results

6. Conclusions

Acknowledgments

Author Contributions

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI