## Fluctuating Line-of-Sight Fading Distribution: Statistical Characterization and Applications

∗ ∗ † ,

Thanh Luan Nguyen , Georges Kaddoum ∗ Department of Electrical Engineering, Ecole de Technologie Sup´ erieure ( ´ TS), Montr´ eal, QC, Canada ´ E † Artificial Intelligence &amp; Cyber Systems Research Center, Lebanese American University Emails: ∗ thanh-luan.nguyen.1@ens.etsmtl.ca, ∗ † , georges.kaddoum@etsmtl.ca,

Abstract -We introduce the fluctuating Line-of-Sight (fLoS) fading model, characterized by parameters K k λ , , , and Ω . The fLoS fading distribution is expressed in terms of the multivariate confluent hypergeometric functions Ψ 2 , Φ ( n ) 3 , and Φ 3 = Φ (2) 3 and encompasses well-known distributions, such as the Nakagamim , Hoyt, Rice, and Rician shadowed fading distributions as special cases. An efficient method to numerically compute the fLoS fading distribution is also addressed. Notably, for a positive integer k , the fLoS fading distribution simplifies to a finite mixture of κ µ -distributions. Additionally, we analyze the outage probability and Ergodic capacity, presenting a tailored Prony's approximation method for the latter. Numerical results are presented to show the impact of the fading parameters and verify the accuracy of the proposed approximation. Moreover, we illustrate an application of the proposed fLoS fading distribution for characterizing wireless systems affected by channel aging.

Index Terms - -κ µ fading, confluent hypergeometric functions, outage probability, channel capacity, Prony's method

## I. INTRODUCTION

T HE advent of sixth-generation (6G) wireless networks is expected to bring forth various applications for ubiquitous, high-speed connectivity [1], [2]. To name a few, the use reconfigurable intelligent surfaces (RIS) can configure radio environments for enhanced performance. Moreover, highmobility communication (HMC) ensures broadband services for users traveling at speeds of up to 1000 km/h. In addition, the utilization of terahertz (THz) bands, ranging from 0 1 . THz to 10 THz, offer transmission rates reaching multiple terabits per second. However, despite these advancements, significant challenges persist within each of these technologies. For instance, RIS-aided communication requires compelling applications that demonstrate substantial performance gains over existing technologies like mmWave and massive multipleinput-multiple-output (mMIMO) systems [3]. HMC faces issues like channel aging, where the system relies on outdated channel state information (CSI) for information processing [4]. THz communications encounter blockages and lack convincing theoretical frameworks and practical channel models [5], [6]. Therefore, to fully realize the potential of these emerging technologies and address the associated challenges, the integration of novel fading channel models that accurately simulate real-world 6G environments is essential.

The literature presents numerous generalized fading models that effectively depict the statistical behavior of mobile radio signals. Various distributions, including the Rayleigh,

Rice, Nakagamim , Hoyt (Nakagamiq ), and κ µ -distributions, describe short-term signal fluctuations [7]. Long-term variations, usually referred to as shadowing, are represented by the lognormal distribution, assuming independent attenuation from obstacles [7], [8]. However, practical scenarios introduce dependencies among obstacle attenuations due to phenomena like diffraction and refraction, challenging the independence assumption. While the lognormal model is widely accepted, alternative approaches like the Nakagamim distribution provide tractable substitute [9]. Moreover, in the context of dynamic environments, such as in vehicle-to-vehicle (V2V) communications, Gaussian, Rayleigh, and Rician distributions are favored for their ability to account for non-Gaussian shadowing with non-zero skewness in the dB scale [8], [10]. Despite these models, significant theoretical gaps persist in characterizing generalized shadowing and other Line-of-Sight (LoS) variations due to various Radio Frequency (RF) propagation phenomena, leading to the need of comprehensive models for capturing real-world signal propagation phenomena.

This paper seeks to address these gaps by providing a novel approach for characterizing the fluctuations of the dominant specular component, i.e., the LoS component, using generalized distributions. Specifically, we propose a novel fluctuating LoS (fLoS) fading distribution, where LoS fluctuations are modeled by a noncentral chi-square distribution. This fLoS fading distribution encompasses classical models, like the Hoyt, Rice, Nakagamim , Rician shadowed, and mixture of κ µ -distributions, as special cases. Our main contributions are summarized as follows 1 :

- · We introduce a novel probability density function (PDF) and cumulative distribution (CDF) for the proposed fLoS fading characterized by 4 parameters K k λ , , , and Ω .
- · We offer numerical methods for tackling the derived distributions due to their high computational complexity. Notably, for integer k , the fLoS fading distribution can be expressed as a mixture of κ µ -distributions.
- · We derive the outage probability (OP) and Ergodic capacity (EC) by introducing a Prony's method to approximate various functions and compute the EC for wireless systems.

1 Our results are reproducible with MATLAB code, uploaded at https://github.com/thanhluannguyen/fluctuatingLoS

- · We present a use case for the fLoS fading distribution and we analyze the wireless system under the influence of channel aging.

The paper is organized as follows: Section II introduces the physical model of the fLoS fading channel. Subsequent analysis of the fLoS fading distribution, including the PDF, CDF, moments, and moment generating function (MGF), is presented in Section III. Section IV presents the performance analysis and a potential application for the fLoS fading distribution. Section V details the numerical results. The paper is concluded with a summary of findings in Section VI.

## II. PHYSICAL MODEL

We consider the first-order scattering model, introduced in [11], [12], to present the received signal under the fLoS fading model with LoS fluctuations as follows:

$$S = \omega _ { 0 } \xi e ^ { j \phi _ { 0 } } + \sigma G,$$

