Next Article in Journal
Integrated Smart Gas Tracking Device with Artificially Tailored Selectivity for Real-Time Monitoring Food Freshness
Next Article in Special Issue
A Compact Broadside Coupled Stripline 2-D Beamforming Network and Its Application to a 2-D Beam Scanning Array Antenna Using Panasonic Megtron 6 Substrate
Previous Article in Journal
Cybersecurity in Internet of Medical Vehicles: State-of-the-Art Analysis, Research Challenges and Future Perspectives
Previous Article in Special Issue
A Novel Monopole Ultra-Wide-Band Multiple-Input Multiple-Output Antenna with Triple-Notched Characteristics for Enhanced Wireless Communication and Portable Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

SDR-Based 28 GHz mmWave Channel Modeling of Railway Marshaling Yard

1
Postgraduate Department, China Academy of Railway Sciences, Beijing 100081, China
2
Signal and Communication Research Institute, China Academy of Railway Sciences Corporation Limited, Beijing 100081, China
3
National Research Center of Railway Intelligence Transportation System Engineering Technology, Beijing 100081, China
4
National Railway Track Test Center, Beijing 100015, China
*
Author to whom correspondence should be addressed.
Submission received: 16 August 2023 / Revised: 17 September 2023 / Accepted: 20 September 2023 / Published: 27 September 2023

Abstract

:
Compared with railway communication service requirements on the mainline, requirements in hotspots such as stations and yards are more complicated in terms of service types as well as bandwidth, of which railway-dedicated mobile communication systems such as 5G-R facilitated with dedicated frequency support cannot meet the entire communication requirements. Therefore, other radio-communication technologies need to be adopted as a supplement, among which the mmWave communication system is a promising technology, especially for large bandwidth communication between train and trackside. However, there is a lack of evaluation of the 28 GHz mmWave channel characteristics for the railway marshaling yard scenario. In this paper, the railway marshaling yard mmWave propagation scenario is deeply analyzed and classified into three typical categories, based on which, a measurement campaign is conducted using an SDR channel sounding system equipped with a 28 GHz mmWave phased-array antenna. A self-developed software under the LabVIEW platform is used to derive the channel parameters. Conclusions on the relationship between the parameters of MPC numbers, time-spread, and received power and position, as well as the impact of typical obstructions such as the Catenary, adjacent locomotives, and buildings are drawn. The statistical results and conclusions of this paper are helpful for facilitating the design and performance evaluation of future mmWave communication systems for railway marshaling yards and can also be further extended and applied to the research of mmWave utilization in 6G and other future communication technologies for more scenarios.

1. Introduction

Railway transportation, as an important pillar of national economics, is directly related to the safety of people’s lives and property and is hailed as the “artery of the country”. As a common consensus, safety is the lifeline of railway transportation. Railway radio communication plays a very important role in railway transportation with respect to transportation production as well as safety management; as a means of communication and information technology, railway radio communication reliably carries the dispatching communication service, train control information, etc., and it is the foundation of promoting production efficiency and the high-quality development of railways.
Radio frequency is the foundation and prerequisite for the innovation and development of railway radio-communication technology. In view of the high reliability requirements of railway radio-communication systems, it is widely recommended by the International Union of Railways (UIC) and other countries in the world to allocate dedicated frequency resources for railway radio-communication systems. The International Telecommunication Union (ITU), UIC, and other countries have paid close attention to the coordination of railway-dedicated frequencies and frequency protection. In 2015, the World Radio-communication Conference (WRC-15) decided to set up agenda item 1.11 to study spectrum harmonization for railway radio-communication systems between train and trackside within the existing mobile service allocations, and countries conducted in-depth research on railway radio-communication systems at the ITU level. In UIC, under the umbrella of the future railway mobile communication system (FRMCS) project, the Frequency Working Group (UGFA) was established and is now responsible for studying the related topics of frequency with respect to railway radio-communication systems.
The dedicated frequency for railway radio communication has played a significant role in railway transportation production, ensuring safety and improving efficiency, as well as achieving huge economic and social benefits. The European Union has allocated 876–880 MHz (train to ground) and 921–925 MHz (ground to train) as dedicated operating frequency bands for GSM-R in Europe, with over 100,000 km of lines currently deployed. In China, GSM-R uses 885–889 MHz (train to ground) and 930–934 MHz (train to ground), and covers over 85,000 km of railway lines.
With the trend of digitalization and intelligence, the contradiction between the increasing service demand and the shortage of frequency resources as well as the limited carrying capacity of railway radio-communication systems has become increasingly prominent. Facing the rapid evolution of radio-communication technology, the broadband of railway radio-communication systems has become a future development direction. Considering the limited frequency resources in the low-frequency band and fragmented allocation, it is an inevitable trend to migrate to the middle- and high-frequency band. In Europe, the existing 2 × 4 MHz R-GSM frequency resources of GSM-R are expanded to the 2 × 5.6 MHz (874.4–880/919.4–925 MHz) ER-GSM frequency band; in addition, 1900–1910 MHz is used as a supplementary frequency band for FRMCS(5G-R). In the meantime, China Railway has submitted an application to the Ministry of Industry and Information Technology(MIIT) for using the 2100 MHz 2 × 10 MHz frequency band as a 5G-R system field test frequency.
However, in addition to meeting the requirements of services on the main lines, such as dispatching communication, train control, monitoring, etc., service models in hotspots such as stations and yards are more complicated in terms of service types as well as bandwidth requirements. The 5G-R system facilitated with dedicated frequency support cannot meet the entire communication requirements of hotspots, so other radio-communication technologies and methods need to be adopted as a supplement to the 5G-R communication system. Among these, the mmWave communication system is a promising technology, especially for large bandwidth communication between train and trackside.
There are fundamental differences between mmWave communications and existing communication systems, in terms of high propagation loss, directivity, and sensitivity to blockage. These characteristics of mmWave communications pose several challenges to fully exploiting the potential of mmWave communications, including integrated circuits and system design, interference management, spatial reuse, anti-blockage, and dynamics control [1]. Though mmWave has the disadvantages mentioned above, rich continuous and underutilized frequency resources are available in mmWave bands. As a result, mmWave has become a key enabler for future radio-communication technologies. The 3GPP defined frequency range designation (FR2) is 24.25–52.6 GHz, and corresponding operating bands are defined in Table 1.
With the digitalization of railways, there are more and more types of sensors installed on trains, and a large number of monitoring data are collected during train operation, which need to be transmitted back to the ground center for comprehensive analysis. The traditional method is for these data to be copies by staff after the train is put into the yard or depot, which has low working efficiency and poor real-time performance. With the wide spectrum resources of mmWave and the application of Massive MIMO, beamforming, spatial multiplexing, and other technologies, the data throughput of mmWave communication systems between train and trackside can reach more than Gbps [2] during the process of trains entering the depot or passing through hotspot areas, such as train stations or marshaling yards, at very low speeds, normally 5–10 km/h, so large amounts of monitoring data such as videos collected during train operation can be transmitted during this timespan to the ground via the mmWave communication system in a few seconds for further analysis. Thus, mmWave is a promising solution for this scenario as a supplementary means of 5G-R and can significantly improve the efficiency and real-time performance of train–ground data transmission.

1.1. Related Works

