# A Novel 3D Nonstationary Channel Model Based on the von Mises-Fisher Scattering Distribution

Authors: Yuming Bi; Jianhua Zhang; Ming Zeng; Mengmeng Liu; Xiaodong Xu
Journal: Mobile Information Systems
 (2016)

---

## Abstract

In the last decade, the nonstationary properties of channel models have attracted more and more attention for many scenarios, that is, vehicle-to-vehicle (V2V), mobile-to-mobile (M2M), and high-speed train (HST). However, little research has been done on the real-physical channel model. In this paper, we propose a generalized three-dimensional (3D) nonstationary channel model, in which the scatterers are assumed to be distributed around the transmitter (Tx) and receiver (Rx) on a two-sphere model. By employing the von Mises-Fisher distribution, the mean values of the azimuth angle of departure (AAoD) and elevation angle of departure (EAoD) and the azimuth angle of arrival (AAoA) and elevation angle of arrival (EAoA) are tracked by time-variant (TV) Brownian Markov (BM) motion paths, which ensure the nonstationarity of the proposed channel model. Moreover, the TV autocorrelation function (ACF) and Doppler power spectrum density (DPSD) of the proposed nonstationary channel model are calculated by using signal processing tools, for example, fast Fourier transform (FFT) and short-time Fourier transform (STFT). In addition, the simulation results show that the TV scatterer distribution results in a nonstationary nonisotropic channel model, and the proposed model can be employed to simulate the 3D nonstationary channel model.


---

## Body

## 1. Introduction

The transmission channel is one of the most crucial parts in mobile communication systems, while the generation of an accurate and effective channel model for testing of communication systems represents a particular challenge. The majority of existing channel models in the literature rely on the assumption of wide-sense stationary (WSS) and uncorrelated scattering (US) [1, 2]. Under such an assumption, different transmission delays stay uncorrelated and the second-order channel moments are stationary. However, the physical channel in the real world is nonstationary [2–4]. Furthermore, plenty of empirical and analytical reports reveal that the WSSUS assumption is valid only for a short travelling distance. Therefore, in order to develop a more accurate and suitable channel model for the real-physical channel, it is essential to study nonstationarity channel modeling, which is a newly explored research field.The channel modeling can be classified into three categories, that is, deterministic, geometry, and correlation statistics [4]. In terms of deterministic channel modeling, Ai and He have applied their measurement data to V2V and HST scenarios employing the nonstationary channel modeling. In [2, 3], the stationary region especially has been extracted according to the correlated matrix distance. This is beneficial to elicit the extent of the channel variation if one could reconstruct the nonstationary channel simulators to approximate the real-physical channel. For geometry methods, [5] proposed a two-dimensional (2D) evolutionary spectrum approach to model nonstationary broadband mobile fading channels. In [6, 7], the nonstationary properties of massive MIMO channels have been characterized by cluster evolution (birth-death process) on both the array and time axes. Besides, in [8, 9], Pätzold suggested a nonstationary one-ring scattering channel model by utilizing Brownian random movement paths of the mobile station and deriving the angle of arrival (AoAs) and angle of departure (AoDs). In addition, [10] provided a generic method to find the MIMO correlation channel matrices under nonstationary interferences. However, on the one hand, all methods mentioned above only took into account the nonstationarity of the vertical plane while they ignored the impact of elevation plane; on the other hand, their modeling process required a high computational complexity and lacked flexibility. Therefore, it is essential to develop a simple and feasible 3D channel model for nonstationary research and to extend the 3D nonstationary channel model to capture the real-physical channel characteristics.In this paper, we propose a more general 3D geometric situation where the scatterers are assumed to be distributed around the user in a sphere area. By employing the von Mises-Fisher scattering distributions, the centers of the AoDs and AoAs scattering distribution are tracked by TV Brownian motion paths [11, 12], which could effectively reflect the movement process of the MS. Furthermore, the ACF of the complex channel gain and local PSD of the Doppler frequencies are calculated by nonstationary signal processing tools, for example, STFT and FFT [13, 14]. In addition, the results show that the time-variation of the scatterer distribution results in a nonstationary nonisotropic channel model and that the proposed model has the outstanding ability to describe vital statistical characteristics of nonstationary channel, such as the envelope distribution and the TV ACF and PSD [15]. Moreover, it also offers the advantages of low computational complexity and easy realization [8].The remainder of this paper is organized as follows: Section2 discusses a novel 3D channel model based on the geometrical two-sphere model. In Section 3, the von Mises-Fisher (VMF) distribution is employed to describe the 3D scatter distribution. A time-variant Brownian motion process is provided in Section 4. And then Section 5 proposes a novel 3D nonstationary channel model. The related simulation and numerical analysis are presented in Section 6. At last, the conclusion is drawn in Section 7.


## 2. A Novel 3D Channel Model

Let us consider a narrowband V2V communication channel with omnidirectional antennas. The radio propagation environment around the transmitter and the receiver is characterized by 3D nonisotropic scattering under line-of-sight (LoS) and none line-of-sight (NLoS) conditions. Figure1 illustrates the novel 3D geometrical channel model, in which the local scatterers are modeled on the surface of two spheres of radii R
T and R
R [16]. In order to limit the computational complexity of this model, both local scatterers around the Tx and Rx are considered. Therefore we assume that there are N
1 scatterers on the sphere (around the Tx) of radius R
T and the n
1th (n
1
=
1,2
,
…
,
N
1) scatterer is denoted by S
n
1. Similarly, around the Rx, N
2 scatterers lie on a surface of sphere of radius R
R, and the n
2th scatterer is denoted by S
n
2. The geometric symbols in Figure 1 have the following meaning: α
R
(
n
i
) and β
R
(
n
i
) denote the main AAoA and main EAoA of the n
ith scatterer, respectively. In the same way, α
T
(
n
i
) and β
T
(
n
i
) denote the main AAoD and main EAoD of the n
ith scatterer, respectively, and D is the distance between the Tx and Rx. It is worth noting that the distance D is assumed to be much larger than the radius R
R or R
T (i.e., R
T or R
R
≪
D). Therefore, the channel impulse response (CIR) of the two-sphere model at the carrier frequency f
c can be expressed as [17](1)
h
t
=
h
L
O
S
t
+
h
N
L
O
S
t
,where(2)
h
L
O
S
t
=
K
K
+
1
e
-
j
2
π
f
c
τ
1
×
e
j
2
π
f
T
m
a
x
t
cos
⁡
α
T
L
O
S
-
γ
T
·
cos
⁡
β
T
L
O
S
×
e
j
2
π
f
R
m
a
x
t
cos
⁡
α
R
L
O
S
-
γ
R
·
cos
⁡
β
R
L
O
S
.The CIR especially of two-sphere model under NLoS conditions is constructed as a sum of the single- and double-bounced rays with different energies. So it can be rewritten as(3)
h
N
L
O
S
t
=
∑
i
=
1
2
h
S
B
i
t
+
h
D
B
t
h
S
B
i
t
=
η
S
B
i
K
+
1
l
i
m
N
i
→
∞
∑
n
i
=
1
N
i
1
N
i
e
j
φ
n
i
-
2
π
f
c
τ
2
×
e
j
2
π
f
T
m
a
x
t
cos
⁡
α
T
n
i
-
γ
T
·
cos
⁡
β
T
n
i
×
e
j
2
π
f
R
m
a
x
t
cos
⁡
α
R
n
i
-
γ
R
·
cos
⁡
β
R
n
i
,
h
D
B
t
=
η
D
B
K
+
1
lim
N
1
,
N
2
→
∞
⁡
∑
n
1
,
n
2
=
1
N
1
,
N
2
1
N
1
N
2
e
j
φ
n
1
,
n
2
-
2
π
f
c
τ
3
×
e
j
2
π
f
T
m
a
x
t
cos
⁡
α
T
n
1
-
γ
T
·
cos
⁡
β
T
n
1
×
e
j
2
π
f
R
m
a
x
t
cos
⁡
α
R
n
2
-
γ
R
·
cos
⁡
β
R
n
2with(4)
α
T
L
O
S
≈
β
T
L
O
S
≈
β
R
L
O
S
≈
0
,
α
R
L
O
S
≈
π
,
τ
1
≈
D
c
,
τ
2
≈
ξ
n
i
+
R
T
/
R
c
,
τ
3
≈
ξ
n
1
,
n
2
+
R
T
+
R
R
c
,where S
B
1 and S
B
2 stand for the subcomponents of single-bounced rays from the Tx sphere and Rx sphere, respectively. DB means the double-bounced rays. η
S
B
i and η
D
B are the factors of normalized power, which specify the amount powers of single- and double-bounced rays contribute to the total normalized power. Herein, K is the Rician factor and c denotes the speed of light, and φ
n
i and φ
n
1
,
n
2 are independent and identically distributed (i.i.d.) random variables with uniform distribution over [
-
π
,
π
) [17]. In addition, f
T
m
a
x and f
R
m
a
x represent the maximum Doppler frequencies of the transmitter and the receiver, respectively. γ
T and γ
R stand for moving direction of the Tx and Rx, respectively. The statistical properties and performances of the proposed 3D V2V channel model will be analyzed as follows.Figure 1
The proposed 3D two-sphere channel model.


### 2.1. Envelope and Phase

Multiple uncorrelated fading processes bring the complex envelope, which can be written as(5)
h
t
=
h
i
t
+
j
h
q
t
.Considering the distribution of the channel envelope ξ
(
t
)
=
h
(
t
) and phase ϕ
(
t
)
=
a
r
g
h
(
t
), their corresponding probability density functions (PDF) can be written as [6](6)
p
ξ
z
=
z
δ
0
2
e
-
z
2
+
k
0
2
/
2
δ
0
2
·
I
0
z
k
0
δ
0
2
,
(7)
p
ϕ
ϕ
=
e
-
k
0
2
/
2
δ
0
2
2
π
·
1
+
k
0
δ
0
·
π
2
cos
⁡
ϕ
-
arg
⁡
h
L
O
S
t
,where z represents the amplitude variable, the symbol k
0 is defined as k
0
=
K
/
K
+
1, and I
0
(
·
) is the zeroth-order modified Bessel function of the first kind, while δ
0
2 stands for the mean power of the channel.


### 2.2. ACF and PSD

The normalized ACF between any two complex fading envelopes is defined as(8)
r
h
h
τ
=
E
h
t
·
h
∗
t
-
τ
E
h
t
2
·
E
h
t
-
τ
2
,where ·
∗ is the complex conjugate operation and E
· represents the statistical expectation operator. Then, we obtain the ACF of the proposed channel model. Taking the LoS and NLoS components into account, the equation of the ACF can be written as(9)
r
h
h
=
r
h
h
L
O
S
τ
+
∑
i
=
1
2
r
S
B
i
N
L
O
S
τ
+
r
D
B
N
L
O
S
τ
,where(10)
r
h
h
L
O
S
τ
=
K
K
+
1
e
j
2
π
/
λ
·
2
D
·
e
j
2
π
τ
f
T
m
a
x
cos
⁡
γ
T
-
f
R
m
a
x
cos
⁡
γ
R
,
r
S
B
i
N
L
O
S
τ
=
η
S
B
i
K
+
1
∬
-
π
π
e
j
2
π
τ
f
T
m
a
x
U
i
+
f
R
m
a
x
V
i
f
α
T
/
R
n
i
,
β
T
/
R
n
i
d
α
T
/
R
n
i
,
β
T
/
R
n
i
,
r
D
B
N
L
O
S
τ
=
η
D
B
K
+
1
⨌
-
π
π
⁡
e
j
2
π
τ
f
T
m
a
x
U
D
B
+
f
R
m
a
x
V
D
B
f
α
T
n
1
,
β
T
n
1
f
α
R
n
2
,
β
R
n
2
d
α
T
n
1
,
β
T
n
1
d
α
R
n
2
,
β
R
n
2
,where i means there are two subcomponents for single-bounced rays. For example, i
=
1 means the single-bounced rays from the Tx sphere. Similarly, i
=
2 means the single-bounced rays from the Rx sphere. And f
(
α
T
/
R
(
n
i
)
,
β
T
/
R
(
n
i
)
) stands for the joint azimuth and elevation angle distribution function of the n
ith scatterers on Tx or Rx. According to the mathematical law of triangle functions, the geometric relationships between each parameter are defined as(11)
U
1
=
cos
⁡
α
T
n
1
-
γ
T
·
cos
⁡
β
T
n
1
,
V
1
=
cos
⁡
α
R
n
1
-
γ
R
·
cos
⁡
β
R
n
1
,
α
R
n
1
≈
π
-
R
T
·
cos
⁡
β
T
n
1
D
·
sin
⁡
α
T
n
1
,
β
R
n
1
=
arccos
⁡
D
-
R
T
cos
⁡
β
T
n
1
·
cos
⁡
α
T
n
1
ξ
n
1
,
ξ
n
1
=
Q
n
1
2
+
R
T
2
sin
2
⁡
β
T
n
1
,
Q
n
1
=
D
+
R
·
cos
⁡
β
R
n
1
cos
⁡
α
R
n
1
,
U
2
=
cos
⁡
α
T
n
2
·
cos
⁡
β
T
n
2
,
V
2
=
cos
⁡
α
R
n
2
·
cos
⁡
β
R
n
2
,
(12)
α
T
n
2
≈
R
R
·
cos
⁡
β
R
n
2
D
·
sin
⁡
α
R
n
2
,
(13)
β
T
n
2
=
arccos
⁡
D
+
R
R
cos
⁡
β
R
n
2
·
cos
⁡
α
R
n
2
ξ
n
2
,
(14)
ξ
n
2
=
Q
n
2
2
+
R
R
2
sin
2
⁡
β
R
n
2
,
(15)
Q
n
2
≈
D
+
R
R
cos
⁡
β
R
n
2
cos
⁡
α
R
n
2
,
(16)
U
D
B
=
cos
⁡
α
T
n
1
·
cos
⁡
β
T
n
1
,
(17)
V
D
B
=
cos
⁡
α
R
n
2
·
cos
⁡
β
R
n
2
.The DPSD has been used for nonisotropic scattering channels and the theoretical expression DPSD is(18)
S
h
h
f
=
∫
-
∞
∞
r
h
h
τ
e
-
j
2
π
f
τ
d
τ
.


## 2.1. Envelope and Phase

Multiple uncorrelated fading processes bring the complex envelope, which can be written as(5)
h
t
=
h
i
t
+
j
h
q
t
.Considering the distribution of the channel envelope ξ
(
t
)
=
h
(
t
) and phase ϕ
(
t
)
=
a
r
g
h
(
t
), their corresponding probability density functions (PDF) can be written as [6](6)
p
ξ
z
=
z
δ
0
2
e
-
z
2
+
k
0
2
/
2
δ
0
2
·
I
0
z
k
0
δ
0
2
,
(7)
p
ϕ
ϕ
=
e
-
k
0
2
/
2
δ
0
2
2
π
·
1
+
k
0
δ
0
·
π
2
cos
⁡
ϕ
-
arg
⁡
h
L
O
S
t
,where z represents the amplitude variable, the symbol k
0 is defined as k
0
=
K
/
K
+
1, and I
0
(
·
) is the zeroth-order modified Bessel function of the first kind, while δ
0
2 stands for the mean power of the channel.


## 2.2. ACF and PSD

The normalized ACF between any two complex fading envelopes is defined as(8)
r
h
h
τ
=
E
h
t
·
h
∗
t
-
τ
E
h
t
2
·
E
h
t
-
τ
2
,where ·
∗ is the complex conjugate operation and E
· represents the statistical expectation operator. Then, we obtain the ACF of the proposed channel model. Taking the LoS and NLoS components into account, the equation of the ACF can be written as(9)
r
h
h
=
r
h
h
L
O
S
τ
+
∑
i
=
1
2
r
S
B
i
N
L
O
S
τ
+
r
D
B
N
L
O
S
τ
,where(10)
r
h
h
L
O
S
τ
=
K
K
+
1
e
j
2
π
/
λ
·
2
D
·
e
j
2
π
τ
f
T
m
a
x
cos
⁡
γ
T
-
f
R
m
a
x
cos
⁡
γ
R
,
r
S
B
i
N
L
O
S
τ
=
η
S
B
i
K
+
1
∬
-
π
π
e
j
2
π
τ
f
T
m
a
x
U
i
+
f
R
m
a
x
V
i
f
α
T
/
R
n
i
,
β
T
/
R
n
i
d
α
T
/
R
n
i
,
β
T
/
R
n
i
,
r
D
B
N
L
O
S
τ
=
η
D
B
K
+
1
⨌
-
π
π
⁡
e
j
2
π
τ
f
T
m
a
x
U
D
B
+
f
R
m
a
x
V
D
B
f
α
T
n
1
,
β
T
n
1
f
α
R
n
2
,
β
R
n
2
d
α
T
n
1
,
β
T
n
1
d
α
R
n
2
,
β
R
n
2
,where i means there are two subcomponents for single-bounced rays. For example, i
=
1 means the single-bounced rays from the Tx sphere. Similarly, i
=
2 means the single-bounced rays from the Rx sphere. And f
(
α
T
/
R
(
n
i
)
,
β
T
/
R
(
n
i
)
) stands for the joint azimuth and elevation angle distribution function of the n
ith scatterers on Tx or Rx. According to the mathematical law of triangle functions, the geometric relationships between each parameter are defined as(11)
U
1
=
cos
⁡
α
T
n
1
-
γ
T
·
cos
⁡
β
T
n
1
,
V
1
=
cos
⁡
α
R
n
1
-
γ
R
·
cos
⁡
β
R
n
1
,
α
R
n
1
≈
π
-
R
T
·
cos
⁡
β
T
n
1
D
·
sin
⁡
α
T
n
1
,
β
R
n
1
=
arccos
⁡
D
-
R
T
cos
⁡
β
T
n
1
·
cos
⁡
α
T
n
1
ξ
n
1
,
ξ
n
1
=
Q
n
1
2
+
R
T
2
sin
2
⁡
β
T
n
1
,
Q
n
1
=
D
+
R
·
cos
⁡
β
R
n
1
cos
⁡
α
R
n
1
,
U
2
=
cos
⁡
α
T
n
2
·
cos
⁡
β
T
n
2
,
V
2
=
cos
⁡
α
R
n
2
·
cos
⁡
β
R
n
2
,
(12)
α
T
n
2
≈
R
R
·
cos
⁡
β
R
n
2
D
·
sin
⁡
α
R
n
2
,
(13)
β
T
n
2
=
arccos
⁡
D
+
R
R
cos
⁡
β
R
n
2
·
cos
⁡
α
R
n
2
ξ
n
2
,
(14)
ξ
n
2
=
Q
n
2
2
+
R
R
2
sin
2
⁡
β
R
n
2
,
(15)
Q
n
2
≈
D
+
R
R
cos
⁡
β
R
n
2
cos
⁡
α
R
n
2
,
(16)
U
D
B
=
cos
⁡
α
T
n
1
·
cos
⁡
β
T
n
1
,
(17)
V
D
B
=
cos
⁡
α
R
n
2
·
cos
⁡
β
R
n
2
.The DPSD has been used for nonisotropic scattering channels and the theoretical expression DPSD is(18)
S
h
h
f
=
∫
-
∞
∞
r
h
h
τ
e
-
j
2
π
f
τ
d
τ
.


## 3. von Mises-Fisher Distribution

The VMF distribution [18, 19] is a close model for directional data distributed uniformly with rotational symmetry on the unit hypersphere S
p
-
1. When p is 3, the ordinary sphere usually corresponds to the set of all points embedded in the Euclidean space R
3. The vector Ω stands for any direction of scatterers on the unit sphere surface and it can be described in Cartesian coordinates as(19)
Ω
=
cos
⁡
β
cos
⁡
α
,
cos
⁡
β
sin
⁡
α
,
sin
⁡
β
T
,where [
·
]
T represents the transpose operation. It is obvious that β and α denote the coelevation and azimuth angles of vector Ω, respectively. Then, the general form of the VMF distribution can be written as(20)
f
p
Ω
;
μ
,
k
=
k
/
2
d
Γ
d
+
1
·
I
d
k
exp
⁡
k
μ
T
Ω
,where k controls the concentration of the distribution about the mean direction vector. In particular, when k
=
0, the distribution is isotropic, and when k
→
∞, the scattering becomes a point source on the surface. Furthermore, parameter d
=
p
/
2
-
1. From a physical point of view, the direction spread represents the degree of direction from the arrival subpaths of a cluster to the mean center of the AoA. Therefore, the mean center value of the AoA indicates the mean direction vector from the MS to the cluster and a subpath shows a multiple within a cluster of scatterers. For the ordinary sphere, the VMF PDF to characterize the distribution of effective scatterers can be defined as(21)
f
α
,
β
=
k
4
π
sinh
⁡
k
·
exp
⁡
k
cos
⁡
β
0
cos
⁡
β
cos
⁡
α
-
α
0
+
sin
⁡
β
0
sin
⁡
β
,where α
0 and β
0 represent the mean values of the azimuth and elevation angle, respectively, -
π
≤
α
0
≤
π, 0
≤
β
0
≤
π
.


## 4. TV Brownian Random Process

In this section, we will bring the temporal Brownian random process [17] to simulate a nonstationary TV 3D channel model. A standard BM process {
B
(
t
)
:
t
∈
[
0
,
T
]
} is a Wiener process, in which the increments satisfy a normal and independent distribution. Besides, it needs to satisfy the following conditions:(1)
B
(
0
)
=
0.
(2)
When0
≤
i
≤
j
≤
T, the increment B
(
j
)
-
B
(
i
) obeys the normal distribution with zero mean and variance j
-
i; that is, B
(
j
)
-
B
(
i
)
~
N
(
0
,
j
-
i
).
(3)
When0
≤
m
≤
n
≤
i
≤
j
≤
T, the increments B
(
m
)
-
B
(
n
) and B
(
j
)
-
B
(
i
) are statistically independent.Due to the fact that the same statistical properties impact the receiver between the AoA distribution motion and the MS random movement, this paper proposes that the variation of the AoA distribution can be used for simulating a 3D nonstationary channel. To model the motion process of the AoA distribution in the 3D plane, we provide a path model with BM movement components along the elevation and azimuth plane. The fluctuations of the motion path especially are modeled by two independent temporal BM processes:B
t
(
α
) and B
t
(
β
). Therefore, the motion path can be modeled as(22)
P
α
t
,
β
t
=
α
t
t
=
α
s
+
w
t
·
Δ
t
+
Δ
t
2
·
B
t
α
β
t
t
=
β
s
+
w
t
·
Δ
t
+
Δ
t
2
·
B
t
β
,where w stands for the random drift parameter used to control the behavior of the deterministic drift degrees along each axis and Δ
t represents the degree of randomness of the path. Meanwhile, Δ
t
=
1
/
w; α
s and β
s are the start points of the process along each axis. According to the proof of the BM paths model in [17], it is not difficult to find that the mean value of each angle always depends on time t and the ACF is also determined by a function of the time difference τ. Consequently, we can conclude that the BM path model is a nonstationary process. More details about Brownian random process can be found in [17].


## 5. The Nonstationary Properties of the Proposed Model

In this section, we consider the TV CIRs, ACF, and DPSD of the nonstationary process. At first, we introduce the BM process to the VMF scattering distribution, and the TV-VMF PDF can be rewritten as(23)
f
α
T
/
R
n
i
t
,
β
T
/
R
n
i
t
=
k
4
π
sinh
⁡
k
·
exp
⁡
k
·
cos
⁡
β
T
/
R
0
n
i
t
·
cos
⁡
β
T
/
R
n
i
cos
⁡
α
T
/
R
n
i
-
α
T
/
R
0
n
i
t
+
sin
⁡
β
T
/
R
0
n
i
t
·
sin
⁡
β
T
/
R
n
i
.Then the TV CIRs of the proposed model in (1)–(3) are obtained when the AoA and AoD become TV variables due to the random motion of the scatterers. Furthermore, the envelope PDF of the modified simulation model is the same as that in (6) and the phase PDF is also presented by (7) because both of them are unrelated with TV frequencies. We consider the TV ACF and DPSD of the nonstationary process [20] as(24)
r
h
~
h
~
t
,
τ
=
r
h
~
h
~
L
O
S
t
,
τ
+
∑
i
=
1
2
r
S
B
i
N
L
O
S
t
,
τ
+
r
D
B
N
L
O
S
t
,
τ
,where(25)
r
h
~
h
~
L
O
S
t
,
τ
=
K
K
+
1
e
j
2
π
/
λ
·
2
D
·
e
j
2
π
τ
f
T
m
a
x
cos
⁡
γ
T
-
f
R
m
a
x
cos
⁡
γ
R
,
r
S
B
i
N
L
O
S
t
,
τ
=
η
S
B
i
K
+
1
∬
-
π
π
e
j
2
π
τ
f
T
m
a
x
U
i
+
f
R
m
a
x
V
i
f
α
T
/
R
n
i
t
,
β
T
/
R
n
i
t
d
α
T
/
R
n
i
t
,
β
T
/
R
n
i
t
,
r
D
B
N
L
O
S
t
,
τ
=
η
D
B
K
+
1
⨌
-
π
π
⁡
e
j
2
π
τ
f
T
m
a
x
U
D
B
+
f
R
m
a
x
V
D
B
f
α
T
n
1
t
,
β
T
n
1
t
f
α
R
n
2
t
,
β
R
n
2
t
d
α
T
n
1
t
,
β
T
n
1
t
d
α
R
n
2
t
,
β
R
n
2
t
.In particular, the parameters in (24)-(25) are given in Section 2. Note that the AoAs (
α
R
(
t
)
,
β
R
(
t
)
) and AoDs (
α
T
(
t
)
,
β
T
(
t
)
) are applied together with the TV variables [18].It is noteworthy that the definition of the TV ACF keeps the symmetrical delay characteristicτ, r
h
~
h
~
(
t
,
τ
)
=
r
h
~
h
~
(
t
-
τ
), and r
h
~
h
~
(
t
,
τ
) is the real function if h
~
(
t
) is a real process, which results in the Fourier transform having the symmetrical and real characteristics. Beyond the TV DPSD, the short-time DPSD is widely used to analyze the nonstationary process [17], which is the amplitude’s square of the signal’s STFT [21, 22]:(26)
S
T
h
~
h
~
f
,
t
=
∫
-
∞
∞
r
h
~
h
~
t
,
τ
e
-
j
2
π
f
τ
·
w
t
-
τ
d
τ
,where (
t
-
τ
) is the analytic window slipping with time. The nonstationary process can be viewed as stationary within the analysis window.


## 6. Simulation Results and Numerical Analysis

In this section, simulations are carried out to illustrate the nonstationary properties of our proposed 3D TV channel model based on the VMF scatterer distribution. The impact of the model parameters on the VMF scatterer distribution is investigated first. Then, the nonstationary properties of the proposed two-sphere channel model are evaluated and analyzed in terms of the AoAs and AoDs motion path at the Tx and Rx, the TV ACF of the complex channel gain, and the TV DPSD.Firstly, we consider the performance of the VMF scatterers distribution. From Figure2, it is not difficult to find that the parameter k controls the shape of the distribution. Meanwhile, Figure 4 shows the VMF PDF in 3D coordinates by setting the mean angles α
0
=
9
0
° and β
0
=
4
5
°. It is clear that the two-dimensional (2D) von Mises (VM) PDF is derived from the 3D VMF PDF for azimuth angle α with β
=
0
°. Furthermore, the scatterer distribution based on the VMF PDF has been shown in Figure 3. Assuming scatterer distribution is dependent on the axis of symmetry specified by the z-axis, it can be seen that the larger the value of k, the higher the density near the z-axis. Especially when k
=
0, the distribution is isotropic.Figure 2
2D VM PDF.Figure 3
Effect of the concentration parameter.Figure 4
3D VMF PDF.Secondly, Figure5 illustrates a realization of the temporal BM motion path where k
=
1, w
=
10, and Δ
t
=
0.01. As seen from this figure, the trajectories of α
0
(
i
)
(
t
) and β
0
(
i
)
(
t
) are varying with time, respectively. The properties of the BM process prove that the proposed path results in a nonstationary nonisotropic channel model.Figure 5
The simulation of Brownian paths.Thirdly, Figures6–8 show the proposed TV 3D models in detail for each statistical value. The simulation settings are as follows: f
c
=
5.9 GHz, D
=
300 m, f
T
m
a
x
=
f
R
m
a
x
=
570 Hz, R
=
25 m, γ
T
=
γ
R
=
0
°, α
R
0
(
n
1
)
=
147
.
8
°, β
R
0
(
n
1
)
=
17
.
2
°, α
T
0
(
n
2
)
=
31
.
6
°, and β
T
0
(
n
2
)
=
171
.
6
°. Considering the vehicular traffic density (VTD) based on [4], we have K
=
3.876, k
=
3.6, η
S
B
1
=
0.625, η
S
B
2
=
0.225, and η
D
B
=
0.15 for low and K
=
0.156, k
=
0.6, η
S
B
1
=
0.15, η
S
B
2
=
0.225, and η
D
B
=
0.625 for high VTD scenarios because both parameters K and k are related to the distribution of scatterers (normally, larger values of k corresponding to less dense moving vehicles.).Figure 6
The absolute values of the temporal ACF for the proposed 3D TV channel model.Figure 7
The PDF of amplitude.Figure 8
The PDF of phase.Figure6 validates the absolute values of the temporal ACF for the proposed 3D channel model in both low and high VTD scenarios. It can be observed that the temporal ACF is affected by the VTD; namely, the temporal ACF in low VTD scenarios is always higher than that in high VTD scenario. Additionally, Figures 7 and 8 depict the envelope and the phase PDF of 3D channel model determined by the parameters K and k, respectively.Finally, we assume that the mean elevation and azimuth angles are moving along the TV-BM motion paths shown in Figure5 and we set the maximum Doppler frequency of the transmitter and receiver as 570 Hz. Meanwhile, the effect of the proposed TV-BM motion path on the scattering model is the same as the one in which the vehicle moves randomly in different directions on the AoAs and AoDs. Therefore, as shown in Figures 9 and 12, the absolute value of the resulting local ACF ((24) and (25)) is illustrated. According to these figures, it is manifest that the shapes of the ACF change with different values of time which is due to the nonstationarity of the TV AoAs and AoDs motion processes. It has been proven that the AoAs and AoDs motion processes based on the TV-BM path model are first-order nonstationary in [17]. For the same reason, the TV PSD, which is derived from the local ACF by utilizing the FFT with respect to time delay τ, is illustrated in Figures 10 and 13 for low and high VTD scenarios, respectively. Both of them show that the similar U-shape of the PSD varies with time t. And this is another way to display the nonstationary properties of the proposed channel model. In addition, the short-time (ST) DPSD S
T
h
~
h
~
(
f
,
t
) (see (26)) is shown in Figures 11 and 14. Compared with the FFT DPSD and ST DPSD, we can conclude that the nonstationary properties of them are very similar to each other. This means the STFT has the ability to process the nonstationary channel model. In addition, by comparing the nonstationary properties of the channel model in the low and high VTD scenarios, it can be found that the DPSD derived from the low VTD scenarios is closer to the ideal Jakes DPSD (U-shape spectrum). That is because the effect of nonstationary properties of channel model is more obvious in the high VTD assumption scenarios.Figure 9
Time-variant ACF.Figure 10
Time-variant DPSD (FFT for high VTD).Figure 11
Time-variant DPSD (STFT for high VTD).Figure 12
Time-variant ACF (low VTD).Figure 13
Time-variant DPSD (FFT for low VTD).Figure 14
Time-variant DPSD (STFT for low VTD).


## 7. Conclusion

In this paper, a novel 3D nonstationary channel model for V2V is proposed. By employing the proposed temporal BM process with the VMF PDF scatterer distribution, we derive a nonstationary nonisotropic channel model which can be applied to simulate the 3D channels in real time. Moreover, the dynamic changes of the local ACF and PSD have been provided. Lastly, the nonstationary properties of the proposed channel model are verified by simulation. In particular, we have proven that the STFT is also valid for analyzing nonstationary channel models.



---

*Source: 2161460-2016-05-18.xml*