where ω e 0 jϕ 0 is the LoS (specular) component with constant amplitude ω 0 , ϕ 0 is a random variable (RV) uniformly distributed in [0 , 2 π ) . The term σG , where G ∼ CN (0 1) , , represents the diffuse scattering, or non-LoS component, with total power σ 2 . Here, K = | ω 0 | 2 σ 2 ∈ [0 , ∞ ) , such that ω 0 = √ K √ K +1 and σ = 1 √ K +1 , is the ratio between the power of the dominant LoS component and the average power of the scattering components. In various researches, the component ξ is often assumed as the lognormal shadowing, or normal (Gaussian) shadowing in decibels. This model assumes independent attenuations from independent obstacles, i.e., scatterers and reflectors, within local clusters, common in heavy shadowing environments. Specifically, if each obstacle (reflector or scatterer) independently attenuates the signal by A m along a path, the total attenuation due to M independent obstacles is ξ = ∏ M m A m or ξ | dB = ∑ M m ξ m | dB in dB. With sufficiently large M , the total attenuation (in dB) follows a Gaussian distribution due to the Central Limit Theorem (CLT) regardless of the distribution of A m distribution [8]. Yet, in reality, the obstacles' attenuations are often correlated due diffraction and receiver body vibrations. In high frequency bands (e.g., mmWave, THz), molecular absorption loss significantly influences obstacle attenuation [2]. Consequently, the characterization of the shadowing distribution is useful when M is sufficiently small to prevent the probability of ξ = 0 approaching one. However, when M is small, the CLT cannot be applied to characterize ξ with a lognormal distribution. Despite the widely accepted lognormal assumption, there exist other tractable shadowing models [9]. In this manner, we consider ξ 2 follows a (scaled) noncentral chi-squared distribution with 2 k degrees of freedom (d.o.f.), 2 λ noncentrality parameter, and scale factor Ω 2 , denoted as ξ 2 ∼ ˜ χ 2 k (Ω; λ ) . The PDF of ξ 2 is:

$$f _ { \xi } ( \xi ) = \frac { e ^ { - \frac { \tilde { \eta } } { \Omega } - \lambda } } { \Omega } \left ( \frac { x } { \lambda \Omega } \right ) ^ { \frac { k - 1 } { 2 } } I _ { k - 1 } \left ( 2 \sqrt { \frac { \lambda x } { \Omega } } \right ), x > 0. \quad ( 2 )$$

where I v ( x ) is the v th order modified Bessel function of the first kind. It is noted that for Ω = 1 k + λ with λ → 0 , (2)

reduces to the Gamma distribution with unit mean and shape factor k , as f ξ ( x ) = k k Γ( k ) x k -1 e -kx , x &gt; 0 , and the model (1) is simplified to the well-known Rician shadowed fading, where k depicts the shadowing severity [13]. We must note here that while the term 'shadowing' is used, the fLoS model also includes variations in the direct LoS path due to other factors, such as movement or dynamic environments.

Remark 1: The noncentral chi-square distribution is closely related to the κ µ -fading power distribution with parameters (¯; γ κ, µ ) . These two distributions can be deduced from one another by substituting suitable parameters. Specifically, let F κµ (¯; γ κ, µ ) be a RV following the κ µ -fading power distribution, we have the following relations:

$$\widetilde { \chi } _ { k } ^ { 2 } ( \Omega ; \lambda ) \overset { d } { = } \mathcal { F } _ { \kappa \mu } ( \Omega ( k + \lambda ) ; \lambda / k, k ), \quad \ \ ( 3 )$$

$$\mathcal { F } _ { \kappa \mu } ( \bar { \gamma } ; \kappa, \mu ) \stackrel { d } { = } \widetilde { \chi } _ { \mu } ^ { 2 } ( \bar { \gamma } [ \mu ( 1 + \kappa ) ] ^ { - 1 } ; \mu \kappa ). \quad \ \ ( 4 )$$

## III. STATISTICAL CHARACTERIZATION

In this section, we derive exact analytical formulas for the Moment Generating Function (MGF), the PDF, the CDF, and raw moments of the instantaneous SNR under fLoS fading conditions. For normalization purposes, let Ω = 1 k + λ , so that E {| S | 2 } = 1 and γ = ¯ γ S | | 2 , where ¯ γ is the average SNR.

## A. Probability Functions

Lemma 1: Considering γ ∼ F LoS (¯; γ K k,λ, ; Ω) with nonnegative real values K k λ , , , and Ω , to characterize the SNR under fLoS fading, where ¯ = γ E { γ } denotes the average SNR. Then, its MGF is obtained as:

$$\mathbf 2 _ { \mathbf e } \quad \mathbf M _ { \gamma } ( s ) = \frac { 1 } { 1 - \sigma ^ { 2 } \bar { \gamma } s } \mathbf M _ { \xi ^ { 2 } } \left ( \frac { | \omega _ { 0 } | ^ { 2 } \bar { \gamma } s } { 1 - \sigma ^ { 2 } \bar { \gamma } s } \right ) \quad \ \ ( 5 )$$

$$& = [ 1 - \sigma ^ { 2 } \bar { \gamma } s ] ^ { k - 1 } \frac { \supset ^ { \mu ^ { \prime } } _ { \lambda \Omega | \omega _ { 0 } | ^ { 2 } \bar { \gamma } s } _ { \bar { \gamma } s } } { [ 1 - ( \sigma ^ { 2 } + \Omega | \omega _ { 0 } | ^ { 2 } ) \bar { \gamma } s ] ^ { k } }, \quad ( 6 )$$

for ¯ γs &lt; ( σ 2 +Ω | ω 0 | 2 ) -1 .

Proof: We can characterize the SNR based on (7) as:

$$\gamma \left | \xi \sim \bar { \gamma } \left | \mathcal { C N } \left ( \omega _ { 0 } \xi e ^ { j \phi _ { 0 } }, \sigma ^ { 2 } \right ) \right | ^ { 2 }$$

$$\gamma | \xi & \sim \bar { \gamma } \left | \mathcal { C N } \left ( \omega _ { 0 } \xi e ^ { j \phi _ { 0 } }, \sigma ^ { 2 } \right ) \right | ^ { 2 } \\ & \sim \frac { \sigma ^ { 2 } \bar { \gamma } } { 2 } \left [ \mathcal { N } ^ { 2 } \left ( \sqrt { 2 } \omega _ { 0 } \sigma ^ { - 1 } \xi \cos \phi _ { 0 }, 1 \right ) \\ & \quad + \mathcal { N } ^ { 2 } \left ( \sqrt { 2 } \omega _ { 0 } \sigma ^ { - 1 } \xi \sin \phi _ { 0 }, 1 \right ) \right ]$$

$$\stackrel { ( a ) } { \sim } \widetilde { \chi } _ { 1 } ^ { 2 } \left ( \sigma ^ { 2 } \bar { \gamma }, | \omega _ { 0 } | ^ { 2 } \sigma ^ { - 2 } \xi ^ { 2 } \right ),$$

where ( a ) is due to the fact that the sum of two statistically independent squared normal (Gaussian) RVs with unit variances follows the noncentral chi-square distribution. Hence, conditioned on ξ , the SNR follows the noncentral chi-square distribution with 2 d.o.f., noncentrality parameter 2 | ω 0 | 2 σ 2 ξ 2 , and scale factor σ 2 2 . The conditional MGF of γ on ξ is:

$$\mathbf M _ { \gamma | \xi } ( s | \xi ) = \frac { 1 } { 1 - \sigma ^ { 2 } \bar { \gamma } s } \exp \left ( \frac { | \omega _ { 0 } | ^ { 2 } \bar { \gamma } s } { 1 - \sigma ^ { 2 } \bar { \gamma } s } \xi ^ { 2 } \right ), \quad ( 1 0 )$$

for ¯ γσ s &lt; 2 1 . The marginal MGF of the SNR is obtained by taking the expectation of M γ ξ | ( s ξ | ) over ξ , which results in (5). Finally, applying M ξ 2 ( s ) = 1 (1 -Ω ) s k e Ω λs 1 -Ω s yields (6), which completes the proof of Lemma 1.

Lemma 2: Considering γ ∼ F LoS (¯; γ K k,λ, ; Ω) with nonnegative real values K k λ , , , and Ω , where ¯ = γ E { γ } denotes the average SNR. Then, the n th moment of γ is

$$\mathbf E \left \{ \gamma ^ { n } \right \} = n! ( \sigma ^ { 2 } \bar { \gamma } ) ^ { n } \sum _ { i = 0 } ^ { n } \binom { n } { i } \left ( \frac { \Omega | \omega _ { 0 } | ^ { 2 } } { \sigma ^ { 2 } } \right ) ^ { i } L _ { i } ^ { k - 1 } ( - \lambda ), \ ( 1 1 )$$

where L α n ( x ) denotes the generalized Laguerre polynomial.

Theorem 1: Considering γ ∼ F LoS (¯; γ K k,λ, ; Ω) with nonnegative real values K k λ , , , and Ω , where ¯ = γ E { γ } denotes the average SNR. Then, the PDF of γ is:

$$f _ { \gamma } ( x ) & = \frac { 1 } { \sigma ^ { 2 } \bar { \gamma } [ 1 + \Omega | \omega _ { 0 } | ^ { 2 } \sigma ^ { - 2 } ] ^ { k } } \exp \left ( - \frac { x } { \sigma ^ { 2 } \bar { \gamma } } - \lambda \right ) \quad \text{Henc} \\ & \quad \times \widetilde { \Psi } _ { 2 } \left ( k ; \frac { \lambda } { 1 + \frac { \Omega | \omega _ { 0 } | ^ { 2 } } { \sigma ^ { 2 } } }, \frac { \Omega | \omega _ { 0 } | ^ { 2 } } { \sigma ^ { 2 } + \Omega | \omega _ { 0 } | ^ { 2 } } \frac { x } { \sigma ^ { 2 } \bar { \gamma } } \right ), \quad ( 1 2 ) \quad \text{L} \Big \{$$

for x &gt; 0 , where ˜ Ψ ( ; 2 a w, z ) ≜ Ψ ( ; 1 2 a , a ; z, w ) with Ψ ( ) 2 · being the confluent hypergeometric function of two variables [14, Eq. (1.4.9)].

Proof: The PDF of γ is obtained by taking the inverse Laplace transform of M γ ( -s ) from to x -domain as:

$$f _ { \gamma } ( x ) & = \mathbf L ^ { - 1 } \left \{ \mathbf M _ { \gamma } ( - s ) ; s, x \right \} & \text{Upon} \\ & = \exp \left ( - \frac { x \bar { \gamma } ^ { - 1 } } { \sigma ^ { 2 } + \Omega | \omega _ { 0 } | ^ { 2 } } \right ) \\ & \quad \times \mathbf L ^ { - 1 } \left \{ \mathbf M _ { \gamma } \left ( - s + \frac { \bar { \gamma } ^ { - 1 } } { \sigma ^ { 2 } + \Omega | \omega _ { 0 } | ^ { 2 } } \right ) ; s, x \right \}, \, ( 1 4 ) \\ & =$$

where the last equality is due to the frequency shifting property of the Laplace transform. After that, utilizing [15, Eq. (7)] and [15, Eq. (16)], we get:

$$\mathbf L ^ { - 1 } \left \{ \frac { ( a + s ) ^ { k - 1 } } { s ^ { k } } e ^ { \frac { a b } { s } } ; s, x \right \} & = \Phi _ { 3 } ( 1 - k ; 1 ; - a x, a b x ) \ ( 1 5 ) \quad \ \ - \\ & = e ^ { - a x - b } \widetilde { \Psi } _ { 2 } ( k ; b, a x ), \quad ( 1 6 ) \quad \ \ -$$

and after some mathematical manipulations, we obtain (12). This completes the proof of Theorem 1.

Theorem 2: Considering γ ∼ F LoS (¯; γ K k,λ, ; Ω) with nonnegative real values K k λ , , , and Ω , where ¯ = γ E { γ } denotes the average SNR. Then, the CDF of γ is:

$$F _ { \gamma } ( x ) & = \frac { x } { \sigma ^ { 2 } \bar { \gamma } } A ^ { k } \exp \left ( - \frac { A } { \sigma ^ { 2 } \bar { \gamma } } x - B \lambda \right ) & \text{$\underset { \text{introduc} } { \text{function} } \\ & \quad \times \Phi _ { 3 } ^ { ( 3 ) } \left ( 1, 1 - k ; 2 ; \frac { A } { \sigma ^ { 2 } \bar { \gamma } } x, - \frac { B } { \sigma ^ { 2 } \bar { \gamma } } x, \frac { A B \lambda } { \sigma ^ { 2 } \bar { \gamma } } x \right ), \ ( 1 7 ) & \quad \Omega = \frac { 1 } { k + } \\ & \quad \text{$B.\ Inteq$}$$

for x &gt; 0 , where A ≜ σ 2 σ 2 +Ω | ω 0 | 2 , B ≜ Ω | ω 0 | 2 σ 2 +Ω | ω 0 | 2 , and Φ ( n ) 3 ( ) · is the n -variate confluent hypergeometric function [14, Eq. (1.4.17)].

Proof: The CDF of γ is obtained by taking the inverse Laplace transform of s -1 M γ ( -s ) as:

$$\Lambda _ { \ } & \text{Laplace transform or } s ^ { \ } * \text{IM} _ { \gamma ( - s ) } \text{ as} \colon \\ \intertext { n o n - e n o t e s } F _ { \gamma } ( x ) & = \text{L} ^ { - 1 } \left \{ s ^ { - 1 } \text{M} _ { \gamma ( - s ) } ; s, x \right \} & ( 1 8 ) \\ & = \exp \left ( - \frac { x \bar { \gamma } ^ { - 1 } } { \sigma ^ { 2 } + \Omega | \omega _ { 0 } | ^ { 2 } } \right ) \\ ( 1 1 ) & \quad \times \text{L} ^ { - 1 } \left \{ \frac { 1 } { s } \frac { \text{M} _ { \gamma } \left ( - s + \frac { \bar { \gamma } ^ { - 1 } } { \sigma ^ { 2 } + \Omega | \omega _ { 0 } | ^ { 2 } } \right ) } { 1 - \frac { \bar { \gamma } ^ { - 1 } } { \sigma ^ { 2 } + \Omega | \omega _ { 0 } | ^ { 2 } } \frac { 1 } { s } } ; s, x \right \}. \ ( 1 9 ) \\ \text{nial.} \quad \text{sred } \Lambda _ { \ } & \text{sred } \Lambda _ { \ } & \text{sred } \Sigma ( n ) / \vee \vdots \text{sred} \Lambda _ { \ } & \text{sred } \Lambda _ { \ } & \text{sred } \Sigma ( n ).$$

Note that Φ ( n ) 3 ( ) · is expressed as [14, Eq. (1.4.17)]:

Φ

(

b

(

n

)

3

1

, . . . , b

;

c

;

d t, . . . , d

t

)

=

lim

n

-

Φ

1

(

n

)

-

1

b

2

(

b

1

n

1

, . . . , b

n

;

c

;

d t, . . . , d

1

n

-

1

t, d n n

b

t

)

.

n

→∞

Hence, using the Laplace transform of Φ ( n ) 2 ( ) · in [14, Eq. (9.4.55)], the Laplace transform of Φ ( n ) 3 ( ) · is obtained as:

$$( 9. 4. 5 5 ) & \text{], the Laplace transform of } \Phi _ { 3 } ^ { ( n ) } ( \cdot ) \text{ is obtained as} \colon \\ \mathfrak { L } & \left \{ \frac { t ^ { c - 1 } } { \Gamma ( c ) } \Phi _ { 3 } ^ { ( n ) } ( b _ { 1 }, \dots, b _ { n - 1 } ; c ; d _ { 1 } t, \dots, d _ { n } t ) \right \} \\ \mathfrak { L } & = \lim _ { b _ { n } \to \infty } \frac { 1 } { s ^ { c } } \left ( 1 - \frac { d _ { n } } { s } \frac { 1 } { b _ { n } } \right ) ^ { - b _ { n } } \prod _ { i = 1 } ^ { n - 1 } \left ( 1 - \frac { d _ { i } } { s } \right ) ^ { - b _ { i } } \\ \mathfrak { r } & = \frac { 1 } { s ^ { c } } e ^ { \frac { d _ { n } } { s } } \prod _ { i = 1 } ^ { n - 1 } \left ( 1 - \frac { d _ { i } } { s } \right ) ^ { - b _ { i } }. \\ \text{Upon applying the above identity, we obtain (17). This}$$

Upon applying the above identity, we obtain (17). This completes the proof of Theorem 2.

TABLE I RELATION BETWEEN THE FLOS FADING MODEL AND OTHER FADING MODELS.

| Channels            | fLoS Fading Parameters   | fLoS Fading Parameters   | fLoS Fading Parameters   |
|---------------------|--------------------------|--------------------------|--------------------------|
| One-sided normal    | K →∞ ,                   | k = 1 / 2 ,              | λ = 0                    |
| Nakagami- m         | K →∞ ,                   | k = m ,                  | λ = 0                    |
| Rayleigh            | K = 0 ,                  | ∀ k ,                    | ∀ λ                      |
| Rice                | K = K ,                  | k →∞ ,                   | λ = 0                    |
| Rician shadowed     | K = K ,                  | k = m ,                  | λ = 0                    |
| Hoyt (Nakagami- q ) | K = 1 - q 2 2 q 2 ,      | k = 1 / 2 .              | λ = 0                    |

The expressions in (12), (17), and (11) are novel, and their accuracy is validated through comprehensive Monte Carlo simulations. However, the newly proposed fading model's expressions are quite challenging from a computational perspective as compared to the κ µ -fading case. In section V, we introduce an efficient method to numerically compute these functions. The relationship of the fLoS fading distribution with other well-known distribution is summarized in Table I, where Ω = 1 k + λ .

## B. Integer value of d.o.f. parameter, k

When the parameter k in the fLoS model is a positive integer, the PDF expression is simplified. Specifically, using the property of the Φ ( ) 3 · function in [15, Eq. (8)], we can derive a more tractable expression for the PDF, as detailed in the subsequent corollary.

Corollary 1: Considering γ ∼ F LoS (¯; γ K k,λ, ; Ω) with nonnegative real values K λ , , Ω , and k ∈ Z + is a positive integer, where ¯ = γ E { γ } denotes the average SNR. Then, the PDF of γ is:

$$f _ { \gamma } ( x ) & = \frac { A ^ { k } } { \bar { \gamma } \sigma ^ { 2 } } \exp \left ( - A \frac { x } { \sigma ^ { 2 } \bar { \gamma } } - B \lambda \right ) \sum _ { j = 0 } ^ { k - 1 } \frac { ( - 1 ) ^ { j } } { j! } ( 1 - k ) _ { j } \quad \begin{array} { c } \text{The t} \\ \text{Algorit} \end{array} \\ & \quad \times \left ( \frac { B } { A \lambda } \frac { x } { \sigma ^ { 2 } \bar { \gamma } } \right ) ^ { \frac { j } { 2 } } I _ { j } \left ( 2 \sqrt { A B \lambda \frac { x } { \sigma ^ { 2 } \bar { \gamma } } } \right ), \ x > 0. \ ( 2 1 ) \quad \overline { \text{Algorit} } \quad \begin{array} { c } \overline { \text{Algorit} } \end{array} \end{array}$$

We observe that the PDF is expressed as a finite mixture distribution, greatly simplifying the computational complexity of the fLoS fading distribution. Examining Corollary 1 and utilizing the noncentral chi-square distribution definition in (2) and its relation to the κ µ -fading model, we observe the fLoS fading distribution as a finite mixture of κ µ -fading power distributions as:

$$f _ { \gamma } ( x ) = \sum _ { j = 0 } ^ { k - 1 } C _ { j } f _ { \kappa \mu } ( \bar { \gamma } ; \underline { \kappa }, \underline { \mu } ; x ), \ x > 0. \quad ( 2 2 ) \ \frac { \text{Our} } { \text{$\mathbb{ }$} }$$

$$F _ { \gamma } ( x ) = \sum _ { j = 0 } ^ { k - 1 } C _ { j } F _ { \kappa \mu } ( \bar { \underline { \gamma } } ; \underline { \kappa }, \underline { \mu } ; x ), \ x > 0. \quad \ ( 2 3 ) \quad \text{$\text{$\text{$froi}$} \\ -$$

where ¯ := ( γ Bλ + j +1) Aσ γ 2 ¯ , µ := j +1 , κ := Bλ j +1 , C j ≜ (1 -k ) j j ! ( -B/A ) j A k with ∑ k -1 j =0 C j = 1 . Here, f κµ (¯; γ κ, µ ; x ) and F κµ (¯; γ κ, µ ; x ) are the PDF and CDF of the κ µ -fading power as introduced in [16, Eq. (2)] and [16, Eq. (3)], respectively. Due to the finite mixture of κ µ -distributions, performance metrics for the fLoS fading channel (with integer k ) can be readily derived from existing results for the κ µ -fading case, without resorting to further mathematical manipulations.

## IV. PERFORMANCE ANALYSIS AND APPLICATION

## A. Outage Probability

The outage probability (OP), P out , indicates the probability that the instantaneous SNR drops below a set threshold, γ th , given by P out = Pr { γ &lt; γ th } [17]. Consequently, the OP under fLoS fading is directly obtained by computing the CDF of the SNR, provided in (17), as P out = F γ ( γ th ) . In this context, when ¯ γ → ∞ (high transmission SNR), we neglect Φ (3) 3 ( ) · in (17) to obtain the asymptotic OP as:

$$P _ { \text{out} } | _ { \uparrow } = A ^ { k } e ^ { - B \lambda } \gamma _ { \text{th} } ( \sigma ^ { 2 } \bar { \gamma } ) ^ { - 1 }. \quad \quad ( 2 4 ) \quad \text{chara}$$

A close observation of (24) indicates that the fLoS fading model's diversity order, which is the rate at which the OP decreases in the high-SNR regime, is d = 1 .

## B. Unified Ergodic Capacity

The Ergodic capacity (EC) per unit bandwidth, in bps/Hz, is given by C = 1 ln 2 ∫ ∞ 0 ln(1 + γ f ) γ ( γ )d γ [18]. However, the equation is difficult to derive due to the SNR PDF, f γ ( γ ) , involving confluent hypergeometric functions. To address this limitation, we propose a unified based method to approximate C [bps/Hz]. Our method tailors the Prony's algorithm presented in [19, Algorithm 2.1] to approximate f ( x ) = ln(1 + x ) as:

$$\overset { \text{per} } { \text{or} }, \quad f ( x ) = \sum _ { s = 1 } ^ { S } \sum _ { k = 1 } ^ { M _ { s } } c _ { s, k } e ^ { - T _ { s, k } x } \overset { \triangle } { = } \hat { f } ( x ), \ x \in [ l _ { s }, u _ { s } ]. \quad ( 2 5 )$$

The tailored Prony's approximation method is detailed in Algorithm 1.

## Algorithm 1: Tailored Prony's approximation method

Input: f , S , M s , L s ≫ 2 M s , h s , l s , and u s , s = 1 , . . . , S , u s +1 = l s ;

1 for s = 1 to S do

- 2 Compute [ f ] i = f ( x i ) ≜ f i , i = 1 , . . . , L s +1 and x i = l s + h s ( i -1) with u s = l s + h L s s ;
- 3 Perform Prony's method to approximate f ( x ) for x ∈ [ u , l s s ] and obtain the Prony's coefficients T k and c k , k = 1 , . . . , M s ;

4

end

Output: Prony's coefficients;

From (25) and the fact that E { e sγ } = M γ ( s ) , the EC per unit bandwidth is approximated as:

$$j \stackrel { \vartriangle } { = } \sum _ { \substack { s = 1 } { \cdot \pi \rangle } } ^ { S } \sum _ { k = 1 } ^ { M } c _ { s, k } \left [ F ( T _ { s, k }, u _ { s } ) - F ( T _ { s, k }, l _ { s } ) \right ], \quad ( 2 6 )$$

where F T,x ( ) = ∫ x 0 e -Ts f γ ( s )d s , or [20, Eq. (4)]

$$F ( T, x ) = \mathbf L ^ { - 1 } \left \{ \frac { \mathbf M _ { \gamma } ( - T - p ) } { p } ; p, x \right \}, \quad ( 2 7 )$$

which is derived in a similar manner as F γ ( x ) in (17).

Observing that ln(1 + x ) → ln( x ) as x →∞ , we directly derive the asymptotic EC per unit bandwidth as:

$$\overline { C } | _ { \uparrow } = \sum _ { s = 1 } ^ { S } \sum _ { k = 1 } ^ { M } \frac { c _ { s, k } e ^ { T _ { s, k } } } { \ln 2 } \left [ F ( T _ { s, k }, u _ { s } ) - F ( T _ { s, k }, l _ { s } ) \right ]. \ \ ( 2 8 )$$

C. Multiple-input single-output (MISO) systems under channel aging

The SNR distribution over fLoS fading channel finds application in many real-world scenarios. Within these scenarios, the intended physical interpretation of the fLoS fading distribution may vary, yet the distribution itself remains useful in characterizing SNR distributions effectively.

In downlink MISO systems, with a high-speed moving receiver Rx, channel aging poses a significant challenge, causing discrepancies between current and estimated channel state information (CSI) [4]. This mismatch degrades the wireless system performance and reliability. Specifically, the true CSI h conditioned on the estimated CSI ˆ h is given by h h | ˆ = ρ p h + ¯ ρ z , where z ∼ CN ( 0 N , I N ) , with N being the number of transmit antennas, ρ being the temporal correlation coefficient, and ¯ = ρ √ 1 -ρ 2 [4]. The effective fading from transmitter to receiver, h eff , under Maximal Ratio Transmission (MRT) beamforming, is h eff | ˆ h ∼ CN ( ρ ∥ ˆ h ∥ , ρ ¯ ) 2 , where ∥ ˆ h ∥

is the L 2 -norm of the estimated CSI. Assuming Rician fading: ˆ h ∼ CN ( √ κ √ κ +1 ˇ h , I N κ +1 ) ∈ C N , where κ is the Rician factor and ˇ h represent the LoS components. Hence, the SNR at Rx is γ ≜ P N W 0 L eff | h eff | 2 , where P [W] is the Rx's transmission power, N 0 [W/Hz] is the noise spectral density, W [H] is the channel bandwidth, and L eff denotes the effective path loss between Tx and Rx. Based on (7) and noting that ∥ ˆ h ∥ 2 ∼ ˜ χ 2 N ( 1 κ +1 , Nκ ) , the true distribution of the received SNR, γ , is given by:

$$\gamma \sim \mathcal { F } _ { \text{LoS} } \left ( \bar { \gamma } ; \rho ^ { 2 } / \bar { \rho } ^ { 2 }, N, \kappa N, [ \kappa + 1 ] ^ { - 1 } \right ), \quad ( 2 9 )$$

where ¯ γ ≜ PL eff N W 0 . Thus, when Rx moves with high speed and causes channel aging, the effective channel's true diversity order between Tx and Rx remains 1 regardless the number of transmit antennas, as discussed in the previous section.

## V. NUMERICAL RESULTS

Although the confluent hypergeometric functions Φ ( n ) 3 ( ) · , Φ ( ) = Φ 3 · (2) 3 ( ) · , and Ψ ( ) 2 · are not available in MATLAB nor in Mathematica, we can numerically compute them as:

$$& \Phi _ { 3 } ^ { ( n ) } ( b _ { 1 }, \dots, b _ { n - 1 } ; c ; d _ { 1 } t, \dots, d _ { n } t ) & \cdot \cdot \cdot \\ & \quad = \frac { \Gamma ( c ) } { t ^ { c - 1 } } \int _ { L - j T } ^ { L + j T } s ^ { - c } \frac { e ^ { s t + \frac { d _ { n } } { s } } } { j 2 \pi } \prod _ { i = 1 } ^ { n - 1 } \left ( 1 - \frac { d _ { i } } { s } \right ) ^ { - b _ { i } } d s, \ ( 3 0 ) & \quad \text{ability 1} \\ & \mathfrak { m } \, l \, \omega \cdot d \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \alpha \cdot \beta \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \,.$$

$$\Psi _ { 2 } ( a ; d ; a ; y, x t ) = e ^ { y + x t } \Phi _ { 3 } ( d - a ; d, a ; - x t, x y t ), \quad ( 3 1 ) \quad \text{on} \ \text{inc}$$

where L = max( d , . . . , d 1 n -1 , 0) + ϵ with ϵ &gt; 0 and T ≫ 0 are manually chosen to ensure the computation accuracy. For example, we set ϵ = 1 2 and T = 10 3 in all the simulations.

## A. Effect of fading parameters

Fig. 1. PDF of the SNR over the fLoS fading with different values for the d.o.f., k , where ¯ = 0 γ dB, λ = 5 and K = 10 dB.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00575v1-figure-85ac6bd7-d183-4a11-b740-e2c6785ebebf.jpg)

We now demonstrate how the fading parameters of the fLoS fading model influence the distribution shape. Due to space constraints, we concentrate on the SNR PDF and use Ω = 1 k + λ to ensure E { ξ 2 } = 1 and E { γ } = ¯ γ .

Fig. 1 illustrates the changes in the PDF as the d.o.f. parameter varies. Increasing k reduces the probability for having very low SNR values by concentrating the SNR around

Fig. 2. SNR PDF over fLoS fading versus λ with k = 1 and K = 10 dB. Solid lines correspond to the ¯ = 1 γ case. Dashed lines correspond to the ¯ = 2 γ case.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00575v1-figure-f119b322-4c23-4ab8-9352-7aea61fed69d.jpg)

