The Slashed-Rayleigh Fading Channel Distribution

. We propose an alternative distribution for modelling fading-shadowing wireless channels. This distribution presents certain advantages over the Rayleigh-lognormal distribution and the K distribution and has proved useful in the setting described. We obtain closed-form expressions for the average channel capacity and for the average bit error rate of differential phase-shift keying and of minimum shift keying when the new distribution is used. This distribution can be obtained exactly as the sum of mutual independent Gaussian stochastic processes, because it must represent the simulation of the fading channel; that is, it simulates the signal envelope. Finally, we describe practical applications of this distribution, comparing it with the Rayleigh-lognormal and K distributions.


Introduction
Mobile communications systems must address various challenges that seriously degrade signal strength. A major problem in this respect is that of the fading signal, i.e., interference to the multiple scattered radio paths between the base station and the vicinity of the mobile receptor. When this occurs, the received signal exhibits rapid signal level fluctuations which are generally Rayleigh distributed.
The multipath channel resembles a physical communication channel and is characterised by bandwidth and gain. For the case of a mobile radio channel with constant gain and a linear phase response over the bandwidth that is greater than the bandwidth of the transmitted signal, the signal received at the terminal will undergo flat fading. This is the most common type of fading and the most widely studied (Rappaport [1]). In this situation, the strength of the received signal will oscillate rapidly due to multipath effects. Figure 1 illustrates a simplified fading model between a stationary source (emitter) and a mobile vehicle (for a complete description of fading modelling, see Barts and Stutzman [2], Rappaport [1], Borhani et al. [3], and Lopez-Fernandez et al. [4]). Several components are involved; first, in the case of clear line of sight between the stationary source (emitter) and the mobile vehicle, no scattering mechanism is involved, although Doppler effects will be considered. However, a multipath or diffuse component (phase-incoherent wave) is created by the multiple (random) reflections (scattering processes) of the signal from scattered elements such as mountains and buildings. This component has little directivity and its magnitude is assumed to be Rayleigh distributed while its phase is uniformly distributed.
The specular component is a phase-coherent groundreflected wave that is related to points close to where the receptor (i.e., the vehicle) is dynamically located. This component is responsible for strong fades, with an amplitude comparable to that of the direct component, although its phase is opposite Suh [5]. Figure 1 also shows the case of blocked line of sight between the stationary source and the mobile vehicle. In this case, the diffuse component is as before (and can be modelled by a Rayleigh distribution), but there is a new element, the shadowed direct component, produced by the scattering of the signal by the leaves, branches, and limbs of trees and by vegetation in general. In consequence, the signal is attenuated, to a degree that depends on the length of its path through the scattering element. The fading produced by this process is called fading-shadowing and can be modelled by the Rayleigh-lognormal (RLN) distribution (Hansen and Meno [6]), although this has a complicated integral form, or by the K distribution Abdi and Kaveh [7]. The latter distribution may be viewed as a Rayleigh distribution with a gamma distribution and it has a simpler form than the RLN.
This study examines fading-shadowing mechanisms in particular. However, as our proposal also contains the Rayleigh as a general case, it also applies to modelling the diffuse component.
To achieve a probability distribution that will efficiently model fading effects, a distribution should be expressed through simple mathematical expressions and subsume the Rayleigh distribution. Given this, parameter estimation is straightforward and second-order statistics such as fading descriptors (level crossing-rate (LCR), average fade duration (AFD), average bit error rate (BER), differential phase-shift keying (DFSK), and minimum shift keying (MSK)) can easily be computed (see, for instance, Proackis and Salehi [8], Adachi et al. [9], and Subbarayan [10], for a comprehensive description of these estimators).
The above conditions are met by the Slashed-Rayleigh (SR) distribution, a two-parameter distribution that was first proposed by Iriarte et al. [11]. However, the latter did not address the question of fading applications. As we discuss below, this distribution is well suited to modelling flat fading effects in wireless communications.
Among the benefits offered by the SR distribution, it includes the Rayleigh distribution, thus facilitating the physical modelling of multipath signal propagation through phasors (complex signal representation).
We complete the description of the SR distribution as follows. First, we simulate it by Monte Carlo analysis (i.e., taking a statistical approach). Then, as required for a fading distribution, we simulate the distribution from a summation of phasors while accounting for Doppler effects (i.e., a physical approach), by embedding the SR distribution within Clarkes's model Rappaport [1], which is known to be appropriate for flat fading.
Although fading effects are apparent in general for mobile communications, they are especially noticeable in moving vehicles, where the signal received commonly presents multipath components and not just the direct component. However, for a stationary receiver, too, if the surrounding objects are moving, Doppler shift on the multipath components will affect the quality of the transmitted signal. The models we discuss can be applied to both stationary and moving receptors.
The outline of this paper is as follows. The SR distribution is presented in Section 2. Section 3 then introduces the SR phasor, shows some simulation plots, and discusses specific mathematical measures for the new distribution, together with metrics related to modelling fading effects in wireless communication channels. In Section 4, the SR distribution is compared with other distributions that are commonly used to account for the statistics of mobile radio signals. The simulation of the SR distribution through Monte Carlo analysis and the physical modelling of the multipath fading channel are discussed in Section 5. Finally, the main conclusions drawn are summarised in Section 6.

The Slashed-Rayleigh Distribution
The Rayleigh fading model can be used to simulate the situation in which a radio signal is scattered before it arrives at the receiver due to the presence of multiple objects in the environment. According to the central limit theorem, given sufficient scatter, the channel impulse response will be well modelled as a Gaussian process irrespective of the distribution of the individual components. If there is no dominant component to the scatter, this process will have a zero mean and its phase will be uniformly distributed between 0 and 2 radians. Therefore, the envelope of the channel response will be Rayleigh distributed, with the following probability density function (pdf) where ( 2 ) = 2 is the expected value of 2 . In this case we write ∼ R( ). As Tse ([12], p. 49) pointed out, the model based on the Rayleigh fading is quite reasonable for scattering mechanisms when there are many small reflectors, but it is adopted primarily for its simplicity in typical cellular situations with a relatively small number of reflectors. For this reason some alternatives to the Rayleigh fading model have been proposed, such as the RLN distribution with pdf, described by (Hansen and Meno [6] and Stuber [13])

Mathematical Problems in Engineering
The K distribution (Abdi and Kaveh [7]), obtained by compounding a Rayleigh distribution with a gamma distribution, is similar to the RLN distribution but it has a simpler structure and its pdf admits a closed form, although due to the Bessel function the estimates of the parameters are not direct. Its pdf is given by where ] ( ) denotes the modified Bessel function of the second kind of order ] and argument and which has a complicated integral form. The main advantage of this special function is that it is included in most of the statistical software currently available, such as R, Matlab, and Mathematica (see Ruskeepaa [14]). Apart from the RLN distribution, other alternatives based on the lognormal distribution include the Rayleigh-inverse Gaussian (RIG) distribution (Karmeshu [15]) which is subject to the same restriction as the above distribution and the generalisation of the Rayleigh distribution that was recently proposed by Gómez-Déniz and Gómez [16], which overcomes some of the disadvantages of the above distributions.
The SR distribution, proposed by Iriarte et al. [11], has the following pdf: where > 0, > 0, and ( , ) is the cumulative distribution function of the gamma distribution. Observe that the pdf given in (1) is obtained from (4) when tends to ∞. Furthermore, this distribution tends to the Dirac delta function at 0 when tends to infinity. Henceforth, when a random variable follows this distribution it will be represented as ∼ ( , ).
If ∼ ( , ) with scale parameter and kurtosis parameter , then (4) can be represented as where ( ) and (0, 1) (the uniform distribution) are independent and > 0. Using the following expression, which relates the incomplete gamma function with the Kummer confluent hypergeometric function the pdf given in (4) can be rewritten as ) .
Now, from (7) and using Kummer's first theorem we have ) . (8) Figure 2 shows slopes of the pdf of the ( , ) distribution, revealing the dependency of the scale parameter and the shape parameters (the fading).
which can be used to obtain the hazard function of the random variable ∼ ( , ), given by Let ∼ ( , ). Then, for = 1, 2, . . . and > , it follows that the -th moment of the proposed distribution is given by Therefore, from (11) we obtain the mean and the variance of the distribution, which are given by respectively. Moment and maximum likelihood estimation is studied in detail in Iriarte et al. [11], where a simulation study is also performed. The pdf of the SR distribution can be written as an infinite convex sum of Nakagami distributions. To do so, observe that from the series expansion of the Kummer confluent hypergeometric function in (8) this pdf can be rewritten as Now, by performing on (13) the change of variable = √ , where ∈ N * and taking +1 = , (13) can be written as where ( ; Ω, ) = 2 which is the Nakagami distribution with parameters ∈ N * and Ω = 2 > 0. Hence Thus, we conclude that the SR distribution studied here can be written as an infinite convex sum of Nakagami distributions.

The Slashed Channel Rayleigh Phasor
In this section, we demonstrate that the SR distribution can be obtained in an exact form as the sum of mutual independent Gaussian stochastic processes, as is required in order to simulate the fading channel, i.e., the signal envelope.
It is known that Rayleigh fading envelopes can be generated from zero-mean complex Gaussian random variables. Other fading distributions (see, for instance, Yacoub et al. [17] for the Nakagami-case) and the generalised Rayleigh distribution in Gómez-Déniz and Gémez [16] are obtained in a similar manner after some mathematical considerations. Hence, in line with these precedents, we must now prove that the phase and the amplitude of a given propagating signal are distributed according to a uniform pdf in the [0, 2 ] interval and the SR distribution, respectively. Following Beckmann ([18], p. 118), consider the sum where = √ −1, the terms (the in-phase phasor) and (quadrature phasor) are independent uniformly distributed phasors (UDP), and the are all distributed identically. When is large, and assuming that is not correlated with the Φ , both and will be distributed normally with mean 0 and variance (1/2) = ∑ =1 2 . Goldsmith [19], see page 69, pointed out that under some conditions this is also true for small . Now, let (1/2) = ∑ =1 2 = −2/ , where represents the uniform distribution in (0, 1). Then, the joint distribution of and is Then, expressing (18) as polar coordinates In the following, we obtain the phase and the amplitude distributions. Proposition 1. The SR distribution satisfies the following: Mathematical Problems in Engineering 5 (ii) The amplitude distribution is given by (4).
Proof. From (19) it is straightforward to obtain the conditional distributions of and given = , which are given by Then, the unconditional distribution (independent of ) for the phase is The unconditional distribution for the amplitude is given by Now, in (22) we perform the change of variable = ( 2 /2 ) 2/ and obtain Hence, the result follows after identifying a gamma distribution within the integral.
According to result (ii) in Proposition 1, if | = ∼ ( −2/ ) and ∼ (0, 1), then ∼ ( , ), i.e., pdf (4) can be represented as a scale mixture (compound) of the Rayleigh distribution and the uniform distribution on the unit interval. Moreover, this representation of the distribution facilitates parameter estimation via the expectationmaximisation (EM) algorithm.

Measures of Interest in the Fading
Channel. The SR distribution can readily be obtained as a scale mixture of the Rayleigh distribution and the uniform distribution, which facilitates the computation of some measures of interest in the framework of the fading channel, such as the amount of fading (AF) (also known as the strength of intensity fluctuations), and the BER for DPSK and MSK when the SR distribution is employed as that of the fading channel.
For a single-input/single-output system, the amount of fading (see Abdi and Kaveh [21]), given by = var( 2 )/ 2 ( 2 ) for the ( , ) distribution proposed here is given by which varies in the interval (4 − 1, ∞). We recall that the value of AF for the Rayleigh, ( ), distribution is 4 − 1.
It is well known (see Abdi and Kaveh [7]) that, for the standard Rayleigh, ( −2/ ), the BERs for DPSK and MSK are given by Now, from (ii) in Proposition 1 we obtain, by compounding, the corresponding average BERs of DPSK and MSK for the SR distribution, as follows.

Proposition 2. The average BERs of DPSK and MSK for the SR distribution are given by
respectively. Here = / 0 , where is the transmitted energy per bit and 0 is the noise power spectral density and 2 1 ( , ; , ) is the hypergeometric function.
Proof. By applying the composite rule we obtain Now, with the change of variable = 2/ expression (29) reduces to , which is immediately identified with (28) after simple algebraic manipulation. The expression for , is obtained in a similar way.
The average channel capacity for fading channel is a useful metric, in that it provides an estimation of the information rate that the channel can support, with little probability of error. Channel capacity, (see, for instance, Li et al. [22] Singh and Rawat [23], among others), is defined as where is the received signal bandwidth. Following Li et al. [22] the Shannon capacity of the Rayleigh fading channel is given by where ( ) = ∫ ∞ 1 exp(− )/ , = /( 0 ), is the average transmit power receiving bandwidth , and the mean channel gain is √ /2.

Proposition 3. The average Shannon capacity of the SR fading channel is given by
where , and ≈ 0.577216 is Euler's constant.
Proof. Again, from the composite rule and using (32) with ≡ −2/ we obtain It is known (see Gautschi et al. [24] and Lin et al. [25]) that Thus, we have from which the result follows, after some algebraic manipulation.  Figure 3 shows the average BERs for DPSK and MSK for the RLN, K, and SR distributions, for the three sets of parameter values given in Table 1. The BERs for the RLN and RIG distributions were computed numerically for DPSK and MSK using the same parameter values for the K distribution as in Abdi and Kaveh [21], in three settings: Setting A for = 0.35, Setting B for = −0.37, and Setting C for = −0.65. In every case = 1. The corresponding parameter values for the RLN and SR distributions were computed by equating the population moments to the moments of the K distribution obtained using the above values of ( , ). We recall that the population moments of the K distribution are given by ( ) = ( 2 ) = 2 exp( 2 /2 + ), respectively. The moments of the SR distribution are given in (11). The estimated parameter values of the RLN distribution, and , differ from those obtained by Abdi and Kaveh [21] because the latter study used the approximation given in Abdi and Kaveh [7].

