# Analytical Nonstationary 3D MIMO Channel Model for Vehicle-to-Vehicle Communication on Slope

Authors: Kaizhen Liu; Zaixue Wei; Sibo Chen
Journal: International Journal of Antennas and Propagation
 (2021)

---

## Abstract

Vehicle-to-vehicle communication plays a strong role in modern wireless communication systems, appropriate channel models are of great importance in future research, and propagation environment with slope is one special kind. In this study, a novel three-dimensional nonstationary multiple-input multiple-output channel model for the sub-6 GHz band is proposed. This model is a regular-shaped multicluster geometry-based analytical model, and it combines the line-of-sight component and multicluster scattering rays as the nonline-of-sight components. Each cluster of scatterers represents the influence of different moving vehicles on or near a slope, and scatterers are, respectively, distributed within two spheres around the transmitter and the receiver. In this model, it is considered that the azimuth and elevation angles of departure and arrival are jointly distributed and conform to the von Mises–Fisher distribution, which can easily control the range and concentration of the scatterers within spheres to mimic the real-world situation well. Moreover, the impulse response and the autocorrelation function of the corresponding channel is derived and proposed; then, the Doppler power spectrum density of the channel is simulated and analyzed. In addition, the nonstationary characteristics of the presented channel model are observed through simulations. Finally, the simulation results are compared with measurement data in order to validate the utility of the proposed model.


---

## Body

## 1. Introduction

In recent years, the applications of vehicle-to-vehicle (V2V) communications have become more extensive and valuable, such as vehicular networking [1–3], the cooperative vehicular system [4], and intelligent transportation [5]. Besides, the accuracy of performance evaluation for V2V communications is to some extent determined by precise modeling of the channel. Reference [6] analyzed vehicular communications from the physical layer. With the mature multiple-input multiple-output (MIMO) technology [7, 8], the performance of V2V communications has been greatly improved. Compared to other communication scenarios, V2V communications have their unique characteristics, i.e., both the transmitter (Tx) and the receiver (Rx) are in high-speed movements, the shapes and heights of different vehicles are similar, and the heights of vehicles are generally between one and two meters, so vehicles usually have antennas with low elevation angles. Thus, the obstruction of the slope will affect the performance of V2V communication seriously. With these special features, the channel modeling of V2V communications is essential and necessary, especially the propagation model on slope terrain. The path loss for V2V scenario with slope has been studied in several works. References [9, 10] proposed path loss prediction models for slope terrain in small-sized cells. However, research on the modeling of scatterers still needs to be finished.The modeling of V2V channels can be of great help in testing the performance and properties of V2V channels. Reference [11] presented propagation models for urban environment. Besides, several precise channel models for other scenarios have been proposed by other researchers. According to this literature, channel models can generally be classified into deterministic models and stochastic models. Stochastic models are further categorized as geometry-based stochastic models (GBSMs) [12] and nongeometry-based stochastic models (NGBSMs) [13]. GBSMs assume that scatterers are distributed on a certain geometry and use methods similar to ray-tracing to represent scattering rays. In further classification, GBSMs can be categorized as regular-shaped GBSMs (RS-GBSMs) [14] and irregular-shaped GBSMs (IS-GBSMs) [15]. RS-GBSMs assume that these geometries are regular shapes, e.g., sphere, ring, or cylinder. In other ways, GBSMs can be divided into two-dimensional (2D) GBSMs [16] and three-dimensional (3D) GBSMs [17] according to the dimension of models. Considering the pattern of antenna arrays, the models can be categorized as single-input single-output (SISO) models [18] and MIMO models [19]. Additionally, in [20], a 3D model adopting massive-MIMO technology was proposed.Comparisons with measurements [21–23] showed that GBSMs could be applied to V2V channel modeling suitably. V2V channels were modeled in 2D due to the very narrow elevation angle of the vehicle antennas at the early time. In [24], the authors proposed a nonstationary 2D RS-GBSM consisting of multiple confocal ellipses for high-mobility intelligent transportation systems. In [25], a 2D RS-GBSM with two-ring for a nonisotropic scattering environment was proposed. Furthermore, 2D models were extended to 3D models because they could better simulate the realistic propagation environment, e.g., the 2D two-circular RS-GBSM [25] was updated to the 3D two-cylinder RS-GBSM in [26], and the 2D two-circular and one-ellipse RS-GBSM [27] was extended to the 3D two-sphere and one-ellipse-cylinder RS-GBSM in [28]. Including the above listed, 3D GBSMs can be classified according to the geometries of scattering, such as spheres, cylinders, one/two semiellipsoid [29–31], and one/two hollow semiellipsoid around the Tx/Rx [32, 33]. In addition, GBSMs can be divided into stationary GBSMs and nonstationary GBSMs [34–42]. The nonstationary GBSMs are closer to reality due to the time-varying parameters of the V2V channel. According to the abovementioned literature, it can be noticed that in the RS-GBSM based channel modeling, the modeling of the scatterer distribution is crucial. In RS-GBSM-based models, isotropic and nonisotropic distributions of scatterers are considered, and the latter is more reasonable. In [27, 35, 36, 43], von Mises distributions were adopted to model the distributions of the scatterers in 2D modeling as nonisotropic models. In [28, 37, 44–48], 3D nonisotropic distributions of scatterers were modeled with von Mises–Fisher (VMF) distributions. By now, the VMF distribution is perhaps the most suitable distribution for 3D nonisotropic scattering modeling.In the aforementioned literature, although nonisotropic 3D scatterers have been modeled, scatterers were always assumed to be concentrated in the direction of the line-of-sight (LoS) path. Further discussions on multidirection scatterers modeling should be conducted. Besides, limited specific models are for vehicles near or on a slope as well as overtaking cases.The motivation of this research is to present a flexible 3D nonstationary V2V channel model with multidirection distributing scatterers on or near slope terrain. Based on the two-sphere 3D RS-GBSM V2V nonstationary models in [28, 37, 48, 49], we abandoned the additional elliptical cylinder and introduced multidirection scattering clusters to model the multidirectional isotropic scatterers. More distinctly, in the model we proposed, rays are scattered by other vehicles in several directions in the form of cluster around Tx and Rx. This is the major difference between our model and [28, 49, 50] and basic theories found in [51]. For the modeling of nonisotropic scattering, the elevation and azimuth angles of arrival (AoAs) and the elevation and azimuth angles of departure (AoDs) are jointly distributed and subjected to VMF distribution, as in [45, 47]. Finally, we compare the simulation results with the measurement data.In summary, the main contributions of this study are as follows:(1)
A nonstationary 3D V2V MIMO RS-GBSM is proposed, which is applied to the scenario where Tx and Rx are surrounded by multiple vehicles on slope, and the vehicle can change lanes or overtake during the movement. The statistical properties of this MIMO channel model are derived and analyzed.(2)
The nonstationary characteristics of the proposed model are derived and simulated. The simulation results support the utility of the model.(3)
The proposed model is compared with the measurement data in four different communication scenarios. The comparisons include root mean square (RMS) delay spread, RMS Doppler spread, and root mean squared error (RMSE).The rest of this study is arranged as follows. Section2 gives detailed definition of the model. The expression of channel impulse response (CIR) is given in Section 3. In Section 4, the statistical properties are calculated. Section 5 is the simulation results and the comparison with measurement data. The conclusion is presented in Section 6.Notation: the superscript⋅T denotes the transpose operator. The superscript ⋅∗ denotes the conjugate operator. x is the norm of a vector x. E⋅ denotes the expectation operator. pα,β is the probability density function (PDF) of the VMF distribution with parameters α and β.


## 2. Nonstationary 3D RS-GBSM for V2V MIMO Channels

The presented 3D RS-GBSM is for MIMO V2V communication systems, and the system structure is shown in Figures1 and 2. We assume that both of Tx and Rx are surrounded by several vehicles. For convenience on building the coordinates, we assume that Tx follows Rx, and it can be very easily adapted to other cases. The coordinates are as shown in Figure 2. According to the different positions of Tx and Rx on the slope, it can be divided into 2 scenarios: (1) all vehicles are on the slope; (2) Tx and surrounding vehicles are off the slope and moving towards the bottom end of the slope while Rx and surrounding vehicles are on the slope.Figure 1
The proposed 3D RS-GBSM for V2V MIMO channels att = 0 s.Figure 2
The real scenario of the 3D RS-GBSM. (a) Scenario 1 and (b) scenario 2.
(a)(b)


### 2.1. Model Specific Description

The distribution of scatterers in the multicluster form at the Tx side is shown in Figure3. Since the distribution of scatterers at the Rx side is the same as the Tx side, we only illustrate the Tx side. Taking cluster S1 as an example, the scatterers are distributed around the main path L1 and obeying the VMF distribution. Since vehicles have a certain volume, the radii of different scatterers to Tx and Rx are different. Limited by the radii, they are distributed between s11 and s12, i.e., all scatterers are distributed within the outermost and innermost two spheres. For clearer description of the structure, the concerned paths and parameters are shown in Figure 1. Major parameters in this study are listed in Table 1. According to the realistic situation, both the LoS and non-LoS (NLoS) components are included. A two-sphere model is developed for the NLoS components. As shown in Figures 1 and 3, effective scatterers are distributed between two spheres, and they are distributed in the form of clusters. Taking the four-cluster scenario as an example, there are four clusters of scatterers at both ends of Tx and Rx, representing, respectively, four vehicles around Tx and Rx, and scatterers are distributed on the surface and in these vehicles. As mentioned above, each effective scatterer has a different radius, and the specific way of modeling is shown in Figure 3. In Figure 1, we assume that both of Tx and Rx are equipped with multiple antennas. The antenna elements are omnidirectional, and the numbers of antennas at Tx and Rx are denoted as nt and nr, respectively. In this model, taking 2 × 2 MIMO as an example, the two antennas of the Tx are labeled p1 and p2, and the two antennas of the Rx are labeled l1 and l2. In addition, we consider that both Tx and Rx are moving at high speeds, and the directions of the movements of the Tx and Rx are jointly defined by elevation angles and azimuth angles. The azimuth and elevation angles of the movement directions are recorded as γTt,γRt,ϕTt,andϕRt, and the velocities of the movements are recorded with vTtandvRt. As mentioned before, the effective scatterers are distributed over a section of the sphere and are distributed in the form of clusters, denoted by Sinij (i=1,…,Nc and j = 1, 2), where i represents that the scatterers are located in the ith cluster, j = 1 represents that the scatterers are on the Tx sphere, and j = 2 represents that the scatterers are on the Rx sphere. In the ith cluster, nij represents that the scatterer is the nijth scatterer of this cluster. αTnijt and βTnijt are the azimuth AoD (AAoD) and elevation AoD (EAoD) of the path. αRnijt and βRnijt are the azimuth AoA (AAoA) and elevation AoA (EAoA) of the path. The scattering paths can be divided into three types: (1) rays that are only reflected once at the Tx end, and these paths are called LSB1. (2) Rays that are only reflected once at the Rx end, and these paths are called LSB2. (3) Rays that are reflected twice at both of the Tx and Rx ends, and these paths are called LDB. With the distances between adjacent antennas δTandδR and the radii of the two spheres RTandRR, we have δT,δR≪RT,RR.Figure 3
The detailed distribution of scatterers on the Tx sphere. The distribution of scatterers near Rx is consistent with Tx.Table 1
Definition and unit of key parameters in model.
ParametersDefinitionsUnitsfcCentral frequencyHzDtDistance between Tx and Rx in scenario 1mD1tDistance between Tx and the bottom end of slope in scenario 2mD2tDistance between the bottom end of slope and Rx in scenario 2mRTt,RRtThe radius of Tx and Rx spheresmδT,δRThe antenna elements spacing at Tx and RxmvTt,vRtVelocity vectors of Tx and Rxm/svscatterertVelocity vectors of vehicles around Tx and Rxm/sϕTt,ϕRtThe elevation angles of movement direction of Tx and RxradγTt,γRtThe azimuth angles of movement direction of Tx and RxradfTmaxt,fRmaxtThe maximum Doppler shifts caused by Tx and RxHzfTDopplert,fRDopplertThe Doppler shifts caused by Tx and RxHzθTt,θRtThe horizontal orientation of antenna array of Tx and RxradαTnijtThe AAoD atSinijradαRnijtThe AAoA fromSinijradβTnijtThe EAoD atSinijradβRnijtThe EAoA fromSinijradαTLoSt,βTLoStThe AAoD and EAoD of the LoS pathradαRLoSt,βRLoStThe AAoA and EAoA of the LoS pathraddS1,S2tThe distance betweenS1 and S2mNcThe number of vehicles (clusters) around Tx and RxNi1,Ni2The number of scatterers in theith cluster on the Tx/Rx sphereθslopeThe tilt angle of the sloperad


### 2.2. Vehicle Trajectory

We have considered not only the situation where vehicles move in a straight line but also the situation where one vehicle changes lane or overtakes. Vehicle’s trajectory is shown in Figure4. The initial time and final time are recorded as tini and tfin. The location, velocity, and acceleration of the vehicle in the X-axis and Y-axis directions at the initial time and the final time are recorded as (xini, x˙ini, x¨ini, yini, y˙ini,y¨ini) and xfin,x˙fin,x¨fin,yfin,y˙fin,y¨fin; x, x˙, and x¨ denote distance, velocity, and acceleration, respectively.Figure 4
Vehicle lane-changing trajectory.We define the location of the vehicle in theX-axis and the Y-axis direction as(1a)xt=a5t5+a4t4+a3t3+a2t2+a1t+a0,(1b)yt=b5t5+b4t4+b3t3+b2t2+b1t+b0.Then, we letA=a5a4a3a2a1a0, B=b5b4b3b2b1b0, and we define the time parameter matrix:(2)T6×6=tini5tini4tini3tini2t15tini44tini33tini22tini1020tini312tini26tini200tfin5tfin4tfin3tfin2tfin15tfin44tfin33tfin2tfin1020tfin312tfin26tfin200.Now, we obtain(3a)xinix˙inix¨inixfinx˙finx¨finT=T6×6⋅AT,(3b)yiniy˙iniy¨iniyfiny˙finy¨finT=T6×6⋅BT.Finally, the solution of the inhomogeneous linear equations (3a) and (3b) can be derived through the boundary conditions. Using this method, we can quickly obtain the track of vehicle’s overtaking or lane-changing. However, the trajectory does not consider the dynamic characteristics of the vehicle, and the generated acceleration should be limited by the condition of vehicle engine, tire, ground friction, and other restrictions. Commonly, we have −2.5m/s2<x¨<2.5m/s2 and 2m/s2<y¨<2m/s2.


## 2.1. Model Specific Description

The distribution of scatterers in the multicluster form at the Tx side is shown in Figure3. Since the distribution of scatterers at the Rx side is the same as the Tx side, we only illustrate the Tx side. Taking cluster S1 as an example, the scatterers are distributed around the main path L1 and obeying the VMF distribution. Since vehicles have a certain volume, the radii of different scatterers to Tx and Rx are different. Limited by the radii, they are distributed between s11 and s12, i.e., all scatterers are distributed within the outermost and innermost two spheres. For clearer description of the structure, the concerned paths and parameters are shown in Figure 1. Major parameters in this study are listed in Table 1. According to the realistic situation, both the LoS and non-LoS (NLoS) components are included. A two-sphere model is developed for the NLoS components. As shown in Figures 1 and 3, effective scatterers are distributed between two spheres, and they are distributed in the form of clusters. Taking the four-cluster scenario as an example, there are four clusters of scatterers at both ends of Tx and Rx, representing, respectively, four vehicles around Tx and Rx, and scatterers are distributed on the surface and in these vehicles. As mentioned above, each effective scatterer has a different radius, and the specific way of modeling is shown in Figure 3. In Figure 1, we assume that both of Tx and Rx are equipped with multiple antennas. The antenna elements are omnidirectional, and the numbers of antennas at Tx and Rx are denoted as nt and nr, respectively. In this model, taking 2 × 2 MIMO as an example, the two antennas of the Tx are labeled p1 and p2, and the two antennas of the Rx are labeled l1 and l2. In addition, we consider that both Tx and Rx are moving at high speeds, and the directions of the movements of the Tx and Rx are jointly defined by elevation angles and azimuth angles. The azimuth and elevation angles of the movement directions are recorded as γTt,γRt,ϕTt,andϕRt, and the velocities of the movements are recorded with vTtandvRt. As mentioned before, the effective scatterers are distributed over a section of the sphere and are distributed in the form of clusters, denoted by Sinij (i=1,…,Nc and j = 1, 2), where i represents that the scatterers are located in the ith cluster, j = 1 represents that the scatterers are on the Tx sphere, and j = 2 represents that the scatterers are on the Rx sphere. In the ith cluster, nij represents that the scatterer is the nijth scatterer of this cluster. αTnijt and βTnijt are the azimuth AoD (AAoD) and elevation AoD (EAoD) of the path. αRnijt and βRnijt are the azimuth AoA (AAoA) and elevation AoA (EAoA) of the path. The scattering paths can be divided into three types: (1) rays that are only reflected once at the Tx end, and these paths are called LSB1. (2) Rays that are only reflected once at the Rx end, and these paths are called LSB2. (3) Rays that are reflected twice at both of the Tx and Rx ends, and these paths are called LDB. With the distances between adjacent antennas δTandδR and the radii of the two spheres RTandRR, we have δT,δR≪RT,RR.Figure 3
The detailed distribution of scatterers on the Tx sphere. The distribution of scatterers near Rx is consistent with Tx.Table 1
Definition and unit of key parameters in model.
ParametersDefinitionsUnitsfcCentral frequencyHzDtDistance between Tx and Rx in scenario 1mD1tDistance between Tx and the bottom end of slope in scenario 2mD2tDistance between the bottom end of slope and Rx in scenario 2mRTt,RRtThe radius of Tx and Rx spheresmδT,δRThe antenna elements spacing at Tx and RxmvTt,vRtVelocity vectors of Tx and Rxm/svscatterertVelocity vectors of vehicles around Tx and Rxm/sϕTt,ϕRtThe elevation angles of movement direction of Tx and RxradγTt,γRtThe azimuth angles of movement direction of Tx and RxradfTmaxt,fRmaxtThe maximum Doppler shifts caused by Tx and RxHzfTDopplert,fRDopplertThe Doppler shifts caused by Tx and RxHzθTt,θRtThe horizontal orientation of antenna array of Tx and RxradαTnijtThe AAoD atSinijradαRnijtThe AAoA fromSinijradβTnijtThe EAoD atSinijradβRnijtThe EAoA fromSinijradαTLoSt,βTLoStThe AAoD and EAoD of the LoS pathradαRLoSt,βRLoStThe AAoA and EAoA of the LoS pathraddS1,S2tThe distance betweenS1 and S2mNcThe number of vehicles (clusters) around Tx and RxNi1,Ni2The number of scatterers in theith cluster on the Tx/Rx sphereθslopeThe tilt angle of the sloperad


## 2.2. Vehicle Trajectory

We have considered not only the situation where vehicles move in a straight line but also the situation where one vehicle changes lane or overtakes. Vehicle’s trajectory is shown in Figure4. The initial time and final time are recorded as tini and tfin. The location, velocity, and acceleration of the vehicle in the X-axis and Y-axis directions at the initial time and the final time are recorded as (xini, x˙ini, x¨ini, yini, y˙ini,y¨ini) and xfin,x˙fin,x¨fin,yfin,y˙fin,y¨fin; x, x˙, and x¨ denote distance, velocity, and acceleration, respectively.Figure 4
Vehicle lane-changing trajectory.We define the location of the vehicle in theX-axis and the Y-axis direction as(1a)xt=a5t5+a4t4+a3t3+a2t2+a1t+a0,(1b)yt=b5t5+b4t4+b3t3+b2t2+b1t+b0.Then, we letA=a5a4a3a2a1a0, B=b5b4b3b2b1b0, and we define the time parameter matrix:(2)T6×6=tini5tini4tini3tini2t15tini44tini33tini22tini1020tini312tini26tini200tfin5tfin4tfin3tfin2tfin15tfin44tfin33tfin2tfin1020tfin312tfin26tfin200.Now, we obtain(3a)xinix˙inix¨inixfinx˙finx¨finT=T6×6⋅AT,(3b)yiniy˙iniy¨iniyfiny˙finy¨finT=T6×6⋅BT.Finally, the solution of the inhomogeneous linear equations (3a) and (3b) can be derived through the boundary conditions. Using this method, we can quickly obtain the track of vehicle’s overtaking or lane-changing. However, the trajectory does not consider the dynamic characteristics of the vehicle, and the generated acceleration should be limited by the condition of vehicle engine, tire, ground friction, and other restrictions. Commonly, we have −2.5m/s2<x¨<2.5m/s2 and 2m/s2<y¨<2m/s2.


## 3. Channel Impulse Response

In this section, the detailed deduction of the CIR is given. Take pathp1⟶l1 as an example. Since scenario 1 and scenario 2 are very similar in derivation, we only give the derivation process for scenario 1. The CIR can be divided into three parts: LoS, single-bounced (SB), and double-bounced (DB) components. Then, the CIR can be expressed as(4)hp1l1t=hLoSt+∑i=12hp1l1SBit+hp1l1DBt.The Doppler shift caused by the Tx’s motion is(5a)fTDopplert=fTmaxtcosβTtcosαTtcosϕTtcosγTt+cosβTtsinαTtcosϕTtsinγTt+sinβTtsinϕTt.Correspondingly, the Doppler shift caused by the Rx’s motion can be expressed as(5b)fRDopplert=fRmaxtcosβRtcosαRtcosϕRtcosγRt+cosβRtsinαRtcosϕRtsinγRt+sinβRtsinϕRt,where fTmaxt=vTt/λ and fRmaxt=vRt/λ are the maximum Doppler shifts at time t, λ is the carrier wavelength, and vTt and vRt are the velocity vectors of the Tx and Rx.Then, the three parts of the CIR can be rewritten as(6a)hp1l1LoSt=KK+1e−j2π/λdp1,l1t×ej2πfTDopplerLoSttej2πfRDopplerLoStt,(6b)hp1l1SBit=1K+1εSBilimNi⟶∞∑k=1Nc∑nki=1Nki1Ni×ejψnki−2π/λdp1,Sknkit+dSknki,l1t×ej2πfTDopplerSBittej2πfRDopplerSBitt,(6c)hp1l1DBt=1K+1εDBlimN1,N2⟶∞∑i,k=1Nc∑ni1,nk2=1Ni1,Nk21N1N2×ejψni1,nk2−2π/λdDBt×ej2πfTDopplerDBttej2πfRDopplerDBtt,where Ni1 is the number of effective scatterers in the ith cluster on the Tx sphere, Ni2 is the number of effective scatterers in the ith cluster on the Rx sphere, N1=∑i=1NcNi1, N2=∑i=1NcNi2, dDBt=dp1,Sini1t+dSini1,Sknk2t+dSknk2,l1t, j2=−1, K is the Ricean factor, and SB1, SB2, and DB, respectively, represent LSB1, LSB2, and LDB which have been introduced in Section 2. Parameters εSBiandεDB are related to the power of different paths, they control the proportions of power among paths LSB1, LSB2, and LDB, and they satisfy εSB12+εSB22+εDB2=1. ψnki and ψni1,nk2 are the random phases, and they are uniformly distributed within −π,π. fT/RDopplerLoS/SBi/DBt, respectively, represent the Doppler shifts at Tx and Rx of the LoS, SB, and DB components. The expressions are derived as (5a) and (5b). For LoS component, αT/Rt=αT/RLoSt and βT/Rt=βT/RLoSt. In terms of SB and DB components, we have αT/Rt=αT/Rnijt and βT/Rt=βT/Rnijt.We consider that the locations of Tx and Rx are known through the global positioning system; the location ofp1 is xp1t,yp1t,zp1t, and the location of l1 is xl1t,yl1t,zl1t. In scenario 1, Tx, Rx, and surrounding vehicles are constantly moving on the slope. At initial time, xp1t=0, yp1t=δTcosθTt/2, zp1t=δTsinθTt/2, xl1=Dtcosθslope, yl1t=δRcosθRt/2, and zl1t=δRsinθRtsinθslope/2. The location of Sini1 can be denoted as [xs1t, ys1t, zs1t]:(7a)xs1t=RTtcosβTni1tcosαTni1t,(7b)ys1t=RTtcosβTni1tsinαTni1t,(7c)zs1t=RTtsinβTni1t.The location ofSini2 is denoted as [xs2t, ys2t, zs2t] and(8a)xs2t=RRtcosβRni2tcosαRni2t+Dtcosθslope,(8b)ys2t=RRtcosβRni2tsinαRni2t,(8c)zs2t≈Dtcosθslope+RRtcosαRni2t×tanθslope+RRtsinβRni2t.Then, we can easily calculatedp1,l1t, dp1,Sini1t, dSjnj2,l1t, and dSini1,Sjnj2t; take dp1,l1t, for example, we have(9)dp1,l1t=Δxl1p12t+Δyl1p12t+Δzl1p12t=Dtcosθslope2+δRcosθRt−δTcosθTt22+δRsinθRtsinθslope−δTsinθTt22.InDB paths, EAoD/AAoD and EAoA/AAoA are independent, whereas in SB1 and SB2 paths, they are related. The generation of angles is introduced in Section 4.


## 4. Space-Time Correlation Function and Doppler Power Spectrum Density

### 4.1. Space-Time Correlation Function

The normalized space-time correlation function (STCF) between pathp1⟶l1 and p2⟶l2 is defined as(10)Rp1l1,p2l2δT,δR,t,τ=Ehp1l1thp2l2∗t+τEhp1l1t2Ehp2l2t2=K+1Ehp1l1thp2l2∗t+τ.The STCF can be written into(11)Rp1l1,p2l2δT,δR,t,τ=Rp1l1LoS,p2l2LoSδT,δR,t,τ+∑i=12Rp1l1SBi,p2l2SBiδT,δR,t,τ+Rp1l1DB,p2l2DBδT,δR,t,τ.(a)
STCF of LoS component: we also only give the derivation for scenario 1 because it is similar to the scenario 2. The STCF of the LoS component is(12)Rl1p1LoS,l2p2LoSδT,δR,t,τ=Ke−j2π/λΔdt×ej2πτfTDopplerLoSt+fRDopplerLoSt=Ke−j2π/λdp1,l1t+dp2,l2t×ej2πτfTDopplerLoSt+fRDopplerLoSt,wheredp1,l1t and dp2,l2t can be calculated as (9). fTDopplerLoSt and fRDopplerLoSt have been expressed in Section 3.(b)
STCF of SB components: we consider that the EAoD and the AAoD inSB1 and EAoA and AAoA in SB2 are contributed to VMF contribution; thus, the general PDF of these angles is(13)pα,β=kcosβ04πsinhk×ekcosβ0cosβcosα−α0+sinβ0sinβ,whereα,β∈−π,π, and α0,β0∈−π,π. Since there are Nc clusters of scatterers, there are also Nc groups of α0 and β0 corresponding to piα,β (i = 1, …, Nc). Next, if the number of the scatterers tends to infinity, the discrete random variables αT/Rnij and βT/Rnij can be replaced by continuous random variables αT/R and βT/R. Then, the STCF of SBi component becomes(14)Rp1l1SB1,p2l2SB1δT,δR,t,τ=εSBi2∑k=1Nc∫−ππ∫−ππe−j2π/λΔdt×ej2πτfTDopplerSBit×ej2πτfRDopplerSBit×pkαT,βTdαTdβT,whereΔdt=dp1,Skt+dSk,l1t−dp2,Skt−dSk,l2t, and the calculations are similar to (9).(c)
STCF of DB components:(15)Rp1l1DB,p2l2DBδT,δR,t,τ=εDB2∑i=1Nc∑k=1Nc∫−ππ∫−ππ∫−ππ∫−ππe−j2π/λΔdej2πτfTDopplerDBtej2πτfTDopplerDBt×piαT,βTpkαR,βRdαTdαRdβTdβR,where Δdt=dp1,Sii1t+dSii1,Skk2t+dSkk2,l1t−dp2,Sii1t−dSii1,Skk2t−dSkk2,l2t.


### 4.2. Doppler Power Spectrum Density

If we apply the Fourier transform to the space-time correlation function, then we can get the corresponding Doppler power spectrum density (PSD).(16)Sp1l1,p2l2t,f=FτRp1l1LoS,p2l2LoSt,τ+∑i=12FτRp1l1SBi,p2l2SBit,τ+FτRp1l1DB,p2l2DBt,τ.


## 4.1. Space-Time Correlation Function

The normalized space-time correlation function (STCF) between pathp1⟶l1 and p2⟶l2 is defined as(10)Rp1l1,p2l2δT,δR,t,τ=Ehp1l1thp2l2∗t+τEhp1l1t2Ehp2l2t2=K+1Ehp1l1thp2l2∗t+τ.The STCF can be written into(11)Rp1l1,p2l2δT,δR,t,τ=Rp1l1LoS,p2l2LoSδT,δR,t,τ+∑i=12Rp1l1SBi,p2l2SBiδT,δR,t,τ+Rp1l1DB,p2l2DBδT,δR,t,τ.(a)
STCF of LoS component: we also only give the derivation for scenario 1 because it is similar to the scenario 2. The STCF of the LoS component is(12)Rl1p1LoS,l2p2LoSδT,δR,t,τ=Ke−j2π/λΔdt×ej2πτfTDopplerLoSt+fRDopplerLoSt=Ke−j2π/λdp1,l1t+dp2,l2t×ej2πτfTDopplerLoSt+fRDopplerLoSt,wheredp1,l1t and dp2,l2t can be calculated as (9). fTDopplerLoSt and fRDopplerLoSt have been expressed in Section 3.(b)
STCF of SB components: we consider that the EAoD and the AAoD inSB1 and EAoA and AAoA in SB2 are contributed to VMF contribution; thus, the general PDF of these angles is(13)pα,β=kcosβ04πsinhk×ekcosβ0cosβcosα−α0+sinβ0sinβ,whereα,β∈−π,π, and α0,β0∈−π,π. Since there are Nc clusters of scatterers, there are also Nc groups of α0 and β0 corresponding to piα,β (i = 1, …, Nc). Next, if the number of the scatterers tends to infinity, the discrete random variables αT/Rnij and βT/Rnij can be replaced by continuous random variables αT/R and βT/R. Then, the STCF of SBi component becomes(14)Rp1l1SB1,p2l2SB1δT,δR,t,τ=εSBi2∑k=1Nc∫−ππ∫−ππe−j2π/λΔdt×ej2πτfTDopplerSBit×ej2πτfRDopplerSBit×pkαT,βTdαTdβT,whereΔdt=dp1,Skt+dSk,l1t−dp2,Skt−dSk,l2t, and the calculations are similar to (9).(c)
STCF of DB components:(15)Rp1l1DB,p2l2DBδT,δR,t,τ=εDB2∑i=1Nc∑k=1Nc∫−ππ∫−ππ∫−ππ∫−ππe−j2π/λΔdej2πτfTDopplerDBtej2πτfTDopplerDBt×piαT,βTpkαR,βRdαTdαRdβTdβR,where Δdt=dp1,Sii1t+dSii1,Skk2t+dSkk2,l1t−dp2,Sii1t−dSii1,Skk2t−dSkk2,l2t.


