Next Article in Journal
An Intrusion Detection System Based on Multi-Level Clustering for Hierarchical Wireless Sensor Networks
Next Article in Special Issue
A Passive Testing Approach for Protocols in Wireless Sensor Networks
Previous Article in Journal
Correlating the Integral Sensing Properties of Zeolites with Molecular Processes by Combining Broadband Impedance and DRIFT Spectroscopy—A New Approach for Bridging the Scales
Previous Article in Special Issue
A Novel Scheme for an Energy Efficient Internet of Things Based on Wireless Sensor Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Timing Estimation Method Based-on Skewness Analysis in Vehicular Wireless Networks

Department of Computer and Communication Engineering, China University of Petroleum (East China), Qingdao 266580, China
*
Author to whom correspondence should be addressed.
Sensors 2015, 15(11), 28942-28959; https://doi.org/10.3390/s151128942
Submission received: 11 October 2015 / Revised: 9 November 2015 / Accepted: 11 November 2015 / Published: 13 November 2015
(This article belongs to the Special Issue Identification, Information & Knowledge in the Internet of Things)

Abstract

:
Vehicle positioning technology has drawn more and more attention in vehicular wireless networks to reduce transportation time and traffic accidents. Nowadays, global navigation satellite systems (GNSS) are widely used in land vehicle positioning, but most of them are lack precision and reliability in situations where their signals are blocked. Positioning systems base-on short range wireless communication are another effective way that can be used in vehicle positioning or vehicle ranging. IEEE 802.11p is a new real-time short range wireless communication standard for vehicles, so a new method is proposed to estimate the time delay or ranges between vehicles based on the IEEE 802.11p standard which includes three main steps: cross-correlation between the received signal and the short preamble, summing up the correlated results in groups, and finding the maximum peak using a dynamic threshold based on the skewness analysis. With the range between each vehicle or road-side infrastructure, the position of neighboring vehicles can be estimated correctly. Simulation results were presented in the International Telecommunications Union (ITU) vehicular multipath channel, which show that the proposed method provides better precision than some well-known timing estimation techniques, especially in low signal to noise ratio (SNR) environments.

1. Introduction

Currently, vehicle collision accident rates remain high and cause great damages every year. The World Health Organization (WHO) has indicated that globally road traffic injuries are the eighth leading cause of death [1]. More than a million people die each year on the World’s roads and current trends suggest that by 2030 road traffic deaths will become the fifth leading cause of death unless more action is taken. One effective action is the application of cooperative collision avoidance (CCA), by which vehicles can predict potential hazards and remind the drivers in advance about hazards by exchanging real-time information such as location, direction, speed, acceleration and so on. For a CCA system, the vehicle’s location is one of the most important parameters to be estimated. If drivers could receive a collision warning 0.5 s before an incident, 60% of collisions could be avoided, but how to accurately and reliably estimate the vehicle’s location is still a challenge. Without an accurate position estimation, the CCA systems cannot produce timely alerts about potential dangers ahead or may provide some unwanted warnings when there is no danger ahead.
Currently, global navigation satellite systems (GNSS) such as the US global positioning system (GPS) and the BeiDou satellite positioning system (BDS) of China have been widely used in vehicular wireless networks to reduce transportation time and traffic accidents. If the satellite signals are blocked in a non-line of sight (NLOS) environment, such as in urban areas where tall buildings or overpass bridges or dense foliage may block the view of the sky they will fail to obtain accurate position estimations. GNSS must therefore be integrated with some other positioning technique such as short range wireless position, Inertial Navigation System (INS), digital maps, and so on. The Chinese Area Positioning System (CAPS) is a regional navigation satellite system. In [2] ultra-wide-band (UWB) pseudo-satellite signals were used in the CAPS, and the integrated method could provide better performance than CAPS itself. A three layer framework was presented in [3] to monitor the performance of real-time relative positioning in vehicular wireless networks that support CCA applications. The IEEE 802.11p [4] standard is just employed to translate the GPS real-time kinematic (RTK) information. Unfortunately, these systems also perform poorly in NLOS environments.
Nowadays, wireless mobile communications systems or cellular networks [5] are widely used in many areas, including in vehicular wireless networks, but they cannot be used for positioning, because they have no real-time communication ability. This is because each message will be delivered through different cellular base stations, which will introduce insufferable time delays compared with the rapid velocity of vehicles. In addition, the positioning precision of cellular signals is about several hundred meters so it can't meet the demands of CCA systems. Traditional wireless sensor networks [6] or indoor positioning methods [7,8] are not applicable in vehicle environments, because of the severe Doppler frequency shift caused by the high velocity of vehicle movement. Fortunately, IEEE 802.11p was approved in 2010 to add wireless access in vehicular environments, which made some amendments to the IEEE 802.11 in the physical (PHY) and media access control (MAC) layers to achieve a robust connection and a fast setup, such as the longer symbol duration, guard time and FFT period.
In order to improve the timing precision for vehicular nodes, a new method with a dynamic threshold using the skewness analysis of IEEE 802.11p short preamble is proposed in this paper. Against this background, the OFDM preamble structure is analyzed in Section 2. Section 3 introduces the related work on timing estimation, and the proposed method is put forward in Section 4. Simulation setup and result analysis are discussed in Section 5, and finally Section 6 concludes the paper. The main contributions of this paper include:
(1)
A new cross-correlation, summing up and skewness analysis (denoted CSS) method is proposed which can be used to estimate the time offset in positioning sensors of vehicular wireless networks.
(2)
The CSS method provides better precision than some well-known time estimation techniques, especially in low signal to noise ratio (SNR), and multi-path environments, which is very important for vehicular wireless networks.
(3)
The CSS method can also be used in other wireless communications which use OFDM with repeated preamble symbols, such as IEEE 802.11a.

2. OFDM Preamble Structure

IEEE 802.11p specifies the extensions to the PHY and MAC layers of IEEE 802.11 for wireless local area networks providing wireless communications while in a vehicular environment, e.g., at high speeds of up to 200 km/h, multi-path fading and long range of range of 1000 m. In IEEE 802.11p, the same OFDM modulation technique and preamble structure as in IEEE 802.11a [9] are used to mitigate the effects of fading, but the larger symbol duration is specified, which results in a longer guard interval (GI) and can provide more tolerance to greater delay spreads and inter-symbol interference (ISI). Thus, IEEE 802.11p can deal with the severe Doppler frequency shifts, and rapidly changing multi-path conditions.

2.1. Modulation Technique

With the modulation technique of OFDM, the inverse fast Fourier transform (IFFT) is used to generate the OFDM symbols. In each OFDM symbol, there are multiple orthogonal subcarriers. In IEEE 802.11p, the channel bandwidth of each subcarrier is 10 MHz and the carrier frequency is 5.9 GHz. Therefore, the frequency spacing between each subcarrier is F = 10/64 = 0.15625 MHz, which is half of the channel bandwidth of IEEE 802.11a. The duration of IFFT (TFFT) in IEEE 802.11p is 6.4 μs, which is twice that of IEEE 802.11a. The duration of GI is 1.6 μs.

2.2. Preamble Symbols

The PHY layer format of IEEE802.11p is shown in Figure 1. The first two part are the preamble fields, which are divided into two parts, 10 short preamble symbols and two long preamble symbols. In the frequency-domain, each short preamble symbol consists of 12 nonzero subcarriers located at positions ±(4, 8, 12, 16, 20, and 24) which gives:
S 26 , 26 = 13 / 6 { 0 , 0 , + 1 + j , 0 , 0 , 0 , 1 j , 0 , 0 , 0 , + 1 + j , 0 , 0 , 0 , 1 j , 0 , 0 , 0 , 1 j , 0 , 0 , 0 , + 1 + j , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 j , 0 , 0 , 0 , 1 j , 0 , 0 , 0 , + 1 + j , 0 , 0 , 0 , + 1 + j , 0 , 0 , 0 , + 1 + j , 0 , 0 , 0 , + 1 + j , 0 , 0 ,   }
The OFDM symbol is defined as:
s ( t ) = w ( t ) k = N / 2 k = N / 2 S k e x p ( j 2 π k F t )
where N = 52 and w ( t ) is the window function. Because there are only 12 nonzero subcarriers, the factor 13 / 6 is used to normalize the average power of the OFDM symbol. Because subcarriers s has nonzero values only for indices that are a multiple of 4, the short preamble symbol duration is TFFT/4 = 1.6 μs. Therefore, the duration of the ten short preamble symbols is 10 × 1.6 = 16 μs.
There are 52 nonzero subcarriers and 12 null subcarriers with a duration of 6.4 μs in each long preamble symbol. Between short preamble and long preamble, there are two guard interval symbols (G12) with the duration of 1.6 μs. Therefore, totally, the duration of the preamble is 10 × 1.6 + 2 × 1.6 +2× 6.4 = 32 μs.
The short preamble symbols are typically used for coarse timing estimation and the long preamble symbols are used for fine timing estimation and frequency shift estimation. The former is critical to the success of the fine estimation so it is considered in this paper.
Figure 1. Format of IEEE802.11p PHY layer.
Figure 1. Format of IEEE802.11p PHY layer.
Sensors 15 28942 g001

3. Related Work

Most of the wireless positioning methods are based on ranging or timing delay estimation from the base station to the target positioning sensor (the vehicular node to be positioned), for example time of arrival (TOA) [10], time difference of arrival (TDOA) and received signal strength (RSS) [11,12]. However, there is a very challenging problem for ranging estimation due to the severe multi-path, reflection, and inter-symbol interference environments encountered [13,14]. Therefore, the ranging or timing delay estimation problem has been extensively studied. Ranging or positioning methods in the traditional application sense have been studied more extensively [15,16,17], but in the vehicle environment less so [18]. Because orthogonal frequency division multiplexing (OFDM) is used in the IEEE802.11p standards, we focus on some timing methods for OFDM.

3.1. SC Method

Schmidl and Cox developed the well-known timing estimation method (denoted as SC) [19] based on OFDM, which is recommended in the IEEE802.11a standard [20]. However, the timing metric will reach a plateau which leads to uncertainty in the timing estimation. In the SC method [19], the preamble symbols can be descripted as the pattern of [GI A A] where GI is the guard interval symbols which are followed by two groups of identical symbols. The metric for timing estimation in SC method is given by:
M ( d ) = | P ( d ) | 2 R 2 ( d )
where P(d) is the auto-correlation of the received samples and R(d) is the energy of the received samples, given by:
P ( d ) = m = 0 L 1 r * ( d + m ) · r ( d + m + L )
R ( d ) = m = 0 L 1 | r ( d + m + L ) | 2
where L is length of the received samples, d is the time index of the received samples and * is the complex conjugate.
Figure 2 presents |P(d)|2, R2(d) and M(d) with L = 32 in an ideal channel with a time delay of 191 samples. For the IEEE802.11p protocol, there are 160 short preamble samples, so the first 160 − 2L = 96 samples are considered to be the GI. The timing metric M(d) forms a plateau which has a length equaling to the length of the GI (Ng = 96 samples) minus the length of channel delay spread. In an AWGN channel, there will be no delay spread, so the length of the plateau is 96 samples.
Figure 2. |P(d)|2, R2(d) and M(d) of SC method with L = 32 in ideal environment.
Figure 2. |P(d)|2, R2(d) and M(d) of SC method with L = 32 in ideal environment.
Sensors 15 28942 g002

3.2. RMB Method

In [21], Reddy, Mahanta, and Bora have modified the timing metric to reduce the plateau, however the single peak will be affected by noises, so the timing estimation is still uncertain. The new timing metric K(d) is shown as:
K ( d ) = M ( d ) M ( d L )
where K(d) will reach the maximum value at the start of training sequence and fall subsequently, thus, there is no plateau. As shown in Figure 3, the M(d) of SC method and the K ( d ) of RMB method are simulated in an AWGN channel with Eb/N0 (the energy per bit to noise power spectral density ratio) = 30 dB, time delay = 191 samples and L = 16.
Figure 3. M(d) of SC method and K(d) of RMB method.
Figure 3. M(d) of SC method and K(d) of RMB method.
Sensors 15 28942 g003

3.3. MathWorks Method

The SC and RMB methods use auto-correlation, but in the method of Matlab [22], MathWorks employs cross-correlation between the received signal and the short preamble given in Equation (2). The timing metric M(d) is the same as the SC method in Equation (3), but P(d) is defined as:
P ( d ) = m = 0 L 1 r ( d + m ) · t * ( 1 + m )
where t(1 + m) is the short preamble sequence, i.e., the reference template, and R(d) is the received energy defined as:
R ( d ) = m = 0 T 1 | r ( d + m ) | 2    
with length T = 16. Because of the repeat of the short preamble samples, there will be nine peaks in the timing metric M(d). In the MathWorks method, during the timing period, the value of 60% × maximum of M(d) is used as a threshold to determine the time offset. If the peaks exceeding the threshold are located in the spacing of 16 times (16, 32, 48,…, 128) and the count of these peaks is more than 5, then the location of the first peak ( L f i r s t ) is deemed as the peak of the first path, and the timing delay ( D ^ est ) is estimated as:
D ^ est = L f i r s t ( L 2 + 1 )    
Two timing metric examples with Eb/N0 = 0 dB and 30 dB are simulated in Figure 4 and the time delay is 191 samples. In these cases, the first peak should be located at L f i r s t =   191 + 16/2 + 1 = 200th sample. As shown in Figure 4, When Eb/N0 = 30 dB, the nine peaks are very sharp and the first peak appears at the 200th sample. When Eb/N0 = 30 dB, all nine peaks exceed the threshold, but when Eb/N0 = 0 dB, there are only two peaks exceeding the threshold, so the timing offset cannot be estimated accurately. In order to solve this problem, the following two steps are put forward in the following Section.
Figure 4. M(d) of MathWorks method in different Eb/N0.
Figure 4. M(d) of MathWorks method in different Eb/N0.
Sensors 15 28942 g004

3.4. Liu et al. Method

Auto-correlation based methods, such as the SC method or RMB method, are immune to carrier frequency offset (CFO), but their timing metric is not sharp enough. On the other hand, a cross-correlation based method such as the MathWorks method has sharp timing metrics, but its performance is sensitive to CFO. In [23], Liu et al. proposed a new timing method which used both auto-correlation and cross-correlation, so this method has the advantages of both approaches. Figure 5 shows the timing metric M(d) for the Liu et al. method in an ideal channel with a time delay of 191 samples. Because of the 10 identical short preamble symbols in IEEE 802.11p, there will be 10 peaks spaced 16 samples apart in an ideal environment, which may lead to uncertainty in the timing estimation just like the MathWorks method.
Figure 5. M(d) of Liu method in ideal environment.
Figure 5. M(d) of Liu method in ideal environment.
Sensors 15 28942 g005

4. Proposed Timing Estimation Method

In this paper, we will pay attention to the correlation receive, that is to say, the received signal is correlated with the reference template (the short preamble symbols), which is given by:
r x s ( τ d ) = 1 T 0 T 0 p ( t ) r ( t τ d ) d t
where T0 is the correlation duration, p(t) is the reference template, the short preamble symbols, and r(t) is the received signal.
The receiver can estimate the time offset according to the output samples from the correlator. Then, the range between the transmitter and receiver nodes can be calculated according to this estimated time offset. Several techniques have been developed for OFDM-based timing estimation. If the reference template signal is selected from the received signal, this method is called auto-correlation, otherwise it is a cross-correlation. The main steps of the proposed method contain cross-correlation, summing-up, and skewness-analysis, so it is denoted as the cross-correlation, summing-up, and skewness-analysis (CSS) method.

4.1. Cross-Correlation

The correlation is the same as in the MathWorks method which includes the Equations (7), (8), and (3). As shown in Figure 4, in an ideal environment, there will be nine peaks in the timing metric M(d) and the spacing between each adjacent desired peak is 16 samples. However, if Eb/N0 is low, the noise can affect the number of peaks that exceed a specific threshold, and it is difficult to determine the optimal threshold value, this is why the performance of the MathWorks method is relatively low.

4.2. Summing-Up

The summing-up of nine correlated samples M(d) with a pace of T = 16 samples apart are considered, as shown in Equation (11):
G ( t ) = d = 1 9 M ( t + ( d 1 ) * L )
and the sample index ( D ^ est ) responding to the time delay is obtained as:
D ^ est = argmin ( t | G [ t ] ξ * max ( G ( t ) ) ) ( T 2 + 1 )    
where ξ is a threshold between 0 and 1.
In ideal environment with no noise, the new timing metric G(t) is simulated in Figure 6 with time delay = 191 samples. Now, there are 17 symmetrical peaks and the location of the maximum peak is just the index to be found, which should be located at the 191 + 16/2 + 1 = 200th sample.
Figure 6. M(d) and G(t) in ideal environment.
Figure 6. M(d) and G(t) in ideal environment.
Sensors 15 28942 g006
As Eb/N0 decreases, the timing metric M(d) of MathWorks method will become more and more hard to be identified, but the new timing metric G(t) is better. As shown in Figure 7, when Eb/N0 = −5 dB, the 9 peaks of M(d) are very hard to be identified, but the maximum of G(t) is very clear, and is located at 200th sample, and then can be estimated correctly.
Figure 7. M(d) and G(t) with Eb/N0 = −5 dB.
Figure 7. M(d) and G(t) with Eb/N0 = −5 dB.
Sensors 15 28942 g007
However, as Eb/N0 decreases more and more, the G(t) will become more difficult to identify by a fixed threshold. As shown in Figure 8, when Eb/N0 = −10 dB, the maximum of G(t) is very hard to locate in the 200th sample. That is to say, the maximum of G(t) is not the proper sample index. At the same time, it is very hard to determine the threshold ξ to find the proper index. Therefore, in the following part of this Section, a dynamic threshold method is put forward.
Figure 8. G(t) with Eb/N0 = −10 dB.
Figure 8. G(t) with Eb/N0 = −10 dB.
Sensors 15 28942 g008

4.3. Skewness-Analysis and Threshold

In order to improve the precision in low Eb/N0, a dynamic threshold using the skewness analysis of IEEE 802.11p short preamble is proposed. In this Section, the skewness and standard deviation (STD) of the G(t) are analyzed. Then, in order to determine the best threshold (ξbest), the relationship between estimation error and the skewness of G(t) was investigated. Finally the relationship between skewness and dynamic threshold is set up to estimate TOA.

4.3.1 Statistical Characteristics of the Summing-Up Samples

(1)
Standard Deviation
STD is used to measure how much variation there is from the average of G(t), which is given by:
δ = i = 1 N b ( x i x ¯ ) 2 ( N b 1 )
where Nb is the number of G(t) samples and xi is the ith value of G(t).
(2)
Skewness
The skewness of G(t) is given by:
S = 1 ( N b 1 ) δ 3 i = 1 N b ( x i x ¯ ) 3
where x ¯ is the mean of G(t), and δ is the standard deviation of G(t). If G(t) is a normal distribution, S will be zero. If S < 0, G(t) will skew left which indicates that the left tail is longer than the right one, while if S > 0, G(t) will skew right and the left tail is shorter than the right one. As the Eb/N0 increases, S will tend to increase. If there are no signals (or Eb/N0 is too low), S will be zero.
In order to examine the characteristics of the two statistical parameters, the ITU channel model is employed. For each Eb/N0 from −10 dB to 20 dB, 2000 channel realizations were generated. The average results of STD and skewness are normalized and shown in Figure 9 which showa that the skewness increases as the Eb/N0 increases but the STD increases at first and then decreases as the Eb/N0 increases. Therefore, skewness is a monotonic function for a large range of Eb/N0 values, and it is more suitable for TOA estimation than STD, because it can better reflect changes in Eb/N0.
Figure 9. Two normalized statistical parameters change with Eb/N0.
Figure 9. Two normalized statistical parameters change with Eb/N0.
Sensors 15 28942 g009

4.3.2. Relationship between Estimation Error, Skewness and Threshold

In order to determine the best threshold (ξbest) based on skewness, the relationship between estimation error, skewness and threshold was investigated. 2000 ITU channel realizations with each Eb/N0 = {−10, 5,…, 20} dB were simulated. With each channel realizations, the thresholds of {0.05, 0.10, 0.15, 0.20,…,1.0} are compared with G(t) to find the first threshold crossing sample index, as shown in Equation (12).
To illustrate the results, Figure 10 shows the estimation error for S = {0.1, 0.8, 1.8} in the ITU channels respectively. The relationship is that for a specific S, as the thresholds increasing, the estimation error decreases at first and then increases. For example, the threshold responding to the lowest errors for S = {0.1, 0.8, 1.8} are 0.85, 0.9 and 0.95. Thus, the threshold with respect to the minimum estimation error is selected as the best threshold ξbest.
Figure 10. Relationship between MAE, threshold and skewness.
Figure 10. Relationship between MAE, threshold and skewness.
Sensors 15 28942 g010

4.3.3. Relationship between Skewness and Threshold

A polynomial with degree 3 was fitted to the best threshold ξbest for each value of S by using the method of least-squares where S is the x-coordinate and ξbest is the y-coordinate. The result of the fitting will generate the coefficients of the polynomial with the minimized summed square of residuals. The ith residual ri for the ith pair of (S, ξbest) is defined as:
r i = y i y ^ i
where yi is the best threshold and ŷi is the fitted threshold value for the ith S, so the summed square of residuals SS is given by:
S s = i = 1 n r i 2 = i = 1 n ( y i y ^ i ) 2
where n is the number of (S, ξbest).
The least-squares fitting result is shown in Equation (17) and Figure 11:
ξ = 0.016085S3−0.088866S2 + 0.17697S + 0.82368
Figure 11. Relationship between skewness and threshold.
Figure 11. Relationship between skewness and threshold.
Sensors 15 28942 g011

5. Simulation Setup and Result Analysis

5.1. Simulation Setup

5.1.1. System Model

Because vehicles are typically in motion, there will be a Doppler shift given by:
F d = v * f c / c
where v is the relative vehicle speed which is typically set to v = 100 km/h, fc = 5.9 GHz is the center frequency of the transmitted signal, and c is the speed of light.
The SC method with L = 32 has a plateau as shown in Figure 2, so the timing estimation is difficult. To mitigate this problem, Ng should equal to 0, so L = 80 is used (denoted by SC-80). The timing metric of SC-80 is presented in Figure 12 for an ideal environment with a time delay of 191 samples, and this indicates that there is no plateau. The other parameters for the simulation are shown in Table 1.
Figure 12. |P(d)|2, R2(d) and M(d) of SC-80 method in ideal environment.
Figure 12. |P(d)|2, R2(d) and M(d) of SC-80 method in ideal environment.
Sensors 15 28942 g012
Table 1. Parameters for the simulation.
Table 1. Parameters for the simulation.
ParametersValue
Modulation modeBPSK
Number of subcarriers52
Symbol duration8 μs
Guard time1.6 μs
FFT period6.4 μs
Preamble duration32 μs
Subcarrier spacing0.15625 MHz
Vehicle speed100 km/h
Channel modeITU-A
The block of coarse time offset estimation using a short preamble is important not only to the ranging of vehicle nodes but also the data communication. Only if the coarse time offset is relatively accurate can the following blocks be correct. In this article we just pay attention to the time offset estimation method using the short preamble in IEEE 802.11p.

5.1.2. Transmission Channel

In an actual application environment, the sent signals may be scattered, reflected, or diffracted by surrounding objects, which may cause some different paths and attenuation. The International Telecommunications Union (ITU) vehicular test environment [24] is a well-known channel model, which can be used to model the multipath delays and channel gains. There are two different delay spreads specified, that is, low delay spread (A), and medium delay spread (B) which are shown in Table 2. The simulation in this paper uses the Channel A.
Table 2. ITU’s vehicular test environment parameters.
Table 2. ITU’s vehicular test environment parameters.
TapChannel AChannel BDoppler
Relative Delay (ns)Average Power (dB)Relative Delay (ns)Average Power (dB)Spectrum
100.00–2.5Classic
2310–1.03000Classic
3710–9.08900–12.8Classic
41090–10.012,900–10.0Classic
51730–15.017,100–25.2Classic
62510–20.020,000–16.0Classic

5.1.3. Performance Metric

In this section, the timing estimation performances are examined in the ITU Channel A parameters. The performance metrics are mean absolute errors (MAE) and the correct percentage of timing estimation error = 0:
M A E = 1 K n = 1 K | D n D ^ n |
where Dn is the peak index corresponding to the nth actual time delay, D ^ n is the peak index corresponding to the nth estimated time delay, and K is the number of simulation iterations.

5.2. Performance Results and Analysis

The performance in ITU channel was simulated K = 2000 times for values of Eb/N0 from −15 dB to 20 dB. In the remainder of this paper, six methods are considered, those are the CSS method with dynamic threshold (denoted by CSS-dynamic), the CSS method with fixed threshold = 1 (denoted by CSS-fixed), the MathWorks method, the SC-80 method, the RMB method, and the Liu et al. method (denoted by Liu).

5.2.1. Estimation Error

Figure 13 presents MAE for the six methods. These results show that the CSS-dynamic method provides the best performance, particularly at low Eb/N0 values.
Figure 13. MAE of different methods.
Figure 13. MAE of different methods.
Sensors 15 28942 g013
For example, when Eb/N0 = −10 dB, the MAE of the CSS-dynamic method is 66 samples, while MAE are only 101 samples for the CSS-fixed method, 102 samples for the SC-80 method, 136 samples for the RMB method, 162 samples for the Liu method, and 192 samples for the MathWorks method. The MAE of the MathWorks method is the highest when Eb/N0 is lower than 9 dB. When Eb/N0 is less than −2 dB the CSS-fixed method is almost the same as the SC-80 method. As shown in Figure 14, because of the multipath, when Eb/N0 is greater than 13 dB, the MAEs of the CSS-dynamic, CSS-fixed, and MathWorks methods are close to 0. Nevertheless, the MAE of the SC-80 method is close to three samples and the MAE of RMB is close to four samples. For the Liu method when Eb/N0 is greater than 18 dB, so its MAEs drop down to two samples.
Figure 14. MAE of different methods with Eb/N0 > 10 dB.
Figure 14. MAE of different methods with Eb/N0 > 10 dB.
Sensors 15 28942 g014

5.2.2. Correct Percentage

Figure 15 presents the percentage of timing estimation error = 0 for the six methods. These results also show that the CSS-dynamic and CSS-fixed methods provide the best performance. For example, when Eb/N0 = 0 dB, the percentage of the CSS-dynamic and CSS-fixed methods having no error are nearly 59%, whereas for the same percentage with MathWorks, Eb/N0 should exceed 8 dB, but for the RMB, Liu and SC-80 methods the percentages are too low for all values of Eb/N0. When Eb/N0 is greater than 14 dB, the percentages of the CSS-dynamic, CSS-fixed and MathWorks methods are close to 100%, but the others are very low which can also be explained using Figure 14. This is because when Eb/N0 is greater than 14 dB the MAE of the other three methods are not close to 0 samples. At the same time, for the methods of SC-80 and RMB, the percentage increases at first as Eb/N0 increases, and then decreases to three or four samples. The reason for the behavior can be explained using Figure 16 and Figure 17 which show the performance in a channel with no noise.
(1)
When Eb/N0 is high, the multipath in the ITU channel causes the peaks to move to the right, making the timing estimation incorrect. Take the SC-80 method as an example, which is illustrated in Figure 16. The percentage will be close to 0 but will not reach 0. On the other hand, Figure 17 gives the corresponding results for the CSS method. This shows that the maximum peak with the CSS method does not move in the multipath channel, and only some smaller peaks are generated, so the percentage will be nearly equal to 100% for a high Eb/N0.
(2)
When Eb/N0 is very low (close to −10 dB), the noise levels are too high, so it is too difficult to locate the peak, so the percentage is again close to 0.
(3)
For Eb/N0 around 0 dB, the signal energy is close to that of the noise, so it is easier to locate the peak. Now because of the randomization of noise, it’s likely that the estimated time happens to be the true time delayed.
Figure 15. Percentage that the timing error is zero in ITU channel.
Figure 15. Percentage that the timing error is zero in ITU channel.
Sensors 15 28942 g015
Figure 16. The effect of multipath channel on the timing metric of the SC-80 method.
Figure 16. The effect of multipath channel on the timing metric of the SC-80 method.
Sensors 15 28942 g016
Figure 17. The effect of multipath channel on the timing metric of the CSS method.
Figure 17. The effect of multipath channel on the timing metric of the CSS method.
Sensors 15 28942 g017

6. Conclusions

In this paper, a high precise ranging estimation method based-on the short preamble has been developed for vehicular node positioning using the IEEE802.11p standard. The correlation-summing-skewness analysis method with a dynamic threshold was developed for timing estimation which may be a valuable reference for the study of positioning in vehicular wireless networks. The performance of the CSS-dynamic method is shown to be the best among several well-known methods such as the SC method, MathWorks method, and RMB method, especially in low SNRs, and multi-path environments which is very important for vehicular wireless networks.
The CSS method can also be used in other wireless communications systems which use OFDM with repeated preamble symbols, such as IEEE 802.11a. On the basis of this coarse time estimation method, we will study on how to use this method in the following fine timing estimation and frequency shift estimation and then the high ranging or position method can be used in the application of cooperative collision avoidance to reduce traffic accidents.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant No. 61301139, the Nature Science Foundation of Shandong Province No. ZR2014FL014 and ZR2014FM017, Distinguished Middle-Aged and Young Scientist Encourage and Reward Foundation of Shandong Province No. 2014BSE28032, the Science and Technology Project in Colleges and Universities of Shandong Province No. J14LN53, and the Project of Basic Research Application of Qingdao City No. 14-2-4-83-jch.

Author Contributions

Xuerong Cui conceived and designed the algorithm, and wrote the paper; Juan Li designed and performed the simulations. Chunlei Wu and Jian-Hang Liu analyzed the data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. Global Status Report on Road Safety 2013: Supporting a Decade of Action. Available online: http://www.who.int/iris/bitstream/10665/78256/1/9789241564564_eng.pdf (accessed on 30 August 2015).
  2. Xiao, Z.; Bai, J.; Ma, G.Y.; Fan, J.T.; Yi, K.C. Research on positioning enhancement scheme of CAPS via UWB pseudolite. Sci. China Phys. Mech. Astron. 2012, 55, 733–737. [Google Scholar] [CrossRef]
  3. Ansari, K.; Feng, Y.M.; Tang, M.L. A runtime integrity monitoring framework for real-time relative positioning systems based on GPS and DSRC. IEEE Trans. Intell. Transp. Syst. 2015, 16, 980–992. [Google Scholar] [CrossRef]
  4. IEEE 802.11p-2010 Standard. Available online: http://standards.ieee.org/getieee802/download/802.11p-2010.pdf (accessed on 30 August 2015).
  5. Haijun, Z.; Chunxiao, J.; Julian, C. Cooperative interference mitigation and handover management for heterogeneous cloud small cell networks. IEEE Wirel. Commun. 2015, 22, 92–99. [Google Scholar] [CrossRef]
  6. Wei, W.; Qi, Y. Information potential fields navigation in wireless Ad-Hoc sensor networks. Sensors 2011, 11, 4794–4807. [Google Scholar] [CrossRef] [PubMed]
  7. Sanchez-Rodriguez, D.; Hernandez-Morera, P.; Quinteiro, J.M.; Alonso-Gonzalez, I. A low complexity system based on multiple weighted decision trees for indoor localization. Sensors 2015, 15, 14809–14829. [Google Scholar] [CrossRef] [PubMed]
  8. Vera, R.; Ochoa, S.F.; Aldunate, R.G. EDIPS: an easy to deploy indoor positioning system to support loosely coupled mobile work. Pers. Ubiquit. Comput. 2011, 15, 365–376. [Google Scholar] [CrossRef]
  9. IEEE 802.11a-1999 Standard. Available online: http://standards.ieee.org/getieee802/download/802.11a-1999.pdf (accessed on 30 August 2015).
  10. Cui, X.-R.; Zhang, H.; Aaron Gulliver, T. Threshold selection for ultra-wideband TOA estimation based on neural networks. J. Netw. 2012, 7, 1311–1318. [Google Scholar] [CrossRef]
  11. Sun, Y.L.; Xu, Y.B.; Li, C.; Ma, L. Kalman/map filtering-aided fast normalized cross correlation-based Wi-Fi fingerprinting location sensing. Sensors 2013, 13, 15513–15531. [Google Scholar] [CrossRef] [PubMed]
  12. Tomic, S.; Beko, M.; Dinis, R. Distributed RSS-based localization in wireless sensor networks based on second-order cone programming. Sensors 2014, 14, 18410–18432. [Google Scholar] [CrossRef] [PubMed]
  13. Miao, C.Y.; Dai, G.Y.; Ying, K.Z.; Chen, Q.Z. Collaborative localization and location verification in WSNs. Sensors 2015, 15, 10631–10649. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, Q.; Foh, C.H.; Seet, B.C.; Fong, A.C.M. Variable elasticity spring-relaxation: improving the accuracy of localization for WSNs with unknown path loss exponent. Pers. Ubiquit.Comput. 2012, 16, 929–941. [Google Scholar] [CrossRef]
  15. Adnan, T.; Datta, S.; MacLean, S. Efficient and accurate sensor network localization. Pers. Ubiquit. Comput. 2014, 18, 821–833. [Google Scholar] [CrossRef]
  16. Chen, Z.H.; Zou, H.; Jiang, H.; Zhu, Q.C.; Soh, Y.C.; Xie, L.H. Fusion of WiFi, smartphone sensors and landmarks using the Kalman filter for indoor localization. Sensors 2015, 15, 715–732. [Google Scholar] [CrossRef] [PubMed]
  17. Guo, J.Q.; Zhang, H.Y.; Sun, Y.C.; Bie, R.F. Square-root unscented Kalman filtering-based localization and tracking in the Internet of Things. Pers. Ubiquit. Comput. 2014, 18, 987–996. [Google Scholar] [CrossRef]
  18. Brida, P.; Machaj, J.; Benikovsky, J. A modular localization system as a positioning service for road transport. Sensors 2014, 14, 20274–20296. [Google Scholar] [CrossRef] [PubMed]
  19. Schmidl, T.M.; Cox, D.C. Robust frequency and timing synchronization for OFDM. IEEE Trans. Commun. 1997, 45, 1613–1621. [Google Scholar] [CrossRef]
  20. Oni, R.R.A. Preamble design problematic with 802.11a IEEE standard (Minn’s training sequence). Radioengineering 2008, 17, 87–91. [Google Scholar]
  21. Reddy, A.K.; Mahanta, A.A.; Bora, P.K. On timing and frequency offset estimation in OFDM systems. In Proceedings of the TENCON 2004. 2004 IEEE Region 10 Conference, 21–24 November 2004; Volume 153, pp. 153–156.
  22. MathWorks. OFDM Synchronization. Available online: http://www.mathworks.com/help/comm/examples/ofdm-synchronization.html (accessed on 30 August 2015).
  23. Yun, L.; Hua, Y.; Fei, J.; Fangjiong, C.; Weiqiang, P. Robust timing estimation method for OFDM systems with reduced complexity. IEEE Commun. Lett. 2014, 18, 1959–1962. [Google Scholar] [CrossRef]
  24. Guidelines for Evaluation of Radio Transmission Technologies for IMT-2000. Available online: https://www.itu.int/dms_pub/itu-r/oth/0A/0E/R0A0E00000C0001MSWE.doc (accessed on 30 August 2015).

Share and Cite

MDPI and ACS Style

Cui, X.; Li, J.; Wu, C.; Liu, J.-H. A Timing Estimation Method Based-on Skewness Analysis in Vehicular Wireless Networks. Sensors 2015, 15, 28942-28959. https://doi.org/10.3390/s151128942

AMA Style

Cui X, Li J, Wu C, Liu J-H. A Timing Estimation Method Based-on Skewness Analysis in Vehicular Wireless Networks. Sensors. 2015; 15(11):28942-28959. https://doi.org/10.3390/s151128942

Chicago/Turabian Style

Cui, Xuerong, Juan Li, Chunlei Wu, and Jian-Hang Liu. 2015. "A Timing Estimation Method Based-on Skewness Analysis in Vehicular Wireless Networks" Sensors 15, no. 11: 28942-28959. https://doi.org/10.3390/s151128942

APA Style

Cui, X., Li, J., Wu, C., & Liu, J. -H. (2015). A Timing Estimation Method Based-on Skewness Analysis in Vehicular Wireless Networks. Sensors, 15(11), 28942-28959. https://doi.org/10.3390/s151128942

Article Metrics

Back to TopTop