Abstract

With the rapid development of high-mobility wireless communication systems, e.g., high-speed train (HST) and metro wireless communication systems, more and more attention has been paid to the wireless communication technology in tunnel-like scenarios. In this paper, we propose a three-dimensional (3D) nonstationary multiple-input multiple-output (MIMO) channel model with high-mobility wireless communication systems using leaky coaxial cable (LCX) inside a rectangular tunnel over the 1.8 GHz band. Taking into account single-bounce scattering under line-of-sight (LoS) and non-line-of-sight (NLoS) propagations condition, the analytical expressions of the channel impulse response (CIR) and temporal correlation function (T-CF) are derived. In the proposed channel model, it is assumed that a large number of scatterers are randomly distributed on the sidewall of the tunnel and the roof of the tunnel. We analyze the impact of various model parameters, including LCX spacing, time separation, movement velocity of Rx, and K-factor, on the T-CF of the MIMO channel model. For HST, the results of some further studies on the maximum speed of 360 km/h are given. By comparing the T-CF between the dipole MIMO system and the LCX-MIMO system, we can see that the performance of the LCX-MIMO system is better than that of the dipole MIMO system.

1. Introduction

Recently, high-mobility wireless communication systems have developed rapidly, for example, high-speed train (HST) and metro wireless communication systems. The speed of HST is getting faster and faster; for example, the trail maximum speed of HST reaches nearly 575 km/h by the French National Rail Corporation, the maximum speed of Japan’s seven-car “magnetic levitation train” reaches 603 km/h, China’s CIT500 HST hits a train wheel of 605 km/h, and the current operating speed of China’s HST is 350 km/h. Leaky coaxial cable (LCX) is widely used in urban subway tunnel environments. LCX has many advantages; for example, its radiation characteristics provide uniform coverage, longitudinal attenuation along the axis of LCX is very small, and its installation is very simple [112]. Therefore, LCX is very suitable for the wireless communication system in a tunnel scenario. In addition, the multiple-input multiple-output (MIMO) technology is widely used in improving the spectrum efficiency and the propagation characteristics of wireless communication systems. In order to improve the stability and reliability of wireless communication systems in the tunnel environment, the MIMO using LCX technology has become an inevitable trend.

The previous research on LCX mainly includes the following: the characteristics analysis of LCX, for example, radiation characteristic [3]; the performance analysis of 2 × 2 MIMO system using single LCX in free space, e.g., channel capacity and correlation [1]; the study of MIMO system performance using LCX in tunnel scenario [2]; the performance of LCX-MIMO in an indoor office environment, e.g., channel capacity and channel power spread [4]; the statistics characterization of 2 × 2 MIMO system using single LCX in linear-cell environment, including capacity using different power allocation (i.e., equal power allocation and water-filling power allocation) over the 2.4 GHz band and the 5 GHz band [5], condition number (CN) with different polarization types (i.e., vertical polarization and horizontal polarization) over wireless local area network (WLAN) frequency band (i.e., 2.4 GHz) [6, 7]; and the performance of 4 × 4 MIMO system using two LCX in linear-cell environment, including capacity using different power allocation (i.e., equal power allocation and water-filling power allocation) with various LCX spacing over the 2.4 GHz band and the 5 GHz band [8], condition number (CN) with different LCX spacing over the 5 GHz band [9], and the 2.4 GHz band [10]. In [11], the author proposed a LCX-MIMO channel model based on electric field distribution of LCX in tunnel, and capacity, CN, and correlation are analyzed by simulation and measurement. In [12], the author proposed a nonstationary LCX-MIMO channel model based on electric field distribution of LCX and geometry-based stochastic model (GBSM) in the tunnel. The space-time correlation is analyzed by simulation. The velocities of movement are 90 km/h, 180 km/h, and 360 km/h. According to the theory of electromagnetic field, with the increase of electromagnetic wave propagation distance, the electric field strength decreases gradually. The channel model in [12] is mainly based on the electric field distribution around the LCX, and only considering the change of the amplitude of the electromagnetic wave signal with the increase of the propagation path distance, but not considering that when the electromagnetic wave signal is scattered by the scatterer, it will not only introduce the phase shift by the scatterer but also attenuate the amplitude of the signal by the scatterer.

With the rapid development of wireless communication, more and more attention has been paid to the nonstationary MIMO fading channel in various scenarios. According to various scenarios, there are different channel models, including spatially and temporally correlated channel model for MIMO system in mobile fading channels [1315] and mobile-to-mobile channels [16], 3D nonstationary MIMO channel model with three-dimensional (3D) antennas array and 3D arbitrary trajectories with mobile station [17], 3D nonstationary GBSM unmanned aerial vehicle (UAV) MIMO channel model for air-to-ground (A2G) [18] and wideband [19] wireless communication system, 3D nonstationary geometry-based MIMO channel model for vehicle-to-vehicle (V2V) communication in the tunnel (i.e., semicircular tunnel and rectangular tunnel) environment [2022], and geometry-based single-bounce (GBSB) channel model for MIMO system in tunnel [23]. In these channel models, the influence of different channel parameters, including the transmitter and receiver antennas spacing, the number of transmitter and receiver antennas, angle of arrival (AOA) and angle of departure (AOD), movement time of transmitter/receiver antennas, movement velocity of transmitter/receiver antennas, and time separation, on channel capacity and correlation is analyzed.

In summary, so much literature has analyzed the performance of the LCX-MIMO system. At the same time, much literature has proposed a large number of nonstationary channel models for various scenarios. However, there is a few of high-mobility channel model for LCX-MIMO wireless communication system in a tunnel environment.

In order to further study the performance of a high-mobility wireless communication system using LCX-MIMO in a tunnel-like scenario, this paper proposes a 3D nonstationary LCX-MIMO channel model with a high-mobility wireless communication system in a rectangular tunnel based on GBSM, and we are considering that when the electromagnetic wave signal is scattered by the scatterer, it will not only introduce the phase shift by the scatterer but also attenuate the amplitude of the signal by the scatterer. The analytical expressions of the channel impulse response (CIR) and the temporal correlation function (T-CF) are derived. The utility and reasonability of the proposed channel model are verified by experiment and simulation. Meanwhile, we analyze the influence of various model parameters, including LCX spacing, time separation, movement velocity of Rx and K-factor, on the T-CF of the MIMO channel model by simulation. For HST, the results of some further studies on the maximum speed of 360 km/h are given.

The remainder of this paper is organized as follows: In Section 2, the LCX-MIMO channel model is derived. In Section 3, the channel T-CF is derived theoretically, and the formula of capacity and K-factor are described. Section 4 is experiment results and analysis. Section 5 is simulation results and validation. Section 6 describes the conclusion of this paper.

2. LCX-MIMO Channel Model

In this paper, for the rectangular tunnel scenario, the transmitters (Tx) and receivers (Rx) of the LCX-MIMO communication system are consisting of two LCXs and two dipole antennas, respectively. In Figure 1, o represents the origin of the coordinate. The parallel direction to the axis of the tunnel is considered as the z-axis of the coordinate system. The cross section of the tunnel is in the xoy plane, the vertical direction of the plane is along the y-axis of the coordinate system, and, finally, the horizontal direction of the plane is along the x-axis of the coordinate system; the S1 and S2 are input signals of LCXs. Figure 1(a) describes the structure of LCX; Figure 1(b) describes the propagation scenario of MIMO using two LCXs in the tunnel; Figure 1(c) describes the nonstationary 3D LCX-MIMO channel model in the tunnel, (−∙∙−) denotes single-bounce non-line-of-sight (NLoS) path, (−∙−) denotes line-of-sight (LoS) paths; and Figure 1(d) describes the randomly distributed scatters in the tunnel.

Here, the slot period of LCX is set as P, L is the length of LCX, and the number of slots is N  = floor (L/P)−1. According to the properties of LCX, the longitudinal attenuation factor of the i-th (i = 1 ,…, N) slot of LCX is , where the longitudinal amplitude attenuation of i-th slot is αi = , α0 is the longitudinal amplitude attenuation per unit period of slots, and βi is the phase variation of i-th slot. The LCX used in this paper refers to the LCX with a working frequency band of 1.8 GHz provided by the Zhongtian technology company (ZTT), and the longitudinal power attenuation constant of LCX is α = 4.0 dB/100 m. It can be seen that α0 = 10αP/10/2/100 ≈ 1, so the influence of αi on CIR is ignored in this paper. The phase variation of i-th slot of LCX and the propagation constant in LCX are βi and kr, respectively, and it can be expressed as follows:where β0 = krP is the phase variation per unit slot in LCX; εr is the dielectric constant in LCX; k0 = 2πf/c is the propagation constant in free space; f is the center frequency of input signal of LCX; c is the velocity of light; and λ is the wavelength.

For Figure 1, the definition of LCX-MIMO model parameters is shown in Table 1. Because the surfaces of the sidewalls of the tunnel and the top of the tunnel are rough, we assume that there are a large number of scatterers at the sidewalls and the roof of the tunnel, and the number of effective and independent scatterers is M. So, there are both line-of-sight (LoS) and non-line-of-sight (NLoS) propagation paths in the propagation channel. The definition of LCX-MIMO model parameters is described in Table 1. The propagation characteristics of the LCX-MIMO system can be described by 2 × 2 CIR matrix; that is, H(t) = [hqp(t)]2×2. Therefore, the CIR hqp(t) from the p-th LCX to the q-th receiving antenna is formed by vector superposition of components (LoS and NLoS), and it can be expressed as follows:

In the expression of CIR, Ωqp represents the transmission power from the p-th transmitting antenna to q-th receiving antenna, that is, Ωlp = E[|hlp(t)|2]; Kqp is Rice factor, that is, ; and is the amplitude of the wave scattered by m-th scatterer Sm to receiving antenna so that as N ⟶ ∞ and M ⟶ ∞. In this paper, it is assumed that is a zero-mean Rayleigh random distribution variable. φim denotes the independent and identically distributed (i.i.d.) random phase shift introduced by m-th scatterer Sm, which is the uniform distribution over [0, 2π). In order to guarantee that Ωqp = E[|hqp (t)|2], the CIR in the above formula needs to be satisfied as N ⟶ ∞ and M ⟶ ∞, where fD denotes the maximum Doppler shift, that is, fD = /λ. The azimuth angle of velocity is . In this paper, we assumed that the receiving antenna moves forward along the z-axis. Therefore, according to Figure 1(c),  = 0.

3. Channel T-CF and K-Factor

In this paper, it is assumed that LoS and NLoS propagation paths are independent of each other, and the mean amplitude of CIR is zero. According to [13, 2023], we consider the transmitter correlation and receiver correlation are independent of each other in this paper. Because the receiver correlation is influenced by the receive antennas spacing, we focus on the correlation between two LCXs in the Tx and the first receive antenna. In his paper, we set the independence between LoS path and NLoS path. Then, the channel T-CF between LCXs (i.e., T-CF between CIR h11(t) and CIR h12(t)) can be expressed as follows:where (•) is the complex conjugate, τ is time separation, is the LoS component of T-CF between LCXs, and is the NLoS component of T-CF between LCXs.

3.1. LoS Channel T-CF

According to formula (2a), the LoS component of T-CF between LCXs can be expressed as follows:

In order to simplify the T-CF of LoS, based on the geometric relationship described in Figure 1, it can be seen that . In this paper, we consider di,11 >> d and di,12 >> d, so  ≈ . Then, the LoS T-CF can be expressed as follows:where di,11 is the distance between i-th slot of LCX1 and Rx1, di,12 is the distance between i-th slot of LCX2 and Rx1, and they are given as

3.2. NLoS Channel T-CF

According to (2b), the NLoS component of T-CF between LCXs can be expressed as follows:where dim,Tx1 is the distance between i-th slot of LCX1 (i.e., Tx1) and m-th scatterer Sm, dim,Tx2 is the distance between i-th slot of LCX2 (i.e., Tx2) and m-th scatterer Sm, dim,Rx1 is the distance between m-th scatterer Sm and Rx1, and they are given as

In order to simplify the T-CF of LoS, based on the geometric relationship described in Figure 1, it can be seen that . Then, the NLoS T-CF can be expressed as follows:

3.3. K-Factor

The Rice K-factor is usually used in evaluating the communication link quality, such as capacity and performance of the MIMO system. Therefore, the estimation of the K-factor is very important in different wireless scenarios. With fixed SNR level, the higher K-factor means larger spatial correlation and smaller channel capacity. Usually, we can estimate K-factor based on the measurement data for the MIMO channel. For a complex channel gain matrix H (t), the K-factor can be estimated as [24]where |•| denotes the operation of matrix determinant, ε denotes the mean value of |H (t)|, var denotes the variance of |H(t)|.

4. Experiment Results and Analysis

4.1. Parameter Setup

In this paper, an experiment was carried out in a tunnel-like scenario at Zhongtian Technology Company (ZTT), Nantong, China. In Figure 2, the overall length of the tunnel is 100m, and the rectangular part and arched part of the tunnel are 50 m, respectively. In this paper, the rectangular part has been considered in the experiment. As shown in Figure 2, it is the experiment scenario description for Nantong tunnel; Figure 2(a) is a photograph of the tunnel; Figure 2(b) is the front view of the tunnel; Figure 2(c) is experiment configuration of LCX-MIMO system; and Figure 2(d) is the top view of the tunnel. And the description of other parameters in the experiment is shown in Table 2.

4.2. Results and Analysis

In this section, we estimated the K-factor using formula (10) to guarantee that the K-factor of the theory is consistent with the experiment. Figure 3 describes the cumulative distribution function (CDF) of the K-factor in the experiment, and we can see that the average value of the K-factor is 4.032 dB. Therefore, we set the theoretical value of the K-factor as 4 dB. The statistics of the K-factor in measurement is shown in Table 3.

5. Simulation Results and Validation

In this section, firstly, we compare the theoretical and simulation results of T-CF to further verify the reasonability of the channel model proposed in this paper. Secondly, we will show that the performance of the LCX-MIMO system is better by comparing the T-CF of the dipole MIMO system with the T-CF of the LCX-MIMO system. According to (3)–(9), the Monte–Carlo method is used in simulating the T-CF (LoS + NLoS) of different time separations τ, different K-factors, and different maximum Doppler shifts fD (i.e., fD = /λ). Because the LoS component of CIR is deterministic, and NLoS component of CIR is a random process, we can multiply the iteration NLoS component of CIR by the Monte–Carlo method. In this paper, we set the number of iteration as 2 × 104.

5.1. Parameter Setup

Here, we assume that the scatterers are randomly distributed on the sidewall of the tunnel and the roof of the tunnel as shown in Figure 1(d). The simulation scenario of LCX-MIMO is the same as Figure 1. The simulation propagation scheme of the MIMO system is shown in Figure 4. For the LCX-MIMO system, Tx is two LCXs and Rx is two dipole antennas. For the dipole MIMO system, the Tx and Rx are both two dipole antennas. These two systems have the same parameters configuration at the same scenario. For the LCX-MIMO system and dipole MIMO system, the Tx antennas have the same x coordinate and y coordinate, but the z coordinate of Tx antennas for dipole MIMO is 0 m. The Rx array moves along the z-axis. The number of tunnel wall scatters and tunnel roof scatters is assumed M = 500. The other parameters are the same as Section 4.

5.2. Results and Analysis

Here, based on (3), (5), and (7), the numerical results are obtained by simulation for the T-CF of LCX-MIMO channel model in rectangular tunnel scenario.

Figure 5 presents the T-CF for various LCX spacing d with different K-factor (here, we assume K11 = K12) when time separation τ = 1 ns, fD = 600 Hz (i.e.,  = 360 km/h), l = 0.5 m. As we can see in Figure 5, the T-CF decreases with the LCX spacing increasing, and the T-CF is higher when K-factor is larger. Meanwhile, the theoretical value of T-CF matches very well with the simulation value of T-CF. It shows that the model proposed in this paper is reasonable. In the MIMO system, the LoS path has a negative effect on the channel capacity. The K-factor is the ratio between the power of LoS path and the power of NLoS path. Therefore, the K-factor is larger when the power of LoS path is higher. Finally, the T-CF is higher when K-factor is larger.

Figure 6 shows the T-CF for different time separation τ with various Doppler shifts fD [150 Hz (90 km/h), 300 Hz (180 km/h), 600 Hz (360 km/h)] of Rx when the LCX spacing d = 4.8λ, Rice factor K11 = K12 = 4 dB, and element spacing of Rx array l = 0.5 m. As shown in Figure 6, the T-CF of theory matches with the T-CF of simulation very well. It is obvious that the T-CF decreases with increasing time separation τ. The T-CF decreases more rapidly when Doppler shift fD increases. As shown in [16, 18], their results show that the movement velocity of the moving vehicles will affect the V2V channel characteristics of the MIMO system.

In Figure 7, it shows the T-CF of dipole MIMO system and the T-CF of LCX-MIMO system for different Doppler shifts fD [150 Hz (90 km/h), 300 Hz (180 km/h), 600 Hz (360 km/h)] of Rx when the LCX spacing d = 4.8λ, Rice factor K11 = K12 = 4 dB, and element spacing of Rx array l = 0.5 m. By comparing the T-CF between dipole MIMO system and LCX-MIMO system, we can see that the T-CF of LCX-MIMO decreases more rapidly than T-CF of dipole MIMO with different time separation τ when Doppler shift fD increases. It means that the performance of LCX-MIMO is better than the performance of dipole MIMO in the tunnel.

6. Conclusion

In this paper, we have proposed a 3D nonstationary LCX-MIMO channel model with a high-mobility wireless communication system using LCX in a rectangular tunnel environment. The effective scatterers of the tunnel are considered to be randomly distributed on the roof and the sidewalls of the tunnel in this model. Taking into account single-bounce scattering under LoS and NLoS propagation condition, we derived the analytical expressions of CIR and T-CF. In order to verify the channel model proposed in this paper, the simulation is carried out in a tunnel-like scenario.

From the numerical results, we conclude that the T-CFs is significantly affected by time separation, K-factor, movement velocity of Rx, and LCX spacing. Finally, by comparing the T-CF between dipole MIMO system and LCX-MIMO system, we can see that the T-CF of LCX-MIMO decreases more rapidly than that of dipole MIMO with different time separation τ when Doppler shift fD increases, which proves that the performance of LCX-MIMO system is better than that of the dipole MIMO system.

Data Availability

The measured data used to support the findings of this study have not been made available because the research findings will be commercialized. The simulation data is obtained by the special software purchased by our laboratory. This paper uses MATLAB software to process the data results.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Foundation of China under Grant 61871261 and Project of Shanghai Science and Technology Commission 18511105704. The author particularly appreciates Mr. Asad Saleem for his help during checks and modifications of this paper.