its mean value ¯ γ . Additionally, higher values of k alleviate the fading severity, with k = 0 indicating the most severe fading conditions.

Fig. 2 depicts the changes in the PDF as the noncentrality parameter λ varies. Similarly, increasing λ reduces the probability for very low SNR values by shifting the distribution to the right causing the peak to move away from zero, while the PDF is not concentrated around its mean as compared to increasing the d.o.f.. Smaller values of λ also result in distributions with heavier tails, reducing the probability for having low SNRs.

## B. Performance Analysis

Fig. 3. OP as a function of the normalized ¯ γ for fLoS fading with k = 1 5 . , K = 13 dB and different values of λ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00575v1-figure-f6be1868-7499-4256-a1f3-ebd4e55e9719.jpg)

Fig. 3 presents the OP as a function of the normalized average SNR, γ/γ th , with different instances of the noncentrality parameter, λ . Analytical and asymptotic results for the OP are obtained by applying (17) and (24), respectively. Simulation results are presented validate the expressions under the current settings. We observe that the OP decreases as λ is increased, especially with increased K , and that the asymptotic OP results are highly accurate for sufficiently large ¯ γ/γ th .

Fig. 4. EC per unit bandwidth, C , as a function of ¯ γ for fLoS fading with K = 13 dB and λ = 0

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00575v1-figure-1b5f634e-7f0b-49a5-9d24-0aa3ca4b6759.jpg)

