Next Article in Journal
Cost-Effective Open-Ended Coaxial Technique for Liquid Food Characterization by Using the Reflection Method for Industrial Applications
Next Article in Special Issue
Resource Allocation in Multi-Carrier Multiplexed NOMA Cooperative System
Previous Article in Journal
Intelligent Reflecting Surfaces Assisted UAV Communications for Massive Networks: Current Trends, Challenges, and Research Directions
Previous Article in Special Issue
Enhancing PAPR and Throughput for DFT-s-OFDM System Using FTN and IOTA Filtering
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Guaranteeing QoS for NOMA-Enabled URLLC Based on κμ Shadowed Fading Model

1
School of Communication and Information Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China
2
Beijing National Research Center for Information Science and Technology, Tsinghua University, Beijing 100084, China
3
School of Information and Communication Engineering, Beijing University of Posts and Telecommunications, Beijing 100876, China
4
Department of Electronic Engineering, Tsinghua University, Beijing 100084, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(14), 5279; https://doi.org/10.3390/s22145279
Submission received: 22 May 2022 / Revised: 30 June 2022 / Accepted: 3 July 2022 / Published: 14 July 2022

Abstract

:
Sixth-generation (6G) wireless communication scenarios are complex and diverse. Small-scale fading is a key part of wireless channels and its impact on performance in scenarios with time sensitivity and 6G ultrareliable and low latency communications (URLLC) quality-of-service requirements cannot be ignored. Therefore, it is necessary to accurately characterize small-scale fading when designing wireless communication systems. In this paper, we derive approximate closed form expressions for the probability density function, cumulative distribution function and moment-generating function of the postprocessing signal-to-noise ratio following the zero-forcing detector in a cell-free massive multiple-input multiple-output (CF mMIMO) system. CF mMIMO system is a nonorthogonal multiple access (NOMA) system that enables users to share all channel uses and can ensure the fairness of the communication quality experienced by different users. Our key contributions include the extension of the κ μ shadowed fading model to a CF mMIMO system and the proposal of theoretical tools (the derived closed-form expression) to improve its mathematical tractability. By exploiting the statistical characterizations of the arrival and service processes, another important contribution is the exploitation of the upper bound of the queuing delay violation probability (UB-QDVP) over the Mellin transforms of the arrival and service processes in the proposed CF mMIMO system under the κ μ shadowed fading model. Corroborated by extensive simulations, our analyses validate that the CF mMIMO system outperforms the orthogonal multiple access and power-domain NOMA systems and reveal the relationships among different small-scale fading types, energy efficiency, delay and the UB-QDVP, as well as the accuracy and effectiveness of the proposed theoretical tools based on the κ μ shadowed fading model.

1. Introduction

1.1. Background and Motivation

Due to the highly time-varying nature of wireless channels, it is difficult to design communication systems that support ultrareliable and low latency communications (URLLC) with reductions in energy consumption. For example, the explosively increasing volumes of delay-sensitive and bandwidth-intensive applications under stringent quality-of-service (QoS) requirements have required strict delay (<0.1 ms), ultrareliability (>99.99999%), and extra-high energy efficiency (EE) demands [1]. Notably, achieving higher data rates and lowering power consumption seem to be contradictory objectives. In particular, when the delay boundary approaches the channel coherence time, the transmit power required to obtain the QoS requirements becomes unbounded [2]. Although both academia and industry have paid attention to the EE of cellular networks in the past few years [3], the related works are still unable to address the increasing complexity of devices, antennas and many frequency bands in the future [4].
Orthogonal multiple access (OMA) systems (e.g., multiple-input multiple-output orthogonal frequency division multiple access (MIMO-OFDMA)) apply the orthogonal resource allocation method, in which different uplink users occupy different channel uses (CUs). This type of system distinguishes different users through different orthogonal time–frequency resource blocks. However, all users share time–frequency resources in the non-OMA (NOMA) system. As a promising technology, NOMA makes full use of the nonorthogonal nature over limited orthogonal resources so as to provide a higher transmission accuracy with lower computational complexity [5], and to meet the massive connectivity requirements in an Internet of things scenario [6,7] and even satellite-based communications [8]. For instance, power-domain NOMA (PD-NOMA), considered as the representative of single carrier NOMA technology [9], distinguishes different groups of users through the spatial domain, and different users in the same group share spatial domain resources. Under PD-NOMA systems, power control is required to meet the QoS requirements of URLLC for energy efficiency, latency and reliability [10]. However, fairness among the users in terms of their throughput is an important goal of NOMA power allocation [11]. For the upcoming sixth generation (6G) mobile networks, a beneficial embodiment of a distributed massive MIMO system is the cell-free massive MIMO (CF mMIMO) system proposed in [12]. As a type of NOMA system, the CF mMIMO system is a better solution between PD-NOMA and MIMO-OFDMA. All users share all CUs, and different users are distinguished by the spatial domain without considering the fairness experienced by uplink users. The CF mMIMO system contains a large number of randomly distributed access points (APs) to coherently serve massive mobile devices; its implementation is more feasible for practical issues than network MIMO with low overhead that exploits channel hardening phenomena and favorable propagation properties [4,13].Relative to traditional cellular mMIMO networks, the notable characteristics of CF mMIMO system are its superior macro diversity and interference suppression capabilities. Macro diversity offers orders of magnitude of coverage probability [12,14], reliability and delay optimization through suitable signal processing. Intercell interference is effectively mitigated by the user-centric approach and coordination [15,16], thereby improving the accuracy of multiuser detection. To the best of our knowledge, a CF mMIMO system that efficiently integrates the unified fading model to provide substantial performance improvements with respect to different types of small-scale fading has still not been well studied.
On the other hand, most current papers concerning the performance analysis in CF mMIMO systems [4,17,18,19] perform analyses in the context of small-scale fading models such as the Rayleigh and Nakagami-m models, which ignore the random nature of wireless fading channels. Such classic fading models usually do not cover all types of fading scenarios encountered in practical cases. More importantly, the analyses performed for practical channel estimations depend on the channel’s statistical properties, which are affected by the selected channel models [20]. Consequently, many results are not reasonably universal enough due to simplified assumptions, such as in many studies of CF mMIMO systems performance based on the Rayleigh model [12,14,21] or Rician model [22]. The design of reliable wireless communication systems is highly dependent on the accurate characterization of the wireless fading channel. However, for the complex communication scenarios in future 6G URLLC scenarios [23,24], the classic fading model is no longer applicable [25,26]. With the research on fading models in recent years, a new fading model that can losslessly portray a variety of fading types has the potential to become a unified fading model [27,28,29,30], which can be used to accurately characterize fading channels in complex future 6G URLLC scenarios. This is a critical issue regarding the theoretical analysis and performance optimization of URLLC, given that the reliability of modeling is of capital importance. Specifically, the κ μ shadowed fading model, which unifies the κ μ fading distribution and the η μ fading distribution, was introduced in [27]. First, by adjusting the parameters κ , μ and m, almost all classical fading models can be unified by the κ μ shadowed fading model [28] with similar mathematical flexibility levels. Second, this shadowed model provides much better tractability and flexibility [31]. These two features allow the κ μ shadowed fading model to accommodate different propagation conditions, such as in practical communication environments that experience the effects of nonhomogeneous fading conditions. However, although the κ μ shadowed fading model can characterize a variety of fading channels, adapting the κ μ shadowed fading model to the CF mMIMO system and specifically analyzing the system performance while supporting URLLC becomes a challenging research issue.
With respect to mathematical tools, towards the exponentially growing demands of time-sensitive services, the traditional Shannon information theory based on the infinite blocklength assumption and a random coding scheme provides a poor benchmark in this scenario [32]. Finite-blocklength (FBL) information theory was proposed in [33] as a short-packet data communication technique for guaranteeing the stringent QoS requirements of URLLC [34]. Well-known results regarding the maximum achievable data rate of FBL packets were derived in [33] as a function of the blocklength and error probability. Moreover, as suggested by Bennis et al. [35] and Arnau et al. [36], the design of URLLC should capture the tail behavior of reliability (i.e., queuing delay violation probability) and delay instead of average metrics. Because the delay and overall reliability consist of multiple components, the corresponding performance indicators should be constrained by a delay boundary and a delay violation probability boundary for URLLC [36]. In this sense, the stochastic network calculus (SNC) approach proposed in [37] provides a nonasymptotic boundary on the delay violation probability as an effective analysis tool for connecting the delay and overall reliability.
Therefore, the CF mMIMO system under the κ μ shadowed fading model is sufficient for analyzing performance in practical communication environments to meet the strict QoS requirement of URLLC. However, a major problem with this kind of performance analysis is appropriately characterizing stochastic queuing behaviors when they are integrated with the CF mMIMO system under the κ μ shadowed fading model with SNC and FBL information theory. Furthermore, it is still an urgent issue to exploit the effect of small-scale fading characteristics on the performance of the CF mMIMO system under the κ μ shadowed fading model to achieve NOMA-enabled URLLC with high resource efficiency.

1.2. Related Works

