Spatial and temporal characteristics of the propagation channel have a significant influence on multiantenna method applicability for fifth-generation- (5G-) enabled Internet of Things (IoT). In this paper, the statistical characteristics of a novel three-dimensional (3D) geometric-based stochastic model for next-generation vehicle-to-vehicle (V2V) multiple-input multiple-output (MIMO) communications under the nonisotropic scattering environment are investigated. In both line-of-sight (LoS) and non-line-of-sight (NLoS) conditions, the proposed model investigates the spatial, frequency, and temporal domain statistical distribution of multipath received signals by using the time-variant transfer function for indoor environments. The probability density function (PDF) of separation distance between the transceiver antennas, angle-of-arrival (AoA), and angle-of-departure (AoD) in the azimuth and elevation planes is derived by using closed-form expressions. For the space, time, and frequency correlation function (STF-CF), a precise analytical expression is derived based on MIMO antenna system. We further determine the effects of several model parameters on the V2V channel performance, such as tunnel width, antenna array spacing, Ricean -factor, and moving velocity. The statistical characteristics of the MIMO channel model are validated by simulation results, confirming the flexibility and effectiveness of our proposed model in the tunnel scenario.

1. Introduction

Recently, the wireless communication technologies like multiple input multiple output (MIMO), channel coding, cooperative communication, and internet of vehicles (IoV) have played a crucial role towards the development of future communication networks [110]. MIMO channel, which uses hundreds or even more antennas, has received a lot of attention. It has been demonstrated that MIMO technique is capable to improve the energy and spectrum efficiency, especially in rich scattering environment. The rapid advancement of the 5G-enabled Internet-of-Things (IoT) communication systems for vehicle-to-vehicle (V2V) radio communications has acquired much more attentions [1114]. Unlike traditional fixed to mobile cellular networks, the V2V channel uses low elevation for multiple antennas when both the transmitting and receiving antennas are conceived in motion. The radio propagation parameters between the transmitting () and receiving () antennas must be designed to facilitate the thorough investigation of V2V radio communications [15]. Therefore, the realistic propagation channel model is required which can provide a quick and easy algorithm to approximate the statistical characteristics of V2V radio channel [16, 17].

To determine the optimized performance of the V2V radio communications in tunnel environments, precise channel models are under consideration. The V2V channel models can be classified as stochastic models and deterministic models (mostly based on the ray tracing (RT) method). Stochastic models are classified into two types: regular-shaped geometrical-based stochastic models (RS-GBSMs) [1822] and non-regular-shaped geometrical-based stochastic models (NGBSMs). The latter, also named as parametric models, are established from channel measurements, whereas the former are constructed from the regular geometrical-shape of the scattering objects. In [18], it is proved that the line-of-sight (LoS) path between the transceiver antennas is likely to be interrupted by buildings and barriers. As a result, Ricean channels must be used to investigate V2V communication systems. In [19], a precise visual scattering model is proposed for the multibounced V2V propagation channels in crowded urban communication environment; however, they did not address the influence of moving car velocity on the V2V channel propagation characteristics. In addition, [20, 21] introduce the 3D RS-GBSMs for microcell and macrocell communication systems, respectively.

However, the majority of the RSGBSMs models are based on narrowband channel assumptions, where the rays have almost the same propagation delay [22], which is not an accurate portrayal of V2V wireless communication system. In [23], the narrowband and wideband V2V channel characteristics based on propagation delays have been investigated through the real-time measurement campaign. The wide-sense stationary uncorrelated scattering (WSSUS) assumptions are used mostly in previous channel models; however, it is only applicable for short time intervals [24]. The WSSUS means actually first- and second-order stationarity in both time and frequency domains (wide-sense stationarity in time, uncorrelated scattering in delay) and holds in successive time and frequency interval. According to the measurement campaign in [25, 26], the stationary interval, in which the WSSUS assumption is still valid for V2V environments, is substantially smaller than the observation interval. This implies that in V2V channels, the WSSUS assumptions are not being fulfilled in most of the cases. Therefore, the stochastic channel models must need to consider the channel’s nonstationarity properties to compensate these shortcomings.

Similarly, in [27], the algebraic framework of the transceiver antennas in the Doppler spectrum based on moving velocity and arbitrary delays has been investigated and velocity vectors are considered to characterize the arbitrary moving direction in V2V radio communication environment. In [28], a novel 3D-GBSM model based on an ellipsoid Gaussian scattering distribution is proposed that can jointly characterize the elevation angles, azimuth angles, and distances from the transmitter (receiver) to the first (last)-bounce scatterers. The spatial consistency of the channel is supported by tracking the locations of the transmitter, receiver, and scatterers. Reference [29] proposes a statistical wideband V2V MIMO channel model based on the Unitary matrix transformation algorithm. The V2V MIMO channel model for other taps is estimated by examining the propagation delays of the first tap with a Unitary matrix transformation algorithm. In [30], a nonstationary stochastic model has been presented between the vehicle and terminal using 3D antenna arrays and random trajectories with fixed velocity consideration. In [31], the delay-based Doppler PDF for V2V radio channels is computed by using a prolate spheroidal system. In [32], the authors proposed a V2V scattering model with the assumption of randomly distributed scatterers across the semicircular shaped tunnel sidewalls. In the context of randomly moving clusters, [33] proposes a nonstationary GBSM channel model for future intelligent V2V communications. Because of the nonstationarity induced by moving clusters, the clusters’ movements are described using a random walk procedure.

Also, [34] presents a regular-shaped 3D-GBSM model for V2V MIMO system for nonisotropic fading channel. However, the channel statistics of PDF, time delays, and tunnel width impact at different locations, which is critical for V2V MIMO systems, have not been investigated in the given model. Therefore, in this paper, the 3D MIMO V2V channel is modeled using the rectangular single-bounce scattering tunnel (RSBST) model, rather than the geometrical semicircular or circular tunnel propagation models. The inspiration came from the idea that both of the geometrical rectangular and semicircular tunnels are generally used in MIMO systems. As arched and rectangular tunnel environments have different propagation characteristics, therefore, we feel that our proposed model is simple and more suitable for indoor rectangular tunnel environment. To our knowledge, no investigation of radio wave propagation in rectangular tunnel from the perspective of scatterer distribution and movement of V2V MIMO systems along tunnel width has been performed.

In this paper, we investigate the new geometrical RSBST model and compute its geometrical statistics with other scattering models to examine the influence of size and cross-sectional shape of tunnels on the V2V radio channel configurations. The key contribution of this study is to analyze the effect of different RSBST model parameters on V2V channel statistics and to obtain an optimized expression for scatterer distribution. The distribution of scatterers has been derived by considering the -axis influence on the system performance and computed the 3D STF-CF of the V2V propagation channel in LoS and NLoS propagation environments. The impacts of numerous model parameters, such as tunnel diameter, antenna array spacing, Ricean -factor, and moving velocity, on the V2V channel statistics are investigated as well.

The remainder of this paper is arranged in the following manner. Section 2 explains the proposed RSBST system model. In Section 3, different characteristics of the 3D-GBSM channel, such as TVTF, STF-CF, and power angular spread, are derived. In Section 4, the analytical results are presented and critically analyzed. Finally, Section 5 summarizes the paper.

2. Analysis of Geometric RSBST Model

The geometrical RSBST propagation model for the V2V MIMO system inside a rectangular-shaped tunnel is described in this section. Figure 1 represents a typical propagation scenario in an -length tunnel. As illustrated in Figure 2, we assume that the scatterers are distributed randomly beside the tunnel sidewalls and ceiling. Due to cross-sectional size, shape, and geometry, the characteristics of a rectangular tunnel are substantially different from those of a circular or arch tunnel [35]. For example, in an arch tunnel, the tunnel width () is dependent on the tunnel radius due to semicircular geometry, whereas in a rectangular-shaped tunnel, the width () of the tunnel is independent from the height () of tunnel. As a result, when the scatterers are distributed on the tunnel ceiling and sidewalls, the -axis affects the statistics of the V2V radio channel in a rectangular tunnel. Thus, we consider the impact of tunnels’ -axis when computing the STF-CF for next-generation V2V MIMO channel configurations in our proposed 3D RSBST model. In addition, Figure 3 shows the geometrical representation of RSBST V2V scattering model.

As illustrated in Figure 3, we employ a Cartesian coordinate approach to investigate the distribution of random scatterers on the tunnel sidewalls and ceiling for , , , , and , . On the rectangle tunnel sidewalls and ceiling, we suppose that the infinite number of scattering components (, , , ) is distributed randomly [36].

However, the WSSUS assumptions can be still fulfilled for V2V radio channels over the small time periods [26]. Therefore, the V2V radio channel inside the tunnel is assumed to be WSSUS in our analysis. The locations of the receiver () and transmitter () are defined by cartesian coordinates and , when and are both moving in the opposite or same directions with velocity . and determine the angle of motion for and antenna elements, respectively. Additionally, the and are equipped with a multiple and antenna elements, and the antenna element spacing is represented by and , respectively. The path length between the th transmitting antenna and the scatterer is denoted by , whereas the distance between the th receiving antenna and the scatterer is denoted by . Moreover, the specifies the distance between the th element of transmitting array and the th element of receiving array in terms of LoS path. In our work, we assume single-bounce scattering (SB) approach and further analyze the scattering effects due to randomly distributed scatterers between transceiver antennas, also known as effective scattering elements. The elevation angles and orientation of the () antenna array elements with respect to the plane can be represented by and , respectively. The symbols and signify the elevation angle of departure (EAoD) and azimuth angle of departure (AAoD) of the signals that are impinging on the scattering elements, whereas the angles and represent the elevation angle of arrival (EAoA) and azimuth angle of arrival (AAoA) of the signals which are scattered from scattering elements, respectively. The EAoD, EAoA, AAoD, and AAoA in coordinate system of scattering elements are given as follows: where or and, the index relates to the receiver or transmitter. The functions and are defined as follows:

3. Statistical Characteristics of RSBST V2V Model

The reference model for the V2V MIMO channel under direct (LoS) and diffuse (NLoS) propagation conditions is presented here. The time-variant transfer function (TVTF), STF-CF, and power angular spread will be calculated from the geometrical RSBST scattering model and provided as a sum of LoS and diffuse components.

3.1. TVTF

The TVTF is denoted as , which is the sum of LoS and diffuse components as where and represent the TVTF of LoS and diffuse components, respectively. The of LoS components is written as follows: where

In (8), the Doppler shifts of the LoS components generated by the movement of the and antennas are denoted by and , respectively. The angles , , , and in (9) and (10) represent the EAoD, AAoD, EAoA, and AAoA of the LoS component, respectively. The in (5) provides the Ricean -factor, which is defined as the mean power ratio of LoS component with the diffuse components. In (5), represents the propagation delays of LoS components and given by , where is speed of light. The of diffuse components for the transmission link can be determined as where where the symbols , , , and represent the maximum Doppler frequency of and , Ricean -factor, and the delays of diffuse components, respectively. In addition, and denote the distances between scatterer to th transmitting antenna array and from the scatterer to the th receiving antenna array components , respectively.

3.2. STF-CF

The STF-CF is the correlation of and in (4), which can be calculated as where represents the conjugate value operator and denotes the expected value operator. and represent the STF-CF of the LoS and the diffuse components, respectively. Thus, the STF-CF can be calculated as follows: where

Moreover, can be estimated as

Note that the propagation delays , and the Doppler frequencies are based on and coordinates. Due to the infinite number of scattering elements in proposed V2V model, the random coordinates and can be replaced by independent variables and . Furthermore, the and are assumed to be uniformly distributed, and their PDF can be estimated as follows:

The expression can be used to calculate the distribution of variable , where is randomly distributed in the range of . By considering the transformation of given variables, the PDF of variable can be derived as follows:

The of variable for different widths is shown in Figure 4. We can see that the of in (23) drops down the uniform value within the range of , as the width () of tunnel increases.

Thus, by considering PDFs of randomly distributed variables and , the joint PDF can be calculated as

Similarly, by exploiting the joint PDF , the STF-CF of diffuse components in (20) can be rewritten as where

It should be noted that the propagation angles , , , and in (27) and (28) are dependent on the , and coordinates of .

3.3. Power Angular Spread

The power angular spread (PAS) can be configured by exploiting the Bartlett method [37]. For electromagnetic analysis, the Bartlett beamforming technique is used to estimate the angle of departure and angle of arrival in azimuth and elevation planes. The Bartlett algorithm is predicated based on the steering vector . The Bartlett beamforming algorithm is commonly known as the Fourier spectrum analyzing approach. The objective is to discover the weighting vector that maximizes the power of the incoming signal. The receiving antenna array elements can receive the signals from different spatially distributed users. Furthermore, the received power includes both LoS and reflected (NLoS) path signals, all of which appear to be emanating from different directions and angles. The steering vector of the transmitted signals can be estimated as where , is the wavelength, is the width of antenna array, and represents the length of antenna arrays, respectively. Also, , , , and denote the antenna array spacing along the -axis, antenna array spacing along the -axis, elevation angle along the positive -axis, and azimuth angle along the positive axis. Furthermore, the signal vector of the arrays can be written as where denotes the noise power. Also, if sources have the same time slot and frequency, the received signals can be written as follows:

By assuming that the received signals are coming through the azimuth direction, so the output of the array by Bartlett can be calculated as follows: where represents the noise power variance and specifies the transpose conjugate of the given vector. The significant solution of (29) is as follows: where denotes the normalized vector; thus, the weighting vector of the Bartlett algorithm becomes

This means that the weighting vector is the same as the incident wave spatial characteristics. The final expression for estimating the PAS is given here only. Reference [37] contains further information about the Bartlett algorithm. The PAS of can be estimated as where is the covariance matrix of received signals, and it may be calculated as where denotes the received signals at the th delay bin and is the linear transformation operator, which sequentially turns the matrix columns into a single column vector. If the normalized vector is used, the total PAS can be stated as

In V2V scenario, the geometric path length is usually time-variant due to the relative movement between the and antenna arrays. Compared to several existing 3D and 2D fixed-to-mobile (F2M) and V2V channel models, the proposed channel model can be utilized more efficiently in tunnel environment; for example, as the relative time is equal to zero, the proposed model tends to be wide-sense-stationary (WSS) channel, which is similar to the F2M channel [38]. However, the proposed model represents nonstationary V2V channels where the relative time is not equal to zero, as described in [39]. Moreover, the proposed 3D model can reflects a wide range of communication scenarios by modifying the geometrical model parameters, as indicated in Sections 2 and 3. For example, the proposed channel model can represent the 3D semiellipsoid channel when the relative time is set to zero, as demonstrated in [40].

4. Numerical Result Analysis

4.1. Simulation Scenario Description and Parameter Adjustment

The numerical findings obtained by assessing the STF-CF for the LoS and the NLoS propagation environments are presented in this section. Simulations are used to verify the accuracy of the proposed analytical results. To determine the optimized performance of the V2V MIMO channel, the STF-CF under different separating distances of transceivers, different moving speed, and different -factor values are compared with similar dimensions of an indoor tunnel between Qihua road to University station in Shanghai, China. We have considered the MIMO system in a rectangular-shaped tunnel with a height of  m, a width of  m, and a height of  m using our proposed geometrical propagation model. In our simulation results, we have considered the SNR as 20 dB. Basically, the SNR is influenced by the path loss, so we have fixed it to check the influence of multipath components effect. The real-time inner view of indoor Shanghai tunnel and geometrical representation is given in Figure 5. The tunnel is divided into two sections, the first of which is a 20 m long tunnel platform or entrance with a 6.5 m high rectangular cross-section. A nonplatform structure with a 5.96 m high arch cross-section makes up the second section. Throughout this 200 m distance, 101 different measuring locations are taken into account, each of which is 2 m apart. The following parameters are considered to produce numerical results: , , , , , and . The is fixed to the tunnel sidewalls at the entrance, whereas the moves along the tunnel axis. The tunnel walls are reinforced with concrete material. The propagation angles , , , and are estimated according to (27) and (28). The first 50 m distance in tunnel is conceived as LoS path, where the rest of 150 m is considered as NLoS path. Table 1 also lists the additional parameters that are used to generate our simulation findings. The Monte Carlo approach [41] is used to simulate the STF-CF under various -factors, time separations , and Doppler shifts . We can use the Monte Carlo approach to multiply the iteration of the TVTF’s NLoS components since the LoS components are deterministic while the NLoS components are based on a random process. The number of iterations for the Monte Carlo technique is fixed to in this paper. In our simulations, we have considered the SNR as 20 dB. Basically, the SNR is influenced by the path loss, so we have fixed it to check the influence of multipath component effect. In this paper, our aim is to check the channel correlation, where the channel correlation is influenced by the SNR and channel matrix simultaneously. In our simulation results, we pay more attention to the influence of wireless channel matrix on correlation in different environments and different - distances, which requires removing the influence of path loss and signal-to-noise ratio so that channel correlation can be compared in a fair way. After fixing the SNR, channel matrix can express the multipath effects in a proper way. In our simulations, the bandwidth is 160 MHz, which means the sampling interval of two adjacent multipath components can be as low as 6.25 ns (1/160 MHz). Because the coherent bandwidth of the channel is lower than the long-term-evolution for metro (LTE-M) signal bandwidth, we used wideband channel assumptions for the calibrations. In this paper, we have considered the reflections from the tunnel sidewalls and ceiling only for simulation results.

4.2. Angular Distribution and Propagation Delays