Fig. 4 depicts the EC considering various values for k . Simulation results (markers) verify the accuracy of (26) and (28). Moreover, the capacity of the AWGN case is included as a benchmark, as it always upper bounds the capacity of a fading channel due to Jensen's inequality. As observed, the performance gap with respect to the AWGN case is reduced as k is increased.

## C. Application

Fig. 5. Coverage probability as a function of the time index, n , for fLoS fading with N = 4 , κ = 10 dB, ¯ γ/γ th = 10 dB, and different values of the normalized Doppler shift, f d T s .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00575v1-figure-49766573-f7c5-4b3c-b407-afebb3cbb34d.jpg)

Fig. 5 depicts the coverage probability, P cov = 1 -P out , versus the time instant index n for varying normalized Doppler shifts, where the data transmission starts at n = 0 . The downward trend illustrates a decreasing coverage probability with increasing n . Additionally, as f T d s increases, the first zero position shifts leftward, indicating a reduction in the coherence time. The lower limit, e -¯ γ/γ th , signifies an effective fading channel behavior similar to Rayleigh fading, one of the most unfavorable fading scenarios.

## VI. CONCLUSION

We introduced the new fLoS fading model, motivated primarily by its ability to encompass the lognormal shadowing and shadowed Rician fading models. This fLoS fading distribution simplifies into a finite combination of the more straightforward κ µ -distribution when k is a positive integer.