Our paper builds on results derived from several research areas. On the physical layer, we consider the uplink CF mMIMO system with imperfect channel state information (CSI) under the κ μ shadowed fading model as well as FBL information theory. On top of a service description in the physical-layer fading channel model, we then investigate service processes. We determine that the effect of different types of small-scale fading on the performance in the design, modeling and analysis of several wireless communication systems cannot be neglected.
The CF mMIMO system has been gradually considered in the URLLC scenario. Nasir et al. [38] first investigated a CF mMIMO system for downlink URLLC in the FBL regime and designed a special class of conjugate beamforming to achieve a lower computational complexity and improved URLLC rate performance. Elwekeil et al. [19] introduced power control schemes in a CF mMIMO system that supported URLLC applications for both traditional ground users and unmanned aerial vehicles, and they evaluated the superior rate performance of this CF mMIMO system. Due to the spatial domain, the CF mMIMO method can be a key to boosting the reliability of URLLC applications. The study in [39] compared different transmission modes in a factory automation scenario and demonstrated that the user-centric transmission mode and power control brought substantial performance improvements in terms of reliability and latency. Zhang et al. [40] integrated the CF mMIMO approach with simultaneous wireless information and power transfer for statistical delay and error-rate bounded QoS to support URLLC and boost the data rate and energy efficiency. Hence, a CF mMIMO network can be viewed as a promising system to support URLLC with a high resource efficiency and low latency.
The κ μ shadowed fading model [27,28] was proposed as a more general yet equally tractable model to capture a wide range of propagation conditions with clear physical interpretation and good analytical properties. Furthermore, Lopez-Martinez et al. [31] derived the probability density function (PDF) and cumulative distribution function (CDF) of the κ μ shadowed fading model with integer fading parameters. Ramirez-Espinosa et al. [25] provided a statistical characterization of the κ μ shadowed fading model and asserted that these results could be used to derive some performance metrics of wireless communication systems over fading channels. With respect to the characterization of the κ μ shadowed fading model, its applications have been studied in several works. ElHalawany et al. [41] carried out a performance (i.e., ergodic capacity, outage probability and average bit error rate) analysis of downlink NOMA systems subject to the κ μ shadowed fading model. Chun et al. [42] proposed an analytic framework to evaluate the average of an arbitrary signal-to-noise-plus-interference ratio (SINR) function over the κ μ shadowed fading model and evaluate the spectral efficiency, moments of the SINR and outage probability of heterogeneous cellular systems. The authors in [43] characterized the impacts of realistic propagation conditions on the achievable secrecy performance of MIMO systems over the κ μ shadowed fading model. Fully absorbing the attractive nature of the κ μ shadowed fading model is seen as one of the potential keys to characterizing the propagation media in emerging practical scenarios.
Several works have also examined SNC-based network performance analysis and theoretical boundary calculation for URLLC. Schiessl et al. used SNC to investigate the delay performance of a multiuser MIMO system with zero-forcing beamforming [44] and NOMA with joint decoding [45] under imperfect CSI and FBL. Xiao et al. [46] derived a closed-form of the upper bound of the queuing delay violation probability (UB-QDVP) in a downlink NOMA system and designed an optimal power allocation scheme for guaranteeing the delay violation probability. Furthermore, by utilizing moment generating function (MGF)-based SNC, probabilistic delay bounds for traffic dispersion, network densification and a hybrid scheme were obtained in millimeter-wave communications over Nakagami-m fading [47]. A study close to our work is the research by Zhang et al. [22], where the stochastic QoS performance in terms of both the delay and error rate metrics was obtained for a downlink CF mMIMO based system model across Rician fading in the FBL regime. However, only traditional fading channels have been considered for discussion in the related literature, while no investigation is available with respect to the performance analysis for CF mMIMO systems with the juxtaposition of shadowing and small-scale fading.

1.3. Main Contributions

For 6G URLLC requirements, we combine the potential key technology of future 6G, CF mMIMO [48], to build the uplink wireless communication system, and extend the κ μ shadowed fading model to the CF mMIMO system, for lossless characterization of its fading channel in the complex communication scenario of 6G [23,24]. The key contributions of the paper can be summarized as follows.
  • We derive approximate closed-form expressions for statistical characteristics (the PDF, CDF and MGF) of the sum of independent and nonidentically distributed (i.n.i.d.) κ μ shadowed random variables (RVs). The analysis is nontrivial, as we optimize the theoretical analysis performance of the κ μ shadowed fading model and improve its mathematical tractability.
  • Based on the κ μ shadowed fading model, we derive approximate closed-form expressions for the PDF, CDF and MGF of the postprocessing signal-to-noise ratio (PPSNR) after the zero-forcing detector in the proposed CF mMIMO system, and extend the κ μ shadowed fading model to 6G NOMA wireless communication systems.
  • By utilizing FBL information theory, SNC and the Mellin transform on the service process, we exploit the UB-QDVP in the proposed CF mMIMO system under the κ μ shadowed fading model. Furthermore, based on extensive simulations, we analyze the system performance with the delay and UB-QDVP indicators and validate the necessity of analyses performed under the κ μ shadowed fading model; the CF mMIMO system outperforms the OMA system and PD-NOMA system. (For the OMA and PD-NOMA systems, there are L antennas equipped in their base station and K users. For illustration convenience, in the PD-NOMA system, we divide K users into N pairs based on the channel gains of the users ( K = 2 N ). Specifically, we divide K users into a group of N “strong users” and a group of N “weak users” according to their channel gains and sort the users in the descending order of their channel gains within each of the groups. The “strong user” and “weak user” with the same label are grouped in the same pair.)

1.4. Organization

The remainder of this paper is organized as follows. In Section 2, the system model for a CF mMIMO system under the κ μ shadowed fading model is described. Section 3 introduces SNC and presents the Mellin transform over both arrival and service processes as well as the UB-QDVP. Simulation results and an analysis are given in Section 4, followed by the conclusions in Section 5.
Notations: C K × N is a complex matrix with K rows and N columns. The modulus and expectation operators are denoted by and E ( ) , respectively. ( ) T and ( ) H denote the transpose and conjugate transpose, respectively, of their arguments. ( ) l k denotes the element in the lth row and kth column of a matrix. C N ( 0 , 1 ) denotes a complex Gaussian distribution with a mean of zero and unit variance.

2. System Model

We consider an uplink CF mMIMO system in which L single-antenna APs and K single-antenna users are randomly uniformly distributed in a circular area with L K , as shown in Figure 1. Some conventional fading can be obtained by adjusting the parameters of the κ μ shadowed fading model: Rayleigh fading ( κ 0 , μ = 1 , m ), Rician fading ( κ = κ _ , μ = 1, m ), Rician shadowed fading ( κ = κ _ , μ = 1, m = m _ ), and other fading ( κ = κ _ , μ = μ _ , m = m _ , where κ _ , μ _ and m _ are nonnegative real numbers. When m _ , it denotes that no shadow effect is observed). All APs simultaneously serve all users at the same time–frequency resource and are connected to a baseband unit (BBU) for centralized signal processing via ideal backhaul links [19,49] (There is no delay in the ideal backhaul links between APs and the BBU [50]). The transmission from the APs to the users (downlink transmission) and the transmission from the users to the APs (uplink transmission) proceed by a time-division-duplexing operation. Each coherence interval is divided into three phases: uplink training, downlink payload data transmission, and uplink payload data transmission. In the uplink training phase, the users send pilot sequences to the APs and the BBU estimates the channel for all users. In this paper, the QoS delay for URLLC is characterized by the queuing delay (Considering ideal backhaul links and tail behaviors based on [35,36], we consider the queuing delay as the delay of the URLLC QoS requirement). Each user transmits a short packet with D information bits through N CUs, which are spread over B MHz of bandwidth and t f milliseconds of time ( N = B t f ).
The channel coefficient g l k t between the lth AP and the kth user can be characterized as
g l k t = β l k t h l k t ,
where β l k t is a large-scale fading coefficient that changes very slowly with time and is assumed to be constant over many coherence conditions. h l k t is the small-scale fading coefficient of the corresponding links, which is assumed to be static during a finite-sized time–frequency coherence interval, and it changes independently from one coherence interval to the next [51]. The channel response is assumed to follow the κ μ shadowed fading model. The gain of small-scale fading between the kth ( k 1 , K ) user and the lth AP ( l 1 , L ) can be modeled as
h l k 2 = i = 1 μ l k X l k , i + ξ l k p ^ l k , i 2 + Y l k , i + ξ l k q ^ l k , i 2 ,
where μ l k is a natural number. Each multipath cluster is modeled by one term of the sum; thus, μ l k is the number of multipath clusters. Each cluster has scattered components with the same power and dominant components with a certain arbitrary power. X l k , i and Y l k , i are independent and identically distributed (i.i.d.) Gaussian RVs with means of zero and variances of σ 2 , and X l k , i + j Y l k , i denotes the scattered components of the ith cluster. p ^ l k , i and q ^ l k , i are real numbers and the mean values of the in-phase and quadrature parts of the multipath components of the ith cluster, respectively. ξ l k is a Nakagami-m RV with shaping parameter and E ξ l k 2 = 1 , where ξ l k 2 Γ m l k , 1 1 m m l k ; m l k represents the fluctuation degree of dominant components caused by the shadow [28]. For each cluster, the total power of the scattered components is 2 σ 2 . Furthermore, the ratio of dominant components to scattered components in small-scale fading is κ l k = i = 1 μ p ^ l k , i 2 + q ^ l k , i 2 2 μ l k σ 2 . In addition, h l k 2 statistically follows the κ μ shadowed distribution, whose PDF can be expressed approximately by the Gamma distribution, namely, h l k 2 Γ K l k , θ l k , where K l k = m l k μ l k 1 + κ l k 2 m l k + μ l k κ l k 2 + 2 m l k κ l k = μ l k 1 + κ l k 2 1 + μ l k κ l k 2 μ l k κ l k 2 m m l k + 2 κ l k ; θ l k = Ω l k K l k ; Ω l k is the mean value of h l k 2 .
The received signal at the BBU under imperfect CSI can be denoted as in the following equation.
Y p t , Y d t = p u G t X p t , X d t + Z p t , Z d t ,
where p u is the average transmit power of each user. Y p t C L × n and Y d t C L × n ^ are received pilot signal and received data signal transmitted by the users at time slot t, respectively. n and n ^ are the length of the pilots (LoP) and the length of the data signal, respectively, which satisfy n + n ^ = N = B t f . X p t = x 1 p t , , x K p t T C K × n is the transmitted pilot signal at slot t, which satisfies X p t X p t H = n I K . X d t = x 1 d t , , x K d t T C K × n ^ is the transmitted data signal at time slot t, obtained by assuming that the elements of X d t are not related to each other, namely E x k d t 2 = n , E x k d t T x j d t = 0 , k , j 1 , K and k j . G t = g 1 t , g K t C L × K is the CSI of K users and L APs, where g k t = g 1 k t , g L k t T , k 1 , K . As we mentioned above, the expression for the channel coefficient is g l k t = β l k t h l k t , k 1 , K , and l 1 , L , which can be simplified as g l k = β l k h l k , g k t g k and G t G because the large-scale fading coefficient β l k t and small-scale fading coefficient h l k t are generally assumed to be RVs that are independent of time. Z p t and Z d t represent additive white Gaussian noise (AWGN), whose elements are i.i.d.; this can be expressed as C N ( 0 , 1 ) . In a practical communication system, serving users with a large number of APs requires accurate CSI to be available. The estimated channel matrix G ^ obtained via the least squares estimation method which relies on the received signal with the known pilot sequence is expressed as [52]
G ^ = Δ G + E ^ ,
where G ^ = g ^ 1 , g ^ K C L × K , Δ = Λ 1 , , Λ K T C K × L , Λ k = Λ 1 k , , Λ L K T and Λ l k = n p u β l k Ω l k 1 + n p u β l k Ω l k . β l k and Ω l k are the large-scale fading coefficient and mean square value of the small-scale fading coefficient, respectively. E ^ = e ^ 1 , , e ^ K C L × K is the channel estimation error matrix, and its elements e ^ k = e ^ 1 k , , e ^ L k T follow e ^ l k C N 0 , σ e ^ l k 2 , where σ e ^ l k 2 = Λ l k 2 σ e l k 2 = Λ l k 2 n p u = β l k Ω l k 1 + n p u β l k Ω l k . The elements of G ^ can be expressed as g ^ l k = Λ l k g l k + e ^ l k = β l k h ^ l k ; thus, h ^ l k = Λ l k h l k + e ^ l k β l k .
Adopting the standard linear zero-forcing detector A ^ = G ^ ( G ^ H G ^ ) 1 C L × K , the received signal of the kth user is
R ( t ) = p u A ^ H G t X d t + A ^ H Z d t .
After completing zero-forcing detection, the PPSNR of the kth user is [53]
γ ^ k = p u p u l = 1 L a ^ l k 2 n = 1 K σ e ^ ln 2 + a ^ k 2 ( a ) p u K n a ^ k 2 + a ^ k 2 = p u K n + 1 a ^ k 2 ( b ) p u K n + 1 l ^ L k g l ^ k 2 = p ^ u l ^ L k η l ^ k = p ^ u ζ k ,
where L k = L / B k = l ^ l ^ = 1 , . . . , L ^ , L = l l = 1 , . . . , L , B k is defined as B k = Δ Unique l n * = arg max l β l k n k and Unique T returns the same values as those in the set T but with no repetitions. a ^ l k is the lth entry of a ^ k , and a ^ k is the kth column of A ^ . η l ^ k is the channel gain of the lth AP and the kth user and is an RV with both large-scale fading coefficient and a small-scale fading coefficient. (a) is based on the consideration that when p u is high enough, p u σ e ^ lk 2 1 n . (b) occurs because 1 a ^ k 2 l ^ L k g l ^ k 2 [53]. Because the gain of small-scale fading h l k 2 statistically follows the κ μ shadowed distribution, the distribution of η l ^ k is obtained in Lemma 1 as follows.
Lemma 1.
Based on a κ–μ shadowed fading model, the channel gain η l ^ k follows the Gamma distribution, which is
η l ^ k Γ K l ^ k , Q l ^ k ,
where K l ^ k = m l ^ k μ l k 1 + κ l ^ k 2 m l l ^ k + μ l ^ k κ l ^ k 2 + 2 m l ^ k κ l ^ k = μ l ^ k 1 + κ l ^ k 2 1 + μ l ^ k κ l ^ k 2 μ l ^ k κ l ^ k 2 m m l ^ k + 2 κ l ^ k , and Q l ^ k = β l ^ k θ l ^ k = β l ^ k Ω l ^ k K l ^ k = Ω ^ l ^ k K l ^ k . β l ^ k Ω l ^ k is the mean value of η l ^ k , as well as the mean value of the channel gain.
Proof. 
Please refer to Appendix A. □
Based on the κ μ shadowed fading model, ζ k = l ^ L k η l ^ k can be regarded to the sum of i.n.i.d. Gamma RVs. Then, the PDF, CDF and MGF of ζ k are obtained from Theorem 1.
Theorem 1.
Considering a κ–μ distribution with parameters ( A k , K k , Ω ^ k ), the PDF can be approximately obtained by a V-order generalized Puiseux series expansion process as follows.
f ζ k z ; A k , K k , Ω ^ k v = 1 V c v A k D k A k K k ω v Ω k A k K k Γ K k z A k K k 1 exp D k z ω v Ω ^ k A k ,
where A k = A l ^ k = 1 , K k = l ^ L k K l ^ k , K k = l ^ L k K l ^ k and D k = Γ K k + 1 1 A k A k Γ K k .
We can obtain approximate parameters c v and ω v via the following formula:
v = 1 V c v = 1 , v = 1 V c v ω v = 1 , v = 1 V c v ω v 2 Ω k , 2 2 Γ K k Γ K k + 1 1 A k A k 2 Γ K k + 2 2 A k A k Γ K k = E ζ k 2 , v = 1 V c v ω v 2 V 2 Ω k , 2 2 V 2 Γ K k Γ K k + 1 1 A k A k 2 V 2 Γ K k + 2 M k 2 2 M k 2 A k A k Γ K k = E ζ k 2 V 2 , v = 1 V c v ω v A k K k A k D k A k K k Ω k , 2 A k K k Γ K k = A L k ,
where A L k = l ^ L k A D l ^ k A K l ^ k Ω l ^ k , 2 A K l ^ k Γ K l ^ k l ^ L k Γ A K l ^ k Γ A l ^ L k K l ^ k and
E ζ k n = n 1 = 0 n n 2 = 0 n 1 n L ^ 1 = 0 n L ^ 2 n n 1 n 1 n 2 n L ^ 2 n L ^ 1 E X 1 n n 1 E X 2 n 1 n 2 E X L ^ n L 1 .
The CDF of ζ k is given by
F ζ k z ; A k , K k , Ω ^ k v = 1 V c v γ K k , D k z ω v Ω ^ k A k Γ K k ,
where γ · , · is the lower incomplete Gamma function and is expressed by
γ s , x = 0 x t s 1 exp t d t .
The MGF of ζ k is
M ζ k s = A k D k A k K k Ω ^ k A k K k Γ K k v = 1 V c v ω v A k K k q 1 / 2 p A k K k 1 / 2 2 π p + q p + q 2 2 1 s A k K k G p , q q , p D k ω v Ω ^ k A k q p p s p q q Δ p , 1 A k K k Δ q , 0
where G p , q q , p · is the Meijer G-function and p q = A k = A n , k = 1 . According to gcd p , q = 1 , both p and q can be 1.
Then, the MGF of ζ k can be simplified to
M ζ k s = A k D k A k K k Ω ^ k A k K k Γ K k v = 1 V c v ω v A k K k 1 s A k K k G 1 , 1 1 , 1 D k ω v Ω ^ k A k 1 s Δ 1 , 1 A k K k Δ 1 , 0 .
Proof. 
Please refer to Appendix B. □
Based on Equation (6) and Theorem 1, the statistical characteristics of the PPSNR in the CF mMIMO system are provided in the following corollary.
Corollary 1.
Due to γ ^ k p ^ u ζ k , the PDF of the PPSNR γ ^ k at the BBU is derived as
f γ ^ k z ; A k , K k , p ^ u Ω ^ k v = 1 V c v A k D k A k K k ω v p ^ u Ω ^ k A k K k Γ K k z A k K k 1 exp D k z ω v p ^ u Ω ^ k A k .
The CDF of γ ^ k is given by
F γ ^ k z ; A k , K k , p ^ u Ω ^ k v = 1 V c v γ K k , D k z ω v p ^ u Ω ^ k A k Γ K k .
The corresponding MGF of γ ^ k can be expressed as
M ζ k s = A k D k A k K k Ω ^ k A k K k Γ K k v = 1 V c v ω v A k K k 1 s A k K k G 1 , 1 1 , 1 D k ω v Ω ^ k A k 1 s Δ 1 , 1 A k K k Δ 1 , 0 .

3. Analysis of the Delay Performance

Because each packet is very short, both imperfect CSI and FBL channel coding cause transmission errors. The maximum achievable data rate is subject to transmission errors and an instantaneous signal-to-noise ratio (SNR), as described in Section 3.1. Due to transmission errors and dynamic data rates in the physical layer, the APs need to keep their packets in a buffer until successful completion, which leads to a random queuing delay, as discussed in Section 3.2. Finally, in Section 3.3, we derive the Mellin transform over the arrival and service processes and investigate the UB-QDVP.

3.1. Service Transmission in FBL Regime

Based on FBL information theory, the maximum achievable data rate (bit/CU) with an instantaneous SNR and a decoding error probability for the kth user can be obtained as follows. [22,33]
r k γ ^ k , ε k d = C γ ^ k V γ ^ k n ^ Q 1 ε k d = log 2 1 + γ ^ k log 2 e 1 1 1 + γ ^ k 2 n ^ Q 1 ε k d = log 2 f r k γ ^ k , ε k d ,
where f r k γ ^ k , ε k d = 1 + γ ^ k exp γ ^ k 2 + 2 γ ^ k 1 + γ ^ k 2 n ^ Q 1 ε k d , n ^ is the number of CUs occupied by data transmission and ε d is the decoding error probability. Q x = 1 2 π x e y 2 y 2 2 2 d y is the tail distribution function of the standard normal distribution. C γ ^ k = log 2 1 + γ ^ k and V γ ^ k = 1 1 1 + γ ^ k 2 log 2 e 2 are the Shannon capacity and channel dispersion, respectively. Given a decoding error probability of ε k d , the achievable data rate r k γ ^ k , ε k d is less than zero when the PPSNR γ ^ k is below a threshold γ 0 . As a consequence, we redefine f r k γ ^ k , ε k d as g r k γ ^ k , ε k d
g r k γ ^ k , ε k d = f r k γ ^ k , ε k d , γ ^ k > γ 0 1 , γ ^ k γ 0 .
Furthermore, the maximum achievable data rate can be simplified as follows:
r k γ ^ k , ε k d = log 2 f r k γ ^ k , ε k d , γ ^ k > γ 0 0 , γ ^ k γ 0 .

3.2. Queuing Model

For our queuing system analysis, we consider a usual statistical model such as that in [37], which is widely used in the SNC methodology. Then, the cumulative arrival, achievable service, and departure processes of user k during time interval [ τ , t ) are defined as follows.
A k τ , t = i = τ t 1 a k i , R k τ , t = i = τ t 1 r k i , U k τ , t = i = τ t 1 u k i .
The instantaneous traffic arrival a k i models the number of bits that arrive at the queue within a discrete time slot i. The increment of service process (service rate) r k i is independent of time and equal to the maximum achievable data rate r k γ ^ k , ε k d . The departure process u k i denotes the number of bits that arrive successfully at the destination and is subject to both the service process and the amount of data waiting in the queue. In addition, we assume that all queues are work-conserving first-come-first-served queues. Then, the queuing delay w k t of user k at time slot t is defined as the number of time slots required for all data that arrive before time t to depart from the transmit buffer and is expressed as
w k t = inf u 0 : A k ( 0 , t ) U k ( 0 , t + u ) .

3.3. SNC in the SNR Domain

A notable aspect of SNC is that it can obtain the delay distribution and the UB-QDVP based on statistical characterizations of the arrival and service processes in terms of their Mellin transforms. However, it is still difficult to capture the statistics of the random arrival and service processes in practical scenarios. To facilitate the analysis, we convert these processes from the bit domain through the exponential domain to the SNR domain. These cumulative processes in the SNR domain are denoted by calligraphic letters [54] as follows:
A k τ , t = e A k τ , t , R k τ , t = e R k τ , t , U k τ , t = e U k τ , t .
M X ( s ) = E X s 1 is the Mellin transform of any nonnegative random variable X for a parameter s R [37]. Both a k ( i ) and r k ( i ) are independent of each other at different time slots, so we use two time-independent random variables a k and r k to denote the number of arrival and service bits at either time slot, respectively. Then, the Mellin transform over the accumulated arrival process is denoted by M A k ( s , τ , t ) = M A k τ , t s = E e A k τ , t s 1 = M φ k s t τ , where φ k = e a k . The Mellin transform over the accumulated achievable service process is M R k ( s , τ , t ) = M R k τ , t s = E e R k τ , t s 1 = M ϕ k s t τ , where ϕ k = e r k .
If the stability condition M φ k 1 + s M ϕ k 1 s < 1 holds, denoting t T as a target delay (target number of time slots), the kernel function K s , t T is defined as [37]
K s , t T = M ϕ k 1 s t T 1 M φ k 1 + s M ϕ k 1 s .
Thus, the UB-QDVP denoted by ε k v , UB is [22,37]
ε k v , UB = Pr w k t > t T inf 0 < s < s 0 K s , t T ,
where s 0 = sup s : M φ k 1 + s M ϕ k 1 s < 1 .

3.3.1. Mellin Transform over the Arrival Process in the SNR Domain

As shown in Equation (24), the UB-QDVP is related to the arrival process on the Mellin transform in the SNR domain. Consequently, deriving the Mellin transform of a k is essential for evaluating the UB-QDVP. We consider the traffic class of δ s , λ ( s ) -bounded arrivals, whose arrival processes characterized in the SNR domain for some s > 0 on the Mellin transform are bounded by [37]
M A k ( s , τ , t ) = M φ k s t τ e ( s 1 ) · ( λ ( s 1 ) · ( t τ ) + δ ( s 1 ) ) .
To simplify the notation, we consider the case where δ ( s ) = 0 and λ ( s ) λ ; namely, δ s , λ ( s ) are independent of s, which is true for constant arrivals [54]. Thus,
M φ k 1 + s e λ s ,
where λ is the constant number of arrival bits that arrive successfully at one slot, as well as the constant bit arrival rate.

3.3.2. Mellin Transform over the Service Process in the SNR Domain

Equation (24) also shows that deriving the closed-form expression for the Mellin transform over the service process in the SNR domain at user k is a key component when analyzing the UB-QDVP, which motivates the following Theorem 2.
Theorem 2.
The Mellin transform over the service process in the SNR domain of user k at one time slot is
M ϕ k 1 s I k s , ε k d + G k γ 0 , ε k d ,
where I k s , ε k d = 1 ε k d M γ ^ k s ln 2 and G k γ 0 , ε k d = ε k d + 1 ε k d F γ ^ k γ 0 . ε k d is the decoding error probability. γ 0 is the upper bound of the PPSNR when the maximum achievable data rate of user k is 0.
Proof. 
Please refer to Appendix C. □
Combining and extending Theorem 2, Equations (24) and (26) yield the following corollary.
Corollary 2.
The approximate UB-QDVP of the kth user is
ε k v , UB I k s , ε k d + G k γ 0 , ε k d t T 1 e λ s I k s , ε k d + G k γ 0 , ε k d .
Proof. 
Please refer to Appendix D. □

4. Simulation Results and Analyses

