Radio Resource Dimensioning for Low Delay Access in Licensed OFDMA IoT Networks †

In this paper, we focus on the radio resource planning in the uplink of licensed Orthogonal Frequency Division Multiple Access (OFDMA) based Internet of Things (IoT) networks. The average behavior of the network is considered by assuming that active sensors and collectors are distributed according to independent random Poisson Point Process (PPP) marked by channel randomness. Our objective is to statistically determine the optimal total number of Radio Resources (RRs) required for a typical cell. On one hand, the allocated bandwidth should be sufficiently large to support the traffic of the devices and to guarantee a low access delay. On the other hand, the over-dimensioning is costly from an operator point of view and induces spectrum wastage. For this sake, we propose statistical tools derived from stochastic geometry to evaluate, adjust and adapt the allocated bandwidth according to the network parameters, namely the required Quality of Service (QoS) in terms of rate and access delay, the density of the active sensors, the collector intensities, the antenna configurations and the transmission modes. The optimal total number of RRs required for a typical cell is then calculated by jointly considering the constraints of low access delay, limited power per RR, target data rate and network outage probability. Different types of networks are considered including Single Input Single Output (SISO) systems, Single Input Multiple Output (SIMO) systems using antenna selection or Maximum Ratio Combiner (MRC), and Multiuser Multiple Input Multiple Output (MU-MIMO) systems using Zero-Forcing decoder.


Introduction
Cellular licensed IoT technology has been an emerging and evolving Low Power Wide Area (LPWA) technology which provides long range, low power and low cost connectivity for IoT devices [1]. It can be deployed in existing cellular networks from which it inherits many of the features that determine its behavior [2], such as Long Term Evolution (LTE) networks. With the booming development of Fifth Generation (5G) technology, evolutions of cellular IoT standards have been proposed and put into practice by 3rd Generation Partnership Project (3GPP) [3]. It will further promote the development of IoT technology. 5G technologies primarily include non-standalone and standalone technologies. The former is the early-drop technology that operators plan to connect new radio 5G base stations to the Fourth Generation (4G) core network until 5G connectivity migrates λ a and λ b . We assume that the sensors are equipped with a single antenna and that the collector is equipped with n r receiver antennas. We assume that a sensor is active n a times per day, the mean service time is ν −1 (s) and the interarrival rate is n a 24×60×60 ρ per second and per km 2 . The active sensors nodes form then a spatial PPP Φ a with intensity λ a = n a 24 × 60 × 60 ρν −1 .
The frequency reuse pattern in the network is equal to 1. At a given collector situated at y 0 ∈ Φ b , the received power on a given RR from a sensor x ∈ Φ a transmitting with power P RR (mW per RR) is computed as, P r (y 0 , x) = P RR α |y 0 − x| −β A f where α and β are respectively the attenuation factor (that includes the average shadowing) and the path-loss exponent that are computed from the Okumura-Hata model, A f is the fading coefficient with distribution depending on the antenna configuration and the used transmission mode (single-user or multiuser). The sensor nodes are considered as static and the fading channel is considered as flat during the transmission. The Channel State Information (CSI) is only available at the receiver side and not at the transmitter side. Each sensor is connected to the collector on which the average received power (by averaging over the fading and the shadowing) is the highest. This is equivalent to connect the sensor to the nearest collector. Assuming a collector y 0 , the set of the sensors connected to this latter is defined as Φ c (y 0 ) = {x ∈ Φ : ∀y ∈ Φ b − {y 0 } : |y 0 − x| < |y − x|}. Nodes transmitting in the same frequency band generate additive interference with power with Φ I (y 0 ) = ∪ y∈Φ b ,y =y 0 {x i : x i randomly selected in Φ c (y)} .
The received SINR at the given collector y 0 is, SINR(y 0 , x) = P RR α |y 0 − x| −β A f P n + I(y 0 ) , (4) with P n being the random exponential noise power with mean power ofP n = KTB where K is the Boltzmann's constant K = 1.379 × 10 −23 W Hz −1 K −1 , T the absolute temperature in kelvins T = 290 K and B the bandwidth. The power of the random exponential noise power is characterized by its Laplace transform as, L P n (s) = E e −sN = 1 sP n + 1 .
Due to Slivnyak-Mecke Theorem in [8], the statistical behavior in the PPP remains unchanged when adding a collector at the center of this region. This defines the typical cell centered at the origin Φ c (0) as illustrated in Figure 1. The typical cell properties are reviewed in Proposition 1.
Proposition 1 (Typical cell properties). The typical cell average number of sensors is, The probability distribution function of r = |x| the distances between x ∈ Φ c (0) and the o collector is, The average radius of the typical cell is 1/(2 √ λ b ).
Proof. We provide for completeness the proof in Appendix A.

Proposed Statistical Dimensioning Model
Given the density of collectors in the network, we propose a statistical method for resource planning in the uplink of cellular IoT networks. We assume the network access is an OFDMA and consider only the Narrowband Physical Uplink Shared Channel (NPUSCH) which carries the uplink user data and uplink control information. In an NB-IoT network or LTE-M, identical to the LTE network, the smallest block of radio resource elements that can be invoked is called the Resource Block (RB) or Radio Resource Block (RRB). It contains several OFDM symbols. Each RR corresponds to 180 kHz in the frequency domain and 0.5 ms in the time domain. In the 5G, different numerologies were defined to reduce the RR duration and hence the latency by increasing the RR bandwidth. In the following, the RR corresponds to the LTE resource block and the RR with numerology 0 in the 5G terminology. For the single-user and multiuser cases, our main objective here is to statistically determine the total number of required RRs to minimize the occurrence of the network outage event. When the total number of RRs required is larger than the number of RRs available on the collector side, this is considered as an outage event.

Dimensioning Objectives in a Typical Cell
In the typical cell illustrated in Figure 1, the collector o allocates according to the level of the received SINR, 1 to N max RRs to the sensor node in order to achieve its target data rate C 0 . The required number of RRs is, with γ k being the threshold SINR required to achieve a target rate of C 0 /k within a single RR and γ 0 > γ 1 > . . . > γ N max and γ 0 → ∞. Note that γ 0 (respectively γ N max ) is the threshold SINR to decode the highest (respectively lowest) Modulation and Coding (MCS) scheme. If SINR(0, x) < γ N max , the sensor will not be able to decode the lowest MCS and the collector does not attribute any RR to this user. This event of having insufficient SINR occurs with probability, we intentionally do not refer to this probability as outage to avoid confusion with the network outage probability.
By setting with A f ,0 = ∞ the threshold fading to achieve a rate C k ∈ [C 0 /k ; C 0 /(k − 1)] over one RR, the number of RR can be rewritten as, The total of required RR in this typical cell depicted in Figure 1 is, The network is in outage if, In order to ensure an optimized network dimensioning, the number of total radio resources N t that ensure a network outage probability of p th,n , should be found. Using the concentration inequality, the typical cell outage probability is upper-bounded by, where with the function g(t) = (1 + t) log(1 + t) − t and N t > m N . By setting a threshold p th,n on the network outage,

Average Delay and Choice of the Network Threshold
In this subsection, we compute the average delay to connect the sensor device to the network independently of its transmitting rate. This average delay should not exceed a maximal delay τ max that we assume proportional to the Time Transmission Interval (TTI). As long as the sensor is not accepted by the network due to the lack of resource, a new trial will be performed after TTI. We assume that the probabilities of being rejected after i × TTI are independent and are equal to P out,c (N t ) i . The average delay to access to the network is hence, This expression can be simplified to,τ ≤ TTI p th,n (1 − p th,n ) 2 . (19) By choosing p th,n such that we make sure that the average delay in the network does not exceed τ max . By letting κ = τ max /TTI, the network threshold is,

Expressions of m N and v N
In this subsection, we derive the expressions of m N and v N in (16) and (17) considering the single-user and multiuser cases. We recall that all sensor devices are equipped with a single antenna and the number of antennas at the collector is n r . The distribution of the fading coefficient in (4) will be specified in the next section depending on the antenna configuration and transmission mode.

Single-User Case
Proposition 2 (Single-user case). The expressions of m N and v N are, with N s = λ a /λ b the average number of sensors as shown in (6), and E r [(.)] = ∞ 0 (.) f (r)dr with f (r) defined in (7).
Proof. Please refer to Appendix B.
Note that (22) can be interpreted as following: N s is the average number of sensor in the typical cell; each sensor requires each 1 to N max RRs depending on their distance to the collector, the channel conditions and the additive interference and noise.

Multiuser Case
In the multiuser case, we assume that each RR is shared by n u users simultaneously as depicted in Figure 2. To cancel the interuser interference, the collectors uses a Zero-Forcing (ZF) decoder.
The equivalent fading observed by each sensor j device is denoted by A (j) f having a distribution that will be specified in the next section. Distance based multiuser scheduling scheme: In this scheme, the active users are sorted according to their proximity to the collector into different groups of n u users, each independently of their fading. As shown in Figure 3, the first group contains the n u nearest neighbor to the typical collector, the second group the (n u + 1) th to 2n th u collector neighbor and so on. We adjust the number of allocated RRs to the furthest user from the collector in each group independently of the fading coefficient experienced by other users. The furthest node x in each group i is the i × n th u neighbor of the collector and it corresponds to the radius of the ball B(|x|) with radius |x| containing i × n u nodes.
where B(r) is the ball with radius r, and (26) with N * being the set of nonzero positive integers. The distribution of r in the typical cell is given by f (r) defined in (7).
Proof. Please refer to Appendix C.
Note that Proposition 3 is a generalization of Proposition 2 with n u = 1. For n u = 1, and the expressions of m N in (24) becomes equal to (22).

Distribution of the Sensor Power Consumption
We assume that the transmitted power is fixed to P RR per RR and the total power over N max RRs does not exceed P max = N max P RR . The power is then a discrete random variable and its Probability Mass Function (PMF) is given in Proposition 4.

Proposition 4 (Total sensor power PMF).
For the general case, the sensor power PMF is, with 1 ≤ k ≤ N max , N g the number of the groups of n u sensors in the typical cell such that, and N s = λ a /λ b .
Proof. The proof is similar to Proposition 3 and the PMF can be obtained by averaging over the whole power in the typical cell normalized by the size of this latter. In each ring delimited by the (i − 1) × n th u and i × n th u neighbor, the n u sensors adjust their power to the the i × n th u furthest neighbor.

Dimensioning Tools: Interference and Fading Characterization
In this section, we characterize the fading and interference distribution required to compute the average number of RR considering different antenna configuration using Single-User and Multiuser (MU-MIMO) communications.

Single-User: Case of Single Antenna Receiver
In this subsection, we review first from [9] the characterization of the Laplace transform of interference and the distribution of the fading in the SISO case averaged on the random noise plus interference.

Interference Laplace Transform
Considering a SISO antenna configuration, the Laplace transform of the interference is, in which Φ I is the set of the interfering nodes that transmit on the same RR as shown in (3). The set of interfering nodes Φ I cannot be considered as homogeneous PPP and its exact statistical distribution is generally a difficult problem to model. For this, we use the same approximation proposed in [7] that is shown to almost capture the statistical behavior of the interference field. Let d i be the distance between the interferer and the intended collector. The probability of finding at least one collector in the ball with radius d i centered at the interferer is (1 − e −λ b πd 2 i ). In this case, a node at distance d i from the intended collector is considered as interfering with probability of (1 − e −λ b πd 2 i ). It was shown in [7] that the effective interference observed at the tagged collector can be approximately modeled as a nonhomogeneous PPP with intensity Considering exponential power fading attenuation coefficients, the interference Laplace transform from (30) for the SISO case is, Using the Probability Generating Functional (PGFL) property, the Laplace transform is then,

Average Fading Distribution
Proposition 5 (SISO case). Given a node position r = |x|, the SISO fading distribution averaged on the random noise and interference, Proof. Please refer to Appendix D.

Single-User Case with Multiantenna Receiver
We assume a SIMO configuration depicted in Figure 4 which n r antennas are used at the receiver and only a single antenna at each sensor node. Let F 1 , . . . , F n r denote the fading coefficients between the sensor device and the n r receive antennas. Assume that the receiver antennas are sufficiently separated to assume that the random exponential fading coefficients are independent.

Antenna Selection
Given a node position r = |x|, the SIMO fading distribution with antenna selection averaged on the random noise and interference, A f = max(F 1 , . . . , F n r ) is distributed as the maximum between n r i.i.d exponential variables. The CDF function is then, Proposition 6 (Antenna selection). The interference Laplace transform with antenna selection is identical to the SISO case in (32), Given a node position r = |x|, the SIMO fading distribution with MRC averaged on the random noise and interference, Proof. The interfering signal is received on a random fading coefficient that is exponentially distributed as in the SISO case. The identify in (36) can be deduced from the distribution of A f given in (34).

Maximum Ratio Combiner (MRC)
Using a MRC, the equivalent fading coefficient A f = F 1 + . . . + F n r is the sum of n r random exponential i.i.d. random variables and it is distributed as chi-squared random variable with 2n r degrees of freedom. The CDF of A f is, Given a node position r = |x|, the SIMO fading distribution with MRC averaged on the random noise and interference, where d j f (s) ds j is the derivative of order j of the function f (s) for j ∈ N * and d 0 f (s) ds 0 = f (s).
Proof. The proof of this proposition is detailed in Appendix E.

Multiuser with Multiantenna Receiver
In order to cancel the multiuser interference, a Zero-Forcing decoder is used by projecting the n r × 1 received signal on the kernel of the space formed by the n r × 1 fading vectors of the (n u − 1) other scheduled users. As stated by [25], the orthogonality constraints consume (n r − n u + 1) degrees of freedom. The receiver diversity for each sensor device is limited then to (n r − n u + 1).
To extract this diversity, we assume that the ZF decoder is followed by a MRC that combines the (n r − n u + 1) received observation. The fading coefficient A f observed by each sensor in Figure 2 is then chi-squared distributed with 2(n r − n u + 1) degrees of freedom.
Proposition 8 (MU case with ZF-MRC decoder). For the MU-MIMO case with ZF decoder followed by a MRC, the equivalent fading distribution averaged on the random noise and interference, The interference Laplace transform is approximated by, Proof. For the multiuser case, n u users in each cell transmit in the same RR. The interference becomes higher than the SISO case. As all neighboring users are multiplexed on the same RR, we approximate the sum of interference coming from the group of n u users by an interfering signal generated by a random interferer transmitting with n u P RR . We have validated this approximation using numerical results as it will be shown in the next section. The rest of the proof is similar to Proposition 7 considering 2(n r − n u + 1) degrees of freedom.

Numerical Results
In this section, we consider the uplink of a sensor network corresponding to the licensed IoT with parameters summarized in Table 1. We assume that the density of nodes ρ = 500 sensors per km 2 transmitting on average once each half-an-hour, n a = 48 during 20 seconds. The active node density is then λ a = 5.5 sensors per km 2 . We consider a cellular network with ranges between 500 m to 1.5 km corresponding to a collector density of λ b = 0.9 nodes per km 2 down to 0.1 nodes per km 2 (the ranges are obtained with a confidence margin of 95%). The maximal power is limited to 14 dBm and is uniformly distributed among the allocated RRs. Table 2 gives the matching between the SINR range with the required number of RRs to achieve a target rate of C 0 = 500 bps. This data are derived from the Link Layer Simulation (LLS) provided in [26] on the Physical Uplink Shared Channel (PUSCH) of LTE-Cat M. To evaluate network performances, we consider a path-loss in an urban/suburban environment with α = 10 −14.1 and β = 3.5. We consider the following antenna configurations: SISO, single SIMO with n r = 8, multiuser SIMO with n r = 8 and n u = 2 or n u = 4.  Table 3 indicates the percentage difference between the mean number of RRs required in the typical cell derived by using our statistical dimensioning tools and the empirical network simulation results. The mean number of required RRs using former method is denoted as m N , the simulation one is denoted as m N,s . Consequently, the percentage difference is ∆m N /m N,s = |m N − m N,s |/m N,s . We consider different collector intensities (from 0.1 to 0.9 nodes/km 2 ) and different antenna configurations as well as transmission modes. As presented in the Table 3, the comparison shows that the difference between the results obtained by our statistical tools and the simulated values is very small, which verifies the accuracy of our model and approach.  Figure 5 further illustrates the specific value of m N considering different network configurations. As shown in the figure, for the given active sensor and collector intensities λ a and λ b , respectively, assuming a maximum delay of 1 ms, the SISO system always requires the highest mean number of RRs.