Scholars have carried out a series of related work with respect to mmWave. Corresponding research will be reviewed from four categories involving communication scenarios, frequency bands, and research methodology, as well as the object of the research.
From the communication scenario point of view, railway, metro, residential house, indoor commercial, tropical area, and many other scenarios are covered by the open literature. With respect to railway scenarios, in [3], channel characteristics are studied in the 5G mmWave band centered at 25.25 GHz for typical urban HSR scenarios, including straight and curved route shapes. In [4], the 28 GHz mmWave channel characteristics of rural railway scenarios are studied via a calibrated ray-tracing simulator. Ref. [5] also focuses on the 28 GHz mmWave channel characteristics of rural railway scenarios, in which the dominant propagation mechanisms (direct, penetration, reflection, scattering, etc.) are determined and the three-dimensional environment model and the electromagnetic parameters of different objects are calibrated. The automatic coupling application scenario is studied in [6], in which a dynamic train-to-train (T2T) mmWave propagation measurement campaign was conducted, and the received signal power in the open field area was analyzed and modeled using a two-ray path loss model. The results indicated that the received power next to a platform was higher compared to the open field counterpart due to strong contributions from the signals reflected by the platform. The work in [7] studied the wireless coverage for intra-wagon scenarios at the 60 GHz band; by considering the balance of the wireless coverage and the cost of the transmitters, the deployment of two transmitters in one wagon is suggested. Typical objects in the railway propagation environment are also important research areas. In [8], a significant analysis methodology, which analyzes and distinguishes the contributions of different railway objects to the mmWave propagation channel, is proposed. In [9], the influence of typical objects to the mmWave propagation channel is analyzed for “train-to-infrastructure” and “intra-wagon” railway scenarios with various configurations. In [10], a statistical mmWave channel modeling for railway communications backhaul in 5G networks was proposed, involving 28 GHz and 60 GHz, by simulating in a novel software, NYUSIM. Viaducts and tunnels are common scenarios of railway terrain, especially in high-speed railways. The work in [11] considered measurement-based ray-tracer calibration and channel analysis for high-speed railway viaduct scenarios at 93.2 GHz, and the conclusion that the typical structure of viaducts and the application of horn antenna lead to a small value of Rician K-Factor and RMS delay spread is drawn. In [12], the authors investigate the mmWave propagation characteristics of a high-speed moving train based on field measurements in tunnel and viaduct scenarios. The measurements were carried out at 28 GHza and path loss (PL) and other channel parameters, including delay spread and Doppler shift, were investigated. In [13], smart-rail mobility was discussed, and the authors identified the main technical challenges and the corresponding chances concerning the reference scenario modules, accurate and efficient simulation platforms, beamforming strategies, and handover design. A mmWave beamforming scheme for disaster detection in high-speed railways was proposed in [14]; the key point of the study was that the antenna array generates multiple beams with different beam widths in different frequency bands simultaneously and the multiple beams are responsible for different detection areas. Another way to improve the reliability of railway communication is interlaced redundant coverage; in [15], a location-fair-based mmWave stable beamforming scheme under interlaced redundant coverage architecture is proposed to improve the stability and reliability of high-speed railway communications. Hotspots, such as depots and shunting yards, are also important railway scenarios. The work in [16] presented a comprehensive study on millimeter-wave-based mobile hotspot network (MHN) systems for high-speed train communications, including system design, field trial, channel modeling based on measurement campaign, simulation, and validation. In [17], mmWave channel measurements are reported for a railway depot environment using a wideband channel sounder operating at 60 GHz; path loss, RMS delay, as well as K-factors were extracted. The vacuum tube ultra-high-speed train (UHST) is a hot topic worldwide, and is a potential development direction of future transportation. In [18], a three-dimensional non-stationary mmWave geometry-based stochastic model (GBSM) is proposed to investigate the channel characteristics of UHST channels in vacuum tube scenarios.
Concerning the frequency bands, 26 GHz, 28 GHz, 30 GHz, 38 GHz, 40 GHz, 60 GHz, 73 GHz, 90 GHz, 93.2 GHz, and up to 300 GHz have all been studied by scholars, among which, 28 GHz and 60 GHz seem to have a higher level of attention. In [19], the authors conducted a 28 GHz-band 5G experimental trial on the actual Shinkansen that runs at a maximum speed of 283 km/h, and the experimental trial could achieve throughput exceeding 1.0 Gbps and consecutive handovers among the three BSs. In [20], field experiments of a 28 GHz-band 5G system at indoor train station platforms were carried out, and through the actual measurements, the results confirmed that the prototype 5G system could achieve mobile broadband capacity (more than 1 Gbps), even when the UE was located anywhere at the indoor train station platform. Ref [4] focuses on channel characteristics in rural railway environments at 28 GHz, and the mmWave channel characteristics of rural railway scenarios are studied via a calibrated ray-tracing simulator. In [5], influence analysis of typical objects in rural railway environments at 28 GHz is performed. In [21], millimeter-wave (mmWave) high-speed train measurements were conducted at 28 GHz with a speed up to 170 km/h in two different HST scenarios, viaduct and tunnel, in which large- and small-scale fading characteristics were extracted.
Moving on to research methodology, measure campaigns using measurement results published in the open literature, ray-tracing simulations, ultrawideband (UWB) channel sounding, deep learning, deep reinforcement learning (DRL), machine learning(ML), the space alternating generalized expectation-maximization (SAGE) algorithm, simulation platforms such as MATLAB, NYUSIM, mathematical methods such as non-convex problem with a near-optimal solution, coalition game, as well as the Bayesian approach, have all been used by scholars. Among all these methodologies, measure campaigns are basically very costly and time-intensive [22]; thus, ray-tracing simulations calibrated with measured data are commonly used, which means only a relatively small number of channel samples can be obtained, and then measurements are thus used to quantify the accuracy of a ray tracer. In [9], propagation measurements are conducted in the mmWave band for the 12 most common railway materials, the corresponding electromagnetic parameters are obtained, and a 3D ray-tracing (RT) simulator is calibrated. In [23], a series of horizontal directional scan-sounding measurements are performed inside a real high-speed train wagon at 60 GHz and 300 GHz frequency bands, and the channel characteristics are extracted, based on which, a self-developed RT simulator is validated through the reconstruction of the three-dimensional wagon model and the calibration of the electromagnetic properties of the main materials. In [5], based on the channel measurement conducted at 28 GHz with train-to-infrastructure deployment in the rural railway environment, a ray tracing (RT) simulator is calibrated in terms of environment modeling and electromagnetic calculation. In [11], ray tracer (RT) is employed to simulate the propagation in a high-speed railway viaduct scenario at a frequency of 93.2 GHz. In the article, by comparing the path loss between simulation and measurement, the permittivity and loss tangent of concrete, which is the construction material of viaducts, are calibrated in the ray tracer. In [24], based on the wideband measurements conducted in the tunnel scenario by using the “mobile hotspot network” system, 3D ray tracing (RT) is calibrated and validated to explore more channel characteristics in different HSR scenarios. In [25], since the measurement of a mmWave channel is deficient for high mobility, the proposed channel model is developed using a ray-tracing (RT) simulator that is validated with the channel measurements performed in the HST scenario at 93.2 GHz. Ref. [26] focuses on the analysis of propagation characteristics for train–ground communication systems in tunnel scenarios at both low-frequency and mmWave bands, based on ray-tracing (RT) simulation, and the material parameters in the RT simulation are calibrated by measurement data collected in realistic tunnel environments.
As for the research object, except for channel parameters such as power delay profiles (PDP), path loss, Rician K-factor, root-mean-square (RMS) delay spread, azimuth spread of arrival, and azimuth spread of departure, reconfigurable intelligent surface (RIS), beamforming, Doppler shift, energy efficient, user association, resource allocation and computation offloading, the influence of meteorological attenuation, sensing, minimal distance between base stations, as well as mmWave network architectures for railway have all been involved. In [27], the authors proposed a RIS-assisted scheduling scheme for scheduling interrupt flows and improving quality of service (QoS), and in the proposed scheme, an RIS is deployed between the BS and multiple mobile relays (MRs). In [28], multiple beams with different beamwidths are formed by the base station (BS) simultaneously to improve the system capacity, and the mobile relays (MRs) are provided with the ability to adjust the receiver (RX) beams automatically to enhance the received signal-to-noise ratio (SNR). In [29], control information and user information are transmitted through the high-frequency micro base station (BS) and the low-frequency macro BS, respectively, and an optimized beam width and power allocation scheme is proposed, which is combined with mobile relay (MR) technology to utilize the architecture of large-scale antenna beamforming. In [30], dynamic and fixed beamforming is evaluated based on the generalized model and the measured data, and the results show that the average throughput of dynamic beamforming is only 4% higher than that of fixed beamforming in the HSR tunnel, but 21% higher in the train station when severe beam misalignment is present. In [31], the authors proposed a new Doppler shift estimator for mmWave communication systems on HSRs: an equally-divided structure-based estimator (ESBE) that divides the effective orthogonal frequency-division multiplexing (OFDM) symbol into multiple equal fragments. In [32], the modeling of the Doppler effect for mmWave in HSR communications is conducted, and data-aided Doppler estimation and compensation algorithms are designed based on the new model. In [33], the authors proposed a new machine learning-based Doppler shift estimator (MLDSE), which estimates the Doppler shift by using the reference signal received power (RSRP) values measured by the mobile receiver at all times.

