Time- and Frequency-Varying K-Factor of Non-Stationary Vehicular Channels for Safety Relevant Scenarios

Time- and Frequency-Varying -Factor of Non-Stationary Vehicular Channels for Safety Relevant Scenarios

Laura Bernadó, Member, IEEE, Thomas Zemen, Senior Member, IEEE, Fredrik Tufvesson, Senior Member, IEEE, Andreas F. Molisch, Fellow, IEEE, Christoph F. Mecklenbräuker, Senior Member, IEEE Manuscript received yyyyyyy.2014. L. Bernadó and T. Zemen (corresponding author) are with Forschungszentrum Telekommunikation Wien (FTW), Vienna, Austria (e-mail: bernado@ftw.at, thomas.zemen@ftw.at). F. Tufvesson is with the Department of Electrical and Information Technology, Lund University, Lund, Sweden (e-mail: Fredrik.Tufvesson@eit.lth.se). C. F. Mecklenbräuker is with the Institut of Telecommunications, Technische Universität Wien, Vienna, Austria (e-mail: cfm@nt.tuwien.ac.at). A. F. Molisch is with the Department of Electrical Engineering, University of Southern California, Los Angeles, CA, USA (e-mail: molisch@usc.edu). This research was supported by the project NOWIRE funded by the Vienna Science and Technology Fund (WWTF), the project NFN SISE (S10607) funded by the Austrian Science Fund (FWF) as well as the strategic FTW project I-0. FTW is supported by the Austrian Government and the City of Vienna within the competence center program COMET. The research work described in this paper was carried out in cooperation with the COST Action IC1004 on Cooperative Radio Communications for Green Smart Environments. Part of this research has been presented in the IEEE 21st International Symposium on Personal Indoor and Mobile Radio Communications (PIMRC), 2010.

Vehicular communication channels are characterized by a non-stationary time- and frequency-selective fading process due to fast changes in the environment. We characterize the distribution of the envelope of the first delay bin in vehicle-to-vehicle channels by means of its Rician -factor. We analyze the time-frequency variability of this channel parameter using vehicular channel measurements at 5.6 GHz with a bandwidth of 240 MHz for safety-relevant scenarios in intelligent transportation systems (ITS). This data enables a frequency-variability analysis from an IEEE 802.11p system point of view, which uses 10 MHz channels. We show that the small-scale fading of the envelope of the first delay bin is Rician distributed with a varying -factor. The later delay bins are Rayleigh distributed. We demonstrate that the -factor cannot be assumed to be constant in time and frequency. The causes of these variations are the frequency-varying antenna radiation patterns as well as the time-varying number of active scatterers, and the effects of vegetation. We also present a simple but accurate bi-modal Gaussian mixture model, that allows to capture the -factor variability in time for safety-relevant ITS scenarios.

I Introduction

Intelligent transportation systems (ITS) have gained much interest in the last years. They have the potential to strongly reduce the rate of accidents and environmental pollution with the help of reliable wireless vehicle-to-vehicle (V2V) and vehicle-to-infrastructure (V2I) communications.

New receiver algorithms for vehicular communications are first evaluated using numeric link level simulations which rely on accurate channel models. Key factors influencing the wave propagation in vehicular channels are the low position of the antenna on the rooftop of the vehicle, the larger number of metallic objects close to the communication link, and the high mobility of the transmitter (TX) and the receiver (RX). The combination of all these aspects, as in vehicular communications, give rise to a non-stationary fading process. Therefore, it is necessary to study and characterize the small-scale fading statistics in vehicular channels, see e.g. [1, 2, 3, 4, 5, 6, 7, 8, 9], since non-stationary vehicular channel properties have a large impact on the performance of a communication system [10, 11].

In this paper we particularly focus on the characteristic of the Rician -factor of non-stationary vehicular channels. The -factor


is defined as the ratio of the energy of the specular part and the diffuse part of the received signal [12, 13]. Often, the specular part consists only of the line of sight (LOS) component, but specular components can also stem from flat good reflecting surfaces (e.g. traffic signs) where the impinging wave is reflected into a single direction.

The use of the Rician -factor is widespread for describing the small-scale fading characteristics. Therefore the -factor is frequently used as input parameter in numerical simulations to test receiver algorithms, coding schemes, or cooperative strategies [14, 15]. Often the -factor is assumed to be deterministic and is used as constant value for these simulations. However, due to the non-stationary nature of vehicular communication channels this assumption does not hold true, as we will demonstrate in this paper.

In the literature there are only few contributions that characterize the -factor in vehicular environments at GHz. In [7] and [16] the authors show by means of channel measurements that the envelope of the first delay bin is Rician distributed, but they do not give concrete values. -factors between dB and dB are reported in [17] at GHz for a measurement bandwidth of MHz, with an average value of dB. Results in the GHz band are also presented in [18], where the -factor is obtained for a MHz bandwidth system and a value of about dB is reported for different measurement conditions.

