Mixed FSO/RF Based Multiple HAPs Assisted Multiuser Multiantenna Terrestrial Communication

Abstract

In this work, a mixed free-space optics (FSO)/radio-frequency (RF) based multiple serial high altitude platforms (HAPs) assisted multiuser multiantenna terrestrial communication system is considered. For the considered multi-hop system, earth station to HAP and HAP to HAP links are assumed as FSO links, and HAP to terrestrial mobile users (MUs) link is assumed as RF link. At the FSO detector both the heterodyne detection and intensity modulation direct detection techniques are considered. Atmospheric turbulence of the ES to HAP FSO link is modeled with Gamma-Gamma fading along with the pointing error impairments. As HAP situates in the stratosphere, negligible atmospheric turbulence exists at HAP altitude (20 km), hence, only pointing error is considered for the inter HAP links. Further, the HAP consists of multiantenna array to provide high data-rates to the terrestrial users via RF links, considered to be Nakagami-m distributed. The intermediate HAPs perform selective decode-and-forward relaying and opportunistic user scheduling is performed for the best user selection. For the performance analysis, analytical expression of overall outage probability is obtained and the impact of various selection parameters like pointing error, FSO detection type, MU selection, number of antennas, and RF fading severity are observed on the outage performance. Asymptotic outage probability is also derived to obtain the diversity order of the considered communication system. Further, considering various modulation schemes, a generalized average bit-error-rate expression is derived and the impact of pointing error, FSO detection type, MU selection, number of antennas, and RF fading severity are observed on its performance. Finally, derived results are validated through Monte-Carlo simulations.

1 Introduction

For reliable and robust high data-rate communication for various wireless applications including the emergency disastrous situations, temporary mass events, and even in military applications, promptly deployable flying base stations such as unmanned aerial vehicles (UAVs) have gained significant research attention as they are capable to provide high data-rates with improved capacity and coverage, targeted for beyond fifth generation (5G) communications(Zhan et al., 2011; Mozaffari et al., 2019). UAVs are broadly categorized into low altitude platform (LAP) and high altitude platform (HAP) based on their size and operating range. The LAPs are situated only few hundred meters above the earth surface, however, the HAPs are usually situated in stratosphere around 17–25 km (km) altitude from the earth surface. As the HAP situates in stratosphere, the wind velocity is low and atmospheric turbulence is negligible, hence, the HAP can be in quasi-static position and is capable to provide ubiquitous connectivity (Kurt et al., 2021). Due to the enhancement in communications technologies, lightweight composite materials, improved efficiency of solar panels and antenna, HAP became an viable aerial network component for various communication applications. The HAPs can be deployed easily and rapidly, their maintenance is easy, and having moderate operational cost than satellite which makes them economically feasible (Fidler et al., 2010; Kurt et al., 2021). Various technological trends of HAPs and their applications in wireless communications have been seen in (d’Oliveira et al., 2016). On the other hand, free-space optics (FSO) has been proven as an efficient replacement to radio-frequency (RF) due to its high operational bandwidth with improved capacity to fulfill the increasing demand of high data-rates for present and future wireless communication networks. FSO operates in unlicensed band, having low cost and easy deployment, and provides high-speed line-of-sight (LoS) communication (Singya et al., 2020; Trichili et al., 2020). The FSO links provide high data-rate secured transmission than the RF links, however requires LoS communication. Further, atmospheric turbulence and weather issues are more severe in terrestrial mobile communication which limits the FSO link’s performance. In such situations, RF links provide more robust communication to mobile users (MUs) than the FSO links. Thus, a hybrid system which can utilize the high data-rate capability of FSO links and the reliable transmission through RF links to terrestrial MUs has gained significant research attention (Trichili et al., 2021).

Various HAP networks like HAPCOS and CAPANINA in Europe, Helicos in USA, SkyNet is Asia, etc. have been deployed successfully for broadband data transmission to a larger distance (in hundreds of km). A considerable work on HAP networks have been seen in the literature. In (Vaiopoulos et al., 2013), authors proposed a HAP assisted setup for WiMAX orthogonal frequency division multiplexing (OFDM) transmission and derived the outage probability. In (Fidler et al., 2010), authors have shown various field trials to and from HAP using optical links to provide high data-rates. In (David et al., 2004), various considerations for inter-HAP optical link design are shown. In (Parthasarathy et al., 2014), channel modeling of inter-HAP optical link is proposed. In (Shibata et al., 2020), a gigabit HAP system for mobile communications is proposed. An HAP assisted ground to satellite uplink FSO/RF system is proposed in (Swaminathan et al., 2021) with Gamma-Gamma and shadowed-Rician fading for FSO and RF links, respectively, and various performance measures are derived. In (Singya and Alouini, 2021), an HAP assisted downlink hybrid FSO/RF based multiuser multiantenna terrestrial communication system is proposed, where the FSO link is characterized with the Gamma-Gamma distribution with pointing error impairments and Nakagami-m distribution is preferred for RF links. For the performance analysis, analytical expressions of outage probability, asymptotic outage probability, ergodic capacity, effective capacity, and generalized average symbol-error-rate (ASER) expressions of hexagonal-quadrature amplitude modulation (QAM), cross-QAM, and rectangular QAM are derived. (Liu et al., 2020) focuses on amplify-and-forward (AF) based HAP assisted satellite-terrestrial communication with Gamma-Gamma and Rayleigh distributions for the FSO and RF links, respectively, and obtained the outage probability. In (Antonini et al., 2006), capacity and efficiency of FSO link between the satellite and HAP are obtained for its feasibility. In (Gao et al., 2021), optimization of an intelligent reconfigurable surface assisted HAP network is proposed for downlink communication to terrestrial MUs.

