DOI QR코드

DOI QR Code

Low Complexity Hybrid Precoding in Millimeter Wave Massive MIMO Systems

  • Received : 2021.07.19
  • Accepted : 2022.03.18
  • Published : 2022.04.30

Abstract

As a preprocessing operation of transmitter antennas, the hybrid precoding is restricted by the limited computing resources of the transmitter. Therefore, this paper proposes a novel hybrid precoding that guarantees the communication efficiency with low complexity and a fast computational speed. First, the analog and digital precoding matrix is derived from the maximum eigenvectors of the channel matrix in the sub-connected architecture to maximize the communication rate. Second, the extended power iteration (EPI) is utilized to obtain the maximum eigenvalues and their eigenvectors of the channel matrix, which reduces the computational complexity caused by the singular value decomposition (SVD). Third, the Aitken acceleration method is utilized to further improve the convergence rate of the EPI algorithm. Finally, the hybrid precoding based on the EPI method and the Aitken acceleration algorithm is evaluated in millimeter-wave (mmWave) massive multiple-input and multiple-output (MIMO) systems. The experimental results show that the proposed method can reduce the computational complexity with the high performance in mmWave massive MIMO systems. The method has the wide application prospect in future wireless communication systems.

Keywords

1. INTRODUCTION

As the key technology of fifth generation mobile communication (5G), the massive multiple-input and multiple-output (MIMO) systems have hundreds of antennas to improve the spectral efficiency [1][2]. The millimeter wave (mmWave) ensures the feasibility of the massive MIMO systems in practical communication systems. Due to the limited wavelength of the millimeter waves, the Antenna in Package (AIP) technology can integrate hundreds or even thousands of system antennas into smaller chips or circuit boards [3]. Therefore, the combination of mmWave and massive MIMO systems is a promising technology.

In mmWave massive MIMO systems, the precoding technology is used to enhance the strength of the desired directional signal through the spatial processing to increase the system capacity [4]. The precoding technology adjusts the phase or amplitude of each element in the antenna array to obtain high transmission efficiency. However, the transmitter antennas of mmWave massive MIMO systems require many RF chains [5][6], which greatly increases the system complexity and the power consumption cost of the equipment [7][8]. Therefore, the hybrid precoding technology is the main solution in mmWave massive MIMO systems.

Based on the connection method, the hybrid precoding technology has three main categories. The first one is the fully-connected architectures-based hybrid precoding [9-13], and many scholars have made some researches about it. Such as, Ayach proposed sparse precoding, which transforms the achievable rate optimization problem into a sparse approximation problem and solves it by the matching pursuit algorithm, and its computation complexity is high [9]. Wonil et al. exploited codebooks to optimize the hybrid precoding by an iterative algorithm. The fully-connected architecture leads to a high computational complexity for iterative algorithms [10]. An alternative minimization algorithm based on the manifold optimization was introduced in [11]. The algorithm has the approximately optimal performances, but due to the manifold optimization, the digital and analog precoding matrixes are updated iteratively, so the computation complexity is still very high. All designed precoding for hybrid precoding employ a fully-connected architecture causing three additional limitations: 1) more phase shifters are required; 2) more energy is required to compensate for the energy consumption of the phase shifters; and 3) the computational complexity is higher, which leads to more energy consumption [9-14].

The second one is the sub-connected architecture-based hybrid precoding [15-16]. The each radio frequency (RF) chain is fixed to a certain portion of the antennas and can consume less power than the fully-connected structure. The successive interference cancelation (SIC) precoding has been proposed in sub-connected architectures [15]. The hybrid precoding optimization problem with non-convex constraints is transferred into a series of simple sub-antenna array optimization problems, each of which only considers one sub antenna array. However, the SIC precoding is a local-optimum for the sub-antenna array. The transmission efficiency is very limited. The hybrid precoding based on the alternating minimization algorithm has the similar drawback [16]. A hybrid precoding algorithm based on particle swarm optimization (PSO) has been introduced. The hybrid precoding uses the PSO algorithm to search the elements in the analog precoding matrix column by column. However, the method reduces the system performance to some extent [17]. Maximizing the spectral efficiency of the system is regarded as the objective function under the sub-connected architecture, and the optimization process of the analog precoding and the digital precoding are separated in two stages. However, the performance is still low [17-19]. Thus, the main concern of the sub-connected architecture is to improve the transmission performance.

The last one is the adaptive-connected structures-based hybrid precoding [20]. The adaptive-connected structure exploits antenna-selected technology to connect subset RF chains, which can consume less power. Some studies have exploited the dynamic connection technology to select RF chain and antenna arrays [21-24]. [23-25] proposed the hybrid precoding based on a sub-connected architecture, which uses a dynamic connection algorithm to select sub-arrays. The adaptive-connected structure can provide a trade-off between the transmission speed and the power consumption by additional switched and antenna selection algorithms, which require more system resources. In addition, some previous studies on the hybrid precoding only investigated how to maximize the spectral efficiency. However, the computational complexity is the key in the hybrid precoding of the wireless communication systems, especially in 5G systems.

To address the problems of high complexity and limited performance, a sub-connected architecture with the low complexity hybrid precoding is proposed in this paper. To obtain well transmission performance and avoid the local optimization problem, the extended power iteration (EPI) method is utilized to calculate the maximum eigenvalues and their eigenvectors of the channel matrix. The eigenvectors of the EPI method are used to obtain the hybrid precoding matrix with low computation complexity. The contributions of this paper can be summarized as follows.

1) The EPI method is proposed to calculate the maximum eigenvalues and their eigenvectors. The power iteration (PI) is used to calculate the maximum eigenvalue and its eigenvectors of the channel matrix. Then the deflation method is used to obtain the maximum eigenvalues and their eigenvectors. Based on the EPI method, the calculation speed of the hybrid precoding is improved.

2) The Aitken method is used to accelerate the convergence of the EPI method. When the maximum eigenvalue is close to the second eigenvalue, the convergence is very slow. Therefore, the Aitken method is used to solve this problem. Additionally, the Aitken method can increase the convergence speed of the EPI method.

