Next Article in Journal
Patent Portfolio Analysis of the Synergy between Machine Learning and Photonics
Previous Article in Journal
High Sensitivity Fiber Gas Pressure Sensor with Two Separated Fabry–Pérot Interferometers Based on the Vernier Effect
Previous Article in Special Issue
Hybrid Ring- and Tree-Topology RoF Transmission System with Disconnection Protection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mixed FSO/RF Integrated Satellite-High Altitude Platform Relaying Networks for Multiple Terrestrial Users with Presence of Eavesdropper: A Secrecy Performance

1
Department of Electrical and Electronic Engineering, Faculty of Technology, University of Ibadan, Ibadan 0234, Nigeria
2
Department of Computer Systems Engineering, Tshwane University of Technology, Pretoria 0001, South Africa
*
Author to whom correspondence should be addressed.
Submission received: 27 October 2021 / Revised: 10 December 2021 / Accepted: 23 December 2021 / Published: 4 January 2022
(This article belongs to the Special Issue Optical Wireless Communication (OWC) Systems)

Abstract

:
In this paper, the secrecy performance of a mixed free space optical (FSO)/radio frequency (RF) integrated satellite-high altitude platform (HAP) relaying networks for terrestrial multiusers with the existence of an eavesdropper is investigated. In this network, FSO is adopted to establish the link between the satellite and HAP for which it experiences Gamma-Gamma distributions under different detection schemes (i.e., heterodyne and intensity modulation direct detection). The transmission between the amplify-and-forward (AF) relaying HAP and terrestrial multiusers is through the RF and is modeled as shadowed-Rician fading distribution. Owning to broadcasting nature of RF link, it is assumed that an eavesdropper attempts to intercept the users’ confidential message, and the eavesdropper link is subjected to Rician distributions. Specifically, the closed-form expression for the system equivalent end-to-end cumulative distribution function is derived by exploiting the Meijer’s G and Fox’s H functions. Based on this expression, the exact closed-form expressions of the system connection outage probability, secrecy outage probability, and strictly positive secrecy capacity are obtained under the different detection schemes at HAP. Moreover, the asymptotic analyze of the system secrecy outage probability is provided to obtain more physical insights. Furthermore, the accuracy of all the derived analytical closed-form expressions is verified through the Monte-Carlo simulations. In addition, the impact of atmospheric turbulence, pointing errors, shadowing severity parameters, and Rician factor are thoroughly evaluated. Under the same system conditions, the results depict that heterodyne detection outperforms the intensity modulation direct detection.

1. Introduction

Recently, there is tremendous growth in the deployment of emergency wireless communication supporting services for large scale disasters in order to provide vital information for both rescue teams and survivors [1]. In this scenario, a fast and reliable communication is highly difficult to establish owing to geographical complexity of the disaster area where the communication infrastructures of such base stations and optical fibers have been destroyed [2]. As a result of this, satellite communication has been suggested as most promising solution due to its large coverage footprint and reliable wireless connectivity for various fixed and mobile users [3]. Conventionally, satellite communication systems are RF-based, in which different frequency bands are used for a particular application. However, this RF-based systems suffer from spectrum congested, licensed band, low security, and interference with other frequency bands [4]. To overcome this problems, FSO has been recently proposed for satellite communication link owing to its ability to offer extremely large bandwidth, high security, unlicensed spectrum, low interference, ease of deployment, etc., compared to its RF counterpart [5]. Despite these unique features, the FSO link usually suffers from atmospheric turbulence due to variation in refractive index as the result of inhomogeneity in atmospheric temperature and pressure that can degrade the system performance [6]. In addition, the link is highly prone to pointing error caused by building sways, leading to misalignment between the transmitter and receiver [7]. As a directional link, it requires line-of-sight (LOS) between the transmitter and receiver, which may be stochastically blocked by the moving, surrounding objects, such as birds, airplanes, etc. [8].
High altitude platforms (HAPs) have attracted immense attention in the research community due to their inherit advantages from both the terrestrial and satellite communication systems [9]. They have been suggested as a future candidate for communication system to offer efficient broadband wireless access at low cost [10]. As flying aerial platforms, such as airships, planes, and balloons; they are usually deployed in stratosphere at altitude between 17   to   22   km above the earth’s surface, where there is existence of strong LOS between the satellite and HAP for efficient wireless transmission [11]. According to International Telecommunication Union (ITU), a larger footprint of 100   km radius has been suggested to act as HAPs, and this can replace various terrestrial base stations [12]. Regardless of the communication link, satellite communication system link usually experiences a masking effect due to unavailability of LOS between the satellite and user equipment. This occurs as a result of obstacles and shadowing, leading to poor quality of services and long delay propagation, especially in the case of geostationary satellites [13]. To overcome this effect, an integrated satellite-high altitude platform (ISHAP) mixed FSO/RF relaying scheme has been proposed in literature, where the RF communication can be used in conjunction with the FSO counterpart. This combined technology can significantly improve the performance of wireless communication networks by taking advantage of the complementary benefits of both connections [14] In this concept, therefore, HAPs can, thus, serve as intermediate relay node to enhance the transmission between the satellite and the ground users by using decode-and-forward (DF) and amplify-and-forward (AF) relaying protocols [13].
Today, security is the most critical challenge in wireless communication system due to the broadcasting nature of RF medium [15]. As a result of this, the ISHAP system is highly susceptible to any external eavesdroppers which have the capacity to intercept the confidential information of the legitimate users via the system RF link. In this view, physical layer security has been suggested to offer better security performance by exploiting the channel characteristics and imperfections without data encryption [16]. Scanning through the open literature, it was found that there have been intense studies on the performance of cooperative satellite-terrestrial relaying networks. A dual-hop satellite communication system with AF relay was presented in Reference [17], where maximum ratio transmission and maximum ratio combiner scheme were adopted at the system source and destination, respectively. In Reference [18], the performance of a hybrid satellite-terrestrial free space optical cooperative network was evaluated, where the system error rate was obtained. In addition, the performance of an AF hybrid satellite terrestrial system over generalized fading channels was investigated in Reference [19], and the authors obtained the analytical expression for the average symbol error rate of the network. In Reference [20], the impact of hardware impairment on the reliability performance of hybrid satellite-terrestrial system was evaluated, where the geosynchronous orbit satellite sends its information to the multiple terrestrial relay nodes. Furthermore, the performance of hybrid satellite-terrestrial network having multi-antennas terrestrial relays was investigated in Reference [21]. It was assumed that there was direct link between the satellite and the destination, and the effect of hardware impairments with the interference on both relay and destination were investigated. The outage performance of a partial relay selection-based hybrid satellite-terrestrial network was analyzed in Reference [22]. In Reference [23], the performance of an underlay hybrid cognitive satellite terrestrial network was evaluated, where the primary satellite coexists with the second terrestrial relay system. The performance of joint impact of CSI and hardware impairments on a cognitive satellite terrestrial relay network was presented in Reference [24]. However, all these aforementioned works only focused on the scenario of single terrestrial user. In this view, future satellite communication network is required to provide reliable and efficient services for a large number of users. As a result of this, many research studies are focusing on evaluating the performance of hybrid satellite-terrestrial multiuser networks. Thus, the performance of a multi-antenna multiusers hybrid satellite terrestrial relay network was studied in Reference [25]. The authors employed opportunistic user scheduling with outdated channel state information (CSI) and AF relaying with co-channel interference. The outage performance of a multi-antenna multiple users integrated in a satellite-unmanned aerial vehicle (UAV) cooperative network was presented in Reference [13]. However, all these works were RF-based multiusers hybrid satellite-terrestrial systems. In addition, the outage performance of a mixed FSO/RF satellite-aerial terrestrial network was evaluated in Reference [26], where the satellite sent messages to multiple users with the help of a multi-antenna aerial platform as a relay. In Reference [27], the performance analysis of a mixed FSO/RF satellite-UAV multiuser network was studied, where the UAV relay node adopts DF relaying protocol. Kong et al. [28] evaluated the outage probability performance of an asymmetric FSO/RF of multiusers satellite-UAV-terrestrial system, and the authors considered DF relay UAV node for transmission and RF link exploit transmit beamforming scheme. However, in all these previous research works on a mixed FSO/RF, multiusers integrated satellite-terrestrial systems are DF-based relay node and the secrecy performance of cooperative satellite terrestrial relaying networks was not considered.
Recently, some research studies have examined hybrid satellite-terrestrial systems in the context of physical layer security. The secrecy performance of a hybrid satellite relay network with multiple external eavesdroppers was studied [29] with the assumption that multiple eavesdroppers receive signals from the terrestrial networks. In Reference [30], the secrecy outage performance of a land mobile satellite system was evaluated, where both the legitimate user and eavesdropper utilize the diversity technique with multiple antennas. Moreover, the physical layer secrecy performance of a hybrid satellite and terrestrial free space optical cooperative network was presented in Reference [31], where the eavesdropper was assumed to intercept the satellite signal. In addition, the secrecy performance of RF eavesdropping for a HAP aided satellite system was investigated in Reference [32]. However, the security performance of all aforementioned studies was not based on multiple legitimate users, and DF relay was considered for the system. Motivated by this, the performance of the secrecy performance of a mixed FSO/RF integrated satellite-high altitude platform (HAP) relaying networks for terrestrial multiusers with the existence of an eavesdropper. The HAP relay node is considered to utilize AF relaying with different detection schemes for FSO signal detection. The closed-form expression for the system equivalent end-to-end cumulative distribution function is derived. Thus, to quantify the system performance, the exact closed-form expressions of the system connection outage probability, secrecy outage probability, and strictly positive secrecy capacity are obtained under the different HAP detection schemes. To gain better insight, the asymptotic expression for the system SOP is derived. The main contributions of this work are summarized as follows:
  • The exact closed-form expression for the system equivalent end-to-end cumulative distribution function (CDF) is derived by exploiting the Meijer’s G and Fox’s H functions. To the best of authors’ knowledge, the derived CDF is novel as a mixed Gamma-Gamma and shadowed-Rician structure under AF relaying protocol is not found in the existing literature.
  • The analytical closed-form expressions of the system performance in terms of connection outage probability (COP), secrecy outage probability (SOP), and strictly positive secrecy capacity (SPSC) are obtained.
  • The asymptotic expression for the system secrecy outage probability is derived to obtain physical insight.
  • Relative to Reference [26], where multiple users are considered, the system performance was not based on physical layer secrecy. In addition, the RF link was subjected to Nakagami-m distributions. In this paper, the system performance is based on physical layer secrecy, and the RF links are subjected to Rician fading distributions.
  • Relative to Reference [32], where the HAP relay aided node employed DF relaying protocol, and the system only considered a single terrestrial user. The HAP in this paper is an AF-based relaying protocol, and the multiple legitimate users are considered at the ground station.