However, a significant work is witnessed on various hybrid FSO/RF systems in the literature. (Bhatnagar and Arti, 2013) considers a hybrid RF/FSO based satellite-terrestrial communication system, where Gamma-Gamma and shadowed-Rician distributions are considered for the FSO and RF links, respectively. Further, analytical ASER expression for M-ary phase-shift keying (MPSK) and asymptotic results are obtained. (Zedini et al., 2016) focuses on a hybrid FSO/RF system with Nakagami-m and Gamma-Gamma fading for the RF and FSO links, respectively, and the analytical expressions of ergodic capacity, outage probability, and bit-error-rate (BER) of binary phase-shift keying (BPSK) are derived. (Chen et al., 2019) covers energy harvesting in hybrid FSO/RF system by obtaining the outage probability. An FSO/RF based satellite-terrestrial communication system is considered in (Ahmad et al., 2017), where satellite is connected through the optical feeder link and provides RF connectivity to terrestrial users, and various performance measures are derived for analysis. (Zedini et al., 2020) focuses on an FSO/RF based satellite-terrestrial communication system with Gamma-Gamma distributed uplink optical feeder link and shadowed-Rician distributed downlink RF terrestrial users links, and ergodic capacity, outage probability, and BER expressions of BPSK are derived. In (Sharma et al., 2019), switching-based FSO/RF system is proposed and various performance measures are obtained. (Lee et al., 2020) focuses on the data throughput maximization for UAV based FSO/RF system. Further, in (Singya et al., 2020), ASER performance of various higher order complex QAM schemes is obtained for a mixed RF/FSO system with outdated channel state information. In (Xu and Song, 2021) and (Xu and Song, 2020), hybrid RF/FSO based communication systems are considered, where RF link is characterized with κ − μ distribution and the FSO link is characterized with the - distribution with pointing error impairments. For performance analysis, outage probability, asymptotic outage probability, ergodic capacity, and ABER results are obtained. In (Xu and Zhang, 2021), a hybrid RF/FSO based deep space communication system is considered, where exponentiated Weibull and Nakagami-m distributions are considered respectively for the FSO and RF links; and outage probability, ergodic capacity, and ABER of MPSK are derived.

In the literature, some papers have discussed FSO based multi-hop communication systems. In (Datsikas et al., 2010; Zedini and Alouini, 2015; Ashrafzadeh et al., 2020), authors have analyzed the performance of multi-hop FSO systems. (Altubaishi and Alhamawi, 2019) focuses on the capacity analysis of a multi-hop FSO/RF system. In (Wang et al., 2015), authors have discussed the performance of multi-hop exponentiated Weibull distributed FSO system. However, above discussed multi-hop works focus only on generalized FSO systems without including HAP networks.

On the other hand, HAPs are having a wide 500 km radius footprint as recommended by International Telecommunication Union (ITU) (ITU-Radiocommun, 2000). Further, most of the HAP assisted projects have smaller coverage area and combining multiple HAPs can provide coverage to even entire country. For example, authors in (Miura and Oodo, 2001) and (Milas et al., 2003) have proposed respectively 16 and 18 HAPs configurations to provide coverage to entire Japan and Greece. However, only few works focus on the performance of multi-hop HAP/LAP assisted communication system. In (Yang et al., 2018), outage probability of LAP assisted RF/FSO/RF multi-hop system is derived. In (Sharma et al., 2016), authors have evaluated the performance of inter-HAP FSO communication system. Authors in (Michailidis et al., 2018) have discussed the performance of a three-hop HAP assisted RF/FSO/RF communication system. The above discussed works are summarized in Table 1 in details. The various considerations, channel models, and derived performance metrics are shown in Table 1 which are compared with our work (highlighted in bold letters) to show the novelty of our work.

Note that the proposed system model and the analysis shown in this work is generalized. By considering number of HAPs N = 1, the multi-hop serial HAPs assisted mixed FSO/RF system can be converted into a dual-hop HAP assisted mixed FSO/RF system. By considering number of users U = 1, the performance of Uth terrestrial MU can be obtained. Further, all the parameters can be adjusted and the considered system model can be deployed for various specific applications.

Rest of the work is organized as follows: In Section II, proposed system model and channel model are discussed in detail. Various performance measures like outage probability, asymptotic outage probability, and generalized ABER are discussed in detail in Section III, Section IV, and Section V, respectively. Section VI discusses various theoretical and simulation results of the proposed system and conclusions from the obtained results are finally illustrated in Section VII.

2 System and Channel Models

2.1 System Model