1.2. Contribution and Organization of the Paper

Though a significant amount of research has been carried out from the four categories reviewed above, to the authors’ best knowledge, there is still a lack of evaluation of the 28 GHz mmWave channel characteristics for the railway marshaling yard scenario. However, accurate propagation modeling and parameter analysis are crucial for achieving better performance for the radio-communication system, from the perspective of both link level and system level. To fill the aforementioned research gaps, the novelty of this work is to investigate the railway marshaling yard mmWave propagation scenarios of 28 GHz mmWave and conduct a measurement campaign with the SDR channel sounding system, using self-developed software under the LabVIEW platform to derive channel parameters, based on which, the channel characteristic and influence of typical objects are analyzed. The data collection for this paper was conducted on a sunny morning, and meteorological effects, such as the potential impact of snow or rain, was not in the scope of this study.
The rest of this paper is structured as follows. In Section 2, we will describe the railway marshaling yard propagation scenario, and three testing scenarios are designed. The software-defined radio-based channel modeling system as well as construction of the measurement environment will be presented in Section 3. Then, corresponding analysis results are summarized in Section 4. Finally, Section 5 will present the conclusions.

2. Railway Marshaling Yard Propagation Scenario

We chose the railway marshaling yard located at the National Railway Track Test Center of China Academy Railway Sciences to conduct our measurement campaign. The marshaling yard, which has 8 tracks, is rectangle-shaped with a length and width of 600 m and 40 m, respectively, as shown in Figure 1.
Three typical propagation scenarios were designed in our measurement campaign, which, to the best knowledge of the authors, cover typical scenarios of the railway marshaling yard:
The first is an open field scene with Catenary poles distributed with a distance of 31 to 47 m on our testing track, and we chose the middle track as the testing track, as shown in Figure 2.
As shown in the red circle in Figure 2, the LOS path can be blocked by the horizontal metal beam or the metal net that is used to prevent birds from building nests at some testing points on the selected track; the influence will also be analyze in the following section.
The second is a representation of a locomotive on the adjacent track. To build this environment, a 26.6-m-long locomotive was placed on the adjacent track of our testing track, as shown in Figure 3.
The third is a representation of buildings near our testing track. To collect the data of this scenario, we chose the track next to the building, as shown in Figure 4.
The parameters of the corresponding scenarios are listed in Table 2.

3. Software-Defined Radio-Based Channel Modeling

3.1. SDR mmWave Measurement System

The SDR mmWave measurement system architecture is shown in Figure 5. The receiving and transmitting ends each use a USRP-2974 as the signal receiving and transmitting main body. Both USRP-2974s are connected to the GPS antenna on their side, so the internal GPS disciplined oscillator (GPSDO) can receive the GPS signal, which is used to lock the internal clock, so that the USRP-2974 at each end can be synchronized. By using a phased-array antenna at the receiving and transmitting ends, the intermediate frequency of 2.8 GHz is moved to the 28 GHz frequency band to achieve the transmission and reception of mmWave signals.
The compositions of the transmitting and receiving ends are shown in Figure 6.
At the transmitting end, the REF reference clock locked by the GPSDO will be sent to the RF module of the USRP-2974 and output to the up-converter of the phased-array antenna to ensure system synchronization. At the same time, the GPSDO will send the pulse-per-second signal (PPS) to the FPGA of the USRP-2974 for timing triggering. A ZC sequence is generated by the USRP-2974 host and sent to the USRP-2974 FPGA. After IQ compensation and digital up-conversion, the FPGA sends the digital baseband data to the USRP-2974 RF module, where DA conversion and baseband up-conversion are completed, and the signal is up-converted to 2.8 GHz. At the phased-array antenna end, the 2.8 GHz signal is received and mixed with the local oscillator, whose frequency is 25.2 GHz, to generate the final 28 GHz mmWave RF signal, which is transmitted through the antenna array surface. In this system, the transmitting angle of the phased-array antenna is 0 degrees, which is the normal direction of the array surface.
At the receiving end, the REF reference clock locked by the GPSDO will be sent to the RF module of the USRP-2974 and output to the down-converter of the phased-array antenna to ensure system synchronization. At the same time, the GPSDO will send PPS to the FPGA of the USRP-2974 for timing triggering. The antenna array will receive a 28 GHz mmWave signal, which will be mixed with the local oscillator, whose frequency is 25.2 GHz, and down convert it into a 2.8 GHz intermediate frequency signal. The RF module of the USRP-2974 will complete baseband down-conversion and baseband sampling and transmit the collected IQ data to the FPGA of the USRP-2974. After completing IQ compensation and digital down-conversion at the FPGA end, the data will be transmitted to the Host end. At the Host end, the same ZC waveform as the transmitter side will be generated and correlated with the received IQ data, and the correlated IQ CIR will be obtained. The IQ CIR will be converted from v to dBm to obtain the power delay profile (PDP); in the meantime, the original IQ waveform and PDP data will be recorded and stored on the disk.

3.1.1. Key Hardware Components

