International Journal of Antennas and Propagation
Volume 2019 (2019), Article ID 4196513, 12 pages
https://doi.org/10.1155/2019/4196513

Research Article

Energy-Efficient Hybrid Precoding Scheme Based on Antenna Selection Technology in Massive Multiple-Input Multiple-Output Systems

Jian Jun Ding and Jing Jiang

College of Communication and Information Engineering, Xi’an University of Posts and Telecommunications, Xi’an, Shanxi, China

Correspondence should be addressed to Jing Jiang; jiangjing@xupt.edu.cn

Received 31 August 2018; Accepted 4 March 2019; Published 2 May 2019

Guest Editor: Raed A. Abd-Alhameed

Copyright © 2019 Jian Jun Ding and Jing Jiang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Hybrid precoding is a promising technology for massive multiple-input multiple-output (MIMO) systems. It can reduce the number of radio frequency (RF) chains. However, the power consumption is still very high owing to the large-scale antenna array. In this paper, we propose an energy-efficient precoding scheme based on antenna selection technology. The precoding scheme greatly increases the energy efficiency (EE) of the system. In the first step, we derive an exact closed-form expression of EE. Meanwhile, we further study the relationship between the number of transmit antennas and EE on the basis of the exact closed-form expression of EE. We prove that there exists an optimal value. When the number of transmit antennas equals to the value, the EE of the system can reach the maximum by a proper hybrid precoding scheme. Then, we propose an antenna selection algorithm to select antennas from the transmit antennas. And the number of selected antennas equals to the optimal value. Subsequently, we design the analog precoder based on a codebook to maximize the equivalent channel gain. At last, we further improve the EE by baseband digital precoding. The precoding algorithm we proposed offers a compromise between spectral efficiency (SE) and EE in millimeter wave (mmWave) massive MIMO systems. Finally, simulation results validate our theoretical analysis and show that a substantial EE gain can be obtained over the precoding scheme we proposed without large performance loss.

1. Introduction

Millimeter wave band communication recently acquires more and more attention owing to its great advantages [1–4]. A beneficial feature of mmWave is that many antenna arrays can be packed into a small dimension owing to the small wavelength [5–11]. However, the problem is that the availability of multiantennas causes high interference between different users and high hardware complexity. Fortunately, precoding can eliminate interference between different users and reduce hardware complexity. Precoding is therefore more favored over mmWave MIMO communication systems. Simple linear precoding schemes, for example, zero forcing (ZF) and minimum mean square error precoding (MMSE), are virtually optimal. However, the digital processing in the MIMO system requires a dedicated base-band and radio frequency (RF) chain for every antenna element. Owing to the large amount of antennas deployed in mmWave massive MIMO systems, the costs are very huge and it is impossible to popularize in practice [12]. Researchers have widely studied on the reduction of the hardware cost. A hybrid precoding scheme comprising both digital and analog processing is capable of reducing the number of RF chains greatly. Therefore, the hybrid precoding scheme is widely used in mmWave MIMO systems. In consideration of a single-user scenario, a fully connected architecture-based hybrid precoding scheme is proposed in [13, 14], where each RF chain is connected to all antennas by analog phase shifters (APSs) and RF adders. The number of RF chains is greatly reduced in the fully connected architecture. However, in the fully connected architecture, the number of APSs is equal to the product of the number of RF chains and the number of antennas. The cost of hardware is still excessive owing to the large amount of antenna elements. Different from the fully connected architecture, a partially connected architecture-based hybrid precoding scheme is proposed in [15–19], where each RF chain only connected to an antenna subarray and the number of APSs equals to the number of antennas. Compared with the fully connected architecture, the number of phase converters and energy consumption is greatly reduced. In the partially connected precoding scheme, the mapping relationship between the antennas and RF chains is predetermined, but the channel condition is time-varying. Thus, the partially connected architecture cannot guarantee that the mapping relationship is optimal under different channel conditions. A dynamic subarray architecture-based hybrid precoding scheme is proposed in [20], where an adaptive antenna selection network is added between RF chains and antenna elements to enhance spectral efficiency (SE). All the previous work is aimed at improving the SE of the system. In the fully connected architecture-based and partially connected architecture-based architecture hybrid precoding schemes, all the transmitted antennas are activated. In this paper, we propose an energy-efficient precoding scheme with considerably reduced energy consumption and assume that not all transmitting antennas are activated. We jointly optimize the number of activated antennas, analog precoding matrix, and digital precoding matrix to maximize the of the system. Firstly, instead of designing the digital and analog precoder directly, we perform an antenna selection procedure before digital and analog precoding. Then, we further optimize the performance in terms of by digital and analog precoding.