3) Experimental results obtained with the proposed precoding are given for mmWave massive MIMO systems. The results show that the hybrid precoding based on the EPI method and the Aitken method has low complexity and a faster calculation speed than other methods.

The remainder of this paper is organized as follows. Section 2 briefly describes the mathematical model of mmWave massive MIMO systems. In Section 3, the hybrid precoding based on the EPI method is introduced in detail. The experimental results are provided in Section 4. The conclusions are given in Section 5.

2. MATHEMATICAL MODEL OF HYBRID PRECODING FOR MMWAVE MASSIVE MIMO

2.1 System Model

A sub-connected point-to-point architecture mmWave MIMO systems including one base station (BS) with MNRF antennas transmits N data streams to receivers with K antennas are considered in our paper, as shown in Fig. 1. [8]. The BS has NRF RF chains and N ≤ NRF < K < MNRF. Therefore the hybrid precoding is divided into two parts: a NRF × NRF baseband precoder FBB using N transmit chains and an MNRF × NRF RF precoder FRF using analog circuits. When we transmit N data streams, the vectors f(NRF-N) and f(NRF-N)B are zero vectors. Thus, the received signal can be modelled as

E1KOBZ_2022_v16n4_1330_f0001.png 이미지

Fig. 1. Sub-connected architecture for mmWave massive MIMO systems.

\(\begin{aligned}\mathbf{y}=\sqrt{\gamma} \mathbf{H F}_{R F} \mathbf{F}_{B B} \mathbf{s}+\mathbf{n}\\ \end{aligned}\)       (2)

\(\begin{aligned}\mathbf{F}_{R F}=\left[\begin{array}{lllc}\mathbf{f}_{1} & 0 & 0 & 0 \\ \mathrm{M} & \mathrm{O} & \mathrm{O} & \mathrm{M} \\ 0 & 0 & \mathbf{f}_{N} & \mathrm{M} \\ 0 & 0 & \mathrm{~L} & \mathbf{f}_{\left(N_{R F}-N\right)}\end{array}\right]_{M N_{R F} \times N_{R F}}\\\end{aligned}\)       (3)

where H∈ £K×MNRF is the matrix of the complex channel gain; S is a N×1 symbol vector with E[ssH]=IN/N, F = FRFFBB = diag[f1, f2,..., fN, f(NRF-N)]·diag[f1B, f2B,..., fNB, f(NRF-N)B] is a MNRF x N transmit hybrid precoding matrix in the sub-connected architecture; and FRF is an analog precoding matrix, which controls the direction of the transmitter antennas. FBB is the digital precoding, which controls the power of the transmitting antenna. n is the additive white Gaussian noise and γ = P/∥F∥F2 is a scaling factor to satisfy the transmitting power constraint. Compared to the fully-connected architecture with MNRF2 phase shifters (PS), the sub-connection architecture requires MNRF PSs.

In our system model, the extended Saleh Valenzuela cluster channel model is utilized to describe the mmWave channel characteristics [26-28]. In this channel model, the channel matrix H is assumed to be the sum of the contributions of L scattering clusters. Therefore, the channel matrix can be written as

\(\begin{aligned}\mathbf{H}=\sqrt{\frac{M N_{R F} K}{L}} \sum_{l=1}^{L} \alpha_{l} \Lambda_{r}\left(\phi_{l}^{r}, \theta_{l}^{r}\right) \Lambda_{t}\left(\phi_{l}^{t}, \theta_{l}^{t}\right) \mathbf{a}_{r}\left(\phi_{l}^{r}, \theta_{l}^{r}\right) \mathbf{a}_{t}\left(\phi_{l}^{t}, \theta_{l}^{t}\right)\\\end{aligned}\)       (4)

where αl denotes the complex gain of the lth path; Λrlr, θlr) and Λtlt, θlt) denote the antenna element gain for the arrival and departure angles, respectively [26]; and φt( φr) and θtr) denote the azimuth and elevation departure angles (arrival), respectively. arlr, θlr) and atlt, θlt) represent the receiving and transmitting antenna array response vectors, respectively.

The function of the antenna element gain Λ(φ, θ) can be expressed as follows:

\(\begin{aligned}\Lambda(\phi, \theta)= \begin{cases} 1 \quad \forall \phi \in\left[\phi_{\min }, \phi_{\max }\right], \forall \theta \in\left[\theta_{\min }, \theta_{\max }\right] \\ 0 & \text { otherwise }\end{cases} \end{aligned}\)       (5)

otherwise

where [φmin, φmax] denote the maximum and minimum azimuth angles of the receiving and transmitting antennas respectively. [θmin, θmax] denote the maximum and minimum elevation angles, respectively.

For the array response vector (ARV), the two commonly-used antenna arrays are expressed as follows:

For the uniform linear array (ULA), the ARV with M elements is given by

\(\begin{aligned}\mathbf{a}_{U L A}(\phi)=\frac{1}{\sqrt{M}}\left[1, e^{j \frac{2 \pi}{\lambda} d \sin (\phi)}, \ldots, e^{j(M-1) \frac{2 \pi}{\lambda} d \sin (\phi)}\right]\\\end{aligned}\)       (6)

In practical systems, the azimuth angle is quantized as an array antenna as follows:

\(\begin{aligned}\phi=\left\{\phi_{\min }+\frac{\phi_{\max }-\phi_{\min }}{Q}, \phi_{\min }+2 \frac{\phi_{\max }-\phi_{\min }}{Q}, \ldots, \phi_{\max }-\frac{\phi_{\max }-\phi_{\min }}{Q}\right\}\\\end{aligned}\)       (7)

where Q denotes the number of phase shifter quantization levels.

For a uniform planar array (UPA) with W1 and W2 elements on horizon and vertical, respectively, the ARV is expressed as follows:

\(\begin{aligned}\mathbf{a}_{U P A}(\phi, \theta)=\frac{1}{\sqrt{W_{1} W_{2}}}\left[1, e^{j \frac{2 \pi}{\lambda} d\left(w_{1} \sin (\phi) \sin (\theta)+w_{2} \cos (\theta)\right)}, \ldots, e^{j \frac{2 \pi}{\lambda} d\left(\left(W_{1}-1\right) \sin (\phi) \sin (\theta)+\left(W_{2}-1\right) \cos (\theta)\right)}\right]\\\end{aligned}\)       (8)

where 0 ≤ w1 ≤ W1, 0 ≤ w2 ≤ W2. (φ, θ) is quantized as follows:

\(\begin{aligned}\phi=\left\{\phi_{\min }+\frac{\phi_{\max }-\phi_{\min }}{Q_{1}}, \phi_{\min }+2 \frac{\phi_{\max }-\phi_{\min }}{Q_{1}}, \ldots, \phi_{\max }-\frac{\phi_{\max }-\phi_{\min }}{Q_{1}}\right\}\\\end{aligned}\)       (7)

\(\begin{aligned}\theta=\left\{\theta_{\min }+\frac{\theta_{\max }-\theta_{\min }}{Q_{2}}, \phi_{\min }+2 \frac{\theta_{\max }-\theta_{\min }}{Q_{2}}, \ldots, \theta_{\max }-\frac{\theta_{\max }-\theta_{\min }}{Q_{2}}\right\}\\\end{aligned}\)       (10)

where Q1 and Q2 denote the number of phase shifter quantization levels.

2.2 Hybrid Precoding Design

We consider the design of the hybrid precoding for point-to-point mmWave massive MIMO systems. For the sub-connected architecture, the achievable rate can be expressed as follows:

\(\begin{aligned}R=\log _{2}\left(\left|\mathbf{I}+\frac{\gamma}{N} \mathbf{H F}_{R F} \mathbf{F}_{B B} \times \mathbf{F}_{B B}^{H} \mathbf{F}_{R F}^{H} \mathbf{H}^{H}\right|\right)\\\end{aligned}\)       (11)

To obtain the maximum achievable rate, FRF, FBB are designed under the following constraints:

A) The total transmit power constraint should satisfy ∥FRFFBB∥≤ N

B) The PS requires that the analogy vectors have the same amplitude: fi ∈ A.

where A = [a(φ1, θ1), a(φ2, θ2), ..., a(φi, θi), a(φL, θL)], a(φi, θi) is an M×1 ARV, and all ARVs are set as a similar matrix for every RF chain.

However, the non-convex constraints on FRFFBB to maximize the achievable rates are complicated. Therefore, in some studies, the achievable rate has been maximized by reducing the “distance” between FRFFBB and the optimal unconstrained precoding matrix. For hybrid architectures, many studies have proved that the design of the analog precoding matrix and the digital precoding matrix is developed from the optimal unconstrained precoding matrix by using sparse signal processing techniques. Ayach proved that the right-singular vectors of the channel are the optimal unconstrained vectors due to the anti-interference ability of SVD precoding [9], which inspired us to make FRFFBB approach the digital SVD precoding matrix in a sub-connected architecture. The optimal unconstrained precoding matrix can divide the channel into an independent channel to reduce the interference between the transmit antennas. In real mmWave communication (N ≤ NRF ≤ MNRF), when we transmit N data streams by N RF-chain, the hybrid precoding design can be decomposed as follows:

H = UΣVH       (12)

\(\begin{aligned}\boldsymbol{\Sigma}=\left[\begin{array}{cc}\boldsymbol{\Sigma}_{1} & 0 \\ 0 & \boldsymbol{\Sigma}_{2}\end{array}\right], \mathbf{V}=\left[\mathbf{V}_{1}, \mathbf{V}_{2}\right]\\\end{aligned}\)       (13)

where Σ1 is N × N and V1 = [v1, v2, ...vi, ...vN].

Due to the special sub-connected architecture and the independence of the vectors in V1, we can split the right-singular matrix V1 into N vectors. Then we exploit the N RF chain and N antenna sub-arrays to approach the independence vector vi from V1.

When fiBfi is close to vi,

∥vi - fiBfi22 = (vi - fiBfi)H (vi - fiBfi)

= viHvi + (fiB)2fiHfi - 2fiB Re(viHfi)

= 1 + (fiB)2 - 2fiB Re(viHfi)

= 1 + (fiB)2 - 2fiB Re(viHfi)

= (fiB - Re(viHfi))2 + (1 + Re2(viHfi))       (14)

Because Vi is a unitary matrix and fi has the same power allocation, we obtain FRF, FBB

\(\begin{aligned}\mathbf{f}_{i}=\frac{1}{\sqrt{M}} jangle \left(\mathbf{v}_{1}\right)\\\end{aligned}\)       (15)

\(\begin{aligned}f_{i}^{B}=\frac{\left\|\mathbf{v}_{i}\right\|}{\sqrt{M}}\\\end{aligned}\)       (16)

In practical systems, the ARV is quantized. Therefore, the elements of FRF are derived from A ; thus,

\(\begin{aligned}\mathbf{f}_{i}=\underset{\mathbf{a}^{(m)} \in \mathbf{A}}{\arg \max }\left\|\operatorname{angle}\left(\mathbf{v}_{i}\right) \operatorname{angle}\left(\mathbf{a}^{(m)}\right)\right\|\\\end{aligned}\)       (17)

From (13), (16), and (17), the new hybrid precoding must calculate the N largest eigenvalues and their eigenvectors. It is necessary to design, FRF, FBB. In addition, the rank of the mmWave channel is equal to that of L scattering clusters; if L < N, we can only transmit L data streams to ensure that the independent eigenvectors guarantee the communication performance. It should be noted that the complexity of SVD procedures is substantial when MNRF is large, and so only the N largest eigenvalues and their eigenvectors are used. Therefore, a low computation complexity method is used to solve these constraints.

3. EXTENDED POWER ITERATION PRECODING

