Abstract

We study the performance of the secondary relay system in a power-beacon (PB) assisted energy harvesting cognitive relay wireless network. In our system model, a secondary source node and a relay node first harvest energy from distributed PBs. Then, the source node transmits its data to the destination node with the help of the relay node. Also, fading coefficients of the links from the PBs to the source node and relay node are assumed independent but not necessarily identically distributed (i.n.i.d) Nakagami- random variables. We derive exact expressions for the power outage probability and the channel outage probability. Based on that, we analyze the total outage probability of the secondary relay system. Asymptotic analysis is also performed, which provides insights into the system behavior. Moreover, we evaluate impacts of the primary network on the performance of the secondary network with respect to the tolerant interference threshold at the primary receiver as well as the interference introduced by the primary transmitter at the secondary source and relay nodes. Simulation results are provided to validate the analysis.

1. Introduction

Energy harvesting (EH) has recently emerged as a promising approach to enable self-sustaining operation of energy-constrained wireless systems [1]. In EH-based communication systems, wireless devices scavenge energy from ambient energy sources, such as solar, wind, geothermal, and vibration sources [2]. In addition to these alternative energy sources, radio frequency (RF) signals have been considered as another promising energy source since it owns several favorable properties including availability and controllability. As RF signal is capable of carrying both information and energy at the same time, the authors in [3] devised receiver architectures for simultaneous wireless information and power transfer (SWIPT) systems. However, due to a huge gap regarding the operational sensitivity level between an information decoder and an energy harvester, SWIPT systems are only suitable for short-range transmission. To overcome this disadvantage, the authors in [4] proposed a novel wireless-powered system architecture in which power-beacons (PBs) are deployed to power wireless devices. Several PB-assisted energy harvesting systems have been then proposed and investigated in the literature; see, for example, [510].

Since EH-based point-to-point systems have been well investigated in the literature, many research works have recently focused on investigating energy harvesting in complicated system architectures, such as relay/cooperative systems, cognitive radios, and cognitive relaying networks. For relay systems, two EH relaying protocols, namely, time switching-based relaying (TSR) and power splitting-based relaying (PSR), are frequently used for energy harvesting and signal processing at a relay node [11]. In particular, in the TSR protocol, the relay node harvests energy from the received signal in a portion of time and spends the rest of the time period for information receiving and processing. Meanwhile, for the PSR protocol, the relay splits the received signal from the source into two streams, one of which is used for energy harvesting operation and the other for information decoding. Throughput analysis of amplify-and-forward (AF) and decode-and-forward (DF) relaying systems was performed in [11] and [9], [12], respectively. In addition, several aspects associated with EH relay systems have been examined, including EH-based relaying networks in the presence of interference [13], EH-based cooperative systems with spatially random relays [14], relay selection for trade-off between quality of information transfer and average transferred energy [15], and two-way EH-based relay systems [16].

With respect to cognitive radio (CR), it is well established that CR networks could offer a significant improvement in terms of spectrum utilization efficiency by allowing unlicensed secondary network (SN) users to share the same spectrum licensed to the primary network [17]. Among different cognitive paradigms, an underlay mode in which SN users transmit in the same band with the primary transmitters simultaneously subject to an interference constraint at the primary receiver is the simplest and is easy to implement in practice. Therefore, many studies have considered energy harvesting for underlay CR networks; see, for example, [1822]. In [18], the authors analyzed the outage probability and spatial throughput of the SN users. They also derived the optimal transmission power and density of SN nodes, which provides insights into optimal network design. In [19], a harvested energy-throughput trade-off optimization problem was formulated and a closed-form solution was obtained. Optimal resource allocation (i.e., time sharing and power) for EH-based CR networks was considered in [20, 21]. Also, an extension to joint time and power allocation for maximal sum-throughput in multiuser EH-based CR networks was studied in [22].

