Distributed Pilot Assignment For Distributed Massive-MIMO Networks
Distributed Pilot Assignment For Distributed Massive-MIMO Networks
Distributed Pilot Assignment For Distributed Massive-MIMO Networks
Abstract
Pilot contamination is a critical issue in distributed massive MIMO networks, where the reuse
of pilot sequences due to limited availability of orthogonal pilots for channel estimation leads to
performance degradation. In this work, we propose a novel distributed pilot assignment scheme to
effectively mitigate the impact of pilot contamination. Our proposed scheme not only reduces latency,
but it also enhances fault-tolerance. Extensive numerical simulations are conducted to evaluate the
performance of the proposed scheme. Our results establish that the proposed scheme outperforms
existing centralized and distributed schemes in terms of mitigating pilot contamination and significantly
enhancing network throughput.
Index Terms
I. I NTRODUCTION
The rapid growth and adoption of wireless communication based data services has prompted
the development of technologies that can improve the capacity, reliability, and efficiency of cellu-
lar networks. By reaping the benefits of mMIMO and Ultra Dense Network (UDN), Distributed
mMIMO is emerging as a promising technology [1], [2]. It abandons the concept of distinct cells
but instead deploys a large number of distributed access points (APs) throughout the coverage
area [1]. As there are no traditional cells, the inter-cell interference that plagues cellular systems
is eliminated. Further, as the user equipments (UEs) are very close to APs, this may provide high
coverage probability in the network, thus improving quality-of service to the users. However, as
in traditional mMIMO systems, channel estimation is a critical challenge in distributed mMIMO
2
networks. One commonly employed approach is blind convolution techniques, such as matrix
decomposition-based signal detection. However, these methods often have high computational
complexity. To mitigate this complexity, researchers have explored leveraging the channel hard-
ening property of mMIMO systems [3]. However, in distributed mMIMO, the channel hardening
is not always guaranteed [4], so getting the estimates using matrix decomposition-based channel
estimation is not a viable option.
Pilot-based channel estimation presents a simple and low-complexity approach that can be
employed within the distributed mMIMO systems to obtain channel estimates. However, before
deploying such solutions in the networks an associated issue of pilot contamination needs to
be addressed. Pilot contamination arises when the same pilot sequences are used for estimating
channels of more than one user, leading to estimation errors and subsequent performance degra-
dation. In a practical system, as it is almost impossible to assign orthogonal pilots to all users,
the users may reuse the pilots. To minimize the impact of pilot contamination, pilot assignment
(PA) should be carried out so that it minimizes the contamination. A lot of work on pilot
assignment in distributed mMIMO networks exists, but most of it is centralized in nature [1],
[5]–[11]. In [1], the authors have proposed two pilot assignment schemes namely random pilot
assignment and greedy pilot assignment based on throughput improvement. In [5], the authors
have proposed a scalable scheme which performs joint pilot assignment and AP-UE association.
In [6], the authors have proposed a graph coloring-based pilot assignment, where the AP-UE
association takes place initially and then pilots are assigned using graph coloring. If the pilot
assignment fails, the AP-UE associations are updated, and so is the pilot assignment. In [7],
the authors have proposed a pilot assignment scheme based on Hungarian matching algorithm.
In [8], the authors have proposed a clustering based scalable pilot scheme where UEs in the
same cluster are allocated the same pilots. In [9], the authors have again exploited the graph
theory, but instead by graph coloring, they have used weighted graph approach. In [10], the
authors have tried to exploit the interference-aware reward, calculated by treating noise as a
interference, for the pilot assignment and AP-UE association jointly. Along with random pilot
assignment [1], that can also be implemented in distributed manner, there are a few distributed
schemes [12], [13]. In [12], the authors have proposed a survey propagation based distributed
pilot assignment scheme that incurs high signaling overhead due to messages passing among
the feasible groups. Also, for the proposed scheme the computational cost increases drastically
as the number of users-to-pilots ratio increases. In [13], the authors have proposed a distributed
3
where τt is the uplink transmit power of the UE t ∈ Tp and nm ∼NC (0, σ 2 IA ) denotes the noise
at the mth AP. The MMSE estimate of channel ht,m is given by
Lp τt Rt,m Ψ−1
p
ĥt,m = t,m yp,m , (2)
p
Lp τt Rt,m + σ 2 IA is the correlation matrix of yp,m . The pilot sharing not
P
where Ψt,m = t∈Tp
only leads to poorer channel estimation, but also affects ĥt,m , which becomes more correlated
and thus increases the interference among the UEs.
The association of the UE t and the AP m is denoted by a indicator variable dt,m , which
equals ’1’ when the UE t is served by the AP m, and ’0’ when the UE t is not so served. The
received uplink signal at the AP m is given by
T
X
yum = τt ht,m dt,m xt + nm . (3)
t=1
The payload signal xt ∈ C is the unit power complex signal sent by the UE t and the power
associated with it is denoted by τt . The combining vector vt,m ∈ CA is assigned to the UE t by
the AP m. The estimate of xt is x̂t = t∈T dt,m vH u
P
t,m ym .
Lp
SEut = (1 − )E{log2 (1 + SINRut )}, (4)
Lc
where SINRut represents the uplink instantaneous signal to noise plus interference ratio of the
5
PM PT
where zt = m=1 dt,m vH
t,m ( t=1 τt Ct + σ 2 IM A )vt,m , and Ct is the error correlation matrix for
the collective channel of the UE t.
The precoding vector, st,m ∈CA , assigned to the tth UE by the mth AP and the received
downlink signal at the tth UE respectively, are:
√ ĥt,m
st,m = τt,m q , (6)
E||ĥt,m || 2
M
X T
X
ytd = hH
t,m dj,m sj,m xj + nt , (7)
m=1 j=1
where nt ∼NC (0, σ 2 ) is the thermal noise, xj ∈ C is unit power complex signal sent to the tth
UE and τt,m is the downlink power assigned to the tth UE by the mth AP, allocated using the
fraction power allocation [5].
By utilizing the use-and-then-forget bound [15], the downlink spectral efficiency (SE) is
Lp
SEdt = (1 − ) log2 (1 + SINRdt ), (8)
Lc
where SINRdt represents the downlink instantaneous signal to noise plus interference ratio of the
tth UE and is given by
M
dt,m hH 2
P
|E{ t,m st,m }|
m=1
SINRdt = M
, (9)
dt,m hH
P
z − |E{ t,m st,m }|
2 + σ2
m=1
PT PM
where z = j=1 |E{ m=1 dj,mhH 2
j,m sj,m }| .
We propose a distributed pilot assignment scheme to address the limitations on the number
of mutually orthogonal pilots (Lp <T ) and mitigate pilot contamination. A key prerequisite for
distributed mMIMO is that each UE be served by at least one AP. Additionally, each AP must
6
assign a given pilot to at most one UE to minimize pilot contamination. To guarantee that each
UE is served by at least one AP and to manage pilot assignment, a specific AP is designated
as the controller-AP for each UE. Therefore, in the proposed scheme first, a controller-AP is
assigned to each UE in a distributed manner. Then, pilots are assigned distributively to each
UE by its associated controller-AP. Furthermore, to fully leverage the benefits of mMIMO and
make the system scalable, each UE needs to be connected to an adequate number of APs. To
achieve this, a distributed AP-UE clustering scheme is employed in the final step based on the
constraint that each AP serves one UE per pilot1, thus minimizing pilot contamination while
providing significant power to the UEs.
For each AP m, we define a subset Um of UEs that includes, at most, the leading Lp UEs of
the mth AP, as determined by the decreasing order of the parameter βt,m . Let Tmc be the subset
of UEs with the AP m as their controller-AP. We define Lpt as the set of available pilots for UE
t, initially contains all Lp pilots.
A. Controller-AP assignment
For the controller-AP assignment, we propose a distributed algorithm. In the proposed ap-
proach, each UE is responsible for selecting its controller-AP from the set of available APs
using their respective LSFC values, by following the standard random access procedure [16]
with the broadcast signal [5]. However, there is a challenge when more than Lp UEs intend to
be allocated a particular AP as their controller-AP. This challenges our design requirement that
each AP can serve at most one UE per pilot. However, another issue is that, if a UE requests an
AP with a lower LSFC value to be its controller-AP, this could potentially lead to degradation
in the system performance. To address these, we propose a Controller-AP selection algorithm,
where each UE is assigned one AP from the set Mct of APs as its controller-AP. Also, an AP
can serve as the controller-AP of at most Lp UEs. The details of constructing subset Mct are
given in Algorithm 1. In order to identify the APs that are serving more than Lp UEs and require
some of their associated UEs to request a new controller-AP assignment, and the APs that are
just serving Lp UEs and so are unable to accept controller-AP requests for more UEs, the APs
are categorized into two set. The former are called the oversaturated-APs and are denoted as
MD , and the latter are called the inert-APs and are denoted as MI , respectively. The detailed
1
The number of users can grow but up to a certain extent such that T < M A, for the system to behave as mMIMO [1].
Therefore, the above assumption does not hinder the performance when numbers of users is very large.
7
procedure to select controller-AP distributively from the set of potential controller-APs is given
by Algorithm 1. The Algorithm 1 is carried out at the APs in cooperation with the UEs.
Now we propose a distributed pilot assignment scheme. In order to minimize pilot contami-
nation, we introduce a contamination matrix Ad, a binary indicator matrix, indicating whether
any two UEs can share a pilot or not. Every AP assigns distinct pilots to its serving UEs and
informs2 its neighboring APs in Mnb nb
m about its pilot assignment. The construction of Mm is
explained in Algorithm 2. If any two UEs have the same pilot and Ad indicates that these two
UEs cannot share the pilot, then the AP allows a UE to retain the pilot if either the cardinality
of set Lpt for this UE is smaller or its index is smaller than the other UE. The detailed procedure
for assigning pilots to UEs is outlined in Algorithm 2. The Algorithm 2 is carried out at the
APs.
C. UE-AP Clustering
The distributed UE-AP clustering is performed to improve the spectral diversity, as in [5].
Each AP serves at most one UE per pilot in order to minimize the pilot contamination. Initially,
each AP m chooses to serve each UE t ∈ Tmc . Then the AP m, for each pilot (except pilots of
UEs in Tmc ), chooses a UE to serve based on the largest LSFC.
D. Complexity analysis
on execution of the while loop between lines 17-36. Assuming the while loop runs I times, the
complexity of execution of the while loop is O(IML2p ). Therefore, the overall complexity of
Algorithm 2 is O(MLp + IML2p ). The final complexity of overall procedure including AP-UE
association is O(TM2 + MLp + IML2p + MLp ). Putting everything together, the overall complexity
of our proposed scheme is O(TM2 + IML2p ). The actual execution time of the aforementioned
2
The exchange of information about pilots among APs is carried out via the CPU.
9
processes may be substantially lower when compared to the complexities given above, primarily
due to opportunities for extensive parallel execution of various processes within the procedure.
In Algorithm 1, the formation of the sets MD and MI may take place at the APs independently.
Similarly, each AP in MD may run the while loop in Algorithm 1 in parallel. This may further
reduce the execution time significantly. In Algorithm 2, the pilot assignment is being carried out
in parallel at the APs, thus reducing the execution time.
For numerical simulations we have considered an area of 2000 × 2000 square meters in which
T UEs and 100 APs are uniformly distributed. Each AP is equipped with 4 antennas. We have
considered APs to be 10 meters higher than the UEs. We set Lc = 200, Lp = 10, τt =
100 mW and 1000 mW downlink power of each AP. We have considered bandwidth to be
20 MHz. To simulate the large scale propagation model, we have considered the 3GPP Urban
Microcell scenario and the local scattering model for spatial correlation same as in [17]. We
have averaged our results over 50 network instances and for each instance we have considered
500 channel realizations. For all the simulations, P-MMSE precoder [5] is used. For spectral
efficiency comparison we have considered centralized schemes, such as scalable PA [5], IAR-
sum PA [10], graph-coloring-based PA [6]; and a distributed scheme: survey-propagation-based
(survey) PA [12]. For a fair comparison, in scalable and survey-propagation, a AP can serve at
most Lp UEs.
A. Downlink Operations
1) Comparison with centralized schemes: Fig. 1 shows the plot of the 90%-likely downlink
SE for different numbers of UEs. The proposed distributed scheme outperforms the scalable
PA scheme by 2%, the graph-coloring PA scheme by 2%, and the IAR-sum PA scheme by 6%
in scenario involving 100 UEs. Similarly, improvements with 150 UEs are 3%, 18%, and 9%
respectively, and improvements with 200 UEs are 7%, 36%, and 7%, respectively. It can also
be observed that as the number of UEs increases, the 90%-likely downlink SE of all schemes
decreases, but our proposed scheme sees a lesser decrease than all other schemes. These results
demonstrates the superiority of our proposed scheme in terms of 90%-likely downlink SE than
all other schemes for all user scenarios.
11
1.5
0
100 150 200
Similarly Fig. 2 shows the graph between average downlink SE for different numbers of UEs.
For T = 100, our proposed scheme demonstrates superior performance compared to the graph-
coloring PA and IAR-sum PA by 22% and 1% respectively. However, it is slightly outperformed
(by 1%) by the scalable PA. For T = 150, the proposed scheme outperforms both the graph-
coloring based PA and the IAR-sum PA, while achieving comparable performance to the scalable
PA. Finally, for T = 200, the proposed scheme surpasses all other schemes.
2) Comparison with a distributed scheme: Fig. 3 shows the CDF of 90%-likely downlink
SE for distributed PA schemes. To decrease the computational cost of the distributed survey
PA scheme, we have considered the number of UEs to the number of pilots ratio to be 3. The
plot shows that the proposed scheme outperforms the survey-propagation-based scheme, thus
demonstrating its superiority.
12
1
0.12
0.8 0.1
0.6 0.08
1.4 1.6 1.8 2
0.4
0.2
0
0 1 2 3 4 5 6 7 8 9 10
T
Fig. 3. CDF of downlink SE, when Lp
= 3.
B. Uplink Operations
1) Comparison with centralized schemes: Fig. 4 shows the plot of the 90%-likely uplink
spectral efficiency for different numbers of UEs. In a scenario involving 100 UEs, the proposed
approach outperforms the scalable PA technique by 1%, the graph-coloring-based PA scheme by
36%, and the IAR-sum PA scheme by 4%. Similarly, with 150 UEs, improvements are 3%, 72%,
and 6%, respectively, and improvements with 200 UEs are 11%, 105%, and 9%, respectively.
It can also be seen that when the number of UE increases, the 90%-likely downlink SE of all
schemes decreases, however our proposed scheme has a lesser decrease than all other schemes.
1.5
1.17 1.16 1.13
1 0.86
0.67 0.65 0.63
0
100 150 200
Fig. 5 shows the graph between average uplink SE for different numbers of UEs. For T = 100,
our proposed scheme demonstrates superior performance compared to the graph-coloring PA by
approximately 17% and the IAR-sum PA by nearly 2%. It is, however, outperformed (by ≤ 1%)
13
by the scalable PA. The proposed scheme surpasses both the graph-coloring-based PA and the
IAR-sum PA by a significant margin for T = 150, while matching the scalable PA performance.
Finally, for T = 200, our proposed scheme surpasses all other schemes.
4 4.04 3.94
4
3.41
3.25 3.26 3.17
3 2.7 2.68
2.58 2.63
2.04
2
0
100 150 200
2) Comparison with a distributed scheme: Fig. 6 shows the uplink SE plot for distributed PA
schemes. Again in order to decrease the computational costs in distributed survey PA scheme, we
have considered the number of UEs to the number of pilots ratio to be 3. In terms of 90%-likely
SE, the proposed scheme outperforms the survey-propagation-based scheme by 2% for T = 45
and by nearly 4% for T = 60. Thus, indicating the superiority of our proposed scheme.
1
0.12
0.8 0.1
0.6 0.08
1.8 2 2.2 2.4
0.4
0.2
0
0 1 2 3 4 5 6 7 8 9 10
T
Fig. 6. CDF of uplink SE, when Lp
= 3.
C. Discussion
Superior performance of our proposed distributed scheme in terms of both 90%-likely downlink
as well as uplink SE demonstrates its ability to outperform not only distributed survey PA scheme,
14
but all of the aforementioned centralized schemes, and performance continues to improve as the
number of UE grows. This is primarily attributable to our approach’s priority of reducing pilot
contamination for each UE by assigning distinct pilots to UEs that are more susceptible to
contamination, thus improving fairness among UEs in terms of contamination. Although our
proposed scheme does not focus on improving the average SE, it still achieves greater average
downlink as well as uplink SE than IAR-sum and graph-coloring PA schemes in all scenarios.
However, our approach has a modest disadvantage in terms of both average downlink and uplink
SEs over the scalable PA scheme for low user density scenarios. This minor difference in average
SE demonstrates that, despite being distributed, our approach is competitive with centralized
approaches.
V. C ONCLUSION
Our work shows that the proposed distributed pilot assignment technique in distributed mMIMO
systems may substantially enhance spectral efficiency compared to centralized as well as dis-
tributed existing schemes. By using a distributed approach, this strategy not only boosts spectral
efficiency but also reduces latency and increases the network’s fault-tolerance capability. These
findings highlight the importance of using distributed resource allocation algorithms to achieve
optimal performance in distributed mMIMO systems. The positive findings of this study illustrate
the potential benefits of adopting distributed approaches in future distributed mMIMO network
plans and deployments.
R EFERENCES
[1] H. Q. Ngo, A. Ashikhmin, H. Yang, E. G. Larsson, and T. L. Marzetta, “Cell-free massive MIMO versus small cells,”
IEEE Trans. on Wireless Commun., vol. 16, no. 3, pp. 1834–1850, 2017.
[2] Ö. T. Demir, E. Björnson, L. Sanguinetti et al., “Foundations of user-centric cell-free massive MIMO,” Foundations and
Trends® in Signal Processing, vol. 14, no. 3-4, pp. 162–472, 2021.
[3] P. Aswathylakshmi and R. K. Ganti, “Pilotless uplink for massive MIMO systems,” arXiv preprint arXiv:2305.12431,
2023.
[4] Z. Chen and E. Björnson, “Channel hardening and favorable propagation in cell-free massive MIMO with stochastic
geometry,” IEEE Trans. on Commun., vol. 66, no. 11, pp. 5205–5219, 2018.
[5] E. Björnson and L. Sanguinetti, “Scalable cell-free massive MIMO systems,” IEEE Trans. on Commun., vol. 68, no. 7,
pp. 4247–4261, 2020.
[6] H. Liu, J. Zhang, S. Jin, and B. Ai, “Graph coloring based pilot assignment for cell-free massive MIMO systems,” IEEE
Trans. on Veh. Techno., vol. 69, no. 8, pp. 9180–9184, 2020.
15
[7] S. Buzzi, C. D’Andrea, M. Fresia, Y.-P. Zhang, and S. Feng, “Pilot assignment in cell-free massive MIMO based on the
hungarian algorithm,” IEEE Wireless Commun. Letters, vol. 10, no. 1, pp. 34–37, 2020.
[8] S. Chen, J. Zhang, E. Björnson, J. Zhang, and B. Ai, “Structured massive access for scalable cell-free massive MIMO
systems,” IEEE Journal on Selected Areas in Commun., vol. 39, no. 4, pp. 1086–1100, 2020.
[9] W. Zeng, Y. He, B. Li, and S. Wang, “Pilot assignment for cell free massive MIMO systems using a weighted graphic
framework,” IEEE Trans. on Veh. Techno., vol. 70, no. 6, pp. 6190–6194, 2021.
[10] S. Chen, J. Zhang, E. Björnson, and B. Ai, “Improving fairness for cell-free massive MIMO through interference-aware
massive access,” IEEE Trans. on Veh. Techno., 2022.
[11] H. Yu, X. Yi, and G. Caire, “Topological pilot assignment in large-scale distributed MIMO networks,” IEEE Trans. on
Wireless Commun., vol. 21, no. 8, pp. 6141–6155, 2022.
[12] S. Kim, H. K. Kim, and S. H. Lee, “Survey propagation for cell-free massive MIMO pilot assignment,” in 2022 IEEE GC
Wkshp. IEEE, 2022, pp. 1267–1272.
[13] M. Rahmani, M. J. Dehghani, P. Xiao, M. Bashar, and M. Debbah, “Multi-agent reinforcement learning-based pilot
assignment for cell-free massive MIMO systems,” IEEE Access, vol. 10, pp. 120 492–120 502, 2022.
[14] H. Q. Ngo, L.-N. Tran, T. Q. Duong, M. Matthaiou, and E. G. Larsson, “On the total energy efficiency of cell-free massive
MIMO,” IEEE Trans. on Green Commun. and Network., vol. 2, no. 1, pp. 25–39, 2017.
[15] E. G. Larsson, O. Edfors, F. Tufvesson, and T. L. Marzetta, “Massive MIMO for next generation wireless systems,” IEEE
Commun. Magazine, vol. 52, no. 2, pp. 186–195, 2014.
[16] S. Sesia, I. Toufik, and M. Baker, LTE-the UMTS long term evolution: from theory to practice. John Wiley & Sons,
2011.
[17] E. Björnson, J. Hoydis, L. Sanguinetti et al., “Massive MIMO networks: Spectral, energy, and hardware efficiency,”
Foundations and Trends® in Signal Processing, vol. 11, no. 3-4, pp. 154–655, 2017.