All these contributions do not analyze the variability of the small scale fading statistics and only few investigations are available that consider the actual randomness of the -factor. In [19] an analytic model is discussed that takes the randomness of the -factor into account, but no validation with channel measurements is presented. In [20, 21] the non-constant -factor is analyzed based on measurement data centered at GHz considering mobility in a cellular setting. In [20] an analytical model is presented for a dense urban environment using a bandwidth of MHz; and in [21] the authors examine the time and frequency variability of the -factor, and they show a fairly constant -factor in the frequency domain within the measurement bandwidth of MHz.

A first preliminary analysis of the time- and frequency-varying -factor for vehicular communications only for a single rural scenario is presented in [22].

Contributions of the Paper:

  • We focus on the characterization of the small-scale fading statistics for non-stationary V2V channels in a large set of different safety-relevant scenarios for ITS using channel measurement data obtained in the DRIVEWAY’09 measurement campaign [23]. To the best of our knowledge no other publication is available dealing with these important scenarios for ITS.

  • We investigate the variation of the -factor of non-stationary V2V channels in time and frequency. The analysis is done from an IEEE 802.11p system point of view, i.e. the parameters for the analysis are chosen such that they are standard compliant and therefore reliable results for ITS systems can be derived.

  • A simple but accurate bi-modal Gaussian mixture distribution is proposed to model the empiric -factor distribution for each safety-relevant scenario. This approach is related to the one used for the root mean square (RMS)-delay and RMS-Doppler spread modeling in the companion paper [9] for the same scenario set.

Organization of the Paper:

The channel measurements used for carrying out the investigation are described in Sec. II. The methodology for estimating the -factor is presented in Sec. III. The results on the time-frequency variability of the -factor are discussed in Sec. IV for a specifically chosen illustrative exemplary measurement. A statistical modeling of the -factor is proposed in Sec. V, and the results are discussed in detail. We conclude the paper in Sec. VI.

Ii Measurement Data Description

For the analysis of the Rician -factor variability of non-stationary vehicular channels in time and in frequency we use vehicular radio channel measurements collected in the DRIVEWAY’09 measurement campaign [23]. The channel measurements were conducted at a carrier frequency of GHz with a measurement bandwidth of MHz in frequency bins, using the RUSK-Lund channel sounder, based on the switched array sounding principle [24]. The measurements were performed for the following six safety-relevant scenarios in ITS:

  1. road-crossing, in urban and suburban environments;

  2. general LOS obstruction on the highway;

  3. merging lanes in a rural environment;

  4. traffic congestion in two different situations, slow traffic and approaching a traffic jam;

  5. in-tunnel; and

  6. on-bridge.

More details regarding the measurement scenarios can be found in [9, 25].

The vehicles containing the transmitter (TX) and receiver (RX) parts of the channel sounder were equipped with a 4 element linear antenna array. Each antenna element has a main lobe that covers one of the main propagation directions. The transmit and receive antenna indices correspond to the main radiation directions [26], see Fig. 3. With this antenna configuration we obtain 16 spatial links between TX and RX. The mapping between the link index and the TX and RX antenna indices is given by


see also [27, Fig. 5.5 and Tab. 5.2].

For the -factor analysis we use the sampled time-varying frequency response, defined as


where , , and represent discrete time, discrete frequency and the link number while and denote the time and frequency resolution, respectively.

The space variability of the -factor was already demonstrated in [22], therefore we center the analysis performed in the current work on the time-frequency-varying -factor, thus concentrating on a single link , corresponding to the front-front antenna configuration, between the antenna pair and . We therefore drop the link index in all equations for the rest of the manuscript.

The measurement parameters used for the -factor investigation are summarized in Tab. I.

parameter value
channel MIMO

carrier frequency

measurement bandwidth

transmit power

sampling sequence length

frequency bins

time resolution
frequency resolution kHz

recording time
TABLE I: Channel measurement parameters

Iii -factor Estimation

In order to conduct a time-frequency-varying -factor analysis from an IEEE 802.11p system point of view with the available measurement data, we will proceed as follows:

  • Time-variability: For each scenario 3 to 15 measurement runs were recorded, each with a duration of or s. Hence, this allows to characterize the time variation of the -factor.

  • Frequency-variability: The IEEE 802.11p standard for vehicular communications utilizes a communication bandwidth of MHz. Having MHz measured bandwidth available, we can conduct a frequency-dependent investigation of the small scale fading over sub-bands of MHz each. The sampling rate in the delay domain is given by the  MHz bandwidth of each sub-band.

The investigation of the -factor is conducted in the delay domain. For that purpose, we first apply an inverse discrete Fourier transform to the recorded channel transfer function in order to obtain the time-frequency dependent channel impulse response (CIR) for all Q frequency sub-bands