In this work, a mixed FSO/RF based multiuser multiantenna inter-HAP configuration is considered to provide high data-rates to the terrestrial MUs as shown in Figure 1. We have considered a communication scenario, where a source ES is communicating to the terrestrial MUs which are far apart through the use of multiple serial HAPs. Initially, the ES transmits the optical signal s(t) to HAP in first time slot. Hence, the received optical signal at the first HAP can be expressed aswhere η is the optical-to-electrical conversion coefficient and r decides the type of FSO detector, wherein r = 1 corresponds to the heterodyne detection and r = 2 represents the IM/DD. Here I = I_aI_p represents the fading coefficient of the FSO link, where I_a is the fading due to atmospheric turbulence and I_p represents the pointing error associated with the FSO link. Further, represents the additive white Gaussian noise (AWGN) with zero mean and variance associated with the ES to HAP link. The noise variance is , where K = 1.38 × 10^–23 is the Boltzmann constant, T_o represents the noise temperature in kelvin, and B_o is the optical receiver bandwidth. Further, δ_EH represents the path-loss while transmitting the source signal from ES to HAP through the FSO link. , where G_tE is the ES transmitter telescope gain, G_rH is the receiver telescope gain at the HAP, and is the free-space loss, wherein L is the ES-HAP path length and λ_o is the optical signal wavelength. A 1550 nm wavelength optical beam is considered which is eye safe and most suitable due to minimum Doppler shift than other frequencies. Therefore, the instantaneous received SNR at HAP will be

FIGURE 1

N intermediate HAPs with DF relaying are considered for e2e signal transmission. Hence, at the HAP, the received signal is first decoded, re-encoded, and then transmitted to the next HAP. The received signal can only be decoded successfully if the instantaneous received SNR at nth HAP is greater than or equal to the threshold SNR i.e. . Hence, probability of successful detection is defined as

Let, is the decoded signal, then the optical signal received at the nth HAP is represented aswhere is the transmitted power from the (n − 1)^th HAP, I_p is the FSO channel fading coefficient which only consists of pointing error and no atmospheric turbulence is considered. Further, n_HH(t) represents the AWGN associated with the HAP-HAP link with zero mean and variance. Also, is the path-loss of the HAP-HAP link and considered to be same between all the HAPs. Here, G_tH represents the gain of the transmitting telescope of the (n − 1)^th HAP and G_rH is the receiver telescope gain of the nth HAP. Further, is the free-space loss in the HAP-HAP link, where D is the inter-HAP distance. Hence, instantaneous received SNR at the nth HAP is given as .

Finally, the decoded signal from the Nth HAP is transmitted to the terrestrial MUs. To improve the channel capacity of the RF link and to improve its reliability, N_t transmit antennas are deployed at the HAP. However, due to size limit constraint, only single antenna is assumed at the terrestrial MUs. Hence, in the last-hop, the signal received from the Nth HAP at the uth MU after transmit beamforming is given aswhere P_N is the Nth HAP’s transmit power and is the probability of successful detection till the Nth HAP. Here and denote the channel vector of N_t × 1 order and transmit beamforming weight vector between the Nth HAP and the uth UE, respectively. Further, (⋅)^H is a Hermitian operator and maximal ratio transmission (MRT) principle is used for the beamforming weight vector as , wherein (⋅)‖_F being the Frobenius norm. Further, represents the AWGN associated with the RF link of zero mean and variance, where K = 1.38 × 10^–23 is the Boltzmann constant, T_r represents the noise temperature in kelvin, B_r is the RF receiver bandwidth, and N_F is the noise figure. Furthermore, is the path-loss associated with the RF link which can be given as (Swaminathan et al., 2021)where G_tH is the Nth HAP’s transmitter gain, G_ru is the uth user’s receiver gain, L_F = 20 log₁₀f_RF + 20 log₁₀L_u + 92.45 is the free-space loss of the RF link between the HAP and uth terrestrial user of L_u distance (in km) with f_RF transmission frequency, L_A is the gaseous atmospheric loss, L_R is the loss by rain attenuation in dB/km, and L_O represents the miscellaneous losses including antenna degradation, polarization mismatch, antenna mispointing, etc. Hence, e2e instantaneous received SNR at the uth MU is calculated as

2.2 Channel Model

2.2.1 Earth Station to HAP Link

ES-HAP link is considered as an FSO link and its fading coefficient is characterized as I = I_pI_a, where I_p is the pointing error and I_a is the atmospheric turbulence induced fading. Let us assume, a Gaussian beam of beamwidth ω_L is propagating L distance from the transmitter to receiver with a aperture radius, the fraction of the transmitted power collected by the receiver is approximately given as (Farid and Hranilovic, 2007)where A₀ represents the fraction of the collected power at r_d = 0, r_d is the radial displacement between the detector and beam center, and is the equivalent beamwidth which is given as . Further, , where , wherein erf(⋅) represents the error function. It is to be noted that the approximation (8) exists only for ω_L > 6a. Figure 2 shows various pointing errors for an FSO system. Let us consider d_x and d_y are the horizontal and vertical displacements of the beam in the detector plane, then denotes the vector for radial displacement as shown in Figure 2. It is assumed that and are independent and Gaussian distributed as and , respectively. Then the amplitude of the radial vector is Beckmann distributed (Beckmann and Spizzichino, 1987). Thus, the pointing error of the FSO link can be modeled in various ways according to the jitter and boresight values (Jung et al., 2020).

FIGURE 2

Considering the case of identical jitters and zero boresight error, mean values ψ_x and ψ_y of the radial vectors d_x and d_y are zero and variances are identical as . Hence, the radial displacement follows Rayleigh distribution. For such case, the pointing error can be modeled aswhere is the ratio of the equivalent beam radius to the standard deviation of jitter at the receiver.