The two most important components of the hardware of the measurement system are the receiver/transmitter USRP-2974 and the mmWave phased-array antenna mmPSA-TR64MX, as shown in Figure 7.
USRP-2974 is a high-performance software-defined radio (SDR) device produced by NI company, and is widely used for channel measurements in wireless communication systems. There are many advantages to using the USRP-2974 for channel measurements. Firstly, it provides high-performance analog-to-digital converters and processors that can achieve high-precision signal acquisition and analysis. Secondly, it supports a wide range of programming tool and development libraries, including GNU Radio, LabVIEW, and MATLAB, so that users can quickly customize their own channel measurement applications. The main parameters of the USRP-2974 are listed in Table 3.
The mmWave phased-array antenna mmPSA-TR64MX is a 64-unit two-dimensional mmWave phased-array antenna, the duplex mode of which is TDD, and the operating frequency band is 27–29 GHz. It supports the 3GPP 5G FR2 n257 frequency band and can be used for 5G mmWave wireless communication. The mmPSA-TR64MX is equipped with 64 phase-shifting transceiver channels, which can achieve high-performance dynamic beamforming of 4096 states in a two-dimensional ± 60 degrees range, and the EIRP value is 51 dBm. The array can be individually turned on or off per 2 × 2 antenna unit to achieve flexible array reconstruction. The phased array is internally integrated with a bidirectional frequency conversion circuit and a high-performance mmWave local oscillator, realizing a 2–4 GHz S band IF transceiver interface, which is convenient for connection with the low-frequency band digital sampling system. The beam characteristics of the mmWave phased-array antenna are shown in Figure 8:
The main parameters of the mmWave phased-array antenna mmPSA-TR64MX are listed in Table 4.

3.1.2. Self-Developed Sounding Data Process Software

The Laboratory Virtual Instrument Engineering Workbench (LabVIEW) platform was employed to develop our self-developed software in order to analyze the collected data. LabVIEW is a programming environment developed by NI company that uses graphical editing G language to write programs and is perfectly compatible with the USRP-2974 used in this system. The program developed by the authors for analyzing PDP and corresponding channel parameters based on collected data is shown in Figure 9.

3.2. Construction of Measurement Environment

The transmitting end of the measurement system was placed on the bridge over the railway tracks; the height of the bridge from the rail surface is 14 m and the transmitting antenna fixed on a tripod is 1.75 m high from the bridge, as shown in Figure 10.
The antenna of the receiving end is fixed on a tripod, which is placed at the middle of the selected track, and the receiving antenna is 1.6 m high from the rail surface, as shown in Figure 11.
In order to meet the data collection requirements described in Section 2, two tracks are selected for testing, as shown in Figure 12; the track marked with the red circle is selected for the first and second scenarios, and is denoted as T 1 & 2 . The track marked with the orange rectangle is selected for the third scenario, and is denoted as T 3 . The data collection campaign was conduct on the corresponding track, with an interval of 5 m, for further analysis. On T 1 & 2 , data collection points with or without Catenary between the TX/RXs were separated into two categories for performance comparison.

3.3. Corresponding Algorithms

The IQ sampling rate used in the measurement system is 200 MS/s, and the corresponding PDP resolution is 1/200 M = 5 ns. The period of the ZC sequence is 10.24 µs.
In Figure 13, the blue blocks represent the sent ZC sequence and green blocks represent the received ZC sequence.
The process flow is shown in Figure 14.

3.3.1. ZC Sequence Selection

In our measurement system, the sample rate of the ZC waveform is 200 MHz, the bandwidth is 160 MHz, and the length is 2048 points. The expression of the selected ZC sequence is:
x [ n ] = e j π u n 2 N , 0 n N 1 , N = 2048
where u is the root index of the ZC sequence.

3.3.2. Discrete Algorithm of Cross-Correlation

Define h as a discrete sequence, whose index can be negative. N is the number of elements in input sequence X, M is the number of elements in sequence Y, and assuming
x j = 0 , j < 0 o r j N
y j = 0 , j < 0 o r j M
h is calculated as:
h j = k = 0 N 1 x k × y j + k
where j [ ( N 1 ) , ( M 1 ) ] .
The cross-correlation is calculated as:
R x y i = h i ( N 1 )
where i [ 0 , N + M 2 ] .

4. Results

As mentioned in Section 3.2, data samples are collected on the selected Tracks T 1 & 2 and T 3 using the self-developed LabVIEW data processing software mentioned in Section 3.1.2. The channel sounding data are processed and the corresponding channel parameters, including received power and time delay, as well as the phase value of each multipath component, are extracted; at every position, 30 samples of each parameter are obtained, and the 95% statistical value are calculated for each parameter. The 95% statistical value of channel parameters at the position of 12.4 m (horizontal distance from the bridge) on T 1 & 2 are shown in Table 5 and Figure 15; here, we show that the index of the MPC increases as the power decreases and the delay increases
Based on all the channel parameters derived from the collected sounding data, further analyses are made from the following two aspects.

4.1. Relationship between Parameters and Position

4.1.1. Relationship between Number of MPCs and Position

As the index of the path increases, the quality of the sound data received from the path deteriorates, and the probability of extracting channel parameters from it decreases accordingly. If the number of decoded channel parameters exceeds 85% of the total decoding times for a specific path, then the path is determined as an effective path. The number of effective paths at each position is determined based on this principle, and the relationship between number of MPCs and position is shows in Figure 16.
As the distance between Tx and Rx increases, there is no significant change in the number of MPCs, which fluctuates randomly between 7 and 11; the mean value μ = 9.55 and the variance δ = 0.868 .

4.1.2. Relationship between Time-Spread and Position

For a specific sampling location, we denote the 95% statistical value of delay as D f i r s t , 95 % for the first effective MPC and D l a s t , 95 % for the last effective MPC. The time-spread, which is denoted as T S for this location, is defined as:
T S = D l a s t , 95 % D f i r s t , 95 %
Using (6), the time-spread of every sampling position on T 1 is calculated, and the relationship between time-domain parameters and position is shown in Figure 17 and Figure 18.
Figure 17 shows the relationship between the delay of the first effective path and position. The nonlinear curve fitting method is used to fit the curve in Figure 17; the fitting curve is the red line in the figure, and the fitting equation is:
y = a × x b
where a = 18.04322 ± 3.86076 and b = 0.74836 ± 0.0436 .
Figure 18 shows the relationship between time-spread and position; at the positions of 152.4 m and 202.4 m, there is a large time-spread sample in each of these two positions, which will be further analyzed in the section on the impact of typical objects. If we exclude these two samples for now, we get Figure 19.
The nonlinear curve fitting method is used to fit the curve in Figure 19; the fitting curve is the red line in the figure, and the fitting equation is:
y = a 1 + e ( k × ( x x c ) )
where a = 2212.6593 ± 778.51346 , k = 0.02413 ± 0.0063 and x c = 143.79056 ± 36.49551 .

4.1.3. Relationship between Received Power and Position

The 95% statistical received power value of all effective paths in each position are calculated and are shown in Figure 20.
The 95% statistical power values are summed by position, and we get Figure 21.
The nonlinear curve fitting method is used to fit the curve in Figure 21; the fitting curve is the red line in the figure, and the fitting equation is:
y = b × ln ( x a )
where a = 2.45198 ± 3.83339 and b = 7.2857 ± 0.15875 . There is deep fading at the position of 152.2 m as well as 202.4 m, which will be further analyzed in the section on the impact of typical objects.

4.2. Impact of Typical Obstructions

4.2.1. Impact of the Catenary