Comparisons
The graphs show that the SR distribution is a valuable means of predicting the BER in multipath dispersion fades. Compared with the K distribution, SR also achieves a good fit with the RLN distribution. The third scenario is conclusive in this case. We recall that the RLN distribution is commonly used in DPSK and MSK modulation schemes. It should also be noted that, for the RLN distribution, there is no closed-form expression for the average BER, which must be calculated by numerical integration methods (generally, the Gauss-Hermite method). Cygan [26] proposed an exact but complicated formula for estimating the BER in the DPSK case when the RLN distribution is used.
The analytical expressions of the SR and K distributions both include special functions, namely, the incomplete gamma function and the modified Bessel function, respectively. The two are similar and either can be used as a substitute for the RLN distribution. Therefore, the proposed distribution can be applied efficiently to capture the bleached shading aspects of the wireless channels. Figure 4 shows the pdf of the RLN, K and SR distributions for the parameters given in Settings A, B, and C. Note that the SR distribution presents a longer tail than the other distributions.
Mathematical Problems in Engineering

Comparison with the RL and K Distributions.
The distance or relative information between two probability distributions can be determined by the Kullback-Leibler divergence measure (see Hall [27], among others), which is defined as follows. Let and be probability densities on R such that is absolutely continuous with respect to (that is, ( ) = 0 implies ( ) = 0), and then the relative information or Kullback-Leibler divergence, ( ‖ ) of with respect to , is with the convention that 0/0 = 1. When is not absolutely continuous with respect to we define ( ‖ ) = ∞. A disadvantage of (38) is that the Kullback-Leibler divergence is not symmetric and therefore is not a genuine distance metric. To overcome this problem, we use the Jensen-Shanon divergence (see, for instance, Lin [28]) given by where = ( + )/2. The integrated squared error (ISE) (Bowman [29]) is given by (40) Table 2 shows the Jensen-Shannon divergence for the K and SR distributions. Here, we see that, except for Setting A, the SR distribution presents the least distance with respect to the RL distribution, which was taken as the reference distribution.

Simulating the Proposed Distribution
To properly model a fading process, a random variate must be generated according to the proposed density distribution. Moreover, the random variate must be generated at low computational cost.
To achieve this, we proceed as in Gómez-Déniz and Gómez [16], using the well-known inverse transform method. Although there exist a plethora of methods for simulating

Input:
Alpha: parameter ( ∈ R) Q: parameter ( ∈ R) A: range of data (Amplitude ∈ R) NSamples: number of samples (NSamples ∈ N) Output: Z: set of data ∼ ( , ) Begin: v ← CDF(Alpha, Q) (0 → A) ⊳ Obtain the CDF and store it within a vector for ← 1 do u ← rand ⊳ Obtain a uniform random sample ( ) ← max{V ≤ } end for End Algorithm 1: Simulating the SR distribution. a random variate, the inverse transform method is simple to understand and easy to programme. Algorithm 1 shows a simplified version of the algorithm used. Note that the whole cdf of the SR is precalculated to speed up the implementation. A Monte Carlo simulation was coded in Matlab MATLAB [30] and tested on an i7-7700HQ CPU @ 2.80 GHz (16 GB RAM), taking around 0.3 seconds to generate = 100000 random variates. This low computational cost makes the approach sufficiently effective for fading channel analysis.
Two datasets were generated through Monte Carlo simulations with 100 000 samples and then compared with the analytical model. Figure 5 shows the analytical pdf for the SR distribution and the pdf simulated for both datasets. The simulated data provide a good fit with the analytic model for both datasets. The parameters used for the SR distribution were = 0.3 and = 3 (left) and = 6 and = 5 (right). Table 3 shows the analytical mean and variance, in comparison with the values estimated from the Monte Carlo samples. The relative error for both measures and both cases is low (less than 2%).
This kind of simulation is only valid for mathematical analysis. To simulate the fading effects and provide the fading channel characteristics required, the modelling must be performed from the phasor formulation derived above and linked to physical variables (signal carrier frequency, signal sampling, speed of receiver, and Doppler effects). These questions are examined in the next section.