## 4.2. Doppler Power Spectrum Density

If we apply the Fourier transform to the space-time correlation function, then we can get the corresponding Doppler power spectrum density (PSD).(16)Sp1l1,p2l2t,f=FτRp1l1LoS,p2l2LoSt,τ+∑i=12FτRp1l1SBi,p2l2SBit,τ+FτRp1l1DB,p2l2DBt,τ.


## 5. Simulation Results and Analysis

In this section, several statistical characteristics of the proposed channel model are simulated and analyzed, including the spatial ACF, temporal ACF, and Doppler PSD. Based on measurement data in [28, 52] and measurement methods in [53–55], the simulation parameters at the initial time are set and listed in Table 2. As for the scenario with slope scenario 2 is chosen in simulations. Due to the energy of the scattering path is much smaller than that contained in the LoS path, the LoS path is ignored in the spatial and temporal ACFs simulation. Besides, since the above measurement campaigns are completed on common roads without slope, in CIR, time-varying ACF, and Doppler PSD simulations, the propagation scenario is reduced to a level load.Table 2
Simulation parameters.
ParameterValueScenario without slopeScenario with slopefc5.9 GHz5.9 GHzK3.7863.786k1515D200 m(D1,D2)(100, 100) m(RT,RR)((15, 20), (15, 20)) m((15, 20), (15, 20)) m(δT,δR)(λ,λ)(λ,λ)(vT,vR)(30, 40) m/s(30, 40) m/svscatterer35 m/s35 m/s(ϕT,ϕR)(0, 0) rad(0, 0.08) rad(γT,γR)(0, 0) rad(0, 0) rad(fTmax,fRmax)(590, 787) Hz(590, 787) Hz(θT,θR)(1.57, 1.57) rad(1.57, 1.57) radαTnij(−2.36, −0.76, 0.76, 2.36) rad(−2.36, −0.76, 0.76, 2.36) radαRnij(−2.36, −0.76, 0.76, 2.36) rad(−2.36, −0.76, 0.76, 2.36) radβTnij(0, 0, 0, 0) rad(0, 0, 0, 0) radβRnij(0, 0, 0, 0) rad(−0.07, 0.07, 0.07, −0.07) rad(αTLoS,βTLoS)(0, 0) rad(0, 0.06) rad(αRLoS,βRLoS)(3.14, 0) rad(3.14, −0.06) radNc44Ni1,Ni2(10, 10)(10, 10)θslope0.08 rad


### 5.1. Simulation of CIR

As mentioned above, in the realistic scenario, the velocity and position vectors of the vehicle are constantly changing; thus, the nonstationary model is closer to reality. The tapped delay line (TDL) model is adopted in the CIR simulation. In Figure5, it can be clearly seen that the locations of different clusters change with time. It needs to be mentioned that the reality scenario is more complicated, since the power of the ray is constantly changing, and new scatterers will appear that accompany existed scatterers’ disappearance.Figure 5
The normalized TDL model without the LoS component.In Figure6, we give a snapshot of the TDL at t=5s. It can be seen that due to the setting of Tx, Rx, and other vehicles’ velocities and locations, different clusters gather together after a certain period of time.Figure 6
A snapshot of the TDL model att=5s.


### 5.2. Spatial ACF

This part is the simulation result of spatial ACF. Figure7 shows the spatial ACFs of scenarios with and without slope at Tx end. The spacing between the two receiving antennas at Rx is fixed, and the antenna spacing in Figure 7 refers to the spacing between adjacent transmitting antennas at Tx. For example, in our simulation, there are Nc=4 clusters of scatterers around the Tx and Rx, respectively, and each clusters contains 10 effective scatterers (Ni1=Nj2=10). Since there are at total Nc×Ni1×Nc×Nj2DB paths, the spatial diversity is huge. As a result, the spatial ACF of the channel is therefore relatively low.Figure 7
Spatial ACFs for the 3D model with and without slope.Figure7 presents normalized spatial ACFs for scenarios including slope with different tilt angles, i.e., θslope=0,θslope=π/36,θslope=π/6,andθslope=π/3. The figure shows that the tilt angle of the slope has little effect on the spatial ACF. As the tilt angle of slope increases, the main lobe of the spatial ACF becomes narrower.


### 5.3. Temporal ACF

#### 5.3.1. The ACF at Initial Time

When simulating the temporal ACF, it is necessary to, respectively, fix the spacing of adjacent transmitting and receiving antennas at Tx and Rx. We takeδT=δR=λ in the simulations. The generation of the scatterers is the same as that in the simulation of the spatial ACF. The specific results are shown in Figure 8.Figure 8
Temporal ACF for the 3D model with and without slope.In the VMF distribution, the number ofα0 is 2×Nc because there are Nc clusters of scatters near Tx and another Nc clusters of scatters near Rx, the values of them are listed in Table 2, and the values of β0 are set to zero.Figure8 shows the temporal ACFs for scenarios with a slope of different tilt angles, and the angles are the same as above. It can be seen from the figure that as the tilt angle of slope increases, the main lobe of the temporal ACF becomes wider.


#### 5.3.2. Time-Varying ACF

In this simulation, we assume that the speeds of the Tx, Rx, and vehicles around them are 30 m/s, 40 m/s, and 35 m/s, respectively. The results of the simulation are shown in Figures9 and 10.Figure 9
Time-varying ACF for the 3D model.Figure 10
ACF att = 0 s, 1 s, 3 s, and 5 s for the 3D model.From Figure9, we can observe the nonstationary feature of the channel. It also verifies that when the speed directions of Tx and Rx coincide, the channel can be regarded as a stationary channel in a short interval (within 1 second). Figure 10 illustrates more detailed time-variant features of the ACFs at different moments, e.g., at the moments t=0s,1s,3s, and 5s. In the simulated scenario set, different vehicles are moving away from each other as time increases. Therefore, the different scattering paths are constantly moving closer together, as shown in Figure 11, the AAoA of SB1 becomes closer to the LoS path (0 radians). And the difference between the angles of the scattering paths and the delays are continuously decreasing. Eventually, it leads to an increase in the main lobe of the ACF.Figure 11
The pdf of AAoA of (a)cluster1, (b)cluster2, (c)cluster3, and (d)cluster4SB1 path at t = 0, 1, 3, and 5 s.
(a)(b)(c)(d)


### 5.4. Doppler PSD

The simulated Doppler PSDs are shown in Figures12 and 13. The result of Figure 12 is consistent with Figure 5. As time increases, Tx and Rx keep moving away from each other, and different clusters gather with each other. It results in a decrease in the width of Doppler frequency.Figure 12
Doppler PSD for the 3D model.Figure 13
Doppler PSD for the 3D model att=0s.We haveNi1=Nj2=10, so there should be a total of N clusters of impulses in Figure 13, i.e., N=N1+N2+N3, N1=Nc×Ni1 and N2=Nc×Nj2, N3=Nc×Ni1×Nc×Nj2. However, in Figure 13, only N1+N2 clusters of SB path impulses can be seen because the energy of DB path is too low. Furthermore, in each cluster of scatterers, all scatterers are surrounding the same direction, so as it can be seen in Figure 13 that these impulses are in the form of clusters. Besides, the Doppler PSD provides help to verify the correctness of the model. After comparison and verification, the Doppler PSDs shown in Figures 12 and 13 are consistent with the theoretical analysis.


### 5.5. Comparison with Measurements

In this section, we compare the simulation results of the model with measurement data to verify the accuracy of the model. We choose the parameter and fitting cure given in [52] and adjust the vehicle’s distribution and movement in the model to match the measurement scenario. After adjusting the model, we simulate RMS delay spread, RMS Doppler spread, and RMSE to be compared with the measurement data.


#### 5.5.1. Adjustments of Scenario

The adjusted scenario is shown in Figure14. Since the measurement environments are level roads, we set the tilt angle of slope to zero. At the initial time, Tx and Rx are moving on two adjacent lanes, then Rx accelerates and overtakes Tx. The velocity of different vehicles varies according to different scenarios.Figure 14
Adjusted scenario.In [52], measurement campaigns are in four different scenarios, highway, urban, suburb, and municipal lake. In the lake scenario, there is very little scattering caused by the environment because there are almost no buildings in this scenario, and the small-scale fading of the channel mainly comes from the scattering caused by the surrounding vehicles. In other measurement scenarios, the environment causes severe channel fading. Therefore, we added environment scatterers in the simulation. The distribution of environment scatterers varies according to the scenario. In suburb and urban scenarios, the environment scatterers are distributed on both sides of the lane. The main difference is that there are dense trees on both sides of the road in suburb areas, as for urban areas, there are mainly dense large buildings on both sides of the lane. In the highway scenario, the environment scatterers are only distributed on one side, and they are mainly sparse houses.Besides, it can be divided into overtaking (OT) and nonovertaking (NOT) cases according to the distance between Tx and Rx. In the case of OT, the distance between Tx and Rx is very close, less than 11 m. While in the NOT case, they keep a certain distance.


#### 5.5.2. RMS Delay Spread

In wireless communication, RMS delay spread is a statistical parameter that can describe the delay characteristics of the radio channel. It is often used in the evaluation of channel characteristics. In V2V communication, the multipath effect caused by the scattering between vehicles and scattering from the environment caused small-scale fading. RMS delay spread can reflect this multipath effect.The complex time-varying CIR can be derived by an inverse Fourier transform of the channel transfer function (6a)–(6c). Then, we can calculate the time-variant power delay profile (PDP) by averaging the magnitude squared of the CIR.(17)Pt,τ=ht,τ2=∑iait2δτ−τi,where ai is the complex coefficient of each delay path, and τi is the excess delay of the ith path.Then, we can get the RMS delayTrms as follows:(18a)Tm=1∫−∞∞PAτdτ⋅∫−∞∞PAττdτ,(18b)Trms=1∫−∞∞PAτdτ⋅∫−∞∞PAττ2dτ−Tm2,where PAτ is the averaged power delay profile.The comparison results are shown in Figure15. The fitting curves corresponding to different scenarios are given in [52]. We put simulation results, measurement data, and fitting curves together for comparison. In the following section, RMSE will be used for evaluation.Figure 15
The CDF of RMS delay spread for (a) highway, (b) urban, (c) suburb, and (d) municipal lake scenario.
(a)(b)(c)(d)


#### 5.5.3. RMS Doppler Spread

In V2V communication, vehicles are moving at high speed, and the speed direction of vehicles may change, so it causes a strong Doppler effect. Therefore, RMS Doppler spread is also applied to the description of V2V channel delay characteristics.We can get the RMS Doppler spreadDrms as(19a)PBm=∫−∞∞PBνdν,(19b)Dm=∫−∞∞PBννdνPBm,(19c)Drms=∫−∞∞PBνν2dνPBm−Dm2,where PBm represents the integration of scattering function.The comparison results are shown in Figure16. The comparison of simulation results, measurement data, and fitting curves are given in the figure.Figure 16
The CDF of RMS Doppler spread for (a) highway, (b) urban, (c) suburb, and (d) municipal lake scenario.
(a)(b)(c)(d)


#### 5.5.4. RMSE

To verify that the model can fit the measurement data better than the fitting curve, we calculate the RMSE of delay and Doppler spread.(20)RMSE=1n∑i=1nxi−yi2,where yi denotes the measurement data, and xi denotes the simulation result or the fitting curve.RMSE can reflect the fitting degree of simulation results and measurement data. A smaller RMSE means a better fitting. The calculation results are shown in Table3. When the difference of RMSE is negative, it indicates that the model can better fit the measurement data than the fitting curve. As can be seen in the table, the model has a smaller RMSE in most cases, and the model has a similar degree to the fitting curve in the remaining scenarios.Table 3
Data of RMSE.
ScenarioHighwayUrbanSuburbLakeCaseOTNOTOTNOTOTNOTOTNOTDelay RMSE between simulation result and measurement dataRMSE0.09900.07290.14640.08640.06220.15350.06080.0698Delay RMSE between fitting distribution and measurement dataRMSE0.10580.07010.04760.15160.08160.08470.06680.1743Difference of delay RMSERMSE−0.00680.00280.0988−0.0652−0.01930.0688−0.0059−0.1044Doppler RMSE between simulation result and measurement dataRMSE0.10070.09010.14630.06940.13560.08670.10640.1233Doppler RMSE between fitting distribution and measurement dataRMSE0.16600.09880.25450.12380.15030.07200.10610.2210Difference of Doppler RMSERMSE−0.0653−0.0087−0.1082−0.0544−0.01460.01470.0003−0.0977


## 5.1. Simulation of CIR

As mentioned above, in the realistic scenario, the velocity and position vectors of the vehicle are constantly changing; thus, the nonstationary model is closer to reality. The tapped delay line (TDL) model is adopted in the CIR simulation. In Figure5, it can be clearly seen that the locations of different clusters change with time. It needs to be mentioned that the reality scenario is more complicated, since the power of the ray is constantly changing, and new scatterers will appear that accompany existed scatterers’ disappearance.Figure 5
The normalized TDL model without the LoS component.In Figure6, we give a snapshot of the TDL at t=5s. It can be seen that due to the setting of Tx, Rx, and other vehicles’ velocities and locations, different clusters gather together after a certain period of time.Figure 6
A snapshot of the TDL model att=5s.


## 5.2. Spatial ACF

This part is the simulation result of spatial ACF. Figure7 shows the spatial ACFs of scenarios with and without slope at Tx end. The spacing between the two receiving antennas at Rx is fixed, and the antenna spacing in Figure 7 refers to the spacing between adjacent transmitting antennas at Tx. For example, in our simulation, there are Nc=4 clusters of scatterers around the Tx and Rx, respectively, and each clusters contains 10 effective scatterers (Ni1=Nj2=10). Since there are at total Nc×Ni1×Nc×Nj2DB paths, the spatial diversity is huge. As a result, the spatial ACF of the channel is therefore relatively low.Figure 7
Spatial ACFs for the 3D model with and without slope.Figure7 presents normalized spatial ACFs for scenarios including slope with different tilt angles, i.e., θslope=0,θslope=π/36,θslope=π/6,andθslope=π/3. The figure shows that the tilt angle of the slope has little effect on the spatial ACF. As the tilt angle of slope increases, the main lobe of the spatial ACF becomes narrower.


## 5.3. Temporal ACF

### 5.3.1. The ACF at Initial Time

When simulating the temporal ACF, it is necessary to, respectively, fix the spacing of adjacent transmitting and receiving antennas at Tx and Rx. We takeδT=δR=λ in the simulations. The generation of the scatterers is the same as that in the simulation of the spatial ACF. The specific results are shown in Figure 8.Figure 8
Temporal ACF for the 3D model with and without slope.In the VMF distribution, the number ofα0 is 2×Nc because there are Nc clusters of scatters near Tx and another Nc clusters of scatters near Rx, the values of them are listed in Table 2, and the values of β0 are set to zero.Figure8 shows the temporal ACFs for scenarios with a slope of different tilt angles, and the angles are the same as above. It can be seen from the figure that as the tilt angle of slope increases, the main lobe of the temporal ACF becomes wider.


### 5.3.2. Time-Varying ACF

In this simulation, we assume that the speeds of the Tx, Rx, and vehicles around them are 30 m/s, 40 m/s, and 35 m/s, respectively. The results of the simulation are shown in Figures9 and 10.Figure 9
Time-varying ACF for the 3D model.Figure 10
ACF att = 0 s, 1 s, 3 s, and 5 s for the 3D model.From Figure9, we can observe the nonstationary feature of the channel. It also verifies that when the speed directions of Tx and Rx coincide, the channel can be regarded as a stationary channel in a short interval (within 1 second). Figure 10 illustrates more detailed time-variant features of the ACFs at different moments, e.g., at the moments t=0s,1s,3s, and 5s. In the simulated scenario set, different vehicles are moving away from each other as time increases. Therefore, the different scattering paths are constantly moving closer together, as shown in Figure 11, the AAoA of SB1 becomes closer to the LoS path (0 radians). And the difference between the angles of the scattering paths and the delays are continuously decreasing. Eventually, it leads to an increase in the main lobe of the ACF.Figure 11
The pdf of AAoA of (a)cluster1, (b)cluster2, (c)cluster3, and (d)cluster4SB1 path at t = 0, 1, 3, and 5 s.
(a)(b)(c)(d)


## 5.3.1. The ACF at Initial Time

When simulating the temporal ACF, it is necessary to, respectively, fix the spacing of adjacent transmitting and receiving antennas at Tx and Rx. We takeδT=δR=λ in the simulations. The generation of the scatterers is the same as that in the simulation of the spatial ACF. The specific results are shown in Figure 8.Figure 8
Temporal ACF for the 3D model with and without slope.In the VMF distribution, the number ofα0 is 2×Nc because there are Nc clusters of scatters near Tx and another Nc clusters of scatters near Rx, the values of them are listed in Table 2, and the values of β0 are set to zero.Figure8 shows the temporal ACFs for scenarios with a slope of different tilt angles, and the angles are the same as above. It can be seen from the figure that as the tilt angle of slope increases, the main lobe of the temporal ACF becomes wider.


## 5.3.2. Time-Varying ACF

In this simulation, we assume that the speeds of the Tx, Rx, and vehicles around them are 30 m/s, 40 m/s, and 35 m/s, respectively. The results of the simulation are shown in Figures9 and 10.Figure 9
Time-varying ACF for the 3D model.Figure 10
ACF att = 0 s, 1 s, 3 s, and 5 s for the 3D model.From Figure9, we can observe the nonstationary feature of the channel. It also verifies that when the speed directions of Tx and Rx coincide, the channel can be regarded as a stationary channel in a short interval (within 1 second). Figure 10 illustrates more detailed time-variant features of the ACFs at different moments, e.g., at the moments t=0s,1s,3s, and 5s. In the simulated scenario set, different vehicles are moving away from each other as time increases. Therefore, the different scattering paths are constantly moving closer together, as shown in Figure 11, the AAoA of SB1 becomes closer to the LoS path (0 radians). And the difference between the angles of the scattering paths and the delays are continuously decreasing. Eventually, it leads to an increase in the main lobe of the ACF.Figure 11
The pdf of AAoA of (a)cluster1, (b)cluster2, (c)cluster3, and (d)cluster4SB1 path at t = 0, 1, 3, and 5 s.
(a)(b)(c)(d)


## 5.4. Doppler PSD

The simulated Doppler PSDs are shown in Figures12 and 13. The result of Figure 12 is consistent with Figure 5. As time increases, Tx and Rx keep moving away from each other, and different clusters gather with each other. It results in a decrease in the width of Doppler frequency.Figure 12
Doppler PSD for the 3D model.Figure 13
Doppler PSD for the 3D model att=0s.We haveNi1=Nj2=10, so there should be a total of N clusters of impulses in Figure 13, i.e., N=N1+N2+N3, N1=Nc×Ni1 and N2=Nc×Nj2, N3=Nc×Ni1×Nc×Nj2. However, in Figure 13, only N1+N2 clusters of SB path impulses can be seen because the energy of DB path is too low. Furthermore, in each cluster of scatterers, all scatterers are surrounding the same direction, so as it can be seen in Figure 13 that these impulses are in the form of clusters. Besides, the Doppler PSD provides help to verify the correctness of the model. After comparison and verification, the Doppler PSDs shown in Figures 12 and 13 are consistent with the theoretical analysis.


## 5.5. Comparison with Measurements

In this section, we compare the simulation results of the model with measurement data to verify the accuracy of the model. We choose the parameter and fitting cure given in [52] and adjust the vehicle’s distribution and movement in the model to match the measurement scenario. After adjusting the model, we simulate RMS delay spread, RMS Doppler spread, and RMSE to be compared with the measurement data.


### 5.5.1. Adjustments of Scenario

The adjusted scenario is shown in Figure14. Since the measurement environments are level roads, we set the tilt angle of slope to zero. At the initial time, Tx and Rx are moving on two adjacent lanes, then Rx accelerates and overtakes Tx. The velocity of different vehicles varies according to different scenarios.Figure 14
Adjusted scenario.In [52], measurement campaigns are in four different scenarios, highway, urban, suburb, and municipal lake. In the lake scenario, there is very little scattering caused by the environment because there are almost no buildings in this scenario, and the small-scale fading of the channel mainly comes from the scattering caused by the surrounding vehicles. In other measurement scenarios, the environment causes severe channel fading. Therefore, we added environment scatterers in the simulation. The distribution of environment scatterers varies according to the scenario. In suburb and urban scenarios, the environment scatterers are distributed on both sides of the lane. The main difference is that there are dense trees on both sides of the road in suburb areas, as for urban areas, there are mainly dense large buildings on both sides of the lane. In the highway scenario, the environment scatterers are only distributed on one side, and they are mainly sparse houses.Besides, it can be divided into overtaking (OT) and nonovertaking (NOT) cases according to the distance between Tx and Rx. In the case of OT, the distance between Tx and Rx is very close, less than 11 m. While in the NOT case, they keep a certain distance.


### 5.5.2. RMS Delay Spread

In wireless communication, RMS delay spread is a statistical parameter that can describe the delay characteristics of the radio channel. It is often used in the evaluation of channel characteristics. In V2V communication, the multipath effect caused by the scattering between vehicles and scattering from the environment caused small-scale fading. RMS delay spread can reflect this multipath effect.The complex time-varying CIR can be derived by an inverse Fourier transform of the channel transfer function (6a)–(6c). Then, we can calculate the time-variant power delay profile (PDP) by averaging the magnitude squared of the CIR.(17)Pt,τ=ht,τ2=∑iait2δτ−τi,where ai is the complex coefficient of each delay path, and τi is the excess delay of the ith path.Then, we can get the RMS delayTrms as follows:(18a)Tm=1∫−∞∞PAτdτ⋅∫−∞∞PAττdτ,(18b)Trms=1∫−∞∞PAτdτ⋅∫−∞∞PAττ2dτ−Tm2,where PAτ is the averaged power delay profile.The comparison results are shown in Figure15. The fitting curves corresponding to different scenarios are given in [52]. We put simulation results, measurement data, and fitting curves together for comparison. In the following section, RMSE will be used for evaluation.Figure 15
The CDF of RMS delay spread for (a) highway, (b) urban, (c) suburb, and (d) municipal lake scenario.
(a)(b)(c)(d)


### 5.5.3. RMS Doppler Spread

In V2V communication, vehicles are moving at high speed, and the speed direction of vehicles may change, so it causes a strong Doppler effect. Therefore, RMS Doppler spread is also applied to the description of V2V channel delay characteristics.We can get the RMS Doppler spreadDrms as(19a)PBm=∫−∞∞PBνdν,(19b)Dm=∫−∞∞PBννdνPBm,(19c)Drms=∫−∞∞PBνν2dνPBm−Dm2,where PBm represents the integration of scattering function.The comparison results are shown in Figure16. The comparison of simulation results, measurement data, and fitting curves are given in the figure.Figure 16
The CDF of RMS Doppler spread for (a) highway, (b) urban, (c) suburb, and (d) municipal lake scenario.
(a)(b)(c)(d)


### 5.5.4. RMSE

To verify that the model can fit the measurement data better than the fitting curve, we calculate the RMSE of delay and Doppler spread.(20)RMSE=1n∑i=1nxi−yi2,where yi denotes the measurement data, and xi denotes the simulation result or the fitting curve.RMSE can reflect the fitting degree of simulation results and measurement data. A smaller RMSE means a better fitting. The calculation results are shown in Table3. When the difference of RMSE is negative, it indicates that the model can better fit the measurement data than the fitting curve. As can be seen in the table, the model has a smaller RMSE in most cases, and the model has a similar degree to the fitting curve in the remaining scenarios.Table 3
Data of RMSE.
ScenarioHighwayUrbanSuburbLakeCaseOTNOTOTNOTOTNOTOTNOTDelay RMSE between simulation result and measurement dataRMSE0.09900.07290.14640.08640.06220.15350.06080.0698Delay RMSE between fitting distribution and measurement dataRMSE0.10580.07010.04760.15160.08160.08470.06680.1743Difference of delay RMSERMSE−0.00680.00280.0988−0.0652−0.01930.0688−0.0059−0.1044Doppler RMSE between simulation result and measurement dataRMSE0.10070.09010.14630.06940.13560.08670.10640.1233Doppler RMSE between fitting distribution and measurement dataRMSE0.16600.09880.25450.12380.15030.07200.10610.2210Difference of Doppler RMSERMSE−0.0653−0.0087−0.1082−0.0544−0.01460.01470.0003−0.0977


## 5.5.1. Adjustments of Scenario

The adjusted scenario is shown in Figure14. Since the measurement environments are level roads, we set the tilt angle of slope to zero. At the initial time, Tx and Rx are moving on two adjacent lanes, then Rx accelerates and overtakes Tx. The velocity of different vehicles varies according to different scenarios.Figure 14
Adjusted scenario.In [52], measurement campaigns are in four different scenarios, highway, urban, suburb, and municipal lake. In the lake scenario, there is very little scattering caused by the environment because there are almost no buildings in this scenario, and the small-scale fading of the channel mainly comes from the scattering caused by the surrounding vehicles. In other measurement scenarios, the environment causes severe channel fading. Therefore, we added environment scatterers in the simulation. The distribution of environment scatterers varies according to the scenario. In suburb and urban scenarios, the environment scatterers are distributed on both sides of the lane. The main difference is that there are dense trees on both sides of the road in suburb areas, as for urban areas, there are mainly dense large buildings on both sides of the lane. In the highway scenario, the environment scatterers are only distributed on one side, and they are mainly sparse houses.Besides, it can be divided into overtaking (OT) and nonovertaking (NOT) cases according to the distance between Tx and Rx. In the case of OT, the distance between Tx and Rx is very close, less than 11 m. While in the NOT case, they keep a certain distance.


## 5.5.2. RMS Delay Spread

In wireless communication, RMS delay spread is a statistical parameter that can describe the delay characteristics of the radio channel. It is often used in the evaluation of channel characteristics. In V2V communication, the multipath effect caused by the scattering between vehicles and scattering from the environment caused small-scale fading. RMS delay spread can reflect this multipath effect.The complex time-varying CIR can be derived by an inverse Fourier transform of the channel transfer function (6a)–(6c). Then, we can calculate the time-variant power delay profile (PDP) by averaging the magnitude squared of the CIR.(17)Pt,τ=ht,τ2=∑iait2δτ−τi,where ai is the complex coefficient of each delay path, and τi is the excess delay of the ith path.Then, we can get the RMS delayTrms as follows:(18a)Tm=1∫−∞∞PAτdτ⋅∫−∞∞PAττdτ,(18b)Trms=1∫−∞∞PAτdτ⋅∫−∞∞PAττ2dτ−Tm2,where PAτ is the averaged power delay profile.The comparison results are shown in Figure15. The fitting curves corresponding to different scenarios are given in [52]. We put simulation results, measurement data, and fitting curves together for comparison. In the following section, RMSE will be used for evaluation.Figure 15
The CDF of RMS delay spread for (a) highway, (b) urban, (c) suburb, and (d) municipal lake scenario.
(a)(b)(c)(d)


## 5.5.3. RMS Doppler Spread

In V2V communication, vehicles are moving at high speed, and the speed direction of vehicles may change, so it causes a strong Doppler effect. Therefore, RMS Doppler spread is also applied to the description of V2V channel delay characteristics.We can get the RMS Doppler spreadDrms as(19a)PBm=∫−∞∞PBνdν,(19b)Dm=∫−∞∞PBννdνPBm,(19c)Drms=∫−∞∞PBνν2dνPBm−Dm2,where PBm represents the integration of scattering function.The comparison results are shown in Figure16. The comparison of simulation results, measurement data, and fitting curves are given in the figure.Figure 16
The CDF of RMS Doppler spread for (a) highway, (b) urban, (c) suburb, and (d) municipal lake scenario.
(a)(b)(c)(d)


## 5.5.4. RMSE

To verify that the model can fit the measurement data better than the fitting curve, we calculate the RMSE of delay and Doppler spread.(20)RMSE=1n∑i=1nxi−yi2,where yi denotes the measurement data, and xi denotes the simulation result or the fitting curve.RMSE can reflect the fitting degree of simulation results and measurement data. A smaller RMSE means a better fitting. The calculation results are shown in Table3. When the difference of RMSE is negative, it indicates that the model can better fit the measurement data than the fitting curve. As can be seen in the table, the model has a smaller RMSE in most cases, and the model has a similar degree to the fitting curve in the remaining scenarios.Table 3
Data of RMSE.
ScenarioHighwayUrbanSuburbLakeCaseOTNOTOTNOTOTNOTOTNOTDelay RMSE between simulation result and measurement dataRMSE0.09900.07290.14640.08640.06220.15350.06080.0698Delay RMSE between fitting distribution and measurement dataRMSE0.10580.07010.04760.15160.08160.08470.06680.1743Difference of delay RMSERMSE−0.00680.00280.0988−0.0652−0.01930.0688−0.0059−0.1044Doppler RMSE between simulation result and measurement dataRMSE0.10070.09010.14630.06940.13560.08670.10640.1233Doppler RMSE between fitting distribution and measurement dataRMSE0.16600.09880.25450.12380.15030.07200.10610.2210Difference of Doppler RMSERMSE−0.0653−0.0087−0.1082−0.0544−0.01460.01470.0003−0.0977


## 6. Conclusion

In this study, we developed an analytical multicluster nonstationary 3D MIMO channel model for the V2V communication system, which is suitable for the scenario where there are a certain number of vehicles near Tx and Rx on or near a slope. After simplification, the model can be applied to the normal road without slope. In addition, we also considered the impact of vehicle changing lanes. Furthermore, we derived and simulated statistical properties of the proposed channel such as the spatial ACF, temporal ACF, time-varying ACF, and Doppler PSD. Finally, we compared the simulation results with the measurement data in four scenarios. The comparison results show that the proposed model can mimic the measurement data better than the fitting curve. These simulation results have shown the utility of the proposed model.



---

*Source: 8871255-2021-01-13.xml*