We use the parameters derived from the data collected on T 1 & 2 to analyze the impact of the Catenary. As we mentioned in Section 2, the LOS path can be blocked by the horizontal metal beam or the metal net at some testing points on the selected track, as shown in Figure 2, and these testing points were recorded in detail during the test and are listed in Table 6.
From Figure 16 and Figure 17, we can tell that there are no obvious changes with respect to the number of MPCs and the delay of the first effective path at the blocked positions, and the LOS path, which is the first effective path, can penetrate the metal net and reach the receiving end, even though penetration loss is introduced, which we will further analyze. However, from Figure 18 and Figure 21, we can see that, in terms of the time-spread and received power, significant changes have taken place in these blocked positions.
In Figure 18, at the position of 37.4 m, we can see the time-spread does not change much. It is only increased by 118 ns compared with the expected value; however, at the positions of 152.4 m and 202.4 m, the time-spread compared with the expected value is increased by 2970 ns and 2853 ns, respectively. This is because, as the distance between Tx and Rx increases, the number of Catenary between the Tx and Rx increases. The delay of the first path, which is the LOS path, is basically stable, but the reflection and scattering conditions of other paths reaching the blocked position are more complex, so the delay spread increases obviously.
In Figure 21, at all the blocked positions, deep fading occurs; we can see at the positions of 37.4 m, 152.4 m, and 202.4 m that the penetration loss is 4.63 dB, 5.31 dB, and 10.33 dB, respectively. As mentioned above, as the distance between Tx and Rx increases, the number of Catenary between Tx and Rx increases, and the density of Catenary increases accordingly; therefore, the probability of being blocked by Catenary increases. We noticed that the penetration loss at the position of 202.4 m doubled compared with the former two positions; this is because, at this position, the path was blocked by Catenary twice, as we observed during the test. We can see that the penetration loss of the Catenary is approximately 5.07 dB.

4.2.2. Impact of Adjacent Locomotive and Building

As we mentioned in Section 2, the locomotive, whose length is 26.6 m, was parked on the adjacent track of T 1 and covered the position from 52.4 m to 79 m.
From Figure 16, we can see the number of MPCs of positions near this area slightly increase by one to two paths. In Figure 18, we can see the time-spread of positions near this area slightly increases; this is because, with the introduction of the locomotive, the reflection and scattering conditions of paths become more complex. We can quantify the impact on the time-spread with M A P E :
M A P E = 1 n i = 1 n Y i Y ^ i Y i × 100 %
We select data from 42.4 m to 102.4 m, which are near the locomotive parking area, for fitting in order to obtain a higher fitting accuracy, and the actual time-spread data in the green circle in Figure 22 are used as well as the fitting data to calculate M A P E , as show in Figure 22.
The linear fitting equation is:
y = a + b × x
where a = 84 and b = 2.5 . The actual data marked in the green circle are used to calculate M A P E , and we obtain M A P E = 38.63 % .
Now we analyze the impact of the building, which much larger than the locomotive, with a 7.5-m-high canopy and a 10-m-high main building, as shown in Figure 4. Since the locomotive and the building are basically set in parallel, we conduct the analysis through comparison of these two different objects, as shown in Figure 23.
In Figure 23a, the difference of MPC numbers derived from samples collected on T 1 & 2 and T 3 is not obvious, with a deviation of only 1–2 paths randomly, and a maximum number of 12 MPCs on T 3 .
In Figure 23b, the lines represent the 95% statistical delay value of the first effective path and indicate that the delay parameters derived from samples collected on T 1 & 2 gradually increase as the distance increases, while the delay parameters derived from samples collected on T 3 exhibit an irregular fluctuation state as the distances increases.
In Figure 23c, the lines show the 95% statistical power values summed by position, and we see that the received power of T 1 & 2 gradually decrease as the distance increases, while the delay of T 3 shows an irregular fluctuation state as the distances increases.
In Figure 23d, the lines show the time-spread of the first and the last effective path, and we see that the time-spread of T 1 & 2 remains stable within a small range of 30 m, while the time-spread of T 3 shows an irregular fluctuation state as the distances increases.
Based on the above analysis, we can draw the conclusion that, due to the larger size and more irregular shape of the building, the received power, delay, and time-spread channel parameters within a small position range show dramatic random variations. Since normally there are several tracks in the railway marshaling yard, we recommend that locomotives with unloading data tasks travel on the tracks away from buildings, in order to obtain a more stable wireless transmission channel. In addition, considering the penetration loss caused by the Catenary, an additional signal-level protection margin needs to be added to achieve redundant protection and improve the reliability of the data transmission.

5. Conclusions

In this work, the railway marshaling yard mmWave propagation scenario is deeply analyzed and classified into three typical categories, involving open scene with Catenary poles, a locomotive parked on the adjacent track, and a building near the testing track.
A measurement campaign is conducted at the National Railway Track Test Center of CARS according to the three propagation scenarios, using an SDR channel sounding system, the hardware of which mainly consists of a USRP-2974 transceiver and a mmPSA-TR64MX phased-array antenna. A self-developed software under the LabVIEW platform is used to derive the channel parameters; to achieve this, a ZC sequence of 200 MHz sample rate and 2048 points is selected, and a discrete algorithm of cross-correlation is designed. Parameters involving received power, delay, and time-spread are successfully derived for every sampling position on our selected tracks of T 1 & 2 and T 3 , and the following four conclusions can be drawn based on these parameters:
(1) As the distance between Tx and Rx increase, there is no significant change in the number of MPCs, which fluctuates randomly between 7 and 11, with a mean value of μ = 9.55 and variance of δ = 0.868 .
(2) The relationship between parameter and position is deeply analyzed and fitted with proper equations. The relationship between the delay of the first effective path and position conforms to the equation y = a × x b , the relationship between the time-spread and the position conforms to the equation y = a 1 + e ( k × ( x x c ) ) , and the relation between received power and position conforms to the equation y = b × ln ( x a ) .
(3) The impact of the Catenary is an important factor. As the distance between Tx and Rx increases, the number of Catenary between the Tx and Rx increases, and the density of Catenary increases accordingly; therefore, the probability of being blocked by Catenary increases. This phenomenon affects mmWave propagation in two ways: With respect to received power, at all the blocked positions, deep fading occurs, with a penetration loss of 5.07 dB. As for time-spread, at the blocked position of 37.4 m, the time-spread does not change much, only increasing by 118 ns compared with the expected value; however, at the blocked positions of 152.4 m and 202.4 m, the time-spread compared with the expected value is increased by 2970 ns and 2853 ns, respectively, which means that as the distance between Tx and Rx increases, the time-spread increases accordingly at blocked positions.
(4) The impact of the locomotive parked on the adjacent track is not obvious; the time-spread slightly increases, with a M A P E of 38.63%. However, due to the larger size and more irregular shape of the building, the received power, delay, and time-spread channel parameters within a small position range show dramatic random variations, so the recommendation that locomotives with unloading data tasks travel on the tracks away from buildings is proposed.
The conclusions of this paper are helpful for facilitating the design of mmWave communication systems for railway marshaling yards, to fully exploit the potential of mmWave communications. In the meantime, 6G networks are expected to utilize ultra-high-frequencies such as mmWave in typical scenarios involving smart factory and cell-free communication, as well as non-terrestrial communication [34]. The relevant analysis methods employed in this paper can also be further extended and applied to research of mmWave utilization in 6G and future communication technologies for more scenarios. Moreover, meteorological effects are of great importance in the research of mmWave. In the future, we will conduct further data collection and analysis based on different weather conditions.

Author Contributions

Conceptualization, Y.L. (Yiqun Liang) and H.L.; methodology, Y.L. (Yiqun Liang), Y.T. and Y.L. (Yi Li); software, Y.L. (Yiqun Liang); validation, Y.L. (Yiqun Liang) and W.W.; formal analysis, Y.L. (Yiqun Liang) and H.L.; writing—original draft preparation, Y.L. (Yiqun Liang); writing—review and editing, Y.L. (Yiqun Liang); visualization, Y.L. (Yiqun Liang); supervision, Y.L. (Yiqun Liang) and H.L.; project administration, Y.L. (Yiqun Liang) and H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Technology R&D Program of China State Railway Group Co., Ltd., grant number SY2021G001, and the Science and Technology R&D Program of China Academy of Railway Sciences Co., Ltd., grant number of 2022YJ033.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