where denotes the time index, the delay index, the frequency sub-band index, and is a Hanning window of length . The values used in the measurements in this publication are or for a s or s measurement run respectively, , and .

Iii-a Data Preprocessing

We characterize the distribution of the envelope of the first strong delay bin. For that, we need to first conduct some data preprocessing, which consists of two steps:

  1. Search for the delay bins corresponding to the first strong multipath component (MPC) in the CIR, and shift them to the origin (ns). This is done on a per link, per frequency sub-band, and per time instance basis. We first perform a maximum search in data chunks of time samples length, where is the number of samples that will be used for removing the large scale fading effects. We then shift the CIR towards the origin by a number of delay bins corresponding to the mode where the maximum is placed within the chunk. This approach does not eliminate the small-scale fading effects.

  2. Remove the large-scale fading from the measured data by applying a moving average filter of length samples in the time domain. We assume the large-scale fading to be stationary during samples [28, 20]. The large-scale fading changes slower than the small-scale fading effects, as a result, the averaging window has to be considered at least as large as the number of samples used later for the -factor estimation. For this reason, we set , where denotes the length of the observation window used for -factor estimation.

    The CIR without large-scale fading is calculated as




    is the average power, represents absolute time, and denotes the local time interval used to calculate the average power.

    After large scale fading removal, the variance of a delay component is an indicator of its -factor. When the received signal is composed of one specular paths its power does not vary. Whereas when the received signal is composed of additional diffuse paths, the superposition of them results in constructive or destructive interference, thus, the envelope varies significantly from one snapshot to the other. The more diffuse components the signal contains, the more variability, and also the smaller the -factor.

    As an example, we show in Fig. 1 the envelope of versus time , for delay , and sub-band . In this figure we depict a setting where the -factor is expected to be low up to and afterwards the -factor is expected to be high, corresponding to periods where the variance of the envelope is large or small, respectively. This will be corroborated in the following sections.

Fig. 1: Magnitude square of the envelope of the CIR of the first delay bin without large-scale fading.

For conducting a first analysis of the -factor, we use the general LOS obstruction scenario on the highway as an illustrative measurement example. Figure 2 shows the time-varying power delay profile (PDP) for this measurement. It will be useful for the reader to refer to this picture as a support for the -factor results discussed in the following sections. The TX and RX are driving in the same direction on the highway at around  km/h (m/s) each, with big trucks driving in both directions beside them, and sometimes obstructing the LOS between the TX and the RX cars.

Fig. 2: Time-varying PDP for a general LOS obstruction on highway scenario: convoy measurement with temporarily obstructed LOS (oLOS), constant velocities of  km/h (m/s).

In this section we focus on the time-frequency variability of the -factor using a single exemplary measurement from link between the antenna pair and , as indicated previously. Later in Section V we perform a statistical analysis of all performed measurement runs in all scenarios.

Details on the radiation pattern of antenna element can be found in Fig. 4, see also [22]. Figure 3 shows an schematic view of the orientation of the radiation patterns for the TX and the RX.

Fig. 3: Scenario layout - general LOS obstruction on highway. In blue the TX car, in red the RX car. The green truck is the element intermittently obstructing the LOS during the measurement.
(a) TX antenna .
(b) RX antenna .
Fig. 4: Radiation pattern of antenna element , points into driving direction [29].

Iii-B Sample Size Selection

It is important to properly select the number of samples used for estimating the -factor. The sample size has to be chosen such that the stationarity length of the process is not exceeded. On the other hand, a large number of samples is needed in order to obtain meaningful statistical results. Furthermore, the condition must be fulfilled.

We choose the length such that it corresponds to travelled distance. Since the snapshot repetition time is fixed to , the number of samples varies depending on the measured scenario, i.e. depending on the speed of the TX and RX. Table II lists the average speed and the selected for the different scenarios 111Due to unreliable GPS data, the window length is not adapted to the instantaneous velocity. Instead, we use average speed values to compute .. The number of samples used for removing the large-scale fading effects is then selected as , being added samples at the beginning and end of the window .

Scenario Average speed Number of samples
road crossing  m/s ( km/h)
general LOS obstruction  m/s ( km/h)
merging lanes  m/s ( km/h)
slow traffic  m/s ( km/h)
approaching traffic jam  m/s ( km/h)
in-tunnel  m/s ( km/h)
on-bridge  m/s ( km/h)
TABLE II: : Number of samples used for estimating the -factor

Iii-C Envelope Distribution per Tap

Even though it is well accepted that the Rician distribution is the best for describing the small-scale fading effects, before further proceeding with the -factor analysis, we evaluate the fitting of the envelope of the first delay bins to different distributions. In order to decide which distribution fits best, we apply the Kolmogorov-Smirnov (KS) test as a goodness-of-fit (GoF) indicator [30], calculated as