The rest of the paper is organized as follows. The system and channel models are provided in Section 2. In Section 3, the exact closed-form expression for the system equivalent end-to-end cumulative distribution function is presented We derive the analytical expressions of the COP, SOP, and SPSC subject to different detection schemes in Section 4. Numerical results and discussions are depicted in Section 5. Finally, concluding remarks are detailed in Section 6.

2. System and Channel Models

A mixed FSO/RF integrated satellite-HAP relaying network is illustrated in Figure 1, where the satellite (S) transmits confidential information to the legitimate terrestrial users (D) through an AF HAP relay (R) in the presence of an eavesdropper attempting to wiretap confidential information. It is assumed that there is no direct link between satellite and the terrestrial users due to masking effects as a result of weather, environmental obstacles, long distance, etc. It is assumed that the FSO link follows Gamma-Gamma fading distribution, while the RF link between the HAP relay and the terrestrial users follows Shadowed-Rician fading distribution. In addition, the eavesdroppers RF link between the E and R is considered to undergo Rician distribution. In addition, the HAP relay in the network employs AF relaying protocol, which will introduce a fixed gain into the receive signal irrespective of the fading amplitude on the FSO link. Thus, the network end-to-end signal-to-noise ratio (SNR) can be defined as [33]:
γ e q = γ 1 γ 2 γ 2 + C ,
where C denotes the fixed relay gain parameter, γ 1 = γ ¯ 1 | h 1 | 2 represents the instantaneous SNR at the HAP with γ ¯ 1 , and | h 1 | 2 , respectively, signifies the average SNR and channel power gain of satellite-to-HAP relay link, while the γ 2 = γ ¯ 2 h 2 indicates the instantaneous SNR at the destination, with γ ¯ 2 and h 2 representing the average SNR and channel power gain of HAP-to-terrestrial users link, respectively.

2.1. FSO Link Statistical Distributions

Since the FSO link experiences Gamma-Gamma fading distribution, the probability density function (PDF) of the instantaneous SNR γ 1 under different detection techniques can be defined as [5]:
f γ 1 ( γ 1 ) = ω γ 1 G 1 , 3 3 , 0 ( M ( γ 1 μ r ) 1 / r | ξ 2 + 1 ξ 2 , α , β ) ,
where ω = ξ 2 r ( Γ ( α ) Γ ( α ) ) , M = ξ 2 α β ( 1 + ξ 2 ) , Γ ( . ) signifies the Gamma function, ξ signifies the pointing error, α and β are the scintillation parameters which are specified in references [34,35], the r parameter shows the type of detection at the HAP relay node (i.e., r = 1 for heterodyne detection, and r = 2 for IM/DD), μ r indicates the average electrical SNR of the link. Specifically, for heterodyne detection, μ 1 = γ ¯ 1 with γ ¯ 1 represents the average instantaneous SNR of the link, and, for IM/DD, μ 2 = γ ¯ 1 α β ξ 2 ( ξ 2 + 1 ) ( ξ 2 + 1 ) ( α + 1 ) ( β + 1 ) .
By integrating (2) through the integral identity detailed in Reference [36] (Equation (26)), the cumulative density function (CDF) of the instantaneous SNR γ 1   can be expressed as:
F γ 1 ( γ 1 ) = ξ 2 r α + β 2 2 π r 1 ( Γ ( α ) Γ ( α ) ) G r + 1 , 3 r + 1 3 r , 1 ( ( α β ) r r 2 r μ r γ 1 | 1 ,   ψ 1 ψ 2 , 0 ) ,
where ψ 1 = Δ ( r ,   ξ 2 + 1 ) and ψ 2 = Δ ( r ,   ξ 2 ) ,   Δ ( r ,   α ) ,   Δ ( r , β ) with Δ ( ψ , a ) = a ψ , a + 1 ψ , ,   a + ψ 1 ψ .

2.2. RF Link Statistical Distributions