This research was supported by Liu Liang from Beijing Zengyi Technology Company, in terms of validation of the SDR measurement system.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
3GPPThe 3rd Generation Partnership Project
5G-R5G for Railway
ADCAnalog-to-Digital Converter
BSBase Station
CARSChina Academy of Railway Sciences
CIRChannel Impulse Response
DRLDeep Reinforcement Learning
EIRPEffective Isotropic Radiated Power
EMElectromagnetic
FRMCSFuture Railway Mobile Communication System
G/TGain-to-Noise Temperature Ratio
GBSMGeometry-Based Stochastic Model
GPSDOGPS Disciplined Oscillator
GSM-RGSM for Railway
HSRHigh-speed Railway
IQIn-phase & Quadrature
ITUInternational Telecommunication Union
LabVIEWLaboratory Virtual Instrument Engineering Workbench
LOSLine of Sight
MHNMobile Hotspot Network
MRMobile Relay
PDPPower Delay Profile
PLPath Loss
QoSQuality of Service
RISReconfigurable Intelligent Surface
RMSRoot Mean Square
RSRPReference Signal Received Power
RTRay Tracing
SAGESpace Alternating Generalized Expectation-Maximization
SDRSoftware Defined Radio
sFDRspurious-Free Dynamic Range
SNRSignal-to-Noise Ratio
T2TTrain-to-Train
UEUser Equipment
UGFAUIC Group for Frequency Aspects
UHSTUltra-High-Speed Train
UICInternational Union of Railway
UWBUltrawideband
WRCWorld Radiocommunication Conference
ZC SequenceZadoff-Chu Sequence