To avoid the SVD procedure, we must identify an approximation method to obtain the N largest eigenvalues and their eigenvectors of the channel matrix. This paper proposes a novel EPI method to solve this problem. Before the EPI method, the PI method is given as follows.

3.1 Power Iteration Method

Let A∈£nxn be a matrix with eigenvalues [29][30]:

1| > |λ2| ≥...≥|λn|       (18)

where λ1 is a maximum eigenvalue, then we assume that A has a set of independent eigenvectors x1, x2, ..., xn of unit length associated with the eigenvalue.

Assume an initial vector p0 = p. The power method forms the sequence of vectors p1 = Ap, p2 = A2p, ..., and the recursion can be expressed as follows:

pk = Apk-1 = A2pk-2 = A3pk-3 = Akp0       (19)

Then, by expanding the initial vector along the eigenvectors p0 = ∑j=1nαjxj and assuming that α1 ≠ 0, we obtain pk for k =1,2,3, ... .

\(\begin{aligned} \mathbf{p}_{k} & =\mathbf{A}^{k} \mathbf{p}_{0}=\sum_{j=1}^{n} \lambda_{j}^{k} \alpha_{j} \mathbf{x}_{j} \\ & =\lambda_{1}^{k}\left(\alpha_{1} \mathbf{x}_{1}+\sum_{j=2}^{n}\left(\frac{\lambda_{j}}{\lambda_{1}}\right)^{k} \alpha_{j} \mathbf{x}_{j}\right) \\ & =\lambda_{1}^{k} \alpha_{1} \mathbf{x}_{1}+\mathrm{O}\left(\left|\frac{\lambda_{2}}{\lambda_{1}}\right|^{k}\right)\end{aligned}\)       (20)

In practice, the l2-norm of p0 is set as 1 to avoid overflow or underflow, and the initial recursion is modified as follows:

JAKO202213841062561_21.png 이미지       (21)

Due to |λ1| > |λ2|, when k → ∞,

\(\begin{aligned} \begin{cases} \begin{array}{l}\lambda_{1} \approx m_{k} \\ \mathbf{x}_{1} \approx \mathbf{p}_{k} \end{array} \end{cases} \end{aligned}\)       (22)

The PI method yields the maximum singular value λ1 and the first right eigenvectors x1 by using (22). However, this PI method computes only the largest eigenvalue and its corresponding eigenvector. Moreover, the sparse precoding needs the N largest eigenvalues and their eigenvectors. Therefore, the EPI method is proposed to calculate the maximum N eigenvalues and their eigenvectors.

3.2 Hybrid Precoding Based on the EPI Method

Because PI only can calculate the maximum value, it is difficult to meet the hybrid precoding requirements. Therefore, we derive the EPI method from the PI method by using the deflation method. Based on the PI method, we delete the maximum eigenvalues and their eigenvectors from the matrix through the deflation method. Then, we obtain new sets of maximum eigenvalues and their eigenvectors from the new matrix. The EPI method for the channel matrix is presented as follows:

Step 1: Calculate the first maximum eigenvalue and eigenvector of H(i) as follows:

z(n-1) = H(i)u(n-1)       (23)

λ(n-1) = argmax|zi(n-1)|       (24)