The atmospheric turbulence I_a can be modeled with Gamma-Gamma distribution, its probability density function (PDF) is given aswhere K_v(⋅) represents the vth order modified Bessel function of second kind. Further, α and β are the fading coefficients of the atmospheric turbulence which are given aswhere is the Rytov variance corresponds to the turbulence strength metric. A slant path exists between the ES and the HAP, hence, the Rytov variance for this path can be given aswhere θ is the zenith angle, k₀ = 2π/λ_o denotes the wave number, H is the HAP altitude, and h₀ is the ES height. Further, is the refractive index structure parameter, depends upon the strength of the atmospheric turbulence and is modeled by Hufnagel-Valley model as a function of altitude h which is given aswhere shows the turbulence strength at ground in m^−2/3 and w₀ is the rms wind speed in m/s. The PDF of the FSO channel fading coefficient (I) can be given as

Considering the identical jitters and zero boresight error case, the resulting conditional distribution is expressed as

Substituting (10) and (15) in (14), we get

Representing the modified Bessel function of the second kind into Meijer-G form as , solving (16) with the help of (Wolframe, 1998, (07.34.21.0085.01)), and after some mathematical simplifications, we get

Instantaneous received SNR of the ES-HAP FSO link is calculated as γ_H = μ_r(I)^r, where r represents the type of FSO detector and μ_r is the average power of the FSO link. For heterodyne detection, r = 1, and hence, . For the IM/DD detection, r = 2, and hence, . After change of variable in (17), the PDF of the e2e SNR of the FSO link is given as

Substituting (18) in and applying (Wolframe, 1998, (07.34.21.0084.01)) with some mathematical computations, the cumulative distribution function (CDF) of the SNR of the FSO link for both the detection techniques is given aswhere , , , and .

2.2.2 Inter HAP Links

The HAP is situated in the stratosphere above the clouds and there is almost negligible atmospheric turbulence. Hence, only the pointing error is considered between inter-HAP links. As identical jitters with zero boresight error is considered, the pointing error can be modeled by Rayleigh distribution and its PDF is given in (9). Now instantaneous received SNR at the nth HAP for the FSO link is calculated as , where is the average received power at the FSO detector of nth HAP. After change of variable in (17), the PDF of the instantaneous received SNR of the HAP-HAP link is given as¹

Performing , the CDF of the HAP-HAP link can be given as

2.2.3 HAP to MUs RF Links

To connect the terrestrial MUs, communication between the HAP and the MUs is established through the RF links. Considering the size constraint of the MUs, single antenna is assumed at the MUs. Further, multiple transmit antennas (N_t) are considered at the HAP to avail the antenna diversity. The Nakagami-m distribution is considered for the RF link i.e. Nak(m, Ω), where m is the fading severity and Ω is the average power of the RF link. Hence, the CDF and PDF of the uth MU’s SNR can be given asrespectively, where Γ(⋅) represents the complete gamma function and . It is considered that the terrestrial MUs are independent and identically distributed (i.i.d.) as they are close to each others. Hence, opportunistic user scheduling is performed to achive the multiuser diversity. Therefore, the instantaneous received SNR of the HAP-MU link after opportunistic scheduling is . Thus, the PDF of γ_HM after the order statistics is given as

Invoking the PDF and CDF values from (22) in (23) and using the series representation of the incomplete Gamma function with multinomial expansion using (Gradshteyn and Ryzhik, 2000, (0.314)), the PDF of γ_HM can be derived aswhere is recursively calculated as , , for 2 ≤ l ≤ (mN_t − 1), for mN_t ≤ l < j(mN_t − 1), and , wherein (Singya et al., 2018, 2019). Using , the CDF of (24) can be obtained aswhere , , , and .

3 Outage Probability

Outage probability is one of the important performance measures and defines the probability of reaching the instantaneous e2e SNR below a fixed threshold (γ_th). For the considered system, e2e instantaneous received SNR is shown in (7), hence the overall outage probability is given as

Substituting the respective CDF values of the individual links in (26), we get

Using the binomial series expansion and after some mathematical computations, (27) can be solved as

4 Asymptotic Outage Probability

Various insights on the system’s performance can be highlighted through the obtained outage probability (28). However, (28) is quite complex to decide the diversity order of the considered system. Hence, in this Section, we have conducted asymptotic analysis on the outage probability to obtain the diversity order. For this, the transmit SNR is assumed to tend to infinity. Therefore, at high SNR, (26) can be approximated as

5 Average BER Analysis

For various modulation schemes, a generalized ABER expression for the considered communication system can be given as (Zedini et al., 2020)where n₀, κ, and g_p are the selection parameters for various modulation schemes and their values are shown in Table 2.

Please note that by using the IM/DD, we can obtain the ABER performance of on-off keying (OOK), preferred commonly in practical FSO systems as it is flexible to laser nonlinearity with ease of use. Further, we have analyzed the ABER performance of various MPSK (including the BPSK) and MQAM schemes using the coherent heterodyne detection at FSO receivers. Now plugin (28) in (31), and performing various mathematical computations by using (Gradshteyn and Ryzhik, 2000, (3.351),(6.455),(7.813), (8.352)) to solve the integrals, the final generalized analytical ABER expression is derived aswherewherein and .

6 Theoretical and Simulation Results

