A Simplified Multipath Component Modeling Approach for High-Speed Train Channel Based on Ray Tracing

High-speed train (HST) communications atmillimeter-wave (mmWave) band have received a lot of attention due to their numerous high-data-rate applications enabling smart railmobility. Accurate and effective channelmodels are always critical to theHST system design, assessment, and optimization. A distinctive feature of the mmWave HST channel is that it is rapidly time-varying. To depict this feature, a geometry-based multipath model is established for the dominant multipath behavior in delay and Doppler domains. Because of insufficient mmWave HST channel measurement with high mobility, the model is developed by a measurementvalidated ray tracing (RT) simulator. Different from conventional models, the temporal evolution of dominant multipath behavior is characterized by its geometry factor that represents the geometrical relationship of the dominant multipath component (MPC) to HST environment. Actually, during each dominant multipath lifetime, its geometry factor is fixed. To statistically model the geometry factor and its lifetime, the dominant MPCs are extracted within each local wide-sense stationary (WSS) region and are tracked over different WSS regions to identify its “birth” and “death” regions. Then, complex attenuation of dominant MPC is jointly modeled by its delay and Doppler shift both which are derived from its geometry factor. Finally, the model implementation is verified by comparison between RT simulated and modeled delay and Doppler spreads.


Introduction
Characterized by high mobility, the high-speed train (HST) scenario is expected to be a typical scenario in the fifthgeneration (5G) wireless communication systems [1].To enable intelligent transportation system, various high-datarate services are required in this scenario [2][3][4], for example, high definition (HD) video surveillance, and onboard realtime connectivity.For satisfying these services, applying advanced technologies in HST communications, including millimeter-wave (mmWave), mobile relay station (MRS), and smart antennas, has been studied in the literature [5][6][7][8].Beforehand, to design these advanced technologies and evaluate their performance, accurate channel models are necessary.However, existing channel models at mmWave band mostly focus on indoor and static scenarios [9][10][11][12][13], without considering the rapidly time-varying HST channels between transmitter () and MRS in open space scenario.
Due to these facts, mmWave HST channel modeling has been underresearched in the literatures.
For modeling approaches, the current HST channel models can be roughly classified into three categories [5]: deterministic, stochastic, and geometry-based stochastic modeling.
Deterministic modeling is mainly based on geometrical optics or, more generally, on solving Maxwell's equations.Scenario reconstruction is the first and key step, which gives detailed geometries, specific locations of  and receiver (), and so on.Then, combining with antenna configuration the complex attenuation, delay, angle, and polarization information of the propagation waves are obtained through intensive electromagnetic simulation.Therefore, the deterministic models are approximately exact.For HST channels deterministic modeling is proposed in [14][15][16][17][18].In [14], the temporal evolutions of delay spread and Doppler spread at 5.2 GHz in orthogonal frequency division multiplexing (OFDM) HST communication system are evaluated based on ray tracing (RT).Similarly, [15,16] using RT analyze HST channel characteristics, for example, path loss (PL), delay spread (DS), and angle spread (AS), at 930 MHz and 5.25 GHz, respectively.As for these studies, time-varying characteristics in HST channels are analyzed in the non-mmWave band.Moreover, in [17] a linear HST environment is conducted and RT simulation at 60 GHz for analysis of DS and AS changing over time is performed.In [18], adopting RT at 60 GHz as well, six HST scenario modules are modeled and the DS for different modules are compared.However, in [17,18] the validation by measurement in the HST deterministic models at mmWave band, which is crucial in channel modeling, is not performed.Furthermore, such modeling approach represented by RT fails in the situation as it is too complicated to reconstruct the exact scenario.In addition, the extensive computational resources are needed, which is not suitable for communication link-level and system-level simulations.Thus, stochastic modeling which is void of exact geometrical assumption and is computationally less demanding has been proposed.
Stochastic modeling provides the statistical manner and conditional dependencies between different channel parameters.Although such modeling approach gives very little insight in an instantaneous channel structure, it is useful in presenting the general characteristic of a typical scenario [19].From this point of view the stochastic modeling is practically useful.The existing stochastic models for HST channels mainly focus on the narrowband [20][21][22][23][24].For example, at 930 MHz, the PL, shadow fading, Rician -factor, and decorrelation distance are stochastically modeled based on various HST measurements [20][21][22][23].Their time-varying characteristics are described by Markov-chains [24].Moreover, for wideband HST channels, the rough position-based tapped delay line channel models with specified Doppler shifts for open space scenario and hilly scenario are proposed in [25,26], respectively.Nevertheless, comparing with deterministic modeling the stochastic models cannot keep the spatial consistency which is crucial in time-varying HST channel modeling [27].Furthermore, the nonstationarity caused by high mobility of MRS is not characterized in these HST stochastic models.
To address the nonstationarity in time-varying channel, Bello proposed the wide-sense stationary (WSS) concept [28].Under the WSS assumption the channel can be stationary over a limited time region and therein the stochastic modeling becomes physically meaningful.Particularly, the power delay profile (PDP) in delay domain and Doppler power spectral density (DSD) in Doppler domain completely describe the second-order statistics of time-varying channels [29,30].
The third modeling approach, the geometry-based stochastic channel model (GSCM), is useful in nonstationary channel [31].This is the type of modeling approach we aim for in this paper.Within the GSCMs specific locations of , , and scatterers are incorporated, where the distribution for scatterers is predefined by the statistical law.The standard channel models, like rural macrocell (RMa) scenario in the WINNER II channel model [32] and moving networks scenario in the IMT-A channel model [33], have considered the HST scenario.However, they are for the non-mmWave band.In [34], the narrowband GSCM is modeled and the scattering components are represented by the measurement-defining cluster.In [35], the regular shaped GSCM is developed by considering time-varying model parameters, that is, angle of departure (AoD) and angle of arrival (AoA).It provides the insights into the relationship of channel characteristics and certain geometry.However, due to oversimplified propagation environment the parametrization work from measurements would be hard [36].In [37] the GSCM with random-cluster is proposed in delay and Doppler domains at 2.14 GHz.And yet for the limited bandwidth, multipath components (MPCs) are unresolvable and only a LoS cluster is modeled according to LoS geometry.While these GSCMs might not be scalable to the mmWave HST channels, they provide some geometrical insights for the HST channel modeling.
In this paper, we propose a geometry-based dominant multipath model for the mmWave HST channel.It fully depicts the temporal evolution of dominant multipath in the delay and Doppler domains.The modeling approach is simple and effective.Without using intensive RT or numerical simulations, the proposed modeling approach is based on the geometry factor extracted from the propagation environment.Through the right triangle constructed by the geometry factor,  location, and instantaneous location of MRS, the time-varying delay and Doppler shift can be derived.The amplitude of dominant multipath is modeled by the delay and Doppler shift which is presented first.The amplitude modeling implies the correlation between the delay and angle domains in the mmWave HST channel.
The proposed modeling approach shows insight into the relationship between the delay and Doppler shift for dominant MPCs in the mmWave HST channels.A similar modeling approach is introduced in [19,38], where only single bounced rays were included.The proposed approach is capable of modeling first-order reflected paths (FRPs) and second-order reflected paths (SRPs), and the reflected paths strictly obey specular reflection law.Moreover, the number of newly dominant paths, their lifetimes, geometry factors (for the delay and Doppler shift calculation), and complex attenuation are statistically modeled and parameterized.A detailed implementation of the proposed model is presented, which is verified by comparison with measurement-validated RT.This simple and effective modeling approach is practical for system-level simulation and communication network evaluation.
The remainder of this paper is structured as follows.In Section 2, the RT simulator is validated through the measured data obtained in the mmWave HST channel and then the temporal evolution of PDP and DSD simulated by the validated RT in the mmWave HST scenario is discussed.Section 3 describes the structure of the proposed model, and the dominant paths are detected and tracked using geometry factors.Section 4 shows the model parameters.Section 5 presents the model implementation and validation.This paper concludes in Section 6.

Ray Tracing Validation and Simulation
The typical mmWave HST scenario is shown is Figure 1, whose construct is proposed in the 3rd Generation Partnership Project (3GPP) [6].The remote radio head (RRH) is placed 5 m from the track.The MRS is deployed on the top of the train.In addition, with the longitudinal separation of 580 m three RRHs are organized as a group connected to the baseband processing unit (BBU), which is a new cellular architecture for the mmWave HST communication system.Due to deficient mmWave HST channel measurement with high mobility, development of the modeling approach is based on RT simulator.Therefore, first we will validate the RT simulator by the channel measurement which is performed in the mmWave HST scenario.In detail, the measurement campaign conducted in Japan was performed at 93.2 GHz and with 2 GHz bandwidth [39].The frequency band was chosen in this measurement for two primary reasons.One is to consider the coexistence with active and passive services in the HST scenarios [40], and the other is for the high data transmission rate.The authors in [41] demonstrate that the attenuation of the 90 GHz band millimeter-wave in the atmosphere is small compared with millimeter-waves of the other frequency bands.The MRS moved along the track and the RRH-MRS distance ranged from 20 m to 350 m.  and  antennas were installed on RRH and MRS with the same half power beam width (HPBW), respectively.The antenna settings are summarized in Table 1.Note in RT simulation the antenna pattern is inserted by the Jones vector at all the angles of arrival or departure.These Jones vectors cover the gain and the polarization of the antennas.The same operation was done and reported in [15].Moreover, the front view of typical structure in the measured scenario is shown in Figure 2.

Ray Tracing Validation.
RT simulation was executed in the measured scenario with settings listed in Table 1.Specifically, the LoS path, reflected paths (up to the third order), and scattering rays are simulated.Since there is no  obstacle between  and , the diffraction mechanism is not included in the RT simulation.Since the mmWave HST channel was not measured dynamically, the measured channel characteristic in Doppler domain cannot be obtained.In addition, only the received power has been presented in [39] with no information about the delay domain.Due to the limitations for the measurement campaign, just the large-scale and small-scale parameters between measured and simulated channel are compared for validation.We separate large-scale fading and small-scale fading by a sliding/overlapped window (size of 20 wavelengths) with an interval of ten wavelengths [25,42].

Large-Scale Fading
Validation.The decibel-valued path loss (PL) model is given by a linear Least Squares (LS) regression of the large-scale fading where  is the distance between the  and  in meters and PL( 0 ) is the frequency-dependent free space path loss at reference distance  0 , generally  0 = 1. refers to the average path loss exponents (PLE) over the observed distance.Moreover, for modeling the shadow fading,   denotes the zero-mean lognormal random value with standard deviation , that is, shadow factor (SF).If  0 and frequency are known, PL( 0 ) in ( 1) is fixed.Hence, for large-scale fading validation, we compare the two parameters, that is, PLE and SF.
The measured and simulated large-scale fading are shown in Figure 3, where their respective PLE and SF are given and the free space path loss at 93.2 GHz is added for comparison.The same PLE, that is, 1.6, is estimated from the two channels.The SF for the measured channel is 2.7 dB, which is just slightly higher than the simulated SF of 2.0 dB.The good match of the PLE and SF indicates reasonable agreement in large-scale fading.In addition, due to the superposition of MPCs, the PLE in the HST scenario is lower than the PLE from free space.

Small-Scale Fading Validation.
The empirical cumulative distribution functions (CDFs) of measured and RT simulated amplitudes for small-scale fading are shown in Figure 4.There is a good match between measured and RT simulated statistical properties in small-scale fading.Moreover, the mean vales of the measured and RT simulated amplitudes are just slightly different; that is,  meas = 0.98 and  simu = 0.97.Due to existent LoS path, the Rician distribution is adopted to fit the distribution.Figure 4 shows that the Rician distribution can accurately characterize the amplitude of small-scale fading.The Rician fitted curve is derived from the measurement data.
The validation confirms that the RT simulator is able to accurately reproduce mmWave channel characteristics in HST scenario.It is worth noting that this scenario is characterized by its effective reflectors, most of which are trackside walls.Actually, the trackside walls are not necessarily from the barriers next to the track.Assuming that in the scenario the main effective reflectors are perpendicular to the ground, such as buildings and tunnel walls, the RT simulator should work adequately.Due to the measurement-validated RT we have the opportunity to extend the limited HST measurement scenario to a general mmWave HST scenario.

Ray Tracing Simulation.
The simulated scenario is a typical mmWave HST environment shown in Figure 1.For the general RT simulation, the antenna heights and patterns in measured scenario have been changed and different antenna heights are set, that is, 3.2 m at  and 3 m at .The antenna heights for the RT simulation are raised compared with the measurement settings and exceed the heights of the trackside barriers as shown in Figure 2. Specifically, to raise the antenna heights is aimed at breaking the semiclosed space constructed by trackside barriers.Then, the impact of trackside buildings (blue blocks in Figure 1) can be increased, just like in the general propagation environment.To characterize the propagation channel in the mmWave HST scenario, the omnidirectional, vertically polarized pattern is used at  and , instead of the directional patterns.The movement range of MRS is from RRH2 to RRH1 with length of 580 m, shown in Figure 1.The MRS speed is set to 500 km/h that is the top speed discussed in 3GPP [6].In addition, the parameter setting for the RT simulator is the same as above.
In this paper, we focus on the channel characteristics in delay and Doppler domains which are crucial in frame structure design (time domain) [43] and subcarrier spacing selection (frequency domain) [44].Particularly for Doppler analysis, the sampling interval should be less than halfwavelength, that is, 1.6 mm (carrier frequency is 93.2 GHz), which indicates the huge data volume for RT simulation.
During the MRS movement 720747 snapshots are simulated with a snapshot repetition time of 5.79 s.Note that, for reducing the volume of simulated data, only the IRs in time and frequency domains for each snapshot are stored and the multipath property, that is, rays absolute delay, angle, electrical field, and so on, has not been exported from RT.
The simulated data are stored in the same format as the data recorded in the channel measurement campaign, and then needed channel characteristics are extracted from the IRs, which makes the proposed approach general.
Analyses of the RT simulated data showed that the observed fading process in mmWave HST channel is nonstationary.Due to the propagation environment changing with a finite rate, the nonstationarity can be overcome by consecutive local WSS regions [45], that is, WSS time windows Δ.The local scattering function (LSF) introduced in [46] is used to determine the size of Δ.To ensure the stationarity analysis throughout, the minimum stationary time interval 3.6 ms is chosen for Δ.Therefore, within each WSS time window the statistics of delay and Doppler domains are locally stationary and we can average the local PDP and DSD within each WSS time window as introduced in [47].
For the mmWave HST channel the temporal evolutions of DSD and PDP are shown in Figure 5. Within each WSS time window, we normalize the instantaneous DSD and PDP to their own maximum values and set the cutoff threshold of 50 dB below the strongest peak.The continuous Doppler and delay trajectories of the time-varying LoS path can be easily distinguished in Figures 5(a Since the continuous or discontinuous trajectories seem to follow certain curves, the underlying geometry for them will be found hereafter.Moreover, it is obvious that the multipath resolution is higher in Doppler domain than in delay domain.Thus, for clarity, the following illustrations are mostly in Doppler domain.

Geometry-Based Dominant Multipath
Channel Model and Data Processing

Geometry-Based Dominant Multipath Channel Model.
The channel model is based on the channel spreading function [37].At time   , the channel spreading function is expressed as the superposition of (  ) paths where (⋅) is the Dirac delta function and   (  ),   (  ), and ]  (  ) represent the complex attenuation, delay, and Doppler shift of the th path observed at time   , respectively.For the HST channels, the behavior of each path, that is, complex attenuation, delay, and Doppler shift, changes over time.
Particularly, due to the high-speed MRS, the MPCs, except LoS path, frequently appear in the "birth" and "death" process, shown in Figure 5.For convenience, we define three sets of the path in (2) as follows [45]: (i) L → is the set of dominant paths that are firstly detected at time   , whose path index is denoted by  → .
(ii) L → is the set of dominant paths that are firstly detected at time   (∀ = 1, 2, . . .,  − 1) and still exist at time   , whose path index is denoted by  → .
(iii) L  is the union set of the above two sets, which contains all the dominant paths that exist at time   .
By detecting and tracking dominant paths in time-varying PDP and DSD, L → and L → can be constructed and the dominant multipath channel model can be reconstructed at time   .The dominant multipath number at time   can be expressed as where ( → ) represents the number of dominant paths being newly detected at time   , that is, the number of elements in L → ; the last term of (3) refers to the dominant multipath number which was detected before   and still exists at time   .Therefore, during the MRS movement, the proposed channel model is able to describe the birth and death process of dominant MPCs.Furthermore, by modeling the temporal evolution of dominant multipath behaviors in (2) during their lifetime, the proposed model can generate the exactly time-varying HST channels.Note, when a dominant path vanishes, its lifetime ends.This assumption simplifies the channel modeling and does not significantly affect the model accuracy [48].
For time-varying HST channels, parameters in (2) are generated by random geometry factors.To statistically model the geometry factor, a WSS time window should be used in a way that the statistics of geometry factor can be characterized within each WSS time window [49][50][51] and thereby channel modeling becomes physically meaningful.The size of WSS time window has to be defined previously.As shown in Figure 5, dominant MPCs exist in multiple consecutive WSS windows.Therefore, the lifetime of dominant path can be characterized by the WSS window number.In the proposed channel model, we use two distinct timescales to describe each path: (i) Large-scale: the geometry factor is indexed by the WSS window index, which indicates that the propagation environment remains changeless within a WSS time window.
(ii) Small-scale: the delay, Doppler shift, and complex attenuation of dominant path are derived at each time instant using the geometry factor.
To summarize, we statistically model the following parameters: (i) The size of dominant multipath lifetime in terms of WSS windows (ii) The number of dominant path being newly detected in each WSS time window (iii) The geometry factor for the newly dominant path in each WSS time window.
The developed dominant multipath model contains LoS path, FRPs, and SRPs.In practical terms, the channel model includes FRPs, SRPs, and third-or higher-order reflected paths.However, during MRS movement the energy ratio of third-or higher-order reflected paths is generally lower than 2% as shown in Figure 6.At each time instant the energy ratio is calculated as the sum power of certain order reflected path divided by the total power of all existent reflected paths.Thus, in the proposed channel model the third-or higher-order reflected path contributions are ignored.Note that this figure is obtained from the downsampling RT simulation with the output of multipath properties.The downsampling interval is 3.6 ms and the corresponding distance sampling interval is 0.5 m (the speed of 500 km/h).Considering that the distances among the RRH, physical objects, and the MRS are much larger than 0.5 m, it is reasonable to assume that during the interval no reflected path slipped away.Moreover, due to lower power contribution of diffuse multipath components (DMCs) at mmWave band DMCs are not modeled [52][53][54].To sum up, for the dominant multipath channel model, modeling LoS path, FRP, and SRP is sufficient.

Data Processing.
Generally, in HST channel the FRPs and SRPs are reflected from the trackside walls whose normal is perpendicular to the track.As illustrated in Figure 7 the potentially geometrical characteristic is fully explored for the dominant MPCs, and the geometry factor is explicitly presented.For each dominant multipath type there is a specific right triangle, where the length of the right-angle side (which is parallel to -axis) is defined as the geometry factor, that is, Δ LoS , Δ FRP , and Δ SRP .With the moving MRS (  ), the lengths of the other two sides are time-varying.
It is worth noting that for LoS path the hypotenuse is equal to the path distance, while, for the FRP and SRP, the length of the hypotenuse is, according to the specular reflection law, also equal to the path distance.Hence, the multipath delays  LoS (  ),  FRP (  ), and  SRP (  ) are able to be derived from the right triangles.Additionally, the respective angle between the path arrival at MRS and the MRS moving direction refers to  LoS (  ),  FRP (  ), and  SRP (  ), by which the Doppler shifts for dominant MPCs can be estimated.Therefore, through the right triangles, defined mainly by the geometry factors, the delay and Doppler shift of dominant MPCs can be derived at time   .
Based on the analysis, the dominant multipath model can be fully characterized by the geometry factor.Moreover, by corresponding geometry factor, the temporal evolution of multipath delay and Doppler shift can be derived.Hence, we detect and track dominant MPCs by geometry factor and its geometrical relationship.

Dominant Multipath Detection.
Before detecting the information of each dominant path (e.g.,   ,   , and ]  ), the analytical expressions for the delay and Doppler shift involving the geometry factor have to be clarified.For uniformity and convenience, the subscript indicating propagation mechanism in Figure 7 is replaced with the path index , where  = [1, . . ., (  )].Note that the physical significance where V denotes the MRS speed 500 km/h.As the hypotenuse length is equal to the travel distance of dominant path, the delay   (  ) is derived as where  refers to the speed of light and Δ  denotes the geometry factor for th path.Furthermore, the Doppler shift ]  (  ) can be given through   (  ) where   represents the carrier frequency of 93.2 GHz.Due to little height difference between  and  the elevation angle is not considered.The analytical form of cos( Afterwards, by involving the geometry factor, the dominant path is detected in delay and Doppler domains.Benefitting from the dilation and erosion mathematical morphology of binary images, the trajectories in Figure 5 are refined as being more clear, smooth, and coherent, as shown in Figure 8, which is beneficial for peak search.Specifically, for the image processing the time-varying DSD and PDP in Figure 5 are firstly transformed to binary image with threshold 25 dB [55].The pixel with power 25 dB lower than the maximum power is set to 0; otherwise, the pixel is set to 1.In the th WSS time window including time instant   , by identifying the peaks in the refined DSD and PDP, Doppler shift and delay as well as the associated power are obtained.Each piece of peak information is stored in a peak set []  (  ), τ (  ), p (  )].
Subsequently, using the geometrical relation determines whether the detected peak belongs to dominant MPCs.By combining ( 5), (6), and (7), Δ in (4) can be rearranged as To pick the dominant from the detected sets, the detected peak set in (8) has to be incorporated and the threshold  has to be used where | ⋅ | denotes absolute value and the value Δ(  ) is obtained from (4), which is according to time-varying position of the MRS, the truth value.If the inequality is satisfied, the peak set is stored in L  for dominant multipath modeling; otherwise, the peak set is discarded.In this analysis, we selected  = 0.9 m and the choices have to be done on a rather empirical basis.Afterwards, the geometry factor can be calculated as After the detecting, for the th path within th WSS time window, there is a set [Δ ŷ (  ), ] (  ), τ (  ), p (  )] stored in extended set L  .For each WSS time window, the detection process will be repeated.According to the detecting, the marginal values of geometry factors are obtained, that is, the minimum and maximum values.Figure 8 shows two Doppler trajectories for marginal values (derived by ( 6)) and the area confined by the two which includes all the Doppler trajectories of dominant MPCs.Moreover, the example trajectory (dash line) in Figure 8 is generated by the geometry factor observed in multiple WSS time windows.Its trajectory is discontinuous with birth and death process several times, which indicates the interruptedly distributed trackside walls.Furthermore, due to applying mathematical morphology, the temporal evolution of DSD in Figure 8 is clearer and with less thin protrusions compared to the DSD in Figure 5(a).

Dominant Multipath Tracking.
To capture the dominant multipath evolution within individual lifetimes, the tracking based on the geometry factor is performed.In two successive WSS windows, if the same geometry factor appears, its associated paths are considered to be the same path.This can be used to track the dominant path and estimate its lifetime.
For two arbitrary dominant paths   and  V , from the two successive WSS time windows, the sets of the two dominant paths are given as To enhance the tracking, the aforementioned threshold  is adopted.The tracking algorithm is described as follows.
Step 1. Calculate the absolute differences between any geometry factors from L  and L +1 and obtain the matrix  with dimension (  ) × ( +1 ).Step 2. If the condition is satisfied, the th MPC at  +1 is considered as evolving from the th MPC at   .For matching, a unique path ID is assigned to them.Examine all other paths in the two successive WSS windows and match the paths evolving over time.Afterwards, if the path at  +1 has no ID, a new path ID is also assigned for the newly observed path.
Step 3. Repeat Steps 1 and 2 for L +1 and L +2 , and inherit the path ID for the path evolving over multiple WSS time windows according to (12).Likewise, if new path has been observed, the new path ID is given.
Step 4. Repeat the steps for the WSS time windows after  +2 , and coherently match in every two consecutive indices of WSS windows.Assign path IDs for all dominant paths.
Step 5. Average the geometry factors collected from the same path ID and assign the mean value as the unique geometry factor to the path ID.
Step 6. Cope with small temporal "gaps," where the information of the dominant path is missed in one or a few (smaller than  gap ) consecutive WSS windows.The set of the missed dominant paths can be covered through the geometry factor in adjacent WSS time windows, that is, involving the geometry factor in ( 5) and ( 6) to calculate the delay and Doppler shifts and then deriving the power by the power model presented in the next section.
The tracking algorithm starts from  = 1 and the dominant paths at initial time are assigned to unique path IDs.In our analysis  gap = 5 is selected.By grouping dominant paths with the path ID, we obtain the lifetime and the temporal evolution for each dominant path.Figure 9 shows the outcome of the tracking algorithm in Doppler domain for the mmWave HST channel.The tracking trajectories calculated by the geometry factors (within their lifetimes) are in agreement with the time-varying DSD from RT simulation.The good match validates the geometry-based detecting and tracking algorithm.

Model Parameters
The channel parameters are first estimated within each WSS time window Δ and then statistically modeled.

Dominant Multipath Lifetime.
If the dominant path is first observed within the  1 th WSS window and is last observed within the  2 th WSS window ( 1 <  2 ), its lifetime is expressed as where the index difference Δ is used to indirectly model lifetime .Note that the LoS path is not included in Δ modeling due to its existing throughout.Figure 10 shows the empirical CDF of Δ and it is well fitted by the truncated Gaussian distribution N(71.5, 48.5), where 71.5 and 48.5 represent the mean value and stand deviation, respectively.In Figure 10Δ𝑠 ranges from 6 to 180 and thus the truncated Gaussian distribution is bounded within [6,180].The fit is implemented by a nonlinear LS regression.To evaluate the fit, the Kolmogorov-Smirnov (KS) test at 95% confidence level is used and it fails to reject the null hypothesis, which validates the fit.In addition, the two blue dash lines shown in Figure 10 are the confidence bounds which are constructed by the given confidence level.Obviously, the empirical CDF of Δ is within the area limited by the two boundaries.the number of old dominant paths.The old dominant paths have been observed before time   2 and their disappearance is determined by their lifetime.Hence, only the number of new dominant paths should be modeled within each WSS time window.

Newly
For accurate modeling, the number of newly emerging dominant paths within the  2 th time window is estimated based on dominant multipath number within the ( 2 − 1)th time window.In detail, according to the conditional probabilities ((  2 → 2 ) | (  2 −1 )) the number of newly emerging dominant paths (  2 → 2 ) can be given.In Table 2 the conditional probabilities are provided.It is evident form the table that one dominant path can be newly observed within each time window at most.Moreover, the probability of generating new path decreases when the dominant multipath number at previous time window is increased.

Newly Dominant Multipath Geometry
Factor.The geometry factor is fixed during dominant multipath lifetime, and, through it, the temporal evolution of delay and Doppler shifts can be calculated by ( 5) and (6).Hence, when a dominant path is newly observed in a WSS time window, the geometry factor should be generated for it.It is found that the PDF of the geometry factor for the newly emerging dominant path can be modeled by Extreme Value distribution: where (⋅) denotes the PDF function of Extreme Value distribution and location parameter  and scale parameter  are model parameters that can be estimated by nonlinear LS regression method.Figure 11 shows the empirical CDF of geometry factor for newly dominant path, that is, Δ, and the LS fit for it, where the good match indicates the reasonable fit.Parameters  and  are estimated as 10.5 and 3.7, respectively.Moreover, the upper and lower confidence bounds (dash lines) obtained by 95% confidence level are presented for comparison.Furthermore, the KS test at the confidence level to evaluate the fit is adopted and the null hypothesis is accepted, which validates the fitted Extreme Value distribution for Δ.

Dominant Multipath Complex Attenuation.
Recalling the parameters in (2), the number of dominant MPCs, delay, and Doppler shift have been modeled above while the complex attenuation of dominant MPCs is investigated here.Note that the parameters, that is, delay, Doppler shift, phase, and complex attenuation, are modeled using the small time-scale, that is,   .Figure 12 shows the 3D scatterplot of delay, Doppler shift, and power of dominant MPCs, which are obtained from RT simulation.Obviously, the dominant multipath path power is dependent on delay as well as Doppler shift.Generally, for the attenuation of MPCs the travel distance and angle information are the major factors.The travel distance can be obtained by the geometrical relationship, which is characterized by the delay obtained from (5).The angle information is included in Doppler shift that is derived from (6).Specifically, when the carrier frequency is given, the Doppler shift is mainly determined by the angle information.Hence, in this paper the complex attenuation is molded with combination of delay and Doppler shift.
To have a good fit to the RT simulated results, the amplitude |  (  )| is modeled in decibel scale and can be divided into two parts, that is, delay-dependent part and Doppler-dependent part.Particularly, amplitude of LoS path only contains delay-dependent part The amplitude model for the other MPCs is expressed as (16), where Δ bp denotes breakpoint for the geometry factor and the parameters  1 ,  2 ,  3 , and  4 for modeling Dopplerdependent part are fitted using nonlinear LS regression.Note that the Doppler shift in ( 16) is in decibel scale.After removing the delay-dependent part from the raw data, the power of Doppler-dependent part obviously presents piecewise behavior, shown in Figure 13, which is the reason for applying piecewise fitting.We postulate that this behavior is due to the different orders of reflection as the higher order with higher path loss.Through the geometrical analysis, the two RT simulated curves comprised of black points can be separated well from each other by the geometry factor.It is reasonable that the higher order reflected paths shall be separated with larger geometry factors, as illustrated in Figure 7. Hence, when  < Δ bp indicates the FRPs, the other in ( 16) is for SRPs.As shown in Figure 13, the piecewise model using exponential functions offers a reasonable fit.Moreover, we note that the fitting curves for Doppler-dependent power have the same shape but different "positions."It can be considered as the exponential function with same shape parameters, that is,  1 ,  2 , moving to different "positions." The "positions" are determined by  3 (along -axis) and  4 (along -axis).In ( 16), the -value of "positions" for SRPs is the double of the FRPs' .The fitting result gives  1 = 1.2,  2 = 3.0,  3 = 45.2, and  4 = −16.0 and the breakpoint Δ bp is set as 65 m in HST scenario.
The phase for dominant paths is expressed as To sum up, the complex attenuations of dominant paths can be calculated as   (  ) = |  (  )| ⋅ exp(  (  )).

Model Implementation and Validation
5.1.Implementation.The implementation steps of the proposed model are summarized as follows.Note that the LoS path throughout exists and thus its generation is deterministic; that is, its lifetime is equal to the duration of MRS moving; geometry factor is equal to Δ LoS (see Figure 7).
(1) Initialization: set the initial dominant multipath number (range from 1 to 5) and the total number of modeled WSS time windows.The lifetime and geometry factor of each dominant path (except LoS path) are initialized by the results generated in Sections 4.1 and 4.3, respectively.Then, by the geometry factor, using ( 5) and ( 6), the dominant multipath delay and Doppler shift are calculated.The complex attenuation of each initial dominant path is derived by its delay and Doppler shift by ( 16) and ( 17).
(2) Modeling for "old" MPCs: in the second WSS window the lifetime for the old dominant multipath is first checked to determine whether to be deleted.For the remaining dominant multipath, involved with the inherited geometry factors from the previous WSS window, their instantaneous delays, Doppler shifts, and complex attenuations are derived the same way as the above step.
(3) Modeling for "new" MPCs: in the second WSS window, a new dominant path generation by the procedure in Section 4.2 is checked.If new one is generated, its lifetime and geometry factor are initialized following the initialization process and then instantaneous delay, Doppler shift, and complex attenuation are derived.
(4) Temporal evolution: by repeating steps (2) and ( 3), the time-varying delays, Doppler shifts, and complex attenuations in spreading function (2) are derived within the consecutive WSS time windows.

Validation.
The validity of the proposed model is examined by comparing the extensive model simulation with the measurement-validated RT simulation.In order to resemble the RT simulation environment, the same system parameters are set, that is, frequency, locations of  and , MRS speed, and moving range.We modeled the spreading functions for 1160 consecutive WSS time windows.Initial dominant multipath number is according to RT simulated results set to 2. Here, within each WSS window, we drop the dominant MPCs with power 25 dB lower than the maximum power [55].
Similarly the RMS Doppler spread  ] (  ) can be calculated by (18) with the variable   (  ) replaced by ]  (  ). Figure 14 depicts the empirical CDFs of   and  ] , which are from RT simulated and modeled channels in Figures 14(a) and 14(b), respectively.It is shown that the model is fairly close to the RT simulation.Note that when   < 5 ns some differences are found between the RT simulated and modeled delay spread, which could emerge because of the model error of the geometry factor that is difficult to accurately model real scenario.Anyway, the overall performance of the mmWave HST channel model is satisfactory.

Conclusion
In this paper, a geometry-based multipath model for mmWave HST channels has been proposed.

Figure 1 :
Figure 1: Top view of mmWave HST scenario deploying the MRS.

Figure 2 :
Figure 2: Front view of typical structure in the measured mmWave HST scenario.

2 Figure 3 :
Figure 3: Directional path loss models estimated by measured and RT simulated channels, respectively, and compared to the free space propagation model.

Figure 4 :
Figure 4: Comparison of CDFs between measured and simulated amplitude of small-scale fading.The Rician distribution is given for modeling.

Figure 5 :
Figure 5: The temporal evolution of power spectral density in Doppler and delay domains for the mmWave HST channel, respectively.
) and 5(b), respectively.Obviously, the LoS path always exists.The other discontinuous Doppler trajectories in Figure 5(a) present at certain consecutive WSS time windows and then vanish.Their corresponding delay trajectories are shown in Figure 5(b).

Figure 6 :
Figure 6: Energy ratio statistical distribution for different order reflected paths (order up to 5th) during the MRS movement.

Figure 8 :
Figure 8: Temporal evolution of refined DSD and the Doppler trajectories for marginal geometry factors obtained (solid line).Detected geometry factor's availability is illustrated by the example trajectory (dashed line).

Figure 9 :
Figure 9: Dominant path tracking in Doppler domain.

Figure 12 :
Figure 12: 3D scatterplot of power in dB versus delay and Doppler shift using the data extracted from RT simulated results.

Figure 13 :
Figure 13: Power of Doppler-dependent part for dominant reflected path and its piecewise model obtained by nonlinear LS regression method.

Figure 14 :
Figure 14: Comparisons of RMS delay spread   and Doppler spread  ] between RT simulated and modeled channels.

Table 2 :
Conditional probabilities for emerging dominant multipath number.Empirical CDF of Δ and its fitted curve with analytical expression.For the fit, two bounds at 95% confidence level are constructed (dash lines).
Dominant Multipath Number.Recall (3); at time   2 the totally dominant multipath number (  2 ) is determined by the number of new dominant paths (  2 → 2 ) and Empirical CDF of Δ for newly dominant path and its fitted curve derived by Extreme Value distribution.