\(\begin{aligned}\mathbf{u}^{(n)}=\frac{\mathbf{z}^{(n-1)}}{\lambda^{(n-1)}}, n  < n'\end{aligned}\)       (25)

Step 2: Deflate the channel matrix H(i) by the deflation method:

\(\begin{aligned}\mathbf{H}^{(i+1)}=\mathbf{H}^{(i)}-\frac{\lambda_{i}^{(n)}}{\left\|\mathbf{u}_{i}^{(n)}\right\|_{2}^{2}}\left(\mathbf{u}_{i}^{(n)}\right)^{T} \mathbf{u}_{i}^{(n)}\\\end{aligned}\)       (26)

When a new matrix H(i+1) is obtained, skip back to step 1. The eigenvalue and eigenvector of H(i+1) are calculated using (24) and (26). The proofs of the deflation method in (26) are presented in the APPENDIX I.

3.3 EPI Method Optimization Based on the Aitken Method

Based on (20), although the EPI method is convergent, its convergence rate may be slow when |λ1| ≈ |λ2|. Therefore, the Aitken method [31][32] is used to accelerate the convergence rate. We can describe the overall procedure of the Aitken iterations as follows:

JAKO202213841062561_27.png 이미지       (27)

JAKO202213841062561_28.png 이미지       (28)

JAKO202213841062561_29.png 이미지       (29)

The geometric interpretation of the Aitken acceleration method is shown in Fig. 2. Suppose that x0 is an approximate root of equation y = g(x), where x1(1) = g(x0), x1(2) = g(x1(1)). The point pair can give P0(x0, x1(1))and p1(x1(1), x1(2)), and the cross point of the chain P0P1 and function y = x is x1, which comes from the Aitken method. The proofs of the Aitken method are presented in APPRODIX II.

By the Aitken method, the (24) can be transformed to

\(\begin{aligned}\lambda^{n}=\begin{cases} \begin{array}{lr}m^{(n)}=\arg \max \left|z_{i}^{(n)}\right|, & \text { for } 1 \leq n \leq 2 \\ \frac{m^{(n)} m^{(n-2)}-\left(m^{(n-1)}\right)^{2}}{m^{(n)}-2 m^{(n-1)}+m^{(n-2)},} & \text { elsewise }\end{array} \end{cases} \end{aligned}\)       (30)

E1KOBZ_2022_v16n4_1330_f0002.png 이미지

Fig. 2. Geometric interpretation of the Aitken method.

3.4 Procedure of Hybrid Precoding Based on EPI Method

Through N iterations, we can obtain the largest eigenvalue and the corresponding eigenvector of the channel matrix H, FBB and FRF are accepted from the EPI method. Those specific steps in the EPI method are as follows in Algorithm 1 and Fig. 3.

E1KOBZ_2022_v16n4_1330_f0003.png 이미지

Fig. 3. A block diagram of EPI hybrid precoding

Algorithm 1. Expansion power iteration method

Require: H(1) = H, u0, n', N, M A

1 : Do :

2 : For 1 ≤ n ≤ n′

3 : z(n) = H(i)u(n-1), m(n) = arg max|z(n)|

4 : If 1 ≤ n ≤ 2 : λi(n) = m(n).

5 : Else: \(\begin{aligned}\lambda_{i}^{(n)}=\frac{m^{(n)} m^{(n-2)}-\left(m^{(n-1)}\right)^{2}}{m^{(n)}-2 m^{(n-1)}+m^{(n-2)}}, \mathbf{u}_{i}^{(n)}=\frac{\mathbf{z}^{(n)}}{\lambda_{i}^{(n)}}\\\end{aligned}\).

6 : End for

7 : \(\begin{aligned}f_{i}^{B}=\left\|\mathbf{u}_{i}^{(n)}\right\| / \sqrt{M}\\\end{aligned}\)

8 : \(\begin{aligned}\mathbf{f}_{i}=\underset{\mathbf{A}^{(m)} \in \mathbf{A}}{\arg \max }\left\|\operatorname{angle}\left(\mathbf{u}_{i}^{(n)}\right) \operatorname{angle}\left(\mathbf{A}^{(m)}\right)\right\|\\\end{aligned}\)

9 : \(\begin{aligned}\mathbf{H}^{(i+1)}=\mathbf{H}^{(i)}-\frac{\lambda_{i}^{(n)}}{\left\|\mathbf{u}_{i}^{(n)}\right\|_{2}}\left(\mathbf{u}_{i}^{(n)}\right)^{T} \mathbf{u}_{i}^{(n)}\\\end{aligned}\)

10: Until i > N > or |λi(n)| ≤ ε?

11: Return FBB, FRF

3.5 Complexity Analysis

In this section, we compare the computational complexity of the SVD, the EPI, the dynamic sub-array (DS)[23], and the sparse precoding. The number of floating point operations (flops) needed to implement these complex operations varies greatly according to the used hardware and the complex number representation [33]. Our EPI precoding needs to compute the optimal victors u(n) through the EPI method according to [31]. The matrix multiplication operation of z(n) = H(i)u(n-1) and the Aitken method is O (MNRFK). In addition, the division operation of the EPI method is O(2n′ − 2). Therefore, the total operation of the EPI method is O(NNRFMKn') + O(2n'N). the SVD and the DS precoding have similar computation complexity, i.e., O((NRFM)3). The computation complexity of the sparse precoding method is O(NRF4M + NRF2L2 + NRF2M2L). All hybrid precoding complexities are shown in Table 1.

Table 1. Complexity of five hybrid precoding in mmWave massive MIMO systems

E1KOBZ_2022_v16n4_1330_t0001.png 이미지

In mmWave systems, the number of iterations n and the number of RF chains N are less than NM, M, and L. and the number of BS antennas NM is far larger than the other parameters. The critical parameter in the sparse precoding complexity is NM. Therefore, the SVD of the sparse precoding is the main contributor to the computational complexity. As shown in Fig. 4, the complexity of the EPI method is lower than that of the sparse precoding.

E1KOBZ_2022_v16n4_1330_f0004.png 이미지

Fig. 4. The complexity comparisons of sparse, EPI, and SIC precoding for a (L = 32, n=5) mmWave MIMO systems

Fig. 4 shows the complexity of antennas per RF chain and RF chains for mmWave MIMO systems, considering mmWave massive MIMO systems with L=32, n′ = 5, K=16, and NRF=16. Based on Fig. 4, the SIC precoding has the lowest complexity, but the lowest achievable rate cannot satisfy the communication requirement. The EPI precoding has a lower complexity than the sparse precoding. The precoding complexity has a light relationship with M; the main complexity is determined by the number of data streams N. In other words, when we increase the number of antennas per RF chain, the computation complexity has limited influence, but it can improve the achievable rate.

4. SIMULATION AND ANALYSIS

To demonstrate the performance of the hybrid precoding, we consider the classical mmWave massive MIMO systems. We use the software defined ratio (SDR), the channel emulation (CE), and the upconverter for our simulation experiment. The CE is used to generate real communication channel information with 8 XILINX KINTEX7 FPGAs; the modulation and the EPI precoding are deployed in the SDR systems. The data stream is upconverted into CE, and all experimental data is analyzed in MATLAB 2014a on a personal computer with 3.3 GHz quad-core processor with 16 GB of RAM. The simulation parameters are shown in Table 2. The ULA is used for simulation with the angle of departure (AoD) randomly distributes in [−π/6, π/6] and the direction of departure (DoD) randomly distributes in [−π, π]. Therefore, we obtaine the ARV from (5).

Table 2. The experiment simulation parameter

E1KOBZ_2022_v16n4_1330_t0002.png 이미지

4.1 Analysis of Achievable Rate with Antennas per RF Chain M.

In this section, the experimental results of the antennas per RF chain are analyzed and compared with the results of five different precoding techniques. Fig. 5 and Fig. 6 show the achievable rates of the SVD precoding, the sparse precoding, the DS precoding, the EPI precoding, and the SIC precoding with different M and signal-to-noise ratio (SNR) values.

E1KOBZ_2022_v16n4_1330_f0005.png 이미지

Fig. 5. Achievable rate of five precoding methods with the number of transmitter per RF chain antennas and SNR( N=8).

E1KOBZ_2022_v16n4_1330_f0011.png 이미지

Fig. 6. Achievable rate comparison for 128 × 16 mmWave MIMO systems (N=8).

As shown in Fig. 5, the number of RF chains N is set to 8. We obtaine different achievable rates with various M. The fully digital SVD precoding has the best performance comparing to the other precoding methods under all SNRs and BS antennas. The sparse precoding has the next highest performance and the DS precoding shows a performance similar to that of the sparse precoding. The SIC precoding has low performance due to the local-optimum of the sub-connected architecture. The proposed precoding has a similar performance to the sparse precoding. Due to the sparsity of the massive MIMO channel matrix, the maximum eigenvalues and their eigenvectors represent the main components of the channel matrix. Thus, the EPI precoding has a better performance by using the top eigenvalues and their eigenvectors.

Fig. 6 shows the achievable rate of five precodings with MN=128, N=8, and SNR=0 dB. Similar to Fig. 5. the performance in descending order is the SVD precoding, the sparse precoding, the DS precoding, the EPI precoding, and the SIC precoding. However, the EPI precoding is close to the performance of the SVD precoding when the SNR value is large, but the computation complexity of the EPI precoding is lower than that of the sparse precoding.

4.2 Analysis of Achievable Rate with Data Streams N.

In this section, the experimental results of the data streams N are analyzed with the result five different precoding. Fig. 7 and Fig. 8 show the achievable rate of five precoding when M=16. The performances are like that shown in Fig. 5, the EPI precoding is close to the performance of the sparse precoding and the DS precoding when SNR>0 dB. The achievable rate of the sparse precoding and the EPI precoding reach 102 bps/Hz when the SNR is 30 dB. Under a low SNR environment, those five precoding has small performance gaps due to high noise power, which limits the total achievable rate. Especially in Fig. 8, as the SNR increasing, all precoding had a better achievable rate. The “gap” between the sparse precoding and the EPI precoding decreases as the SNR increases. However, the EPI precoding has a lower cost and computation complexity than the sparse precoding.

E1KOBZ_2022_v16n4_1330_f0006.png 이미지

Fig. 7. Achievable rate of five precoding methods with different data stream N and SNR(M=16).

E1KOBZ_2022_v16n4_1330_f0012.png 이미지

Fig. 8. Achievable rate comparison for 256 × 16 mmWave MIMO systems (M=16).

Fig. 9 shows the achievable rate variation against the number of different data streams N, when M =16 and SNR = 0 dB. Based on the Fig. 9, the performance of the EPI precoding is worse than that of the sparse precoding, the DS precoding, and the digital SVD precoding. The fully-connected architecture hybrid precoding has the better performance than the DS and the sub-connected architecture. However, the sub-connected architecture has a lower complexity and cost.

E1KOBZ_2022_v16n4_1330_f0007.png 이미지

Fig. 9. Achievable rate comparison against the number of the data streams (SNR=0dB, L=32).

To compare the cost and complexity of the different precoding, the performance comparisons of those five precoding are given in Table 3. As shown in Table 3, the four hybrid precoding architectures including the full-digital (SVD precoding), the full-connected (sparse precoding), the dynamic subarray (DS precoding) and the sub-connected (EPI and SIC precoding) architectures, the RF chains, the PSs, the complexity, and the achievable rate are compared for the mmWave massive MIMO systems with M=16, NRF=16, L=32, n=5, SNR=10 dB and N=16. As shown in Table 3, the SVD precoding required 256 RF chains, and the RF chain and PS are costly for communication systems; thus the cost of the SVD precoding is immense. Although the fully-connected architecture requires the similar RF chains, it needs more PSs than the dynamic subarray and sub-connected architecture. the DS precoding has a better performance than the EPI precoding, but its computational complexity is enormous. Although the EPI precoding and SIC precoding are based on a sub-connected architecture and the complexity of the EPI precoding is greater than that of the SIC precoding, the achievable rate of the EPI precoding (45.84bps/Hz) is greater than that of the SIC precoding (15.12bps/Hz) with sub-connected architecture. Overall, the EPI precoding provides a trade-off between transmission rate, the computation complexity and the cost. Thus, in a real communication environment, the hybrid precoding based on the EPI method is suitable for mmWave communication systems.

Table 3. The performance comparison of five precodings

E1KOBZ_2022_v16n4_1330_t0003.png 이미지

4.3 Analysis of Energy Efficiency

The energy efficiency of the hybrid precoding is expressed as the ratio of the achievable rate and energy consumption to express the trade-off between these two factors [25]:

\(\begin{aligned}E \stackrel{\Delta}{=} \frac{R}{P_{t o t}}\\\end{aligned}\)       (31)

where Ptot is the total power consumption of the transmitter. The total power consumption for different hybrid precoding architectures is defined in Table 4.

P is the transmit power, PBB is the digital power, PRF is the RF power, PPS is the PS power, and PSW is the switch network power. In our simulation, P is 10w, PBB is 200mW, PRF is 300mW, PPS is 50mW, and PSW is 5mW. The energy efficiency of different precoding is given in Fig. 10. Due to the dynamic subarray architecture, the DS precoding has the best energy efficiency but its complexity is immense compared to that of other precoding. Additionally, the sub-connected architecture has a better energy efficiency than the sparse precoding with a fully-connected architecture. The EPI precoding has a better energy performance than the SIC precoding and the sparse precoding. Therefore, the proposed low complexity precoding is more energy efficient for mmWave massive MIMO systems.

Table 4. Transmit power consumption for different architecture in mmWave massive MIMO systems

E1KOBZ_2022_v16n4_1330_t0004.png 이미지

E1KOBZ_2022_v16n4_1330_f0008.png 이미지

Fig. 10. Average energy efficiency for different architecture in 128 × 16 mmWave MIMO system.

4.4 Analysis of Achievable Rate with Imperfect CSI.

The effect of the imperfect CSI is evaluated with the EPI hybrid precoding. The imperfect CSI \(\begin{aligned}\hat{\mathbf{H}}\\\end{aligned}\) can be modelled as [1]

\(\begin{aligned}\hat{\mathbf{H}}=\xi \mathbf{H}+\sqrt{1-\xi^{2}} \mathbf{E}\\\end{aligned}\)       (32)

where H is the perfect channel matrix, ξ ∈[0,1]is the CSI accuracy, and E is the error matrix with the distribution i.i.d. CN (0, 1). Fig. 11 shows the sum achievable rate comparison for an MN × K = 256 × 16 (N =8) mmWave MIMO system, for which different ξ scenarios are simulated. The EPI hybrid precoding is robust to different ξ. As shown in Fig. 11, the achievable rate of the EPI precoding approaches the perfect CSI scenario (ξ = 1) when ξ = 0.9. Even when the CSI accuracy is quite poor (i.e., ξ = 0.5), the EPI precoding achieves more than 80% sparse precoding performance in the perfect CSI scenario. Therefore, the EPI precoding has the better performance in the imperfect CSI scenario.

E1KOBZ_2022_v16n4_1330_f0009.png 이미지

Fig. 11. Achievable rate comparison of different ξ for an MN × K= 256 × 16 (M = 16) mmWave MIMO system.

4.5 Analysis of Achievable Rate with Scatting Cluster L.

The rank of the mmWave channel is equal to the number of scattering clusters L; if L is small, the number of data streams is N, and the performance of the hybrid precoding decreases, because the independent eigenvectors of the channel matrix cannot support transmitting N data streams. To explain the relationship between the spectral effectiveness and the number of scattering clusters L < N, the simulation results are shown in Fig. 12. The number of data streams is set to 16. When the scattering clusters are small at 16, the achievable rate increases as L increases because when L > N, the performance of the hybrid precoding is not limited by the eigenvectors of the channel matrix. The number of data streams N influences the achievable rate. As shown in Fig. 9, when L=32, the achievable rate increases as N increases.

E1KOBZ_2022_v16n4_1330_f0010.png 이미지

Fig. 12. Achievable rate comparison of different scattering clusters mmWave MIMO systems (M =16, N=16, and SNR=0dB).

5. Conclusion

n this paper, we propose a novel hybrid precoding for mmWave massive MIMO systems. First, we utilize the hybrid precoding to approach the maximum eigenvectors of the channel matrix, with the aim of obtaining the maximum achievable rate in the sub-connected architecture. Then we introduce the EPI method to obtain the maximum eigenvalues and their eigenvectors to avoid the SVD procedure. Then the Aitken method is used to accelerate the convergence of the EPI method. Based on the EPI and Aitken methods, we can apply hybrid precoding in mmWave massive MIMO systems with low computation complexity. Finally, the simulation results show that the EPI algorithm has a near-optimal performance approaching sparse precoding and digital SVD precoding under low levels of computational complexity.

In addition, some further investigations will focus on two directions. 1) The impacts of the SNR on our EPI method will be studied; and we will conduct deeper research on filter technology before the hybrid precoding. 2) Quantitative analysis will be given more attention between lower complexity and faster computational speed in our upcoming work.