References

  1. Niu, Y.; Li, Y.; Jin, D.; Su, L.; Vasilakos, A.V. A survey of millimeter wave communications (mmWave) for 5G: Opportunities and challenges. Wirel. Netw. 2015, 21, 2657–2676. [Google Scholar] [CrossRef]
  2. Zhang, G.; Wang, L.; Gao, H.; Lin, H.; Wang, X.; Cao, H.; Liang, Y. Test of key intelligent technology of Haolebaoji-Ji’an Railway. China Railw. 2022, 11, 23–32. [Google Scholar]
  3. He, D.; Ai, B.; Guan, K.; Zhong, Z.; Tian, L.; Dou, J. Channel analysis for millimeter-wave railway communications in urban environment. In Proceedings of the 2017 XXXIInd General Assembly and Scientific Symposium of the International Union of Radio Science (URSI GASS), Montreal, QC, Canada, 19–26 August 2017; pp. 1–4. [Google Scholar]
  4. Wang, L.; Ai, B.; Guan, K.; He, D.; Zhong, Z.; Kim, J. Channel Characteristics in Rural Railway Environment at 28 GHz. In Proceedings of the 16th International Conference on Intelligent Transportation Systems Telecommunications (ITST), Lisbon, Portugal, 15–17 October 2018; pp. 1–5. [Google Scholar]
  5. He, D.; Guan, K.; Ai, B.; Zhong, Z.; Wang, L.; Kim, J. Influence Analysis of Typical Objects in Rural Railway Environments at 28 GHz. IEEE Trans. Veh. Technol. 2019, 68, 2066–2076. [Google Scholar] [CrossRef]
  6. Soliman, M.; Schneckenburger, N.; Michailow, N.; Deneire, L.; Liénard, M.; De Doncker, P.; Moeneclaey, M. Dynamic Train-to-Train Propagation Measurements in the Millimeter Wave Band—Campaign and First Results. In Proceedings of the 2019 13th European Conference on Antennas and Propagation (EuCAP), Krakow, Poland, 31 March–5 April 2019; pp. 1–5. [Google Scholar]
  7. Yi, H.; Guan, K.; He, D.; Ai, B.; Zhong, Z.; Kim, J. Wireless Coverage Analysis for Intra-Wagon Scenario at 60 GHz Band. In Proceedings of the 2018 12th International Symposium on Antennas, Propagation and EM Theory (ISAPE), Hangzhou, China, 3–6 December 2018; pp. 1–4. [Google Scholar]
  8. He, D.; Ai, B.; Guan, K.; Zhong, Z.; Tian, L.; Dou, J. Significance Analysis for Typical Objects in mmWave Urban Railway Propagation Environment. In Proceedings of the 2017 IEEE Globecom Workshops (GC Wkshps), Singapore, 4–8 December 2017; pp. 1–5. [Google Scholar]
  9. He, D.; Ai, B.; Guan, K.; Zhong, Z.; Wang, L.; Kim, J. Influence of Typical Railway Objects in a mmWave Propagation Channel. IEEE Trans. Veh. Technol. 2018, 67, 2880–2892. [Google Scholar] [CrossRef]
  10. Dutty, H.B.H.; Mowla, M.M.; Mou, M.A. A Statistical mmWave Channel Modeling for Railway Communications Backhaul in 5G Networks. In Proceedings of the 2019 3rd International Conference on Electrical, Computer and Telecommunication Engineering (ICECTE), Rajshahi, Bangladesh, 26–28 December 2019; pp. 121–124. [Google Scholar]
  11. Lin, X.; Xu, J.; Yi, H.; Ai, B.; Zhong, Z.; Kim, J. Measurement based ray tracer calibration and channel analysis for high-speed railway viaduct scenario at 93.2 GHz. In Proceedings of the 2017 IEEE International Symposium on Antennas and Propagation and USNC/URSI National Radio Science Meeting, San Diego, CA, USA, 9–14 July 2017; pp. 617–618. [Google Scholar]
  12. Park, J.-J.; Lee, J.; Kim, K.-W.; Kim, M.-D. 28-GHz High-Speed Train Measurements and Propagation Characteristics Analysis. In Proceedings of the 2020 14th European Conference on Antennas and Propagation (EuCAP), Copenhagen, Denmark, 15–20 March 2020; pp. 1–5. [Google Scholar]
  13. Guan, K.; He, D.; Hrovat, A.; Ai, B.; Zhong, Z.; Kürner, T. Challenges and chances for smart rail mobility at mmWave and THz bands from the channels viewpoint. In Proceedings of the 2017 15th International Conference on ITS Telecommunications (ITST), Warsaw, Poland, 29–31 May 2017; pp. 1–5. [Google Scholar]
  14. He, L.; Fang, X.; Li, H.; Li, C.; Liu, Y. An mmWave beamforming scheme for disaster detection in high speed railway. In Proceedings of the 2016 IEEE/CIC International Conference on Communications in China (ICCC), Chengdu, China, 27–29 July 2016; pp. 1–6. [Google Scholar]
  15. Yin, H.; Jiang, R.; Xu, Y.; Li, D.; Tian, F. Location-Fair Based mmWave Stable Beamforming Scheme for High Speed Railway. In Proceedings of the 2019 IEEE International Conference on Consumer Electronics - Taiwan (ICCE-TW), Yilan, Taiwan, 20–22 May 2019; pp. 1–2. [Google Scholar]
  16. Kim, J.; Yi, H.; Park, J.-J.; Lee, J.; Zhong, Z.; Ai, B. A Comprehensive Study on mmWave-Based Mobile Hotspot Network System for High-Speed Train Communications. IEEE Trans. Veh. Technol. 2019, 68, 2087–2101. [Google Scholar] [CrossRef]
  17. Bulut, B.; Güneş, F.; Akkaya, G.; Karaçalı, A.; Kurnaz, F.; Çimen, S.B.; Güvendi, I.; Çayören, M. Millimeter wave channel measurements in a railway depot. In Proceedings of the 2017 IEEE 28th Annual International Symposium on Personal, Indoor, and Mobile Radio Communications (PIMRC), Montreal, QC, Canada, 8–13 October 2017; pp. 1–5. [Google Scholar]
  18. Xu, Y.; Yu, K.; Li, L.; Lei, X.; Hao, L.; Wang, C.-X. A 3D Non-stationary MmWave Channel Model for Vacuum Tube Ultra-High-Speed Train Channels. In Proceedings of the 2021 IEEE Wireless Communications and Networking Conference (WCNC), Nanjing, China, 29 March–1 April 2021; pp. 1–6. [Google Scholar]
  19. Nonaka, N.; Nakamura, T.; Ito, K.; Kumagai, T.; Kishiyama, Y.; Benjebbour, A.; Li, J.; Takeda, K. 28 GHz-Band Experimental Trial at 283 km/h Using the Shinkansen for 5G Evolution. In Proceedings of the 2020 IEEE 91st Vehicular Technology Conference (VTC2020-Spring), Antwerp, Belgium, 25–28 May 2020; pp. 1–5. [Google Scholar]
  20. Okano, M.; Hasegawa, Y.; Kanai, K.; Wei, B.; Katton, J. Field Experiments of 28 GHz Band 5G System at Indoor Train Station Platform. In Proceedings of the 2020 IEEE 17th Annual Consumer Communications and Networking Conference (CCNC), Las Vegas, NV, USA, 10–13 January 2020; pp. 1–6. [Google Scholar]
  21. Park, J.-J.; Lee, J.; Kim, K.-W.; Kim, M.-D. Large- and Small-Scale Fading Characteristics of mmWave HST Propagation Channel Based on 28-GHz Measurements. In Proceedings of the 2021 15th European Conference on Antennas and Propagation (EuCAP), Dusseldorf, Germany, 22–26 March 2021; pp. 1–5. [Google Scholar]
  22. Hur, S.; Kim, T.; Lee, J.; Kang, J.; Kim, K. Proposal on Millimeter-Wave Channel Modeling for 5G Cellular System. IEEE J. Sel. Top. Signal Process. 2016, 10, 454–469. [Google Scholar] [CrossRef]
  23. Guan, K.; Ai, B.; He, D.; Wang, L.; Zhong, Z.; Kim, J. Channel Sounding and Ray Tracing for Intrawagon Scenario at mmWave and Sub-mmWave Bands. IEEE Trans. Antennas Propag. 2021, 69, 1007–1019. [Google Scholar] [CrossRef]
  24. He, D.; Ai, B.; Guan, K.; Wang, L.; Zhong, Z.; Kim, J. Channel Measurement, Simulation, and Analysis for High-Speed Railway Communications in 5G Millimeter-Wave Band. IEEE Trans. Intell. Transp. Syst. 2018, 19, 3144–3158. [Google Scholar] [CrossRef]
  25. Yang, J.; Wang, Y.; Zhang, Y.; Li, J.; Liu, Y.; Zhang, W. A Geometry-Based Stochastic Channel Model for the Millimeter-Wave Band in a 3GPP High-Speed Train Scenario. IEEE Trans. Veh. Technol. 2018, 67, 3853–3865. [Google Scholar] [CrossRef]
  26. Zhou, T.; Li, H.; Sun, R.; Wang, Y.; Liu, L.; Tao, C. Simulation and Analysis of Propagation Characteristics for Tunnel Train-Ground Communications at 1.4 and 40 GHz. IEEE Access 2019, 7, 105123–105131. [Google Scholar] [CrossRef]
  27. Li, P.; Wang, L.; Ai, B.; Guan, K.; Zhong, Z.; Kim, J. RIS-Assisted Scheduling for High-Speed Railway Secure Communications. IEEE Trans. Veh. Technol. 2023, 72, 3488–3501. [Google Scholar] [CrossRef]
  28. Yin, H.; Jiang, R.; Xu, Y. An mmWave-Based Adaptive Multi-Beamforming Scheme for High Speed Railway Communications. In Proceedings of the 2018 10th International Conference on Wireless Communications and Signal Processing (WCSP), Hangzhou, China, 18–20 October 2018; pp. 1–6. [Google Scholar]
  29. Zhang, Q.; Zhao, J.; Chen, A.; Wang, C. Adaptive Beamforming and Power Allocation for mmWave Communication in High-Speed Railway. Radio Sci. 2021, 56, 1–11. [Google Scholar] [CrossRef]
  30. Yu, D.; Ai, B.; Guan, K.; Wang, L.; Zhong, Z.; Kim, J. Empirical Study on Directional Millimeter-Wave Propagation in Railway Communications Between Train and Trackside. IEEE J. Sel. Areas Commun. 2020, 38, 2931–2945. [Google Scholar] [CrossRef]
  31. Zhang, C.; Wang, G.; He, R.; Zou, Y. A Doppler Shift Estimator for Millimeter-Wave Communication Systems on High-Speed Railways. In Proceedings of the 2018 IEEE/CIC International Conference on Communications in China (ICCC), Beijing, China, 16–18 August 2018; pp. 227–231. [Google Scholar]
  32. Gong, Z.; Li, C.; Jiang, F.; Win, M.Z. Data-Aided Doppler Compensation for High-Speed Railway Communications Over mmWave Bands. IEEE Trans. Wireless Commun. 2021, 20, 520–534. [Google Scholar] [CrossRef]
  33. Kim, T.; Ko, K.; Hwang, I.; Hong, D.; Choi, S.; Wang, H. RSRP-Based Doppler Shift Estimator Using Machine Learning in High-Speed Train Systems. IEEE Trans. Veh. Technol. 2021, 70, 371–380. [Google Scholar] [CrossRef]
  34. Alsabah, M.; Naser, M.A.; Mahmmod, B.M.; Abdulhussain, S.H.; Eissa, M.R.; Al-Baidhani, A.; Noordin, N.K.; Sait, S.M.; Al-Utaibi, K.A.; Hashim, F. 6G wireless communications networks: A comprehensive survey. IEEE Access 2021, 9, 148191–148243. [Google Scholar] [CrossRef]
Figure 1. Railway marshaling yard. (a) Sketch map of the railway marshaling yard. (b) Picture of the railway marshaling yard.
Figure 1. Railway marshaling yard. (a) Sketch map of the railway marshaling yard. (b) Picture of the railway marshaling yard.
Sensors 23 08108 g001
Figure 2. Open field scene with Catenary poles.
Figure 2. Open field scene with Catenary poles.
Sensors 23 08108 g002
Figure 3. Locomotive on adjacent Track.
Figure 3. Locomotive on adjacent Track.
Sensors 23 08108 g003
Figure 4. Building adjacent to the testing track.
Figure 4. Building adjacent to the testing track.
Sensors 23 08108 g004
Figure 5. SDR mmWave measurement system architecture.
Figure 5. SDR mmWave measurement system architecture.
Sensors 23 08108 g005
Figure 6. Composition diagram of transmitting and receiving ends. (a) Transmitting end. (b) Receiving end.
Figure 6. Composition diagram of transmitting and receiving ends. (a) Transmitting end. (b) Receiving end.
Sensors 23 08108 g006
Figure 7. Receiver/transmitter and mmWave phased-array antenna. (a) USRP 2974. (b) mmWave phased-array antenna mmPSA-TR64MX.
Figure 7. Receiver/transmitter and mmWave phased-array antenna. (a) USRP 2974. (b) mmWave phased-array antenna mmPSA-TR64MX.
Sensors 23 08108 g007
Figure 8. Beam Characteristics of the mmWave phased-array antenna. (a) Horizontal beam. (b) Vertical beam. (c) Beam Width.
Figure 8. Beam Characteristics of the mmWave phased-array antenna. (a) Horizontal beam. (b) Vertical beam. (c) Beam Width.
Sensors 23 08108 g008
Figure 9. Self-developed software based on LabVIEW.
Figure 9. Self-developed software based on LabVIEW.
Sensors 23 08108 g009
Figure 10. Construction of the transmitting end of the measurement system. (a) Bridge over the tracks. (b) Transmitting end on the bridge.
Figure 10. Construction of the transmitting end of the measurement system. (a) Bridge over the tracks. (b) Transmitting end on the bridge.
Sensors 23 08108 g010
Figure 11. Construction of the receiving end of the measurement system. (a) Receiving antenna on the selected track. (b) Equipment of the receiving end.
Figure 11. Construction of the receiving end of the measurement system. (a) Receiving antenna on the selected track. (b) Equipment of the receiving end.
Sensors 23 08108 g011
Figure 12. Selection of tracks.
Figure 12. Selection of tracks.
Sensors 23 08108 g012
Figure 13. ZC sequence synchronization of the transmitting and receiving ends.
Figure 13. ZC sequence synchronization of the transmitting and receiving ends.
Sensors 23 08108 g013
Figure 14. Process flow of the measurement system.
Figure 14. Process flow of the measurement system.
Sensors 23 08108 g014
Figure 15. The channel parameters at the position of 12.4m. The horizontal axis coordinates 1# to 10# represent the index of MPC, which increases as the power decreases and the delay increases. (a) Delay. (b) Received Power.
Figure 15. The channel parameters at the position of 12.4m. The horizontal axis coordinates 1# to 10# represent the index of MPC, which increases as the power decreases and the delay increases. (a) Delay. (b) Received Power.
Sensors 23 08108 g015
Figure 16. Relationship between number of MPCs and position.
Figure 16. Relationship between number of MPCs and position.
Sensors 23 08108 g016
Figure 17. The relationship between delay of the first effective path and position.
Figure 17. The relationship between delay of the first effective path and position.
Sensors 23 08108 g017
Figure 18. The relationship between time-spread and position.
Figure 18. The relationship between time-spread and position.
Sensors 23 08108 g018
Figure 19. Fitting of time-spread.
Figure 19. Fitting of time-spread.
Sensors 23 08108 g019
Figure 20. The relationship between received power and position.
Figure 20. The relationship between received power and position.
Sensors 23 08108 g020
Figure 21. The relationship between summed received power and position.
Figure 21. The relationship between summed received power and position.
Sensors 23 08108 g021
Figure 22. Fitting of time-spread from 42.4 m to 102.4 m.
Figure 22. Fitting of time-spread from 42.4 m to 102.4 m.
Sensors 23 08108 g022
Figure 23. Comparison of the locomotive and the building. (a) Number of MPCs. (b) Received power. (c) Delay of the first effect path. (d) Time-spread.
Figure 23. Comparison of the locomotive and the building. (a) Number of MPCs. (b) Received power. (c) Delay of the first effect path. (d) Time-spread.
Sensors 23 08108 g023
Table 1. NR operating bands in FR2.
Table 1. NR operating bands in FR2.
Operating BandUplink (GHz)Downlink (GHz)Duplex Mode
n25726.5–29.526.5–29.5TDD
n25824.25–27.524.25–27.5TDD
n26037–4037–40TDD
n26127.5–28.3527.5–28.35TDD
Table 2. Parameters of corresponding scenarios.
Table 2. Parameters of corresponding scenarios.
ParameterDetail Information
Length of the yard600 m
Width of the yard40 m
Number of tracks8
Track gauge1.435 m
Distance between tracks5.05 m
Distance between Catenary poles31–47 m
Hight of Catenary poles10 m
Hight of the buildingcanopy 7.5 m; building 10 m
Distance between the track and the building6 m
Size of the locomotiveLength × Width × Height (m) = 26.6 × 3.105 × 4.433
Table 3. Parameters of USRP-2974.
Table 3. Parameters of USRP-2974.
ParametersDetail Information
ADC resolution14 bit
DAC resolution16 bit
BasebandMaximum I/Q sample rate200 MS/s
ADC sFDR88 dB
DAC sFDR80 dB
Number of channels2
Frequency range10 MHz to 6 GHz
Frequency step<1 kHz
TransmitterMaximum output power5 mW to 100 mW
Gain range0 dB to 31.5 dB
Gain step0.5 dB
Maximum instantaneous real-time bandwidth160 MHz
Number of channels2
Frequency range10 MHz to 6 GHz
Frequency step<1 kHz
ReceiverGain range0 dB to 37.5 dB
Gain step0.5 dB
Maximum input power10 dBm
Noise figure5 dB to 7 dB
Maximum instantaneous real-time bandwidth160 MHz
Table 4. Parameters of the mmWave phased-array antenna mmPSA-TR64MX.
Table 4. Parameters of the mmWave phased-array antenna mmPSA-TR64MX.
ParametersDetail Information
Frequency Band27–29 GHz
Duplex modeTDD
Array size8 × 8
Unit spacing0.5 × λ
G/T≥−7 dB/K
Polarization modeVertical polarization
Sidelobe suppression≥15 dB @ normal-direction
Beam scanning range±60° @ Horizontal&Vertical direction
Beam switching speed<3 μ s
Transceiver switching speed<3 μ s
EIRP51 dBm
Gain dynamic range≥30 dB
Gain adjustment accuracy1 dB
Table 5. The channel parameters at the position of 12.4 m on T 1 & 2 .
Table 5. The channel parameters at the position of 12.4 m on T 1 & 2 .
Multipath ComponentDelay (ns) (95%)Received Power (dBm) (95%)
1165−21.15
2185−36.078
3210−43.687
4235−48.653
5260−52.216
6285−53.432
7365−56.995
8385−58.209
9405−58.767
10410−59.293
Table 6. Blocking positions on T 1 .
Table 6. Blocking positions on T 1 .
Position (m)Block Object
37.4Middle of the medal net
42.4Edge of the medal net
132.4Edge of the medal net
152.4Middle of the medal net
162.4Edge of the medal net
202.4Edge of the medal net
212.4Middle of the medal net (twice)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liang, Y.; Li, H.; Tian, Y.; Li, Y.; Wang, W. SDR-Based 28 GHz mmWave Channel Modeling of Railway Marshaling Yard. Sensors 2023, 23, 8108. https://doi.org/10.3390/s23198108

AMA Style

Liang Y, Li H, Tian Y, Li Y, Wang W. SDR-Based 28 GHz mmWave Channel Modeling of Railway Marshaling Yard. Sensors. 2023; 23(19):8108. https://doi.org/10.3390/s23198108

Chicago/Turabian Style

Liang, Yiqun, Hui Li, Yuan Tian, Yi Li, and Wenhua Wang. 2023. "SDR-Based 28 GHz mmWave Channel Modeling of Railway Marshaling Yard" Sensors 23, no. 19: 8108. https://doi.org/10.3390/s23198108

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop