Research Article  Open Access
Juan Moreno GarcíaLoygorri, César BrisoRodríguez, Israel Arnedo, César Calvo, Miguel A. G. Laso, Danping He, Florentino Jiménez, Vicente González Posadas, "Wideband Channel Modeling for mmWave inside Trains for 5GRelated Applications", Wireless Communications and Mobile Computing, vol. 2018, Article ID 6916954, 11 pages, 2018. https://doi.org/10.1155/2018/6916954
Wideband Channel Modeling for mmWave inside Trains for 5GRelated Applications
Abstract
Passenger trains and especially metro trains have been identified as one of the key scenarios for 5G deployments. The wireless channel inside a train car is reported in the frequency range between 26.5 GHz and 40 GHz. These bands have received a lot of interest for highdensity scenarios with a hightraffic demand, two of the most relevant aspects of a 5G network. In this paper we provide a full description of the wideband channel estimating PowerDelay Profiles (PDP), SalehValenzuela model parameters, timeofarrival (TOA) ranging, and pathloss results. Moreover, the performance of an automatic clustering algorithm is evaluated. The results show a remarkable degree of coherence and general conclusions are obtained.
1. Introduction
Metro lines are one of the most acknowledged highdensity scenarios, as well as stadiums and other large venues. These places are a clear target for 5G mobile [1] deployments, which is intended to provide very large throughputs and ultralow latencies in dense environments and sustain a high number of simultaneous connections, among others. It is noteworthy that it is possible to have more than 1200 people inside a 120meter long train, which means around 5 persons/m^{2}. There are many railwayrelated applications like TraintoGround, TraintoTrain, TraintoCar, TraintoSatellite, intratrain communications, and so on, which are usually grouped under the T2X notation, or perhaps in the more generic V2X notation, where “V” stands for vehicular. In this paper we focus on the intratrain channel, which has an enormous importance to provide connectivity to passengers. Also, this wireless channel could be used for the internal communications of the train, which is one of the objectives of the EUH2020 Roll2Rail Project [2].
There is still some vagueness in the frequency bands to be assigned, but some things are clear: there will be 5G deployments both above and below 6 GHz [3]. This fact has attracted much attention on millimeter waves (mmwave) in order to have larger bandwidths than in lower frequencies, where the entire spectrum has been allocated. For example, in USA both the 28 GHz (27.5–28.35 GHz) and the 37 GHz band (37–40 GHz) have been licensed for mobile applications; in both China and the European Union, the 26 GHz band (24.25–27.5 GHz) is the chosen one. On the other side, there is no clear direction on which bands will be allocated for unlicensed use, but the 37 GHz band (37–37.6 GHz) is perhaps the most promising one. All these spectrum policies lead to the need of a proper channel model for 5Grelated scenarios in the aforementioned bands.
There are many papers related to channel measurements and modeling in railway scenarios but, as far as we know, only very few of them are focused on mmwave for intratrain environments (i.e., the work carried out by ETRI in Korea [4]). Regarding the TraintoGround link for highspeed trains, a very complete reference is [5] but it is focused on frequencies below 6 GHz; for tunnels which are a very common scenario in railways, there are many papers, like [6, 7]; inside metro stations there are many research papers like [8]; and for mmwave research in railways there are some raytracing simulations in mmwave and THz [9]. For intratrain link, there are some recent papers but below 6 GHz [10, 11].
The structure of the paper is the following: in Section 2 we will describe the experiment setup, covering the environment details, the measurement setup, and the processing of the data; in Section 3 we will explain all the results obtained from the measurements, and, finally, in Section 4 a conclusion and some discussion are provided.
2. Experimental Setup
2.1. Environment
The environment for the measurement campaign is the train interior of a metro train. It is depicted in Figures 1 and 2. This vehicle under test is a Madrid Metro train of the s3000 series. It has the capacity to carry 442 persons in the fourcar consist, which amounts to 132 persons on each car.
(a)
(b)
The maximum dimensions of the train are meters, but the interior of the train has a maximum section of 2.2 m 2.2 m. As usual in this type of trains, there are many handholds, seats, and other furniture. The seats are resincoated reinforced with glass fiber and both the handholds and the roof are made of stainless steel covered with a 500 micrometers thick yellow painting. The floor is made of stratified rubber, the windows are made of laminated glass and the doors are of aluminum and glass like the windows’ one. This layout and materials are very common in modern subway trains, with very small differences between them.
2.2. Measurement Setup
The measurements were carried out using an Agilent 8722ES Vector Network Analyzer (VNA). The approach for this VNAbased measurement is the usual one: connecting the transmitter to port 1 and the receiver to port 2, measure the s21 parameter, which is a good approximation for the channel transfer function [12]. The connection of both antennas to the VNA ports was done using phase stable RTK040 wires. The length of this wire was 18 m and we used 2.4 mm connectors.
All the measurements but one were carried out using horn antennas (see Figure 3), model 22240 by FLANN, which is an antenna designed to be used in the 26.5–40 GHz, which is the band where our experiment was focused on. The gain of this antenna is 20 dB, and the variations of this parameter are bounded to 2 dB in the band of interest, which gives us a great stability. In the last measurement an omnidirectional antenna (Vivaldi, see Figure 3) was used in the receiver as a replacement of the horn antennas mentioned before.
(a)
(b)
The measured positions for the receiver (the transmitter was fixed) are depicted in Figure 2. We measured at five different transmitterreceiver distances = , 4, 6, 9.5, m, moving only the receiver (in the last one the receiver is in the next train car). The height of both the transmitting and the receiving antennas was 1 m over the train floor. Both the transmitter and the receiver were placed inside the train at all time, as well as the VNA. In the measurement setup the number of dots in the 13.5 GHz wideband (26.5–40 GHz) was 1601. The processing of the data was done offline with MATLAB. As we explain in a more detailed form later in this paper, we decided to employ an automatic clustering algorithm instead of identifying them by visual inspection.
3. Channel Model
3.1. PathLoss
The pathloss is obtained from the measured data in the VNA. This result is the average in the entire frequency range (26.5–40 GHz) [13]. equals 1601 which is the number of dots
The obtained values are detailed in Table 1.

The pathloss increases together with the distance, both the average pathloss and the LOS (θ = 0°) scenario, as it was expected. Regarding the influence of the angle (θ), for θ = 0° the received power is the maximum one for every measured distance, and we can see a great increase in the pathloss when we misalign the antennas = 45°, 90°, 135°, . This increase is more significant at lower distances. For example, the average pathloss increases 31.46 dB at = 2 m and only 8.84 dB at = 17 m. This is because at higher distances the multipath is less severe and the impact of reflections on handlers and other diffuse scatters is significantly lower (the socalled Waveguide Effect).
The variation of the pathloss for the different angles (given a fixed distance) is more or less random. However, at the highest distances = 9.5, 17 the angle that experiences a lower pathloss is clearly θ = 45°. The reason for this is the same as before: the Waveguide Effect. To properly assess the Waveguide Effect, we have calculated the pathloss exponent () for 0°. The obtained value is 1.6, which is below the freespace one () so, under this premises, we can conclude that we are experiencing this effect inside the train at this frequency band. The whole oneslope model is shown in where is the total losses in dB, is the reference distance, is the losses at the reference distance, is the pathloss exponent, and is the distance.
We also measured the pathloss at a given distance ( = 9.5 m) with six people moving around between the transmitter and the receiver. This impact is an average pathloss 2.28 dB higher and some fading 30 dB deep.
3.2. Ranging
If we measure the time of arrival (TOA) we can easily estimate the distance between the transmitter and the receiverwhere is the time for the first peak in the PDP (i.e., the first path that arrives to the receiver) and is the speed of light. In Table 2 we can see the estimated distances using this TOAbased technique to estimate the distance between transmitter and receiver. The estimation based on the LOS path is the most accurate (5.56% error in the worst case and 0.35% in the best one); for θ = 45° it is slightly worse and for the other angles we have large errors. This method seems to be more accurate when the distance increases (for = 17 m, the error is lower than 3% for θ = 0°, 45°, and 90° and 13.20% for θ = 180°). The reason behind this good performance at higher distances is because the impact of reflections from handlers and other furniture is less significant than at lower distances.

3.3. PowerDelayProfiles and Clustering
Regarding the PowerDelay Profiles (PDP) we have computed them for all the scenarios described in Section 2. In Figures 5–9 we can see the PDPs for = , 4, 6, 9.5, m and an extra one for the omnidirectional antenna at a distance of = 4 m.
Multipath components are clustered according to SalehValenzuela [14] model and the estimation of the clusters that these multipath components form has been done using the automatic algorithm proposed in [15]. This algorithm takes into account both the amplitude and time delay to form the clusters. The usual approach to identify multipath component clusters is the visual inspection [15] but we found preferable to use an automatic algorithm. The main advantage of this algorithm is that it allows us to do a fast estimation of the clusters, and that it takes into account both time and amplitude values. Therefore, the clustering process is very dependent on the thresholds we propose (as it is acknowledged in [15]), but after some iterations, it seems to have a good performance if we look to Figures 5–11. In these figures we have highlighted in red the clusters that this algorithm has identified (dots in red are the first ray from a cluster).
For = 2 m (Figure 5) we can see that the number of clusters increases with θ. For θ = 90° the number of clusters increases significantly, and for θ = 180° it decreases. This is an expectable result, because as we go towards a pure NLOS scenario, with many scatters between transmitter and receiver, this leads to different times of arrival (i.e., clusters).
Regarding the = 4 m scenario (Figure 6) the number of clusters does not increase significantly from the previous scenario, but we have more or less the same dependence to θ. The difference between θ = 0° and θ = 45° is small if we visually inspect both PDPs, but for θ = 180° we have more clusters in the receiver. In our opinion, this is attributable to the layout of the train that makes this PDP possible for this angle. If we employ an omnidirectional antenna in the receiver (Figure 11) we can see that we are able to see a significantly higher number of clusters. This is an expected result because the wider radiation pattern allows more different directions of arrival to the receiver. We still have a strong first cluster, associated with the LOS direct link plus four more clusters for the reflection through the train car.
In Figure 7 we can see the PDPs for = 6 m. The number of clusters remain stable and similar to = , 4 for θ = 0°. For θ = 45° we see more clusters for the same reason as before: losing the LOS link leads to more reflections due to the presence of obstacles inside the train. This is coherent with the fact that for θ = , the number of clusters does not increase.
For = 9.5 m (Figure 8), the estimated number of clusters is stable (and lower than before) for both θ = , because we are experiencing a significant Waveguide Effect and there are many rays that are suppressed. For θ = , we have the same increase in the number of clusters experienced before. So the pattern does not change from the = 6 m scenario. Nevertheless, for this distance we also measured the channel with people moving randomly inside the train (between transmitter and receiver). The result can be seen in Figure 10, and the number of clusters increases significantly from any other angle at the same distance. The reason is the expected one: people act as supplementary scatters, so the PDP changes (we have more clusters, obviously). This result is very important because intratrainbased applications will work on trains with people onboard.
Finally, for = 17 m (Figure 9), if we inspect visually the PDPs for the four angles considered, we can see that it does not change significantly, but the algorithm estimates 2, 6, 3, and 4, clusters for θ = , 45°, 90°, . This is perhaps the most debatable outcome of the algorithm. In this case, the distance is significantly higher than in the other scenarios, where the Waveguide Effect is dominant and the PDP is smoother (which means fewer clusters). It is also the scenario where we have less differences between angles (together with = 2 m, where there were almost no scatters between transmitter and receiver). An expected result is that, at a distance like this, the power is mostly concentrated on the first cluster of rays that arrives to the receiver (this is easy to see in the θ = 0°, obviously). This is another consequence of the “Waveguide Effect.”
3.4. SalehValenzuela Parameters
The classical SalehValenzuela (SV) model is intended for indoor scenarios and we have chosen it for its statistical nature and simplicity: it describes the channel with only four parameters (λ, Λ, γ, and Γ). Λ is the cluster arrival rate, which is assumed to follow a Poisson distribution [14], the same distribution as the ray arrival rate “λ”; “γ” is the exponential decay within a cluster; and “Γ” is the exponential decay for the whole PDP. In Table 3 we have included our estimation of these four SV parameters for all the scenarios that we have measured.

If we look at the average SV parameters in the column on the right of Table 3, we can see that we have a “slow” decay. This is because the distance between the objects that reflect the signal is very low, and all of them are very similar (metallic surfaces). That is for the “overall decay” (parameter Γ). If we look inside the clusters, we can see that they decay rapidly (parameter γ), which means a higher attenuation (due to the high frequencies) and very few reflections (apart from the first one). Anyway, this is coherent because the clustering algorithm that we chose is adequate, because the calculation of SV model parameters is very dependent on the clustering algorithm.
Explaining properly all the estimated SV parameters would take a lot of time, so here we will discuss some of the most significant parameters. First of all, it is noteworthy that λ (ray arrival rate) has small variations no matter in which scenario we are. This parameter is very dependent on the limitations of the measurement devices, so perhaps with a different device able to increase the number of resolvable rays this parameter could experience higher variations. The main difficulty to improve this is that the VNA employed is a stateoftheart device, so it will take some time to achieve a better resolution in this parameter.
The cluster arrival rate (Λ) has a direct relation with the number of clusters and we have gone through this parameter in Section 3.3. Regarding the overall decay rate (Γ), it is generally higher than the cluster decay rate (γ), which is an expected result. Moreover, the clustering algorithm does not allow negative slopes (if it resolves a cluster with a negative slope; that cluster is joined to the next one) which increases—in general terms—the values of γ. If we look at some particular results, we can see that the Γ for the shortest distance ( = 2 m) and LOS scenario has the smallest value calculated. For higher distances and angles, the Γ increases, but the general trend is to decrease (but for = 2 m), which is an expected result. It is hard to find a clear pattern in the angle variation, but a possible reason behind this behavior is that the first ray of the cluster is not always the one with more power [15]. This could lead to an underestimation of the Γ parameter.
We also took two more measurements: with people moving around inside the train ( = 9.5 m) and also using an omnidirectional antenna ( = 4 m). The SalehValenzuela model parameters for the omnidirectional antenna are γ = 0.187; Γ = 0.208; λ = 3.800; Λ = 0.277. This is similar “γ” and “λ” for horn antenna in the receiver at the same distance. The apparent invariability of “λ” was discussed before and for “γ” it could be an impact of the clustering algorithm, which prevents higher variations inside a cluster, which leads to an approximately constant “γ.” The number of clusters is significantly different, as the parameter “Γ” highlights, so this is coherent with the “conservative” nature (regarding “γ”) of the clustering algorithm.
If we consider the difference between having people in motion between transmitter and receiver (the SV parameters are γ = 0.134; Γ = 0.188; λ = 3.745; Λ = 0.277) and not, using horn antennas and the same distance, equal to 9.5 m, we can see that we have a lower “Λ” when we have people in motion, and also a lower “Γ.” It is also noteworthy that the average parameters using horn antennas (this is θ = , 45°, 90°, ) are a good approximation of the scenario with an omnidirectional antenna. This is a coherent result because using an antenna like this is the same as we average the contributions from all directions. However, as we can see in Figure 4, the radiation pattern of the Vivaldi antenna is not perfectly omnidirectional, so this “averaging” process is not as precise as it could be otherwise.
(a)
(b)
(c)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
4. Conclusion
The SV model can be applicable directly to an intratrain scenario, as it is an indoor scenario (the scope of the SV model). It is only necessary to estimate the parameters of the model to each environment. In particular, intratrain communications at high frequency are a high multipath environment, where SV model is intended to work.
The intratrain link is one of the most likely scenarios to require high throughputs, ultradense networks, and very low latencies. These three requirements are the core of 5G technologies. In this paper we have assessed this link in a region of the spectrum related to future 5G developments (26.5–40 GHz), providing estimations for the pathloss, PowerDelay Profile and SalehValenzuela model parameters. We have seen the dependence of all these results on distance and angle; moreover, we have seen the impact of people moving around and also the differences coming from the radiation pattern of the receiving antenna (using an omnidirectional antenna instead of the horn antennas used in all measurements but one). The performance of the clustering algorithm is remarkable (with some limitations, as we have discussed) and the obtained SV parameters are coherent as well.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
Enabling 5G TEC201455735C32R is funded by the Spanish Ministry of Economy and Competitiveness and also is funded by the Chinese Strategic International Cooperative Project of National key R&D Plan, 2016YFE0200200.
References
 5G for Europe: An Action Plan, European Commission, Brussels, 2016, https://ec.europa.eu/transparency/regdoc/rep/1/2016/EN/12016588ENF11.PDF.
 Roll2Rail Project, http://roll2rail.eu.
 S.Y. Lien, S.L. Shieh, Y. Huang, B. Su, Y.L. Hsu, and H.Y. Wei, “5G New Radio: Waveform, Frame Structure, Multiple Access, and Initial Access,” IEEE Communications Magazine, vol. 55, no. 6, pp. 64–71, 2017. View at: Publisher Site  Google Scholar
 G. Noh, J. Kim, H. S. Chung, B. Hui, Y. Choi, and I. Kim, “mmWavebased mobile backhaul transceiver for high speed train communication systems,” in Proceedings of the IEEE Globecom Workshops (GC Wkshps), pp. 1–5, Singapore, Singapore, December 2017. View at: Publisher Site  Google Scholar
 T. Zhou, C. Tao, S. Salous, L. Liu, and Z. Tan, “Implementation of an LTEbased channel measurement method for highspeed railway scenarios,” IEEE Transactions on Instrumentation and Measurement, vol. 65, no. 1, pp. 25–36, 2016. View at: Publisher Site  Google Scholar
 J.M. MolinaGarcíaPardo, M. Lienard, P. Degauque, C. GarcíaPardo, and L. JuanLlácer, “MIMO Channel capacity with polarization diversity in arched tunnels,” IEEE Antennas and Wireless Propagation Letters, vol. 8, pp. 1186–1189, 2009. View at: Publisher Site  Google Scholar
 S. Hairoud, P. Combeau, Y. Pousset, Y. Cocheril, and M. Berbineau, “WINNER model for subway tunnel at 5.8 GHz,” in Proceedings of the 12th International Conference on ITS Telecommunications (ITST '12), pp. 743–747, Taiwan, November 2012. View at: Publisher Site  Google Scholar
 et. Al. Ke Guan, “Broadband Channel Measurements inside Metro Station,” in Proceedings of the Progress In Electromagnetics Research Symposium Proceedings, Guangzhou, China, 2014. View at: Google Scholar
 K. Guan, X. Lin, D. He et al., “Scenario modules and raytracing simulations of millimeter wave and terahertz channels for smart rail mobility,” in Proceedings of the 11th European Conference on Antennas and Propagation (EUCAP), pp. 113–117, Paris, France, March 2017. View at: Publisher Site  Google Scholar
 I. Val, A. Arriola, P. M. Rodriguez et al., “Wireless channel measurements and modeling for TCMS communications in metro environments,” in Proceedings of the 11th European Conference on Antennas and Propagation (EUCAP '17), pp. 108–112, Paris, France, March 2017. View at: Publisher Site  Google Scholar
 A. GonzalezPlaza, J. Moreno, I. Val et al., “5G communications in high speed and metropolitan railways,” in Proceedings of the 11th European Conference on Antennas and Propagation (EUCAP '17), pp. 658–660, Paris, France, March 2017. View at: Publisher Site  Google Scholar
 A. Chandra, P. Kukolev, A. Prokes, T. Mikulasek, and C. F. Mecklenbrauker, “UWB measurements for spatial variability and ranging: parked car in underground garage,” IEEE Antennas and Wireless Propagation Letters, vol. 16, pp. 1859–1862, 2017. View at: Publisher Site  Google Scholar
 A. Chandra, A. Prokeš, T. Mikulášek et al., “FrequencyDomain InVehicle UWB Channel Modeling,” IEEE Transactions on Vehicular Technology, vol. 65, no. 6, pp. 3929–3940, 2016. View at: Publisher Site  Google Scholar
 A. A. Saleh and R. A. Valenzuela, “A statistical model for indoor multipath propagation,” IEEE Journal on Selected Areas in Communications, vol. 5, no. 2, pp. 128–137, 1987. View at: Publisher Site  Google Scholar
 M. Corrigan, A. Walton, W. Niu, J. Li, and T. Talty, “Automatic UWB clusters identification,” in Proceedings of the IEEE Radio and Wireless Symposium (RWS '09), pp. 376–379, USA, January 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Juan Moreno GarcíaLoygorri et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.