Recently, some research works have considered energy harvesting for cognitive relay networks (CRN). The motivation behind this approach is to exploit the advantages of cognitive radio and relay techniques [23]. In [24], the authors investigated EH-based underlay CRN with DF relaying and the TSR protocol. In their system model, secondary source and relay nodes harvest energy from multiple primary users that are closely located. They performed exact analysis of the outage probability as well as asymptotic analysis when the number of primary nodes goes to infinity. An analysis of outage probability in EH-based CRN with the PSR protocol was investigated in [25]. CRN with the PSR protocol was further investigated in [26] in which outage performance, energy-rate trade-off, and energy efficiency metrics were analyzed and optimized. In addition, analysis of the outage probability in multihop CRN and CRN with opportunistic relay selection was studied in [27] and [28], respectively. It is worth noting that these studies assume that the energy harvesting circuitry at receivers is always activated. Fading channel coefficients are assumed to follow independent and identical distribution (i.i.d). Moreover, instantaneous channel state information (CSI) of the links between SN nodes and a primary receiver (PR) is available. In many practical scenarios, it might be hard to access the instantaneous CSI since the collaboration between SN nodes and primary receiver is very limited.

Unlike the existing works, we consider in this study a distributed power-beacons (PBs) assisted energy harvesting cognitive relay wireless network. Some distinguished features of our system model are as follows: (i) simple PBs equipped with single-antenna are deployed to broadcast energy and they perform no additional information processing tasks; (ii) secondary source node and relay node need to harvest energy from PBs (and a primary transmitter if available) to sustain information transmission operations but they might be inactive if the harvested power is not sufficient, (iii) fading coefficients of the channels from PBs to secondary nodes follow i.n.i.d Nakagami- distributions; and (iv) instantaneous CSI of the links between SN nodes and PR might be available or not (i.e., only statistical CSI). Our aim is to comprehensively analyze the outage probability of the secondary relay system under a spectrum sharing context with energy harvesting capability. The main contributions of this work can be summarized as follows:(i)Derivation of analytical expressions for the power outage probability (i.e., the probability that the harvested power at the source node and/or relay node is below the minimum circuit activation power threshold) and the channel outage probability (i.e., the probability that the signal-to-noise ratio (SNR) of the received signal at the destination node is below the SNR threshold).(ii)Derivation of the total outage probability of the secondary network in both noise-limited scenario (i.e., the primary transmitter is far away from the secondary network) and interference-limited scenario (i.e., the primary transmitter is near to the secondary network).(iii)Derivation of the total outage probability of the secondary network subject to the availability of the instantaneous CSI or statistical CSI of the links between the secondary network and the primary receiver.(iv)Asymptotic analysis of the outage probability of the secondary relay system under different scenarios, which provides important insights into the behavior of the outage performance.(v)Evaluation of impacts of the transmit power of power-beacons, the tolerant interference threshold at the primary receiver, and the transmit power of the primary transmitter on the outage probability of the secondary relay system.

The rest of the paper is organized as follows. In Section 2, we describe the cognitive relay system model with RF energy harvesting. In Section 3, we analyze the power outage probability. Exact and asymptotic analysis of the total outage probability of the secondary relay system are carried out in Section 4. Simulation results and discussions are provided in Section 5. Finally, we conclude the paper in Section 6.

2. System Model

We consider a distributed PB-assisted energy harvesting cognitive relay wireless network as depicted in Figure 1. All nodes are equipped with single-antenna. The secondary network consists of a source node (SS), a decode-and-forward (DF) relay node (SR), and a destination node (SD). The source node and relay node are energy-constrained devices that need to harvest energy from PBs for data transmission. Also, the direct link between SS and SD is assumed not to exist due to poor fading conditions or obstacles. Hence, data transmission from SS to SD is accomplished with the help of SR. In this work, we consider an underlay cognitive paradigm, in which the secondary nodes can perform data transmission in the same spectrum band licensed to the primary network subject to interference constraint at the primary receiver.

We consider one frame of duration , where the fading channels are assumed unchanged. In the first phase of duration , SS and SR harvest energy from PBs. In the second phase , SS transmits its data to SR, and in the third phase , SR forwards the information to SD. For convenience, we denote parameters of the links between nodes in the network as in Table 1. Due to short distance between PBs and SS and SR, the line-of-sight (LoS) path is likely to exist. To this end, we assume that and , follow i.n.i.d Nakagami-m distribution with fading parameters and and squared means and , respectively, where denotes an expectation operation. For the other links, we assume that the fading coefficients follow i.n.i.d Rayleigh distribution.

2.1. Primary Transmitter Is Far Away from the Secondary Network

We first consider a scenario in which the primary transmitter (PT) is located far away from the secondary network (SN). Therefore, PT does not introduce interference to both SR and SD (i.e., a noise-limited scenario). Also, the energy harvested at SS and SD associated with the RF signal radiated by the primary transmitter is negligible.

The total energy harvested at the SS node at the end of the energy harvesting phase can be expressed aswhere is the efficiency of an energy harvester, is the transmit power of PB, is the distance between the th PB and the source node SS, is the path-loss at reference distance of 1 m, and is the path-loss exponent. Also, and . We note that the probability distribution function (PDF) of is given aswhere . Thus, the PDF of can be expressed by [29]where is a weighting factor defined similar to that in [29, Eq. ()]. In particular, the recursive formula for evaluating can be expressed asSimilarly, the energy harvested at the SR node is given aswhere and the PDF of is expressed aswhere and is computed similar to .

For an underlay cognitive paradigm, the spectrum licensed to the primary network is utilized by the secondary network. However, the transmit powers of both SS and SR nodes need to be adjusted subject to the interference constraint at the primary receiver. We note that the transmit powers of SS and SR are governed by the availability of the channel state information (CSI) of the SS-PR and SR-PR links at SS and SR nodes, respectively. Thus, in this work, we consider two cases depending on whether the instantaneous CSI is available or not.

2.1.1. Instantaneous CSI of SS-PR and SR-PR Links Is Available

When the instantaneous CSI of the SS-PR and SR-PR links is available, the SS and SR can adjust their transmit powers to guarantee that the interference at PR is always below the maximum acceptable value . Therefore, the transmit powers at SS and SR can be computed, respectively, aswhere , , and . Given that the primary transmitter is located far away from the secondary network, we can express the SNR at SR and SD asrespectively, where , , and is the noise power of the additive white Gaussian noise (AWGN) at SR and SD. Finally, the end-to-end SNR at the SD node in the DF relay system is obtained as

2.1.2. Statistical CSI of SS-PR and SR-PR Links Is Available

Let us now consider the case when only statistical CSI of the SS-PR and SR-PR links is available. We note that a focus on the statistical CSI case is useful from a practical viewpoint. This is due to the fact that the collaboration between secondary nodes and PR is limited. Thus, it is hard to obtain perfect instantaneous CSI of SS-PR and SR-PR links at SS and SR, respectively. Under this condition, the interference constraint at PR due to SS transmission can be formulated as [30]Where Pr stands for probability, is the transmit power of SS, and is the tolerated error. Given that follows an exponent distribution with mean , the transmit power by SS needs to satisfySimilarly, the transmit power of SR is constrained bywhere .

Thus, the actual transmit powers of SS and SR are obtained, respectively, asAlso, the end-to-end SNR is given aswhere

2.2. Primary Transmitter Is Near to the Secondary Network

We next consider a scenario where the primary transmitter is near to the secondary network. In this case, SS and SR can further harvest energy from the RF signal radiated by PT. However, it is worth noting that PT will introduce interference to both SR and SD during the information transmission period. For analytical simplicity, we assume that the interference dominates the noises at SR and SD nodes (i.e., an interference-limited scenario).

Let denote the transmit power of the primary transmitter. In this scenario, the energy harvested by SS and SR during the energy harvesting phase consists of the energy broadcasted by PBs and the energy radiated by PT. Similar to the previous case, we can express the harvested energy at SS aswhere . It is worth noting that follows the Rayleigh distribution which can be interpreted as a Nakagami- variable with a parameter . Therefore, is a sum of Nakagami- variables and its PDF can be evaluated in a similar manner as (cf. (6)). Also, the harvested energy at SR is given aswhere . Similar to the previous scenario, the end-to-end SNR when PT is near to SN is obtained aswhere, in the presence of the instantaneous CSI,where and . Also, when only statistical CSI is available, we have

3. Analysis of Power Outage Probability

In RF-powered wireless communication systems, the received power at EH-enabled devices must be larger than the minimum power threshold to activate the energy harvesting circuitry. A typical value of is in the range from −30 dBm to −10 dBm in practice [1]. In the considered system, if the received power is below , the source node and relay node are inactive and thus no information transmission can occur. Thus, in this section, we evaluate the power outage probability of the system, that is, the probability that the harvested power is smaller than .

Let us first consider the scenario where PT is far away from SN. The received power at the source node and relay node can be expressed as , where (cf. (1) and (5)). Thus, the power outage probability due to insufficient harvested power at SS is calculated aswhere is the cumulative distribution function (CDF) of that can be obtained based on (2) and (3). Consequently, (25) can be further expressed asThe power outage due to the insufficient harvested power in the system happens when SS or SR or both of them experience the power outage. Consequently, the power outage probability of the SN network can be computed as

In case that PT is near to SN, we obtain a similar expression for the power outage probability; that is,whereNumerical results based on these expressions are provided in Section 5.

4. Analysis of Total Outage Probability

In the considered system, the information outage at the destination node can occur due to (i) the power outage at SS and/or SR node since the harvested power is not sufficient; (ii) SS and SR being active but the channel conditions of SS-SR and SR-SD links are poor so that the SNR of the received signal at the SD node is below the SNR threshold (hereafter we refer to this condition as the channel outage). Therefore, the total outage probability of the secondary network can be computed aswhere and are the power outage probability and the channel outage probability, respectively. In what follows, we derive expressions for the channel outage probability in two scenarios subjecting to the location of PT relative to the SN network.

4.1. Primary Transmitter Is Far Away from the Secondary Network

The channel outage probability of the secondary DF relay system can be expressed as [9]where is the SNR threshold and and are the CDF of and , respectively.

4.1.1. Instantaneous CSI

Let us first consider the case in which the instantaneous CSI of the SS-PR and SR-PR links is available. We have the following results.

Proposition 1. The exact CDF expression of when the instantaneous CSI of SS-PR and SR-PR links is available is shown in (32), where , and is the th order modified Bessel function of the second kind [31, Eq. ()].

Proof. See Appendix A.

Based on the result in Proposition 1, we can easily deduce the CDF of by replacing the parameters in by their counterparts in (i.e., , , , , and ). Then, by substituting the CDF expressions of and into (31), we obtain the closed-form expression of the channel outage probability.

Asymptotic Analysis. In the high SNR regime (i.e., or, equivalently, ) and assuming that is fixed and independent of , and in (8) and (9) can be approximated toBy evaluating and , we obtain the following result.

Corollary 2. The asymptotic CDFs of and at the high SNR regime when is independent of are, respectively, given as

Proof. The proof is straightforward and is thus omitted to save space.

In cases that the interference threshold is proportional to transmit power , that is, , where is a positive constant, we obtain the following asymptotic result.

Corollary 3. The asymptotic expression for the CDF of at the high SNR when is proportional to is given asAlso, the asymptotic CDF expression of is obtained by exchanging the parameters in by their counterparts in .

Proof. See Appendix B.

It is worth noting that the condition of implies that the impact of the interference is nonnegligible regardless of the transmit power of secondary nodes [25].

4.1.2. Statistical CSI

For the case that only statistical CSI is available, we obtain the exact and asymptotic expressions as shown below.

Proposition 4. The exact CDF expression of is shown in (36), where and is the upper incomplete Gamma function [31, Eq. ()]:

Proof. See Appendix C.

Corollary 5. The asymptotic expressions for CDFs of and at the high SNR regime when is independent of are obtained as

Proof. When and is fixed, the above result is obtained based on the exact CDF expression in (36) by using the approximation of , when , and the property of .

Corollary 6. The asymptotic expression for CDF of when is proportional to is obtained asThe asymptotic CDF expression of is obtained by replacing the parameters in by their counterparts in .

Proof. The result is obtained by approximating the exact CDF expression in (36). Please see Appendix D for details.

4.2. Primary Transmitter Is Near to the Secondary Network

The channel outage probability of the system in this scenario can be expressed aswhere and are the CDFs of and , respectively.

4.2.1. Instantaneous CSI

In scenarios that instantaneous CSI is available, we obtain the following results.

Proposition 7. The exact CDF expression of is given in (40), where , , , , , and denotes the Whittaker function [31, Eq. ()].

Proof. See Appendix E.

Corollary 8. The asymptotic CDFs of and at the high SNR regime when is independent of are, respectively, given aswhere denotes the exponential integral function [31, Eq. ()].

Proof. See Appendix F.

Corollary 9. The asymptotic CDF expression of when is proportional to is obtained aswhere is the Psi function [31, Eq. ()].

Proof. See Appendix G.

4.2.2. Statistical CSI

When only statistical CSI is available, we obtain the exact and asymptotic CDF expressions as below.

Proposition 10. The exact CDF of is shown in (44), where ,

Proof. As a sketch of the proof, we first calculate the unconditional CDF of (cf. (23)). Then, the unconditional CDF expression is obtained by taking the expectation of with respect to .