The contributions existing in this paper are summarized as follows:

(1)An exact closed-form expression of was derived in this paper, and we derived an optimum number of antennas for maximizing the EE according to the exact closed-form expression of (2)In previous precoding schemes, all the antennas are activated. In this paper, not all the transmit antennas are activated (an antenna is activated, which means the antenna is used to transmit message). Based on channel state information and the exact closed-form expression of , we select antennas from the transmitting antennas to activate and the remaining antennas will be temporarily closed(3)In the prior hybrid precoding scheme, we design the analog precoder and the digital precoder according to the channel state information which was obtained by channel estimation to improve the sum rate and reduce the interference between different users. The hybrid precoding scheme we proposed differs from the prior hybrid precoding scheme. The transmitter first performs an antenna selection process according to the channel state information and then optimizes system performance in terms of by analog and digital precoding. Through the scheme, the number of antennas and power consumption is reduced, and the is improved(4)An energy-efficient hybrid precoding scheme for a single user in mmWave systems is developed in this paper. First, we calculate the optimal number of transmit antennas and then the antenna selection algorithm is used to select the subset of transmit antennas. Then, we use the analog precoding scheme to maximize the gain of the equivalent channel between BS and objective users. Furthermore, we use the digital precoding scheme to maximize the of the system

The rest of this paper is organized as follows: Section 2, Section 3, and Section 4 introduce the system model, channel model, and power consumption model of the system, respectively. In Section 5, we introduce the optimization problem and propose an antenna selection algorithm. In Section 6, we proposed a hybrid beamforming algorithm. Computer simulation results are shown in Section 7. Finally, conclusions are drawn in Section 8.

Notations: we use the following notation throughout the paper. denotes a matrix; is a vector; denotes a scalar; denotes the transpose of . denotes the conjugate of . denotes the inversion of . denotes the -th element of . We express as . is an identity matrix. The acronym denotes “subject to,” and denotes “independent and identically distributed.” denotes the mathematical expectation of , and denotes complex Gaussian distribution.

2. System Model

In this paper, we consider a downlink SU-MIMO mmWave system where a base station (BS) is equipped with antennas and RF chains to serve a single user which has RF chains and antennas. The number of RF chains satisfies to guarantee multistream transmission.

First, transmit data streams at the BS are processed by a digital precoder in the baseband and then processed by an analog precoder (RF precoder using analog circuitry). of size denotes the transmitting analog beam former. denotes the base band digital precoder satisfying , and denotes the total transmit power. Notably, the RF precoder can realize only phase changes because of phase-only control. For hybrid precoding systems as shown in Figure 1, the received signal vector of the objective user can be expressed by

where denotes the received signal vector. denotes the channel matrix with denoting the channel vector between the BS and the -th receiving antenna. denotes the transmitted signal vector, satisfying . And represents the number of data streams. is the vector of additive complex Gaussian noise with zero mean and variance . denotes the power of noise.

Figure 1: A mmWave single-user MIMO system.

3. Channel Model

Since mmWave channels have high free-space path loss, the mmWave propagation can be perfectly characterized by a clustered channel model [21] in this paper. The Saleh-Valenzuela model was used in this paper.

where the number of cluster and the number of rays can be expressed by and , and the gain of the -th ray in the -th cluster can be expressed as . We assumed that satisfies the complex Gaussian distribution and satisfies , and is a normalization factor, where and denotes the receive and transmit array response vectors, respectively. and denote the azimuth and elevation angles of arrival and departure, respectively.

In this paper, we adopt a uniform square planar array with antenna elements.

denotes the distance between antennas, and denotes the wavelength of the signal. satisfies , the same with . and are the antenna indices in the 2D plane.

4. Power Consumption Model

This section discusses the power consumption model of the downlink single-user MIMO mmWave system. The power consumption model of the system [22–24] can be denoted by

where denotes the inefficiency factor of the power amplifier (PA). denotes the power consumption of the transmitter, and . consists of the dynamic power consumption and static power consumption . can be denoted by

where denotes the power consumption of all RF chains. The number of APSs is ; thus, the power consumption of phase shifters is , where denotes the energy consumption of APS and denotes the power consumption of the RF chain.

5. Problem Formulation and Antenna Selection Algorithm

5.1. Problem Formulation

First, we consider the relationship between system capacity and signal-to-noise ratio. According to Shannon’s theorem, the relationship between system capacity and SNR can be expressed by the following formula:

Once the antenna set is determined, the corresponding channel matrix is determined. We can define the channel matrix as . Equation (6) can be denoted as

The objective of this paper is to select the transmit antenna subset according to channel state information and design the hybrid precoders to maximize the of the system. The problem can be formulated as

where denotes the selected antenna subset. The elements of are the index of selected antennas (e.g., denotes that the selected antennas are the first, second, and fourth antennas of the BS). It is worth pointing out that the constraint in equation (8) is nonconvex, because of the coupling between the analog precoding matrix and digital precoding matrix. This makes equation (8) to be a difficult problem to solve. Fortunately, we can decouple the optimization problem into two convex optimization problems. First, we focus on the antenna selection process. We select the optimal subset of the antennas from the transmit antennas according to channel matrix . Then, with the channel matrix between the selected antenna and objective user, the hybrid precoding scheme can be adopted to improve the of the system. Specially, we adopt the analog precoding to maximize the gain of the equivalent channel and digital precoding to maximize the .

At first, we would briefly introduce the of the system without consideration of hybrid precoding. The is defined as the throughput divided by the total energy consumed, which can be denoted by

where denotes the mathematical expectation of , and denotes the mutual information. If the power consumption of the receiver is taken into account, the circuit power consumption model can be expressed:

where and denote the number of transmitters and receivers, respectively, and , , , , , , and are the power consumption of the digital analog converter, the mixer, the filter at the transmit side, the synthesizer, the low-noise amplifier, the intermediate-frequency amplifier, the filters at the receiver side, and the analog digital converter, respectively. We set

where is a very small value much less than one.

Lemma 1. Suppose there exists a communication system. transmit antennas were equipped by the BS, and antennas are selected from these transmitter antennas; the distribution of mutual information is approximately [25] given by

Thus, . So, the closed-form expression of the EE becomes

Next, we will analyze the effects of the number of transmitted antennas on the EE.

Lemma 2. If the power consumption of the mmWave MIMO system can be modeled as the addition of the transmit power and the circuit power consumption , increases at first and then decreases as the number of selected transmit antennas increases.

Proof 1. can be expressed as a function of :

We can take the first derivative of

setting . Because the denominator of (17) is positive, the denominator of (17) has no effect on the monotonicity of function (17). Setting the numerator of (17) as ,

Then, we take the first derivative of again:

It is obvious that is a negative value. So, is a decreasing function of ; thus, :

and

We do not know whether is a positive or negative value. However, equation (21) can also be viewed as a function of . We set . Then, we take the first derivative of .

We take the first derivative of (22) again in respective of , we can conclude that

As mentioned earlier in this paper, is a very small number much less than one, so . Then, , so is an increasing function of . Besides,

so . So, increases first and then decreases with the increase in . Assume that satisfies . With mathematical knowledge, we can conclude that

(1)if , the of the system increases with the increase in (2)if , the of the system decreases with the increase in

That is to say, there exists an optimal number of antenna that can maximize the of the system.

5.1.1. Discussion

Based on the system model and power consumption model that we proposed, increases at first and then decreases with the increase in the number of antennas. This is because with the increase in the number of selected antennas, the power consumption increases linearly. However, this is not true for the sum rate of the system. When the number of antennas increases to a certain value, the increase in the sum rate is not enough to compensate for the increase in the energy consumption of the system, so the of the system will reduce.

Denote the selected transmit antenna subset as , the number of active transmits is where denotes the cardinality of . could also be also expressed as