We present numerical results to validate and evaluate the stochastic queuing behaviors and the effects of fading characteristics on performance in the proposed CF mMIMO system under the κ μ shadowed fading model. Different from most studies [4,12,13,17,18], which explore the effects of the number of APs, the coverage radio, the number of potential users and other traditional system parameter configurations on performance in communication scenarios, we focus on fading characteristics. Throughout our simulations, we set such parameters constant. For large-scale fading β l k = u p d k d 0 α p , u p , α p , d 0 and d l k are the constant path loss across the minimum distance, the path loss exponent, the minimum distance and the distance between the kth user and the lth AP, respectively. The key simulation parameters are listed in Table 1. The overall reliability for the kth user can be represented by the overall error probability ε k [55], which is ε k ε k d + ε k v . Because ε k d 0 , ε k can be expressed by ε k v , UB .
First, we compared the maximum achievable data rates obtained under the same decoding error probability but different fading models in Figure 2, when the transmit power is 5 dBm and the distance from the user to the center of the coverage area is 200 m. It is obvious that the maximum achievable data rates corresponding to the same decoding error probability are inconsistent under different fading environments. This indicates that the effect of different fading environments cannot be ignored with respect to system performance. In particular, when designing a wireless communication system (e.g., the CF mMIMO system) for strict decoding error probability, it is crucial to design customized coding design and power control schemes by combining the fading characteristics of the wireless communication environment to adjust the service rate to meet the reliability requirements of NOMA-enabled URLLC QoS. Therefore, the κ μ shadowed fading model accommodating different propagation conditions that are applied to design or analyze the performance of a wireless communication system is necessary and significant.
Figure 3 shows the fitting ability of our derived theoretical tool (Theorem 1) for the PDFs of the sum of i.n.i.d. RVs. We used the Monte Carlo method to randomly generate six independent RVs with different distributions and sum some of them. We verify that the cylindrical diagram (obtained by a simulation using the Monte Carlo method for 10 7 iterations) and the curve with an asterisk (obtained by numerical calculation (Equation (8)) have almost the same slope. This indicates that our proposed theoretical tools based on Theorem 1 can accurately characterize the statistical characteristics of the sum of multiple i.n.i.d. RVs. Furthermore, we summarize the second-order approximate parameters obtained under the different numbers of summed elements and different distributions listed in Table 2, which were obtained by using fsolve in MATLAB. From Figure 3, we also observe that the generalized Puiseux series expansion process converges very quickly, so it can effectively fit the PDF without complex approximate orders. All these properties lay the foundation that appropriately characterizing the statistical characteristics of the sum of multiple i.n.i.d. shadowed RVs is the key to unsolved problems. This again verifies the effectiveness and good analytical properties of our proposed theoretical tool (Theorem 1) in complex future communication environments.
For the proposed CF mMIMO system, Figure 4 depicts the PDFs of the PPSNR under different fading environments when the transmit power is 10 dBm and the distance from the user to the center of the coverage area is 200 m. We verify that the cylindrical diagram (obtained by simulation using Monte Carlo method for 10 7 iterations) and the curve with asterisk (obtained by numerical calculation (Corollary 1)) have almost the same slope. This indicates that our proposed theoretical tools based on Corollary 1 have superior flexibility, mathematical tractability and practicality for wireless communication system design. Furthermore, these findings also reveal the κ μ shadowed fading model based on our proposed theoretical tools, which form a reliable and mathematically tractable channel fading model can be extended to other common future 6G wireless communication systems such as CF mMIMO systems, enabling this model as a more applicable approach to complex 6G applications in the future.
Figure 5 plots the UB-QDVPs obtained with different target delays over a κ μ shadowed fading model for the CF mMIMO system when the transmit power is 10 dBm and the distance from the user to the center of the coverage area is 200 m. The UB-QDVP decreases as the target delay increases. For the same target delay, the users have inconsistent UB-QDVPs under different fading environments. In Figure 5c, under Nakagami-m fading ( κ 0 , μ = 3 , m ) and a target delay of 0.05 ms, the UB-QDVP is 10 7 . However, when the UB-QDVP is 10 7 under One-sided Gaussian fading ( κ 0 , μ = 0.5 , m ), the target delay exceeds 0.05 ms. For the performance evaluation, analysis and design of the system, it is necessary to consider the impacts of different fading environments on the system performance. Otherwise, the final result or conclusion obtained will differ greatly from reality, eventually leading to analysis results with less practical value and significance.
Figure 6a compares the UB-QDVPs obtained under different distances from the user to the center of coverage area (200 m or 1000 m) in three communication systems under One-sided Gaussian fading ( κ 0 , μ = 0.5 , m ) when the transmit power is 10 dBm. These three communication systems are the CF mMIMO system we proposed, PD-NOMA and OMA systems equipped with 200 antennas in their base station. We observe that the UB-QDVP of the PD-NOMA system is lower than that of the OMA system, and the UB-QDVP of the CF mMIMO system is lower than that of the PD-NOMA system under the same target delay when the distance from the user to the center of the coverage area is 200 m. In detail, when the target delay is 0.06 ms and the distance from the user to the center of the coverage area is 200 m, the UB-QDVPs of the PD-NOMA and OMA systems exceed 10 7 , and only the UB-QDVP of the CF mMIMO system is lower than 10 7 . Under a target delay of 0.07 ms, the UB-QDVPs of the CF mMIMO system and the PD-NOMA system are below 10 7 . However, under the same parameter settings, the UB-QDVP of the OMA system is still higher than 10 7 . Specifically, we extended the distance from the user to the center of the coverage area to 1000 m and set the target delay to 0.08 ms. Only the CF mMIMO system can achieve a UB-QDVP lower than 10 7 , and both the PD-NOMA and the OMA systems cannot. Figure 6b plots the UB-QDVP for a user at 200 m or 1000 m from the center of the coverage area under Nakagami-m fading ( κ 0 , μ = 3 , m ) for the CF mMIMO system, PD-NOMA system, and OMA system when the transmit power is 10 dBm. It is also clear from Figure 6b that when the user is 200 m from the center of the coverage area, the UB-QDVP of the CF mMIMO system is lower than that of the PD-NOMA system, which in turn is lower than that of the OMA system, for the same target delay. For example, when the target delay is 0.06 ms, the UB-QDVP is lower than 10 7 for the CF mMIMO system and PD-NOMA system, and higher than 10 7 for the OMA system when the user is 200 m away from the center of the coverage area. When the user is 1000 m away from the center of the coverage area and the target delay is 0.07 ms, only the UB-QDVP of the CF mMIMO system is lower than 10 7 , and the UB-QDVPs of the PD-NOMA system and OMA system are higher than 10 7 . Figure 6a,b confirm that the CF mMIMO system not only optimizes the reliability and delay performance, but also expands the coverage area for providing URLLC services in comparison with the PD-NOMA and OMA systems.
Figure 7a shows the relationship between E E dB and the UB-QDVP under different fading environments in the proposed CF mMIMO communication system when the transmit power is 10 dBm, the target delay is 0.08 ms and the distance from the user to the center of the coverage area is 200 m. As seen from Figure 7a, as the UB-QDVP becomes tighter, E E dB decreases. Under different fading environments, the E E dB of the user for the same UB-QDVP is different, and a gap is observed. The E E dB of the user under Nakagami-m fading ( κ 0 , μ = 3 , m ) for controlling the transmit power to meet the UB-QDVP requirement ( 10 7 ) is improved by 16.6% in E E dB over that of users under Rician shadowed fading ( κ = 4 , μ = 1 , m = 3 ) with the same UB-QDVP. Therefore, for delay-sensitive applications, performance analyses for indicators such as delay and reliability, power control and system design, the impacts of different small-scale fading characteristics on the performance of communication systems need to be considered. The accurate characterization of small-scale fading using the κ μ shadowed fading model may become a necessity for future wireless communication performance analysis and system design.
Figure 7b–d compare the E E dB and UB-QDVP values obtained under different distances from the user to the center of the coverage area (200 m or 1000 m) in three communications systems under Rician shadowed fading ( κ = 4 , μ = 1 , m = 3 ), Rayleigh fading ( κ 0 , μ = 1 , m ) and κ μ shadowed fading ( κ = 3 , μ = 3 , m = 4 ) when the transmit power is 10 dBm, respectively. These three communication systems are the proposed CF mMIMO system, PD-NOMA and OMA systems equipped with 200 antennas in their base station. In Figure 7b, we can see that when the target delay is 0.1 ms and the distance to the center of the coverage area is 200 m, the user achieves ε k v , UB = 10 7 with E E dB < 40 bit/J in the OMA system. Under the same conditions, the user in the PD-NOMA system can arrive at the same UB-QDVP with E E dB > 40 bit/J, which is an E E dB improvement of 5.1% over that of the OMA system. Furthermore, the user in the CF mMIMO system achieves a 2.4% increment in E E dB compared with that of the PD-NOMA system. When E E dB = 27.65 bit/J, the target delay is 0.1 ms and the distance to the center of the coverage area expands to 1000 m, the UB-QDVPs of the user in both PD-NOMA and OMA systems are greater than 10 7 . However, the UB-QDVP of the user in the CF mMIMO system is less than 10 7 and is 97.5% lower than that of the PD-NOMA system. In Figure 7c, when the target delay is 0.1 ms and the distance to the center of the coverage area is 200 m, the user achieves ε k v , UB = 10 7 with E E dB < 40 bit/J in the OMA system. Under the same conditions, the user in the PD-NOMA system can arrive at the same UB-QDVP with E E dB = 40 bit/J, which is an E E dB improvement of 11.1% over that of the OMA system. Furthermore, the user in the CF mMIMO system achieves a 13.9% increment in E E dB compared with that of the OMA system. When E E dB = 27.65 bit/J, the target delay is 0.1 ms and the distance to the center of the coverage area expands to 1000 m, the UB-QDVPs of the user in both PD-NOMA and OMA systems are greater than 10 7 . However, the UB-QDVP of the user in the CF mMIMO system is less than 10 7 and is 99.4% lower than that of the PD-NOMA system. In Figure 7d, when the target delay is 0.1 ms and the distance to the center of the coverage area is 200 m, the user achieves ε k v , UB = 10 7 with E E dB = 40 bit/J in the OMA system. Under the same conditions, the user in the PD-NOMA system can arrive at the same UB-QDVP with E E dB > 40 bit/J, which is an E E dB improvement of 4.9% over that of the OMA system. Furthermore, the user in the CF mMIMO system achieves a 6.25% increment in E E dB compared with that of the OMA system. When E E dB = 27.65 bit/J, the target delay is 0.1 ms and the distance to the center of the coverage area expands to 1000 m, the UB-QDVPs of the user in both PD-NOMA and OMA systems are greater than 10 7 . However, the UB-QDVP of the user in the CF mMIMO system is less than 10 7 and is 99.6% lower than that of the PD-NOMA system. Hence, the CF mMIMO approach is a prospective technology that can improve EE while ensuring that the target delay and reliability meet the URLLC QoS requirements. At the same time, it can greatly decrease the UB-QDVP and improve user reliability, thereby ensuring ultrareliable communication for all users covered.

5. Conclusions

In this paper, we developed an analytical model to characterize stochastic queuing behaviors and the effects of fading characteristics on performance in the proposed CF mMIMO system under the κ μ shadowed fading model. We derived approximate closed-form expressions for the PDF, CDF and MGF of the sum of i.n.i.d. κ μ shadowed RVs and improved their mathematical tractability. Then, based on the κ μ shadowed fading model, we derived approximate closed-form expressions for the PDF, CDF and MGF of the PPSNR after the zero-forcing detector in the proposed CF mMIMO system. The κ μ shadowed fading model can be extended to 6G wireless communication systems, and we verified the reasonableness and feasibility of its superior theoretical analysis performance with simulations and theoretical comparisons. In particular, we applied FBL information theory and SNC to model and characterize the arrival, service and departure processes, and we derived a closed-form expression for the UB-QDVP to analyze the effectiveness and reliability of the CF mMIMO system under the κ μ shadowed fading model. Furthermore, by applying the theoretical analysis above, we also conducted a set of simulations to validate and evaluate the UB-QDVP gaps under different systems (the CF mMIMO, PD-NOMA and OMA systems) and fading environments and the relationships between the UB-QDVPs and EE under different fading environments. Making full use of the characteristics of the κ μ shadowed fading model can improve the system EE, reduce the error probability and better satisfy the QoS requirements of 6G URLLC. Applying the κ μ shadowed fading model to accurately describe small-scale fading and to model and analyze unmanned aerial vehicles, underwater acoustics and other wireless communication systems in the future is a promising direction.

Author Contributions

Conceptualization, J.Z.; methodology, J.Z.; software, Y.S. and T.W.; validation, J.Z., Y.S. and T.W.; data curation, T.W.; writing—original draft preparation, Y.S.; writing—review and editing, J.Z. and Y.S.; supervision, J.Z., T.L. and S.Z.; project administration, J.Z.; funding acquisition, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (no. 62001264), the Natural Science Foundation of Beijing (no. L192025), and the National Key R&D Program (no. 2018YFB1801102).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript.
URLLCUltrareliable and low-latency communications
CF mMIMOCell-free massive multiple-input multiple-output
NOMANonorthogonal multiple access
UB-QDVPUpper bound of queuing delay violation probability
QoSQuality-of-service
EEEnergy efficiency
OMAOrthogonal multiple access
CUChannel use
PD-NOMAPower-domain NOMA
APAccess point
FBLFinite blocklength
SNCStochastic network calculus
CSIChannel state information
PDFProbability density function
CDFCumulative distribution function
SINRSignal-to-noise-plus-interference ratio
MGFMoment-generating function
i.n.i.d.Independent and nonidentically distributed
BBUBaseband unit
i.i.d.Independent and identically distributed
AWGNAdditive white Gaussian noise

Appendix A. Proof of Lemma 1

According to h l k 2 Γ K l k , θ l k , η l ^ k = β l ^ k h l ^ k 2 and the characteristics of the Gamma distribution, we can obtain η l ^ k Γ K l ^ k , β l ^ k θ l ^ k ; thus, η l ^ k Γ K l ^ k , Q l ^ k , where K l ^ k = m μ 1 + κ 2 m + μ κ 2 + 2 m κ = μ 1 + κ 2 1 + μ κ 2 μ κ 2 m m + 2 κ and Q l ^ k = β l ^ k θ l ^ k = β l ^ k Ω l ^ k K l ^ k = Ω ^ l ^ k K l ^ k . Lemma 1 is proven.

Appendix B. Proof of Theorem 1

In reference [57], the authors derived the PDF of the sum of i.n.i.d. generalized Gamma RVs. Therefore, we transform the PDF of η l ^ k into the PDF of generalized Gamma RVs, and its expression is as follows.
f η l ^ k x ; α , κ , μ , m , Ω ^ l ^ k = f ζ m , k x ; A l ^ k , K l ^ k , Ω ^ l ^ k = A l ^ k D l ^ k A l ^ k K l ^ k Ω ^ l ^ k A l ^ k K l ^ k Γ K l ^ k x A l ^ k K l ^ k 1 exp D l ^ k x Ω ^ l ^ k A l ^ k
where D l ^ k = Γ K l ^ k + 1 1 A l ^ k A l ^ k Γ K l ^ k . For l ^ L k and k 1 , K , A l ^ k = 1 can be obtained.
Therefore, when l ^ L k , k 1 , K and A l ^ k is constant, the sum of L ^ i.n.i.d. generalized Gamma RVs can be approximated by the generalized Puiseux series expansion. The approximation for the PDF of the sum of L i.n.i.d. κ μ shadowed RVs z is given by
f Z z = i = 1 L A D i A K i Ω i , 2 A K i Γ K i z A i = 1 L K i 1 i = 1 L Γ A K i Γ A i = 1 L K i I z ; u 1 , u 2 , , u L 1 = A L z A i = 1 L K i 1 I z ; u 1 , u 2 , , u L 1
where
A L = i = 1 L A D i A K i Ω i , 2 A K i Γ K i i = 1 L Γ A K i Γ A i = 1 L K i
I z ; u 1 , u 2 , , u L 1 = 0 1 0 1 0 1 k = 1 L 1 u k A j = 1 k K j 1 1 u k A K k + 1 1 B A i = 1 k K k , A K k + 1 Ξ ψ u 1 , u 2 , u L 1 d u 1 d u 2 d u L 1
Ξ = exp z A
ψ u 1 , u 2 , u L 1 = D 1 u 1 Ω 1 , 2 A + D 2 1 u 1 Ω 2 , 2 A u 2 A + D 3 1 u 2 Ω 3 , 2 A u 3 A + D 4 1 u 3 Ω 4 , 2 A u L 1 A + D L 1 u L 1 Ω L , 2 A
On the basis of the generalized Puiseux series expansion process, I z ; u 1 , u 2 , , u L 1 can be approximated by the truncation order of any integer V, namely,
I z ; u 1 , u 2 , , u L 1 = v = 1 V c v Ξ σ v
When V , the series can accurately approach I z ; u 1 , u 2 , , u L 1 . Moreover, the series expansion converges very quickly. Equation (A2) can be expressed as
f Z z v = 1 V c v A L z A i = 1 L K i 1 Ξ σ v v = 1 V c v A L z A i = 1 L K i 1 exp σ v z A
To simplify the derivation, we transform the content of the cumulative sum into the PDF of the generalized Gamma distribution as follows.
f Z z v = 1 V c v A L z A i = 1 L K i 1 exp σ v z A v = 1 V c v i = 1 L A D i A K i Ω i , 2 A K i Γ K i i = 1 L Γ A K i Γ A i = 1 L K i z A i = 1 L K i 1 exp σ v z A v = 1 V c v i = 1 L A D i A K i Ω i , 2 A K i Γ K i i = 1 L Γ A K i Γ A i = 1 L K i ω v Ω z , 2 A z K z Γ K z A z D z A z K z A z D z A z K z ω v Ω z , 2 A z K z Γ K z z A K z 1 exp D z ω v Ω z , 2 z A v = 1 V c v A z D z A z K z ω v Ω z , 2 A z K z Γ K z z A z K z 1 exp D z z ω v Ω z , 2 A
Additionally, according to the nature of the Gamma distribution, the corresponding CDF of the approximation can be estimated as
F Z z v = 1 V c v γ K , D z z ω v Ω z , 2 A Γ K z
The coefficients c m and ω m can be obtained by using the moment-matching method [57]. Using this approach, the moments of the approximate PDF are equal to the exact moments. For this purpose, the following equations are formulated, in which the approximate moments appear on the left and the exact moments appear on the right.
v = 1 V c v = 1 , v = 1 V c v ω v Ω z , 2 = E Z = E X 1 + X 2 + + X L = i = 1 L E X i v = 1 V c v ω v 2 Ω z , 2 2 Γ K z Γ K z + 1 1 A z A z 2 Γ K z + 2 2 A z A z Γ K z = E Z 2 v = 1 V c v ω v 2 V 2 Ω z , 2 2 V 2 Γ K z Γ K z + 1 1 A z A z 2 V 2 Γ K z + 2 M z 2 2 M z 2 A z A z Γ K z = E Z 2 V 2 v = 1 V c v ω v A z K z A z D z A z K z Ω z , 2 A z K z Γ K z = v = 1 V c v A L = A L
Theorem 1 is proven.

Appendix C. Proof of Theorem 2

The Mellin transform over the service process of user k at one time slot is
M ϕ k 1 s = E e 1 s 1 r k γ ^ k , ε k d = 0 ε k d + 1 ε k d e s r k γ ^ k , ε k d f γ x d x = ε k d + 1 ε k d γ 0 g r k γ ^ k , ε k d s ln 2 f γ ^ k x d x + 1 ε k d 0 γ 0 g r k γ ^ k , ε k d s ln 2 f γ ^ k x d x = c ε k d + 1 ε k d γ 0 f r k γ ^ k , ε k d s ln 2 f γ ^ k x d x + 1 ε k d 0 γ 0 f γ ^ k x d x = ε k d + 1 ε k d γ 0 f r k γ ^ k , ε k d s ln 2 f γ ^ k x d x + 1 ε k d F γ ^ k γ 0 F γ ^ k 0 ε k d + 1 ε k d 0 f r k γ ^ k , ε k d s ln 2 f γ ^ k x d x + 1 ε k d F γ ^ k γ 0 = ε k d + 1 ε k d 0 1 + γ ^ k exp γ ^ k 2 + 2 γ ^ k 1 + γ ^ k 2 n ^ Q 1 ε k d s ln 2 f γ ^ k x d x + 1 ε k d F γ ^ k γ 0 .
When γ ^ k γ 0 , g r k γ ^ k , ε k d = f r k γ ^ k , ε k d , γ ^ k > γ 0 1 , γ ^ k γ 0 , (c) is established.
When γ ^ k 1 , V γ ^ k = γ ^ k 2 + 2 γ ^ k 1 + γ ^ k 2 1 , we can obtain γ ^ k 2 + 2 γ ^ k 1 + γ ^ k 2 n ^ Q 1 ε 0 and 1 + γ ^ k exp γ ^ k 2 + 2 γ ^ k 1 + γ ^ k 2 n ^ Q 1 ε γ ^ k [22].
The Mellin transform over the service process in the SNR domain of the kth user at one time slot can be rewritten as
M ϕ k 1 s ε k d + 1 ε k d 0 1 + γ ^ k exp γ ^ k 2 + 2 γ ^ k 1 + γ ^ k 2 n ^ Q 1 ε s ln 2 f γ ^ k x d x + 1 ε k d F γ ^ k γ 0 d ε k d + 1 ε k d 0 γ ^ k s ln 2 f γ ^ k x d x + 1 ε k d F γ ^ k γ 0 = ε k d + 1 ε k d 0 e s ln 2 ln γ ^ k f γ ^ k x d x + 1 ε k d F γ ^ k γ 0 e ε k d + 1 ε k d 0 e s ln 2 γ ^ k 1 f γ ^ k x d x + 1 ε k d F γ ^ k γ 0 f ε k d + 1 ε k d M γ k s ln 2 + 1 ε k d F γ ^ k γ 0 = I k s , ε k d + G k γ 0 , ε k d
where I k s , ε k d = 1 ε k d M γ ^ k s ln 2 , G k γ 0 , ε k d = ε k d + 1 ε k d F γ ^ k γ 0 . (d) and (f) are obtained from γ ^ k 1 . (e) is approximately obtained in terms of the first-order expansion of the Taylor formula. Theorem 2 is proven.

Appendix D. Proof of Corollary 2

The approximate UB-QDVP of the kth user is
ε k v , UB = Pr w k t > t T M ϕ k 1 s t T 1 M φ k 1 + s M ϕ k 1 s M ϕ k 1 s t T 1 e λ s M ϕ k 1 s I k s , ε k d + G k γ 0 , ε k d t T 1 e λ s I k s , ε k d + G k γ 0 , ε k d
Corollary 2 is proven.

References

  1. Nguyen, D.C.; Ding, M.; Pathirana, P.N.; Seneviratne, A.; Li, J.; Niyato, D.; Dobre, O.; Poor, H.V. 6G Internet of Things: A comprehensive survey. IEEE Internet Things J. 2022, 9, 359–383. [Google Scholar] [CrossRef]
  2. Berry, R.A. Optimal power-delay tradeoffs in fading channels-Small-delay asymptotics. IEEE Trans. Inf. Theory 2013, 59, 3939–3952. [Google Scholar] [CrossRef]
  3. Buzzi, S.; Chih-Lin, I.; Klein, T.E.; Poor, H.V.; Yang, C.; Zappone, A. A survey of energy-efficient techniques for 5G networks and challenges ahead. IEEE J. Sel. Areas Commun. 2016, 34, 697–709. [Google Scholar] [CrossRef] [Green Version]
  4. Papazafeiropoulos, A.; Ngo, H.Q.; Kourtessis, P.; Chatzinotas, S.; Senior, J.M. Towards optimal energy efficiency in cell-free massive MIMO systems. IEEE Trans. Green Commun. Netw. 2021, 5, 816–831. [Google Scholar] [CrossRef]
  5. Ye, N.; Li, X.; Yu, H.; Zhao, L.; Liu, W.; Hou, X. DeepNOMA: A unified framework for NOMA using deep multi-task learning. IEEE Trans. Wirel. Commun. 2020, 19, 2208–2225. [Google Scholar] [CrossRef]
  6. Pan, J.; Ye, N.; Yu, H.; Hong, T.; Al-Rubaye, S.; Mumtaz, S.; Al-Dulaimi, A.; Chih-Lin, I. AI-Driven Blind Signature Classification for IoT Connectivity: A Deep Learning Approach. IEEE Trans. Wirel. Commun. 2022. early access. [Google Scholar] [CrossRef]
  7. Ye, N.; An, J.; Yu, J. Deep-Learning-Enhanced NOMA Transceiver Design for Massive MTC: Challenges, State of the Art, and Future Directions. IEEE Wirel. Commun. 2021, 28, 66–73. [Google Scholar] [CrossRef]
  8. Ye, N.; Yu, J.; Wang, A.; Zhang, R. Help from space: Grant-free massive access for satellite-based IoT in the 6G era. Digit. Commun. Netw. 2021, 8, 215–224. [Google Scholar] [CrossRef]
  9. Ding, Z.; Liu, Y.; Choi, J.; Sun, Q.; Elkashlan, M.; Chih-Lin, I.; Poor, H.V. Application of non-orthogonal multiple access in LTE and 5G networks. IEEE Commun. Mag. 2017, 55, 185–191. [Google Scholar] [CrossRef] [Green Version]
  10. Ren, H.; Pan, C.; Deng, Y.; Elkashlan, M.; Nallanathan, A. Joint power and blocklength optimization for URLLC in a factory automation scenario. IEEE Trans. Wirel. Commun. 2019, 19, 1786–1801. [Google Scholar] [CrossRef] [Green Version]
  11. Liu, Y.; Qin, Z.; Elkashlan, M.; Ding, Z.; Nallanathan, A.; Hanzo, L. Non-orthogonal multiple access for 5G and beyond. Proc. IEEE 2017, 105, 2347–2381. [Google Scholar] [CrossRef] [Green Version]
  12. Ngo, H.Q.; Ashikhmin, A.; Yang, H.; Larsson, E.G.; Marzetta, T.L. Cell-free massive MIMO versus small cells. IEEE Trans. Wirel. Commun. 2017, 16, 1834–1850. [Google Scholar] [CrossRef] [Green Version]
  13. Zhang, J.; Chen, S.; Lin, Y.; Zheng, J.; Ai, B.; Hanzo, L. Cell-free massive MIMO: A new next-generation paradigm. IEEE Access 2019, 7, 99878–99888. [Google Scholar] [CrossRef]
  14. Björnson, E.; Sanguinetti, L. Making cell-free massive MIMO competitive with MMSE processing and centralized implementation. IEEE Trans. Wirel. Commun. 2019, 19, 77–90. [Google Scholar] [CrossRef] [Green Version]
  15. Interdonato, G.; Frenger, P.; Larsson, E.G. Scalability aspects of cell-free massive MIMO. In Proceedings of the IEEE International Conference on Communications (ICC), Shanghai, China, 20–24 May 2019. [Google Scholar]
  16. Ammar, H.A.; Adve, R.; Shahbazpanahi, S.; Boudreau, G.; Srinivas, K.V. User-centric cell-free massive MIMO networks: A survey of opportunities, challenges and solutions. IEEE Commun. Surv. Tutor. 2022, 24, 611–652. [Google Scholar] [CrossRef]
  17. Jin, S.N.; Yue, D.W.; Nguyen, H.H. Spectral and energy efficiency in cell-free massive MIMO systems over correlated Rician fading. IEEE Syst. J. 2020, 15, 2822–2833. [Google Scholar] [CrossRef]
  18. Wang, X.; Ashikhmin, A.; Wang, X. Wirelessly powered cell-free IoT: Analysis and optimization. IEEE Internet Things J. 2020, 7, 8384–8396. [Google Scholar] [CrossRef]
  19. Elwekeil, M.; Zappone, A.; Buzzi, S. Power control in cell-free massive MIMO networks for UAVs URLLC under the finite blocklength regime. arXiv 2021, arXiv:2111.10613. [Google Scholar]
  20. Zhang, J.; Fan, J.; Zhang, J.; Ng, D.W.K.; Sun, Q.; Ai, B. Performance analysis and optimization of NOMA-based cell-free massive MIMO for IoT. IEEE Internet Things J. 2021, 9, 9625–9639. [Google Scholar] [CrossRef]
  21. Abdallah, A.; Mansour, M.M. Efficient angle-domain processing for FDD-based cell-free massive MIMO systems. IEEE Trans. Commun. 2020, 68, 2188–2203. [Google Scholar] [CrossRef] [Green Version]
  22. Zhang, X.; Wang, J.; Poor, H.V. Statistical delay and error-rate bounded QoS provisioning for mURLLC over 6G CF M-MIMO mobile networks in the finite blocklength regime. IEEE J. Sel. Areas Commun. 2020, 39, 652–667. [Google Scholar] [CrossRef]
  23. Saad, W.; Bennis, M.; Chen, M. A vision of 6G wireless systems: Applications, trends, technologies, and open research problems. IEEE Netw. 2020, 34, 134–142. [Google Scholar] [CrossRef] [Green Version]
  24. Dang, S.; Amin, O.; Shihada, B.; Alouini, M.S. What should 6G be? Nat. Electron. 2020, 3, 20–29. [Google Scholar] [CrossRef] [Green Version]
  25. Ramirez-Espinosa, P.; Lopez-Martinez, F.J.; Paris, J.F.; Yacoub, M.D.; Martos-Naya, E. An extension of the κ-μ shadowed fading model: Statistical characterization and applications. IEEE Trans. Veh. Technol. 2018, 67, 3826–3837. [Google Scholar] [CrossRef] [Green Version]
  26. Chen, Z.; Wang, L.; Han, D.; Zeng, H. A unified performance analysis of relaying communication system for IoT application with hybrid fading. IEEE Internet Things J. 2020, 7, 570–583. [Google Scholar] [CrossRef]
  27. Moreno-Pozas, L.; Lopez-Martinez, F.J.; Paris, J.F.; Martos-Naya, E. The κ-μ shadowed fading model: Unifying the κ-μ and η-μ distributions. IEEE Trans. Veh. Technol. 2016, 65, 9630–9641. [Google Scholar] [CrossRef] [Green Version]
  28. Paris, J.F. Statistical characterization of κ-μ shadowed fading. IEEE Trans. Veh. Technol. 2013, 63, 518–526. [Google Scholar] [CrossRef] [Green Version]
  29. Yacoub, M.D. The αηκμ Fading Model. IEEE Trans. Antennas Propag. 2016, 64, 3597–3610. [Google Scholar] [CrossRef]
  30. Al-Hmood, H.; Al-Raweshidy, H.S. On the effective rate and energy detection based spectrum sensing over αηκμ fading channels. IEEE Trans. Veh. Technol. 2020, 69, 9112–9116. [Google Scholar] [CrossRef]
  31. Lopez-Martinez, F.J.; Paris, J.F.; Romero-Jerez, J.M. The κ-μ shadowed fading model with integer fading parameters. IEEE Trans. Veh. Technol. 2017, 66, 7653–7662. [Google Scholar] [CrossRef] [Green Version]
  32. Sun, C.; She, C.; Yang, C.; Quek, T.Q.; Li, Y.; Vucetic, B. Optimizing resource allocation in the short blocklength regime for ultra-reliable and low-latency communications. IEEE Trans. Wirel. Commun. 2018, 18, 402–415. [Google Scholar] [CrossRef]
  33. Polyanskiy, Y.; Poor, H.V.; Verdú, S. Channel coding rate in the finite blocklength regime. IEEE Trans. Inf. Theory 2010, 56, 2307–2359. [Google Scholar] [CrossRef]
  34. Zhang, X.; Wang, J.; Poor, H.V. AoI-driven statistical delay and error-rate bounded QoS provisioning for mURLLC over UAV-multimedia 6G mobile networks using FBC. IEEE J. Sel. Areas Commun. 2021, 39, 3425–3443. [Google Scholar] [CrossRef]
  35. Bennis, M.; Debbah, M.; Poor, H.V. Ultrareliable and low-latency wireless communication: Tail, risk, and scale. Proc. IEEE 2018, 106, 1834–1853. [Google Scholar] [CrossRef] [Green Version]
  36. She, C.; Yang, C.; Quek, T.Q. Radio resource management for ultra-reliable and low-latency communications. IEEE Commun. Mag. 2017, 55, 72–78. [Google Scholar] [CrossRef]
  37. Al-Zubaidy, H.; Liebeherr, J.; Burchard, A. Network-layer performance analysis of multihop fading channels. IEEE/ACM Trans. Netw. 2014, 24, 204–217. [Google Scholar] [CrossRef]
  38. Nasir, A.A.; Tuan, H.D.; Ngo, H.Q.; Duong, T.Q.; Poor, H.V. Cell-free massive MIMO in the short blocklength regime for URLLC. IEEE Trans. Wirel. Commun. 2021, 20, 5861–5871. [Google Scholar] [CrossRef]
  39. Alonzo, M.; Baracca, P.; Khosravirad, S.R.; Buzzi, S. Cell-free and user-centric massive MIMO architectures for reliable communications in indoor factory environments. IEEE Open J. Commun. Soc. 2021, 2, 1390–1404. [Google Scholar] [CrossRef]
  40. Zhang, X.; Wang, J.; Poor, H.V. Statistical delay and error-rate bounded QoS provisioning for SWIPT over CF M-MIMO 6G mobile wireless networks using FBC. IEEE J. Sel. Top. Signal Process. 2021, 15, 1272–1287. [Google Scholar] [CrossRef]
  41. ElHalawany, B.M.; Jameel, F.; Da Costa, D.B.; Dias, U.S.; Wu, K. Performance analysis of downlink NOMA systems over κ-μ shadowed fading channels. IEEE Trans. Veh. Technol. 2019, 69, 1046–1050. [Google Scholar] [CrossRef]
  42. Chun, Y.J.; Cotton, S.L.; Dhillon, H.S.; Lopez-Martinez, F.J.; Paris, J.F.; Yoo, S.K. A comprehensive analysis of 5G heterogeneous cellular systems operating over κ-μ shadowed fading channels. IEEE Trans. Wirel. Commun. 2017, 16, 6995–7010. [Google Scholar] [CrossRef] [Green Version]
  43. Sánchez, J.D.V.; Osorio, D.P.M.; López-Martínez, F.J.; Paredes, M.C.P.; Urquiza-Aguiar, L.F. Information-theoretic security of MIMO networks under κ-μ shadowed fading channels. IEEE Trans. Veh. Technol. 2021, 70, 6302–6318. [Google Scholar] [CrossRef]
  44. Schiessl, S.; Gross, J.; Skoglund, M.; Caire, G. Delay performance of the multiuser MISO downlink under imperfect CSI and finite-length coding. IEEE J. Sel. Areas Commun. 2019, 37, 765–779. [Google Scholar] [CrossRef] [Green Version]
  45. Schiessl, S.; Skoglund, M.; Gross, J. NOMA in the uplink: Delay analysis with imperfect CSI and finite-length coding. IEEE Trans. Wirel. Commun. 2020, 19, 3879–3893. [Google Scholar] [CrossRef] [Green Version]
  46. Xiao, C.; Zeng, J.; Ni, W.; Su, X.; Liu, R.P.; Lv, T.; Wang, J. Downlink MIMO-NOMA for ultra-reliable low-latency communications. IEEE J. Sel. Areas Commun. 2019, 37, 780–794. [Google Scholar] [CrossRef]
  47. Yang, G.; Xiao, M.; Poor, H.V. Low-latency millimeter-wave communications: Traffic dispersion or network densification? IEEE Trans. Commun. 2018, 66, 3526–3539. [Google Scholar] [CrossRef] [Green Version]
  48. Elhoushy, S.; Ibrahim, M.; Hamouda, W. Cell-free massive MIMO: A survey. IEEE Commun. Surv. Tutor. 2022, 24, 492–523. [Google Scholar] [CrossRef]
  49. Femenias, G.; Riera-Palou, F. Fronthaul-constrained cell-free massive MIMO with low resolution ADCs. IEEE Access 2020, 8, 116195–116215. [Google Scholar] [CrossRef]
  50. Becvar, Z.; Mach, P.; Elfiky, M.; Sakamoto, M. Hierarchical scheduling for suppression of fronthaul delay in C-RAN with dynamic functional split. IEEE Commun. Mag. 2021, 59, 95–101. [Google Scholar] [CrossRef]
  51. Van Chien, T.; Björnson, E.; Larsson, E.G. Joint power allocation and load balancing optimization for energy-efficient cell-free massive MIMO networks. IEEE Trans. Wirel. Commun. 2020, 19, 6798–6812. [Google Scholar] [CrossRef]
  52. Fodor, G.; Di Marco, P.; Telek, M. On minimizing the MSE in the presence of channel state information errors. IEEE Commun. Lett. 2015, 19, 1604–1607. [Google Scholar] [CrossRef]
  53. Liu, P.; Luo, K.; Chen, D.; Jiang, T. Spectral efficiency analysis of cell-free massive MIMO systems with zero-forcing detector. IEEE Trans. Wirel. Commun. 2019, 19, 795–807. [Google Scholar] [CrossRef] [Green Version]
  54. Arnau, J.; Kountouris, M. Delay performance of MISO wireless communications. In Proceedings of the 2018 16th International Symposium on Modeling and Optimization in Mobile, Ad Hoc, and Wireless Networks (WiOpt), Shanghai, China, 7–11 May 2018. [Google Scholar]
  55. She, C.; Yang, C.; Quek, T.Q. Cross-layer optimization for ultra-reliable and low-latency radio access networks. IEEE Trans. Wirel. Commun. 2017, 17, 127–141. [Google Scholar]
  56. Ngo, H.Q.; Larsson, E.G.; Marzetta, T.L. Energy and spectral efficiency of very large multiuser MIMO systems. IEEE Trans. Commun. 2013, 61, 1436–1449. [Google Scholar]
  57. Payami, M.; Falahati, A. Accurate variable-order approximations to the sum of αμ variates with application to MIMO systems. IEEE Trans. Wirel. Commun. 2020, 20, 1612–1623. [Google Scholar] [CrossRef]
Figure 1. CF mMIMO system under the κ μ shadowed fading model.
Figure 1. CF mMIMO system under the κ μ shadowed fading model.
Sensors 22 05279 g001
Figure 2. Relationship between the decoding error probability and the maximum achievable data rate.
Figure 2. Relationship between the decoding error probability and the maximum achievable data rate.
Sensors 22 05279 g002
Figure 3. PDFs of the sum of i.n.i.d. κ μ shadowed RVs. (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4.
Figure 3. PDFs of the sum of i.n.i.d. κ μ shadowed RVs. (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4.
Sensors 22 05279 g003
Figure 4. PDFs of the PPSNR for CF mMIMO system. (a) under the Rician fading channel; (b) under the Rayleigh fading channel; (c) under the One-sided fading channel; (d) under the Nakagami-m fading channel; (e) under the Rician shadwoed fading channel; (f) under the κ μ shadowed fading channel.
Figure 4. PDFs of the PPSNR for CF mMIMO system. (a) under the Rician fading channel; (b) under the Rayleigh fading channel; (c) under the One-sided fading channel; (d) under the Nakagami-m fading channel; (e) under the Rician shadwoed fading channel; (f) under the κ μ shadowed fading channel.
Sensors 22 05279 g004
Figure 5. Performance comparison conducted for the CF mMIMO system under the κ μ shadowed fading model. (a) Comparison between Rician and Rayleigh fading channel; (b) Comparison between Rayleigh and One-sided fading channel; (c) Comparison between One-sided and Nakagami-m fading channel; (d) Comparison between Nakagami-m and Rician shadowed fading channel; (e) Comparison between Rician shadowed and the κ μ shadowed fading channel; (f) Comparison between the κ μ shadowed and Rician fading channel.
Figure 5. Performance comparison conducted for the CF mMIMO system under the κ μ shadowed fading model. (a) Comparison between Rician and Rayleigh fading channel; (b) Comparison between Rayleigh and One-sided fading channel; (c) Comparison between One-sided and Nakagami-m fading channel; (d) Comparison between Nakagami-m and Rician shadowed fading channel; (e) Comparison between Rician shadowed and the κ μ shadowed fading channel; (f) Comparison between the κ μ shadowed and Rician fading channel.
Sensors 22 05279 g005
Figure 6. Comparison of the UB-QDVPs in the CF mMIMO, PD-NOMA and OMA systems. (a) Comparison of different communication systems under One-sided Gaussian fading; (b) Comparison of different communication systems under Nakagami-m fading.
Figure 6. Comparison of the UB-QDVPs in the CF mMIMO, PD-NOMA and OMA systems. (a) Comparison of different communication systems under One-sided Gaussian fading; (b) Comparison of different communication systems under Nakagami-m fading.
Sensors 22 05279 g006
Figure 7. The UB-QDVP vs. E E dB under the κ μ shadowed fading model. (a) Comparison of different small-scale fading types; (b) Comparison of different communication systems under Rician shadowed fading; (c) Comparison of different communication systems under Rayleigh fading. (d) Comparison of different communication systems under κ μ shadowed fading.
Figure 7. The UB-QDVP vs. E E dB under the κ μ shadowed fading model. (a) Comparison of different small-scale fading types; (b) Comparison of different communication systems under Rician shadowed fading; (c) Comparison of different communication systems under Rayleigh fading. (d) Comparison of different communication systems under κ μ shadowed fading.
Sensors 22 05279 g007
Table 1. Parameters Values for Numerical Results.
Table 1. Parameters Values for Numerical Results.
ParameterValue
Farthest coverage radius (m)1000
Number of potential users2000
Activation rate of potential users5%
Number of APs200
Noise power spectral density (dBm/Hz)−106
Path loss exponent, α p 3
Minimum distance, d 0 (m)10
Constant path loss, u p (dB)−10
Bandwidth, B (MHz)20
Carrier frequency (GHz)3.5
Arrival bit per slot, λ (bit/slot)50
Length of time slot, t f (ms)0.01
Length of pilots, n Max ( K , 0.2 N )  1
Target delay (target number of time slots), t T [1,10]
Decoding error probability, ε k d [ 10 9 , 10 8 ]
1 This is the fixed ratio pilot allocation method [56].
Table 2. The second-order approximate parameters in different cases shown in Figure 3.
Table 2. The second-order approximate parameters in different cases shown in Figure 3.
c 1 c 2 ω 1 ω 2
Case 1:
κ = 4 , μ = 1 , m
κ = 4 , μ = 1 , m = 3
0.99980.00020.997616.6538
Case 2:
κ = 4 , μ = 1 , m
κ = 4 , μ = 1 , m = 3
κ = 6 , μ = 1 , m = 3
0.99850.00150.99235.9651
Case 3:
κ = 4 , μ = 1 , m
κ = 4 , μ = 1 , m = 3
κ = 6 , μ = 1 , m = 3
κ = 4 , μ = 2 , m
0.86050.13950.93551.3979
Case 4:
κ = 4 , μ = 1 , m
κ = 4 , μ = 1 , m = 3
κ = 6 , μ = 1 , m = 3
κ = 4 , μ = 2 , m
κ = 4 , μ = 2 , m = 3
0.89020.10980.94861.4166
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zeng, J.; Song, Y.; Wu, T.; Lv, T.; Zhou, S. Guaranteeing QoS for NOMA-Enabled URLLC Based on κμ Shadowed Fading Model. Sensors 2022, 22, 5279. https://doi.org/10.3390/s22145279

AMA Style

Zeng J, Song Y, Wu T, Lv T, Zhou S. Guaranteeing QoS for NOMA-Enabled URLLC Based on κμ Shadowed Fading Model. Sensors. 2022; 22(14):5279. https://doi.org/10.3390/s22145279

Chicago/Turabian Style

Zeng, Jie, Yuxin Song, Teng Wu, Tiejun Lv, and Shidong Zhou. 2022. "Guaranteeing QoS for NOMA-Enabled URLLC Based on κμ Shadowed Fading Model" Sensors 22, no. 14: 5279. https://doi.org/10.3390/s22145279

APA Style

Zeng, J., Song, Y., Wu, T., Lv, T., & Zhou, S. (2022). Guaranteeing QoS for NOMA-Enabled URLLC Based on κμ Shadowed Fading Model. Sensors, 22(14), 5279. https://doi.org/10.3390/s22145279

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop