#### Abstract

Intelligent Transportation System (ITS) is more and more crucial in the modern transportation field, such as the applications of autonomous vehicles, dynamic traffic light sequences, and automatic road enforcement. As the upcoming fifth-generation mobile network (5G) is entering the deployment phase, the idea of cellular vehicle-to-everything (C-V2X) is proposed. The same 5G networks, coming to mobile phones, will also allow vehicles to communicate wirelessly with each other. Hence, 3.5 GHz, as the main sub-6 GHz band licensed in 5G, is focused in our study. In this paper, a comprehensive study on the channel characteristics for vehicle-to-infrastructure (V2I) link at 3.5 GHz frequency band is conducted through channel measurements and ray-tracing (RT) simulations. Firstly, the channel parameters of the V2I link are characterized based on the measurements, including power delay profile (PDP), path loss, root-mean-square (RMS) delay spread, and coherence bandwidth. Then, the measurement-validated RT simulator is utilized to conduct the simulations in order to supplement other channel parameters, in terms of the Ricean *K*-factor, angular spreads, the cross-correlations of abovementioned parameters, and the autocorrelation of each parameter itself. This work is aimed at helping the researchers understand the channel characteristics of the V2I link at 3.5 GHz and support the link-level and system level design for future vehicular communications of 5G.

#### 1. Introduction

Intelligent Transportation System (ITS) is more and more crucial in the modern transportation field, such as the applications of autonomous vehicles, dynamic traffic light sequences, and automatic road enforcement. Since the vehicles are at high speed, the real-time and fast exchange of dynamic information should be transmitted in a short time from vehicle-to-everything (V2X) [1]. The V2X provides wireless services for vehicle-to-vehicle (V2V), vehicle-to-infrastructure (V2I), and vehicle-to-pedestrian (V2P) communications. Such transferred information includes not only small data such as the speed, the location, and the directions of neighboring vehicles but also the big data such as video of surrounding environments and three-dimensional (3D) high-resolution maps. The transportation control system can use these applications for congestion avoidance, control performance improvement [2], overall traffic efficiency improvement [3], and cyberattack mitigation [4].

In order to provide the reliable and low latency for the V2X services, efforts have been made in recent years to develop V2X communications using IEEE 802.11p [5] and ITS-G5 [6]. Cellular vehicle-to-everything (C-V2X) communication is now regarded as another promising and feasible solution for the fifth-generation mobile networks- (5G-) enabled vehicular communications [7]. C-V2X will allow vehicles to communicate wirelessly with each other, with traffic signs, and with other roadside infrastructures by using the same 5G networks coming to mobile phones.

In [8], the authors propose an architecture that incorporates V2V communications into C-V2X and IEEE 802.11p-based the vehicular ad hoc networks (VANETs). For the frequency bands, 5.9 GHz band has been standardized in the IEEE 802.11p, while ITS-G5 is operating at 5 GHz frequency band. C-V2X will operate on the cellular networks using the same frequency bands as 5G new radio (NR). Frequency band for 5G NR [9] is being separated into two different ranges, including the sub-6 GHz frequency band [10] and millimeter wave (mmWave) frequency band [11]. The very high 5G frequency bands at mmWave, such as 24.25-27.5 GHz schemed in Europe, 24.25-27.5 GHz and 37-43.5 GHz schemed in China, 27.5-28.35 GHz and 37-40 GHz schemed in the United States, will allow the deployment of hotspots providing very high throughput thanks to the large available bandwidth for operators. On the other hand, the sub-6 GHz spectrum is less complicated in the development of infrastructure, deployment, and future network enhancements. Thus, it can be quickly cleared for the early deployment of 5G cellular networks across the globe. Europe has awarded the trial licenses at 3.4-3.8 GHz band [12]. China is ongoing at 2.5-2.6 GHz, 3.4-3.6 GHz, and 4.8-4.9 GHz frequency bands [13], while Korea is planning at 3.4-3.7 GHz. The United States schemes the frequency bands of 3.1-3.55 GHz and 3.7-4.2 GHz. Overall, the frequency from 2.5 GHz to 5.9 GHz is the main sub-6 GHz frequency band for V2X communications. Therefore, it is necessary to thoroughly investigate the V2X channel characteristics at this frequency band in different scenarios.