where denotes the supremum, and is the empirical cumulative distribution function (CDF), and is the analytical CDF, of the random variable . We consider to achieve a good fit if the outcome of a KS-test is below , which indicates that the maximum distance between empirical and fitted CDF is lower than . Figure 5 shows the GoF result for the first taps for different distributions for the exemplary measurement run for link , and , corresponding to MHz. Different markers correspond to different distributions, the GoF plotted in the figure is the result of averaging the GoF values for time instances (). We corroborate from this figure, that the Rician distribution fits best the data, thus we carry on with the -factor analysis.

Fig. 5: Goodness of fit for the first taps using different distributions to fit their envelope.

Figure 6 (a), (b), and (c) show the CDF of the first five delay bins from . The measured data is represented by a solid line and the fitted CDF uses a dashed line. Results are shown for frequency sub-band , corresponding to a frequency of MHz, at time instances  s,  s, and  s respectively. At  s the LOS between TX and RX is obstructed, therefore the distribution of the first tap (denoted ’tap 0’) is closer to being Rayleigh (note the shift towards the left hand side of the plot). At s, the truck blocking the LOS has moved to another lane allowing a LOS communication. As a result, the -factor of the first delay bin increases and its CDF is shifted to the right in Fig. 6 (b). The same effect can be noticed at s (Fig. 6 (c)).

(a) s
(b) s
(c) s
Fig. 6: Empirical and fitted CDF for delay bins and sub-band (MHz), using a sample size of samples at different time instances for a general LOS obstruction on highway scenario.
(a) s
Fig. 7: Weibull plot for delay bins and sub-band (MHz), using a sample size of samples for for a general LOS obstruction on highway scenario.

Furthermore, we observe that delay bins (denoted ’tap 1…4’) follow a Rayleigh distribution throughout all time instances. On the other hand, the first delay bin is clearly Rician distributed with a varying -factor for the three different time instances, dB, dB, and dB.

For we also show an examplary Weibull plot [31] of the first 5 tap distributions in Fig. 7. The Weibull distribution


where is the shape parameter and is the scale parameter. The CDF is given as


for , and for . For the Weibull distribution is identical with the Rayleigh distribution and for it is analogous to a Rician distribution [4], [32]. For the Weibull plot the empirical cumulative distribution function is shown using axes vs. resulting in a linearization of the Weibull CDF


From Fig. 7 we can see that the taps have a Rayleigh distribution with shape parameter while tap is clearly Rician distributed with . We choose the time instance  s for exemplification because of the clear values of the shape parameter , however, this is also valid for the other time instances, being their shape parameters listed in Tab. III.

shape parameter  s  s  s
TABLE III: Shape parameters for the Weibull distribution

The analysis in this section for the the delay bins clearly demonstrates, that the delay bin is strongly Rician distributed with a factor larger than zero while the delay bin are mostly Rayleigh distributed. Hence, we will focus our analysis on delay bin in this paper.

Iii-D -Factor Estimation

The -factor estimate


is calculated for the first delay bins , where denotes the -factor estimation from a finite set of observation at time indices . For estimating the -factor we use a technique introduced in [33] which is based on the method of moments (MoM) [34].

Iv Empirical Results

The goal of this investigation is to characterize the variation of the -factor in time and frequency. The analysis is only carried out for the first delay bin .

The time- and frequency-varying -factor estimate is plotted in Fig. 8, next to its corresponding time- and frequency varying power of the first delay bin. Here, the time-frequency variation of the -factor can be appreciated already.

Noteworthy is that there is not necessarily a correspondence between received power and -factor (see the scale in the figures). This is due to the fact that the received power is calculated as the sum of the power of all specular and diffuse components, , whereas the -factor is an indicator of the ratio between them, see (1).

(a) Time-frequency dependent -factor.
(b) Time-frequency dependent power.
Fig. 8: General LOS obstruction on highway: Time-frequency dependent -factor and power for tap of link .

On what follows, we discuss the time and frequency variability of the -factor of the first delay bin independently. We remind here that the chosen illustrative measurement corresponds to a general LOS obstruction on highway scenario. In the studied measurement run, the LOS between the TX and the RX cars is obstructed for the first s, afterwards, the truck blocking the LOS moves away.

Due to the large variation of the factor all following figure displaying the -factor will be given in logarithmic scale. We would like to point out, that for regions where the -factor is close to , the logarithmic scale might overemphasize its variation since the range in linear scale is mapped to the range in logarithmic scale. Please keep this in mind when interpreting the figures.

Iv-a Time-Varying -Factor

Figure 9 shows the time evolution of the -factor for three different frequency sub-bands, together with their corresponding power. The three curves are obtained by making a cut in Fig. 8 at the three different frequency sub-bands. The -factors obtained at the three frequency bands are comparable although they do not have the same numerical value. Furthermore, they change in time following a similar tendency.

Fig. 9: General LOS obstruction on highway: Time-varying -factor and power of the first delay tap without large-scale fading for three different frequency sub-bands for a sample size of .

Large -factors up to dB are observed between and  s. However, the estimated -factor is around dB smaller for the period from to  s. In this case, there are more significant diffuse components in the received signal, and thus it is better described by a lower -factor. We call obstructed LOS (oLOS) the situation when the LOS between TX and RX is obstructed, but the waves still propagate by diffraction over the roof-top of the obstructing vehicle, as it occurs between seconds and .

We observe estimated -factor values for the oLOS to be between and  dB, being those lower than the -factor in pure LOS communication, with values around  dB. Lower -factors are obtained in non-LOS (nLOS) situations (shown in results in Sec. V). Thus, we can establish the following relationship: .

Iv-B Frequency-Varying -Factor

For the frequency-variation analysis we select four different time instances and plot their estimated -factor as a function of frequency sub-bands in Fig. 10, together with their corresponding power. The three curves here are obtained by making a cut in Fig. 8 at three different time instances . Looking at the figure, it can be concluded that the assumption of a constant -factor in frequency (within the MHz measured frequency band) is not realistic.

Fig. 10: General LOS obstruction on highway: Frequency-varying -factor and power of the first delay tap without large-scale fading for three different time instances for a sample size of .

Variations of up to dB are observed throughout the frequency sub-bands, mainly due to the frequency dependent antenna radiation pattern. The variation of the antenna gains for antenna elements and shown in Fig. 4 are significant, when comparing the gains at the lower, central, and upper frequency band. This effect is mostly observed for element , where the gain experiences a variation of  dB at the bandwidth edges with respect to the carrier frequency of MHz. The directionality of the pattern is affected by the frequency variations, but notches are significant, such as at , and , giving rise to drops at some frequencies, as observed in Fig. 8 (a). The -factor also varies in frequency due to the effect of foliage and vegetation in the surrounding environment.

The results for are shown in Fig. 10 to provide a link to Fig. 9 and Fig. 8.a where a very low -factor near zero, or in the logarithmic domain of , is obtained for this time instant.

V Statistical Modeling

V-a Distribution Fitting

We are interested in characterizing the distribution of the -factor for the sub-band bandwidth of . In order to see whether we can use the whole data ensemble of the link for the statistical characterization, we first have a look at the CDF of the envelope of the first tap for the frequency sub-bands individually, as done in Sec. III-C for a single sub-band. The result of applying the KS-test shows values of for all sub-bands, which we consider low enough and therefore conclude that the first tap is described by a Rician distribution in all sub-bands. Therefore we will use the estimated -factors from all available sub-bands jointly for fitting a distribution.

We show in Fig. 11 the normalized histogram (histogram divided by the total number of samples) for three different scenarios. Two Gaussian-shaped distributions can be envisioned from these histograms, which motivates us to select a bi-modal Gaussian distribution for fitting purposes, similar as in the companion paper for the root-mean-square (RMS) Delay and Doppler Spread [9].

Fig. 11: Normalized histogram (histogram/number of samples) for three different scenarios.

One could argue that the pdf of the K-factor is different depending on the selected sample ensemble, which is time-dependent, for each measurement run. For instance, for the exemplary general LOS obstruction scenario, if we separate the CIRs of our measurement runs into two ensembles, one set containing the CIRs when the LOS is not obstructed and the other one when the LOS is obstructed, the resulting pdf of the K-factor for the two sets of CIRs will be different. However, we consider adequate to characterize both sets jointly with a single bi-modal distribution, in order to describe the transition probability between them as well.

The probability density function (pdf) of a bi-modal Gaussian mixture distribution reads


where and are the mean and the standard deviation of each Gaussian component with index , and and are the weighting factors of the first and second component respectively. The first Gaussian function describes the nLOS situation (small -factor), the other one describes the LOS (large -factors) situation. The weighting factors of the two Gaussians and indicate the probability of having nLOS or LOS for a particular scenario.

In Fig. 12, we plot the empirical CDF of the estimated -factor as solid line, and the fitted CDF


as dashed line.

We apply the KS test as a GoF indicator


to assess the accuracy of the bi-modal Gaussian distribution [9]. For all safety relevant measurement scenarios we obtain , thus showing a very good fit, as listed in Tab. IV. Hence, the bi-modal Gaussian distribution (12), based on only five parameters, is a simple, but at the same time accurate model for the -factor distribution in safety-relevant vehicular scenarios.

Fig. 12: General LOS obstruction on highway: Joint CDF - time-frequency dependent -factor.

V-B Discussion

V-B1 Weighting factor

Looking at the results, in dB, in Tab. IV, we see that in some cases the weighting factor of the one of the component is very small, namely for the general LOS obstruction, merging lanes, and approaching traffic jam, where one of the weighting factors is below . Comparable weighting factors for the both components of the bi-modal Gaussian distribution are obtained for road crossing - urban scenarios, for both single lane and multiple lanes. Precisely in these two scenarios the cars experience nLOS conditions more often, since there are buildings in the four quadrants of the crossing.

A particular case exhibiting similar weighting factors is the on-bridge scenario. Nevertheless, here we do not have two clear Gaussian distributions, since the two mean values, and are less than  dB apart from each other. This indicates that the actual shape of the histogram for this particular scenario is close to a single Gaussian.

The other scenarios, road crossing in suburban areas, slow traffic, and in-tunnel, have different weighting factors, but with very dis-similar mean values for each component.

V-B2 Mean value

In general, the two mean components, and are well separated being . The lower component corresponds to the nLOS period of the measurement, as discussed previously, with values around  dB in most of the cases (we have already commented on the special case of on-bridge scenario).

If we now concentrate on , the highest value is obtained for the merging lanes (14.2 dB) and on-bridge (14.6 dB) scenarios, where there is a clear LOS (the two cars drive one after the other without obstacles between them).

On the other extreme, the lowest are obtained in the two scenarios in urban environment, road crossing - single lane (-0.6 dB) and road crossing - multiple lane (0.11 dB), with comparable values. In these measurements (urban) the number of scatterers close to the TX-RX link is large, leading to low -factors [35].

In general, scatterers are closer in road crossing scenarios. Nevertheless, in suburban environments, their number is lower, and for the road crossing - suburban scenarios are a bit higher than in urban environments. The mean value of the -factor, , for situations with traffic (3.7 dB) is sightly lower than for situations without traffic (4.5 dB), as the more traffic (i.e. more vehicles around), the more scattering. The rich scattering leading to low values is also present in in-tunnel (dB), and slow traffic (dB) scenarios. Finally, the values for the two scenarios on the highway are also comparable, being 7.6 dB for the general LOS obstruction scenario, and 8.1 dB for the approaching traffic jam scenario.

To summarize, we can list the values from low to high and also depending on the richness of the scattering environment as: crossing-urban (-0.6, 0.1 dB) crossing-suburban (3.7, 4.5 dB) slow-traffic (4.4 dB) in-tunnel (4.7 dB) highway (general LOS obstruction and approaching traffic jam) (7.6, 8.1 dB) merging lanes (14.2 dB).

V-B3 Standard deviation

Nearly all have comparable values among themselves, as well as all among themselves. The standard deviation of the first component is around 7 dB, except for the merging lane scenario, which is 21.7 dB. We do not give importance to this value, since the the first component in this particular scenario is almost negligible ().

The standard deviation of the second component, , is around 5 dB for all scenarios, except for the general LOS obstruction (7.5 dB) and the two traffic congestion (6.5 dB) scenarios. In these two cases, oLOS occurs with higher frequency, with values close to the LOS case (discussed in Sec. IV), thus being reflected in larger standard deviations in our modeling.

 [dB]  [dB]  [dB]  [dB] runs
road crossing - suburban with traffic
road crossing - suburban without traffic
road crossing - urban single lane
road crossing - urban multiple lane
general LOS obstruction - highway
merging lanes - rural
traffic congestion - slow traffic
traffic congestion - approaching traffic jam
TABLE IV: Modeling the time-varying -factor

Vi Conclusions

In this contribution, we analyzed and characterized the small scale fading from vehicular channel measurements. We found that the first delay bin is Rician distributed with a varying -factor, and the following delay bins are mostly Rayleigh distributed. We also showed that the -factor varies in the time- and frequency-domain. These variations are significant and must be taken into account for modeling the propagation conditions for safety relevant scenarios in intelligent transportation systems (ITS).

The -factor variation depends mainly on the following factors: (i) The number of illuminated scatterers due to the antenna radiation pattern and objects in between TX and RX, (ii) the antenna radiation pattern variation over the measurement bandwidth, (iii) vegetation and foliage, and (iv) the presence of good reflecting objects in the vicinity of the TX-RX link (e.g. traffic signs).

Furthermore, we provided a statistical model for the -factor of the first delay bin in a MHz communication bandwidth for important safety critical vehicular scenarios matched to the IEEE 802.11p standard. For that, we used a bi-modal Gaussian distribution, where each one of the components models the -factor in LOS and nLOS propagation conditions.

The weighting factors can be interpreted as the transition probabilities between LOS and nLOS situations. They are similar for road-crossing scenarios in urban environments. Otherwise, one component is more important than the other (with a higher weighting factor), except for the general LOS obstruction, merging lanes, and traffic congestion scenarios, where the obstructed LOS (oLOS) situation is more common than nLOS.

We also saw that the mean values corresponding to the nLOS Gaussian component are around dB. The mean values corresponding to LOS and oLOS situations are different depending on the scattering richness of the surroundings, being  dB the lowest for road-crossing in urban environment, and dB for merging lanes in rural environment.

Finally, we also analyzed the standard deviation of the two Gaussian components showing that the largest values are obtained for situations with LOS and oLOS mixed, since their mean values are close enough.

We consider that the work presented in this manuscript is the basis for more complex modeling of a variable -factor, where time correlation can be included. Furthermore, the bi-modal Gaussian fitting can open the doors to a Markov-chain model of the -factor for non-stationary channels.


  • [1] G. Matz, “On non-WSSUS wireless fading channels lololo,” IEEE Trans. Wireless Commun., vol. 4, no. 5, pp. 2465–2478, September 2005.
  • [2] T. Willink, “Wide-sense stationarity of mobile MIMO radio channels,” IEEE Trans. Veh. Commun., vol. 57, no. 2, pp. 704–714, March 2008.
  • [3] L. Bernadó, T. Zemen, A. Paier, G. Matz, J. Karedal, N. Czink, F. Tufvesson, M. Hagenauer, A. F. Molisch, and C. F. Mecklenbräuker, “Non-WSSUS Vehicular Channel Characterization at 5.2 GHz - Spectral Divergence and Time-Variant Coherence Parameters,” in Assembly of the International Union of Radio Science (URSI), August 2008, pp. 9–15.
  • [4] O. Renaudin, V.-M. Kolmonen, P. Vainikainen, and C. Oestges, “Non-stationary narrowband MIMO inter-vehicle channel characterization in the 5-GHz band,” IEEE Trans. Veh. Technol., vol. 59, no. 4, pp. 2007–2015, May 2010.
  • [5] C. F. Mecklenbräuker, A. F. Molisch, J. Karedal, F. Tufvesson, A. Paier, L. Bernadó, T. Zemen, O. Klemp, and N. Czink, “Vehicular channel characterization and its implications for wireless system design and performance,” Proc. IEEE, vol. 99, no. 7, pp. 1189–1212, July 2011.
  • [6] A. F. Molisch, F. Tufvesson, J. Karedal, and C. F. Mecklenbrauker, “A survey on vehicle-to-vehicle propagation channels,” IEEE Wireless Commun. Mag., vol. 16, no. 6, pp. 12–22, December 2009.
  • [7] J. Maurer, T. Fugen, and W. Wiesbeck, “Narrow-band measurement and analysis of the inter-vehicle transmission channel at 5.2 GHz,” in IEEE 55th Vehicular Technology Conference (VTC Spring), vol. 3, 2002, pp. 1274–1278.
  • [8] G. Acosta, K. Tokuda, and M. Ingram, “Measured joint Doppler-delay power profiles for vehicle-to-vehicle communications at 2.4 GHz,” in IEEE Global Telecommunications Conference (GLOBECOM), vol. 6, December 2004, pp. 3813–3817.
  • [9] L. Bernadó, T. Zemen, F. Tufvesson, A. Molisch, and C. Mecklenbräuker, “Delay and Doppler spreads of non-stationary vehicular channels for safety relevant scenarios,” IEEE Trans. Veh. Technol., vol. 63, no. 1, pp. 82 – 93, January 2014.
  • [10] T. Zemen, L. Bernado, N. Czink, and A. F. Molisch, “Iterative time-variant channel estimation for 802.11p using generalized discrete prolate spheroidal sequences,” IEEE Trans. Veh. Technol., vol. 61, no. 3, pp. 1222–1233, March 2012.
  • [11] T. Zemen and A. F. Molisch, “Adaptive reduced-rank estimation of non-stationary time-variant channels using subspace selection,” IEEE Trans. Veh. Technol., vol. 61, no. 9, pp. 4042–4056, November 2012.
  • [12] S. Rice, Statistical properties of a sine wave plus random noise, ser. Technical publications: Monograph.   Bell Telephone Laboratories, 1948, vol. 27.
  • [13] G. Durgin, Space-time wireless channels, 1st ed.   Upper Saddle River, NJ, USA: Prentice Hall Press, 2002.
  • [14] P. Castiglione, S. Savazzi, M. Nicoli, and T. Zemen, “Impact of fading statistics on partner selection in indoor-to-outdoor cooperative networks,” in IEEE International Conference on Communications (ICC), May 2010.
  • [15] K.-Y. Lin, R.-T. Juang, H.-P. Lin, W.-J. Shyu, and P. Ting, “Link adaptation of MIMO-OFDM transmission exploiting the Rician channel K-factor,” in IEEE Mobile WiMAX Symposium (MWS), July 2009, pp. 184–188.
  • [16] J. Karedal, “Measurement-based modeling of wireless propagation channels - MIMO and UWB,” Ph.D. dissertation, Lund University, 2009.
  • [17] L. Thiele, M. Peter, and V. Jungnickel, “Statistics of the Ricean K-factor at 5.2 GHz in an urban macro-cell scenario,” in IEEE 17th International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC), September 2006.
  • [18] O. Renaudin, V.-M. Kolmonen, P. Vainikainen, and C. Oestges, “Wideband MIMO car-to-car radio channel measurements at 5.3 GHz,” in 68th IEEE Vehicular Technology Conference (VTC Fall), September 2008.
  • [19] M. Ilic-Delibasic and M. Pejanovic-Djurisic, “Impact of random K factor on Ricean fading wireless system performance,” in 16th IEEE Mediterranean Electrotechnical Conference (MELECON), March 2012, pp. 233–236.
  • [20] G. G. Messier and J. A. Hartwell, “An empirical model for nonstationary Ricean fading,” IEEE Trans. Veh. Technol., vol. 58, no. 1, pp. 14–20, January 2009.
  • [21] L. J. Greenstein, S. Ghassemzadeh, V. Erceg, and D. G. Michelson, “Ricean K-factors in narrow-band fixed wireless channels: Theory, experiments, and statistical models,” IEEE Trans. Veh. Technol., vol. 58, no. 8, pp. 4000–4012, October 2009.
  • [22] L. Bernadó, T. Zemen, J. Karedal, A. Paier, A. Thiel, O. Klemp, N. Czink, F. Tufvesson, A. F. Molisch, and C. F. Mecklenbräuker, “Multi-dimensional K-factor analysis for V2V radio channels in open sub-urban street crossings,” in 21st Annual IEEE International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC), Istanbul, Turkey, September 2010.
  • [23] A. Paier, L. Bernadó and, J. Karedal, O. Klemp, and A. Kwoczek, “Overview of vehicle-to-vehicle radio channel measurements for collision avoidance applications,” in 71st IEEE Vehicular Technology Conference (VTC Spring), May 2010.
  • [24] R. Thomä, D. Hampicke, A. Richter, G. Sommerkorn, A. Schneider, U. Trautwein, and W. Wirnitzer, “Identification of time-variant directional mobile radio channels,” IEEE Trans. Instrum. Meas., vol. 49, no. 2, pp. 357–364, April 2000.
  • [25] L. Bernadó, “Non-stationarity in vehicular wireless channels,” Ph.D. dissertation, Vienna University of Technology, Vienna, Austria, April 2012.
  • [26] A. Thiel, O. Klemp, A. Paiera, L. Bernado, J. Karedal, and A. Kwoczek, “In-situ vehicular antenna integration and design aspects for vehicle-to-vehicle communications,” in 4th European Conference on Antennas and Propagation (EuCAP), April 2010.
  • [27] L. Bernadó, “Non-Stationarity in Vehicular Wireless Channels,” Ph.D. dissertation, Vienna University of Technology, 2012.
  • [28] R. H. Clarke, A statistical theory of mobile-radio reception.   Bell Telephone Laboratories, 1968, vol. 47.
  • [29] L. Bernadó, A. Roma, A. Paier, T. Zemen, N. Czink, J. Karedal, A. Thiel, F. Tufvesson, A. F. Molisch, and C. F. Mecklenbräuker, “In-tunnel vehicular radio channel characterization,” in 73rd IEEE Vehicular Technology Conference (VTC Spring), May 2011.
  • [30] F. J. Massey, “The Kolmogorov-Smirnov test for goodness of fit,” Journal of the American Statistical Association, vol. 46, no. 253, pp. 68–78, 1951. [Online]. Available: http://www.jstor.org/stable/2280095
  • [31] J. I. McCool, Using the Weibull Distribution: Reliability, Modeling and Inference.   Wiley, 2012.
  • [32] D. Matolak and J. Frolik, “Worse-than-Rayleigh fading: Experimental results and theoretical models,” Communications Magazine, IEEE, vol. 49, no. 4, pp. 140–146, April 2011.
  • [33] C. Oestges, N. Czink, B. Bandemer, P. Castiglione, F. Kaltenberger, and A. J. Paulraj, “Experimental characterization and modeling of outdoor-to-indoor and indoor-to-indoor distributed channels,” IEEE Trans. Veh. Technol., vol. 59, no. 5, pp. 2253–2265, June 2010.
  • [34] L. Greenstein, D. Michelson, and V. Erceg, “Moment-method estimation of the Ricean K-factor,” IEEE Communications Letters, vol. 3, no. 6, pp. 175–176, 1999.
  • [35] T. Abbas, J. Karedal, F. Tufvesson, A. Paier, L. Bernadó, and A. F. Molisch, “Directional analysis of vehicle-to-vehicle propagation channels,” in Vehicular Technology Conference (VTC Spring), 2011 IEEE 73rd, May 2011.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description