The performance analysis of wireless communication systems under the influence of the proposed fLoS fading model was examined using the OP and EC, where a unified Prony's method is proposed to approximate the EC. Finally, we presented an application for the fLoS fading distribution in characterizing wireless systems under the influence of channel aging.

## REFERENCES

- [1] W. Jiang, B. Han, M. A. Habibi, and H. D. Schotten, 'The road towards 6G: A comprehensive survey,' IEEE Open J. Commun. Soc. , vol. 2, pp. 334-366, Feb. 2021.
- [2] C.-X. Wang, X. You, Gao et al. , 'On the road to 6G: Visions, requirements, key technologies, and testbeds,' IEEE Commun. Surveys Tuts. , vol. 25, no. 2, pp. 905-974, Feb. 2023.
- [3] E. Bj¨rnson, O. Ozdogan, and E. G. Larsson, 'Reconfigurable intelligent o ¨ ¨ surfaces: Three myths and two critical questions,' IEEE Commun. Mag. , vol. 58, no. 12, pp. 90-96, Dec. 2020.
- [4] R. Chopra, C. R. Murthy, H. A. Suraweera, and E. G. Larsson, 'Performance analysis of FDD massive MIMO systems under channel aging,' IEEE Trans. Wirel. Commun. , vol. 17, no. 2, pp. 1094-1108, Feb. 2018.
- [5] G. Geraci, D. L´ opez-P´ erez, M. Benzaghta, and S. Chatzinotas, 'Integrating terrestrial and non-terrestrial networks: 3D opportunities and challenges,' IEEE Commun. Mag. , vol. 61, no. 4, pp. 42-48, April 2023.
- [6] M. M. Azari et al. , 'Evolution of non-terrestrial networks from 5G to 6G: A survey,' IEEE Commun. Surv. Tut. , vol. 24, no. 4, pp. 2633-2672, Aug. 2022.
- [7] C.-X. Wang, J. Huang, H. Wang, X. Gao, X. You, and Y. Hao, '6G wireless channel measurements and models: Trends and challenges,' IEEE Veh. Technol. Mag. , vol. 15, no. 4, pp. 22-32, Dec. 2020.
- [8] J. D. V. Gutierrez and C. Oestges, 'Measurement-based large-scale propagation characterization in 5G micro-cells at 3.8 GHz,' IEEE Trans. on Veh. Technol. , vol. 71, no. 11, pp. 11 325-11 335, Nov. 2022.
- [9] S. Dang, J. Ye, M. A. Beach, and M. Chafii, 'On optimality of Gamma approximation for Lognormal shadowing models,' IEEE Antennas Wirel. Propag. Lett. , vol. 22, no. 5, pp. 1099-1103, May 2023.
- [10] H. Wang, X. Yin, J. Rodr´ ıguez-Pi˜eiro, J. Lee, and M.-D. Kim, 'Shadn owing and multipath-fading statistics at 2.4 GHz and 39 GHz in vehicleto-vehicle scenarios,' in Proc. 14th Eur. Conf. Antennas Propag. , Mar. 2020, pp. 1-5.
- [11] J. L´ opez-Fern´ andez, P. R. Espinosa et al. , 'A fluctuating Line-of-Sight fading model with double-Rayleigh diffuse scattering,' IEEE Trans. Veh. Technol. , vol. 71, no. 1, pp. 1000-1003, Jan. 2022.
- [12] J. Salo, H. M. El-Sallabi, and P. Vainikainen, 'Statistical analysis of the multiple scattering radio channel,' IEEE Trans. Antennas Propag. , vol. 54, no. 11, pp. 3114-3124, Nov. 2006.
- [13] A. Abdi, W. Lau, M.-S. Alouini, and M. Kaveh, 'A new simple model for land mobile satellite channels: first- and second-order statistics,' IEEE Trans. Wirel. Commun. , vol. 2, no. 3, pp. 519-528, May 2003.
- [14] H. M. Srivastava and P. W. Karlsson, 'Multiple gaussian hypergeometric series,' USA:Wiley , 1985.
- [15] N. Y. Ermolova and O. Tirkkonen, 'Laplace transform of product of generalized Marcum Q, Bessel I, and power functions with applications,' IEEE Trans. Sig. Proc. , vol. 62, no. 11, pp. 2938-2944, June 2014.
- [16] M. D. Yacoub, 'The κ µ -distribution and the η µ -distribution,' IEEE Antennas and Propag. Mag. , vol. 49, no. 1, pp. 68-81, Feb. 2007.
- [17] P. Bhardwaj, E. S. Karnawat, and S. Zafaruddin, 'An exact statistical representation of α η κ µ ---fading model for THz wireless communication,' in 2023 IEEE GC Wkshps . IEEE, 2023, pp. 661-666.
- [18] F. Yilmaz and M.-S. Alouini, 'A novel ergodic capacity analysis of diversity combining and multihop transmission systems over generalized composite fading channels,' in 2012 IEEE ICC , June 2012, pp. 46054610.
- [19] R. Zhang and G. Plonka, 'Optimal approximation with exponential sums by a maximum likelihood modification of Prony's method,' Advances in Computational Mathematics , vol. 45, no. 3, pp. 1657-1687, 2019.
- [20] F. J. Lopez-Martinez, J. M. Romero-Jerez, and J. F. Paris, 'On the calculation of the incomplete MGF with applications to wireless communications,' IEEE Trans. Commun. , vol. 65, no. 1, pp. 458-469, Jan. 2017.
Thanh Luan Nguyen Georges Kaddoum
Fluctuating Line-of-Sight Fading Distribution: Statistical Characterization and Applications
We introduce the fluctuating Line-of-Sight (fLoS) fading model, characterized
by parameters $K$, $k$, $\lambda$, and $\Omega$. The fLoS fading distribution
is expressed in terms of the multivariate confluent hypergeometric functions
$\Psi_2$, $\Phi_3^{(n)}$, and $\Phi_3 = \Phi_3^{(2)}$ and encompasses
well-known distributions, such as the Nakagami-$m$, Hoyt, Rice, and Rician
shadowed fading distributions as special cases. An efficient method to
numerically compute the fLoS fading distribution is also addressed. Notably,
for a positive integer $k$, the fLoS fading distribution simplifies to a finite
mixture of $\kappa$-$\mu$ distributions. Additionally, we analyze the outage
probability and Ergodic capacity, presenting a tailored Prony's approximation
method for the latter. Numerical results are presented to show the impact of
the fading parameters and verify the accuracy of the proposed approximation.
Moreover, we illustrate an application of the proposed fLoS fading distribution
for characterizing wireless systems affected by channel aging.