Appendix

References

  1. F. Rusek et al., "Scaling Up MIMO: Opportunities and Challenges with Very Large Arrays," IEEE Signal Proce. Mag., vol. 30, no. 1, pp. 40-60, Jan. 2013. https://doi.org/10.1109/MSP.2011.2178495
  2. R. Zi, X. Ge, J. Thompson, C.-X. Wang, H. Wang, and T. Han, "Energy effciency optimization of 5G radio frequency chain systems," IEEE J. Sel.Areas Commun., vol. 34, no. 4, pp. 758-771, Apr. 2016. https://doi.org/10.1109/JSAC.2016.2544579
  3. J. Park, D. Choi and W. Hong, "Millimeter-Wave Phased-Array Antenna-in-Package (AiP) Using Stamped Metal Process for Enhanced Heat Dissipation," IEEE Antennas and Wireless Propagation Letters, vol. 18, no. 11, pp. 2355-2359, Nov. 2019. https://doi.org/10.1109/lawp.2019.2938229
  4. A. Alkhateeb, J. Mo, N. Gonzalez-Prelcic, and R. Heath, "MIMO precoding and combining solutions for millimeter-wave systems," IEEE Commun. Mag., vol. 52, no. 12, pp. 122-131, Dec. 2014. https://doi.org/10.1109/MCOM.2014.6979963
  5. W. Roh et al., "Millimeter-wave beamforming as an enabling technology for 5G cellular communications: Theoretical feasibility and prototype results," IEEE Commun. Mag., vol. 52, no. 2, pp. 106-113, Feb. 2014. https://doi.org/10.1109/MCOM.2014.6736750
  6. Z. Pi and F. Khan, "An introduction to millimeter-wave mobile broadband systems," IEEE Commun. Mag., vol. 49, no. 6, pp.101-107, June. 2011. https://doi.org/10.1109/MCOM.2011.5783993
  7. W. Zhang, X. Xia, Y. Fu and X. Bao, "Hybrid and full-digital beamforming in mmWave Massive MIMO systems: A comparison considering low-resolution ADCs," China Commun., vol. 16, no. 6, pp. 91-102, June 2019. https://doi.org/10.23919/j.cc.2019.06.008
  8. Q. Ding, Y. Deng, X. Gao and M. Liu, "Hybrid precoding for mmWave massive MIMO systems with different antenna arrays," China Commun., vol. 16, no. 10, pp. 45-55, Oct. 2019. https://doi.org/10.23919/jcc.2019.10.003
  9. O. El Ayach, S. Rajagopal, S. Abu-Surra, Z. Pi, and R. Heath, "Spatially sparse precoding in millimeter wave MIMO systems," IEEE Trans. Wire. Commun., vol. 13, no. 3, pp. 1499-1513, Mar. 2014. https://doi.org/10.1109/TWC.2014.011714.130846
  10. C. Chen, "An Iterative Hybrid Transceiver Design Algorithm for Millimeter Wave MIMO Systems," IEEE Wire. Commun. Let., vol. 4, no. 3, pp. 285-288, June 2015. https://doi.org/10.1109/LWC.2015.2409268
  11. X. Yu, J. Shen, J. Zhang and K. B. Letaief, "Alternating Minimization Algorithms for Hybrid Precoding in Millimeter Wave MIMO Systems," IEEE Journal of Selected Topics in Signal Processing, vol. 10, no. 3, pp. 485-500, April 2016. https://doi.org/10.1109/JSTSP.2016.2523903
  12. J. Jin, Y. R. Zheng, W. Chen and C. Xiao, "Hybrid Precoding for Millimeter Wave MIMO Systems: A Matrix Factorization Approach," IEEE Trans. Wire. Commun., vol. 17, no. 5, pp. 3327-3339, May 2018. https://doi.org/10.1109/twc.2018.2810072
  13. L. Liang, W. Xu and X. Dong, "Low-Complexity Hybrid Precoding in Massive Multiuser MIMO Systems," IEEE Wire. Commun. Lett., vol. 3, no. 6, pp. 653-656, Dec. 2014. https://doi.org/10.1109/LWC.2014.2363831
  14. W.U. Khan, X. Li, A. Ihsan, et al., "NOMA-enabled Optimization Framework for Next-generation Small-cell IoV Networks under Imperfect SIC Decoding," IEEE Trans. Intel. Trans. Sys., 2021.
  15. X. Gao, L. Dai, S. Han, C. I and R. W. Heath, "Energy-Efficient Hybrid Analog and Digital Precoding for MmWave MIMO Systems With Large Antenna Arrays," IEEE J. Sel.Areas Commun., vol. 34, no. 4, pp. 998-1009, April 2016. https://doi.org/10.1109/JSAC.2016.2549418
  16. J. Du, W. Xu, H. Shen, X. Dong and C. Zhao, "Hybrid Precoding Architecture for Massive Multiuser MIMO With Dissipation: Sub-Connected or Fully Connected Structures?," IEEE Trans. Wire. Commun., vol. 17, no. 8, pp. 5465-5479, Aug. 2018. https://doi.org/10.1109/twc.2018.2844207
  17. O. Alluhaibi, Q. Z. Ahmed, C. Pan and H. Zhu, "Capacity Maximization for Hybrid Digital-to-Analog Beamforming mm-Wave Systems," in Proc. of 2016 IEEE Global Communications Conference (GLOBECOM), Washington, DC, pp. 1-6, 2016.
  18. W.U. Khan, F. Jameel, X. Li, et al., "Joint Spectrum and Energy Optimization of NOMA-enabled Small-Cell Networks with QoS Guarantee," IEEE Trans. Vehic. Techn., vol. 70, no. 8, pp. 8337-8342, 2021. https://doi.org/10.1109/TVT.2021.3095955
  19. W. Ye, Z. Weixia, "Hybrid precoding design for millimetre wave systems with the partially-connected structure," IET Commun., vol.14, no.4, pp.561-567, 2020. https://doi.org/10.1049/iet-com.2019.0188
  20. X. Xue, Y. Wang, L. Yang, et al., "Energy-Efficient Hybrid Precoding for Massive MIMO mmWave Systems with A Fully-Adaptive-Connected Structure," IEEE Transactions on Communications, vol. 68, no. 6, pp. 3521-3535, June 2020. https://doi.org/10.1109/tcomm.2020.2979139
  21. V. V. Ratnam, A. F. Molisch, O. Y. Bursalioglu, and H. C. Papadopoulos, "Hybrid beamforming with selection for multiuser massive MIMO systems," IEEE Trans. Signal Process., vol. 66, no. 15, pp. 4105-4120, Aug. 2018. https://doi.org/10.1109/TSP.2018.2838557
  22. Y. Gao, H. Vinck, and T. Kaiser, "Massive MIMO antenna selection: Switching architectures, capacity bounds, and optimal antenna selection algorithms," IEEE Trans. Signal Process., vol. 66, no. 5, pp. 1346-1360, Mar. 2018. https://doi.org/10.1109/tsp.2017.2786220
  23. S. Park, A. Alkhateeb and R. W. Heath, "Dynamic Subarrays for Hybrid Precoding in Wideband mmWave MIMO Systems," IEEE Trans. Wirel. Commun., vol. 16, no. 5, pp. 2907-2920, May 2017. https://doi.org/10.1109/TWC.2017.2671869
  24. L. Yan, C. Han and J. Yuan, "A Dynamic Array-of-Subarrays Architecture and Hybrid Precoding Algorithms for Terahertz Wireless Communications," IEEE J. Sel.Areas Commun., vol. 38, no. 9, pp. 2041-2056, Sept. 2020. https://doi.org/10.1109/JSAC.2020.3000876
  25. H. Li, M. Li, Q. Liu and A. L. Swindlehurst, "Dynamic Hybrid Beamforming With Low-Resolution PSs for Wideband mmWave MIMO-OFDM Systems," IEEE J. Sel.Areas Commun., vol. 38, no. 9, pp. 2168-2181, Sept. 2020. https://doi.org/10.1109/JSAC.2020.3000878
  26. Z. Pi and F. Khan, "An introduction to millimeter-wave mobile broadband systems," IEEE Commun. Mag., vol. 49, no. 6, pp. 101-107, Jun. 2011. https://doi.org/10.1109/MCOM.2011.5783993
  27. H. S.F., I. C.L., X. Z.K., and R. C., "Large-scale antenna systems with hybrid precoding analog and digital beamforming for millimeter wave 5G," IEEE Commun. Mag., vol. 53, no. 1, pp. 186-194, Jan. 2015. https://doi.org/10.1109/MCOM.2015.7010533
  28. C.A. Balanis, Antenna Theory: Analysis and Design, Hoboken, NJ, USA: Wiley, 2012
  29. G. H. Golub and C. F. Van Loan, Matrix Computations, MD, USA: JHU Press, 2012.
  30. A. Bjorck, Numerical methods in matrix computations, Berlin, German: Springer, 2015.
  31. Chen J., Gan M., Zhu Q., et al., "Robust Standard Gradient Descent Algorithm for ARX Models Using Aitken Acceleration Technique," IEEE Trans. Cybern., (early access), 2021.
  32. T. X.Q., W. L.J., and F. W., Numerical Calculation, Beijing, China: Sci. Press, 2015.
  33. S. Boyd, L. Vandenberghe, "Numerical linear algebra background,". Accessed 12 Feb 2016. [Online]. Available: num-lin-alg2.pdf (stanford.edu)