In Figure 6, we have shown power angular distribution by considering the elevation angle of departure (EAoD), azimuth angle of departure (AAoD), elevation angle of arrival (EAoA), and azimuth angle of arrival (AAoA) at 100 m separation distances between the transceiver antenna arrays ( m). Here, the number of scatterers is conceived as 2000 [42]. We can see that EAoA and AAoA both have negative magnitude which is because of and , respectively. The power and angular spread is configured by using the Bartlett method [37]. From Figure 6, we can observe that the propagative angles reproduce two main beams at of elevation and azimuth of and for EoA and AoA (see Figure 6(a)) and at elevation and azimuths of and for EAoD and AAoD (see Figure 6(b)), respectively. Moreover, the propagation delays are shown in Figure 7 by using RSBST approach at 5.6 GHz operating frequency. The results represent the cumulative distribution function (CDF) of propagation delays in the presence of different tunnel widths ( m, 7 m, and 9 m). It can be seen that propagation delays are increased by increasing the tunnel width. A modal theory can explain this pretty well [43]. Multiple propagation modes are triggered at the shorter to distances. At a greater distance, particularly in the case of NLoS propagation, high-order model attenuation becomes critical, and only the low-order primary model is retained. As a result, the reflected environment’s delay spread is relatively steady. The tunnel has a smaller path loss exponent in the given frequency band, but it still displays the waveguide effect. Furthermore, the signal’s shorter wavelength increases the specular reflections, which boosts the waveguide effect in the tunnel for higher frequency transmissions. When the multipath components are obstructed by the tunnel sidewalls and ceiling, the waveguide effect is considerably decreased. Also, due to the enormous 1500 m radius of curvature , there are only a few path loss exponents in the NLoS region.

4.3. PDF and Correlation Function

In Figure 8, the theoretical vs. the simulation analysis of joint PDF of proposed RSBST model for different tunnel widths () is shown by using (24). Here, antenna array spacing is set as . It can be shown that joint PDF grows as the separation distance between the and the reduces and decreases as the width () of the tunnel increases. Furthermore, the simulated results are in good agreement with the theoretical values of , verifying the validity of the proposed model. In a multipath channel, the initial geometric path length has a significant impact on the PDF of the transceiver separation distance. It is worth noting that the and are installed in two semiellipsoid scenarios in [44], and the PDF statistics are compatible with our proposed 3D model. Due to the symmetry of the geometric radio channel model, each joint PDF curve in Figure 8 is symmetric, which is also compatible with [44]. Thus, the proposed 3D MIMO model can be used to characterize the real tunnel environment.

Figure 9 presents the STT-CF for different -factor values, when  Hz, coherence time  ns, and separating distance between transceiver arrays is  m. As shown in Figure 9, the STF-CF decreases as the tunnel width or and antenna array spacing increases, and the STF-CF increases as the -factor increases. Meanwhile, the STF-CF simulation values are in good agreement with the theoretical values, indicating that the proposed model is suitable. The LoS components have a negative impact on MIMO channel performance. Figure 10 represents the STF-CF for different normalized time lag and Doppler shift for when , and distance between and is set to  m. The STF-CF of the simulation results fits very well with the theoretical results, as illustrated in Figure 10. It is worth noting that when temporal separation increases, the STF-CF decreases. In our proposed model, we have considered the LoS and diffused components which provide different STF-CF values instead of 1 under different antenna spacing and normalized time lag, as given in [45]. When the Doppler shift raises, the STF-CF reduces rapidly. The moving velocity of vehicles affects the V2V channel statistics of the MIMO system, as shown in [46], and this behavior can also be shown in Figure 10.

5. Conclusions

In this paper, we have presented a novel 3D scattering model for V2V MIMO radio communications in geometrical rectangle tunnel environment. Moreover, the statistical characteristics and analytical results of the proposed model allow one to understand the real tunnel propagation environment for 5G-enabled IoT communications, as it assumes the randomly distributed scatterers along the tunnel sidewalls and ceiling. A generic analytical expression for the STF-CF has been obtained by considering the nonisotropic scattering conditions in both LoS and NLoS propagation environments. Near- and far-field effects are taken into account in this RSBST model to characterize the spherical wavefronts, resulting in AAoA, EAoA, EAoD, and AAoD statistics and PDF and STF-CF responses on the MIMO antenna arrays. It is proved that the nonstationarity, including moving direction and propagation delays, affects the MIMO channel performance. From the analytical results, we can deduce that the STT-CF is seriously affected by propagation delay, -factor, movement velocity of receiving antenna array, and antenna array spacing. The proposed model results are aligned very well with previous propagation models, indicating that the proposed massive MIMO channel model can be generalized.

Data Availability

The data that support the research findings are available on request. The data are not publicly available due to the privacy of research participants.

Conflicts of Interest

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