where is the channel matrix of size between the transmit antenna and the receive antennas. The objective of antenna selection can be expressed as

where (26) is a nonconvex optimization problem [26–29]. It is often intractable to directly solve, because it is a joint optimization problem over three variables. We propose a suboptimal algorithm to decouple the nonconvex optimization problem into a convex optimization problem that can be solved directly. First, just as it is proved above, there exists an optimal value . When the number of transmitter antennas equals to the value, the of the system can achieve an optimal value. Thus, the proposed algorithm will be divided into three parts.

At first, as we know, there exists an optimal value that maximizes the . Thus, we can select antennas from transmit antennas to optimize the of the system before hybrid precoding.

Then, an analog precoder is designed to maximize the gains of the equivalent channel, where the equivalent channel is defined as . At this step, it is assumed that the digital precoding matrix is fixed.

The design of digital precoder will become an easily solved optimization problem with a single variable after the analog precoding. The digital precoding matrix will be designed based on the maximized criterion.

The proposed antenna select algorithm and hybrid precoding scheme will be discussed in the next section.

Input: the number of transmit antenna , the power consumption of digital analog, converter, filter, synthesizer, low-noise amplifier, the intermediate-frequency amplifier, and the analog digital amplifier: , , , , , , the transmit power p, channel matrix .

Initialization:

    1: for do

    2:      Calculate:

    3:      Update ;

    4:      ;

    5: end for

    6: Sort in a descend order as:

    7: Setting

    8: Setting

    9: Sort in a descend order as

    10:

    11: Output:

Algorithm 1: The proposed antenna selection scheme.

5.2. Antenna Selection Algorithm

In this section, a low-complexity antenna select algorithm is proposed which enables the BS to reduce the number of transmit antennas and without large performance loss. The antenna selection algorithm is described as below.

As discussed above, increases first and then decreases with the increase in transmit antenna in mmWave massive MIMO systems. So, we can select antennas from transmit antennas to optimize the of the system before hybrid precoding at first.

The pseudocode of the proposed antenna selection algorithm is summarized in Algorithm 1, which can be explained as follows. At the beginning, we formulate the of the system as the sum rate divided by energy consumption and set the denoted in (13). During the antenna selection, in step 1, step 2, and step 3, we calculate the of the mmWave MIMO system with different numbers of antennas. In step 4, we find the optimal number of antennas by comparing under different conditions. Then, in step 6, step 7, and step 8, we choose optimal antennas from transmitted antennas by comparing the Frobenius norm of the channel vector of each antenna. Finally, the antennas are used to transmit information. Note that the antenna selection algorithm enables the system to adaptively choose the transmit antenna under different channel conditions to maximize the of the system.

6. The Design of the Hybrid Precoder

In this section, we will discuss the design of the hybrid precoder for the mmWave massive MIMO system under the condition that the transmit antenna subset has been obtained.

As the antenna selection process effectively reduces the number of transmit antennas, the power consumption of the system is reduced. After the antenna subset has been obtained, analog precoding and base band digital precoding will be utilized to improve the of the system. The optimization problem (26) can be written as

where denotes the sum rate of the proposed system and denotes the energy consumption of the whole system.

At first, an analog precoding algorithm is proposed to maximize the power of the received signal for the objective user [30–32]. Further, a digital precoder will be designed to maximize the of the system. However, it is difficult to solve problem (27) because of the coupling between two matrix variables: and . We solve this problem by decoupling the joint optimization problem into two convex optimization problems.

Firstly, an analog precoder is designed that is aimed at maximizing the equivalent channel gain and assuming that the base-band precoding matrix is fixed at this stage.

After the analog precoding matrix is obtained, the original optimization problem (27) becomes an easily solved convex optimization problem with one matrix variable , which can be expressed as

Then, the optimal baseband precoding matrix can be obtained by solving (29).

6.1. The Design of the Analog Precoder

In this section, we will discuss the optimal analog precoding matrix , where ; the equivalent channel is defined as . We rewrite the receive signal of the objective user

The gain of the equivalent channel can be expressed as

We aim to find an analog precoding matrix to maximize the equivalent channel gain based on codebook . is a predefined RF beamforming codebook. can be specified as a quantized matrix . For , it takes the value 0,1,2,3,…,, and for , it takes the value 0,1,2,3,…,. and denote the quantized precision of azimuth and elevation angles, respectively. The most optimal precoding matrix is designed based on the following metric.

The pseudocode of the proposed analog precoding algorithm is summarized in Algorithm 2.

Input: Channel matrix H, analog beamforming codebook , the number of RF chains

Initialization:

    1: for do

    2:      while do

    3:            User feed back the index of codebook vector ri to BS. BS finds corresponding vector from codebook and setting as the ith column of FRF, that is

    4:            

    5:      end while

    6:      

    7: end for

Algorithm 2: Analog precoding scheme.

At the beginning, we set , and we search in the codebook to find out an optimal vector which has the largest inner product with channel matrix and make it as the first column of the analog precoding matrix .

Then, is updated by , and the procedure above is repeated to find optimal vector and is set as the ()-th column of analog precoding matrix .

The procedure is repeated until . Finally, could be obtained, i.e., .

6.2. Digital Precoding Design to Improve EE

In this section, we discuss the design of digital beamformer to further improve the of the system. After the analog precoding matrix is obtained, the original nonconvex problem is conveyed to a univariate optimization [33] problem which is given by

There exist many approaches to solve this classical fractional programming problem [34]. At first, we transform it into an equivalent parametric programming problem, , where can be expressed as

, where satisfies and denotes the max sum rate of the system. The optimal digital precoding matrix which is the solution of (35) is also the optimal solution of (33). The general method to solve (35) is the half-interval search method. It is summarized as in Algorithm 3.

Initialization: , and

    1: while do

    2:        

    3:        Solve the problem (35) for a given ,

    4:        update and .

    5:        If let . Otherwise, .

    6: end while

    7:

Algorithm 3: Half-interval search method.

In Algorithm 3, the key step is solving (35). In this section, we will discuss how to solve (35). Problem (35) can be expressed as

According to the knowledge of the determinant, it is obvious that ; (35) could be conveyed into

where and ; it is obvious that for any matrix which is nonnegative, . Thus, only if is a diagonal matrix will (37) achieve its maximum value. Then, (38) will further expressed as

The problem in (38) was a convex optimization problem, and its general solution method is to convert it into an unconstrained optimization problem with the Lagrange multiplier method [35, 36]. Equation (38) could be converted into (39) by the Lagrange multiplier method:

By solving the KKT conditions of (39), the optimal diagonal entries of are obtained for a given via the water-filling algorithm [37, 38].

is a decreasing function of ; thus, could be found by a one-dimensional search algorithm. We proposed a binary water-filling algorithm in this paper. The binary water-filling algorithm is summarized in Algorithm 4.

Initialization: ,

    1: for do

    2:      

    3:      calculate based on (39), (40)

    4:      update and

    5:      if ,

    6:      otherwise

    7:      

    8: end for

    9:

    10: Calculate according to (40)

    11:

Algorithm 4: Binary iterative water-filling algorithm.

7. Simulation Results

In this section, we present simulation results to evaluate the performance of the proposed algorithm in terms of and the sum rate. Specially, the simulation parameters are summarized in Table 1. In this section, we first compare the performance of different algorithms in terms of with different antenna array sizes. Then, we compare the performance of the proposed algorithm with a conventional algorithm under different power constraints in terms of . At last, we compare the performance of different precoding schemes in terms of the sum rate.

Table 1: Simulation parameters.

Parameters
Value

Number of antennas
0~180

Channel model
IEEE 802.11ad living room

Number of users
1

Carrier frequency (GHz)
60

Bandwidth
2.56

Power amplifier efficiency (%)
26

Modulation coding scheme
16PSK, 1/2 turbo codes

Simulation frames
20000

Number of RF chains
2~20

Figure 2 compares the performance between the energy-efficient algorithm and spectral efficient algorithm under different SNRs in terms of SE. As Figure 2 shows, with the increase in SNR, the system performance of different algorithms increases, but the performance of the pure digital precoding system is optimal. The performance of the spectral efficient hybrid precoding scheme is close to that of the digital precoding scheme, and the performance of the energy-efficient algorithm is lower than that of the spectral efficient algorithm in terms of SE. This is because the energy-efficient algorithm sacrifices part of the spectrum efficiency to improve the energy efficiency of the system. And, for the same algorithm (energy-efficient algorithm/spectral efficient algorithm), using a different number of RFs will lead to a different performance in terms of spectral efficiency. As can be seen from the figure, under the same SNR condition, the higher the number of radio frequency chains, the better the system performance.