In this Section, numerical values from the derived results are obtained and are validated through the simulations performed through Matlab by conducting 10⁸ realizations. The atmospheric turbulence of ES-HAP link is generated by the product of two generalized Gamma distributed random variables. The radial displacement parameter is considered to be Rayleigh distributed and is applied to (8) for pointing error impairments. Throughout the analysis, η = 1 is considered. Unity transmit power is considered at all the nodes and γ_th = 1 dB is considered for the analysis. Further, parameters related to various links are shown in Table 3 unless otherwise stated.

In Figure 3, analytical and simulation results of outage probability against transmit SNR are shown for all three links. In Figure 3A, analytical and simulation results of ES-HAP link are compared for both the IM/DD and heterodyne detection (Het. in figures) techniques for different atmospheric turbulence conditions based on various ground turbulence i.e. , , and , (in m^−2/3) respectively for weak, moderate, and strong ground turbulence conditions. Simulation results match well with the analytical results and validate the accuracy. From Figure 3A, we observe a significant improvement in outage performance while going from strong to weak turbulence for both the detection techniques. Also, heterodyne detection provides significant performance improvement than IM/DD due to its capability in handling turbulence effects more efficiently despite having implementation complexity (Zedini et al., 2016). For the performance analysis, a = 10 cm, ω_L = 50 cm, and σ = 10 cm are considered which corresponds to ϵ = 2.553. The same can be conducted for different values of ϵ by varying these selection parameters to show the impact of pointing error. In Figure 3B, analytical and simulation results of HAP-HAP link are compared for various values of ω_L/a which corresponds to different values of ϵ_o. Simulation results match well with the analytical results and validate the accuracy. Further, selection of appropriate values of a, ω_L, and σ results in an appropriate ϵ_o for desired system performance. Figure 3C compares the analytical and simulation results of HAP to terrestrial MUs link for different values of fading severity m, number of antenna N_t, and number of terrestrial users U. U = 1 gives the performance of uth terrestrial MU while increasing U gives the flexibility of user selection. Considering U = 1, m = 1, N_t = 1 as a reference case, we observe significant performance improvement while increasing m or N_t or both. However, increase in N_t gives better performance than increase in m. Increase in U provides the user selection flexibility and improves the outage performance further.

FIGURE 3

In Figure 4, analytical, simulation, and asymptotic results of overall system outage probability are compared for various values of U, m, and N_t while heterodyne detection is considered at FSO receiver. For the analysis, is considered and ω_L/a = 5 for moderate pointing error (Figure 4A) and ω_L/a = 2 for severe pointing error (Figure 4B) are considered. Rest of the parameters are shown in Table 3. From Figure 4A, considering U = 1, m = 1, N_t = 1 as reference case, for an outage probability of 10^–2, approximately 8.7 and 11.7 dB gains are achieved while increasing m from one to two and N_t from 1 to 2, respectively. This indicates better performance with an increase in N_t than m by approximately 3 dB. Further, the gain is same throughout the SNR range except for U = 1, m = 2, N_t = 2 and U = 2, m = 2, N_t = 2 cases. For U = 2, m = 2, N_t = 2, there must be a significant performance improvement as shown in Figure 3C, however, only slight improvement is observed over U = 1, m = 2, N_t = 2. This is because the ES-HAP FSO link’s performance reaches to its saturation which limits the overall outage performance. Hence, no further performance improvement is achieved with the increase in U, m, and N_t beyond U = 2, m = 2, N_t = 2 case. In Figure 4B, ω_L/a = 2 is maintained for the FSO links for severe pointing error case. From Figure 4B, by considering U = 1, m = 1, N_t = 1 as reference case, we can observe a clear performance improvement for U = 1, m = 2, N_t = 1 or U = 1, m = 1, N_t = 2. For further increase in U, m, or N_t, we can observe slight improvement only in the low and medium SNR range. However, no such improvement can be seen at high SNRs and all the curves merge with each other. This is because the ES-HAP FSO link’s performance saturates early due to severe pointing error which limits the overall outage performance. Simulation results validate the analytical results presented in Figure 4. Further, asymptotic results on outage probability are also shown in Figure 4 which match with the outage results at high SNRs for all the considered cases and validates the high SNR approximation. Further, the diversity order of the overall system is obtained as .

FIGURE 4

In Figure 5, analytical and simulation results of ABER for 4-QAM are compared for various values of U, m, and N_t while heterodyne detection is considered at the FSO receiver. For the analysis, is considered and ω_L/a = 5 for moderate pointing error (Figure 5A) and ω_L/a = 2 for severe pointing error (Figure 5B) are considered. Rest of the parameters are shown in Table 3. From Figure 5A, a significant improvement in ABER performance is observed with the increase in m, N_t, U or all as compared to the U = 1, m = 1, N_t = 1 reference case. Further, simulation results validate the derived analytical results. Considering U = 1, m = 1, N_t = 1 as a reference case, for an ABER of 10^–2, approximately 5.2 and 8.1 dB gains achieved while increasing m from one to two and N_t from 1 to 2, respectively. Hence, increase in N_t gives nearly 2.9 dB gain than the increase in m. Also for U = 2, m = 2, N_t = 2, only marginal improvement is observed over the U = 1, m = 2, N_t = 2 case. This is due to the fact that the ES-HAP FSO link’s performance reaches to its saturation which limits the overall ABER performance. Hence, no further performance improvement is achieved with the increase in U, m, and N_t beyond U = 2, m = 2, N_t = 2 case. In Figure 5B, ω_L/a = 2 is maintained for the FSO links for severe pointing error case. From Figure 5B, by considering U = 1, m = 1, N_t = 1 as reference case, we can observe a clear ABER performance improvement for U = 1, m = 2, N_t = 1 or U = 1, m = 1, N_t = 2. However, for further increase in U, m, or N_t, we can observe slight improvement only in the low and medium SNR range. However, no so such improvement can be seen at high SNRs and all the curves merge with each other. This is because the ES-HAP FSO link’s performance saturates early due to severe pointing error which limits the overall ABER performance. Further, simulation results validate the analytical results presented in Figure 5 for all the cases.

