1. Introduction
The use of satellite communication systems nowadays is widespread in many diverse applications, such as navigation, mobile communication, broadcasting and disaster relief. Therefore, their proper working in the above diverse practical scenarios is very important for both the users and service providers. One highly probable event in land mobile satellite communications is the difficulty in maintaining line-of-sight (LOS) communications [
1,
2], because of the following propagation impairments: the blocking of signals caused by large obstacles (shadowing), severe attenuation (path loss) and the multipath channel exhibiting frequency selective fading may cause intersymbol interference [
3,
4]. This situation is also commonly called masking when LOS is no longer available between satellite and terrestrial users, and it severely affects the indoor users or in the case of low satellite elevation angles.
A number of works have studied the performance of hybrid satellite-terrestrial system (HSTS) [
1,
2] in the presence of the masking effect, e.g., [
5,
6,
7,
8,
9,
10]. The average symbol error rate (SER) of the fixed gain amplify-and-forward (AF) hybrid satellite terrestrial-system with generalized fading channels was derived in [
6]. The analysis of cooperative relaying strategies combined with the delay diversity technique in a digital video broadcasting-satellite services to handhelds (DVB-SH)-compliant hybrid satellite/terrestrial network, for reducing the impairments caused by the masking, was given in [
7]. A cooperative diversity scheme for mobile satellite multimedia broadcasting systems, utilizing space-time block codes and rate-compatible turbo codes to achieve diversity gains and additional coding gains, respectively, was analyzed in [
8]. A hybrid satellite-terrestrial networking approach, where land mobile users have both satellite communication and cooperative-networking capabilities, to solve the masking problem was proposed in [
9]. The use of different cooperative transmission schemes for the delivery of satellite services, in which mobile terminals forward the received signal from the satellite, was studied in [
10]. The outage performance of an HSTS was analytically evaluated in [
11]. The performance of a selective decode-and-forward (DF) protocol-based HSTS, employing the selection of the best relay terminal, was investigated in terms of outage probability by [
12]. The symbol error probability of the hybrid/integrated satellite-terrestrial cooperative network described in [
12] was evaluated in [
13]. In [
14], the authors considered the performance of AF-based relaying in a hybrid satellite-terrestrial free space optical (FSO) cooperative link, where a masked destination node was made to receive the relayed signals through an FSO link. The use of multiple antenna techniques in HSTS with masked ground node was considered in [
15].
All of the above cited papers have contributed to our understanding of the performance analysis of HSTS; however, they have focused on the ideal case without co-channel interference (CCI). The assumption of no CCI is impractical nowadays due to the deployment of many wireless standards in the same spatial location and the increased practice of reusing the spectrum resource in traditional wireless networks [
16,
17,
18]. The effect of CCI on the average SER of a particular HSTS was investigated in [
19], where the DF protocol was used at the terrestrial relay, and the destination was corrupted by multiple Rayleigh faded interferers. The LOS path between satellite and destination node was assumed to be available in [
19].
As is evident from the above literature survey, with the exception of [
19], the effect of interference on the performance characteristics of dual-hop HSTS has not been investigated before in the open literature. However, we notice that the work in [
19] was limited in the sense of practical deployment, due to the following reasons: DF was used at the relay instead of more general AF; a direct path between the satellite and destination was present (which is uncommon in urban canyons and indoor environments); and Rayleigh fading was used instead of generalized Nakagami-
m fading. A more recent and practically oriented study was conducted in [
20], where exact performance analysis of AF-based HSTS with co-channel interference was presented, considering multiple independent interferers and the effects of different channel parameters in a network. However, we note that the cumbersome performance analysis framework of [
20] can be simplified with reasonable system model simplifications without any loss in practical insight, e.g., by considering single dominant co-channel interferer. Furthermore, extensions to other digital modulation techniques are possible in addition to M-ary phase shift keying (M-PSK) studied therein, such as commonly-used M-ary pulse amplitude modulation (M-PAM) and M-ary quadrature amplitude modulation (M-QAM). For M-QAM with chunk-based resource allocation in orthogonal frequency-division multiple access (OFDMA) systems, the papers [
3,
4] have given its performance in Rayleigh fading. Furthermore, the mathematical expressions for a simple, yet accurate high signal-to-noise ratio (SNR) approximation (asymptotic SER) of the exact SER, not provided in [
20], could also be provided. As multiple co-channel-independent interferers were considered in [
20], therefore, the resulting performance analysis developed was too cumbersome and tedious for the purpose of system-level analysis and simulations.
Specifically, we point out that computationally-simple and holistic performance characterization of three node AF-based HSTS for generalized fading channels, with noise and CCI at the relay, is not reported in the open literature despite the presence of some loosely-related studies. Therefore, in this work, we deal with the many open problems/issues (explained in the last paragraph) regarding the performance analyses of HSTS and develop the holistic and computationally-tractable performance analysis framework, which also addresses all of the shortcomings in [
19,
20] explained previously. As a result, the novelty and main contributions of this work for HSTS can be listed as follows: (1) We analyze the average SER of a digitally-modulated dual-hop fixed gain AF-based relay network with interference and noise at the relay, while the destination only experiences noise. We derive the average SER for the following constellations: M-PSK, M-QAM and M-PAM; (2) We consider a network where a direct connection between source node (satellite) and destination node (terrestrial receiver) is absent, so a terrestrial relay forwards the source symbol to the destination; (3) We assume generalized fading channels where the source-relay link follows the shadowed Rice LMS model [
21]; and the channels of the relay-destination and interferer-relay links follow the Nakagami-
m model; (4) Using the moment generating function (MGF)-based approach [
22], the exact MGF of the proposed HSTS is derived, based on the derived MGF, the average SER of the considered network is given; (5) Finally, and importantly, we develop a simple (in terms of computational complexity), yet appropriately accurate high SNR approximation (asymptotic SER) for the exact SER, which provides important insight in the high SNR regime. Extensive numerical results in terms of the average SER and the asymptotic SER of the considered system model are provided in this paper. According to the authors’ best knowledge, this practically significant performance characterization of HSTS with CCI over generalized fading channels, which is analytically holistic while being of low computational complexity, shows the novelty of the presented work. Furthermore, for the considered HSTS, all of the mathematical expressions derived and the methods used, as well as the corresponding analyses presented here are reported in the literature for the first time and, so, show the novelty of this study.
The remainder of this article is organized as follows.
Section 2 gives the detailed description of the considered dual-hop relay network with co-channel interference.
Section 3 develops the performance analysis framework of the proposed system model.
Section 4 derives the computationally-efficient asymptotic average SER.
Section 5 presents the detailed numerical results. Finally,
Section 6 concludes this paper.
2. System Model
We consider an HSTS, where a satellite communicates with a destination node at the ground with the help of a relay node located at the ground, as shown in
Figure 1. Therefore, there is no direct connection between the satellite and destination node. The overall communication is divided into two orthogonal phases. In the first phase, the satellite transmits its signal to the relay. At the relay, the received signal in the presence of a single dominant interferer will be:
where
is the channel gain between the satellite and the relay;
x is the satellite’s transmitted symbol with
power;
is the channel gain between the interferer and the relay;
y is the interferer’s transmitted symbol with
power;
is the zero-mean additive white Gaussian noise (AWGN) at the relay with
variance.
In the second phase, the relay multiplies the received signal
with a fixed gain
and forwards the amplified signal to the destination. The received signal at the destination, by using Equation (
1), is given by:
where
is the channel gain between the relay and the destination and
is the zero-mean AWGN at the destination with
variance.
The satellite-relay link is assumed to follow the shadowed Rice fading channel with the following probability density function (PDF) [
21]:
where
,
,
is the average power of the LOS component,
is the average power of the multipath component,
is the confluent hypergeometric function [
23] and
is the Nakagami parameter. The relay-destination channel is assumed to follow the Nakagami-
m distribution; hence,
follows the Gamma distribution [
6] as:
where
,
; and
and
denote the shape and scale parameters, respectively, of the relay-destination channel. The interferer-relay channel is also assumed to follow the Nakagami-
m distribution; hence,
follows the Gamma distribution [
6] as:
where
,
; and
and
denote the shape and scale parameters, respectively, of the interferer-relay channel.
Here, we comment on the choice of the selected channel models for different links, i.e., the shadowed Rice and the Nakagami-
m. The satellite-relay link (satellite-to-terrestrial wireless communication channel) is assumed to follow the shadowed Rice model, which is an example of the LMS channel, while the relay-destination link (terrestrial wireless communication channel) and interferer-relay link (terrestrial wireless communication channel) are assumed to follow the Nakagami-
m model. The satellite-to-terrestrial propagation channels are inherently different from the terrestrial propagation channels because of the differences in the following underlying properties of the wireless channels: multipath fading (diffraction, reflection and scattering of the transmitted signal), LOS obstruction, path loss and shadowing. The above listed factors affect the signal level, while the channel additionally causes temporal dispersion and the Doppler shifts. The resulting first- and second-order statistics of the models are different due to these two different propagation environments. Therefore, we have selected the most widely-used channel models in their respective categories, i.e., Nakagami-
m for the traditional wireless channel and shadowed Rice for the LMS channel. The evaluation of the distribution of the sum of squared shadowed-Rice random variables and its application to the LMS channel was carried out in [
24]. The performance evaluation of an energy detector under multi-path fading/shadowing effects in a Gamma-shadowed Rician fading condition was done in [
25].
Thus, we can see that the HSTS model under consideration alleviates the need of statistics for the sum of squared Nakagami-
m random variables, required for multiple independent interferers [
20]. As previously published results are either in the form of infinite sums or higher order derivatives of the fading parameter
m, this makes the resulting modeling and analysis impractical and quite difficult to realize computationally [
20]. Our work captures the essence of HSTS in [
20] with the dominant single interferer, reduces the resulting approximations assumed therein and extends it to multiple dimensions.
4. Asymptotic SER
In this section, we present an accurate approximation for the average SER, which provides a reasonable insight into the performance of the system under consideration in the high SNR regime. Let us assume that when SNR takes a very large value, then
is approximated by its closest integer [
6,
21]; therefore, by using the Binomial expansion in Equation (
11), we have:
By substituting Equation (
4) in Equation (
12), we get:
By substituting Equation (
29) in Equation (
30) and after interchanging summation and integration, we get:
After multiplication and the simplification of terms inside the integral in Equation (
31), we obtain:
By using [
28] (Equation (2.3.6.9)), we can solve the integrals in Equation (
32) and get:
where
is defined in [
23] (Equation (9.210.2)). By substituting Equations (
5) and (
33) in Equation (
20), using the fact that
takes very large value and employing [
23] (Equation (9.211.4)), we finally get the MGF for the high SNR case as:
where:
and:
By considering the terms corresponding to
in Equation (
34), the asymptotic average SER of the considered HSTS employing M-PSK/M-QAM/M-PAM, can finally be obtained from Equation (
34) and expressions for
in
Section 3.2.
5. Numerical Results
This section presents the analytical and simulated results of the considered HSTS model using M-PSK, M-QAM and M-PAM modulation formats over generalized fading channels. We demonstrate the usefulness of the expressions derived in
Section 3 and
Section 4 using numerical examples and study the effects of interference on the system’s performance. The simulated results are obtained by generating
channel realizations for the M-PSK, M-QAM and M-PAM symbols. It is assumed that relay-destination and interferer-relay channels follow the Nakagami-
m fading with parameters adopted from [
29]. The satellite-relay LMS channel is changed according to different practical shadowing conditions. The parameters of the shadowed Rice LMS model are shown in
Table 1; see
Section 2 for more details on these channel parameters.
Figure 2 shows the average SER versus SNR of the considered HSTS, for infrequent light shadowing (in satellite-relay LMS channel), with multiple values of CCI (−5 dB, 0 dB and 5 dB) using different M-PSK modulation schemes: BPSK, QPSK and 8-PSK. It is assumed that
; and on the x-axis of
Figure 2, SNR denotes
; in the plots, CCI represents
. We consider the situation when relay is interfered by a single dominant interferer; we assume that the interferer can operate in a range of different transmit power levels. The value of total CCI
experienced by the source to relay link is varied in increasing order as: −5 dB, 0 dB and 5 dB. An arrow shown in
Figure 2 traversing the plots represents the direction of increasing CCI, i.e., from −5 dB to +5 dB. The analytical average SER for M-PSK is plotted by using Equation (
28), and that for M-QAM and M-PAM is evaluated by employing the corresponding expressions given in
Section 3.
Figure 3 shows the average SER versus SNR of the considered HSTS, for average shadowing (in satellite-relay LMS channel), with multiple values of CCI (−5 dB, 0 dB and 5 dB) using different M-PSK modulation schemes: BPSK, QPSK and 8-PSK.
Figure 4 shows the average SER versus SNR of the considered HSTS, for frequent heavy shadowing (satellite-relay LMS channel), with multiple values of CCI (−5 dB, 0 dB and 5 dB) using different M-PSK modulation schemes: BPSK, QPSK and 8-PSK.
Figure 5 shows the average SER versus SNR of the considered HSTS, for infrequent light shadowing (in satellite-relay LMS channel), with multiple values of CCI (−5 dB, 0 dB and 5 dB) using different M-QAM modulation schemes: 8-QAM and 16-QAM.
Figure 6 shows the average SER versus SNR of the considered HSTS, for average shadowing (satellite-relay LMS channel), with multiple values of CCI (−5 dB, 0 dB and 5 dB) using different M-QAM modulation schemes: 8-QAM and 16-QAM.
Figure 7 shows the average SER versus SNR of the considered HSTS, for infrequent light shadowing (in satellite-relay LMS channel), with multiple values of CCI (−5 dB, 0 dB and 5 dB) using different M-PAM modulation schemes: 4-PAM and 8-PAM. The above discussion about the values of network and interference parameters for
Figure 2 is also applicable for
Figure 3,
Figure 4,
Figure 5,
Figure 6 and
Figure 7. We observe from
Figure 3,
Figure 4,
Figure 5,
Figure 6 and
Figure 7 that the simulated SER very closely follows the analytical SER for all shadowing conditions of the LMS channel and modulations considered in the figures; indicating the correctness of the approximations taken in deriving the performance analysis framework of
Section 3.
As we can notice from
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6 and
Figure 7, when CCI at the relay increases from −5 dB to +5 dB, there is clearly a notable increase in average SER of the considered HSTS for a given modulation scheme. We see that the increase in SER is more prominent for higher-order modulation, such as 8-PSK/16-QAM, than that for the lower-order modulation scheme of BPSK/8-QAM/4-PAM. This can be seen from
Figure 2,
Figure 3 and
Figure 4, e.g., by comparing the plots for 8-PSK and BPSK for the same value of given CCI. The same line of reasoning is also applicable for comparing combinations of 8-PSK/QPSK and QPSK/BPSK modulations. When source-relay LMS channel suffers an increase in the amount of shadowing, as shown by the sequence of
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6 and
Figure 7, respectively, we notice that the average SER of the HSTS also increases correspondingly. The reader can view the effect of shadowing (in the source-relay LMS channel) on the considered system by comparing the curves for particular modulation (M-PSK) with the same given CCI from
Figure 2,
Figure 3 and
Figure 4, and also, the same observation is valid for M-QAM from
Figure 5 and
Figure 6. We also comment here about the computational complexity of the expression for average SER in Equation (
28), since it contains multiple integrals. We tested in MATLAB that for different SNRs (dB), e.g., 10, 20, 30; a modern personal computer takes approximately 0.1 s to calculate Equation (
28).
The difference between exact average SER and asymptotic average SER (high SNR approximation) is plotted in
Figure 8. The analysis is performed for BPSK-modulated HSTS with different values of CCI over the average shadowed LMS source-relay channel. It can be seen clearly from
Figure 8 that with the increasing value of SNR, the gap between the exact SER curve and the asymptotic SER curve keeps on decreasing. The difference between the exact and asymptotic values is very small in the high SNR region of 25–30 dB, and when SNR is over 30 dB, the magnitude of the difference between exact and asymptotic values is unnoticeable. Another example of the difference between exact average SER and asymptotic average SER (high SNR approximation) is presented in
Figure 9. The analysis is done for QPSK-modulated HSTS with different values of CCI over average shadowed LMS source-relay channel. It can be noticed clearly from
Figure 9 that with the increasing value of SNR, the gap between the exact SER curve and asymptotic SER curve keeps on decreasing. The difference between the exact and asymptotic values is very small in the high SNR region of 25–30 dB, and when SNR is over 30 dB, the magnitude of the difference between exact and asymptotic values is unnoticeable. As discussed above,
Figure 8 and
Figure 9 indicate the correctness of the assumptions and approximations taken for deriving the asymptotic SER. The so-derived approximation is a useful tool for substituting the computationally-intricate formulation inherent to the exact expression of Equation (
28).