In the literature, the sub-6 GHz vehicular communication channels for 5G are characterized by academia and industry. Channel measurement is regarded as a valid approach to describe realistic propagation information. Many measurements have been taken in vehicular communications to characterize radio channel [14]. In vehicle communications, measurements are even more challenging due to high mobility scenarios and more dynamic network conditions [15]. The authors of [16] investigate channels under the conditions from Line-of-Sight (LoS) to non-Line-of-Sight (NLoS) in terms of delay spread, Doppler spread, and the RMS delay spreads at 5.9 GHz. The author of [17] present thorough characterizations for intracar propagation based on the practical measurement on-board an HSR at 2.4 and 5.8 GHz. By considering the traffic flows and the velocity of the vehicles, the wireless communications at 2.48 GHz under the NLoS condition are compared between the measurements and proposed time-variant temporal correlation function (TCF) model in [18]. The cross-polarized channel measurements are conducted in a small-cell street canyon scenario at 2.6 GHz [19]. It is found that the polarization does not influence the angles of arrival. In [20], channel measurements in an urban scenario are performed at 2.6 GHz using a virtual uniform linear array (ULA). In [21], the mobility measurements are conducted in an urban macrocellular (UMa) scenario to analyze the path loss at 2.54 GHz, 3.5 GHz, 4.9 GHz, and 5.4 GHz, respectively. Path loss, Ricean *K*-factor, and angular power spectrum are analyzed. The path loss of the measurements in an overpass scenario is fitted by two-ray model at 5 GHz [22]. The distributions of time-varying delay spread and Doppler spread follow a bimodal Gaussian mixture distribution in the street cross scenario at 5.5 GHz [23]. The authors of [24] survey the propagation characteristics at sub-6 GHz and mmWave bands, in which the preferred bands for initial deployment are 3.3-4.2 GHz and 26/28 GHz for 5G communications. The authors of [25] summarize the advanced channel measurements for vehicular communications at sub-6 GHz frequency band.

Channel models can be classified as stochastic and deterministic ones. The geometry-based stochastic channel models (GBSMs) with wide-sense stationary (WSS) assumption for vehicular channels have been widely accepted in the past studies [26, 27]. In [28], the GBSM combined with a two-ray model is proposed for nonisotropic multiple-input multiple-output (MIMO) Ricean fading channels in regular shaped environments and is validated for V2V channel at 5.9 GHz. Moreover, several measurements show that the WSS assumption validates only in a short distance interval and time interval [29]. For the nonstationary vehicular channel modeling, most works focus on the surrounding scatterers, e.g., vehicles, pedestrians [22], and traffic signs [30]. Deterministic channel models, such as ray-tracing (RT), can provide accurate channel information in propagation scenarios. The LoS, reflected, scattered, diffracted, and penetrated multipaths can be traced as rays, which contain the amplitude, angle, and delay information. In [31], the authors integrate the radar cross sections of the small-scale structures (e.g., lampposts and traffic signs) into RT simulator and conduct extensive RT simulations for V2V channel in urban and open space environments at 5 GHz. In [32], the authors propose that RT can be utilized to detect hidden obstacles in an external environment of a vehicle with the help of a light detection and ranging (LIDAR) sensor. In [33], the authors compare the RT simulations against V2V channel measurements using a channel sounder in an urban scenario at 5.9 GHz. The received power, delay spreads, and Doppler spreads are characterized in LoS and NLoS conditions, respectively.

Nevertheless, a comprehensive characterization of V2I channel based on cellular networks at 3.5 GHz frequency band is rarely found. Due to the limits of measurement equipments, propagation scenarios, and traveling speeds of the vehicles, a tendency of vehicular channel characterization is to combine the RT simulations. Thus, in this paper, we characterize the V2I channel at 3.5 GHz based on channel measurements and simulations with the following contributions:
(i)Dynamic channel measurements are performed in a suburban scenario at 3.5 GHz frequency with the bandwidth of 100 MHz, because 3.5 GHz is the common spectrum band in Europe [12], China [13], and other countries of the world for 5G. The transmitter (Tx) is fixed, while the moving receiver (Rx) is carried by a vehicle during the measurements. Double-cone omnidirectional antennas for both Tx and Rx are utilized in the measurements. The power delay profiles (PDPs), path loss (PL), shadow fading (SF), root-mean-square (RMS) delay spread (DS), and the coherence bandwidth are analyzed based on the measured results(ii)We reconstruct the suburban scenario and conduct the simulations utilizing the measurement-validated RT simulator with the same configurations as the measurement. The RT simulations breakthrough the limits of the measurement data and expand the measured two-dimensional (2D) channel into a 3D channel. The electromagnetic (EM) properties of relevant materials are calibrated. Then, the channel is characterized more comprehensively with the aid of Ricean *K*-factor (KF), angular spreads, and cross-polarization ratio (XPR). Combined with the channel parameters extracted from measurements, the cross-correlations of the abovementioned parameters are listed, and the autocorrelation of each parameter itself is calculated as well. The provided key parameters of this paper help to understand the V2I channel characteristics and to support the link-level and system-level design. These will enable the C-V2X communications for 5G

The remainder of this paper is organized as follows. The measurement systems and campaigns for the V2I channel are detailed in Section 2. The analysis of the measurement results is described in Section 3. In Section 4, the RT simulations are conducted in the reconstructed measurement scenario. Finally, the conclusions are drawn in Section 5.

#### 2. Channel Measurement Campaign

In order to measure the V2I channel, a dynamic frequency domain channel sounding method is applied in a suburban scenario. As shown in Figure 1, this measurement system is supported by a controlling laptop, a radio frequency (RF) signal source by NI PXIe-5673E module, a power amplifier, and an antenna at the Tx side. At the Rx side, it is composed of a controlling laptop, a wideband vector signal analyzer (VSA) by NI PXIe-5668 module, a disk array, and an antenna. The XHTF3311 global positioning system (GPS) and a rubidium clock are utilized to ensure the synchronization of the signal for Tx and Rx during the measurements. The carrier frequency is 3.5 GHz with a bandwidth of 100 MHz, and the number of frequency points is 821. The transmitting power is 40.77 dBm by the power amplifier. The transmit signal is a Zadoff-Chu (ZC) sequence and can provide an almost constant magnitude in the frequency domain [34].

The measurement campaigns are conducted in Haidian District, Beijing, China. We select the measurement area far from the main road of the city and only conduct the measurements in the evening so that by-passing traffics and pedestrians do not affect the measurements. We measure the LoS scenario, which exists mostly in V2I channels. As shown in Figures 2 and 3, the fixed Tx, with a height of 2.55 m, is at the intersection of the two roads. In the measurements, the Tx transmits signal sequences continuously. In order to emulate a typical vehicular environment, the two crossroads are selected in the measurements. The Rx is carried by a vehicle with a height of 1.88 m (as shown in Figure 3) forward or reversely following the track of the red arrows. Route 1 from north to south (approximately 210 m) represents the vehicle moving towards the fixed Tx, while Route 2 from west to east (approximately 340 m) represents the vehicle moving away from the Tx. Either Tx or Rx is a double-cone omnidirectional antenna with 0 dBi gain. The 3D radiation patterns and the 2D radiation patterns in the horizontal plane (H-plane) and elevation plane (E-plane) of utilized antennas are shown in Figure 4.

#### 3. Analysis of Measurement Results

In this section, the measured data are preprocessed in the frequency domain and time domain. The channel characteristics, including PDP, path loss, shadow fading, RMS delay spread, and coherence bandwidth, are calculated and analyzed for the V2I channel as follows.

##### 3.1. Preprocessing

The received signal in the frequency domain can be expressed as follows [34]: where and are the transmitted and received signals, respectively. is the distance between Tx and Rx, and is the operating frequency. and are the transfer functions of the Tx and Rx equipment, such as antennas, VSA, and cables, while is the transfer function of the wireless channel. In order to remove the influence of the transceivers, the reference signal measurements with the same setting as in the channel measurements are conducted in an anechoic chamber. The reference received signal can be expressed as follows: where is the free space transfer function.

Thus, the channel transfer function (CTF) can be calculated as follows:

Then, the channel impulse response (CIR) can be calculated by the inverse Fourier transform: where is the time delay of the received signal and is the effective number of the measured frequency points. A Hann window is used to suppress side lobes, and the sliding window is with the size of 20 wavelengths. A total of 4080 effective CIRs are extracted for the entire travel trajectory.

In addition, at the largest distance between Tx and Rx, the lowest average signal-to-noise-ratio (SNR) is about 20 dB, and at most of the measurement locations, the SNR is over 25 dB. Hence, the accurate estimation of channel parameters is possible in the following part.

##### 3.2. Power Delay Profile

The signal power on each multipath against their respective propagation delay is defined by the PDP, which is defined as:

A time-varying PDP of the measurement is shown in Figure 5. The calculated optimal threshold after filtering noise is -95.69 dBm, which is 6 dB higher than the mean value of the noise power and effectively separates the signal from noise. The dynamic range is larger than 40 dB. Then, the local maximums as the “peaks” higher than the noise threshold are identified as the “paths.” It can be seen clearly that in the propagation environment, there are always two energy-intensive paths. One is the LoS path, and the other is the strongly reflected path from the ground. Besides, there are rich multipaths in the whole scenario.

##### 3.3. Path Loss Exponent and Shadow Fading

The local wideband path gain can be calculated directly from the measured CTFs in the frequency domain as follows: where is the sampling frequency and is the number of the measured frequency points. As can be seen from Figure 2, the moving track of Route 1 is Rx towards the fixed Tx, while the moving track of Route 2 is Rx away from the fixed Tx. The environments and the length of these two roads are not similar. In Route 1, there are buildings along the two sides of the road, while only one side of Route 2 has buildings, and the other side is a fence with advertisements. Thus, two routes should be separated to analyze. To remove the effect of small scale fading, in this paper, we use a sliding window with the size of 20 wavelengths for averaging.

The log-distance path loss model is to fit the opposite number of measured path gain. The change of the path loss, along with the distance, is depicted by the path loss exponent. This exponent and shadow fading are extracted from the measured results by using the following expression [35]:
where is the path loss and is the distance between the Tx and Rx in . *d*_{0} is the reference distance; we choose the measured PL at as a reference in our study. is the media path loss at , and is the path loss exponent. By using the least-square criterion, the path loss exponent can be obtained. is the shadow fading, which can be expressed as a Gaussian variable with zero mean value and a standard deviation of . The path loss exponent and the shadow fading standard deviations based on the measurements are summarized in Table 1. The fact of higher path loss exponent and higher for Route 1 indicates a more complicated environmental condition of the Route 1. There are a large number of trees, two rows of iron fences, and concrete walls on the side of Route 1, which cause the scattering during the wireless propagation.

##### 3.4. RMS Delay Spread and Coherence Bandwidth

The RMS delay spread is one of the key parameters for a multipath channel. It is defined in [36], namely, the square root of the second central moment of the PDP as follows:
where denotes the RMS delay spread and and denote the power and the excess delay of the -th multipaths, respectively. The distribution of RMS delay spreads for the whole trajectory is following a Gaussian distribution. The mean value () of the measurement routes is 670.34 ns with the standard deviation () 355.69 ns. It is in line with the outdoor suburban environment, where the RMS delay spread value is 0.2-2.0 *μ*s. In [23], the RMS delay spreads around 47 ns and 75 ns are measured in urban and tunnel scenario, respectively. However, sometimes, although the LoS component exists in the channel, several significant reflected multipaths and scattered multipaths can still result in a large delay spread, such as in our study.

The inversely proportional relation between RMS delay spread and the coherence bandwidth can be expressed: where is a constant that depends on the environment or/and on how is defined. can be defined as the bandwidth at which the complex correlation function has a value of 0.5, 0.9, or even 0.95. For example, means the coherence bandwidth can be specifically denoted as with subscript 0.5 as the threshold [37].

The frequency correlation function of the channel can be obtained by setting a reasonable threshold [38]: where is the CIR at a specific time delay . The statistical estimation of the mean value and the standard deviation of coherence bandwidths can be obtained according to the threshold of 0.5 in our study. The values of RMS delay spread and coherence bandwidths are summarized in Table 2. The coherence bandwidth is narrower when using the onmi-directional antennas. It is consistent with the RMS delay spread results and the inverse relation between these two parameters. As happens with the RMS delay spread, the multipath propagation causes the coherence bandwidth.

Moreover, the study on Doppler spectrum of these measurement campaigns has been published in [39]. It can be concluded that the faster the vehicle (Rx) moves, the larger the Doppler spreads. The most significant scatterers during the measurements are trees and billboards, whereas other cars parked on the roadside do not significantly contribute to the multipath propagation.

#### 4. RT Simulation for V2I Channels at 3.5 GHZ

In order to physically interpret the measurement results, a self-developed dynamic RT simulator of Beijing Jiaotong University is used in this study. Supported by high-performance computing (HPC), this RT simulator is now furnished to an HPC cloud-based platform (CloudRT), which contains 1600 CPUs and 10 GPUs. Parallel processing can be implemented by a compute node or monitored by a management node. The first premise of accurate ray-tracing simulation is a thorough and factual description of the scenarios. The finer the scenarios modeled, the more accurate the corresponding prediction is. Also, simulation time is proportional to the accuracy. The simulation time for a single snapshot in this article is about 10 seconds, and the total simulation time depends on the number of nodes selected for parallel computing. More details of this platform can be found in [40] and the website http://raytracer.cloud.

As shown in Figure 6, the 3D model of the suburban scenario is reconstructed through SketchUp [41] and OpenStreet Map [42]. The asphalt road of Route 2 is with a width of 14 m and with a length of 340 m. The height of the office buildings on the one side is about 15-20 m. In order to simulate the channel for 100 MHz bandwidth, we use the subband RT approach. The CIR for one snapshot is generated by a number of subbands at multiple center frequency points [43]. Based on the measurement configurations in the same scenario, 821 frequency points are determined. Since the relevant materials are not frequency sensitive, the constant EM properties of them are used for all subbands at 3.5 GHz. The LoS ray, reflected rays (up to third order), scattered rays, and diffracted rays are traced in the simulations for the trade-off between simulation accuracy and computational complexity. The directive scattering model [44] and the uniform theory of diffraction (UTD) model [45] are used in the simulations. In order to get the pure propagation channel without the influence of certain antenna pattern, the antenna of Tx or Rx is omnidirectional antenna, which is in line with the measurements. In addition, in order to get the fully polarimetric information of the channel, the antennas of Tx and Rx are with four cases of polarization: vertical-vertical (VV) polarization, vertical-horizontal (VH) polarization, horizontal-vertical (HV) polarization, and horizontal-horizontal (HH) polarization in the extensive RT simulations. Table 3 summarizes the simulation configurations.

##### 4.1. Validation and Calibration of RT Simulator

Several snapshots are stochastically chosen to validate the geometry and calibration of the EM properties of the suburban scenario, where the validation and calibration processes are according to our previous work [46]. Figure 7, as an example, shows the comparison of path loss between the measurement and the validated and calibrated RT simulation. The antenna influence can be neglected with the omnidirectional antennas both at Tx and Rx sides. Compared with the measurement results, the difference of path loss exponent between simulation and measurement is 0.12, with the relative error 6.4%. The difference between the shadow fading is 0.29 dB, and the relative error is 5.6%, which proves the accuracy of the RT simulator in channel simulation. The calibrated EM properties are summarized in Table 4, which will be used in the following extensive simulations.

##### 4.2. Channel Characterization Based on Extensive RT Simulations

The extensive RT simulations are done for more comprehensive channel characterization. Route 2 selected as an example, the key parameters can be supplemented to fully describe for the V2I channel, in terms of the Ricean *K*-factor, XPR, and four angular spreads (namely, azimuth angular spread of arrival (ASA), azimuth angular spread of departure (ASD), elevation angular spread of arrival (ESA), and elevation spread of departure (ESD)). Combined with the channel parameters extracted from the measurements, the cross-correlations of the abovementioned parameters are listed, and the autocorrelation of each parameter itself is calculated as well.

###### 4.2.1. Ricean *K*-Factor

The Ricean *K*-factor is defined as the ratio of the signal power of the LoS path to the signal power of other multipath components, which can be expressed as:
where is the energy of all the paths of the signal during transmission and is the energy of the LoS path. As shown in Figure 8, most values of the *K*-factor are larger than 0 dB. When the probability is 0.5, the corresponding value of KF is 5.98 dB; when it reaches 0.9, the value is 18.16 dB. Such high values indicate that the LoS component plays the leading role in the propagation scenario, while the multipaths reflected from the surfaces of the buildings and the ground are relatively weak.

###### 4.2.2. Angular Spreads

The angular spread is the distribution of the angle of arrival/departure of each multipath in 3D environments for the Tx and Rx, respectively. Each multipath is assigned an azimuth angle in the horizontal plane and an elevation angle in the vertical plane. Hence, four values of angular spreads can be calculated as follows: where denotes the angular spread (AS), denotes the power of the -th multipath, and is defined by: where is the azimuth angle of arrival (AoA)/azimuth angle of departure (AoD)/elevation angle of arrival (EoA)/elevation angle of departure (EoD) of the -th ray. is:

Figure 9 shows the CDFs of the four angular spreads, including ASA, ESA, ASD, and ESD, respectively. Most elevation angular spreads are below 10°, which is much smaller than azimuth angular spreads in the horizontal plane. In general, the scattering effect in the horizontal direction is stronger than it in the vertical direction, because of the multipath components reflected from the infrastructures along the roadside. We further analyze the variation of the angular spreads during the vehicle moving trajectory, as shown in Figure 10. The values of ESA and ESD are more significant than those of ASA and ASD when the distance between Tx and Rx is within 10 m. This observation implies that multipaths reflected from the ground in the vertical direction are stronger than those in the horizontal direction at the beginning. As the distance increases, the values of ASA and ASD are larger than those of ESA and ESD. Concerning the ASA, it increases sharply to 173° due to energetic scattered paths from the infrastructures far away from the Tx but close to Rx.

**(a) ASA**

**(b) ESA**

**(c) ASD**

**(d) ESD**

###### 4.2.3. XPR

The XPR is to describe the variety of the polarization for the multipath components and is calculated by its definition as follows: where is copolarization defined by the received power when Tx and Rx antennas polarize in the same manner. Conversely, is the power received in vertical (horizontal) polarization and transmitted in horizontal (vertical) polarization. As shown in Figure 11, the high values imply little depolarization in such a suburban scenario. These values are much larger than the higher frequency band in similar scenarios, such as at 28 GHz [47] and at THz band [48]. Generally speaking, the linear polarization of the antennas is sufficient in such a scenario.

###### 4.2.4. Autocorrelation and Decorrelation Distance

The autocorrelation function of a parameter is to describe how fast a certain variable changes with the distance interval (or time delay ), which can be calculated as follows: where is the expected value operator. Correspondingly, the decorrelation distance in (m) is expressed as: where is the threshold.

The larger the threshold is set, the stricter the correlation requirements are, and the smaller the decorrelation distance is. The threshold is generally a practical value. The common value of 0.9 is set in our study. Table 5 summarizes the decorrelation distance of each parameter. It can be preconcluded that the shadow fading is not self-correlated. The KF and DS are highly self-related. The value of the decorrelation distances of the parameters mentioned above is in line with the values provided by 3GPP 38.901 [49].

###### 4.2.5. Cross-Correlation

We consider a wideband V2I channel, study the key parameters of the channels across 100 MHz at 3.5 GHz frequency band, and try to capture the channel correlation across different parameters. The cross-correlation coefficient of the two variables is a measure of their linear dependence. If each variable has scalar observations, then the Pearson cross-correlation coefficient is defined as: where and represent two different aforementioned channel parameters and and denote the mean value and the standard deviation of corresponding parameters, respectively.

Combined with the measurements and simulation results, the cross-correlations of the key parameters are summarized in Table 6. It can be found that shadow fading is not cross-correlated with other key parameters since the correlations of it with other key parameters are approximately equal to 0. Remarkably, KF is negatively correlated to the angular spreads, since with a large value of KF, a significant component comes from a single direction leading to smaller angular spreads. In the angular domain, the parameters show relatively strong correlations between each other. If the variables of interest are highly correlated, it can be used to make accurate predictions.

###### 4.2.6. Summary of the Channel Parameters by Measurements and Simulations

The V2I channel is comprehensively characterized based on the measurements and simulation together. The group of channel parameters can be summarized in Table 7. The validation results of path loss indicate that RT can practically reflect the target scenarios. Besides, the RT simulations can break the limits of the measurement data. Reliable channel information can be obtained with similar materials. Moreover, other important channel information, such as angles, which are not captured from the measurements easily, can be given by RT simulations.

#### 5. Conclusion

In this paper, we comprehensively characterize the V2I channel through the combination of the measurements and simulations in a suburban scenario at 3.5 GHz band for 5G communications. To begin with, the channel measurements are conducted in a typical vehicular scenario. The track of moving vehicle is separated into two routes, where the environment of Route 1 is more complicated than Route 2. The time-varying PDP, path loss, RMS delay spread, and the coherence bandwidth are extracted from the channel measurements. From the time-variant power delay profiles, it can be clearly seen that two significant paths always exist during the measurements. It indicates that the LoS vehicular channel measurement generally follows the classic two-ray models. One is the LoS path, and the other is the reflected path from the ground. The path loss exponent and the are extracted to describe the path loss and shadow fading, respectively. Because of the more complicated environment of Route 1, the values of and are larger than these of Route 2. The coherence bandwidth is to follow the inversely proportional relation to RMS delay spreads. As the threshold is set 0.5, the coherence bandwidth then can be calculated and analyzed.

Moreover, the RT simulations are conducted in order to supply other key parameters. The high-performance self-developed RT simulator of Beijing Jiaotong University is utilized in our study. The RT simulations breakthrough the limits of the measurement data and expand the measured 2D channel into a 3D channel. After validated and calibrated by the measured PDPs and compared with the path loss, the RT simulator is utilized to do extensive simulations in a reconstructed suburban scenario. The large value of Ricean *K*-factor implies the LoS leading a significant role in the propagation. Four cases of the vertical/horizontal polarization of antennas for Tx and Rx are schemed to characterize the fully polarimetric information of the channel. The analysis of the XPR shows little depolarization in such a LoS condition. The angular spreads (ASA, ESA, ASD, and ESD) are analyzed. When the distance between Tx and Rx is within 10 m, the multipaths reflected from the ground are stronger than the multipaths reflected/scattered from the infrastructures, iron fences, and tress. As the distance increases, the azimuth angular spreads become larger than elevation angular spreads. Moreover, the decorrelation distance of each parameter and the cross-correlations of the abovementioned parameters are supplied. SF is neither self-correlated nor cross-correlated with other channel parameters. KF is negatively correlated to the four angular spreads, since a significant component coming from a single direction leads to a large value of KF but small angular spreads.

In summary, the group of channel parameters is thoroughly analyzed, and the V2I channel is comprehensively characterized based on the measurements and simulations. This study on channel characteristics at 3.5 GHz can be applied to the link-level and system-level design for the C-V2X communications of 5G.

#### Data Availability

The raw data required to reproduce these findings cannot be shared at this time as the data also forms part of an ongoing study.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work is supported by the Key-Area Research and Development Program of Guangdong Province, China (2019B010157001), the State Key Laboratory of Rail Traffic Control and Safety (Contract No. RCS2020ZZ005), and the NSFC under Grants 61771036, 61911530260, 61901029, and 61725101.