FIGURE 5

In Figure 6, ABER of various modulation schemes are compared against transmit SNR. For OOK modulation, IM/DD detection is considered at the FSO receiver while heterodyne detection is preferred for the rest of the complex modulation schemes. For the analysis, U = 1, m = 1, and N_t = 2 are considered and other parameters are mentioned in Table 3. It is observed that BPSK provides significant gain in ABER performance than OOK as heterodyne detection provides significant performance improvement over IM/DD. To achieve an ABER of 10^–2, BPSK provides around 15 dB gain over the OOK. Further, it is observed that the ABER performance of 4-QAM is same as 4-PSK. However, with the increase in constellation order, MQAM outperforms the MPSK. For an ABER of 10^–3, 16-QAM provides approximately 3.2 dB gain over the 16-PSK. This gain improves further with the increase in M.

FIGURE 6

7 Conclusion

In this work, a mixed FSO/RF based multiple serial HAPs assisted multiuser multiantenna terrestrial communication system was considered, where ES-HAP and HAP-HAP links were FSO links and HAP to terrestrial MUs link was an RF link. Both the heterodyne detection and IM/DD techniques were considered at FSO receivers. The atmospheric turbulence of ES-HAP link was modeled with gamma-gamma fading with pointing error impairments. The HAP-HAP links were assumed to have pointing errors only. Multiantenna array was considered at the Nth HAP and best user selection was performed to improve the performance of RF links for the terrestrial users which were Nakagami-m distributed. For the performance analysis, analytical expressions of individual link’s outage probability, overall outage probability, asymptotic outage probability, and a generalized ABER expression of various modulation schemes were derived. Finally, the impact of various selection parameters like pointing error, FSO detection type, MU selection, RF fading severity, and the number of antennas were observed on their performance. The proposed system model and the analysis shown in this work is generalized. All the parameters can be adjusted and the considered system model can be deployed for various specific wireless communication applications. By considering number of HAPs equal to one, the multi-hop serial HAPs assisted hybrid FSO/RF system can be converted into a dual-hop HAP assisted hybrid FSO/RF system which can be applied for a short range highly reliable and robust high data-rate communication for various wireless applications including the emergency disastrous situations, temporary mass events, and even in military applications. By considering N HAPs in a serial mode, the coverage for the communication can be improved till several hundred kilometers and the problem of connecting the unconnected can be solved.

References

• 1

  AhmadI.NguyenK. D.LetzepisN. (2017). “Performance Analysis of High Throughput Satellite Systems with Optical Feeder Links,” in IEEE Global Commun. Conf. (GLOBECOM), Singapore, 4–8 Dec 2017 (IEEE), 1–7. 10.1109/glocom.2017.8255105

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AntoniniM.BettiS.CarrozzoV.DucaE.RuggieriM. (2006). “Feasibility Analysis of a HAP-LEO Optical Link for Data Relay Purposes,” in IEEE Aero. Conf., Big Sky, MT, USA, 4–11 March 2006 (IEEE), 1–7.

  • Google Scholar
  • View reference in article

• 3

  AshrafzadehB.ZaimbashiA.Soleimani-NasabE.UysalM. (2020). Unified Performance Analysis of Multi-Hop FSO Systems over Double Generalized Gamma Turbulence Channels with Pointing Errors. IEEE Trans. Wireless Commun.19, 7732–7746. 10.1109/twc.2020.3015780

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BeckmannP.SpizzichinoA. (1987). The Scattering of Electromagnetic Waves from Rough Surfaces. Norwood, MA, USA: Artech House Radar Library; Artech Print on Demand.

  • Google Scholar
  • View reference in article

• 5

  BhatnagarM. R.ArtiM. K. (2013). Performance Analysis of Hybrid Satellite-Terrestrial FSO Cooperative System. IEEE Photon. Technol. Lett.25, 2197–2200. 10.1109/lpt.2013.2282836

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  ChenJ.YangL.WangW.YangH.-C.LiuY.HasnaM. O.et al (2019). A Novel Energy Harvesting Scheme for Mixed FSO-RF Relaying Systems. IEEE Trans. Veh. Technol.68, 8259–8263. 10.1109/tvt.2019.2925736

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  DatsikasC. K.PeppasK. P.SagiasN. C.TombrasG. S. (2010). Serial Free-Space Optical Relaying Communications over Gamma-Gamma Atmospheric Turbulence Channels. J. Opt. Commun. Netw.2, 576–586. 10.1364/jocn.2.000576

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  DavidF.GiggenbachD.HennigerH.HorwathJ.LandrockR.PerlotN. (2004). “Design Considerations for Optical Inter-HAP Links,” in AIAA Int. Commun. Satellite Syst. Conf. Exhibit (ICSSC), 3144. 10.2514/6.2004-3144

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  d’OliveiraF. A.MeloF. C. L. D.DevezasT. C. (2016). High-Altitude Platforms—Present Situation and Technology Trends. J. Aero. Technol. Manage.8, 249–262. 10.5028/jatm.v8i3.699

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  FaridA. A.HranilovicS. (2007). Outage Capacity Optimization for Free-Space Optical Links with Pointing Errors. J. Lightwave Technol.25, 1702–1710. 10.1109/jlt.2007.899174

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  FidlerF.KnapekM.HorwathJ.LeebW. R. (2010). Optical Communications for High-Altitude Platforms. IEEE J. Select. Top. Quan. Electron.16, 1058–1070. 10.1109/jstqe.2010.2047382

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  GaoN.JinS.LiX.MatthaiouM. (2021). Aerial RIS-Assisted High Altitude Platform Communications. IEEE Wireless Commun. Lett.10, 2096–2100. 10.1109/lwc.2021.3091164

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  GradshteynI.RyzhikI. (2000). Table of Integrals, Series and Products. 6th ed. New York, NY, USA: Academic.

  • Google Scholar
  • View reference in article

• 14

  ITU-Radiocommun (2000). Preferred Characteristics of Systems in the Fixed Service Using High Altitude Platforms Operating in the Bands 47.2-47.5 GHz and 47.9-48.2 GHz. Geneva: International Telecommunication Union, Electronic Publication. https://www.itu.int/dms_pubrec/itu-r/rec/f/R-REC-F.1500-0-200005-I!!PDF-E.pdf.

  • Google Scholar
  • View reference in article

• 15

  JungK.-J.NamS. S.AlouiniM.-S.KoY.-C. (2020). Unified Finite Series Approximation of FSO Performance Over Strong Turbulence Combined with Various Pointing Error Conditions. IEEE Trans. Commun.68, 6413–6425. 10.1109/tcomm.2020.3008459

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  KongH.LinM.ZhuW.-P.AmindavarH.AlouiniM.-S. (2020). Multiuser Scheduling for Asymmetric FSO/RF Links in Satellite-UAV-Terrestrial Networks. IEEE Wireless Commun. Lett.9, 1235–1239. 10.1109/lwc.2020.2986750

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  KurtG. K.KhoshkholghM. G.AlfattaniS.IbrahimA.DarwishT. S. J.AlamM. S.et al (2021). A Vision and Framework for the High Altitude Platform Station (HAPS) Networks of the Future. IEEE Commun. Surv. Tutorials23, 729–779. 10.1109/comst.2021.3066905

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  LeeJ.-H.ParkK.-H.KoY.-C.AlouiniM.-S. (2020). Throughput Maximization of Mixed FSO/RF UAV-Aided Mobile Relaying with a Buffer. IEEE Trans. Wireless Commun.20, 683–694. 10.1109/TWC.2020.3028068

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  LiR.ChenT.FanL.DangA. (2019). Performance Analysis of a Multiuser Dual-Hop Amplify-And-Forward Relay System with FSO/RF Links. J. Opt. Commun. Netw.11, 362–370. 10.1364/jocn.11.000362

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  LiuX.LinM.ZhuW.-P.WangJ.-Y.UpadhyayP. K. (2020). Outage Performance for Mixed FSO-RF Transmission in Satellite-Aerial- Terrestrial Networks. IEEE Photon. Technol. Lett.32, 1349–1352. 10.1109/lpt.2020.3025452

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  MichailidisE. T.NomikosN.BithasP.VouyioukasD.KanatasA. G. (2018). “Outage Probability of Triple-Hop Mixed RF/FSO/RF Stratospheric Communication Systems,” in Proc. IEEE Int. Conf. Adv. Satellite Space Commun. (SPACOMM), Athens, Greece, April 22–26, 2018, 1–6.

  • Google Scholar
  • View reference in article

• 22

  MilasV.KolettaM.ConstantinouP. (2003). Interference and Compatibility Studies Between Satellite Service Systems and Systems Using High Altitude Platform Stations. ESA Spec. Publ.541, 34.

  • Google Scholar
  • View reference in article

• 23

  MiuraR.OodoM. (2001). Wireless Communications System Using Stratospheric Platforms: R and D Program on Telecom and Broadcasting System Using High Altitude Platform Stations. J. Commun. Res. Lab.48, 33–48.

  • Google Scholar
  • View reference in article

• 24

  MozaffariM.SaadW.BennisM.NamY.-H.DebbahM. (2019). A Tutorial on UAVs for Wireless Networks: Applications, Challenges, and Open Problems. IEEE Commun. Surv. Tutorials21, 2334–2360. 10.1109/comst.2019.2902862

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  ParthasarathyS.GiggenbachD.KirstädterA. (2014). “Channel Modelling for Free-Space Optical Inter-HAP Links Using Adaptive ARQ Transmission,” in Unmanned/Unattended Sensors and Sensor Networks X (International Society for Optics and Photonics), Amsterdam, Netherlands, 17 October 2014, 9248, 92480Q. 10.1117/12.2067195

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  Saleh AltubaishiE.AlhamawiK. (2019). Capacity Analysis of Hybrid AF Multi-Hop FSO/RF System Under Pointing Errors and Weather Effects. IEEE Photon. Technol. Lett.31, 1304–1307. 10.1109/lpt.2019.2926679

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  SharmaM.ChadhaD.ChandraV. (2016). High-Altitude Platform for Free-Space Optical Communication: Performance Evaluation and Reliability Analysis. J. Opt. Commun. Netw.8, 600–609. 10.1364/jocn.8.000600

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  SharmaS.MadhukumarA. S.R.S. (2019). Switching-Based Cooperative Decode-And-Forward Relaying for Hybrid FSO/RF Networks. J. Opt. Commun. Netw.11, 267–281. 10.1364/jocn.11.000267

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  ShibataY.KanazawaN.KonishiM.HoshinoK.OhtaY.NagateA. (2020). System Design of Gigabit HAPS Mobile Communications. IEEE Access8, 157995–158007. 10.1109/access.2020.3019820

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  SingyaP. K.AlouiniM.-S. (2021). Performance of UAV Assisted Multiuser Terrestrial-Satellite Communication System Over Mixed FSO/RF Channels. IEEE Trans. Aero. Electron. Syst. 10.1109/taes.2021.3111787

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  SingyaP. K.KumarN.BhatiaV.AlouiniM.-S. (2019). On Performance of Hexagonal, Cross, and Rectangular QAM for Multi-Relay Systems. IEEE Access7, 60602–60616. 10.1109/access.2019.2915375

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  SingyaP. K.KumarN.BhatiaV.AlouiniM.-S. (2020). On the Performance Analysis of Higher Order QAM Schemes Over Mixed RF/FSO Systems. IEEE Trans. Veh. Technol.69, 7366–7378. 10.1109/tvt.2020.2990747

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  SingyaP. K.KumarN.BhatiaV. (2018). Impact of Imperfect CSI on ASER of Hexagonal and Rectangular QAM for AF Relaying Network. IEEE Commun. Lett.22, 428–431. 10.1109/lcomm.2017.2778153

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  SwaminathanR.SharmaS.VishwakarmaN.MadhukumarA. (2021). HAPS-Based Relaying for Integrated Space-Air-Ground Networks with Hybrid FSO/RF Communication: A Performance Analysis. IEEE Trans. Aero. Electron. Syst.57, 1581–1599. 10.1109/TAES.2021.3050663

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  TrichiliA.CoxM. A.OoiB. S.AlouiniM.-S. (2020). Roadmap to Free Space Optics. J. Opt. Soc. Am. B37, A184–A201. 10.1364/josab.399168

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  TrichiliA.RaghebA.BriantcevD.EsmailM. A.AltamimiM.AshryI.et al (2021). Retrofitting FSO Systems in Existing RF Infrastructure: A Non-Zero-Sum Game Technology. IEEE Open J. Commun. Soc.2, 2597–2615. 10.1109/ojcoms.2021.3130645

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  VaiopoulosN.SandalidisH. G.VaroutasD. (2013). Using a HAP Network to Transfer WiMAX OFDM Signals: Outage Probability Analysis. J. Opt. Commun. Netw.5, 711–721. 10.1364/jocn.5.000711

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  WangP.CaoT.GuoL.WangR.YangY. (2015). Performance Analysis of Multihop Parallel Free-Space Optical Systems Over Exponentiated Weibull Fading Channels. IEEE Photon. J.7, 1–17. 10.1109/jphot.2015.2396115

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  Wolframe (1998). The Wolfram Function Site. Available at: http://functions.wolfram.com/.

  • Google Scholar
  • View reference in article

• 40

  XuG.SongZ. (2020). Performance Analysis for Mixed κ-μ Fading and M-Distribution Dual-Hop Radio Frequency/Free Space Optical Communication Systems. IEEE Trans. Wireless Commun.20, 1517–1528. 10.1109/TWC.2020.3034104

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  XuG.SongZ. (2021). Performance Analysis of a UAV-Assisted RF/FSO Relaying Systems for Internet of Vehicles. IEEE Internet Things J.10.1109/JIOT.2021.3051211

  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  XuG.ZhangQ. (2021). Mixed RF/FSO Deep Space Communication System Under Solar Scintillation Effect. IEEE Trans. Aero. Electron. Syst.57, 3237–3251. 10.1109/taes.2021.3074130

  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  YangL.YuanJ.LiuX.HasnaM. O. (2018). On the Performance of LAP-Based Multiple-Hop RF/FSO Systems. IEEE Trans. Aero. Electron. Syst.55, 499–505. 10.1109/TAES.2018.2852399

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  ZediniE.AlouiniM.-S. (2015). Multihop Relaying over IM/DD FSO Systems with Pointing Errors. J. Lightwave Technol.33, 5007–5015. 10.1109/jlt.2015.2492244

  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  ZediniE.KammounA.AlouiniM.-S. (2020). Performance of Multibeam Very High Throughput Satellite Systems Based on FSO Feeder Links with HPA Nonlinearity. IEEE Trans. Wireless Commun.19, 5908–5923. 10.1109/twc.2020.2998139

  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  ZediniE.SouryH.AlouiniM.-S. (2016). On the Performance Analysis of Dual-Hop Mixed FSO/RF Systems. IEEE Trans. Wireless Commun.15, 3679–3689. 10.1109/twc.2016.2524685

  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  ZhanP.YuK.SwindlehurstA. L. (2011). Wireless Relay Communications with Unmanned Aerial Vehicles: Performance and Optimization. IEEE Trans. Aerosp. Electron. Syst.47, 2068–2085. 10.1109/taes.2011.5937283

  • CrossRef
  • Google Scholar
  • View reference in article