Corollary 11. The asymptotic CDF expressions of and when is fixed and independent of are, respectively, given as

Proof. The result is obtained based on the exact CDF expression in (44) by using the approximated formula of , when , and the property of .

Corollary 12. The asymptotic CDF expression of when is proportional to is obtained as

Proof. See Appendix H.

Remarks 1. From the derived asymptotic expressions above, we have some remarks in the high SNR regime as follows: (1)From Corollaries 2, 5, 8, and 11, it is observed that when the interference threshold is fixed and independent of the transmit power , an outage probability floor occurs in the high SNR regime. The larger the value is, the lower the floor appears. In addition, this outage probability behavior also implies that no diversity gain is achieved in the high SNR region under these scenarios.(2)In scenarios that PT is near to the secondary network, a smaller transmit power of the primary transmitter results in lower outage probability. When only statistical CSI is available, a higher tolerated error also leads to lower outage probability (see, e.g., Corollary 11).(3)When is proportional to , it can be concluded from Corollaries 6 and 12 that a diversity order equal to one is attained.Numerical results based on the derived expressions will be presented in Section 5.

5. Results and Discussions

In this section, we provide simulation results to validate the derived expressions in the previous section. Also, we evaluate impacts of the transmit power of power-beacons, transmit power of the primary transmitter, and the interference threshold at the primary receiver on the outage probability of the secondary relay system. We examine two scenarios, namely, Scenario ( PBs) and Scenario ( PBs). Without loss of generality, we assume that the coordinates of the SS, SR, SD, and PR nodes in the -plane are (), (), (), and (), respectively. For Scenario , the PBs are located at (), (), and (). In Scenario , the locations of the two additional PBs are at () and (). In case that PT is located near the secondary network, we assume that the coordinate of PT is (). We consider , , for simple presentation given that listing out all the values is lengthy. Other simulation parameters are , , , , , , = 0 dB,  W, , , , and .

In Figure 2, we plot the power outage probability versus the transmit power of power-beacons. The power threshold for EH circuitry activation is set to −30 dBm [1]. It can be seen that when the transmit power of PBs is 20 dBm, the probability that the outage occurs due to insufficient harvested energy is very low. Moreover, when the number of PBs is increased (i.e., in Scenario ), the power outage is reduced significantly since a larger amount of energy is harvested. Another observation is that when the primary transmitter (PT) is near to SN, the outage probability is smaller as additional energy is harvested from PT. However, this benefit becomes less significant when the number of PBs is increased.

In Figures 3 and 4, we plot the channel outage probability versus when the instantaneous CSI is available. When there is only the statistical CSI, the results are shown in Figures 5 and 6. It can be seen from all the figures that the curves based on analysis match exactly with the simulation curves in all cases, which confirms the accuracy of the derived expressions. Moreover, as expected, the asymptotic curves are close to the exact curves in the high SNR regime. We have some other important observations from these results as follows:(1)In Scenario , the outage probability is lower than that in Scenario since higher energy is harvested due to a larger number of employed PBs.(2)When PT is located close to the secondary network, the outage probability is higher compared to its counterparts. This is because, under this condition, the benefit of increased energy harvesting at SS and SR nodes thanks to the signal radiated by PT is not significant compared to the disadvantages due to the interference that this signal introduces at SR and SD nodes.(3)When is fixed and independent of , the outage probability floor appears at the high SNR region. This is due to the fact that, for a fixed value , increasing does not necessarily lead to higher transmit power of SS and SR nodes since the transmit power is now determined by (cf. (7), (14)-(15)). Note that this behavior of the outage probability suggests that cognitive relay systems are not applicable in the high SNR regime.(4)When is proportional to , a full diversity order can be attained. A similar conclusion was made in [25], where a different EH-based CRN system was investigated.

To further investigate impacts of the primary network on the outage probability of the secondary relay network, we plot in Figures 7 and 8 the total outage probability versus the tolerant interference threshold and the transmit power , respectively. It can be observed from Figure 7 that the outage probability is decreased when the interference threshold is higher. This is because larger implies that higher transmit power of SS and SR can be allowed, which in turns leads to lower outage probability. However, the outage probability exhibits an error floor when is sufficiently large. This is due to the fact that, in this region, the maximum transmit power of SS and SR depends on the harvested energy only since the interference constraints are relaxed (cf. (7), (14)-(15)). This also explains why in Figure 7 the error floors in Scenario 2 are lower than those in Scenario 1. In Figure 8, we notice that higher results in higher outage probability in all scenarios. This can be explained similar to Comment 2 regarding the location of PT relative to SN nodes presented in the previous paragraph.

6. Conclusions

In this paper, we have studied distributed power-beacon-assisted energy harvesting cognitive relay wireless networks. We have considered practical scenarios that (i) secondary nodes might be inactive due to insufficient harvested energy; (ii) the fading channels from PBs to secondary nodes undergo i.n.i.d Nakagami- fading distribution, and (iii) the instantaneous CSI of the links between SN nodes and PR might not be available. The exact expressions for the outage probability of the secondary network under difference scenarios have been derived. Asymptotic analysis for the outage probability has also been performed, which provides insights into the system behavior. We have evaluated the impacts of the transmit power of PBs, the interference threshold at the primary receiver, and the transmit power of the primary transmitter on the outage probability. The results show that when is fixed and independent of , the outage probability floor exhibits at the high SNR region. On the other hand, if is proportional to , a full diversity order could be achieved. Our future work would be a consideration of a general system model where nodes are equipped with multiple antennas.

Appendix

A. Proof of Proposition 1

The CDF of can be calculated as (cf. (8))Note that and are exponentially distributed random variables with parameters of and , respectively (i.e., , , and , ). Also, the CDF of can be obtained from the PDF function in (3) asBy substituting these results into (A.1) and applying the integral results in [31, Eq. ()] and [31, Eq. ()], we obtain (32) after some mathematical manipulations.

B. Proof of Corollary 3

When and , we note that and is a constant. Thus, the second and the third terms within the summation in (32) cancel out with each other and are similar to the fourth and the fifth terms. Consequently, the exact CDF expression in (32) can be approximated toOn the other hand, it can be proved from (2) and (3) thatTherefore, by using the series representation of the Bessel function provided in [31, Eq. ()], we can approximate (B.1) to (35) after some manipulations.

C. Proof of Proposition 4

From (17), we can evaluate the CDF of asBy substituting the expressions of , , , and into (C.1) and then using the results in [31, Eq. ()], [31, Eq. ()], and [31, Eq. ()], we finally obtain (36).

D. Proof of Corollary 6

By adopting the approximation of , when , we can approximate the integral defined in Proposition 4 when aswhere is the lower incomplete Gamma function [31, Eq. ()]. Recall that and . Thus, we have that is a constant. Thus, under the condition that and , the exact CDF expression in (36) can be approximated toBy using the equation of and the property of , we can simplify (D.2) to (38) after some manipulations. This completes the proof.

E. Proof of Proposition 7

Let us begin by calculating the conditional CDF of with respect to , that is,We note that there is the similarity in terms of mathematical description between (A.1) and (E.1) (i.e., versus ). Thus, similar to the proof in Appendix A, we can express the closed-form expression for as in (E.2).The unconditional CDF of is then obtained by taking the expectation of with respect to . In particular, we havewhere is the PDF of , that is, (). By applying the integral result of , where denotes the Whittaker function [31, Eq. ()] and is the Gamma function [31, Eq. ()], the integral in (E.3) can be solved. The result in (40) is thus obtained.

F. Proof of Corollary 8

In the high SNR regime (i.e., ), defined in (14) can be approximated toThe CDF of conditioned on can be computed asThe unconditional CDF of is then obtained asBy using the integral result in [31, Eq. ()], we can solve the integral in (F.3), which finally yields (41). In addition, the CDF expression of is obtained by exchanging the parameters in by its counterparts in . This completes the proof.

G. Proof of Corollary 9

We use a similar approach as done in Appendix E. At first, we express the approximated conditional CDF expression of with respect to as (cf. (35))The unconditional CDF expression is then obtained asWith the help of [31, Eq. ()] and [31, Eq. ()], we can calculate the resulting integrals in (G.2). The result in (43) is thus proved.

H. Proof of Corollary 12

Let us start by noting that the Taylor series expansion of at can be expressed as . Therefore, we can approximate the integral defined in Proposition 10 when asWhen and , the exact CDF expression in (44) can be approximated toBy using the equation of , the property of , and the Laurent series expansion of at , that is, , we can simplify (H.2) to (46). This completes the proof.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.