Figure 2: SE comparison of different precoding algorithms under different SNRs.

Figure 3 compares the performance between the energy-efficient algorithm and spectral efficient algorithm under different transmit power constraints. As Figure 3 shows, at low power constraint, the energy-efficient performance of the proposed algorithm is exactly close to that of the spectral efficient algorithm. As the transmit power increases, the of the proposed algorithm remains a constant but shows a decline in the spectral efficient algorithm. The reason is that in the spectral efficient algorithm, the increase in the sum rate cannot keep path with the increase in power consumption.

Figure 3: comparison of different precoding algorithms under different power constraints.

Figure 4 compares the performance between the energy-efficient algorithm and spectral efficient algorithm under different antenna array sizes in terms of . As Figure 4 shows, there exists an optimal antenna array size to maximize the of the system.

Figure 4: comparison of different precoding algorithms under different transmit antennas.

The increases at first and then decreases as the number of transmit antennas increases. The more power constraint there is, the larger the antenna array size is. It is reasonable, because once the number of transmit antennas increases over a certain value, the increase in cannot compensate for the negative influence of power consumption increase .

Figure 5 compares the energy efficiency achieved under a given optimal number of antennas and transmits power but a different number of RF chains. As Figure 5 shows, under the condition that the optimal number of transmitting antennas is 80, with the increase in the number of RF chains, the of the system increases first and then decreases. It can be seen from the figure that the optimum number is 8. Meanwhile, it can be seen from Figure 5 that when the number of transmitting antennas is the same, the optimal number of the RF chain is the same. At this time, the larger the transmitting power of the system, the higher the that the system can achieve.

Figure 5: comparison of the energy-efficient precoding scheme under different chains.

Figure 6 compares the achieved under a given optimal number of antennas and transmit power but a different number of RF chains. As Figure 5 shows, when the optimal number of transmitting antennas is 80, the of the system increases at first and then is kept constant with the increase in the number of radio frequency chains. As can be seen from the figure, when the number of RF chains increases from 1 to 10 under the condition that and , the of the system increases. When the number of chains is greater than 10, the of the system will be constant. However, with the increase in the number of , the power consumption of the system will increase, which will lead to the decrease in the of the system. This is consistent with the result of Figure 6. In Figure 6, we can conclude that the of the system reaches its maximum when the number of chains is 8 under the condition that and . Meanwhile, the larger the number of transmitting antennas , the greater the maximum the system can achieve.

Figure 6: SE comparison of the energy-efficient precoding scheme under different radio frequency chains.

8. Conclusions

In this paper, an energy-efficient hybrid precoding scheme for single-user mmWave massive MIMO systems was considered. First, we calculate the optimal number of transmit antennas according to the channel state information and then design an antenna selection algorithm to select the antenna to maximize the channel gain. Finally, the of the system is optimized by digital precoding. Numerical results show that the near-optimal hybrid precoding algorithm achieved a near-optimal performance. Simulation results show that the algorithm we proposed can achieve higher than the traditional precoding scheme can. In future work, we will study on the joint optimization of transmit antennas, receive antennas, and precoder to improve the of the system.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This paper is sponsored by the National Natural Science Foundation of China under Grant 61871321, National Science and Technology Major Project under Grant 2016ZX03001016, and Innovation Team Project of Shaanxi Province under Grant 2017KCT-30-02.

References

D. W. Bliss, A. M. Chan, and N. B. Chang, “Mimo wireless communication channel phenomenology,” IEEE Transactions on Antennas and Propagation, vol. 52, no. 8, pp. 2073–2082, 2004.

K.-K. Wong, R. D. Murch, and K. B. Letaief, “Performance enhancement of multiuser mimo wireless communication systems,” IEEE Transactions on Communications, vol. 50, no. 12, pp. 1960–1970, 2002.