The SR Distribution for Modelling Fading Effects. The physical model is completed by reformulating the phasors
as described by Rappaport [1] among others, as follows: Here, is the amplitude and satisfies the condition that the ensemble average < ∑ =1 2 ) >= 1. The phase is now replaced by the term = − . The phase is the random phase uniformly distributed in [− , ], = V cos( ) is the Doppler shift, where V is the vehicle speed, = 2 / is the wave number, is the wavelength, and V is the maximum Doppler shift (in radians per second). The angle of arrival of the signal is , also distributed in [− , ], and is the number of sinusoidal waves, which, if large enough, ensures that both ( ) and ( ) can be considered as Gaussian processes. Accordingly, the SR signal envelope is given by 2 = 2 + 2 , and = 1, 2, . . . The amplitude of the envelope is then fitted to the SR( , ) amplitude.
Taking into account the above, we first performed a Monte Carlo simulation to obtain the dataset, using 15 scattered random phasors and 20 000 samples. We then compared the resulting pdf with the analytic SR pdf. This comparison is shown in Figure 6 (left column) for two sets of SR distribution parameters ( = 2, = 10 and = 0.3, = 3). As can be seen, the fit between the analytic and the physical models is reasonably good. Figure 6 (right column) also shows a set of simulated samples, for both cases. As can be seen, these data samples are mostly grouped around the mean value provided by the parameter setting. Note also that    Table 4 shows the analytical means and variances, compared with those estimated from the Monte Carlo samples. The relative errors for both measures and both cases are low (around 5%), but larger than when the data are generated directly from the cdf, because this simulation is more complex and includes physical effects (wave scattering). Figure 7 shows Monte Carlo simulated samples obtained using phasors for the SR ( = 0.3, = 3) and RLN ( = 1.029, = 0.73) distributions, with similar results (note, the RLN distribution is widely accepted for modelling fading effects).
To conclude our demonstration of how the SR distribution efficiently models fading effects, we performed a simulation of the physical channel. Algorithm 2 shows the code used. This implementation reproduces Clarke's model (mapped to the timing domain). A baseband Doppler filter was also included in the simulations, in order to generate a realistic fading spectrum, that is, to produce fading waveforms that are properly time-correlated. Figure 8 shows a simulated fading signal for the SR ( = 0.3, = 3) distribution, corresponding to a signal with a mean value amplitude of -2.2262 dB. Five Rayleigh processes were simulated to obtain the SR( , ) envelope. This simulated envelope pdf of the dataset fits well with the analytical one (not shown). In this case, the envelope includes very deep fading levels (around -40 dB) which, although rather infrequent have been reported in rapid fading in longdistance HF propagation (see M et al. [31]). The results shown are for two vehicle velocities (60 km/h and 120 km/h). The Doppler effect is significantly more apparent at the higher velocity, as expected.
In the present paper, the physical interpretation of the parameters of the SR distribution is not discussed, but it seems clear from expression (13), and Algorithm 2, that there must exist a straightforward relation between these two parameters and the average signal power strength (through the means and variances for the amplitude of scattered waves). However, more analysis of this question is needed for a complete physical interpretation.
A simulator block implementing the SR( , ) signal has been developed (coded in Matlab) and included in a more general mobile radio channel software simulator, which performs Gaussian/uniform processes to simulate the phasors and thus obtain the signal envelope distributed with SR( , ) pdf. Some typical routines are also included, such as apole Tchebicheff filter block and a simple RF combiner (for the equal gain and maximum-ratio cases), together with the possibility of simulating samples for the other commonly used fading distributions (i.e., Rayleigh, GR and RLN).

Final Comments
In this paper, we present a new distribution, the twoparameter SR distribution, which is competitive with the RLN distribution and the K distribution. The Rayleigh distribution is considered as a special case when one of the parameters tends to infinity. Parameter estimation for the new distribution is also discussed. From a mathematical standpoint, the distribution we propose has certain advantages over the RLN distribution within the framework of fading signal strength. Two methods to obtain the simulated envelope are discussed; the first is based strictly on the pdf of the distribution, and the second is based on a physical model constructed from the Rayleigh physical model, using the classical method described by Clarke. In addition, closed-form expressions are derived for the average channel capacity and the bit error rate (BER) for DPSK and MSK modulations for the proposed distribution.

Data Availability
The data used to support the findings of this work can be reproduced by referring to the mathematical formulas and the algorithms described. In addition, the codes used are available from the corresponding author upon request. These data include Mathematica and Matlab codes and therefore users will require the corresponding licenses.

Conflicts of Interest
The authors declare that they have no conflicts of interest.