Average and Total Number of RR
With the use of SIMO and MU-MIMO, multiantenna techniques achieve higher data rates through increased spectral efficiency, with a corresponding reduction in the required RRs. The multiple receive antennas system has the potential to enhance signal robustness and increase system capacity. In particular, the MU-MIMO system significantly reduces the mean number of RRs required in a typical cluster. Meanwhile, as collector intensity λ b increases, the mean number of RRs required for various configurations decreases. When λ b increases, the radius of the typical cell will decrease, while the number of active sensor nodes per unit area remains the same, which means that each collector needs to serve fewer nodes in the cell and thus the number of RRs required decreases.  Figure 6 presents the total number of required RRs in a typical cell versus the collector intensity. It maintains the same trend as in Figure 5, with MU-MIMO being the best performing case, followed by SIMO and then SISO.    For example, according to Figure 6, for the same value of λ a , λ b , we find that the N t value required for the SISO system is 19. In Figure 7, we can see that N t = 19 corresponds to a CDF value of 0.88. This means that there is a network outage probability of 1 − 0.88 = 0.12 meaning that the network cannot attributed RR to 12% of nodes in the typical cell. The corresponding average delay can be then computed using (19). Figure 8 show the actual average access delay corresponding to different antenna configurations and transmission modes for given intensities of active sensors and collectors λ a , λ b and the maximum access delay τ max = 1 ms. It can be seen from the figure that the average access delay τ is less than the maximum access delay τ max . SIMO with MRC transmission mode requires the longest access time because it has the lowest individual off probability, i.e., it serves the largest percentage of active nodes. The second highest access delay is the SIMO with antenna selection scheme. For the SISO configuration, the average access delay is relatively lower because in this mode, many nodes are in off state and not served. For the MU-MIMO system, each time there are n u nodes transmitting information simultaneously and being processed by the collector. The larger the n u is, the shorter the average access delay required.  From the figure, it can be seen that the N t derived by our statistical tool is always larger than the empirical N t , i.e., it is overdimensioning. The excess in the radio-resource provided by the theoretical model ensures that the expected data rate is achieved and the access to the cellular IoT network does not exceed the predefined average access delay.

Tolerated Delay and Overdimensioning
In Figure 10, it shows the relationship between the total number of required RRs and the maximal access delay τ max with λ a = 5.5 nodes/km 2 and λ b = 0.5 nodes/km 2 . Assuming that τ max ranges from 0.5 to 4 ms, a larger τ max means that the node can wait longer for access to the network. The average access delay τ of the active nodes must not exceed the maximal delay τ ≤ τ max . In general, a lower total number of RRs is required to response to the demands of active sensor nodes with a larger τ max .

Individual Sensor off Probability
As mentioned in Section 3.1, when the SINR received on the collector side is below the SINR threshold γ N max , no RR will be assigned to this sensor. In other words, this active sensor is forced to be off in this trial. Table 4 illustrates the individual off probability in the typical cell according to the different antenna configurations, including SISO, 1 × 8 SIMO with antenna selection and MRC, 2 users and 4 users 1 × 8 MU-MIMO. From the table, we can see the 1 × 8 SIMO with MRC outperforms other modes in the individual off probability. Although it does not reduce the number of RRs required in the typical cell significantly compared to the SIMO with antenna selection or the MU-MIMO, it does reduce the sensor individual off probability, ensuring more sensors to be served. This result gives us a new way of thinking about resource planning, i.e., we need to take into account both the total number of RRs assigned to the network and the individual off probability to improve network service.

Power Distribution
In Table 5, for a given pair of λ a and λ b , we elaborate the power distribution with the different transmission modes in the typical cell. The transmission power of each RR is fixed to P RR = 25/6 mW (6.3 dBm), so the maximal transmission power is N max P RR = 25 mW (14 dBm) with N max = 6. For each transmission mode, Table 5 shows the specific percentage share of each transmit power. Where off represents a number of nodes to which no RR is assigned and therefore have a transmission power of 0 mW. From P RR to 6P RR corresponds to the transmission power of the nodes that are assigned the corresponding number of 1 to 6 RRs, respectively. As marked in the light grey part of the table, it can be noticed that SIMO with MRC mode has the least number of off nodes, meanwhile with 99.92% of the nodes being transmitting information at P RR power. Therefore, it has the best power distribution. MU-MIMO with n u = 2 also has a good power distribution, with only 0.05% of nodes off and 99.80% of nodes transmitting at P RR . SISO has the worst performance, with up to 3.97% of nodes unable to transmit information and only 90.76% of nodes transmitting at low power. Table 5. The power distribution with the different transmission modes in the typical cell, P RR = 25/6 mW, λ a = 5.5 nodes/km 2 , λ b = 0.5 nodes/km 2 .

Conclusions
In this paper, we focus on licensed OFDMA-based IoT networks which corresponds to NB-IoT, LTE-M and to some extent 5G network using the numerology 0. We assume a wireless sensor network in which sensor nodes and collectors are distributed according a spatial PPP. A statistical method based on PPP is proposed to model and analyze the average behavior of the cellular IoT networks. This method is used for statistical resource planning to support the dominant uplink communications in the IoT networks and fulfill the network requirements in terms of limited transmission power per RR, low access delay, preset network outage probability and target data rate. Different antenna configurations and transmission modes are taken into account, e.g., SISO, SIMO, and MU-MIMO. Based on the preset tolerated average access delay and target data rate, the specific total number of RRs required for the typical cell is further calculated. Numerical results are finally given to assess our statistical analysis. We compare the theoretical results with the empirical simulation results in detail. The comparison shows that the results obtained from our statistical model are very close to the empirical results. Moreover, the results highlight and quantify the radio resources gain obtained by the receiver diversity techniques and the multiuser gains.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The Consider the set of all active nodes x ∈ Φ a generated by the homogeneous PPP, the node x is in the typical cluster Φ c (0) if the origin is o is the closest point to x among all the other collectors y ∈ Φ b , i.e., |x| < |y − x|, ∀y ∈ Φ b . We can note that the following indicator function The average number of users in this typical cell is, For a given x ∈ Φ a , we use the PGFL property on the PPP of Φ b , Note that |y − x| 2 = |y| 2 + |x| 2 − 2|y||x| cos(θ z − θ y ) where |y| and θ y (resp. r = |y| and θ x ) are the polar coordinates of z (resp. y). The condition |y − x| 2 < r 2 in (A4) is then equivalent to |y| < 2|x| cos(θ y − θ x ) with 0 ≤ (θ y − θ x ) ≤ π to guarantee a positive range of r y . The corresponding integral is then, |y−x| 2 <r 2 dy = π−θ x 0 2r cos(θ y −θ x ) 0 r y dr y dθ y = π−θ x 0 r 2 cos 2 (θ y − θ x )dθ y = πr 2 (A5) independently of x. Replacing in (A3), Equations (A4) and (A5), we have, Finally, the pdf of r can be deduced using the CDF of r,

Appendix B. Proof of Proposition 2
The average number of radio-resources in the typical cluster is, Using Campbell Theorem, with dx = rdrdθ and Using the relationships (A4) and (A5), Notice that 2πλ b r exp(−λ b πr 2 ) is nothing but the pdf of the smallest distance between the sensor and its nearest collector. This implies that,

Appendix C. Proof of Proposition 3
Considering the distance based scheduling algorithm, we adjust the dimensioning of each group with respect to the i × n th u nearest neighbor to the typical collector (i ∈ N * ). We let r i denote the distance between this i × n th u nearest neighbor x i = r i e jθ i and the typical collector. The average number of required RR in the typical cluster is the total number of RR required by this i × n th u nearest neighbors, such that, m N = E r i ,θ i ,A f ,I,P n ∑ r i :i∈N * N RR (r i , A f , I, P n ) ∏ y∈Φ b By average on the collectors distribution as in Proof in Appendix B, The average number of required RR is then, being the distribution of i × n th u nearest neighbor of the typical collector.

Appendix D. Proof of Proposition 5
Using a single antenna receiver, the SISO fading coefficients are Rayleigh distributed and A f has an exponential distribution, and, E A f [1 {A f ,k ≤A f <A f ,k−1 } ] = e −(P n +I)s k − e −(P n +I)s k−1 .
(A16) Proposition 5 is obtained considering the average behavior of the noise and interference.

Appendix E. Proof of Proposition 7
Appendix E. 1

. Laplace Transform of the Interference with MRC
To show that the Laplace transform remains the same, we let h x = [h x,1 . . . h x,n r ] denote the complex Gaussian vector containing the (signal) fading coefficients and h y = [h y,1 . . . h y,n r ] the fading vector between the interferer and the collector. When using a MRC decoder, the n r observations of x are projected into the orthonormal vector u x = h x / h x . The equivalent interferer channel is then the scalar product between u x and h y . Let u y = h y / h y , the equivalent fading seen by each interferer A f ,i = | u x , h y | 2 is exponentially distributed. This is a consequence that the cosinus square of an angle between a Gaussian vector and a given direction cos 2 (u x , h y ) is beta distributed with parameters 1 and (n r − 1) as shown in [27]. This leads to a distribution of β(1, n r − 1) multiplied by a chi-squared with 2n r degrees of freedom which is exponentially distributed. In other words, A f ,i = | u x , h y | 2 = h y 2 cos(u x , u y ) ∼ χ 2 2n r × β(1, n r − 1) = exprnd(1). (A17)

Appendix E.2. Fading Distribution
Using the distribution of A f in (37), By averaging over the noise and the interference level, Proposition 7 is obtained by invoking the following lemma: Lemma A1. Let Ξ i (s) = E s i (P n + I) i e −s(P n +I) Then, Due to the independence of P n and I, Note that, for a given random variable X, we have,