# Modeling of Non-WSSUS Double-Rayleigh Fading Channels for Vehicular Communications.

1. IntroductionTerrestrial vehicle-to-vehicle (V2V) communication systems are emerging as an enabling technology for a variety of new wireless applications and services, such as information relaying for mobile cellular networks [1] and peer-to-peer data transmission for vehicular communications [2]. Some of the most important applications of these systems target the prevention of vehicular accidents and the optimization of traffic flow. Such applications have caught the attention of the automotive industry and different government bodies around the world, who have become major promoters of the V2V communications technology [3].

One of the main challenges in the design of V2V communication systems is to develop a robust air interface that supports delay sensitive applications under the constraints of a rapidly changing propagation environment and a dynamic network topology. To successfully design and optimize such an air interface, a realistic reference model of the time-frequency (TF) dispersive V2V fading channel is required. This is of primary importance, since the performance of the wireless communication systems is highly influenced by the propagation environment. In addition to measurement-based models, proper analytic channel models are needed that provide insights into the physics of V2V radio reception and, at the same time, that lend themselves to mathematical and numerical system performance investigations.

Important advances in the analytical characterization of fixed-to-mobile (F2M) multipath fading channels were prompted by the emergence of the mobile cellular communication systems in the late 1970s. However, the modeling of fading channels for V2V communications not only demands more exhaustive research work, but also requires a shift of paradigm, because some of the assumptions that are often invoked for the characterization of F2M channels are not valid when the terminals at both ends of the radio link are able to move at high speeds. For example, most of the existing statistical models for F2M channels have been formulated assuming the fulfillment of the wide-sense stationary uncorrelated scattering (WSSUS) condition introduced by Bello in [4] (e.g., see [5-7]). This assumption facilitates the mathematical analysis of TF-dispersive channels, as it implies that the channel's statistics are simultaneously wide-sense stationary (WSS) in the time and the frequency domains. Nevertheless, recent empirical investigations carried out in vehicular communication environments suggest that the WSSUS condition is not valid for V2V channels [8]. While the nonstationary features of multipath wireless channels have been a subject of analysis since the early days of the mobile radio communications (e.g., see [4, 9]), the modeling of such nonstationarities has predominantly been addressed from a large-scale propagation perspective. Measured data obtained independently in [10-12] demonstrates that the nonstationary characteristics of V2V channels are also meaningful at a small-scale level.

Empirical investigations have further shown that the signal fading induced by V2V channels is more severe than the one produced by F2M channels [13], which is typically modeled by Rayleigh or Rice distributions. The exacerbation of signal fading is not surprising if one considers that, in a V2V communications system, the mobile terminals are located at ground level. This scenario increases the chances of receiving echoes of the transmitted signal that interact with multiple interfering objects (IOs) on their way to the receiver antenna. As a result of such multiple interactions, the received multipath signal is subject to a form of cascaded fading that is modeled by worse-than-Rayleigh distributions [13,14] (e.g., the double-Rayleigh [15] and double-Rice distributions [16]).

Notable contributions to the analytical characterization of non-WSSUS V2V channels have recently been made in [1719] on the basis of the geometry-based statistical modeling approach. The geometrical channel models proposed in these papers assume the propagation of spherical waves to account for the nonstationarities of V2V channels stemming from small-scale propagation. In the spherical wave propagation (SWP) framework, the angle of departure (AOD) and angle of arrival (AOA) of the received multipath signals are determined by the instantaneous spatial position of the transmitting and receiving mobile stations (MSs). The angular statistics of the resulting channel models are therefore time-dependent. This feature is particularly convenient for the characterization of nonstationary small-scale channels but renders the mathematical analysis of the channel's statistics a cumbersome task.

To facilitate the modeling and analysis of non-WSSUS V2V channels, we recently proposed in [20-22] a novel framework that builds instead on the principles of plane wave propagation (PWP). Our proposal is well suited for the analysis of V2V radio reception over small local areas spanning a few tens of wavelengths, where a plane wave approximation of the more realistic spherical electromagnetic waves can be applied. For such propagation scenarios, the angular statistics of the V2V channel can be modeled by time-invariant distributions, which are more mathematically tractable than their time-varying counterparts. The focus of [20,21] was on the characterization of nonstationary Rayleigh fading channels for single-input single-output (SISO) and multiple-input multiple-output (MIMO) V2V communication systems, respectively. In [21], the MSs are assumed to move at constant speeds and on linear trajectories, whereas the effects of acceleration and nonlinear motion are investigated in [20] following a parallel approach to the TF analysis techniques employed in [23]. On the other hand, in [22], we apply our modeling framework to the characterization of non-WSSUS SISO V2V channels that experience double-Rayleigh fading.

In this paper, we complete our preliminary work presented in [22] by providing a detailed description and a thorough statistical analysis of the proposed geometrical model for non-WSSUS V2V double-Rayleigh fading channels. The scope and depth of the work in [22] are expanded here as follows:

(i) Important statistical quantities of the proposed channel model, such as the autocorrelation functions in the time and frequency domains, as well as the TF-dependent delay and Doppler profiles, were not investigated in [22]. An in-depth analysis of these statistical quantities is presented here.

(ii) Details on the derivations of the four-dimensional (4D) TF correlation function (4D TF-CF) of the proposed channel model were not presented in [22]. An outline of the derivations is given here in the Appendix.

(iii) In [22], our discussion of the proposed channel model's stationary (nonstationary) characteristics was constrained to a single paragraph due to space limitations. In this paper, we complement our discussion with important additional remarks.

(iv) Finally, new numerical examples are presented in this paper to illustrate our findings regarding the autocorrelation, spectral, and stationary (nonstationary) characteristics of the proposed geometry-based statistical model (GBSM) for non-WSSUS V2V double-Rayleigh fading channels.

The remainder of the paper is organized as follows. Our proposal for the geometrical modeling of nonstationary V2V double-Rayleigh fading channels is presented in Section 2. In Section 3, we derive general expressions for the envelope and phase distributions, the 4D TF-CF, and the TF-dependent delay and Doppler profiles of the proposed channel model. It is worth pointing out that such expressions are valid regardless of the underlying geometry of the propagation area. In Section 4, we compute a closed-form solution of the 4D TF-CF by considering the particular case of the geometrical two-ring scattering model. Numerical examples illustrating our theoretical findings are presented in Section 4. Finally, our conclusions are given in Section 5.

Notation. The complex conjugate, the argument, and the absolute value operations are denoted by [(*).sup.*], arg{*}, and [absolute value of *], respectively. Vectors are written in bold face. The transpose operation is denoted by [(*).sup.[dagger]], [parallel] * [parallel] stands for the Euclidean norm, and the scalar product between two vectors [z.sub.1] and [z.sub.2] is represented as <[z.sub.1], [z.sub.2]>. The operator E{*} designates the statistical expectation. The set of positive real numbers is denoted by [R.sup.+], and the operator card(*) indicates set cardinality.

2. The Proposed Geometrical Model for Non-WSSUS V2V Channels

2.1. Geometrical Modeling of the Propagation Scenario. The aim of this paper is to model the nonstationarities of TF-dispersive V2V channels stemming from small-scale propagation. We are particularly interested in characterizing the nonstationarities arising from the time-varying nature of the propagation delays. For that purpose, we consider a SISO V2V communication system and assume that the transmitted signal reaches the receiver antenna through a double interaction with nonmoving IOs randomly located in the propagation environment. Specifically, we assume that the transmitted signal interacts first with a set [S.sub.T] of L IOs (card([S.sub.T]) = L) that are located in the surroundings of the transmitting MS ([T.sub.x]). Then, the L scattered signals that result from such an interaction impinge on a second set [S.sub.R] of M IOs (card([S.sub.R]) = M) that lie in the vicinity of the receiving MS ([R.sub.X]). Thereby, a total of L x M double-scattered waves are produced, which combine with one another at the receiver antenna. Figure 1 shows an illustration of the propagation scenario under consideration at the time instant t = [t.sub.0] when the MSs start to communicate with each other.

In Figure 1, the L IOs in the set [S.sub.T] are represented by black dots, whereas white dots stand for the M IOs in [S.sub.R]. We denote the lth IO in [S.sub.k] by [S.sup.k.sub.l], for l [member of] {1, 2, ..., card([S.sub.k])}, and k [member of] {T, R}. The positions of [T.sub.X] and [R.sub.X] at the time when the communications begin (i.e., at t = [t.sub.0]) are denoted by the time-invariant vectors [O.sub.T] and [O.sub.R]. The distance between [O.sub.T] and [O.sub.R] is given by D. The velocity vectors of [T.sub.X] and [R.sub.X] are represented by [v.sub.T] and [v.sub.R], respectively. We assume that [T.sub.X] and [R.sub.X] move at constant speeds and on linear trajectories.

The time-invariant vectors [[??].sup.k.sub.l], for l [member of] {1, 2, ..., card([S.sub.k])} and k [member of] {T, R}, indicate the position of the lth IO [S.sup.k.sub.l] in [S.sub.k], with respect to the fixed reference point [O.sub.k]. In addition, the instantaneous position of the lth IO [S.sup.T.sub.l] in [S.sub.T], as seen from the perspective of the moving [T.sub.X], can be represented by the time-varying vector [p.sup.T.sub.l](t). Analogously, the time-varying vector [p.sup.R.sub.m](t) indicates the instantaneous position of the mth fixed IO [S.sup.R.sub.m] in [S.sub.R], as seen by the moving [R.sub.X]. Regardless of the geometrical configuration of the propagation scenario, we can express these two vectors as

[P.sup.k.sub.l](t) = [[??].sup.k.sub.l] - t * [v.sub.k] (1)

for k [member of] {T, R} and l [member of] {1, 2, ..., card([S.sub.k])}. On the other hand, the position of the mth nonmoving IO [S.sup.R.sub.m] in [S.sub.R] relative to that of the lth IO [S.sup.T.sub.l] in [S.sub.T] is given by the time-invariant vector

[P.sup.S.sub.l,m] = D - [[??].sup.T.sub.l] + [[??].sup.R.sub.m] (2)

for l [member of] {1, 2, ..., card([S.sub.T])} and m [member of] {1, 2, ..., card([S.sub.R])}, where D = [O.sub.R] - [O.sub.T] and [parallel]D[parallel] = D. The time- invariant vectors [u.sup.T.sub.l], [u.sup.S.sub.l,m], and [u.sup.R.sub.m] are unit vectors which point at the direction of propagation of the waves that travel from [T.sub.X] to [S.sup.T.sub.l], from [S.sup.T.sub.l] to [S.sup.R.sub.m], and from [S.sup.R.sub.m] to [R.sub.X], respectively. Note that these vectors do not need to point exactly towards [S.sup.T.sub.l], [S.sup.R.sub.m], and [O.sub.R], if we assume (as we do in this paper) the propagation of plane waves. If, on the other hand, we consider the propagation of spherical waves, then such vectors should point exactly towards the aforementioned observation points. In fact, in the case of SWP, [u.sup.T.sub.l] and [u.sup.R.sub.m] should be modeled as time- varying vectors, since the position of the observer relative to the source changes over time because [T.sub.X] and [R.sub.X] are moving [18].

Geometrical configurations of the propagation scenario similar to the one shown in Figure 1 have already been considered in other papers for modeling SISO and MIMO V2V fading channels (e.g., see [24-26]). However, the relative positions among the MSs and the IOs are characterized in this paper following a vector framework that allows capturing the temporal dynamics of the channel in a simple and compact manner. In contrast, such relative positions are typically modeled in the state of the art by a framework of time-invariant scalar quantities (angles and distances) that provide a less flexible description of the propagation scenario's temporal dynamics. Capitalizing on the flexibility of the vector framework presented here, we propose in the following subsection a general model for the channel impulse response (CIR) of V2V fading channels that is valid regardless of the geometrical configuration of the propagation area, as long as the plane wave model applies and the transmitted signal arrives at [R.sub.X] via a double interaction with IOs.

2.2. Mathematical Model of the Channel Impulse Response. We define the small-scale CIR in the equivalent baseband at time t due to an impulse applied t seconds in the past by the superposition of Lx M plane waves as follows:

[mathematical expression not reproducible]. (3)

In (3), [j.sup.2] = -1; [[theta].sub.0] is the initial phase of the transmitted signal; [g.sup.k.sub.l] and [[theta].sup.k.sub.l] stand for the gain and phase shift, respectively, introduced by the interaction of the transmitted signal with the lth IO [S.sup.k.sub.l] in [S.sub.k], for l [member of] {1, 2, ..., card([S.sub.k])} and k [member of] {T, R}. The Dirac delta function is denoted by [delta](*), and

[mathematical expression not reproducible] (4)

is a rectangular windowing function introduced as a means to limit the length of the CIR h(t; [tau]) within an interval of length [T.sub.0] inside of which the large-scale variations of the channel are negligible, and the plane wave approximation is reasonably justified. For simplicity, and without loss of generality, we will henceforth assume that [t.sub.0] = 0.

The time-varying parameters [[tau].sub.l,m](t) and [[??].sub.l,m](t) stand for the instantaneous propagation delay and path-length dependent phase rotation, respectively, of the multipath signal that arrives at [R.sub.X] via the interaction with the IOs [S.sup.T.sub.l] and [S.sup.R.sub.m]. These two parameters are related to each other by

[[??].sub.l,m] (t) = 2[pi][f.sub.c][[tau].sub.l,m](t), (5)

where [f.sub.c] = [C.sub.l][lambda] is the carrier frequency, [C.sub.l] denotes the speed of light, and [lambda] is the transmitted signal's wavelength. We model the propagation delays as

[[tau].sub.l,m](t) = <[p.sup.T.sub.l](t), [u.sup.T.sub.l]> + <[p.sup.S.sub.l,m], [u.sup.S.sub.l,m]> - <[p.sup.R.sub.m](t), [u.sup.R.sub.m]>. (6)

The sum of scalar products at the right-hand side of (6) provides the path length of the plane wave that travels from [T.sub.X] to [R.sub.X] via [S.sup.T.sub.l] and [S.sup.R.sub.m] [27]. The path length is computed in (6) by considering the instantaneous position of the moving [T.sub.X] and [R.sub.X], which is accounted for by the time-varying position vectors [p.sup.T.sub.l](t) and [p.sup.R.sub.m](t). Equation (6) therefore allows modeling the time-varying nature of the propagation delays. This is a salient feature of our proposal, because to the best of our knowledge, the GBSMs for V2V channels that have been proposed under similar considerations (i.e., by assuming a double interaction with IOs and the propagation of plane waves) do not take into account the temporal dynamics of the propagation delays (see, e.g., [24-26]). Such geometrical channel models define the propagation delays as time-invariant quantities [[tau].sub.l,m] that depend only on the position of [T.sub.X] and [R.sub.X] at time t = 0. (When the propagation of plane waves is assumed, the propagation delays are modeled in the state of the art as [[tau].sub.l,m] = [<[[??].sup.T.sub.l], [u.sup.T.sub.l]> + <[p.sup.T.sub.l,m], [u.sup.S.sub.l,m]> - <[[??].sup.R.sub.l], [u.sup.R.sub.m]>]/[C.sub.l], [for all]m.)

From (1), and assuming that the unit vectors [u.sup.T.sub.l], [u.sup.S.sub.m], and [u.sup.R.sub.m] are collinear with the time-invariant vectors [[??].sup.T.sub.l], [p.sup.S>sub.l,m], and [[??].sup.R.sub.m], respectively, in such a way that

<[[??].sup.T.sub.l], [u.sup.T.sub.l]> = [parallel][[??].sup.T.sub.l][parallel] (7)

<[[??].sup.S.sub.l,m], [u.sup.S.sub.l,m]> = [parallel][P.sup.S.sub.l,m][parallel] (8)

<[[??].sup.R.sub.m], [u.sup.R.sub.m]> = [parallel][[??].sup.R.sub.m][parallel] (9)

we have

<[p.sup.T.sub.l](t), [u.sup.T.sub.l]> = [parallel][[??].sup.T.sub.l][parallel] - t <[v.sub.T], [u.sup.T.sub.l]> (10)

<[p.sup.S.sub.l,m], [u.sup.S.sub.l,m]> = [parallel][P.sup.S.sub.l,m][parallel] - [parallel]D - [[??].sup.T.sub.l] + [[??].sup.R.sub.m][parallel] (11)

<[p.sup.R.sub.m](t), [u.sup.R.sub.m]> = -[parallel][[??].sup.R.sub.m][parallel] - t<[v.sub.R], [u.sup.R.sub.m]>. (12)

Thereby, the propagation delays [[tau].sub.l,m] in (6) can be written as

[[tau].sub.l,m] = [parallel][[??].sup.T.sub.l][parallel] + [parallel][P.sup.S.sub.l,m][parallel] + [parallel][[??].sup.R.sub.m][parallel]/ [C.sub.l] + t [f.sup.D.sub.l,m]/[f.sub.c], (13)

where

[f.sup.D.sub.l,m] = [f.sup.T.sub.l] + [f.sup.R.sub.m] (14)

is a Doppler frequency shift due to the combined movement of [T.sub.X] and [R.sub.X], and

[f.sup.T.sub.l] = <[v.sup.T], [u.sup.T.sub.l]>/[lambda] (15)

[f.sup.R.sub.m] = <[v.sub.R], [u.sup.R.sub.m]>/[lambda]. (16)

It is worth highlighting that the channel modeling framework defined by (3)-(16) is not restricted to a particular geometrical arrangement of the IOs' positions. Moreover, the proposed framework is valid for both two-dimensional (2D) and three-dimensional (3D) propagation scenarios, as it applies regardless of whether the vectors introduced in Figure 1 are defined in 2D or 3D spaces.

3. Statistical Properties of the Proposed Non-WSSUS V2V Channel Model

3.1. Considerations. In this section, we analyze the first-order (FO) statistics, the correlation properties, and the spectral characteristics of the proposed GBSM for V2V channels. To reduce the number of variables in our analysis, we will restrict our attention to a 2D propagation scenario. We point out, however, that the extension with respect to 3D propagation is straightforward. In addition, we will make the following considerations:

(i) The gains [g.sup.k.sub.l] in (3) are statistically independent, although not necessarily identically distributed, positive random variables (r.v.), each having a probability density function (PDF) [mathematical expression not reproducible] (c), for l [member of] {1, 2, ..., card([S.sub.k])}, k [member of] {T, R}, and c [greater than or equal to] 0.

(ii) The phases [[theta].sup.k.sub.l], l [member of] {1, 2, ..., card([S.sub.k])}, k [member of] {T,R}, are statistically independent r.v. uniformly distributed over [-[pi], [pi]).

(iii) The direction of the time-invariant position vector [[??].sup.T.sub.l] is modeled by a random angle [[phi].sup.T.sub.l] characterized by a circular PDF [p.sup.T.sub.[phi]]([[phi].sub.T]), for all l, and [[phi].sub.T] [member of] [-[pi], [pi]). Due to the collinearity between [[??].sup.T.sub.l] and [u.sub.l] (see (7)), [[phi].sup.T.sub.l] can be identified as the AOD of the plane wave that interacts with the lth IO [S.sup.T.sub.l] in [S.sub.T].

(iv) The r.v. [[phi].sup.R.sub.m] describes the direction of the time-invariant position vector [[??].sup.R.sub.m] and is characterized by a circular PDF [p.sup.R.sub.[phi]]([[phi].sub.R]), for all l, and [[phi].sub.R] [member of] [~[pi], [pi]). Due to the collinearity between [[??].sup.R.sub.m] and [u.sub.m] (see (9)), we can refer to [[phi].sup.R.sub.m] as the AOA of the mth multipath signal that arrives at [R.sub.X] through the interaction with the mth IO [S.sup.R.sub.m] in [S.sub.R].

(v) The time-invariant vectors [[??].sup.T.sub.l] and [[??].sup.R.sub.m] have magnitude modeled by functions of [[phi].sup.T.sub.l] and [[phi].sup.R.sub.m], respectively; that is,

[parallel][[??].sup.k.sub.l][parallel] = [G.sub.k] ([[phi].sub.l]), [G.sub.k] :[-[pi], [pi]) [right arrow] [R.sup.+] (17)

for l [member of] {1, 2, ..., card{[S.sub.k]}} and k [member of] {T, R}.

(vi) The aforementioned gains, phases, AODs, and AOAs are mutually independent r.v.

Assumptions (i)-(iv) and (vi) are customary in the state of the art, but the one introduced in (v) is not. This latter assumption provides a generic description of the IOs' position and is introduced in this paper with the purpose of obtaining insights into the channel model's statistics that are not limited to a specific geometrical configuration of the propagation area. The fact that the magnitude of [[??].sup.T.sub.l] and [[??].sup.R.sub.m] is modeled as functions of [[phi].sup.T.sub.l] and [[phi].sup.R.sub.m] implies that the IOs are located on the contour of two surfaces. The shape of such surfaces will be determined by the functions [G.sub.T] and [G.sub.R].

3.2. First-Order Statistical Analysis. The instantaneous mean value and average power of a doubly dispersive V2V fading channel are TF-varying quantities given as [m.sub.h](t; f) [??] E{H(t; f)} and [P.sub.h](t; f) [??] E{[[absolute value of H(t; f)].sup.2]}, where H(t; f) [??] [[integral].sup.[infinity].sub.-[infinity]] h(t; [tau]) exp{-j2[pi]f[tau]}d[tau] is the channel transfer function. By computing the Fourier transform of the CIR in (3) with respect to r, we find that

[mathematical expression not reproducible]. (18)

Under Assumptions (i)-(vi), we find [m.sub.h](t; f) = 0 and [P.sub.h](t; f) = [[sigma].sub.H] = [[zeta].sub.T][[zeta].sub.R], for t [member of] [0, [T.sub.0]], where

[mathematical expression not reproducible]. (19)

The previous equation indicates that the average power of the channel is just a cumulus of the average powers of the received multipath signal's components (as was to be expected). This equation also shows that the channel's large-scale variations are not taken into account by the model in (18), as the instantaneous average power of H(t; f) is a TF-invariant quantity for t [member of] [0, [T.sub.0]]. We have intentionally neglected such variations to focus the spotlight of our work on the nonstationarities of V2V channels arising from the time-varying nature of the propagation delays. However, the large-scale variations can easily be incorporated, for example, by modeling the gains as random functions of time (i.e., stochastic processes) [g.sup.k.sub.l](t), l [member of] {1, 2, ..., card([S.sub.k])}, k [member of] {T, R}, each having a time-varying average power.

The envelope and the phase of H(t; f) are also TF-varying quantities given by [XI](t,f) [??] [absolute value of H(t; f)] and [PSI](t, f) [??] arg{H(t; f)}. The computation of the distributions of [XI](t, f) and [PSI](t, f) is not a trivial task, because the term [parallel][p.sup.S.sub.l,m][parallel] in (13) is a nonlinear transformation of the random vectors [[??].sup.T.sub.l] and [[??].sup.R.sub.m] (see (11)). However, if the condition D [much greater than] max{[parallel][[??].sup.T.sub.l][parallel], [parallel][[??].sup.R.sub.m][parallel]} holds for all l, m, then this transformation becomes linear and equal to [parallel][p.sup.S.sub.l,m][parallel] [approximately equal to] [D.sup.T.sub.l] + [D.sup.R.sub.m], where = [D.sup.k.sub.l] = D/2 + [q.sub.k] <D, [[??].sup.k.sub.l]>/D, for l [member of] {1, 2, ..., card{[S.sub.k]}}, k [member of] {T, R}, and

[mathematical expression not reproducible]. (20)

With this simplification, the FO PDF of [XI](t, f) can be computed following the procedures presented in [28, 29]. Thereby, the FO PDF of [XI](t, f) is given by

[mathematical expression not reproducible] (21)

for [zi] [greater than or equal to] 0 and t [member of] [0, [T.sub.0]], where [J.sub.0](*) is the Bessel function of the first kind and zeroth order. Analogously, the PDF of [PSI](t, f) is given as

[P.sub.[psi]] (t, f; [psi]) = 1/2 [pi], 0 < [psi] [less than or equal to] 2[pi], t [member of] [0, [T.sub.0]]. (22)

The details on the derivations are omitted for reasons of brevity, but an outline can be found in [30].

We can observe from (21) and (22) that [XI](t, f) and [PSI](t, f) are FO-stationary random processes, as their PDFs [p.sub.[XI]](t, f; [xi]) and [p.sub.[psi]](t, f; [psi]) do not change over time and frequency within the relevant interval t [member of] [0, [T.sub.0]]. This feature indicates that the nonstationary characteristics of the proposed channel model--if any--are not caused by large-scale factors, such as shadowing or path loss, or by gross changes in the location of the IOs, such as the appearance and disappearance of IOs. We note again that the influence of such large-scale factors has intentionally been neglected in this paper to place emphasis on the nonstationarities caused by the time-varying propagation delays.

Unlike the PDF of the phase [PSI](t, f), which does not depend on the number of multipath components L x M, the PDF of the envelope [XI](t, fi) is highly influenced by the values of L and M. Nevertheless, if the gains [g.sup.k.sub.l] are Rayleigh distributed r.v., that is, if

[mathematical expression not reproducible], (23)

for l [member of] {1, 2, ..., card([S.sub.k])} and k [member of] {T, R}, where the parameters [[rho].sup.k.sub.l] control the dispersion of the Rayleigh distribution, and they are not necessarily equal to each other, then the PDF in (21) can be simplified to

[P.sub.[XI]]([xi]) = 4[xi]/[[sigma].sub.H] [K.sub.0] (2[xi]/[[sigma].sub.H]), [xi] [greater than or equal to] 0, (24)

for t [member of] [0, [T.sub.0]], where [K.sub.0] is the modified Bessel function of the second kind and zeroth order (see [30] for details on the derivation of (24)). The PDF in (24) is known as the double-Rayleigh distribution [13]. The same result is obtained regardless of the distribution of the gains [g.sup.k.sub.l] if the value of L is large (infinitely large in theory).

3.3. Four-Dimensional TF-CF. For the analysis of the correlation properties of H(t; f), we will consider the following definition of the 4D TF-CF given in [31, Eq. (3b)]:

[R.sub.H] (t, f; [DELTA]t, [DELTA]f) [??] E {[H.sup.*] (t - [DELTA]t; f) H(t; f + [DELTA]f)}. (25)

Starting from (18) and (25) and assuming that [t.sub.0] = 0, we show in Appendix A that

[mathematical expression not reproducible], (26)

where [mathematical expression not reproducible], while [[phi].sub.k] is an arbitrary r.v. in [mathematical expression not reproducible], and

[f.sub.D] ([[phi].sub.T], [[phi].sub.R]) = [f.sub.T] ([[phi].sub.T]) + [f.sub.R] ([[phi].sub.R]). (27)

The Doppler frequencies [f.sub.T]([[phi].sub.T]) and [f.sub.R]([[phi].sub.R]) are equal to

[f.sub.k] ([[phi].sub.k]) = [f.sup.k.sub.max] cos ([[phi].sub.k] - [[gamma].sub.k]), k [member of] {T, R}, (28)

where [[gamma].sub.k] is the angle of the velocity vector [v.sub.k], for k [member of] {T, R}, [f.sup.k.sub.max] = [v.sub.k]/[lambda], with [v.sub.k] denoting the speed of [T.sub.X] (when k = T) and [R.sub.X] (when k = R). In turn, the function [G.sub.S]([[phi].sub.T], [[phi].sub.R]) in (26) is given as

[G.sub.S] ([[phi].sub.T], [[phi].sub.R])

= [parallel]D - [G.sub.T] ([[phi].sub.T]) [u.sub.T] ([[phi].sub.T]) + [G.sub.R] ([[phi].sub.R]) [u.sub.r] ([[phi].sub.R])[parallel], (29)

where [u.sub.k]([[phi].sub.k]) is a unit vector that points at the direction specified by [[phi].sub.k], for k [member of] {T, R}. Invoking the expected value theorem, we can express [R.sub.H](t, f; [DELTA]t, [DELTA]f) as

[mathematical expression not reproducible]. (30)

The result in (30) can be simplified to the product of two line integrals if the condition D [much greater than] max{[G.sub.T]([[phi].sub.T]), [G.sub.R]([[phi].sub.R])} holds. Under this condition,

[G.sub.S] ([[phi].sub.T], [[phi].sub.R]) [approximately equal to] [D.sub.T] ([[phi].sub.T]) + [D.sub.R] ([[phi].sub.R]), (31)

where

[D.sub.k] ([[phi].sub.k]) = D/2 + [q.sub.k][G.sub.k] ([[phi].sub.k]) <D, [u.sub.k] ([[phi].sub.k])>/D, k [member of] {T, R}. (32)

Thereby, we have

[R.sub.H] (t, f; [DELTA]t, [DELTA]f) [approximately equal to] [R.sup.T.sub.H] (t, f; [DELTA]t, [DELTA]f) [R.sup.R.sub.H] (t, f; [DELTA]t, [DELTA]f) (33)

with

[mathematical expression not reproducible]. (34)

Equations (33) and (34) show that the 4D TF-CF of the proposed channel model can be factorized as the Kronecker product (tensor) of two different TF-CFs, one due to the signal's dispersion at the transmitter side ([R.sup.T.sub.H](t, f; [DELTA]t, [DELTA]f)) and the other to the dispersion at the receiver side ([R.sup.R.sub.H](t, f; [DELTA]t, [DELTA]f). This is not a novel result, nevertheless, as it has been shown in previous papers that the channel's correlation function can be expressed by the Kronecker product if the transmitted signal arrives at the receiver via double interactions with IOs [32, 33]. However, to the best of the authors' knowledge, the GBSMs that have been proposed for doubly and triply selective V2V channels assuming the propagation of plane waves and a double interaction with IOs fulfill the WSSUS condition by design [24-26]. In contrast, we can observe from (26), (30), (33), and (34) that the proposed channel model is a non-WSSUS random process, because its 4D TF-CF is a TF-varying function, meaning that ([R.sub.H]([t.sub.1], [f.sub.1]; [DELTA]t, [DELTA]f) [not equal to] ([R.sub.H]([t.sub.2], [f.sub.2]; [DELTA]t, [DELTA]f) for different observation instants ([t.sub.1], [f.sub.2]) and ([t.sub.2], [f.sub.2]). We recall that a TF-dispersive channel is said to fulfill the WSSUS condition if its mean value is a constant, and its TF-CF is invariant over the time t and the frequency f variables, that is, if ([R.sub.H]([t.sub.1], [f.sub.1]; [DELTA]t, [DELTA]f) = ([R.sub.H]([t.sub.2], [f.sub.2]; [DELTA]t, [DELTA]f), for ([t.sub.1], [f.sub.1]) [not equal to] ([t.sub.2], [f.sub.2]) [5, 34]. In addition, we recall that the fulfillment of the WSSUS condition implies that the channel is simultaneously WSS in the time and the frequency domains [5].

The nonstationary features of our channel model are noteworthy because our modeling framework does not consider macroscopic (large-scale) or microscopic (small-scale) factors that are already well-known sources of nonstationarities, such as path loss, shadowing, the appearance and disappearance of IOs, or time-varying AODs and AOAs. The nonstationarities of the channel model defined by (3) stem from the propagation delays' temporal variations. This microscopic source of nonstationarities has passed practically unnoticed in the literature, and its analysis has therefore received little attention.

Equation (30) shows that the time dependence of [R.sub.H](t, f; [DELTA]t, [DELTA]f) is given by the factor t * [DELTA]f * [f.sub.D]([[phi].sub.T], [[phi].sub.R])/[f.sub.c], whereas its frequency dependence is given by f * [DELTA]t * [f.sub.D]([[phi].sub.T], [[phi].sub.R])/[f.sub.c]. These two factors must be equal to zero in order for the proposed channel model to be WSS in both the time t and the frequency f variables (i.e., to enforce the WSSUS condition). Hence, this condition can only be met if the Doppler frequency [f.sub.D]([[phi].sub.T], [[phi].sub.R]) is equal to zero, meaning that the channel is static. Nevertheless, a static channel model is of little relevance for the analysis of V2V communication systems.

While the proposed channel model does not fulfill the WSSUS condition, it is worth noting that [R.sub.H](t, f; [DELTA]t, [DELTA]f) becomes a time-independent function over the observation region defined by Y(t, [DELTA]t) if we make [DELTA]f = 0. For that particular case, [R.sub.H](t, f; [DELTA]t; [DELTA]f)[|.sub. [DELTA]j=0] = Y(t, [DELTA]t)[T.sub.H](f; [DELTA]t), where [T.sub.H](f; [DELTA]t) is a frequency-varying time correlation function (TCF) given as

[mathematical expression not reproducible]. (35)

This means that, for a constant value of f, say f = f', the channel transfer function H(t; fi') can be deemed WSS in the time domain within the relevant observation interval t [member of] [0, [T.sub.0]]. Analogously, [R.sub.H](t, f; [DELTA]t, [DELTA]f) can be simplified to a frequency-independent function if [DELTA]t = 0. Under such conditions, [mathematical expression not reproducible], where [F.sub.H](t; [DELTA]f) is a time-varying frequency correlation function (FCF) equal to

[mathematical expression not reproducible]. (36)

Thus, for a constant value of t, say t = t', the channel transfer function H(t'; f) can be modeled as a WSS process in the frequency domain. We note that the Kronecker factorization holds also for both the frequency-varying TCF [T.sub.H](f; [DELTA]t) and the time-varying FCF [F.sub.H](t; [DELTA]f).

Even though the proposed channel model is a non-WSSUS random process, the two particular cases discussed above show that H(t;f) fulfills the WSS condition in one dimension (either time or frequency) if the channel's dispersion is neglected in the other dimension. This characteristic makes our proposal compatible with some notable channel models that are widely accepted benchmarks for the performance analysis of wireless communication systems. Consider, for example, the model proposed by Clarke in [9] for WSS frequency-nonselective (narrowband) Rayleigh fading channels.

The mathematical definition of wide-sense stationarity requires the first- and second-order (SO) statistical properties of a random process to be invariant over the corresponding index set (the time or frequency lines, or the TF plane). While difficult to meet, the condition of statistical invariance should always be the yardstick to determine whether a random process is WSS or not. Then, if the process is found to be nonstationary, we should turn our attention to the concept of quasi-stationarity to obtain a less rigid notion of its stationarity properties [35]. Even though a thorough analysis of quasi-stationarity is beyond the scope of this paper, some interesting remarks can be made in that regard from (26). For example, in practice, the values of f are restricted to f [member of] [-B/2,B/2], where B is the signal's bandwidth. Given that B is typically much smaller than the carrier frequency [f.sub.c], the factor f * [DELTA]t * [f.sub.D]([[phi].sub.T], [[phi].sub.R])/[f.sub.c], which determines the dependence of [R.sub.H](t, f; [DELTA]t, [DELTA]f) on the frequency variable f, can be approximated as [DELTA]t * [f.sub.D]([[phi].sub.T], [[phi].sub.R])f/[f.sub.c] [approximately equal to] 0. This means that the channel can be modeled by a quasi-WSS random process in the frequency domain if we consider a frequency observation interval of a length similar to the signal's bandwidth and B [much less than] [f.sub.c]. On the other hand, the dependence of [R.sub.H](t, f; [DELTA]t, [DELTA]f) on the time variable t is given by the term t * [DELTA] * [f.sub.D]([[phi].sub.T], [[phi].sub.R])/[f.sub.c], which is influenced by the degree of mobility of [T.sub.X] and [R.sub.X] (through the Doppler frequency shift [f.sub.D]([[phi].sub.T], [[phi].sub.R])), the signal's bandwidth (through the frequency lag [DELTA]f [member of] [-B, B]), and the carrier frequency [f.sub.c]. Hence, the channel maybe modeled by a quasi-WSS random process in the time domain over the region associated with Y(t, [DELTA]t) if the MSs move at very low speeds (such that [f.sub.D]([[phi].sub.T], [[phi].sub.R]) [approximately equal to] 0), or if the signal's bandwidth B is much smaller than [f.sub.c], such that [f.sub.D]([[phi].sub.T], [[[phi].sub.R])[DELTA]f/[f.sub.c] [approximately equal to] 0, [for all][DELTA]f [member of] [-B, B]. Nonetheless, to properly assess the quasi-stationarity properties of the proposed channel model, a more formal and mathematically rigorous analysis is required.

3.4. TF-Dependent Delay and Doppler Profiles. The spectral properties of the proposed channel model can be analyzed on the grounds of the TF-dependent delay and Doppler profiles. These two functions are given, respectively, as [31]

[mathematical expression not reproducible] (37)

[mathematical expression not reproducible] (38)

[mathematical expression not reproducible] (39)

[mathematical expression not reproducible]. (40)

By substituting (36) and (35) into (38) and (40), we find

[mathematical expression not reproducible], (41)

where sinc(x) [??] sin(x)/x is the sinc function. Equation (41) shows that the delay profile of the proposed channel model is a frequency-invariant function [P.sub.H](t; [tau]), which depends on the channel's angular statistics and the propagation area geometry. The Doppler profile, on the other hand, is a TF-dependent function that depends only on the channel's angular statistics.

4. Particular Case: The Geometrical Two-Ring Scattering Model

4.1. Background. The definite integrals in (30), (33), and (34) can be evaluated numerically without major difficulties if [G.sub.T], [G.sub.R], [p.sup.T.sub.[phi]] and [p.sup.R.sub.[phi]] are known. In fact, the 4D TF-CF can be expressed in a closed form if the abovementioned functions allow writing the integral at the right-hand side of (34) in terms of standard functions. This is demonstrated in this section by considering the particular case of the geometrical two-ring scattering model and the von Mises distribution of the AOD and AOA. The geometrical two-ring model has widely been employed as a reference to analyze the correlation properties of F2M and V2V fading channels assuming the fulfillment of the WSSUS condition [24, 33]. In this geometrical model, [T.sub.X] and [R.sub.X] are each surrounded by a ring of IOs, in such a way that

[G.sub.T] ([[phi].sub.T]) = [r.sub.T], [for all][[phi].sub.T] (43a)

[G.sub.R] ([[phi].sub.R]) = [r.sub.R], [for all][[phi].sub.R], (43b)

where [r.sub.k], k [member of] {T, R}, is the radius of the ring of IOs surrounding [T.sub.X] (when k = T) and [R.sub.X] (when k = R). An illustration of this model is shown in Figure 2.

The von Mises distribution, on the other hand, is a well-known model for the statistics of circular data and is analogous to the Gaussian distribution on the line [36]. Following the von Mises distribution, we define

[p.sup.k.sub.[phi]] ([phi]) = exp {[[kappa].sub.k] cos([phi] - [[mu].sub.k])}/2[pi][I.sub.0] ([[kappa].sub.k]), k [member of] {T, R}, (44)

where [[mu].sub.k] and [[kappa].sub.k] are the distribution's mean value and concentration parameter and [I.sub.0] is the modified Bessel function of the first kind and zeroth order.

4.2. Closed-Form Solution of the 4D TF-CF. Without loss of generality, we will assume that the time-invariant vectors [O.sub.T] and [O.sub.R] introduced in Figure 1 are collinear with the x-axis, such that D = [O.sub.R] - [O.sub.T] = [[D, 0].sup.[dagger]]. In addition, we will assume that D [much greater than or equal to] max{[r.sub.T], [r.sub.R]}. Under these conditions, and substituting (43a), (43b), and (44) into (34), we show in Appendix B that

[mathematical expression not reproducible], (45)

where

[B.sup.k.sub.c] = [f.sup.k.sub.max] cos ([[gamma].sub.k]) Z (t, f, [DELTA]t, [DELTA]f) - [q.sub.k] [DELTA]f[r.sub.k]/[C.sub.l] (46a)

[B.sup.k.sub.s] = [f.sup.k.sub.max] sin ([[gamma].sub.k]) Z (t, f, [DELTA]t, [DELTA]f) (46b)

for k [member of] {T, R}, and

Z(t, f, [DELTA]t, [DELTA]f) = [DELTA]t([f.sub.c] + f/[f.sub.c]) + [DELTA]f t/[f.sub.c]. (47)

4.3. Numerical Examples. In the remainder of this section, we present some graphical examples of the 4D TF-CF [R.sub.H](t, f; [DELTA]t, [DELTA]f) in (45) and the corresponding time-varying delay profile [P.sub.H](t; [tau]) and TF-dependent Doppler profile [D.sub.H](t, f; v). The expressions presented in Section 3 for the FO PDFs of the envelope and phase are not considered here because they have already been analyzed in [30]. For the evaluation of [R.sub.H](t, f; [DELTA]t, [DELTA]f), [P.sub.H](t; [tau]), and [D.sub.H](t, f; v), we have chosen [f.sub.c] = 5.9 GHz, B = 10 MHz, [T.sub.0] = 6.4 ms, [[gamma].sub.T] = 60[degrees], [[gamma].sub.R] = 250[degrees], D = 500 m, [r.sub.T] = 30 m, [r.sub.R] = 40 m, [f.sup.T.sub.max] = 500 Hz (corresponding to a speed of 91.4 km/h), [f.sup.T.sub.max] = 500 Hz, [[sigma].sub.H] = 1, [[kappa].sub.T] = 1 (moderate nonisotropic dispersion), [[kappa].sub.R] = 10 (highly nonisotropic dispersion), [[mu].sub.T] = 60[degrees], and [[mu].sub.R] = 120[degrees]. The system related parameters [f.sub.c], B, and [T.sub.0] are taken from the IEEE 802.11p standard for vehicular DSRC systems [3] (the value of [T.sub.0] equals the duration of a signal frame comprising 800 data symbols, each having a length of 8 [mu]s).

4.3.1. 4D TF-CF. A 3D graph and a contour plot of the absolute value of [R.sub.H](t, f; [DELTA]t, [DELTA]f) are presented in Figure 3 for an arbitrary observation point (t, f) = (0.5[T.sub.0], 0.3B) in the TF plane. These two graphs show that the proposed channel model is a non-WSSUS random process, because its 4D TF-CF has an asymmetrical shape. We recall that a fundamental property of WSS random processes is that their autocorrelation function (ACF) is symmetric around the origin [37, Theorem 10.12]. Particularly, for a complex-valued one-dimensional (1D) random process x(f) defined in either the time or the frequency domain, the WSS condition implies that [R.sup.*.sub.[chi]]([eta]; [DELTA][eta]) = [R.sub.[ci]]([eta]; - [DELTA][eta]), [for all][eta], [DELTA][eta] [member of] R, where [R.sub.[chi]] ([eta]; [DELTA][eta]) [??] E{[[chi].sup.*]([eta])[chi]([eta] - [DELTA][eta])} is the ACF of [chi]([eta]). Analogously, for a complex-valued 2D random process [chi](t; f) defined in the TF plane, it can be shown that the WSSUS condition implies that [R.sub.[chi]](t, f; [DELTA]t, [DELTA]f) = [R.sup.*.sub.[chi]](t, f; [DELTA]t], - [DELTA]f) = [R.sup.*.sub.[chi]](t, f; [DELTA]t, [DELTA]f) = [R.sub.[chi]](t, f; -[DELTA]t, -[DELTA]f), [for all]t, f, where [R.sub.[chi]](t, f; [DELTA]t, [DELTA]f) is the 4D TF-CF defined in (25).

To make the asymmetries of [R.sub.H](t, f; [DELTA]t, [DELTA]f) more evident, Figure 4 shows curves of the absolute value of [R.sub.H](t, f; [DELTA]t, [DELTA]f) generated for fixed values of either the frequency lag [DELTA]f or the time lag [DELTA]t. The curves shown in Figure 4(a) were computed for [DELTA]f [member of] {-2.5 x [10.sup.6], -1.25 x [10.sup.6], 0,1.25 x [10.sup.6], 2.5 x [10.sup.6]}, while the ones presented in Figure 4(b) were generated for [DELTA]t [member of] {-8 x [10.sup.-4], -3 x [10.sup.-4], 0,3 x [10.sup.-4], 8 x [10.sup.-4]}. The graphs obtained for [DELTA]f = 0 and [DELTA]t = 0 correspond to the particular cases of the frequency-varying TCF [T.sub.H](f; [DELTA]t) and the time-varying FCF [F.sub.H](t; [DELTA]f), respectively.

Aside from the particular cases of [T.sub.H](f; [DELTA]t) and [F.sub.H](t; [DELTA]f), we can observe from Figure 4 that [absolute value of [R.sub.H](t, f; [DELTA]t, [DELTA]f)] is an asymmetrical function. However, the asymmetries of [absolute value of [R.sub.H](t, f; [DELTA]t, [DELTA]f)] smoothen off gradually as the channel becomes nondispersive in one dimension. This behavior shows that the proposed channel model's quasi-stationary characteristics are stronger as the 2D TF lag (At, Afi) approaches the ideal pairings that lead to [T.sub.H](f; [DELTA]t) and [F.sub.H](t; [DELTA]f). On the other hand, the symmetrical shape of [T.sub.H](f; [DELTA]t) and [F.sub.H](t; [DELTA]f) can be taken as graphical evidence for the observation we made in the previous section about the proposed channel model becoming a WSS random process in one dimension (time or frequency) if its dispersiveness is neglected in the other dimension. Even though these two ideal cases have limited practical significance, they are of theoretical relevance, as they allow identifying scenarios where the WSS condition is met in one dimension.

4.3.2. Time-Frequency Dependent Delay Profile. A 3D graph and a contour plot of the time-varying delay profile [P.sub.H](T; [tau]) are presented in Figure 5. To highlight the time-varying nature of [P.sub.H](t; [tau]), we have considered a larger observation time window of length [T.sub.0] = 320 ms. In our simulation setup, the MSs are moving towards each other along nearly parallel trajectories. This explains the shortening of the minimum propagation delay, [[tau].sub.min](t), which is defined as the propagation delay of the plane wave that travels to [R.sub.X] over the shortest path. Note that [[tau].sub.min](t) shrinks from 1.6678 [micro]s at t = 0 to 1.6457 [micro]s at t = 320 ms. On the other hand, the maximum propagation delay, [[tau].sub.max](t), which is defined as the propagation delay of the plane wave that travels to [R.sub.X] over the largest path, increases from 2.1127 ms to around 2.15 ms. This increment is due to the fact that as time evolves, the MSs get closer to some IOs but farther off from others.

4.3.3. Time-Frequency Dependent Doppler Profile. Finally, Figure 6 shows a 3D graph and a contour plot of the TF-varying Doppler profile [D.sub.H](t, f; v) evaluated at t = 0.8[T.sub.0]. While the variability of [D.sub.H](t, f; v) in the frequency f variable is not evident in the figure, the shape of [D.sub.H](t, f; v) does vary in the frequency domain, but at a rate that is much smaller than the length of the observation window f [member of] [-B/2, B/2]. This is demonstrated in Figure 7, where we present the same contour plot of Figure 6(b) but for an expanded observation window f [member of] [-100B, 100B]. Regarding the time variations of [D.sub.H](t, f; v), these are the result of convolving the spectra of the time-invariant TCF [T.sub.H](f; [DELTA]t) and the time-dependent window function [mathematical expression not reproducible]. Such temporal variations manifest themselves by a rippling effect in the spectrum of [T.sub.H](f; [DELTA]t). Our numerical experiments indicate that this effect is rather mild, making the time-varying nature of [D.sub.H](t, f; v) almost imperceptible. This is demonstrated by comparing the 3D surface of Figure 6(a) with the one presented in Figure 8, where the value of the observation time t is equal to t = 0.1[T.sub.0]. Hence, based on our numerical results, we conclude that, for most practical purposes, the Doppler profile of the proposed channel model can be considered a time-invariant function [D.sub.H](f; v). However, the frequency variations of [D.sub.H](f; v) should in general not be neglected.

5. Conclusions

In this paper, we have proposed a novel GBSM for small-scale non-WSSUS V2V double-Rayleigh fading channels. The model is quite flexible and can easily be adapted to a wide range of 2D and 3D geometrical propagation scenarios. Based on this model, we derived general expressions for the envelope and phase distributions, the 4D TF-CF, and the TF-dependent delay and Doppler profiles. The obtained expressions show that the 4D TF-CF of the proposed channel model is a TF-varying function. However, such expressions also show that the corresponding delay profile is invariant in the frequency domain, while the Doppler profile can be deemed invariant in the time domain for most practical purposes. Even though the proposed model does not fulfill the WSS condition simultaneously in the time and the frequency domains, it becomes a WSS random process in one dimension (time or frequency) if the channel's dispersiveness is neglected in the other dimension. For the more relevant case of 2D dispersion, our results suggest that the proposed channel model is quasi-WSS over small observation regions in the TF plane. However, a more formal analysis is necessary to assess the model's quasi-stationarity characteristics.

Appendix

A. Derivation of (26)

By a direct evaluation of (25), and taking into account the definition of the propagation delays [[tau].sub.l,m](t) given by (13), we find

[mathematical expression not reproducible] (A.1)

The propagation delays [[tau].sub.l,m] (t) are statistically equivalent random processes for all l [member of] {1, 2, ..., card{[S.sub.T]}} and m [member of] {1, 2, ..., card{[S.sub.T]}}. For this reason, we can take the third expectation out of the double summation in the latest equation. Then, taking account of (19), we get to the result presented in (26).

B. Derivation of (45)

Assuming that [G.sub.T]([[phi].sub.T]) = [r.sub.T], [for all][[phi].sub.T], [G.sub.T] = [r.sub.R], [for all] [[phi].sub.R], and D = [[D, 0].sup.[dagger]], we can express (32) as

[D.sub.k] ([[phi].sub.k]) = D/2 + [q.sub.k][r.sub.k] cos ([[phi].sub.k]), k [member of] {T, R}. (B.1)

Substituting this result into (34), we find

[mathematical expression not reproducible]. (B.2)

After some simple algebraic manipulations, we obtain

[mathematical expression not reproducible], (B.3)

where [B.sup.k.sub.c] and [B.sup.k.sub.s] are defined in (46a) and (46b). Substituting the von Mises PDF into the above equation and making use of [38, Eq. (3.338-4)], we have

[mathematical expression not reproducible]. (B.4)

Finally, from (33), we obtain the expression presented in (45) for [R.sub.H](t, f; [DELTA]t, [DELTA]f).

https://doi.org/10.1155/2017/6394653

Disclosure

The material in this paper is based on "Modeling of Nonstationary Double-Rayleigh Fading Channels for Mobile-to-Mobile Communications," by C. A. Gutierrez, J. M. Luna-Rivera, and D. U. Campos-Delgado, which appeared in the proceedings of the 2016 European Wireless Conference (EW2016), Oulu, Finland, March 2016 ([c] 2016 IEEE).

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors acknowledge the financial support provided by CONACYT for the development of this research work: Project Grants nos. 236188 and 241272.

References

[1] G. Araniti, C. Campolo, M. Condoluci, A. Iera, and A. Molinaro, "LTE for vehicular networking: a survey," IEEE Communications Magazine, vol. 51, no. 5, pp. 148-157, 2013.

[2] G. Karagiannis, O. Altintas, E. Ekici et al., "Vehicular networking: a survey and tutorial on requirements, architectures, challenges, standards and solutions," IEEE Communications Surveys and Tutorials, vol. 13, no. 4, pp. 584-616, 2011.

[3] J. B. Kenney, "Dedicated short-range communications (DSRC) standards in the United States," Proceedings of the IEEE, vol. 99, no. 7, pp. 1162-1182, 2011.

[4] P. Bello, "Characterization of randomly time-variant linear channels," IEEE Transactions on Communications, vol. 11, no. 4, pp. 360-393, 1963.

[5] H. Schulze and C. LUders, Theory and Applications of OFDM and CDMA, John Wiley & Sons, Ltd, Chichester, UK, 2005.

[6] A. Molisch, Wireless Communications, John Wiley and Sons, Chichester, UK, 2005.

[7] N. Costa and S. Haykin, Multiple-Input, Multiple-Output Channel Models: Theory and Practice, Wiley, New York, USA, 2010.

[8] C. F. Mecklenbrauker, A. F. Molisch, J. Karedal et al., "Vehicular channel characterization and its implications for wireless system design and performance," Proceedings of the IEEE, vol. 99, no. 7, pp. 1189-1212, 2011.

[9] R. H. Clarke, "A statistical theory of mobile-radio reception," The Bell System Technical Journal, vol. 47, no. 6, pp. 957-1000, 1968.

[10] O. Renaudin, V.-M. Kolmonen, P. Vainikainen, and C. Oestges, "Non-stationary narrowband MIMO inter-vehicle channel characterization in the 5-GHz band," IEEE Transactions on Vehicular Technology, vol. 59, no. 4, pp. 2007-2015, 2010.

[11] L. Bernado, T. Zemen, F. Tufvesson, A. F. Molisch, and C. F. Mecklenbrauker, "The (in-) validity of the WSSUS assumption in vehicular radio channels," in Proceedings of the IEEE 23rd International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC '12), pp. 1757-1762, Sydney, NSW, Australia, September 2012.

[12] A. Roivainen, P. Jayasinghe, J. Meinilau, V. Hovinen, and M. Latva-Aho, "Vehicle-to-vehicle radio channel characterization in urban environment at 2.3 GHz and 5.25 GHz," in Proceedings of the 2014 25th IEEE Annual International Symposium on Personal, Indoor, and Mobile Radio Communication, IEEE PIMRC 2014, pp. 63-67, Washington, DC, USA, September 2014.

[13] D. W. Matolak and J. Frolik, "Worse-than-rayleigh fading: experimental results and theoretical models," IEEE Communications Magazine, vol. 49, no. 4, pp. 140-146, 2011.

[14] L. Rubio, J. Reig, V. M. Rodrigo-Penarrocha, H. Fernandez, and S. Loredo, "Analysis of small-scale fading distributions in vehicular-to-vehicular communications," Mobile Information Systems, vol. 2016, Article ID 9584815, 7 pages, 2016.

[15] P. S. Bithas, K. Maliatsos, and A. G. Kanatas, "The Bivariate Double Rayleigh Distribution for Multichannel Time-Varying Systems," IEEE Wireless Communications Letters, vol. 5, no. 5, pp. 524-527, 2016.

[16] B. TalhaandM. Patzold, "A geometrical three-ring-based model for MIMO mobile-to-mobile fading channels in cooperative networks," Eurasip Journal on Advances in Signal Processing, vol. 2011, Article ID 892871, 2011.

[17] J. Karedal, F. Tufvesson, N. Czink et al., "A geometry-based stochastic MIMO model for vehicle-to-vehicle communications," IEEE Transactions on Wireless Communications, vol. 8, no. 7, pp. 3646-3657, 2009.

[18] M. Walter, D. Shutin, and U.-C. Fiebig, "Delay-dependent doppler probability density functions for vehicle-to-vehicle scatter channels," IEEE Transactions on Antennas and Propagation, vol. 62, no. 4, pp. 2238-2249, 2014.

[19] Y. Yuan, C.-X. Wang, Y. He, M. M. Alwakeel, and E.-H. M. Aggoune, "3D Wideband Non-Stationary Geometry-Based Stochastic Models for Non-Isotropic MIMO Vehicle-to-Vehicle Channels," IEEE Transactions on Wireless Communications, vol. 14, no. 12, pp. 6883-6895, 2015.

[20] C. A. Gutierrez, M. Patzold, W. Dahech, and N. Youssef, "A Non-WSSUS Mobile-to-Mobile Channel Model Assuming Velocity Variations of the Mobile Stations," in Proceedings of the 2017IEEE Wireless Communications and Networking Conference (WCNC), pp. 1-6, San Francisco, CA, USA, March 2017.

[21] J. T. Gutierrez-Mena, C. A. Gutierrez, J. M. Luna-Rivera, D. U. Campos-Delgado, and J. V. Castillo, "Geometry-based statistical modeling of non-stationary MIMO vehicle-to-vehicle channels," in Proceedings of the 5th ACM Symposium on Development and Analysis of Intelligent Vehicular Networks and Applications, DIVANet 2015, pp. 85-92, Cancun, Mexico, November 2015.

[22] C. A. Gutierrez, M. Luna-Rivera, and D. U. Campos-Delgado, "Modeling of non-stationary double-Rayleigh fading channels for mobile-to-mobile communications," in Proceedings of the 22nd European Wireless Conference, EW2016, pp. 131-136, Oulu, Finland, May 2016.

[23] W. Dahech, M. Patzold, C. A. Gutierrez, and N. Youssef, "A Non-Stationary Mobile-to-Mobile Channel Model Allowing for Velocity and Trajectory Variations of the Mobile Stations," IEEE Transactions on Wireless Communications, vol. 16, no. 3, pp. 1987-2000, 2017.

[24] A. G. Zajic and G. L. Stuber, "Three-dimensional modeling and simulation of wideband MIMO mobile-to-mobile channels," IEEE Transactions on Wireless Communications, vol. 8, no. 3, pp. 1260-1275, 2009.

[25] J. Chen and T. G. Pratt, "A three-dimensional geometry-based statistical model of 2x2 dual-polarized MIMO mobile-to-mobile wideband channels," Modelling and Simulation in Engineering, vol. 2012, Article ID 756508, 2012.

[26] X. Cheng, Q. Yao, M. Wen, C. Wang, L. Song, and B. Jiao, "Wideband channel modeling and intercarrier interference cancellation for vehicle-to-vehicle communication systems," IEEE Journal on Selected Areas in Communications, vol. 31, no. 9, pp. 434-448, 2013.

[27] D. K. Cheng, Field and Wave Electromagnetics, Addison-Wesley, New York, 2nd edition, 1989.

[28] M. Patzold and B. Talha, "On the statistical properties of sum-of-cisoids- based mobile radio channel simulators," in Proceedings of the Proc. 10th International Symposium on Wireless Personal Multimedia Communications (WPMC07, pp. 394-400, Jaipur, India, December 2007.

[29] B. O. Hogstad, C. A. Gutierrez, M. Patzold, and P. M. Crespo, "Classes of sum-of-cisoids processes and their statistics for the modeling and simulation of mobile fading channels," EURASIP Journal on Wireless Communications and Networking, vol. 2013, no. 1, pp. 1-15, 2013.

[30] J. J. Jaime-Rodriguez, C. A. Gutierrez, D. U. Campos-Delgado, and J. M. Luna-Rivera, "First-order statistics analysis of two new geometrical models for non-WSSUS mobile-to-mobile channels," in Proceedings of the 12th IEEE International Conference on Wireless and Mobile Computing, Networking and Communications, WiMob 2016, New York, NY, USA, October 2016.

[31] G. Matz, "On non-WSSUS wireless fading channels," IEEE Transactions on Wireless Communications, vol. 4, no. 5, pp. 2465-2478, 2005.

[32] A. G. Zajic, G. L. Stuber, T. G. Pratt, and S. T. Nguyen, "Wideband MIMO mobile-to-mobile channels: Geometry-based statistical modeling with experimental verification," IEEE Transactions on Vehicular Technology, vol. 58, no. 2, pp. 517-534, 2009.

[33] M. Patzold, B. O. Hogstad, and N. Youssef, "Modeling, analysis, and simulation of MIMO mobile-to-mobile fading channels," IEEE Transactions on Wireless Communications, vol. 7, no. 2, pp. 510-520, 2008.

[34] A. Papoulis and S. Pillai, Probability, Random Variables, and Stochastic Processes, McGraw-Hill, New York, NY, USA, 4th edition, 2002.

[35] R. He, O. Renaudin, V.-M. Kolmonen et al., "Characterization of quasi-stationarity regions for vehicle-to-vehicle radio channels," IEEE Transactions on Antennas and Propagation, vol. 63, no. 5, pp. 2237-2251, 2015.

[36] R. von Mises, "Uber die 'Ganzzahligkeit' der Atomgewichte und verwandte Fragen," Physikalische Zeitschrift, vol. 19, pp. 490-500, 1918.

[37] R. D. Yates and D. J. Goodman, Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers, John Wiley and Sons, New Jersey, 2nd edition, 2005.

[38] I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series and Products, Academic Press, New York, USA, 7th edition, 2007.

Carlos A. Gutierrez, (1) J. J. Jaime-Rodriguez, (1) J. M. Luna-Rivera, (1) Daniel U. Campos-Delgado, (1) and Javier Vazquez Castillo (2)

(1) Faculty of Science, Universidad Autonoma de San Luis Potosi, Av. Salvador Nava Martinezs/n, 78290 San Luis Potosi, SLP, Mexico

(2) Department of Engineering, Universidad de Quintana Roo, Blvd. Bahia Esq. Ignacio Comonfort s/n, 77019 Chetumal, QR, Mexico

Correspondence should be addressed to Carlos A. Gutierrez; cagutierrez@fc.uaslp.mx

Received 29 April 2017; Revised 14 July 2017; Accepted 1 August 2017; Published 3 October 2017

Academic Editor: Xianfu Lei

Caption: Figure 1: The reference propagation scenario at time t = [t.sub.0].

Caption: Figure 2: The geometrical two-ring model.

Caption: Figure 3: Absolute value of the 4D TF-CF [R.sub.H](t, f; [DELTA]t, [DELTA]f) at (t, f) = (0.5 [T.sub.0], 0.3B).

Caption: Figure 4: Absolute value of the 4D TF-CF [R.sub.H](t, f; [DELTA]t, [DELTA]f) at (t, f) = (0.5 [T.sub.0], 0.3B).

Caption: Figure 5: Absolute value of the time-varying delay profile [P.sub.H](t; [tau]) for [T.sub.0] = 320 ms.

Caption: Figure 6: Absolute value of the TF-varying Doppler profile [D.sub.H](t, f; v) for t = 0.8[T.sub.0] and [T.sub.0] = 320 ms.

Caption: Figure 7: Contour plot of the absolute value of the TF-varying Doppler profile [D.sub.H](t, f; v) for t = 0.8[T.sub.0], [T.sub.0] = 320 ms, and an extended observation window in the frequency domain.

Caption: Figure 8: 3D surface of the absolute value of the TF-varying Doppler profile [D.sub.H](t, f; v) for t = 0.1 [T.sub.0] and [T.sub.0] = 320 ms.

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article; wide-sense stationary uncorrelated scattering |
---|---|

Author: | Gutierrez, Carlos A.; Jaime-Rodriguez, J.J.; Luna-Rivera, J.M.; Campos-Delgado, Daniel U.; Castillo, |

Publication: | Wireless Communications and Mobile Computing |

Article Type: | Report |

Date: | Jan 1, 2017 |

Words: | 10438 |

Previous Article: | Compressed RSS Measurement for Communication and Sensing in the Internet of Things. |

Next Article: | REACH: An Efficient MAC Protocol for RF Energy Harvesting in Wireless Sensor Network. |

Topics: |