It is assumed that the RF link between the HAP and terrestrial users link follows Shadowed-Rician fading distribution and the PDF of instantaneous SNR γ 2 can be expressed as [17,25,37]:
f γ 2 ( γ 2 ) = e x p ( a γ 2 ) 1 F 1 ( m h ; 1 ; δ γ 2 ) ,   γ 2 > 0
where = 1 / 2 ( 2 b m h 2 b m h + Ω ) m h , a = 1 / 2 b and δ = Ω / 2 b ( 2 b m h + Ω ) , with Ω , m h , and 2 b being the average power of LOS, Nakagami fading severity parameter with 0 < m h < , and the multipath component, respectively, and   1 F 1 ( x ; y ; z ) denotes the confluent hypergeometric function. By using the identity detail in Reference [38], the confluent hypergeometric function in (4) can be expressed as:
  1 F 1 ( m h ; 1 ; δ γ 2 ) = e x p ( δ γ 2 ) q = 0 m h 1 ( 1 m h ) q ( δ γ S R ) q ( q ! ) 2 ,
where ( . ) q denotes the Pochhammer symbol.
By invoking (5) into (4), the Shadowed-Rician distribution of the RF link can be further simplified as:
f γ 2 ( γ 2 ) = q = 0 m h 1 ( 1 m h ) p ( δ ) q ( p ! ) 2 ( γ ¯ 1 ) q + 1 γ 2 q exp ( Δ H γ S R ) ,
where Δ H = ( a δ ) / γ ¯ 2 .
Since the sum of SNR at the terrestrial users can be expressed as γ 2 = q = 1 N u | h 2 , q | 2   , then, by following the same approach detailed in Reference [17], the PDF of instantaneous SNR γ 2 Shadowed-Rician distribution can be expressed as:
f γ 2 ( γ 2 ) = q 1 = 0 m h 1 q N u = 0 m h 1 Ξ ( N u ) γ 2 λ 1 exp ( Δ H γ 2 ) ,
where
{ Ξ ( N u ) = j = 1 N u ρ ( q j ) N u t = 1 N u 1 B ( k = 1 ε q k + ε , q ε + 1 + 1 ) ρ ( q j ) = ( 1 m h ) q j ( δ ) q j ( q j ! ) 2 ( γ ¯ 2 ) q j + 1 λ = i = 1 N u q j + N u ,
with B ( . , . ) denoting Beta function.
By integration (7) using the integral identity detailed in Reference [38] (Equation (3.351(1), the CDF of the link can be obtained as:
F γ 2 ( γ 2 ) = q 1 = 0 m h 1 q N u = 0 m h 1 Ξ ( N u ) Δ H λ γ ( λ , Δ H γ 2 ) ,
where γ ( . , . ) is the lower incomplete Gamma function.
By converting the incomplete Gamma function to Meijer-G function using the identity detailed in Reference [39] (Equation (8.4.16(1))), then, (9) can be further expressed as:
F γ 2 ( γ 2 ) = q 1 = 0 m h 1 q N u = 0 m h 1 Ξ ( N u ) Δ H λ G 1 , 2 1 , 1 ( Δ H γ 2 | 1 λ ,         0 ) .
Owing to broadcast nature of RF link, the eavesdropper attempts to intercept the secret information of the ground users from the HAP relay node. As a result of this, the eavesdropper link is assumed to follow Rican distribution with the PDF of instantaneous SNR γ e defined as [40,41]:
f e ( γ e ) = ( K + 1 ) exp ( K ) γ ¯ e e x p ( K + 1 γ ¯ e ) I o ( 2 K ( K + 1 ) γ ¯ e γ e ) ,
where K is the Rician fading factor which is defined as the ratio of the power of the line-of-sight component to the scattered components, and I o ( . )   denotes the zero-order modified Bessel function of the first kind. By using the identity detailed in Reference [38], I o = l = 0 x l / 2 2 l ( l ! ) 2 ,
f e ( γ e ) = η e x p ( K ) l = 0 ( K η ) l ( l ! ) 2 γ e l exp ( η γ e ) ,
where η = K + 1 γ ¯ e .

3. Statistics of End-to-End SNR

Since the HAP relay node employs AF relaying protocol, the equivalent end-to-end CDF of the concerned system can expressed as [33,42]:
F e q ( γ ) = 0 P [ γ e q = γ 1 γ 2 γ 2 + C < γ | γ 1 ] f γ 1 ( γ 1 ) d γ 1 .
By variable transformation of y = γ 1 γ , (13) can be further expressed as [33,43]:
F e q ( γ ) = F γ 1 ( γ 1 ) + 0 F γ 2 ( C γ y ) f γ 1 ( y + γ ) d y Υ .
By putting (2), (3), and (10) into (13), the integral term of (14) can be expressed as:
Υ = ω q 1 = 0 m h 1 q N u = 0 m h 1 Ξ ( N u ) Δ H λ 0 G 1 , 2 1 , 1 ( Δ H C γ y | 1 λ ,         0 ) ( y + γ ) 1 G 1 , 3 3 , 0 ( M ( ( y + γ ) μ r ) 1 / r | ξ 2 + 1 ξ 2 , α , β ) d y .
By applying the identity defined in Reference [44] (Equation (1.112)), the Meijer-G functions in (15) can be expressed as follows:
Υ = ω ( 2 π i ) 2 q 1 = 0 m h 1 q N u = 0 m h 1 Ξ ( N u ) Δ H λ L 1 L 2 Γ ( λ + s 1 ) Γ ( s 1 ) Γ ( 1 s 1 ) Γ ( ξ 2 + s 2 ) Γ ( α + s 2 ) Γ ( β + s 2 ) Γ ( ξ 2 + 1 + s 2 )   × ( Δ H C γ ) s 1 ( M μ r 1 / r ) s 2 0 y s 1 ( y + γ ) s 2 r 1 d y d s 1 d s 2 .
By applying the identity in Reference [38] (Equations (3.193(3)) and (8.384(1))), (16) can be solved as:
Υ = ω ( 2 π i ) 2 q 1 = 0 m h 1 q N u = 0 m h 1 Ξ ( N u ) Δ H λ L 1 L 2 Γ ( λ + s 1 ) Γ ( s 1 ) Γ ( 1 s 1 ) Γ ( ξ 2 + s 2 ) Γ ( α + s 2 ) Γ ( β + s 2 ) Γ ( ξ 2 + 1 + s 2 )   × Γ ( s 1 + 1 ) Γ ( s 2 r s 1 ) Γ ( s 2 r + 1 ) ( Δ H C ) s 1 ( γ 1 / r M μ r 1 / r ) s 2 d s 1 d s 2 .
By letting s = s 1 and t = s 2 , then, (17) can be rewritten as:
Υ = ω ( 2 π i ) 2 q 1 = 0 m h 1 q N u = 0 m h 1 Ξ ( N u ) Δ H λ L 1 L 2 Γ ( t r s ) Γ ( λ + s ) Γ ( s ) Γ ( s + 1 ) Γ ( 1 s ) × Γ ( ξ 2 + t ) Γ ( α + t ) Γ ( β + t ) Γ ( ξ 2 + 1 + t ) Γ ( t r + 1 ) ( Δ H C ) s ( γ 1 / r M μ r 1 / r ) t d s d t .
By applying the identity detailed in Reference [45], (18) can be obtained as:
Υ = ω q 1 = 0 m h 1 q N u = 0 m h 1 Ξ ( N u ) Δ H λ × H 1 , 0 : 3 , 1 : 3 , 2 0 , 1 : 1 , 2 : 0 , 3 ( ( 1 , 1 , 1 r ) | ( 1 λ , 1 ) ( 0 , 1 ) ( 1 , 1 ) ( 0 , 1 ) | ( 1 ξ 2 , 1 ) ( 1 α ,   1 ) ( 1 β , 1 ) ( ξ 2 , 1 ) ( 0 ,     1 r ) |   1 Δ H C , μ r 1 / r γ 1 / r M ) ,
where H p 1 ,   q 1 :   p 2 ,   q 2 :   p 3 ,   q 3 m 1 , n 1 : m 2 , n 2 : m 3 , n 3 [.] is the Extended Generalized Bivariate Fox’s H-function (EGBFHF).
By substituting (3) and (19) into (13), the concerned system equivalent end-to-end CDF can be obtained as:
F e q ( γ ) = ξ 2 r α + β 2 2 π r 1 ( Γ ( α ) Γ ( α ) ) G r + 1 , 3 r + 1 3 r , 1 ( ( α β ) r r 2 r μ r γ | 1 ,   ψ 1 ψ 2 , 0 ) + ω q 1 = 0 m h 1 q N u = 0 m h 1 Ξ ( N u ) Δ H λ × H 1 , 0 : 3 , 1 : 3 , 2 0 , 1 : 1 , 2 : 0 , 3 ( ( 1 , 1 , 1 r ) | ( 1 λ , 1 ) ( 0 , 1 ) ( 1 , 1 ) ( 0 , 1 ) | ( 1 ξ 2 , 1 ) ( 1 α ,   1 ) ( 1 β , 1 ) ( ξ 2 , 1 ) ( 0 ,     1 r ) |   1 Δ H C , μ r 1 / r γ 1 / r M ) ) .

4. Performance Analysis

In this section, the exact closed-form expression of the system connection outage probability (COP), secrecy outage probability, and strictly positive secrecy capacity (SPSC) are derived. In addition, to gain more insight about the derived expression, the asymptotic expression of the concerned system is obtained.

4.1. Connection Outage Probability (COP)

This metric characterizes the attainable reliability performance of the proposed system. It describes a situation, whereby the eavesdropper is unable to decode the confidential message of terrestrial users. This usually occurs when the end-to-end instantaneous SNR of the concerned system falls below a predefined threshold value γ t h and can be formulated as [46]:
C O P = P r [ γ γ t h ] .
Thus, by substituting (20) into (21), the COP of the concerned system under different detection schemes can be obtained as:
C O P = ξ 2 r α + β 2 2 π r 1 ( Γ ( α ) Γ ( α ) ) G r + 1 , 3 r + 1 3 r , 1 ( ( α β ) r r 2 r μ r γ t h | 1 ,   ψ 1 ψ 2 , 0 ) + ω q 1 = 0 m h 1 q N s = 0 m h 1 Ξ ( N u ) Δ H λ × H 1 , 0 : 3 , 1 : 3 , 2 0 , 1 : 1 , 2 : 0 , 3 ( ( 1 , 1 , 1 r ) | ( 1 λ , 1 ) ( 0 , 1 ) ( 1 , 1 ) ( 0 , 1 ) | ( 1 ξ 2 , 1 ) ( 1 α ,   1 ) ( 1 β , 1 ) ( ξ 2 , 1 ) ( 0 ,     1 r ) |   1 Δ H C , μ r 1 / r γ t h 1 / r M ) .

4.2. Secrecy Outage Probability (SOP)

4.2.1. Exact Analysis

These are the most vital performance indices which describe the probability that the instantaneous secrecy capacity falls below a predefined threshold rate of confidential information R s and can be defined mathematically as follows [6]:
S O P = 0 F e q ( Θ γ ) f e ( γ ) d γ ,
where Θ = e R s .
By invoking (12) and (22) into (23), the system SOP can be expressed as:
S O P = I 1 + I 2 .
Thus,
I 1 = ξ 2 r α + β 2 η e x p ( K ) 2 π r 1 ( Γ ( α ) Γ ( α ) ) l = 0 ( K η ) l ( l ! ) 2 0 γ e l exp ( η γ ) G r + 1 , 3 r + 1 3 r , 1 ( ( α β ) r r 2 r μ r γ Θ | 1 ,   ψ 1 ψ 2 , 0 ) d γ .
By utilizing the integral identity detailed in Reference [38] (Equation (7.813(1))), the I 1 term of (25) can be expressed as:
I 1 = ξ 2 r α + β 2 e x p ( K ) 2 π r 1 ( Γ ( α ) Γ ( α ) ) l = 0 K l ( l ! ) 2 G r + 2 , 3 r + 1 3 r , 2 ( ( α β ) r Θ η r 2 r μ r | l , 1 ,   ψ 1 ψ 2 , 0 ) .
Similarly,
I 2 = ω η e x p ( K ) l = 0 q 1 = 0 m h 1 q N u = 0 m h 1 ( K η ) l Ξ ( N u ) Δ H λ ( l ! ) 2 0 γ l exp ( η γ ) × H 1 , 0 : 3 , 1 : 3 , 2 0 , 1 : 1 , 2 : 0 , 3 ( ( 1 , 1 , 1 r ) | ( 1 λ , 1 ) ( 0 , 1 ) ( 1 , 1 ) ( 0 , 1 ) | ( 1 ξ 2 , 1 ) ( 1 α ,   1 ) ( 1 β , 1 ) ( ξ 2 , 1 ) ( 0 ,     1 r ) |   1 Δ H C , μ r 1 / r ( γ Θ ) 1 / r M ) d γ .
By converting the H-Fox function to integral form using the identity detailed in Reference [45], (27) can be expressed as:
I 2 = ω η e x p ( K ) l = 0 q 1 = 0 m h 1 q N u = 0 m h 1 ( K η ) l Ξ ( N u ) Δ H λ ( l ! ) 2 ( 2 π i ) 2 L 1   L 2   Γ ( t r s ) Γ ( λ + s ) Γ ( s ) Γ ( s + 1 ) Γ ( 1 s ) × Γ ( ξ 2 + t ) Γ ( α + t ) Γ ( β + t ) Γ ( ξ 2 + 1 + t ) Γ ( t r + 1 ) ( Δ H C ) s ( Θ 1 / r M μ r 1 / r ) t 0 γ l t / r exp ( η γ ) d γ I 3   d s d t .
By applying the integral identity defined in Reference [38] (Equation (3.326(2))), the I 3 term of (28) can be solved as:
I 3 = Γ ( l t r + 1 ) η ( l t r + 1 ) .
Putting (29) into (28) and apply the identity given in Reference [45], the I 2 term of (28) can be obtained as:
I 2 = ω e x p ( K ) l = 0 q 1 = 0 m h 1 q N u = 0 m h 1 K l Ξ ( N u ) Δ H λ ( l ! ) 2 × H 1 , 0 : 3 , 1 : 3 , 3 0 , 1 : 1 , 2 : 1 , 3 ( ( 1 , 1 , 1 r ) | ( 1 λ , 1 ) ( 0 , 1 ) ( 0 , 1 ) ( 0 , 1 ) | ( 1 ξ 2 , 1 ) ( 1 α ,   1 ) ( 1 β , 1 ) ( ξ 2 , 1 ) ( 0 ,     1 r ) |   1 Δ H C , ( η μ r ) 1 / r Θ 1 / r M ) .
By substituting (26) and (30) into (24), the system SOP can be obtained as:
S O P = ξ 2 r α + β 2 e x p ( K ) 2 π r 1 ( Γ ( α ) Γ ( α ) ) l = 0 K l ( l ! ) 2 G r + 2 , 3 r + 1 3 r , 2 ( ( α β ) r Θ η r 2 r μ r |   ψ 3 ψ 4 ) + ω e x p ( K ) l = 0 q 1 = 0 m h 1 q N u = 0 m h 1 K l Ξ ( N u ) Δ H λ ( l ! ) 2 × H 1 , 0 : 3 , 1 : 3 , 3 0 , 1 : 1 , 2 : 1 , 3 ( ( 1 , 1 , 1 r ) | ( 1 λ , 1 ) ( 0 , 1 ) ( 0 , 1 ) ( 0 , 1 ) | ( 1 ξ 2 , 1 ) ( 1 α ,   1 ) ( 1 β , 1 ) ( ξ 2 , 1 ) ( 0 ,     1 r ) |   1 Δ H C , ( η μ r ) 1 / r Θ 1 / r M ) ,
where ψ 3 = l , 1 ,   ψ 1 , and ψ 4 = ψ 2 , 0 .

4.2.2. Asymptotic Analysis

Due to the complexity of exact closed-form expression of the system SOP given in (31), limited physical insight is provided. As a result of this, the analytical asymptotic expression is derived for the system SOP at high SNR by considering a condition at which μ r and can be derived in closed-form as:
S O P = ξ 2 r α + β 2 e x p ( K ) 2 π r 1 ( Γ ( α ) Γ ( α ) ) l = 0 t = 1 3 r K l ( l ! ) 2 ( ( α β ) r Θ η r 2 r μ r ) ψ 4 , t + ω η e x p ( K ) l = 0 q 1 = 0 m h 1 q N u = 0 m h 1 r ( K η ) l Ξ ( N u ) Δ H λ ( l ! ) 2 ( 2 π i ) ( ( η μ r ) 1 r Θ 1 r M ) r H 4 , 6 5 , 2 ( Θ 1 r M ( η μ r ) 1 r | Δ 1 Δ 2 ) ,
where Δ 1 = { ( 1 , 1 ) ( l , 1 / r ) ( ξ 2 + 1 , 1 ) ( 1 , 1 / r ) } and Δ 2 = { ( λ , 1 ) ( 1 , 1 ) ( ξ 2 , r )   ( α , r ) ( β , r ) ( 0 , 1 ) } .
Proof. 
See Appendix A. □

4.3. Strictly Positive Secrecy Capacity (SPSC)

This illustrates the probability of the existence of positive secrecy capacity so as to offer a secure transmission and can be expressed as [32]:
S P S C = 1 0 F e q ( γ ) f e ( γ ) d γ .
Thus,
S P S C = 1 S O P Θ = 1 .

5. Numerical Results and Discussions

In this section, numerical results of the system secrecy performance in terms of COP, SOP and SPSC under different detection schemes are presented. The shadowing severity for the Shadowed-Rician distributions are defined as follows: ( m = 1 ,   b = 0.063 ,   Ω = 0.0007 ) for the frequent heavy shadowing (FHS), and ( m = 5 ,   b = 0.251 ,   Ω = 0.279 ) for the average shadowing (AS) [25]. Regarding the influence of atmospheric turbulence, different turbulence conditions are considered, and these include: weak levels ( α = 3.78 ,   β = 3.74 ) , moderate levels ( α = 2.50 ,   β = 2.06 ) , and strong levels ( α = 2.04 ,   β = 1.10 ) [7]. Except otherwise stated, R s = 0.01   nats / s / Hz , N u = 2 , C = 0.75 , and the system is subjected to strong turbulence and FHS fading conditions.
The COP performance of the system under the variation of threshold SNR γ t h is illustrated in Figure 2 under different atmospheric turbulence. It can be deduced from the results that the increase in atmospheric turbulence significantly deteriorates the system COP performance for the two detection schemes at the HAP. In addition, a similar impact of atmospheric turbulence on the system SOP performance can be observed in Figure 3 as the system SOP gets degraded due to severities of the atmospheric turbulence. In both secrecy performances, it is depicted by the results that heterodyne detection offers the system better performance than IMDD detection. The results also depict that the analytical results perfectly matched with the simulation results, which indicate the accuracy of the derived expressions.
The impact of pointing error on the system SOP performance under the different detection schemes is presented in Figure 4. Under both detections, the lower pointing error value ( ξ = 1.1 ) causes an increase in the effect of pointing errors, which leads to degradation in the system SOP performance. On the other hand, higher values of PE ( ξ = 6.5 ) decrease the effect of PE and improve the system SOP. It is observed from the results that the system performance is better under the heterodyne detection compared to the IMDD detection.
In Figure 5, the effect of γ ¯ e level on the system SOP performance is demonstrated under the different detection schemes. It can be seen that the analytical results collaborate perfectly with the simulation results, which proves the accuracy of the derived expressions. Moreover, it can be deduced from the results that the higher the values of γ ¯ e , the worse the system SOP performance, with heterodyne detection outperforming the IMDD detection. It can also be observed from the result that the asymptotic results follow the SOP results at high SNR. This shows that the quality of eavesdropper channel improves at higher valves of γ ¯ e , which significantly leads to degradation in system SOP performance.
The effect of Rican factor level of the eavesdropper link on the system SOP performance is depicted in Figure 6. It can be observed that the higher value of   K factor has lower system SOP performance. This is because, at large values of K   factor, there is strong LOS for the eavesdropper to wiretap the confidential message from HAP, and this degrades the system SOP performance. The analytical results match well with simulation results, confirming the validity of the derived expressions.
The system SOP performance under different shadowing fading severities is presented in Figure 7. It can be observed that the system becomes degraded as the shadowing fading effect increases from AS to FHS. Under the same channel conditions, the results show that heterodyne detection offers the system better SOP performance compared with IMDD detection. In addition, analytical results of the system SOP closely match with simulation results, showing the accuracy of our derivations.
The effect of the number of terrestrial users on the system SOP performance under strong turbulence condition is presented in Figure 8. It can be deduced that the system SOP performance significantly improves as the number of ground users increases. This is because multiuser diversity gain is achieved, which, in turn, enhances the system SOP performance. The results indicate that the heterodyne detection offers the system better performance compared to the IMDD detection.
The SPSC performance of the system under the different shadowing fading effect is illustrated in Figure 9 under the heterodyne detection at HAP. It can be observed that the system becomes degraded as the shadowing fading effect get severer. The results also demonstrate that the analytical results perfectly agree with the simulation results, which justify the accuracy of the derived expressions.

6. Conclusions

The secrecy performance of a mixed FSO/RF integrated satellite-HAP relaying networks in the presence of an eavesdropper is evaluated. The FSO and RF links are model as Gamma-Gamma distribution and Rican fading distribution, respectively. The exact closed-form expression of the system equivalent end-to-end CDF is derived under different detection schemes. Based on this, the analytical closed form expressions of the concerned system COP, SOP, and SPSC are determined. To obtain more insight about the system performance, the asymptotic analysis of the system SOP is provided. The results shows that the pointing errors, atmospheric turbulence, and shadowing severity significantly degrade the system performance. In addition, the results indicate that heterodyne detection offers the system better performance compared to IM/DD detection. The future work of this paper could be tailored toward the non-orthogonal multiple access (NOMA) for better system spectrum efficiency.

Author Contributions

K.O.O. and P.A.O.; methodology, K.O.O.; software, K.O.O.; validation, K.O.O. and P.A.O.; formal analysis, K.O.O.; writing—original draft preparation, K.O.O.; writing—review and editing, K.O.O. and P.A.O.; supervision, P.A.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Asymptotic Analysis of System SOP

In this section, the asymptotic analysis of the system SOP is derived. Particularly, when γ ¯ r , the asymptotic expression of the SOP can be determined from (31) as:
S O P = S O P 1 + S O P 2 .
Thus, the asymptotic expression of the first term of (A1) can be obtained by applying [47,48] as follows:
S O P 1 = ξ 2 r α + β 2 e x p ( K ) 2 π r 1 ( Γ ( α ) Γ ( α ) ) l = 0 t = 1 3 r K l ( l ! ) 2 ( ( α β ) r Θ η r 2 r γ ¯ r ) ψ 4 , t ,
where
= g = 1 g t 3 r Γ ( ψ 3 , j ψ 3 , t ) t = 1 2 Γ ( 1 ψ 3 , j + ψ 4 , t ) t = 3 r + 2 Γ ( ψ 3 , j ψ 4 , t ) t = 3 r + 1 3 r + 1 Γ ( 1 ψ 3 , j + ψ 4 , t ) .
In addition, the asymptotic expression of the second term of (A1) can be derived by applying [45] as:
S O P 2 = ω η e x p ( K ) l = 0 q 1 = 0 m h 1 q N u = 0 m h 1 ( K η ) l Ξ ( N u ) Δ H λ ( l ! ) 2 ( 2 π i ) 2 L 1   Γ ( λ + s ) Γ ( s ) Γ ( s + 1 ) Γ ( 1 s ) ( Δ H C ) s d s × L 2   Γ ( t r s ) Γ ( ξ 2 + t ) Γ ( α + t ) Γ ( β + t ) Γ ( l t r + 1 ) Γ ( ξ 2 + 1 + t ) Γ ( t r + 1 ) ( Θ 1 / r M ( η γ ¯ r ) 1 / r ) t d t .
By using the identity stated in Reference [44] (Equation (1.2)), (A3) can be further expressed as:
S O P 2 = ω η e x p ( K ) l = 0 q 1 = 0 m h 1 q N u = 0 m h 1 ( K η ) l Ξ ( N u ) Δ H λ ( l ! ) 2 ( 2 π i ) L 1   Γ ( λ + s ) Γ ( s ) Γ ( s + 1 ) Γ ( 1 s ) ( Δ H C ) s × H 3 , 4 4 , 1 ( Θ 1 / r M ( η γ ¯ r ) 1 / r | ( l , 1 / r ) ( ξ 2 + 1 , 1 ) ( 1 , 1 / r ) ( s , 1 / r )   ( ξ 2 , 1 )   ( α , 1 ) ( β , 1 ) ) d s .
According to Reference [49], the asymptotic value of the Fox’s H-function in (A4) can be expressed as the residue of the closest pole to the left of the integration path l when the argument tends to zero. Therefore, by utilizing the identity detailed in Reference [50] (Equation (1.8.4)), (A4) can be expressed as:
S O P 2 = ω η e x p ( K ) l = 0 q 1 = 0 m h 1 q N u = 0 m h 1 r ( K η ) l Ξ ( N u ) Δ H λ ( l ! ) 2 ( 2 π i ) ( ( η γ ¯ r ) 1 / r Θ 1 / r M ) r × L 1   Γ ( λ + s ) Γ ( s ) Γ ( s + 1 ) Γ ( 1 s ) Γ ( ξ 2 + s r ) Γ ( α + t ) Γ ( β + s r ) Γ ( 1 + l s ) Γ ( ξ 2 + 1 + s r ) Γ ( s + 1 ) ( Δ H C Θ 1 / r M ( η γ ¯ r ) 1 / r ) s d s .
By applying the identity stated in Reference [44] (Equation (1.2)), (A5) can be solved as:
S O P 2 = ω η e x p ( K ) l = 0 q 1 = 0 m h 1 q N u = 0 m h 1 r ( K η ) l Ξ ( N u ) Δ H λ ( l ! ) 2 ( 2 π i ) ( ( η γ ¯ r ) 1 / r Θ 1 / r M ) r × H 4 , 6 5 , 2 ( Θ 1 / r M ( η γ ¯ r ) 1 / r | ( 1 , 1 ) ( l , 1 / r ) ( ξ 2 + 1 , 1 ) ( 1 , 1 / r ) ( λ , 1 ) ( 1 , 1 ) ( ξ 2 , r )   ( α , r ) ( β , r ) ( 0 , 1 ) ) .
By invoking (A2) and (A6) into (A1), the asymptotic expression for the system SOP can then be obtained as (32).

References

  1. Wang, W.; Gao, W.; Bai, X.; Peng, T.; Chuai, G.; Wang, W. A framework of wireless emergency communications based on relaying and cognitive radio. In Proceedings of the 2007 IEEE 18th International Symposium on Personal, Indoor and Mobile Radio Communications, Athens, Greece, 3–7 September 2007; IEEE: Piscataway, NJ, USA; pp. 1–5. [Google Scholar]
  2. Zhu, L.-p.; Yan, X.; Zhu, Y.-S. High altitude platform-based two-hop relaying emergency communications schemes. In Proceedings of the 2009 5th International Conference on Wireless Communications, Networking and Mobile Computing, Beijing, China, 24–26 September 2009; IEEE: Piscataway, NJ, USA; pp. 1–4. [Google Scholar]
  3. Berioli, M.; Courville, N.; Werner, M. Emergency communications over satellite: The WISECOM approach. In Proceedings of the 2007 16th IST Mobile and Wireless Communications Summit, Budapest, Hungary, 1–5 July 2007; IEEE: Piscataway, NJ, USA; pp. 1–5. [Google Scholar]
  4. Shah, S.; Siddharth, M.; Vishwakarma, N.; Swaminathan, R.; Madhukumar, A. Adaptive-Combining-based Hybrid FSO/RF Satellite Communication with and without HAPS. IEEE Access 2021, 9, 81492–81511. [Google Scholar] [CrossRef]
  5. Ndjiongue, A.R.; Ngatched, T.; Dobre, O.A.; Armada, A.G.; Haas, H. Performance analysis of RIS-based nT-FSO link over GG turbulence with pointing errors. arXiv 2021, arXiv:2102.03654. [Google Scholar]
  6. Odeyemi, K.O.; Owolawi, P.A.; Olakanmi, O.O. Secrecy performance of cognitive underlay hybrid RF/FSO system under pointing errors and link blockage impairments. Opt. Quantum Electron. 2020, 52, 1–16. [Google Scholar] [CrossRef]
  7. Odeyemi, K.O.; Owolawi, P.A.; Srivastava, V.M. Optical spatial modulation over gamma–gamma turbulence and pointing error induced fading channels. Optik 2017, 147, 214–223. [Google Scholar] [CrossRef]
  8. Petkovic, M.I.; Djordjevic, G.T. Impact of temporary link blockage on ergodic capacity of FSO system. ETRI J. 2018, 40, 330–336. [Google Scholar] [CrossRef]
  9. Nguyen, N.T.; Vu, M.Q.; Pham, H.T.; Dang, B.H.; Ngoc, T.D. Performance enhancement of HAP-based relaying M-PPM FSO system using spatial diversity and heterodyne detection receiver. J. Opt. Commun. 2021, 42, 111–120. [Google Scholar] [CrossRef]
  10. Vu, M.Q.; Nguyen, N.T.; Pham, H.T.; Dang, N.T. Performance enhancement of LEO-to-ground FSO systems using All-optical HAP-based relaying. Phys. Commun. 2018, 31, 218–229. [Google Scholar] [CrossRef]
  11. Swaminathan, R.; Sharma, S.; Vishwakarma, N.; Madhukumar, A. HAPS-based Relaying for Integrated Space-Air-Ground Networks with Hybrid FSO/RF Communication: A Performance Analysis. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 1581–1599. [Google Scholar]
  12. Kandeepan, S.; Rasheed, T.; Reisenfeld, S. Energy efficient cooperative HAP-terrestrial communication systems. In International Conference on Personal Satellite Services; Springer: Berlin/Heidelberg, Germany, 2011; pp. 151–164. [Google Scholar]
  13. Liu, X.; Lin, M.; Huang, Q.; Wang, J.; Ouyang, J. Performance analysis for multi-user integrated satellite and UAV cooperative networks. Phys. Commun. 2019, 36, 100762. [Google Scholar] [CrossRef]
  14. Odeyemi, K.O.; Owolawi, P.A. Wireless energy harvesting based asymmetric RF/FSO system with transmit antenna selection and receive diversity over M-distribution channel and non-zero boresight pointing error. Opt. Commun. 2020, 461, 125219. [Google Scholar] [CrossRef]
  15. Liang, H.; Li, Y.; Miao, M.; Gao, C.; Li, X. Analysis of selection combining hybrid FSO/RF systems considering physical layer security and interference. Opt. Commun. 2021, 497, 127146. [Google Scholar] [CrossRef]
  16. Odeyemi, K.O.; Owolawi, P.A.; Olakanmi, O.O. On secure transmission in hybrid satellite-terrestrial cooperative network with untrusted energy harvesting relay and imperfect channel estimation. Trans. Emerg. Telecommun. Technol. 2021, 32, e4206. [Google Scholar] [CrossRef]
  17. Miridakis, N.I.; Vergados, D.D.; Michalas, A. Dual-hop communication over a satellite relay and shadowed Rician channels. IEEE Trans. Veh. Technol. 2014, 64, 4031–4040. [Google Scholar] [CrossRef]
  18. Bhatnagar, M.R.; Arti, M. Performance analysis of hybrid satellite-terrestrial FSO cooperative system. EEE Photonics Technol. Lett. 2013, 25, 2197–2200. [Google Scholar] [CrossRef]
  19. Bhatnagar, M.R.; Arti, M. Performance analysis of AF based hybrid satellite-terrestrial cooperative network over generalized fading channels. IEEE Commun. Lett. 2013, 17, 1912–1915. [Google Scholar] [CrossRef]
  20. Wu, H.; Zou, Y.; Cao, W.; Chen, Z.; Tsiftsis, T.A.; Bhatnagar, M.R.; De Lamare, R.C. Impact of hardware impairments on outage performance of hybrid satellite-terrestrial relay systems. IEEE Access 2019, 7, 35103–35112. [Google Scholar] [CrossRef]
  21. Guo, K.; Lin, M.; Zhang, B.; Wang, J.B.; Wu, Y.; Zhu, W.P.; Cheng, J. Performance analysis of hybrid satellite-terrestrial cooperative networks with relay selection. IEEE Trans. Veh. Technol. 2020, 69, 9053–9067. [Google Scholar] [CrossRef]
  22. Ruan, Y.; Wang, Y.; Li, Y.; Zhang, R.; Li, T. Outage analysis of partial relay selection based on shadowing side information in hybrid satellite-terrestrial relay networks. Trans. Emerg. Telecommun. Technol. 2020, 31, e3826. [Google Scholar] [CrossRef]
  23. Lai, P.; Bai, H.; Huang, Y.; Chen, Z.; Liu, T. Performance evaluation of underlay cognitive hybrid satellite–terrestrial relay networks with relay selection scheme. IET Commun. 2019, 13, 2550–2557. [Google Scholar] [CrossRef]
  24. Guo, K.; An, K.; Zhang, B.; Huang, Y.; Guo, D. On the performance of cognitive satellite-terrestrial relay networks with channel estimation error and hardware impairments. Sensors 2018, 18, 3292. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Bankey, V.; Upadhyay, P.K.; da Costa, D.B.; Bithas, P.S.; Kanatas, A.G.; Dias, U.S. Performance analysis of multi-antenna multiuser hybrid satellite-terrestrial relay systems for mobile services delivery. IEEE Access 2018, 6, 24729–24745. [Google Scholar] [CrossRef]
  26. Liu, X.; Lin, M.; Kong, H.; Ouyang, J.; Zhu, W.-P. Performance analysis of mixed FSO-RF transmission in multiuser satellite-aerial-terrestrial networks. Opt. Commun. 2021, 496, 127141. [Google Scholar] [CrossRef]
  27. Singya, P.K.; Alouini, M.-S. Performance of UAV assisted Multiuser Terrestrial-Satellite Communication System over Mixed FSO/RF Channels. arXiv 2021, arXiv:2109.05762. [Google Scholar] [CrossRef]
  28. Kong, H.; Lin, M.; Zhu, W.-P.; Amindavar, H.; Alouini, M.-S. Multiuser scheduling for asymmetric FSO/RF links in satellite-UAV-terrestrial networks. IEEE Wirel. Commun. Lett. 2020, 9, 1235–1239. [Google Scholar] [CrossRef]
  29. Xu, R.; Guo, D.; Zhang, B.; Guo, K.; Li, C. Secrecy performance analysis for hybrid satellite terrestrial relay networks with multiple eavesdroppers. In Proceedings of the 2019 28th Wireless and Optical Communications Conference (WOCC), Beijing, China, 9–10 May 2019; IEEE: Piscataway, NJ, USA; pp. 1–5. [Google Scholar]
  30. An, K.; Liang, T.; Yan, X.; Zheng, G. On the secrecy performance of land mobile satellite communication systems. IEEE Access 2018, 6, 39606–39620. [Google Scholar] [CrossRef]
  31. Ai, Y.; Mathur, A.; Cheffena, M.; Bhatnagar, M.R.; Lei, H. Physical layer security of hybrid satellite-FSO cooperative systems. IEEE Photonics J. 2019, 11, 1–14. [Google Scholar] [CrossRef]
  32. Yahia, O.B.; Erdogan, E.; Kurt, G.K. On the Use of HAPS to Increase Secrecy Performance in Satellite Networks. arXiv 2021, arXiv:2106.08180. [Google Scholar]
  33. Li, S.; Yang, L.; da Costa, D.B.; Yu, S. Performance analysis of UAV-based mixed RF-UWOC transmission systems. IEEE Trans. Commun. 2021, 69, 5559–5572. [Google Scholar] [CrossRef]
  34. Benkhelifa, F.; Rezki, Z.; Alouini, M.-S. Low SNR capacity of FSO links over gamma-gamma atmospheric turbulence channels. IEEE Commun. Lett. 2013, 17, 1264–1267. [Google Scholar] [CrossRef] [Green Version]
  35. Sun, Q.; Zhang, Z.; Zhang, Y.; Lopez-Benitez, M.; Zhang, J. Performance analysis of dual-hop wireless systems over mixed FSO/RF fading channel. IEEE Access 2021, 9, 85529–85542. [Google Scholar] [CrossRef]
  36. Adamchik, V.; Marichev, O. The algorithm for calculating integrals of hypergeometric type functions and its realization in REDUCE system. In Proceedings of the International Symposium on Symbolic and Algebraic Computation, Tokyo, Japan, 20–24 August 1990; pp. 212–224. [Google Scholar]
  37. Shuai, H.; Guo, K.; An, K.; Zhu, S. NOMA-Based Integrated Satellite Terrestrial Networks with Relay Selection and Imperfect SIC. IEEE Access 2021, 9, 111346–111357. [Google Scholar] [CrossRef]
  38. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products; Academic Press: Cambridge, MA, USA, 2014. [Google Scholar]
  39. Prudnikov, A.P.; Brychkov, I.U.A.; Marichev, O.I. Integrals and Series: Special Functions; CRC Press: Boca Raton, FL, USA, 1986. [Google Scholar]
  40. Mamaghani, M.T.; Hong, Y. On the performance of low-altitude UAV-enabled secure AF relaying with cooperative jamming and SWIPT. IEEE Access 2019, 7, 153060–153073. [Google Scholar] [CrossRef]
  41. Yang, L.; Meng, F.; Zhang, J.; Hasna, M.O.; di Renzo, M. On the performance of RIS-assisted dual-hop UAV communication systems. IEEE Trans. Veh. Technol. 2020, 69, 10385–10390. [Google Scholar] [CrossRef]
  42. Lou, Y.; Sun, R.; Cheng, J.; Liu, S.; Zhou, F.; Qiao, G. Physical-Layer Security for Two-Hop Air-to-Underwater Communication Systems with Fixed-Gain Amplify-and-Forward Relaying. arXiv 2020, arXiv:2009.09550. [Google Scholar]
  43. Li, S.; Yang, L. Performance Analysis of Dual-Hop THz Transmission Systems over α-μ Fading Channels with Pointing Errors. arXiv 2021, arXiv:2107.13166. [Google Scholar]
  44. Mathai, A.M.; Saxena, R.K.; Haubold, H.J. The H-Function: Theory and Applications; Springer Science & Business Media: Berlin, Germany, 2009. [Google Scholar]
  45. Mittal, P.; Gupta, K. An integral involving generalized function of two variables. In Proceedings of the Indian Academy of Sciences-Section A; Springer: Cham, Switzerland, 1972; Volume 75, pp. 117–123. [Google Scholar]
  46. Ai, Y.; Mathur, A.; Lei, H.; Cheffena, M.; Ansari, I.S. Secrecy enhancement of RF backhaul system with parallel FSO communication link. Opt. Commun. 2020, 475, 126193. [Google Scholar] [CrossRef]
  47. Yang, L.; Hasna, M.O.; Gao, X. Performance of mixed RF/FSO with variable gain over generalized atmospheric turbulence channels. IEEE J. Sel. Areas Commun. 2015, 33, 1913–1924. [Google Scholar] [CrossRef]
  48. Odeyemi, K.O.; Owolawi, P.A.; Olakanmi, O.O. Performance analysis of reconfigurable intelligent surface assisted underwater optical communication system. Prog. Electromagn. Res. M 2020, 98, 101–111. [Google Scholar] [CrossRef]
  49. Alhennawi, H.R.; el Ayadi, M.M.; Ismail, M.H.; Mourad, H.-A.M. Closed-form exact and asymptotic expressions for the symbol error rate and capacity of the $ H $-function fading channel. IEEE Trans. Veh. Technol. 2015, 65, 1957–1974. [Google Scholar] [CrossRef]
  50. Kilbas, A.A. H-Transforms: Theory and Applications; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
Figure 1. A mixed FSO/RF integrated satellite-HAP relaying network model.
Figure 1. A mixed FSO/RF integrated satellite-HAP relaying network model.
Photonics 09 00032 g001
Figure 2. COP performance of the system under atmospheric conditions for different detection scheme at γ ¯ e = 5   dB , γ ¯ 1 = 35   dB , γ ¯ 2 = 15   dB .
Figure 2. COP performance of the system under atmospheric conditions for different detection scheme at γ ¯ e = 5   dB , γ ¯ 1 = 35   dB , γ ¯ 2 = 15   dB .
Photonics 09 00032 g002
Figure 3. Impact of atmospheric turbulence on the system SOP performance under different detection schemes when γ ¯ 2 = 10   dB , γ ¯ e = 5   dB , and ξ = 6.5 .
Figure 3. Impact of atmospheric turbulence on the system SOP performance under different detection schemes when γ ¯ 2 = 10   dB , γ ¯ e = 5   dB , and ξ = 6.5 .
Photonics 09 00032 g003
Figure 4. Effect of pointing errors on the system SOP under different detection schemes at γ ¯ e = 5   dB , γ ¯ 2 = 10   dB .
Figure 4. Effect of pointing errors on the system SOP under different detection schemes at γ ¯ e = 5   dB , γ ¯ 2 = 10   dB .
Photonics 09 00032 g004
Figure 5. Impact of γ ¯ e levels on the system SOP under different detection schemes when γ ¯ 2 = 10   dB and ξ = 1.1 .
Figure 5. Impact of γ ¯ e levels on the system SOP under different detection schemes when γ ¯ 2 = 10   dB and ξ = 1.1 .
Photonics 09 00032 g005
Figure 6. SOP performance of the system under Rician factors level for different detection schemes when γ ¯ e = 15   dB , γ ¯ 2 = 10   dB , and ξ = 1.1 .
Figure 6. SOP performance of the system under Rician factors level for different detection schemes when γ ¯ e = 15   dB , γ ¯ 2 = 10   dB , and ξ = 1.1 .
Photonics 09 00032 g006
Figure 7. Impart of shadowing fading severities on the system SOP performance under different detection schemes when γ ¯ 1 = 45   dB , γ ¯ e = 5   dB , and ξ = 6.5 .
Figure 7. Impart of shadowing fading severities on the system SOP performance under different detection schemes when γ ¯ 1 = 45   dB , γ ¯ e = 5   dB , and ξ = 6.5 .
Photonics 09 00032 g007
Figure 8. Effect of number of terrestrial users on the system SOP performance under the heterodyne detection for different pointing errors when γ ¯ 1 = 35   dB and γ ¯ e = 5   dB .
Figure 8. Effect of number of terrestrial users on the system SOP performance under the heterodyne detection for different pointing errors when γ ¯ 1 = 35   dB and γ ¯ e = 5   dB .
Photonics 09 00032 g008
Figure 9. SPSC performance of the system under shadowing fading severities for heterodyne detection at γ ¯ 1 = 45   dB , γ ¯ e = 5   dB , and ξ = 6.5 .
Figure 9. SPSC performance of the system under shadowing fading severities for heterodyne detection at γ ¯ 1 = 45   dB , γ ¯ e = 5   dB , and ξ = 6.5 .
Photonics 09 00032 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Odeyemi, K.O.; Owolawi, P.A. A Mixed FSO/RF Integrated Satellite-High Altitude Platform Relaying Networks for Multiple Terrestrial Users with Presence of Eavesdropper: A Secrecy Performance. Photonics 2022, 9, 32. https://0-doi-org.brum.beds.ac.uk/10.3390/photonics9010032

AMA Style

Odeyemi KO, Owolawi PA. A Mixed FSO/RF Integrated Satellite-High Altitude Platform Relaying Networks for Multiple Terrestrial Users with Presence of Eavesdropper: A Secrecy Performance. Photonics. 2022; 9(1):32. https://0-doi-org.brum.beds.ac.uk/10.3390/photonics9010032

Chicago/Turabian Style

Odeyemi, Kehinde O., and Pius A. Owolawi. 2022. "A Mixed FSO/RF Integrated Satellite-High Altitude Platform Relaying Networks for Multiple Terrestrial Users with Presence of Eavesdropper: A Secrecy Performance" Photonics 9, no. 1: 32. https://0-doi-org.brum.beds.ac.uk/10.3390/photonics9010032

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

Article Metrics

Back to TopTop