Capacity of κ-μ Shadowed Fading Channels

The ergodic capacity of fading channelsmodeledwith a κ-μ shadowed distribution is investigated to derive closed-form expressions. The κ-μ shadowed distribution is of particular interest because it contains, as special cases, other classical ones like one-side Gaussian, Rayleigh, Rician, Nakagami-m, κ-μ, and Rician shadowed distributions. The paper discusses the physical meaning of the distribution parameter variations and also their impact on the channel capacity. These results can be used to study the behavior of different channels like the ones in underwater acoustic communications, land mobile satellite systems, body centric communications, and other wireless communication applications. The analytical closed-form expression results are validated with numerical simulations.


Introduction
The propagation of radio waves or acoustic waves, used as data carriers, in realistic scenarios where either the transmitter, receiver, or surrounding objects are moving is quite complex.The amplitude and phase instantaneous values of the incident waves exhibit rapid variations according to the interaction patterns with the environment.This interaction produces phenomena like reflection, diffraction, and scattering of the waves on the object surfaces [1].All this makes the problem of modeling the communication channel behavior almost unmanageable with a deterministic approach.An alternative strategy is to formulate mathematical models that reach a reasonable similarity to the expected signal at the receiver from a statistical approach.Such models try to reproduce the received signal according to different time scales.On the one hand, at a large time scale, there is a path loss modeling that estimates the mean signal strength, which is notably influenced by the transmission distance.On the other hand, there is a fading modeling, at a shorter scale, to represent the rapid fluctuations around the mean due essentially to multipath propagation.Even for this fading, two levels of approximation are usual: a longterm fading due to the influence of macroscopic movements in the environment (associated with a channel time-varying behavior) that is slower than a short-term fading, which is caused by variations in the superposition of the different wave components associated with the transmitted signal, even in a quasistatic environment [2,3].There exist many statistical fading models and the - distribution has been used in some of them to represent the received signal level in a flexible way [4][5][6][7][8][9][10].
The present work focused on the short-term, narrowband, fading experienced in multipath propagation scenarios with a line-of-sight condition.Besides the multipath fading, shadowing of the dominant components is also investigated.There are two kinds of models that consider shadowing: in the first one, shadowing is assumed to change the whole signal power (composed of the dominant components and the scattered waves) [11][12][13][14].In the second type of models, the shadowing is considered to affect only the power of the dominant components, like the Rician shadowed model [15].This latter model was initially employed to fit experimental data from land mobile satellite channels and lately has also been applied to actual underwater acoustic communications (UAC) channels [16].
The present work is based on the - shadowed distribution model, first proposed in [17], which belongs to the second type of shadowing models and achieves a better fit to UAC channels measurements.
The paper is organized as follows: first, the system model is addressed in Section 2, explaining its physical interpretation and providing some fundamental statistics.In Section 3, the mathematical analysis of the ergodic capacity of the channel for the assumed model is presented, whereas, in Section 4, numerical simulations that support the analytical results are discussed.Finally, some conclusions are given in the last section.

System Model
The purpose of this section is to define the framework of the problem and to help the reader to follow the remainder of the paper.To this end, a brief summary of the physical model underneath the - shadowed distribution is given and the main statistical functions of this distribution are reviewed.

Physical Model.
In the physical model for the - shadowed distribution [17], the received power  can be expressed in terms of the in-phase and quadrature components of the fading signal by where  is a natural number,   and   are mutually independent Gaussian processes (with , and   and   are real numbers.Each term of the sum represents a multipath cluster and hence  is the number of clusters.The random variable   +   models the scattered components of the th cluster, which is a circularly symmetric complex Gaussian since the number of components is considered large enough so that the centrallimit theorem applies.In each cluster, the total power of the scattered components is 2 2 , while the dominant component of the th cluster is given by the complex random variable   +   , whose power is  2  +  2  as  is a power-normalized random variable.
The following are the main parameters to formulate the - shadowed distribution.
(i) −parameter: all the dominant components exhibit a common shadowing fluctuation that is represented by the random amplitude , that is, a Nakagami- random variable with shaping parameter  and [ 2 ] = 1.
(ii) -parameter: the natural number of clusters  can be replaced in (1) by the nonnegative real extension , which leads to a more general and flexible distribution.
(iii) -parameter: the ratio between the power of the dominant components and the power of the scattered waves is represented by  that is defined as  = (∑ Assume  ∼ S  (; , , ) is a - shadowed random variable with mean  and real nonnegative shaping parameters , , and . represents the instantaneous signal to noise ratio at the receiver and models the fading channel.It can be defined as  ≜ /, where  ≜ [] and  =

Fundamental Statistics.
In this part, the most important statistical functions corresponding to the - shadowed distribution are presented (see details in [17]).
The probability density function (PDF) of  is given by where 1  1 (⋅) is the confluent hypergeometric function defined in [19].
The associated cumulative density function (CDF) of  is given by where Φ 2 (⋅) is the bivariate confluent hypergeometric function defined in [19].
The moment generating function (MGF) of  is given by

Ergodic Capacity Analysis
For the physical model under consideration, the ergodic capacity of the channel (in bit/s/Hz) can be calculated [3] from which, after substituting (2), leads to In order to find the analytical solution to this equation, the Meijer G-function is going to be introduced.This is a very general function, intended to include many of the known special functions as particular cases, which is represented by a Mellin-Barnes type of contour integral (see 9.3 in [19]).In order to be consistent with (10), the variable of integration in the definition shown in [19] has been changed to −, as follows: where 0 ≤  0 ≤  0 , 0 ≤  0 ≤  0 , and the poles of Γ(  + ) must not coincide with the poles of Γ(1−  −) (with  = 1, . . .,  0 and  = 1, . . .,  0 ). is a suitable closed contour in the complex plane which can be chosen among three types of integration paths.
According to the values of the parameters  and , the solutions for the ergodic capacity lead to different close-form expressions.

General Case:
The parameters in the - shadowed model can take any arbitrary values; however ( − ) cannot be a positive integer.
The Meijer G-identities corresponding to each of the special functions in ( 6) are [20] ln The double integral in (10) must satisfy the convergence theorems (see II.13 in [21]).In particular, the property 1  1 ( − , ; ) =    1  1 (, ; −) must be applied to the confluent hypergeometric function so that the argument in the corresponding G-function is positive and the integral converges.
After some algebraic manipulations, the desired closedform expression for the ergodic capacity is given by The above expression is valid for any set of parameters in the - shadowed model, except for particular relationships between  and  that correspond to singularities of the gamma function.For these cases, ( − ) ∈ Z + 0 , a different approach to determine the PDF of  must be followed.This is done in the next section to obtain a new expression for the ergodic capacity.
International Journal of Antennas and Propagation 3.2.Particular Case: (−) ∈ Z + 0 .When (−) is a positive integer, the singularities in (13) are solved by getting the PDF from the MGF of the - shadowed fading model.For this purpose, the inverse Laplace transform can be employed: It is observed in (4) that the exponent of the  function in the numerator of M  () is a natural number, so the following expression (see 5.4(4) in [22]) can be used to calculate the inverse Laplace transform: where  is the generalized Laguerre polynomial.Hence, the PDF is given by Since  − 1 may be a real number and not a natural one, the expression for the generalized Laguerre polynomial must be written as (see 8.97 in [19])    () = ( + 1)  /! 1  1 (−,  + 1; ).Taking into account that For this particular case, the ergodic capacity depends on the sum of weighted Meijer G-functions, which is a simpler expression than the one presented for the general case.It can be proved that when  tends to infinity, ( 16) matches (6) in [5], which is the average channel capacity of the - model, after taking into account the following properties: By setting  = , the MGF expression boils down to and the PDF is directly   () = (  /(  Γ())) −1  −(/) and, by means of (5), the ergodic capacity is

Numerical Results
This section presents numerical simulations to validate the proposed expressions of the ergodic capacity of a - shadowed fading channel.
For the particular cases  =  and ( − ) ∈ Z + , the Meijer G-function can be evaluated by conventional mathematical packages.However, the Meijer G-function of two variables is not readily available.In the present work, the implementation presented in [18] is followed.The EGBMGF has been extended to any values of  1 and  1 in (10) and situations where the contours for the integrals cannot be a straight line parallel to the imaginary axis of the corresponding complex plane.See Appendix for a complete MATHEMATICA implementation of the function that has been used to obtain the numerical results for the ergodic capacity in the general case.
To check the analytical expressions, numerical simulations have been carried out by generating random numbers according to the physical model and the corresponding classical distributions.However, the physical model cannot be used directly when  is not a natural number.In that situation, a sampled version of the CDF in (3) has been calculated for the selected set of physical parameter values.Afterwards, random numbers have been generated with a uniform distribution and then mapped into a - shadowed distribution by using the inverse CDF method.
In Figure 1, the ergodic capacity for three different sets of parameters of the fading model is drawn.The parameters have been selected to showcase the three types of expressions obtained for the ergodic capacity ( 13), (16), and (19).In all cases the simulations and the corresponding analytical results show perfect agreement.
Figures 2-4 highlight the influence of the parameters on the channel capacity, according to the physical model.As mentioned before, the -parameter represents the level of fluctuation due to shadowing in the dominant component and the -parameter indicates the ratio between the dominant components and the scattered waves (and models the fading channel).For large  values, more power is present in the dominant component with respect to the scattered waves, while, for large , the power of the dominant component is more stable and therefore there is less shadowing.In Figures 2-4, the sensitivity of the ergodic capacity on the parameter is shown at different values of  (hence different dominant components importance).For the remainder of the paper, the parameter  has been set to 3, which corresponds to a situation where the scattered waves arrive in the receiver organized in three clusters.
In Figure 2, there is a strong dominant component ( = 5) and the influence of  on the ergodic capacity is very significant.This dependency on  increases with  (the average received signal to noise ratio).For signal to noise ratios above  15-20 dB, the ergodic capacity exhibits a quasilinear behavior, with a constant difference between the minimum capacity (at  = 0.5) and maximum (at  = 100) of approximately 0.8 bit/s/Hz.
In Figure 3, the power of the dominant component is the same as the scattered waves ( = 1) and hence  has a lower influence and there is less than 0.25 bit/s/Hz difference between the capacities at lowest and highest value of .Finally, in Figure 4, the dominant component is very weak compared to the scattered components ( = 0.2) and the impact of  on the ergodic capacity is negligible.
Table 1 summarizes the ergodic capacity and its dependency on both  and  parameters.It can be seen that, under a weak shadowing condition (large ), the ergodic capacity increases with , while, under a strong shadowing condition (low ), the ergodic capacity decreases with  as the contribution of the scattered waves is advantageous.

Conclusion
In the paper, an exact closed-form expression is presented for the ergodic capacity of communications channels that exhibit a fading according to the - shadowed distribution.These communication channels can be found in systems where a line-of-sight component experiences shadowing, for example, in land mobile satellite systems, underwater acoustic communications, or body centric communications.The derived analytical expressions have been validated against numerical simulations for different capacity curves.The influence of parameters on the shadowing level has been also investigated.The results contained in this work can be used for future evaluation of channels communication performance on different applications.

1 Figure 1 :
Figure 1: Representation of mean capacity values for different kinds of analytical solutions.The marks correspond to results obtained by numerical simulation.

Figure 2 :
Figure 2: Influence of shadowing -parameter on the mean capacity values when the dominant component is more powerful than the scattered waves ( = 5).The value of -parameter is set to 3 in all curves.The marks correspond to results obtained by numerical simulation.
to noise ratio (dB) Ergodic capacity per unit bandwidth (b/s/Hz)

Figure 3 :
Figure 3: Influence of shadowing -parameter on the mean capacity values when the dominant component has the same power as the scattered waves ( = 1).The value of -parameter is set to 3 in all curves.The marks correspond to results obtained by numerical simulation.

Figure 4 :
Figure 4: Influence of shadowing -parameter on the mean capacity values when the dominant component is less powerful than the scattered waves ( = 0.2).The value of -parameter is set to 3 in all curves.The marks correspond to results obtained by numerical simulation.