Outage-Constrained Beamforming for Two-Tier Massive MIMO Downlink with Pilot Reuse

Massive multiple-input multiple-output (MIMO) systems and small cell networks are both regarded as promising candidates to meet the exponential growth of mobile data traffic for the next generation (5G) wireless communications. Hence, a new kind of multitier networks which combine massive MIMO macro cells with a secondary tier of small cells is proposed to resolve the contradiction of large network coverage and high data rate. In such multitier networks, it is inevitable to allocate nonorthogonal uplink pilot sequences to user equipment (UE) due to the large number of users. We propose a pilot reuse scheme by exploiting the uniquearchitectureofthisnetworksandanalysethespecialmixedchannelstateinformation(CSI)yieldedbythepilotreusescheme. BasedonthemixedCSI,weformulateadownlinktransmitbeamformingproblemofminimizingthetotalpowerconsumption whilesatisfyingthequalityofservice(QoS)requirementswithoutageconstraints.Afterdecomposingtheoriginalprobleminto simplersubproblems,weprovideanefficientalgorithmtocombinethesesubproblemsandsolvethemiterativelyforgeneratingthe beamformingvectors.MonteCarlosimulationsshowthattheaveragepowerconsumptionoftheproposedpilotreuseschemeand itsassociatedbeamformingalgorithmisclosetothatoftheperfectCSIcase.


Introduction
With the advent of the fifth generation (5G) cellular wireless communications, an ever-increasing demand of substantially higher throughput is a quite pressing task laid in front of us and drives the researchers to find new technology for wireless communications.Massive MIMO is a potential technology for meeting this demand and also taken as an attractive solution for 5G systems [1].By implementing a large number of antennas at the BS, massive MIMO systems offer a high spatial and multiplexing resolution which can drastically improve the communication systems' performance in terms of data rate and reliability [2].However, a critical issue for such massive MIMO systems is the excessive energy cost for their large number of antennas.How to provide higher data rates with lower energy consumption for such massive MIMO systems is a critical problem in front of us.
An innovative solution to this problem is a densified network with a multitier network architecture as stated in [3].Sometimes, such multilayer and multiarchitecture networks are also called heterogeneous networks (HetNets).In [4], the authors have made further discussions about two-tier massive MIMO networks.The fundamental architecture of such networks is based on a deployment of a macro cell with very large antenna arrays in combination with a secondary tier of small cells (SCs) with a few antennas each.Macro cells are deployed for the coverage of large areas and capable of handling low data traffic or the users with a relative high mobility, and SCs with a reduced coverage range of tens of meters are designed for providing localized higher data rate communications.
However, there still exist some implicit obstacles to put such two-tier massive MIMO networks into practice.One problem is how to fulfill the channel estimation.In time division duplexing (TDD) cellular systems, channel estimation is obtained via the uplink training based on the channel reciprocity property.In other words, each user in the cell would be assigned a specific pilot sequence and these pilot sequences would be transmitted to the base station (BS) via the uplink.Acquiring perfect CSI inherently asks for 2 International Journal of Antennas and Propagation sufficient numbers of orthogonal pilot sequences, which may not be possible for massive MIMO cellular systems.Hence, the pilot reuse is inevitable in massive MIMO systems, and the imperfect uplink training due to reusing the same pilot sequences tends to be a critical problem (referred to as "pilot contamination" problem [5]).
It is critical to design an appropriate pilot reuse scheme for mitigating the interference due to the pilot contamination.In [6][7][8], a time-shifted pilot reuse scheme was proposed to combat pilot contamination by symmetrically rearranging the uplink pilot transmission order for different cells and the system performance for such scheme with zeroforcing beamforming and a large number of BS antennas is also studied.Reference [9] studied the optimal pilot reuse factor for sum-rate maximization in massive MIMO systems.In [10], the authors proposed a pilot reuse scheme in homogenous multicell networks based on the degree of the spatial orthogonality in greedy fashion.In [11], the authors partitioned each cell into 3 sectors and assigned the reused pilot sequences in a symmetrical way in order to perform pilot contamination precoding (PCP).Recently, a novel pilot reuse scheme is proposed in [12] which exploits the channel spatial localization property to reduce the number of orthogonal pilots for uplink channel estimation in single cell massive MIMO systems.However, the pilot reuse scheme in [12] requires restrictive assumptions such as high channel spatial correlation and uniform linear array at the BS.In summary, all the works mentioned above focused on either homogenous multicell networks or single cell systems.As a result, they cannot be directly applied to (or are not optimized for) HetNets with asymmetrical system architecture.
Another problem is how to tackle the intercell interference.With limited spectral resources, a cochannel deployment of macro cell and secondary tier SCs is the only viable solution, and this in turn requires a sophisticated interference management scheme across the tiers.The crosstier interference is even thought to be one of the bottlenecks for designing high-performance HetNets [13].Now, a tough problem of proposing a high-performance and low-complexity interference management scheme to guarantee the system's QoS requirements with limited uplink pilot resources is in front of us.
For this problem, robust transmit beamforming methods which take the CSI errors into consideration are in need for the reason that "pilot contamination" causes imperfect CSI.In [14,15], a sum-rate maximization problem with imperfect CSI was investigated.However, these methods mainly focused on the worst case approach under the assumption of normbounded uncertainty which is not suitable for the considered system.Furthermore, a more efficient conservative formulation is presented which involves solving a semidefinite programming (SDP) in [16].By using the Bernstein-type inequality, a chance constrained beamforming problem in cognitive radio networks is also provided in [17].The kernel of these methods is to transform the original problem to a tractable SDP problem conservatively.However, the huge computation complexity of SDP, ( 6.5 ), makes these methods unattractive in massive MIMO systems [18].Compared with the previous works, we propose a pilot reuse scheme by exploiting the unique architecture of this HetNet.A "mixed CSI" is achieved at the BS due to our pilot reuse scheme.Based on the mixed CSI, we formulate the energy-efficient beamforming problem as a chance constraint programming and decompose it into simpler subproblems.Then, an efficient iterative algorithm is provided to combine these subproblems and solve them iteratively for obtaining the beamforming vectors.
The rest of this paper is organized as follows.In Section 2, the system model is provided.How to compress the number of the uplink pilot sequences and how to design the energyefficient beamformer with the mixed CSI are presented and solved in Section 3. Section 4 presents the numerical results and Section 5 concludes this paper.

System Model
We consider the downlink of a TDD orthogonal frequency division multiplexing (OFDM) system.Different from the traditional systems with single base station, one macro base station (MBS) and  low range small cell base stations (SBSs) separately serve their intended single antenna users in the same frequency band shown as in Figure 1.The MBS has  BS antennas and each SBS has  SC antennas, where  BS is assumed to be a large number which is known as massive MIMO.We consider that the total number of the served users is ; since the total number of transmit antennas is  =  BS +  SC , the downlink channel propagation matrix between the  antennas and  users can be described by a  ×  matrix, H.The base stations can exploit channel reciprocity to obtain channel state information at the transmitter (CSIT) based on uplink training.Here we denote the macro cell as cell 0 and the th SC is denoted as cell .The channel of the th user is represented as h  = [h 0, h 1, ⋅ ⋅ ⋅ h , ] ∈ C 1× , where h 0, ∈ C 1× BS and h , ∈ C 1× SC ,  = 1, . . ., .Let ℎ ,, ∈ C, ∀, , , denote the entries of h , , ∀, , and ℎ ,, describes the channel between the th antenna in cell  and the th user.We consider a composite fading channel; that is, ℎ ,, = √ ,,  ,, , where  ,, is the path gain between antenna  in cell  and user , and  ,, ∈ C is the small scale fading with i.i.d.∼ CN(0, 1) entry.Here, we assume a time block fading model.Thus path fading vectors  ,, , ∀, , , stay constant during each coherence interval and these vectors are assumed to be independent in different coherence blocks.Since all the transmit antennas of MBS are considered to be center collected, h 0, ∼ CN(0,  0, I) is satisfied, where  0, , ∀ = 1, . . ., , are independent random variables describing the path loss fading.
We divide the  users into two groups represented by small-cell user equipment (SUE) and macro cell user equipment (MUE).If user  is located in the coverage of the th SC and register at the th SBS for a certain high quality communication service, then user  is a SUE of SBS .The other group is the MUE who do not register at any SC and communicate with the MBS for a common quality communication service.To specify which cell it belongs to, the SUE  located in the th SBS can also be denoted as user   .If user  is a MUE, then it can be also denoted as  0 similarly.For convenience, we suppose that each SC has   registered SUE.Thus, there are    ≤  users registered at all the SCs totally.
We assume that the MBS and SBSs are connected to a backhaul network with limited capacity which enables interference coordination.When the backhaul is capacity limited, the coordination strategy with no data and only limited CSI sharing is preferred via the backhaul [19].Note that each user is considered to be served only by its associated BS in our system; that is, MUE are served by MBS and each SUE is served by its associated SBS.The information symbols for user   from its associated th BS can be denoted as    (or   ).These information symbols are independently precoded by the beamforming vectors, w 0, ∈ C  BS or w , ∈ C  SC , where w , is the beamforming vector at the th BS for the th user.Then the received signal at user  is (note that w , = 0 satisfies when BS  does not serve the th user) where the term   (or    ) is the circularly symmetric complex Gaussian receiver noise which is an i.i.d.CN(0,  2  ) random variable.The signal   is originated from independent Gaussian codebooks and constrained to have a total expected power of one; that is, For a given set of optimized beamforming vectors w , , ∀, , the signal-to-interference and noise ratio (SINR) of the th user can be expressed as Consequently, the information rate of user  is calculated by We model the power consumption of the system as an addition of the transmit power and the circuit power like that in [20] as follows: where  tot denotes the total energy consumption of the system.  is the dynamic power consumption known as the aggregation of the emitted power and   is modeled as the static power consumption cost by the BS components attached to each antenna, for example, A/D (D/A) converters, filters, and mixers, where   and   can be calculated by where   ∈ (0, 1] represents the efficiency of transmitter power amplifiers in the th cell,   is the power consumption in the circuits of each antenna in the th cell, and ‖ ⋅ ‖  stands for the Frobenius norm.

Pilot Reuse and Outage-Constrained Beamforming
In this section, we will introduce our pilot reuse scheme and analyze the special mixed CSI yielded by this scheme.A low complexity algorithm for downlink beamforming is also provided based on the mixed CSI.

Pilot Reuse and System
Operating.In such two-tier HetNet, at least  orthogonal uplink pilot sequences are needed to achieve the global perfect CSI [21] (in this paper, we ignore the issues related to channel estimation and simply assume that the perfect CSI between one UE and certain BS is available at the BS if the UE's allocated pilot sequence is orthogonal to that of the other UEs' .Since pilot power and data power could be assigned individually [22,23], this is an appropriate assumption even when the users' SINR is varying).In practice, the short channel coherence intervals do not allow for such long training sequences [24].
Hence, in such multi-BS networks, nonorthogonal training sequences must be utilized.The pilot contamination problem was studied in [5] which found that pilot contamination would result in corruption of precoding matrices and cause significant intercell interference (ICI).The ICI brings about the saturation effect; that is, the system throughput does not grow with the increasing of the number of BS antennas.
Compressing the number of pilot sequences to a precise value and allocating these pilots in an appropriate way is a hard problem.To the best of the authors' knowledge, there has been no systematic method for this task currently.While, by International Journal of Antennas and Propagation Figure 2: A kind of TDD-based frame structure.The pilot reuse scheme for the two-tier HetNet is presented in (a).The pilot scheme for achieving global perfect CSI is presented in (b).
exploiting the special architecture of the considered HetNet, a simplified scheme for this pilot allocation problem can be obtained.In such HetNet, the SCs are geometrically separated with each other and the coverage regions of different SCs are often nonoverlapped.As a result, the users in a SC hardly receive signal from other SBSs, which means that the interference among the SCs can be ignored.Thus, the crosstier interference is the bottleneck to improve the throughput of such HetNet.Now, we consider how to compress the number of pilot sequences and allocate these pilots to different users in such HetNet.On one hand, the MUE may walk around the macro cell with a relative high mobility; instantaneous CSI is crucial for the MBS to serve these MUE; on the other hand, each SC would have strong channel gain to its associated SUE, instantaneous CSI is also in need for this SBS and its associated SUE.Hence, we allocate orthogonal pilot sequence to each MUE and all the SCs share the same   pilot sequences for their   active users.Considering a TDD frame structure similar to that in [24], we propose the frame structure of this HetNet in Figure 2(a).
For comparison, the frame structure for achieving the global perfect CSI is presented in Figure 2(b).When global perfect CSI is preferred, at least  orthogonal pilot sequences are in need in one TDD frame.This would be intolerable for the considered macro cell massive MIMO systems in which the number of users would be large.Now, we pay attention to the special kind of CSI achieved by our proposed pilot reuse scheme.For any MUE, it will transmit its nonreused pilots to the BSs in the uplink training phase.Both the MBS and the SBSs can receive the pilot signal and then estimate the channels between this MUE and the BSs.Since all the pilot sequences of MUE are nonreused, both the MBS and the SBSs could have the perfect CSI while, for a SUE   , it will transmit its reused pilots to the BSs in the uplink training phase.Since the coverage regions of different SCs are nonoverlapped and the receive antenna gain of MBS is usually much larger than the SBS, it means that the minor channel coefficients between SUE   and other SCs can be ignored, and only the th SBS and the MBS can receive the pilot signal.Then the th SBS and the MBS can estimate the channels between SUE   and them.Note that SBS  cannot receive this pilot sequence from the SUE of other SCs, while the MBS receive this pilot sequence from all the SCs.It means that the th SBS will achieve the perfect CSI of SUE   , and the MBS (M) could just obtain the statistical knowledge of this SUE, for example, its path loss fading.
For simplicity, we rearrange the original channel matrix H as in Figure 3, where the first  BS columns correspond to the channel matrix of the MBS antennas and the following  SC columns correspond to the channel matrix of the antennas belonging to the 1st SBS and so on.It is similar that the first  −   rows correspond to channel matrix of the MUE and the next   rows correspond to the channel matrix of the SUE in the 1st SC and so on.
Besides, we use different colours to distinguish the different kinds of CSI.The perfectly known CSI is marked in green and the contaminated CSI is marked in red.The statistical knowledge of these contaminated CSI is assumed to be known at the BS, for example, the path loss fading of these users.The other space of the channel matrix is marked in white and stands for the minor channel coefficients which can be ignored.We call such kind of CSI as the "mixed CSI" for containing different kinds of CSI.

Outage-Constrained
Beamforming with Mixed CSI.The obtained mixed CSI, that is, some entries of H are perfect known at the BS while the others are statistically known or 0, brings a new technical challenge for the system downlink beamforming.Since some CSI is unknown, the QoS outage is inevitable.To guarantee the QoS and its outage, we impose a probabilistic QoS constraint as follows: where 0 <   < 1 indicates that the system should guarantee the QoS requirements for the th users with probability of at least 1 −   .Then, we can formulate the energy-efficient beamforming problem as a problem of minimizing the total transmit power while satisfying the system chance constraint as Recall that there are two different groups of users, that is, MUE and SUE.Hence, we can denote M ≜ {1, . . .,  −   } as the set of the MUE and let S ≜ { −   + 1, . . ., } stand for the set of all the SUE for convenience.Moreover, there are also two types of QoS constraints for such two types of users.As for the MUE, both the MBS and SBSs can obtain their perfect CSI.It means that the type of QoS for the MUE should be determinate; that is,   = 0, ∀ ∈ M.However, for the SUE, there is uncertainty in their CSI causing occasional QoS outage.Thus, the type of QoS for the SUE should be probabilistic; that is, it should involve outage probability constraint.Hence, the chance constrained beamforming with the achievable mixed CSI is formulated as P 1 : P 1 is to optimize the beamforming vectors {w , , ∀} for a total transmit power minimization.It aims to ensure that each MUE  is served with a minimal information rate of log 2 (1 +   ) and each SUE  is served with a minimal information rate of log 2 (1 +   ) with an outage probability no greater than   .
Usually, by using semidefinite relaxation (SDR), P 1 can be conservative transformed to a SDP.However, time complexity of standard interior-point SDP solvers could be as high as ( 6.5 ).Such intensive computation cost may prohibit SDR applicable to massive MIMO systems.Moreover, the conservatism of such SDR method is also hard to control.
In order to overcome this difficulty, we first decompose the complex problem P 1 into simpler subproblems and then propose an efficient iterative algorithm to combine these subproblems and solve P 1 iteratively.
As for the MUE, only the MBS transmits information symbols to them, and their beamforming vectors are generated by solving the following convex problem: From P 0 and P 1 , we find that the cross-tier interference of SUE and MUE is in relation with each other; that is, it needs the knowledge of the MUE beamforming vectors to calculate the interference power of SUE, and vice versa.Now we are ready to combine P 0 with P 1 and use them to solve P 1 iteratively.For convenience, the cross-tier interference of the th user caused by SBSs or MBS is denoted as    and    , respectively, where Considering the initial interference caused by the SBS and MBS as 0, then the optimization problem P 1 can be independently solved by each SBS and the beamforming vectors of w ,  , ∀ ̸ = 0, can be obtained.Since the CSI between MUE and all SBSs is perfectly known, MBS can broadcast h , , ∀ ̸ = 0,  ∈ M, to each associated SBS beforehand.Hence, the th SBS could compute the cross-tier interference it causes by calculating by and transfer it to the MBS.After getting the total cross-tier interference    = ∑  =1   , , the MBS can generate its beamforming vectors by solving P 0 .Since the CSI between MBS and SUE is statistically known, with the knowledge of {w 0, ,  ∈ M}, the cross-tier interference caused by the MBS to the th SUE can be calculated as where with   ,  = 1, . . .,  being the nonzero eigenvalues of (∑ ∈M w 0, w  0, ) in descending order (a detailed derivation of ( 13) and ( 14) can be found in the appendix).
(2) while Δ  ≥  ∑ Equation ( 13) means that Pr(  ≥    ) ≥ 1 −   ; that is, the cross-tier interference caused by the MBS to the th SUE will be larger than   with a probability no greater than   .
After the MBS broadcasting {   ,   ∈ S} to the th SBSs to update their local    , SBSs repeat the process of solving the problem P 1 mentioned above with the updated    .Such iterative process will stop until the predefined stop criterion is met.
Summarizing the above steps, we thus obtain the distributed algorithm as in Algorithm 1, where Δ  denotes the difference between the beamforming matrices in two successive iterations and  is utilized to constrain the maximum of this difference.The stop criterion means that the iterative process is permissible only when the following two conditions are satisfied simultaneously: (i) Δ  is not smaller than a ratio of the total transmit power; (ii) Δ  is shrinking.

Remark 1.
In Algorithm 1, it requires exchanging the knowledge of cross-tier interference between MBS and SBSs.Note that, in the process of exchanging interference knowledge, only one real value is transferred via the backhaul each time, which results in a total signaling of   + ( −   ) real values per iteration.Some other operations would also affect this total signaling.For example, if SUE broadcast their {w , , ∀ ̸ = 0,  ∈ S} to MBS instead of transferring   , , the total signaling would be   + 2    .Hence, mobile operators should choose an appropriate kind knowledge to exchange so as to yield a minimum backhaul signaling based on the architecture of the HetNet.
Remark 2. It is difficult to bound the exact number of iterations required by the proposed algorithm.However, simulations show that the proposed algorithm always converges within a few number of iterations.For example, in all of our Monte Carlo simulations (over 200,000 random channel realizations), the maximum number of iterations is 25.Most of the time, the proposed algorithm will converge within 5 times of iterations (see Section 4 for details).

Simulation Results
Monte Carlo simulations are performed to verify the conclusions and testify the performance of the proposed pilot reuse scheme and its associated downlink beamforming algorithm.The system performance is evaluated in terms of average total power consumption.When keeping the QoS constraints constant, the energy efficiency of the system is inversely proportional to its average total power consumption.

Simulation Settings.
The simulation parameters are set as follows.We consider a circular cell with one MBS and  = 4 small cells.The radii of the macro cell and small cell are  = 0.5 km and  = 40 m, respectively.Since the user assignment problem is independent of the beamforming problem considered in this paper, the users in the whole systems are scheduled in a round robin fashion for fairness consideration (please refer to [27][28][29] and the references therein for more about the user assignment problem).The goal of fairness user scheduling in this paper is to make the system operate at some point of its ergodic achievable rate region so as to provide useful insights into the impact of model parameters (e.g., number of antennas and outage probability constraint).We set the following: there are  = 10 active users in the whole system, where 6 users are uniformly distributed in the macro cell and each SC has one user uniformly distributed within its small cell coverage.We evaluate the average power consumption performance over user locations and channel realizations.The position of the MBS is set as (0, 0) and the overlaid SBSs are located at (0.5, 0.5), (−0.5, 0.5), (0.5, −0.5), and (−0.5, −0.5), respectively.Though there is no definite standard of the future 5G network, it is commonly thought that the future 5G network is an evolution of the current communication network (3G/4G network) [30].Hence, we take some system parameters of the current communication network for simulation, and the channel model is assumed as follows: (i) there are 600 subcarriers in the 10 MHz total system bandwidth, where the bandwidth of each subcarrier is 15 KHz; (ii) the channel between the MBS and UEs is using the 3GPP LTE pathloss mode: PL = 148.1 + 37.6 log 10 () (in dB) with  (in km) being the MBS-UE distance; (iii) within each SC, the channel is modeled as PL = 127 + 30 log 10 () (in dB) with  (in km) being the SBS-UE distance; (iv) the energy consumption parameters are set according to [31], that is,  0 = 0.388,   = 0.052, ∀ ̸ = 0,  0 = 189 mW, and   = 5.6 mW, ∀ ̸ = 0; and (v) the noise covariance of receivers is set as −127 dBm.

Comparison with the Perfect CSI Case.
It would be interesting to simulate the performance gap between the proposed low complexity algorithm and the optimal beamforming method with chance constraint.However, to the best of the authors' knowledge, there has not been such optimal methods in existing work.Thus, we propose to compare Algorithm 1  with the following baseline which is labeled as "full CSI" in Figures 4 and 5.In this baseline, we assume that the fully coordinated MBS and SBSs are able to share both the data and channel information.Then CoMP-JP can be used to convert the intercell interference into useful signal [32].In other words, this information is available at all BSs, who share their antennas together to serve the users jointly.If the beamforming method in CoMP-JP is optimal, this baseline provides a strictly upper bound for the maximum sum-rate.Since there exists QoS outage constraint in Algorithm 1 which may be helpful to save some transmit power, this baseline is likely to be the upper bound of Algorithm 1 especially when   is tiny.Now, we introduce the beamforming method for CoMP-JP.In this condition, the problem of minimizing the total power consumption while satisfying the QoS constraints is formulated as follows: where Γ  is shown as with After applying SDR by dropping the rank constraints, we can relax (15) to a semidefinite optimization problem on standard form (see details in [33,Chapter 18]) min where W  is positive semidefinite with W  = w  w   .It can be proved that these rank constraints can be relaxed without losing optimality [33,Chapter 18].It means that the original problem (15) can be solved as a convex optimization problem (17).

Performance of the Proposed Algorithm.
In this part, we will present the performance of the proposed algorithm and compare it with the optimal COMP-JP method mentioned above.
First of all, we focus on the impact of the different numbers of antennas at the MBS and SBS. Figure 4 compares the power consumption of the full CSI scenario (or perfect CSI scenario) and the pilot reuse scenario with different numbers of antennas.The solid lines show the case when the full downlink CSI is available at the BSs, while the dotted lines depict the one when mixed CSI is available at BSs.The given QoS requirement of each user is set as {  = 3, ∀} and the QoS outage probability requirement for the mixed CSI scenario is set to be   = 0.1, ∀ ∈ S.
In the mixed CSI scenario, when  SC is fixed, it will be more energy efficient when  is not too large or small.This is because that when  is not too large, the transmit power is much larger than the circuit power consumption, and the consumed power is mainly provided for the data transmission.However, with the increasing of , when the circuit power consumption is comparable to or even dominates the transmit power, too much additional circuit power is consumed by the large number of active antennas and therefore decrease the energy efficiency.In other words, appropriate number of active antennas could improve the energy efficiency of the system.It has the same phenomenon in the full CSI scenario.
Then, we focus on the impact of the number of antennas at each SBS,  SC .When  is not so large, the energy efficiency decreases with the increasing of  SC , while a larger  SC will yield a higher energy efficiency when  is large enough.It tells us that, when  is not large enough, it is not favorable to extract antennas from the macro cell and settle them in the SCs which will leave not enough macro cell antennas to handle the MUE.When  is large enough which is sufficient for the MUE, appropriate extraction antennas from the MBS to the SCs will not only improve the throughput performance of SBSs but also decrease circuit power consumption of the MBS.
Besides, we find that the power consumed by the proposed algorithm is larger than that the full CSI scenario when the number of  is small.With the increasing of , the proposed algorithm is likely to consume less power than the full CSI scenario.This interesting phenomenon is mainly due to the fact the full CSI scenario and the proposed scheme have different QoS requirements.In the full CSI scenario, the SINR constraint in (15) must be satisfied for all channel realizations.On the other hand, in the proposed scheme, SINR outage is allowed; that is, the SINR constraint in ( 8) is satisfied with probability at least 1−  .Hence, the QoS requirement for the full CSI scenario is more stringent, which may lead to larger transmit power than that of the proposed scheme.
Moreover, we find that, when choosing the most energy efficient  (about 35), the average energy consumption of the proposed pilot reuse scheme is larger than that of the full CSI scenario, and the energy consumption gap of them is no larger than 0.3 dB.Though consuming some extra energy than the full CSI scenario, our scheme will greatly compress the number of pilot sequences especially for those highly densified network with large number of small cells.
Hence, we summarize that choosing a suitable number of the total antennas and allocating an appropriate portion of them to the SCs would intensively enhance the energy efficiency of the whole system and the proposed algorithm is especially fit for the massive MIMO systems.
Figure 5 compares the average power consumption for different QoS constraints in the full CSI scenario and the mixed CSI scenario when  = 40 and   = 0.1, ∀ ∈ S. We observe that, when the required QoS constraints are not large, the systems with larger  SC consume less power even than the full CSI cases.This is for the reason that   will relax the QoS requirements to some extent and save some energy consequently while, with the increasing of the required QoS constraints, the gap of the power consumption between the full CSI scenario and the mixed CSI scenario tends to be static, which ensures that the proposed algorithm is widely adaptive to different QoS requirements.
Figure 6 shows the impact of the QoS outage probability constraints on the power consumption when  = 40 and   = 3, ∀.We find that when decreasing   from 0.2 to 0.01 it consumes more power.However, when  SC is fixed, the difference of power consumption is not evident.Hence, we conclude that, over a reasonable range of QoS outage probability, the consumed power of our scheme is likely to be static.Outage probability Total power consumption (dBm) Figure 7 describes the feasible rate of the actual SINR satisfaction probabilities and the presented results are achieved from 40, 000 random channel realizations.When   = 0.05 and the given SINR constraint for each user is   = 3, we obtain the actual SINR satisfaction probabilities of one SUE in Figure 7.We find that Pr{Γ  >   } is 0.9498 ≈ 1 −   which meets our QoS outage probability constraints strictly.
Figure 8 shows the total number of iterations of the proposed algorithm with different .In the presented results, each point is obtained from 10, 000 random channel realizations.We find that the total number of iterations tends to be static with large number of , for example, 3.5 times for  SC = 1 and 3 times for  SC = 2 or 3.In all of our simulations, the maximum number of iterations is 25 when {,  SC } = {16, 1}.With the increasing of , it will take fewer iterations to meet the stop criterion.When  SC = 1, it always costs more iterations to meet the stop criterion.Hence, it is favorable to allocate more antennas at SCs in view of spending fewer times of iterations.
In summary, the above simulations show that the carefully chosen {,  SC } is acquired to maximize the energy efficiency of such two-tier HetNet.When {,  SC } is carefully chosen, the average energy consumption of the proposed pilot reuse scheme is larger than that of the full CSI scenario, and the gap between them is about 0.3 dB.Over a tolerable range of the QoS outage probability, the power consumption of our scheme will stay nearly static.Specifically, when  and  SC are sufficient large or   is not too large, the proposed algorithm consumes less power than the perfect CSI scenario.Hence, the pilot reuse scheme and its associated beamforming algorithm are quite suitable for massive MIMO systems.

Conclusion
In this paper, we focus on a kind of two-tier networks, where massive MIMO is deployed in combination with small cell networks.However, it is not permissible to assign orthogonal pilot sequences for all the supported users due to the large number.We propose a pilot reuse scheme based on the unique architecture of this networks and find that a special kind of mixed CSI is yielded when using our pilot reuse scheme.Since mixed CSI may cause QoS outage the UEs, we formulate the downlink beamforming problem as a chance constraint programming.Then we decompose the original problem into simpler subproblems and provide an algorithm to combine these subproblems and solve them iteratively in a distributed fashion.The numerical results show that the average energy consumption of the proposed pilot reuse scheme and its associated beamforming algorithm is close to that of the perfect CSI scenario and the energy consumption is nearly static over a wide range of the QoS outage probabilities.The QoS requirements and QoS outage probability are also strictly controlled in our algorithm.Our scheme will not only greatly compress the number of pilot sequences, but also owns a high performance with low complexity, which is especially suitable for massive MIMO systems and will be crucially useful to put such multitier massive MIMO systems into practice.

Figure 1 :
Figure 1: One typical architecture of a two-tier network.A macro cell tier is overlaid with a secondary tier of small cells.The interfering signals are marked by dotted lines and the intended signals are marked by solid lines.

Figure 3 :
Figure 3: The rearranged H for the two-tier network.

1 N SC = 2 N SC = 3 Figure 4 :
Figure 4: The average power consumption of the whole system with different numbers of  and  SC .

1 N SC = 2 N SC = 3 Figure 5 :
Figure 5: The average power consumption of the whole cell for different QoS constraints.The QoS constraints are measured in bps/Hz per user.

Figure 6 :Figure 7 :
Figure 6: The average power consumption of the whole cell for different outage probability constraints.

3 Figure 8 :
Figure 8: The number of iterations of the proposed algorithm with different numbers of  and  SC .