Performance Analysis of Hybrid Cellular and Cell-free MIMO Network
Abstract
Cell-free wireless communication is envisioned as one of the most promising network architectures, which can achieve stable and uniform communication performance while improving the system energy and spectrum efficiency. The deployment of cell-free networks is envisioned to be a long-term evolutionary process, in which cell-free access points (APs) will be gradually introduced into the communication network and collaborate with the existing cellular base stations (BSs). To further explore the performance limits of hybrid cellular and cell-free networks, this paper develops a hybrid network model based on stochastic geometric toolkits, which reveals the coupling of the signal and interference from both the cellular and cell-free networks. Specifically, the conjugate beamforming is applied in hybrid cellular and cell-free networks, which enables user equipment (UE) to benefit from both cellular BSs and cell-free APs. The aggregate signal received from the hybrid network is approximated via moment matching, and coverage probability is characterized by deriving the Laplace transform of the interference. The analysis of signal strength and coverage probability is verified by extensive simulations.
I Introduction
The future wireless communication networks will witness a proliferation of mobile applications and unprecedented growth in wireless data. The realization of higher spectrum and energy efficiency with superior costs remains a challenging issue in current research. As one of the most prominent wireless technologies proposed in recent years, cell-free network is considered as a promising network architecture in the beyond fifth-generation (B5G) and sixth generation (6G) mobile communication system [1]. Different from traditional cellular systems, cell-free network is a user-centric coverage architecture that discards the traditional concept of cellular boundaries [2]. The central processing unit (CPU) controls the access points (APs) to cooperate to provide services to user equipment (UE) on the same time-frequency resources, thus realizing higher spatial multiplexing [3, 4]. Cell-free network improves the energy and spectral efficiency of the system, and effectively reduces the performance gaps between UE by ensuring that there are spatially short-range APs that provide stable services to UE.
However, the deployment of cell-free systems in existing commercial mobile networks still faces serious challenges. First, the construction of cell-free systems requires the deployment of the distributed APs throughout the network and the construction of the corresponding fronthaul links, which brings heavy time costs and deployment expenses. Second, simply introducing the cell-free system without cooperation will inevitably cause mutual interference with existing cellular systems, significantly limiting the system performance. Therefore, the deployment of cell-free systems is bound to be a long-term evolutionary process, and hybrid cellular and cell-free cooperation networks are both a necessity and a desirable choice for B5G and 6G.
Some recent works have investigated the performance analysis and resource allocation of hybrid cellular and cell-free networks. A cell-free and legacy cellular coexistence system deployed on the existing system architecture, as well as the corresponding precoding, power control, etc., are outlined in [5]. With appropriate UE association criteria and coordinated beamforming, hybrid cell-free and small cell systems can provide superior downlink rates for static and dynamic UE than the single architecture [6]. However, the above works do not take into account the impact of the spatial distribution of base stations (BSs), APs, and UE on network performance.
Due to the densification and irregularity of wireless node distributions in the network, traditional grid-based deployment models are difficult to reflect the practical system performance. Stochastic geometry models the spatial distribution of wireless nodes with point processes and can effectively characterize the lower bound of the actual system performance. There have been some works using stochastic geometry to analyze the performance of cell-free networks and coordinated multiple points (CoMP) communication in terms of energy efficiency (EE) [7], power control [8] and channel hardening analysis [9], etc. However, there is still a lack of work related to the characterization of hybrid cellular and cell-free networks. In addition, the existing stochastic geometry-based heterogeneous network studies, which separate different network layers from each other [10], are not applicable to the analysis of hybrid cellular and cell-free networks.
To gain some insights of the performance limit of the hybrid network, this paper develops a stochastic geometry-based model for hybrid cellular and cell-free networks, which reveals the coupling of the signal and interference from both the cellular and cell-free networks. However, the aggregate signals from the BSs and APs with conjugate beamforming make it difficult to characterize the distribution of the signal strength and the corresponding signal-to-interference plus noise ratio (SINR). To tackle this issue, we first derive the closed-form expressions for the average signal strength and interference power from the APs, and then the aggregate signal strength distribution is approximated via moment matching. Finally, the coverage probability is characterized on the basis of the Laplace transform of the system interference power. The analysis of network coverage probability is verified by extensive simulations, and it can be used to guide the network deployment and interference management in the hybrid cellular and cell-free networks.
II System Model
As shown in Fig. 1, a hybrid cellular and cell-free network is considered in this paper. The locations of BSs, cell-free APs, and single-antenna UE are modeled by independent homogeneous Poisson point processes (HPPP) , and , with density /, / and /, respectively. Each BS is equipped with antennas, while each AP is equipped with antennas. Considering the different configurations that AP and BS can support, and denote the maximum downlink power of AP and BS with . In the network, and denote the distance between BS and UE , and that between AP and UE , respectively. We consider a typical UE, referred as UE 0, which is jointly served by the closest BS, named BS 0, with the distance , and all the cell-free APs in . The channel vector between BS and UE is denoted by , while the channel vector between AP and UE is . The channel model consisting of distance-dependent large-scale fading and random small-scale fading is considered as
(1) |
(2) |
where and are path loss of the channel with and . and are the sets of all the BSs and APs, respectively. The small-scale fading in both and are independent and identically distributed (i.i.d.) random variables (r.v.s).
The entire network area is represented by .

During the downlink transmission, UE in the same cell will be served by the same BS. In addition, both BSs and APs select conjugate beamforming in order to obtain low computational complexity and good performance, and also to avoid channel state information (CSI) interactions between APs [1]. Therefore, the signal transmitted by BS is
(3) |
where denotes the set of UE that are served by BS , and denotes the information-bearing symbols for UE . denotes the power constraint parameter with . For convenience, is expressed as the average number of users per BS, i.e., , where denotes the average of for any .
Denote the set of UE in the network as . Therefore, the corresponding downlink signal from each cell-free AP is
(4) |
where denotes the power constraint parameter. Note that the averaged number of UE in is , can be denoted as to ensure .
For any UE in the network, is denoted as the index of the associated and nearest BS providing the service. With the collaboration of cell-free APs and cellular BSs, the downlink signal received by the typical UE 0 is
|
(5) |
where the first term and the second term represent the desired signals from BS 0 and cell-free APs, respectively. The total interference is shown in the third term . Each term in consists of the signal sent by cellular BSs and cell-free APs to any other UE. The last term denotes the additive white Gaussian noise (AWGN) with power .
Based on (5), the interference power caused by the signal intended to UE is given by
(6) |
In interference , BS channel vector is independent of channel vector . Meanwhile, the channel gains and from different APs are also independent of each other, . Taking into account the law of large number and the mutual independence of channels as well as beamforming vectors between different BSs and APs, can be approximated as
(7) |
The approximation in (7) indicates that can be expressed in the form of the sum of the power of the signal from the associated BS and each AP. Further, by classifying the interference into intra-cell interference , inter-cell interference , and interference due to the APs, the total interference can be rewritten as
(8) |
where
(9) | ||||
Therefore, the corresponding received SINR at UE 0 is approximately expressed as
(10) |
where
|
(11) |
III Analysis of Signal Strength and Interference
In this section, the statistical distributions of the received signal and interference power are characterized. The signal strength is approximated via moment matching, with its first- and second-order moments derived in closed-form. The intra- and inter-cell interference caused by BSs are approximated as a Gamma r.v. and a weighted sum of Gamma r.v.s, respectively. The average interference caused by cell-free APs is derived in closed-form. In addition, the performance of the network can be obtained by analyzing a typical UE 0 according to Slivnyak’s theorem [11].
III-A Analysis of Channel Distribution
Based on the channel model of BSs and APs, the power of the channel to UE for the th BS and the th AP can be respectively given by
(12) |
(13) |
Since all the entries in both and follow the i.i.d. , and are isotropic vectors in and dimensions respectively [12].
Note that for the isotropic vector with each entry following i.i.d. , is the sum of i.i.d. variables , and thus follows [13]. Therefore, we have and .
Lemma 1.
For the Gamma distributed r.v. and any , [14].
III-B Approximation of the Signal Power Distribution
According to (14), the power of the nearest associated BS channel is the sum of i.i.d. variables following , i.e., . For further analysis of the desired signal in (9), Lemma 2 about the square root of Gamma variable is first introduced.
Lemma 2.
For any Gamma distributed r.v. , the square root of follows the Nakagami distribution as [15], where the parameters are .
Therefore, the distribution of is obtained according to Lemma 2 as
(16) |
while the component of AP channel in has
(17) |
The distribution of the desired signal is composed of the signals from the associated BS 0 together with all APs. Considering the random distribution of the cell-free APs in the network, the following Lemma 3 is introduced.
Lemma 3.
With the law of large number, the desired signal in due to the APs can be approximated by their average when the number of APs is large and , i.e.,
(18) | ||||
where . The detailed derivation is based on the Campbell Theorem [11], and will be shown in an extended journal version.
From Lemma 3, the power expression for the desired signal is simplified as the square of the sum of a Nakagami r.v. and the constant as . Therefore, the following Lemma is introduced.
Lemma 4.
For any Nakagami r.v. , the probability density function (PDF) of the square of the shifted Nakagami r.v. for is
|
(19) |

The verification of the distribution of is shown in Fig. 2. From (19) and Fig. 2, the exact distribution of is difficult to characterize, but the corresponding PDF of has a similar structure to that of the Gamma distribution. Therefore, with given distance of the nearest associated BS 0, can be approximated as a Gamma r.v. based on its first- and second-order moments [16]. The corresponding Lemma (5) is introduced as follows.
Lemma 5.
According to the definition of the Gamma r.v. [17], the desired signal power can be approximated as the Gamma distribution with
(20) | ||||
where the first- and second-order moments of are
(21) |
(22) | ||||
The detailed derivation is based on the raw moments of Gamma distribution, and is omitted here for space.
III-C Analysis of Interference Power Distribution
Similar to the case of the desired signal, the power distribution of the interference is analyzed in this subsection. Considering that conjugate beamforming is applied by both BSs and APs in the network, the Lemma 6 for the projection of isotropic channel vectors is introduced.
Lemma 6.
Denote as an isotropic vector with i.i.d. entries. If is projected onto an s-dimensional beamforming subspace, the power distribution is [18]
(23) |
Based on Lemma 6, the power of intra-cell interference in (9) can be approximated as the sum of i.i.d. variables following . Further, extracting the scale parameter, the power of intra-cell interference can be rewritten as , where .
For the inter-cell interference, since is also independent of , the interference power of each BS in is approximated as . Therefore, as the sum of interference from BSs in , the inter-cell interference can be further expressed as the sum of Gamma variables with the same shape parameters and scale parameters, i.e.,
(24) |
where .
Next, we need to analyze the interference from the APs. According to (3) and (6), the following Lemma about is introduced.
Lemma 7.
By the law of large number, the interference due to the APs is approximated by its average when the number of APs is large and , i.e.,
(25) | ||||
The detailed derivation is based on the Campbell Theorem [11], and will be shown in an extended journal version.
IV Coverage Probability of Hybrid Cellular and Cell-free Network
In this section, the coverage probability of hybrid cellular and cell-free network is analyzed based on the distribution of signal strength and various interference components, derived in the preceding section. In general, the coverage probability is the complementary cumulative distribution function (CCDF) of SINR over the overall network, which can be defined as
(26) |
where denotes the target threshold of the SINR .
IV-A Analysis of Coverage Probability
Taking the distance between the typical UE 0 and its associated and nearest BS 0 as an r.v., the average coverage probability in the network is
(27) |
where denotes the PDF of the nearest point distance in PPP. With the cumulative distribution function (CDF) of r.v. as [11], there is
(28) |
First, the relationship between intra-cell interference and inter-cell interference is analyzed. Clearly, both and are dependent on the distance , i.e., comes from the associated BS 0 with a distance of , and comes from the other BSs with distances . However, there is no interaction between and . Specifically, with the given , depends on the distribution of UE in the cell of BS 0, while depends mainly on the distribution of other BSs with a distance no smaller than . Therefore, and are independent of each other with a given . Therefore, the coverage probability can be further expressed as in Lemma 8.
Lemma 8.
Considering that the desired signal following the Gamma distribution, i.e., the network coverage probability in (26) can be expressed as
(29) | ||||
where the sum of interference from APs and noise is denoted as . The shape parameter is integer. The Laplace transforms of and are
|
(30) |
Based on Lemma 8, the analysis of the coverage probability is transformed into the analysis of the higher-order derivatives, and the coverage of probability is rewritten as
(31) |
where
(32) | ||||
The higher-order derivatives of is derived in the next part.
IV-B Evaluation of Higher-order Derivatives
The objective function of the higher-order derivatives in (31) is rewritten in the form of the exponential function, i.e.,
(33) | ||||
where is applied for convenience.
Since is a composite function of , the special case of Faà di Bruno’s formula with exponential functions can be applied to efficiently derive the th order derivatives of [16, 19], i.e.,
|
(34) |
where denotes the th complete exponential Bell polynomial, whose coefficients can be efficiently obtained according to its definition [20, 21]. The remaining work is to evaluate the higher-order derivatives of , which can be decomposed as
(35) |
For , the derivatives from order 1 to order can be expressed respectively as
(36) |
Additionally, for , there is
(37) |
Finally, the higher-order derivatives of is
(38) | ||||
V Simulation Results
In this section, the analytical results of the coverage probability of the hybrid cell and cell-free network are verified by the comparison with the Monte-Carlo (MC) simulation results. Each result of the MS simulation is averaged from 1000 randomly generated wireless node distributions with 5 realizations per channel. All the wireless nodes are randomly distributed in a circular area of radius 500m and UE 0 is located at the center of the circle. The densities of BSs, APs and UE are , and , respectively. Other relevant parameters are as follows: , , dB, , , , m/s, GHz and .
Fig. 3 shows the coverage probability in cellular networks (), cell-free networks (), and hybrid cellular and cell-free networks for different SINR thresholds . The coverage probability is obtained from the linear weighted probability of the upper and lower integers of , i.e.,
(39) |
From Fig. 3, It is expected that the coverage probability analysis of the hybrid network is generally consistent with the results of MC simulation under different SINR threshold, and can also be applied to the special case where . Compared with traditional cellular networks, hybrid networks effectively improve the communication performance of edge UE and reduce the performance gaps between UE. Compared with cell-free networks, such hybrid networks can achieve higher peak SINR. Therefore, by deploying low-power APs, hybrid cellular and cell-free networks can provide UE with uniformly good communication services while obtaining better peak SINR performance.

VI Conclusion
In this paper, the hybrid cellular and cell-free network is modeled by the stochastic geometry approach, revealing the coupling of the signal and interference from both the cellular and cell-free networks. Moment matching is used to approximate the aggregate signal received from the hybrid network to address the difficulty of distribution analysis due to conjugate beamforming. The coverage probability is then obtained by the Laplace transform for interference. The analysis of the coverage probability of hybrid networks is validated by MC simulation, demonstrating that hybrid networks can reduce the performance gap while improving the peak SINR performance.
VII Acknowledgment
This work was supported by the National Key R&D Program of China with Grant number 2019YFB1803400, and by the National Natural Science Foundation of China with grant number 62071114.
References
- [1] H. Q. Ngo, A. Ashikhmin, H. Yang, E. G. Larsson, and T. L. Marzetta, “Cell-free massive MIMO versus small cells,” IEEE Trans. Wirel. Commun., vol. 16, no. 3, pp. 1834–1850, Mar. 2017.
- [2] D. Wang, M. Wang, P. Zhu, J. Li, J. Wang, and X. You, “Performance of network-assisted full-duplex for cell-free massive MIMO,” IEEE Trans. Commun., vol. 68, no. 3, pp. 1464–1478, Mar. 2020.
- [3] E. Nayebi, A. Ashikhmin, T. L. Marzetta, H. Yang, and B. D. Rao, “Precoding and power optimization in cell-free massive MIMO systems,” IEEE Trans. Wirel. Commun., vol. 16, no. 7, pp. 4445–4459, Jul. 2017.
- [4] Z. Dai, J. Xu, Y. Zeng, S. Jin, and T. Jiang, “Characterizing the rate region of active and passive communications with RIS-based cell-free symbiotic radio,” IEEE Internet Things J., pp. 1–1, 2023, early access, doi:10.1109/JIOT.2023.3308970.
- [5] T. Kim, H. Kim, S. Choi, and D. Hong, “How will cell-free systems be deployed?” IEEE Commun. Mag., vol. 60, no. 4, pp. 46–51, Apr. 2022.
- [6] S. Elhoushy and W. Hamouda, “Towards high data rates in dynamic environments using hybrid cell-free massive MIMO/small-cell system,” IEEE Wireless Commun. Lett., vol. 10, no. 2, pp. 201–205, Feb. 2021.
- [7] A. Papazafeiropoulos, P. Kourtessis, M. D. Renzo, S. Chatzinotas, and J. M. Senior, “Performance analysis of cell-free massive MIMO systems: A stochastic geometry approach,” IEEE Trans. Veh. Technol., vol. 69, no. 4, pp. 3523–3537, Apr. 2020.
- [8] T. M. Hoang, H. Q. Ngo, T. Q. Duong, H. D. Tuan, and A. Marshall, “Cell-free massive MIMO networks: Optimal power control against active eavesdropping,” IEEE Trans. Commun., vol. 66, no. 10, pp. 4724–4737, Oct. 2018.
- [9] Z. Chen and E. Björnson, “Channel hardening and favorable propagation in cell-free massive MIMO with stochastic geometry,” IEEE Trans. Commun., vol. 66, no. 11, pp. 5205–5219, Nov. 2018.
- [10] Y. J. Chun, M. O. Hasna, and A. Ghrayeb, “Modeling heterogeneous cellular networks interference using poisson cluster processes,” IEEE J. Sel. Areas Commun., vol. 33, no. 10, pp. 2182–2195, Oct. 2015.
- [11] J. G. Andrews, F. Baccelli, and R. K. Ganti, “A tractable approach to coverage and rate in cellular networks,” IEEE Trans. Commun., vol. 59, no. 11, pp. 3122–3134, Nov. 2011.
- [12] “Isotropic vectors,” http://encyclopediaofmath.org/index.php?title=Isotropic_vector&oldid=47446, Encyclopedia of Mathematics.
- [13] R. W. Heath Jr, T. Wu, Y. H. Kwon, and A. C. K. Soong, “Multiuser MIMO in distributed antenna systems with out-of-cell interference,” IEEE Trans. Signal Process., vol. 59, no. 10, pp. 4885–4899, Oct. 2011.
- [14] P. G. Moschopoulos, “The distribution of the sum of independent gamma random variables,” Ann. Inst. Stat. Math., vol. 37, no. 1, pp. 541–544, 1985.
- [15] L.-F. Huang, “The nakagami and its related distributions,” WSEAS Trans. Math, vol. 15, no. 44, pp. 477–485, 2016.
- [16] J. Lyu and R. Zhang, “Hybrid active/passive wireless network aided by intelligent reflecting surface: System modeling and performance analysis,” IEEE Trans. Wirel. Commun., vol. 20, no. 11, pp. 7196–7212, Nov. 2021.
- [17] H. Pishro-Nik, Introduction to probability, statistics, and random processes. Kappa Research, LLC Blue Bell, PA, USA, 2014.
- [18] R. J. Muirhead, Aspects of Multivariate Statistical Theory. Wiley Series in Probability and Statistics. New York, NY, USA: Wiley, 1982., 1982.
- [19] W. P. Johnson, “The curious history of Faà di Bruno’s formula,” Amer. Math. Monthly, vol. 109, no. 3, pp. 217–234, 2002. [Online]. Available: http://www.maa.org/news/monthly217-234.pdf
- [20] R. Tanbourgi, H. S. Dhillon, J. G. Andrews, and F. K. Jondral, “Dual-branch MRC receivers under spatial interference correlation and nakagami fading,” IEEE Trans. Commun., vol. 62, no. 6, pp. 1830–1844, Jun. 2014.
- [21] V. F. Ivanoff, “Problems for solution: 4782,” Amer. Math. Monthly, vol. 65, no. 3, p. 212, 1958. [Online]. Available: http://www.jstor.org/stable/2310076.