A low-complexity algorithm for the joint antenna selection and ......lem is to applying antenna selection [11–13] to decide optimal subset of BS transmit antennas for decreasing
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
RESEARCH Open Access
A low-complexity algorithm for the jointantenna selection and user scheduling inmulti-cell multi-user downlink massiveMIMO systemsSaidiwaerdi Maimaiti*, Gang Chuai, Weidong Gao, Kaisa Zhang, Xuewen Liu and Zhiwei Si
Abstract
The massive MIMO (multiple-input multiple-output) technology plays a key role in the next-generation (5G) wirelesscommunication systems, which are equipped with a large number of antennas at the base station (BS) of anetwork to improve cell capacity for network communication systems. However, activating a large number of BSantennas needs a large number of radio-frequency (RF) chains that introduce the high cost of the hardware andhigh power consumption. Our objective is to achieve the optimal combination subset of BS antennas and users toapproach the maximum cell capacity, simultaneously. However, the optimal solution to this problem can beachieved by using an exhaustive search (ES) algorithm by considering all possible combinations of BS antennas andusers, which leads to the exponential growth of the combinatorial complexity with the increasing of the number ofBS antennas and active users. Thus, the ES algorithm cannot be used in massive MIMO systems because of its highcomputational complexity. Hence, considering the trade-off between network performance and computationalcomplexity, we proposed a low-complexity joint antenna selection and user scheduling (JASUS) method based onAdaptive Markov Chain Monte Carlo (AMCMC) algorithm for multi-cell multi-user massive MIMO downlink systems.AMCMC algorithm is helpful for selecting combination subset of antennas and users to approach the maximum cellcapacity with consideration of the multi-cell interference. Performance analysis and simulation results show thatAMCMC algorithm performs extremely closely to ES-based JASUS algorithm. Compared with other algorithms in ourexperiments, the higher cell capacity and near-optimal system performance can be obtained by using the AMCMCalgorithm. At the same time, the computational complexity is reduced significantly by combining with AMCMC.
Keywords: 5G, Massive MIMO systems, Antenna selection, User scheduling, Adaptive Markov chain MonteCarlo algorithm
1 IntroductionIn order to satisfy the rapidly increasing requirements forhigh data rate in current wireless communication systems,a new massive MIMO (multiple-input multiple-output)technology was introduced in [1–3]. Massive MIMO tech-nique plays a key role to enhance the cell capacity withoutincreasing system bandwidth or base station (BS) trans-mission power for the 5G network systems [4]. The keyidea of the massive MIMO technique is to install a large
amount of transmit antennas at the BS of a cellular andprovide services for several users sharing the samespectrum resources. However, as the number of BS anten-nas and users increases, the combination complexity andhardware cost also increase dramatically. Therefore, whenthe numbers of BS transmit antenna and active users areextremely large, the joint antenna selection and userscheduling (JASUS) algorithm [5–8] can be adopted as anapproach to decide the radio frequency (RF) chain config-uration to improve the cell capacity in multi-cell massivemulti-user MIMO systems.In a practical network, one of the key challenges in
* Correspondence: [email protected] Laboratory of Universal Wireless Communications, Ministry of Education,Beijing University of Posts and Telecommunications, Beijing 100876, People’sRepublic of China
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 https://doi.org/10.1186/s13638-019-1529-7
hardware cost and power consumption because theelement of each antenna needs a complete RF chain thatconsists of RF amplifiers and analog-to-digital con-verters, which are very pricey and are the main elementsof the power consumption at the BS [9]. Differentschemes were used in many types of research, such ashybrid precoding and spatial modulation, to reduce thecost of the hardware and the power consumption of thesystem [10]. One of the best schemes to solve this prob-lem is to applying antenna selection [11–13] to decideoptimal subset of BS transmit antennas for decreasingthe required number of high pricey RF chains while de-creasing the resulting network performance loss.However, in multi-cell multi-user massive MIMO sys-
tems, only a limited number of transmit antennas are se-lected to provide services for active users scheduled.Hence, if the number of users exceeds, the number ofselected transmit antennas, user scheduling must be per-formed because different wireless channels have differentproperties. High cell capacity can be obtained by sched-uling users with the high channel quality. Therefore, theresearch of JASUS method for multi-cell multi-usermassive MIMO systems is necessary.Recently, only a few types of research have studied a
low complexity JASUS for downlink massive multi-userMIMO systems. Benmimoune et al. [14] proposed atwo-step JASUS scheme for downlink multi-user massiveMIMO systems. It successively closed unnecessary an-tennas and removes undesired users which contributelittle to system performance. However, due to the highcomputational complexity, this algorithm can only beemployed to the scenarios with a smaller number of can-didate antennas and user sets. Thus, using this algorithmin practical multi-cell multi-user massive MIMO systemscenarios is difficult. Olyaee et al. [15] proposed a JASUSmethod based on zero-forcing (ZF) precoding algorithmfor single-cell multi-user massive MIMO downlink sys-tems. Though the ZF precoding method has a high sys-tem performance, it also has a very high computationalcomplexity. For distributed downlink multi-user massiveMIMO system, a JASUS method was proposed in [16]by Xu. et al. It successively obtains the majority of gainwith limited backhaul capacity. Lee et al. [17] proposed arandom antenna selection algorithm, the algorithm canprovide significant capacity efficiency gain, but it is diffi-cult to use for multi-cell multi-user massive MIMO sys-tems. However, the above researches focused on single-cell multi-user massive MIMO systems. Thus, the re-search of JASUS method for multi-cell multi-usermassive MIMO systems remains a largely open area.Therefore, the novel JASUS algorithm with consideredmulti-cell interference, which causes no or only a fewdecreases of system performance, represents a newpromising research topic.
In this paper, we consider the problem of JASUS inmulti-cell multi-user massive MIMO downlink system op-erating with TDD mode. Considering the trade-off be-tween cell capacity and complexity, we proposed a low-complexity algorithm for JASUS method based onAMCMC. In our proposed method, only a small subset ofBS transmit antennas is selected to serve predeterminedactive users, thus reduces the number of RF chains, avoidsuneconomical hardware costs, and reduces power con-sumption caused by the selection of unnecessary transmitantennas to provide the required services. The main con-tributions of our work are as follows.
1. A low complexity JASUS method based onAMCMC algorithm is proposed for downlinkmulti-cell multi-user massive MIMO systems.AMCMC algorithm is helpful for selectingcombination subset of antennas and users toapproach the maximum cell capacity whiledecreasing the resulting network performanceloss.
2. In this paper, we proposed updating rules for theselection probability of each base station transmitantenna and the scheduling probability for eachuser. In addition, we also proposed a newprojection strategy to satisfy the constraints ofantenna and user selection.
3. Performance analysis and simulation results showthat our proposed algorithm produced promisingresults. Compared with ES-based JASUS algorithm,the proposed algorithm achieved comparableperformance with low complexity. In addition, theAMCMC-based JASUS algorithm outperformsgreedy-based JASUS and norm-based JASUSmethods in terms of cell capacity and SER (symbolerror rate) performance.
Notation: Symbol ℂ denotes the set of complex num-bers, vectors are denoted by using lower-case bold letters,matrices are denoted by using bold letters, |.| denotes theabsolute value of a scalar, ‖·‖F denotes the Frobenius normfunction, and (.) represents the binomial coefficient.The remaining content is organized as follows. In
Section 2, the system model and capacity maximizeproblem formulation are described. In Section 3, we for-mulate the problem of JASUS method based on AMCMCin multi-cell multi-user massive MIMO systems. Section 4presents the simulation setup and assumption. In Section5, we discuss the simulation results and analyze the com-plexity; finally, this work is concluded in Section 6.
2 System model and problem formulationIn this part, we simply give the system model for multi-cell massive multi-user
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 2 of 14
MIMO downlink systems with the system capacityformulation model with consideration of the multi-cellinterference.
2.1 System modelAs shown in Fig. 1, the considered scenario is a multi-cell multi-user massive MIMO downlink system operat-ing in TDD mode, and the cell capacity maximizingproblem is studied with consideration of the inter-cellinterference. The system is composed of Β hexagonalcells. All B BSs, where B = {1, 2,…, Β} are installed withM antennas and serves U (M ≥U ≥ 1) single-antennausers in each cell. The block-fading channel model isassumed. We assume that BS can select N transmit an-tennas among the M transmit antennas and scheduleK(K ≤N) users among the U users within the cell to beserved simultaneously. The channel vector giju ∈ ℂ
M
from the jth BS and user u in cell i can be expressed as
giju ¼ffiffiffiffiffiffiffiβiju
qhiju ð1Þ
where βiju denotes the large scale channel fading be-tween jth BS and user u in cell i, including shadowing
and path loss. hiju is the small-scale fading vector, andhiju = [hiju1, hiju2,…, hijuM]
T ∈ ℂM. Then, the overalldownlink transmission matrix Gij ∈ ℂ
M ×U between theBS in cell j and all users in cell i can be expressed as
Gij ¼ gij1; gij2;…; gijUh i
¼ HijD12ij∈ℂ
M�U ð2Þ
where Hij = [hij1, hij2,…, hijU] ∈ ℂM ×U is the overall
small-scale fading matrix and Dij = diag(βij1, βij2,…, βijU).Our objective is to find optimum combinations subset
of BS antennas and users to approach the maximum cellcapacity while decreasing the resulting network perform-ance loss. Furthermore, we will decrease the number ofexpensive RF chains and avoid the uneconomic costs ofthe hardware and decrease power consumption causedby selecting undesired antennas to provide the require-ment of service.
2.2 Problem formulationIn the downlink system, the signal received by users incell i can be written as
Fig. 1. Illustration of the downlink multi-cell multi-user massive MIMO network transmission model
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 3 of 14
ri ¼ ffiffiffiffipt
p XΒj¼1
GHij W js j þ ni ð3Þ
where pt denotes the transmitted power, sj = [sj1, sj2,…,sjU]
T ∈ℂU is the transmit signals for users in cell j,Wj ∈ ℂM ×U is the precoding matrix of BS in cell j, andni = [ni1, ni2,…, niu]
T ∈ℂU is the noise vector at the uthuser in cell i. The downlink signal received by user u incell i can be written as
riu ¼ ffiffiffiffipt
pgHiiuwiusiu|{z}
desired signal
þ ffiffiffiffipt
p Xb≠u
gHiibwibsib|{z}intra‐cell interference
þ ffiffiffiffipt
p Xj≠i
XUb¼1
gHijuwjbsjb|{z}inter‐cell interference
þniu
ð4Þ
where wjb is the bth column of precoding matrix Wj.Formula (4) including the desired signal for the user u incell i, intra-cell interference signal (comes from otherusers in the same cell), and inter-cell interference signal(comes from other cells), respectively.Since we used the zero-forcing (ZF) methods and sup-
posing that the channel state information (CSI) is per-fectly known at the BS, the intra-cell interference signal(second term) in function (4) drops to zero according toprevenient works [14, 18–20]. Therefore, function (4)becomes
riu ¼ ffiffiffiffipt
pgHiiuwiusiu|{z}
desired signal
þ ffiffiffiffipt
p Xj≠i
XUb¼1
gHijuwjbsjb|{z}inter‐cell interference
þniu ð5Þ
We assume that BS can choose N transmit antennasamong the M transmit antennas, and schedule K(K ≤N)users among the U users within the cell to be servedsimultaneously. For convenience, we give the selectedsubset of the antenna and scheduled subset of user indi-cator functions are ω and ω,
ω ¼ ωmf gMm¼1 ð6Þ
ω ¼ ϖuf gUu¼1 ð7Þwhere ω and ω are binary vectors that include twovalues 0 and 1 to indicate if a given antenna or agiven user is selected. (e.g., 1 → selected, 0 →unselected).For making an easy description, we will define to two
indicator functions, which are Im(ω) ≜ ωm ∈ {0, 1} andIuðωÞ≜ϖu∈f0; 1g , respectively. We use these to indicatewhether the mth BS antenna and the uth user are selectedor not, respectively. s
ifωg∈ℂK denotes transmit signal
vector, sub-block channel matrix of corresponding
denotes by Gijfω;ωg∈ℂN�K and n
ifωg∈ℂK is noise vector,
respectively. Finally, in order to denote the joint antennaand user selection, we employ the 2-tuple Ω≜ðω;ωÞ . Inorder to have an easy explanation, we will interchangeablyuse Ω and ðω;ωÞ for the following part. After using an-tenna selection and user scheduling method, function (5)becomes
ri ϖuf g ¼ ffiffiffiffipt
pgHii ω;ωf gwi ϖuf gsi ϖuf g
þ ffiffiffiffipt
p Xj≠i
XUb¼1
gHij ω;ωf gw j ϖbf gs j ϖbf g
þ ni ϖuf g ð8Þwhere w jfϖbg is bth column of precoding matrixW
jfω;ωg∈ℂN�K and nifϖug is the noise at the uth user
in cell i.
2.3 Capacity of massive MIMOAccording to the aforementioned discussion, the re-ceived signal-to-interference-plus-noise ratio (SINR) forthe user u ∈U (which is connected to cell i) with a se-lected channel vector g
ijfω;ωg can be written as
SINRi ϖuf g ¼pt g
Hii ω;ωf gwi ϖuf g
��� ���2Xj≠i
pt gHij ω;ωf gw j ϖuf g
��� ���2 þ ni ϖuf g�� ��2 ð9Þ
Considering the inter-cell interference, the formula ofsum capacity for cell i can be expressed as
Cisum G ω;ωf g� � ¼ log2 det I þ
pt GHii ω;ωf gWi ω;ωf g
��� ���2Xj≠i
pt GHij ω;ωf gW j ω;ωf g
��� ���2 þ ni ωf gnHi ωf g
0BBB@
1CCCA
0BBB@
1CCCA
ð10ÞOur target is to jointly select the optimal combination
sets of BS transmit antenna and active user to approachthe maximum cell capacity while decreasing the compu-tational complexity. Hence, the problem of JASUS canbe written as
ΦC ¼ maximizeω;ω
Cisum G ω;ωf g� � ð11Þ
subject to
XUu¼1
Iu ωð Þ ¼ K ð12Þ
XMm¼1
Im ωð Þ ¼ N ð13Þ
Addressing the aforementioned problem by employingan ES method needs to evaluating the cell capacity of
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 4 of 14
φ≜CNM � CK
U joint antenna and user combinations, whereCN
M and CKU are the binomial coefficient. This fact indi-
cates that the ES method cannot be used in current themassive MIMO systems where U and M are very numer-ous, it leads to high computational complexity. Thus, alow complexity algorithm for JASUS is needed in order toobtain the best network performance with low computa-tional complexity.Obviously, formula (11) serves as the target function
in this study. Therefore, address the problem (11) inmulti-cell multi-user massive MIMO systems, we needto solve the following three main problems:
1. Inter-cell interference: the first problem is theinter-cell interference coming from other cells. Inorder to solve this problem, cells can adjust theirprecoding matrices, thus can eliminate or decreasethe interference from all users. In order to makebetter use of the coordinated massive MIMOtechnology, a group of users in the cells shouldbe scheduled so that each group of users in thecell has the biggest spatial separation with theinterference channels of users in neighboringcells.
2. Computational complexity: the second problem ishow to obtain the best combination subset of antennaand user in each cell with lower computationalcomplexity so as to decrease or eliminate theinter-cell interference and maximize the sumcapacity of all cells. We know that in multi-cellmassive MIMO systems, the various path lossbetween the antennas coming from neighboringcells and users coming from target cell also bringto much computational complexity same to theantenna selection and user scheduling. Thus, thecomputational complexity of the ES algorithmbecomes very large than of a single-cell JASUS.
3. CSI feedback cost: for massive MIMO systems, theperfect CSI feedback depends mainly on thenumber of active antennas and the users theysupport. Hence, in order to centralize processingfor selecting antennas and scheduling users acrosscells, the BS needs to exchange the overall CSI ofoverall combined subset antennas and users at eachscheduling period, it brings more burden for BS. Inaddition, when the number of BS antennas andusers in each cell increases, the cost of CSIincreases accordingly.
On the base of the aforementioned discussion, alow complexity scheme is needed from the practicalpoint of view for JASUS in multi-cell multi-usermassive MIMO scenarios to reduce the complexity of
function (11) while decreasing the cost of the CSIfeedback.
3 Joint antenna selection and user schedulingalgorithmIn this part, we presented two suboptimal iterative algo-rithms for JASUS before a discussion of the proposedAMCMC method.
3.1 Norm-based JASUS algorithmFirstly, we presented a norm-based JASUS method for ad-dressing the objective function (10). The norm-based JASUSscheme maximizes kGfω;ωgkF
, where ‖.‖F is the Frobenius
norm function. Let CNBsumðGfω;ωgÞ ¼ kGfω;ωgkF
. This
scheme including initialization step and iterative up-dating step, respectively. This both steps use the vec-tor norm as criteria which considerably decrease eachiterative computation complexity. The norm-basedJASUS problem is modeled as
ωNB ¼ arg maxω;ω
CNBsum G ω;ωf g� � ð14Þ
where ωNB is the combination selection indicator.Nevertheless, it has still a problem that when performingtransmit antenna selection and user scheduling simplybased on the F-norm criteria, it would sacrifice some cellcapacity. In sum, the norm-based JASUS method has ex-tremely low complexity, but it cannot guarantee a highsum capacity performance.
3.2 Greedy-based JASUS algorithmIn order to enhance the cell capacity over norm-basedJASUS, we presented a greedy-based JASUS algorithm.Unlike to the norm-based JASUS algorithm, the greedy-based JASUS algorithm maximizes the cell capacity ineach step. This method also includes initialization stepand iterative updating step. The greedy-based JASUSproblem is modeled as
ωGR ¼ arg maxω;ω
CGRsum G ω;ωf g� � ð15Þ
where ωGR is the combination selection indicator. Com-pared with the norm-based JASUS algorithm, thismethod has a good capacity performance, but it has highcomputational complexity.
3.3 AMCMC-based JASUS algorithmAlthough the greedy-based JASUS method improved thecell capacity, it ignores the computational complexity ofall system. For the actual network communication sys-tem, it will not have commercial value or attraction.Hence, considering the trade-off between cell capacity
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 5 of 14
and complexity, we proposed a low-complexity JASUSmethod based on AMCMC algorithm and its descriptionis as follows.MCMC [21] is a method of generating random sam-
ples, which is often used to calculate statistical estima-tion, marginal probability, and conditional probability.MCMC algorithms depend on (Markov) sequences withlimit distributions corresponding to interest distribu-tions. In the past decades, it has been widely used inmany fields such as engineering and statistics [22]. Thekey idea of the MCMC method is that Markov chainsare simulated in state space X, and the stable distribu-tion of the chains is the target distribution π [23].In order to address our objection problem in (11), we
must tackle tow essential issues when applying AMCMCalgorithm. The first problem is how to provide a pro-posal distribution of candidate samples LMCMC. The sec-ond problem is how to design the most suitableupdating rule for the proposal distribution.
3.3.1 Derivation of the candidate sampling distribution forthe MCMC methodThe biggest advantage of the MCMC method is that itcan search the “elite samples” instead of exhaustivelysearching the whole samples. At iteration t, the samples
fΩℓ;t ¼ fωℓ;t;ωℓ;tggLMCMCℓ¼1 from an MCMC method can
be employed to estimate the maximum value of the tar-get function Ci
sumðGfω;ωgÞ
Φ�C ¼ arg max
Ωℓ;t≜ ωℓ;t ;ωℓ;tð Þ; ℓ¼1;2;…;LMCMC
Cisum G ωℓ;t ;ωℓ;tf g� �
ð16Þwhere LMCMC denotes the total number of samples, andΦ�
C is the estimated value of formula (11).Given that scheduled of per user ϖu and selection of
each antenna ωm are binary variables, we use the Boltz-mann distribution of the objective function Ci
sumðGfωℓ;t ;ωℓ;tgÞ with a suitable temperature τ
π Ωℓ;t� �
Ωℓ;t≜ ωℓ;t ;ωℓ;tð Þ¼ exp Ci
sum Ωℓ;t� �
=τ� ��
Γð17Þ
where Γ ¼PΩℓ;t≜ðωℓ;t ;ωℓ;tÞ expfðCisumðΩℓ;tÞ=τÞg is a
normalization constant in the MCMC method that canbe neglected. Thus, maximizing Ci
sumðΩℓ;tÞ is equivalentto maximizing π(Ωℓ, t), and π(Ωℓ, t) is the targetdistribution.In order to prove the MCMC method for searching
the distribution π(Ωℓ, t), we use a MIS (metropolized in-dependence sampler) [23], which is a generic MCMCmethod. The step is as follow. An initial value Ω[0], t ischosen randomly. Given the current sample Ω[ℓ], t,
� A candidate sample Ω[new], t is drawn from proposaldistribution Ψ(Ωℓ, t; Rt − 1).
� Simulate u uniform[0, 1], and according to theaccepting probability α(Ω[ℓ], t,Ω[new], t), let
ð19ÞAfter LMCMC iterations, we can achieve a set of sam-
ples fΩ½0�;t;Ω½1�;t;…;Ω½LMCMC�;tg , which is subjected todistribution π(Ωℓ, t).
3.3.2 Derivation of updating rule for the AMCMC algorithmIn this part, we provide updating rule for the pro-posal distribution. For the AMCMC method, the jointproposal distribution is proportional to the product ofBernoulli distributions, namely
Ψ Ωℓ;t ;Rt−1� �
≜ΠU
u¼1 pIu ωℓ; tð Þu;t−1 1−pu;t−1
� �1−Iu ωℓ; tð Þ
Γ0|{z}
≜Ψ ωℓ; t ;Pt−1ð Þ
�ΠM
m¼1 gIm ωℓ; tð Þm;t−1 1−gm;t−1
� �1−Im ωℓ; tð Þ
Γ0|{z}
≜Ψ ωℓ; t ;gt−1ð Þ
ð20Þ
where pu denotes the probability of the uth user be-ing selected for communicating with the BS. That is,ϖu ∼ Ber(pu) for u = 1, 2, …, U, and gm is the prob-ability of the mth BS antenna being selected. That is,ωm ∼ Ber(gm) for m = 1, 2, …, M. We use the indicatorfunctions Im(ωℓ, t) and Iuðωℓ;tÞ to indicate whether themth BS antenna and the uth user are selected or not, re-
¼ fgm;tgMm¼1, and Γ′ is a normalization constant that
can be ignored in the AMCMC. The adaptationscheme is employed to adjust the parameterized pro-posal distribution Ψ(Ωℓ, t; Rt − 1) and minimize theKullback-Leibler divergence [24, 25] between the tar-get distribution π(Ωℓ, t) and the proposal distributionΨ(Ωℓ, t; Rt − 1), namely
D π Ωℓ;t� �
Ψ Ωℓ;t;Rt−1� ���� ¼ XLMCMC
ℓ¼1
π Ωℓ;t� �� log
π Ωℓ;t� �
Ψ Ωℓ;t;Rt−1� � !
ð21Þ
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 6 of 14
It is observed that ~D ¼ πðΩℓ;tÞ � logπðΩℓ;tÞ−D½πðΩℓ;tÞkΨðΩℓ;t ;Rt−1Þ� is a convex function [26]. Hence, theminimization of the Kullback-Leibler divergenceD[π(Ωℓ,t)‖Ψ(Ωℓ,t; Rt − 1)] w.r.t. R can be achieved when∂~D=∂R ¼ 0. Thus, ~D can be written as
~D ¼XLMCMC
ℓ¼1
π Ωℓ;t� �� logΨ Ωℓ;t;Rt−1
� � ð22Þ
We set the partial derivative of (22) to zero with re-spect to R. Then, the Eq. (22) can be written as
∂∂R
XLMCMC
ℓ¼1
π Ωℓ;t� �� logΨ Ωℓ;t ;Rt−1
� � ¼ 0 ð23Þ
where the partial derivatives of logΨ(Ωℓ,t;Rt − 1) with re-spect to pu and gm are respectively given by
∂∂pu
logΨ Ωℓ;t;Rt−1� � ¼ Iu ωℓ;t
� �−pu;t−1
1−pu;t−1� �
pu;t−1ð24Þ
∂∂gm
logΨ Ωℓ;t;Rt−1� � ¼ Im ωℓ;t
� �−gm;t−1
1−gm;t−1
� �gm;t−1
ð25Þ
By substituting (24) and (25) into (23), we obtain
XLMCMC
ℓ¼1
π Ωℓ;t� �� Iu ωℓ;t
� �−pu;t−1
1−pu;t−1� �
pu;t−1¼ 0 for u ¼ 1; 2;…;U
ð26ÞXLMCMC
ℓ¼1
π Ωℓ;t� �� Im ωℓ;t
� �−gm;t−1
1−gm;t−1
� �gm;t−1
¼ 0 for m ¼ 1; 2;…;M
ð27Þ
Given a number of samples fΩℓ;t ¼ fωℓ;t ;ωℓ;tggLMCMCℓ¼1
drawn from target distribution π(Ωℓ, t), the Monte Carloestimate of ∂~D=∂pu and ∂~D=∂gm are
1LMCMC
1
1−pu;t−1� �
pu;t−1
XLMCMC
ℓ¼1
Iu ωℓ;t� �
−pu;t−1
24
35 ð28Þ
1LMCMC
1
1−gm;t−1
� �gm;t−1
XLMCMC
ℓ¼1
Im ωℓ;t� �
−gm;t−1
24
35
ð29ÞApplying the Robbins-Monro stochastic approxima-
tion scheme [26], we can achieve the recursive up-date function to close to the root of ∂~D=∂pu ¼ 0 and∂~D=∂gm ¼ 0, namely
pu;t ¼ pu;t−1 þrt
1−pu;t−1� �
pu;t−1
� 1LMCMC
XLMCMC
ℓ¼1
Iu ωℓ;t� �
−pu;t−1
!
ð30Þ
gm;t ¼ gm;t−1 þrt
1−gm;t−1
� �gm;t−1
� 1LMCMC
XLMCMC
ℓ¼1
Im ωℓ;t� �
−gm;t−1
!ð31Þ
where rt denotes the sequence of decreasing step sizes[27]. In addition, we can simplify formulas (30) and (31),because (1 − pu, t − 1)pu, t − 1 and (1 − gm, t − 1)gm, t − 1 hasno significant impact on the convergence of (30) and(31). Hence, Eqs. (30) and (31) becomes
pu;t ¼ pu;t−1 þ rt � 1LMCMC
XLMCMC
ℓ¼1
Iu ωℓ;t� �
−pu;t−1
!
ð32Þ
gm;t ¼ gm;t−1 þ rt � 1LMCMC
XLMCMC
ℓ¼1
Im ωℓ;t� �
−gm;t−1
!
ð33ÞThe updated proposal distribution Eqs. (32) and (33)
are iteratively used with the objective to close to the tar-get distribution.
3.4 Constraints for the AMCMC-based JASUS problemFor the JASUS problem, the vectors ω and ω are subjectto constraint functions (12) and (13), respectively. How-ever, employing functions (32) and (33) with MIS togenerate samples, but we cannot ensure that the samplesmeet the constraints (12) and (13). In order to ensurethat the samples drawn from (32) and (33) meet theconstraints (12) and (13), we propose a new projectionstrategy. For convenience, we only introduce the pro-posed projection strategy for ω because the similar pro-jection strategy can be used to ω.Assume the sample ωℓ, t drawn from Ψ(·; gt − 1) at
the tth iteration. We define the two sets, which areϕ0 = {m : Im(ωℓ, t) = 0} and ϕ1 = {m : Im(ωℓ, t) = 1}, re-spectively. We use these to collect the indices for the unse-lected and selected BS antennas, respectively. The following
projection strategy is applied ifPM
m¼1Imðωℓ;tÞ≠N :
� IfPM
m¼1Imðωℓ;tÞ < N , then the proposed projectionstrategy sequentially selects the BS antenna with thebiggest probability from the set ϕ0 to the set ϕ1 until
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 7 of 14
|ϕ1| =N, where ϕ1 is the number of elements of theset ϕ1.
� IfPM
m¼1Imðωℓ;tÞ > N , then the BS antennas withthe smallest probability in the set ϕ1 are closedsequentially according to the proposed projectionstrategy until |ϕ1| =N.
3.5 Convergence analysis for the AMCMC-based JASUSalgorithmIn order to obtain a higher convergence rate, theprobability parameters Rt − 1 of the proposal distribu-tion Ψ(Ωℓ, t; Rt − 1) are adjusted. In this paper, we useliterature [28] to explain the convergence problem,because our proposed method gives a similar descrip-tion of the convergence problem and proves the ef-fectiveness of the method through analysis. Besides,the complexity of the MCMC algorithm has been
proved to be only related to sample size LMCMC in [22].The proposed adaptive strategy requires less sample sizeand iteration times, which can significantly improve theconvergence speed of the MCMC algorithm.
3.6 Constrained AMCMC-based JASUS algorithmOn the base of the aforementioned discussion, we canbe written the proposed AMCMC-based JASUS algo-rithm by the following steps. At iteration t, LMCMC
samples fΩℓ;t ¼ fωℓ;t;ωℓ;tggLMCMCℓ¼1 from the MCMC
method are can be generated by employing MISaccording to proposal distribution Ψ(Ωℓ, t; Rt − 1).Then, the new proposal distribution Ψ(Ωℓ, t; Rt) willbe updated by the Kullback-Leibler divergence untilit approach the target distribution π(Ωℓ, t). The de-tailed AMCMC-based JASUS algorithm is describedas follows.
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 8 of 14
4 Simulation configurationIn this section, the simulation configuration and simula-tion parameters are described. The considered scenariois a multi-cell multi-user massive MIMO downlink sys-tem operating in TDD mode with Β = 7, as shown inFig. 2. The simulation is done with a static networksimulator. The key simulation parameters are summa-rized in Table 1. In this simulation, we assume that theCSI is perfectly known at the transmitter, the totalpower is uniformly allocated among the transmit anten-nas. The system composed of Β hexagonal cells. All BBSs, where B = {1, 2,…, Β} are installed with M antennasand serve U single-antenna users in each cell. Each BS islocated at the cell center while U single-antenna usersare randomly located in the cell. There is no user move-ment and handover during the simulation process.
5 Simulation results and analysisIn this section, we provide numerical results and compu-tational complexity analysis of the proposed algorithmby simulative evaluation.
5.1 Performance evaluationAs can be seen the result from Fig. 3, the multi-cell multi-user massive MIMO downlink system using differentJASUS method at various SINR. We find that the sumcapacity of the AMCMC-based JASUS algorithm is veryclose to the maximum capacity result obtained by the ES-based JASUS algorithm with a wide range of SINRs. Forexample, when the SINR is 20 dB, the achieved values ofcell capacity by using ES and AMCMC algorithms are
50.7 and 49.9 b/s/Hz, respectively. Ninety-eight percent ofthe optimal capacity is obtained by our proposed method.The result shows that the AMCMC-based JASUS methodhas a good performance compared to the greedy-basedJASUS method with the wide range of SINRs, bothAMCMC-based JASUS and greedy-based JASUS methodshave better capacity performances compared to the norm-based JASUS algorithm. This simulation result shows that,when the SINR is 30 dB, ES, AMCMC, and greedy algo-rithms enhanced cell capacity of approximately 9.1, 8.2,and 5.8 b/s/Hz, respectively.Figure 4 shows the increase of sum cell capacity dur-
ing each iteration for AMCMC-based JASUS algorithmwith SINR = 20 dB. As can be seen from Fig. 4, we findthat the AMCMC-based JASUS converges after aboutt = 30 iterations. As expected, the sum cell capacity ob-tained by AMCMC-based JASUS strictly monotonicallyincreased with a number of iteration.Now, the assumed scenario is a multi-cell multi-user
massive MIMO downlink system with 50 active users(U = 50) at SINR = 20 dB, and we assume a differentnumber of transmit antenna, M, from 16 to 60. Sixteenantennas (N = 16) were selected to be used by the trans-mitter and ten users (K = 10) were served. As can beseen in the result from Fig. 5, the cell capacity differencebetween the ES-based JASUS and AMCMC-basedJASUS scheme is relatively small, and the cell capacityachieved by the aforementioned algorithms slightlygrows with M. In summary, when numbers of the se-lected BS antennas (N) and scheduled users (K) are con-firmed, the increase of the number of transmittingantennas (M) has little effect on the system capacity per-formance. Thus, it can be proven from the result thatactivation of more transmit antennas at the BS side isunnecessary.A cell capacity performance comparison of each
JASUS algorithms with various numbers of the selectedantennas at SINR = 20 dB is shown in Fig. 6. The various
-3 -2 -1 0 1 2 3-3
-2
-1
0
1
2
3
BSUser
Fig. 2 Multi-cell multi-user massive MIMO system with B = 7and K = 10
Table 1 Simulation parameters setting
Parameters Values
No. of cells in the simulation Β 7
Average no. of UEs in one cell 2 ≤ U ≤ 50
Number of BS antennas M 10 ≤M ≤ 64
Inter-site distance 500m
Cell radius 295 m
Path loss 128.1 + 37.6×log10(distance(km))[dB]
Each BS transmission power 10 dB
Shadowing standard derivation 7 dB
Noise spectral density − 174 dBm/Hz
Users’ speed 0
System bandwidth 20 MHz
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 9 of 14
numbers of the selected antennas (12 ≤N ≤ 24) corres-pond to the different maximum cell capacity of the net-works. Compared with the norm-based JASUS, theJASUS algorithms, which are ES, AMCMC, and greedy,had more significant enhancement for the system per-formance. For example, the maximum cell capacity en-hancements, which are approximately 2.5, 2.2, and 1.7 b/s/Hz at SINR = 20 dB, are achieved when 18 (N = 18) BSantennas are selected. From the figure, it can be ob-served, when N > 18, the cell capacity achieved by theJASUS algorithms is slightly growing when the numberof selected antenna goes large. Thus, it can be provenfrom the result that when numbers of the scheduledusers K are confirmed, the system capacity sequentiallyincreases until the numbers of selected BS antenna closeto the N = 18, when N > 18, the increasing number of se-lected antennas has no significant effect on the system
capacity performance. Therefore, the results show thatmore antenna selection is unnecessary at the base sta-tion. Thus, the proposed algorithm is demonstrated tobe effective. In addition, we considerably decreased sys-tem cost and power consumption while approach themaximum cell capacity by selected suitable transmit an-tennas at BS side.Figure 7 shows that the cell capacity different K at
SINR = 20 dB, for user scheduling and with transmitantenna selection. It can be observed that the cellsum capacity increases with increasing of user K. Thedifferent numbers of the scheduled users (2 ≤ K ≤ 16)correspond to the different maximum cell capacity ofthe networks. Compared with the norm-based result,ES, AMCMC, and greedy algorithms enhanced thesystem performance. For example, the maximum cellcapacity enhancement, which are approximately 6.2,
Fig. 3 Ergodic capacity versus SINR with B = 7, M = 64, U = 50, N = 16, and K = 10
Fig. 4 Ergodic capacity achieved in iterative by AMCMC-based JASUS with B = 7, M = 64, U = 50, N = 16, K = 10, and SINR = 20 dB
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 10 of 14
5.5, and 4.0 b/s/Hz at SINR = 20 dB, are achievedwhen ten (K = 10) users are scheduled. From the fig-ure, it can be observed that when K > 10, the cellcapacity achieved by the aforementioned algorithmsslightly grows when the number of the scheduled usergoes large. This result proves that the behavior ofJASUS algorithm does not change drastically whenthe scheduled user number becomes large.
Finally, we discuss on SER of different JASUS algo-rithms. A 16-QAM scheme is used with a ZF receiver.The SER performance of the linear ZF receiver system isshown in Fig. 8. Compared with the results of norm-based JASUS algorithm, ES, AMMC, and greedy algo-rithms improved the SER of the system. Same to thecase of the cell capacity performance, the system SER ofthe AMCMC-based JASUS algorithm is close to that of
Fig. 5 Ergodic capacity of different algorithms with versus transmit antenna numbers with B = 7, U = 50, N = 16, K = 10, and SINR = 20 dB
Fig. 6 Ergodic capacity of different algorithms with different numbers of the selected antennas with B = 7, M = 64, U = 50, K = 10,and SINR = 20 dB
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 11 of 14
the ES-based JASUS algorithm. At the same time, wefind that ES-based JASUS and AMMC-based JASUShave better SER performance than greedy-based JASUSand norm-based JASUS, especially the SINR to high.When the SINR is 20 dB, the SER performances of theES-based JASUS, AMCMC-based JASUS, greedy-basedJASUS, and norm-based JASUS algorithms are approxi-mately 3.5 × 10−2, 3.9 × 10−2, 4.9 × 10−2, and 6.8 × 10−2,respectively.
5.2 Computational complexity analysisThe computational complexities of the introduced differ-ent JASUS algorithm are analyzed in this section. Table 2summarizes the computational complexity of our pro-posed algorithm along with the complexity of another al-gorithm. The asymptotic notations, which reflect thecomputational complexity, was used to evaluation howthe scheme responds to changes of parameters which areM, N, U, and K. From Table 2, we can easily observe the
Fig. 7 Ergodic capacity of different algorithms with different scheduled number of users with B = 7, M = 64, U = 50, N = 16, and SINR = 20 dB
Fig. 8 SER versus SINR with B = 7, M = 64, U = 50, N = 16, and K = 10
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 12 of 14
computational complexity gap between the four methods.Note that CN
M denotes the binomial coefficient, and thematrix inverse operation [29–32] makes the computa-tional complexity of per sample up to O(N3). Thus, theoverall complexity of our proposed Algorithm 1 for theproblem of (11) is O(N3tLMCMC), where t × LMCMC is thetotal number of target function evaluations. Ultimately,we can be observed from Table 2 that our proposed algo-rithm has a very low computational complexity comparedto ES-based JASUS and greedy-based JASUS algorithms.However, the norm-based JASUS algorithm has a very lowcomputational complexity compared to our proposed al-gorithm, but it also has a very low cell capacity. This resultshows that our proposed algorithm is suitable for practicalmulti-cell multi-user massive MIMO system.
6 ConclusionIn this paper, we studied the problem of JASUS in amulti-cell multi-user massive MIMO downlink systemoperating with TDD mode. Considering the trade-off be-tween network performance and computational com-plexity, we proposed a low-complexity algorithm forJASUS method based on AMCMC algorithm in thedownlink multi-cell multi-user massive MIMO systems.AMCMC algorithm has been proven helpful for select-ing combination subset of antennas and users to ap-proach the maximum cell capacity with consideration ofthe inter-cell interference. In our algorithm, the updatingrules of the selection probability of each base station an-tenna and scheduling probability for each user are pro-posed. In addition, we proposed a new projection strategyto satisfy the constraints of selection. Performance analysisand simulation results show that our proposed algorithmcan produce promising results and achieve a good trade-offbetween complexity and performance. Compared with ES-based JASUS algorithm, the proposed algorithm achievedcomparable performance with very low complexity. Inaddition, we demonstrate that our proposed algorithm
outperforms greedy-based JASUS and norm-based JASUSmethods in terms of cell capacity and SER performancewith under poorly conditioned channels. At the sametime, the computational complexity is reduced signifi-cantly by combining with the proposed algorithm.
AcknowledgementsThis work described in this paper was supported by the National Scienceand Technology Major Project: No. 2018ZX03001029-004.
Authors’ contributionsWG conceived and designed the study. SM and KZ performed the simulationexperiments. SM and KZ wrote the paper. XL and ZS reviewed and editedthe manuscript. All authors read and approved the final manuscript.
Authors’ informationSaidiwaerdi Maimaiti received the M.Sc. degree in signal and informationprocessing from Southwest Jiaotong University, Chengdu, China, in 2014. Heis currently working towards the Ph.D. degree in Information andCommunications Engineering, Key Laboratory of Universal WirelessCommunications, Ministry of Education, Beijing University of Posts andTelecommunications Beijing, China. His research interests include massiveMIMO, interference management, radio network planning, resourcemanagement, and intelligent network optimization in 5G network systems.
FundingThe funding for the research reported is provided by the National Scienceand Technology Major Project: No. 2018ZX03001029-004. The funds aremainly used for simulation hardware support.
Competing interestsThe authors declare that they have no competing interests.
Received: 10 March 2019 Accepted: 25 July 2019
References1. T.A. Sheikh, J. Bora, A. Hussain, in CCTCEEC. A survey of antenna and user
scheduling techniques for massive MIMO-5G wireless system (Mysore, 2017),pp. 578–583. https://doi.org/10.1109/CTCEEC.2017.8455177
2. T.L. Marzetta, Noncooperative cellular wireless with unlimited numbers ofbase station antennas. IEEE Trans. Wirel. Commun. 9(11), 3590–3600 (2010)
3. E.G. Larsson, O. Edfors, F. Tufvesson, T.L. Marzetta, Massive MIMO for nextgeneration wireless systems. IEEE Comun. Mag. 52(2), 186–195 (2014)
4. S. Han, I. Chih-Lin, Z. Xu, C. Rowell, Large-scale antenna systems with hybridanalog and digital beamforming for millimeter wave 5G. IEEE Comun. Mag.53(1), 186–194 (2015)
5. Y. Dong, Y. Tang, K.Z. Shenzhen, in IEEE/ACIS International Conference onComputer and Information Science IEEE. Improved joint antenna selectionand user scheduling for massive MIMO systems (2017)
6. M. Torabi, D. Haccoun, Performance analysis of joint user scheduling andantenna selection over MIMO fading channels. IEEE Signal Process. Lett.18(4), 235–238 (2011)
7. C. Bandeira, D. Moreira, P. Normando, et al., in Proc. IEEE Vehicular TechnologyConf. (VTC-Spring). Performance of joint antenna and user selection schemes forinterference alignment in MIMO interference channels (2013)
8. S. Xue, Z. Qian, W. Shao, et al., in Proc. IEEE International Conf. ComputationalElectromagnetics (ICCEM). Distributed joint antenna and user selection forMIMO broadcast channels (2016)
9. J.-C. Chen, “Joint Antenna Selection and User Scheduling for MassiveMultiuser MIMO Systems With Low-Resolution ADCs,” IEEE Trans. Veh.Technol. 68(1) 1019–1024 (2019)
Table 2 Computational complexity analysis
Algorithms General case M = 64, N = 16,U = 50, k = 10
ES-based JASUS OðCNMCKUN
3Þ 2.1 × 1028
AMCMC-based JASUS O(N3tLMCMC) 2.8 × 107
Greedy-based JASUS tO(MUN3) 3.9 × 108
Norm-based JASUS O(MU) + tO(N3) 1.3 × 105
Number of iterations t 1 ≤ t ≤ 50 t = 30
Constant δ 0.5/1/1.5/2/2.5 δ = 2
Number of promisingsamples LMCMC
LMCMC = δ × (M + U) LMCMC = 228
Notes: δ = 2 is the best value in our experiments, because when δ is largerthan 2, the cell capacity performance improvement is no longer significant
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 13 of 14
10. X. Gao, L. Dai, A.M. Sayeed, Low RF-complexity technologies to enablemillimeter-wave MIMO with large antenna array for 5G wirelesscommunications. IEEE Commun. Mag. 56(4), 211–217 (2018)
11. T.-H. Tai, W.-H. Chung, T.-S. Lee, "A low complexity antenna selectionalgorithm for energy efficiency in massive MIMO systems", Proc. IEEE Int.Conf. Data Sci. Data Intensive Syst. 284-289 (2015)
12. X. Gao , O. Edfors, F. Tufvesson, and E. G. Larsson, “Massive MIMO in realpropagation environments: Do all antennas contribute equally?,” IEEE Trans.Commun. 63(11) 3917–3928 (2015)
13. B. Makki, A. Ide, T. Svensson, T. Eriksson, M.-S. Alouini, A geneticalgorithm-based antenna selection approach for large-but-finite MIMOnetworks. IEEE Trans. Veh. Technol. 66(7), 6591–6595 (2016)
14. M. Benmimoune, E. Driouch, W. Ajib, D. Massicotte, in Proc. IEEE WCNC. Jointtransmit antenna selection and user scheduling for massive MIMO systems(2015), pp. 381–386
15. M. Olyaee, M. Eslami, J. Haghighat, An energy-efficient joint antennaand user selection algorithm for multi-user massive MIMO downlink.IET Commun. 12(3), 255–260 (2018)
16. G. Xu, A. Liu, W. Jiang, H. Xiang, W. Luo, Joint user scheduling and antennaselection in distributed massive MIMO systems with limited backhaulcapacity. China Commun. 11(5), 17–30 (2014)
17. B.M. Lee, J. Choi, J. Bang, B.-C. Kang, in IEEE Int. Symp. On Circuits andSystems (ISCAS). An energy efficient antenna selection for large scale greenMIMO systems (2013), pp. 950–953
18. Q.H. Spencer, A.L. Swindlehurst, M. Haardt, Zero-forcing methods fordownlink spatial multiplexing in multiuser MIMO channels. IEEE Trans. SignalProcess. 52(2), 461–471 (2004)
19. G. Dimic, N.D. Sidiropoulos, On downlink beamforming with greedy userselection: performance analysis and a simple new algorithm. IEEE Trans.Signal Process. 53(10), 3857–3868 (2005)
20. M. Al-Saedy, M. Al-Imari, M. Al-Shuraifi, et al. Joint User Selection and Multi-Mode Scheduling in Multicell MIMO Cellular Networks[J]. IEEE Transactionson Vehicular Technology, 1-1 (2017)
21. J.C. Spall, Estimation via Markov chain Monte Carlo. Control Syst. IEEE23(2), 34–45 (2003)
22. W.R. Gilks, S. Richardson, D.J. Spiegelhalter, Markov Chain Monte Carlo inPractice (Chapman and Hall, London, 1996)
23. J. Liu, Monte Carlo Strategies in Scientific Computing (Springer-Verlag,New York, 2001)
24. S. Kullback, R. Leibler, On information and sufficiency. Ann. Math. Stat.22(1), 79–86 (1951)
25. R.Y. Rubinstein, D.P. Kroese, The Cross-Entropy Method: A Unified Approach toCombinatorial Optimization, Monte-Carlo Simulation and Machine Learning(Springer-Verlag, Berlin, 2004)
26. H. Robbins, S. Monro, A stochastic approximation method. Ann. Math. Stat.22(3), 400–407 (1951)
27. H. Kushner, G. Yin, Stochastic Approximation Algorithms and Applications,2nd edn. (Springer-Verlag, New York, 2003)
28. Y. Liu, Y.Y. Zhang, C.L. Ji, W.Q. Malik, D.J. Edwards, A low complexityreceive-antenna-selection algorithm for MIMO-OFDM wireless systems.IEEE Trans. Veh. Technol. 58(6), 2793–2802 (2009)
29. X. Gao, L. Dai, Y. Ma, Z. Wang, Low-complexity near-optimal signal detectionfor uplink large-scale MIMO systems. Electron. Lett. 50(18), 1326–1328 (2014)
30. X. Gao, L. Dai, Y. Hu, Y. Zhang, Z. Wang, Low-complexity signal detection forlarge-scale MIMO in optical wireless communications. IEEE J. Sel. AreasCommun. 33(9), 1903–1912 (2015)
31. L. Dai, X. Gao, X. Su, S. Han, I. Chih-Lin, Z. Wang, Low-complexity soft-outputsignal detection based on Gauss-Seidel method for uplink multiuser large-scale MIMO systems. IEEE Trans. Veh. Technol. 64(10), 4839–4845 (2015)
32. Z. Wu, C. Zhang, Y. Xue, S. Xu, Z. You, "Efficient architecture for soft-outputmassive MIMO detection with gauss-seidel method", Proc. IEEE Int. Conf.Circuits Syst. (ISCAS), 1886-1889 (2016)
Publisher’s NoteSpringer Nature remains neutral with regard to jurisdictional claims inpublished maps and institutional affiliations.
Maimaiti et al. EURASIP Journal on Wireless Communications and Networking (2019) 2019:208 Page 14 of 14