T. S. Rappaport, G. R. MacCartney, M. K. Samimi, and S. Sun, “Wide-band millimeter-wave propagation measurements and channel models for future wireless communication system design,” IEEE Transactions on Communications, vol. 63, no. 9, pp. 3029–3056, 2015.

C. Rusu, R. Mendez-Rial, N. Gonzalez-Prelcic, and R. W. Heath, “Low complexity hybrid precoding strategies for millimeter wave communication systems,” IEEE Transactions on Wireless Communications, vol. 15, no. 12, pp. 8380–8393, 2016.

T. E. Bogale and L. B. le, “Massive mimo and mmwave for 5g wireless hetnet: potential benefits and challenges,” IEEE Vehicular Technology Magazine, vol. 11, no. 1, pp. 64–75, 2016.

T. Halinen, A. A. Dowhuszko, and J. Hamalainen, “Performance of distributed beamforming for dense relay deployments in the presence of limited feedback information,” EURASIP Journal on Advances in Signal Processing, vol. 2013, no. 1, 2013.

J. Zhang, X. Huang, V. Dyadyuk, and Y. Guo, “Massive hybrid antenna array for millimeter-wave cellular communications,” IEEE Wireless Communications, vol. 22, no. 1, pp. 79–87, 2015.

J. Singh and S. Ramakrishna, “On the feasibility of beamforming in millimeter wave communication systems with multiple antenna arrays,” IEEE Transactions on Wireless Communications, vol. 14, no. 5, pp. 3802–3808, 2014.

X. Gao, L. Dai, S. Han, C. L. I, and R. W. Heath, “Energy-efficient hybrid analog and digital precoding for mmwave mimo systems with large antenna arrays,” IEEE Journal on Selected Areas in Communications, vol. 34, no. 4, pp. 998–1009, 2016.

D. W. Browne, M. Manteghi, M. P. Fitz, and Y. Rahmat-Samii, “Experiments with compact antenna arrays for mimo radio communications,” IEEE Transactions on Antennas and Propagation, vol. 54, no. 11, pp. 3239–3250, 2006.

A. M. Sayeed and V. Raghavan, “Maximizing mimo capacity in sparse multipath with reconfigurable antenna arrays,” IEEE Journal of Selected Topics in Signal Processing, vol. 1, no. 1, pp. 156–166, 2007.

Z. Pi and F. Khan, “An introduction to millimeter-wave mobile broad-band systems,” IEEE Communications Magazine, vol. 49, no. 6, pp. 101–107, 2011.

X. Zhu, Z. Wang, L. Dai, and Q. Wang, “Adaptive hybrid precoding for multi-user massive mimo,” IEEE Communications Letters, vol. 20, no. 4, pp. 776–779, 2016.

P. Zhao and Z. Wang, “Hybrid precoding for millimeter wave communications with fully connected subarrays,” IEEE Communications Letters, vol. 22, no. 10, pp. 2160–2163, 2018.

C. C. Hu and J. H. Zhang, “Hybrid precoding design for adaptive subconnected structures in millimeter-wave mimo systems,” IEEE Systems Journal, vol. 13, no. 1, pp. 137–146, 2019.

T. Huang, M. Linhua, H. Xing, H. Shaocheng, S. Kangning, and L. Shiping, “Practical hybrid precoding algorithm for millimeter wave massive MIMO,” Journal of Electronics & Information Technology, vol. 39, no. 8, pp. 1788–1795, 2017.

M. Majidzadeh, A. Moilanen, N. Tervo, H. Pennanen, A. Tolli, and M. Latva-Aho, “Hybrid beamforming for single-user mimo with partially connected RF architecture,” in 2017 European Conference on Networks and Communications (EuCNC), Oulu, Finland, 2017.

M. Hanif, H. C. Yang, G. Boudreau, E. Sich, and H. Seyedmehdi, “Practical hybrid precoding for multi-user massive mimo systems,” in 2016 IEEE Canadian Conference on Electrical and Computer Engineering (CCECE), Vancouver, BC, Canada, 2016.

N. Li, Z. Wei, H. Yang et al., “Hybrid Precoding for mmWave Massive MIMO Systems With Partially Connected Structure,” IEEE Access, vol. 5, no. 1, pp. 15142–15151, 2017.

S. Park, A. Alkhateeb, and R. W. Heath, “Dynamic subarrays for hybrid precoding in wideband mmWave MIMO systems,” IEEE Transactions on Wireless Communications, vol. 16, no. 5, pp. 2907–2920, 2017.

Y. Kai and B. E. Ottersten, “Models for MIMO propagation channels: a review,” Wireless Communications and Mobile Computing, vol. 2, no. 7, 666 pages, 2002.

F. Heliot, M. A. Imran, and R. Tafazolli, “On the energy efficiency gain of mimo communication under various power consumption models,” in Future Network & Mobile Summit, Warsaw, Poland, 2011.

K. Wang, T. Nirmalathas, C. Lim, and S. Skafidas, “Power-consumption comparison of high-speed indoor millimeter-wave 60GHz communication system and optical wireless systems,” in Asia Pactific Microwave Photontics Conference in Japan, Kyoto, 2012.

C. Ma, C. Pan, Y. Pan, X. Han, and C. Ming, “Energy-efficient hybrid precoding for millimeter wave MIMO systems,” in International Conference on Communications & Networking in China, Shanghai, China, 2016.

H. Li, L. Song, and M. Debbah, “Energy efficiency of large-scale multiple antenna systems with transmit antenna selection,” IEEE Transactions on Communications, vol. 62, no. 2, pp. 638–647, 2014.

A. Moudafi, “A proximal iterative approach to a non-convex optimization problem,” Nonlinear Analysis: Theory, Methods & Applications, vol. 72, no. 2, pp. 704–709, 2010.

X. Y. Wang, X. W. Jiang, and J. X. Dai, “Homotopy method for nonconvex optimization problem,” Journal of Jilin University, vol. 52, no. 2, pp. 273–276, 2014.

Y. Bian, B. Mirzasoleiman, J. M. Buhmann, and A. Krause, “Guaranteed non-convex optimization: submodular maximization over continuous domains,” 2016, arXiv: 1606.05615.

P. Jain and P. Kar, “Non-convex optimization for machine learning,” Foundations and Trends® in Machine Learning, vol. 10, no. 3-4, pp. 142–336, 2017.

J. Jing, C. Xiaoxue, and X. Yongbin, “Energy-efficiency based downlink multi-user hybrid beamforming for millimeter wave massive mimo system,” The Journal of China Universities of Posts and Telecommunications, vol. 23, no. 4, pp. 53–62, 2016.

J. Jing, K. Deting, L. Guangyue, W. Junxuan, and Y. Wujun, “User scheduling algorithm for mmwave FDMA massive MU-MIMO system with hybrid beamforming,” in Proceedings of the 8th International Conference on Signal Processing Systems - ICSPS 2016, pp. 194–198, Auckland, New Zealand, 2016.

T. E. Bogale and L. B. Le, “Pilot optimization and channel estimation for multiuser massive mimo systems,” in 2014 48th Annual Conference on Information Sciences and Systems (CISS), Princeton, NJ, USA, 2014.

A. Zilinskas and G. Gimbutiene, “On one-step worst-case optimal trisection in univariate bi-objective Lipschitz optimization,” Communications in Nonlinear Science and Numerical Simulation, vol. 35, pp. 123–136, 2016.

J. P. Crouzeix and J. A. Ferland, “Algorithms for generalized fractional programming,” Mathematical Programming, vol. 52, no. 1-3, pp. 191–207, 1991.

D. P. Bertsekas, Constrained Optimization and Lagrange Multiplier Methods, Academic Press, 1982.

G. Haeser and M. L. Schuverdt, “On approximate KKT condition and its extension to continuous variational inequalities,” Journal of Optimization Theory and Applications, vol. 149, no. 3, pp. 528–539, 2011.

M. Kobayashi and G. Caire, “An iterative water-filling algorithm for maximum weighted sum-rate of Gaussian MIMO-BC,” IEEE Journal on Selected Areas in Communications, vol. 24, no. 8, pp. 1640–1646, 2006.

R. S. Prabhu and B. Daneshrad, “An energy-efficient water-filling algorithm for OFDM systems,” in 2010 IEEE International Conference on Communications, Cape Town, South Africa, 2010.