Next Article in Journal
Where There Is Fire There Is SMOKE: A Scalable Edge Computing Framework for Early Fire Detection
Next Article in Special Issue
Fingerprints and Floor Plans Construction for Indoor Localisation Based on Crowdsourcing
Previous Article in Journal
Estimating the Volume of Unknown Inclusions in an Electrically Conducting Body with Voltage Measurements
Previous Article in Special Issue
Lightweight Workload Fingerprinting Localization Using Affinity Propagation Clustering and Gaussian Process Regression
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Virtual Antenna Array and Fractional Fourier Transform-Based TOA Estimation for Wireless Positioning

School of Electronic and Information Engineering, Xi’an Jiaotong University, No. 28 West Xianning Road, Xi’an 710049, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(3), 638; https://doi.org/10.3390/s19030638
Submission received: 7 December 2018 / Revised: 21 January 2019 / Accepted: 30 January 2019 / Published: 2 February 2019
(This article belongs to the Collection Positioning and Navigation)

Abstract

:
In this paper, a novel virtual antenna array and fractional Fourier transform (FRFT)-based 2-dimension super-resolution time-of-arrival (TOA) estimation algorithm for OFDM WLAN systems has been proposed. The proposed algorithm employs channel frequency responses (CFRs) at the equi-spaced positions on a line or quasi-line moving trajectory, i.e., the CFRs of a virtual antenna array, to extract multipaths’ TOA information. Meanwhile, a new chirp-like quadratic function is used to approximate the channel multipaths’ phase variation across the space dimension, which is more reasonable than the traditional linear function, especially for relatively big virtual antenna array sizes. By exploiting the property of chirp-like multipaths’ energy concentration in the FRFT domain, the FRFT can be first used to separate chirp-like multipath components, then the existing TOA estimation methods in frequency domain can be further employed on the separated multipath components to obtain the multipaths’ TOA estimates. Therefore, the proposed algorithm can make more use of the multipaths’ characteristics in the space dimension, thus it can efficiently enhance the multipath resolution and achieve better multipaths’ TOA estimation performance without requiring a real antenna array. Simulation results demonstrate the effectiveness of the proposed algorithm.

1. Introduction

Indoor wireless positioning technology has increasingly attracted research interests [1,2] for its wide and popular applications in the construction industry, health industry [3], people guidance and so on. Geometric location based on time-of-arrival (TOA) is the most popular method for accurate positioning systems [4]. The basic problem in TOA estimation is to extract the TOA of the Line-of-Sight (LOS) path in wireless channels, which is usually interfered by Non-LOS (NLOS) paths in the multipath environment [5]. Due to the serious NLOS multipath effects and the limited system bandwidth, the conventional TOA estimation methods, such as the inverse Fourier transform (IFT)-based method, cannot achieve higher estimation accuracy [6].
To improve the TOA estimation accuracy, several super-resolution TOA estimation methods based on the channel frequency responses (CFRs), such as the multiple signal classification (MUSIC)-based method [1], TLS-ESPRIT [7]-based method, matrix pencil (MP)-based method [8,9] and so on, employ the eigenstructure or eigen-subspace of the CFRs to obtain TOA estimates. In contrast to the conventional methods, these super-resolution methods exploit not only the property of multipaths’ linear phase variation across the frequency dimension but also the sparsity of the main channel multipaths in the time-delay domain, so they can achieve higher accuracy than the conventional methods. Nevertheless, such super-resolution methods are susceptible to the noise since the noise causes serious estimation errors of eigen-subspace or eigenvalues.
Furthermore, the 2-dimension super-resolution TOA estimation methods, such as the 2-dimension MP method [10], joint beamforming and MUSIC-based method [11], and the 2-dimension (2D) MUSIC method [12], employ both the channel frequency responses and the channel space responses to jointly estimate the TOA and AOA. Besides the property of multipaths’ linear phase differences across the frequency dimension and the sparsity of the main channel multipaths, the 2D super-resolution methods also exploits the property of multipaths’ quasi-linear phase variation across the space dimension, hence they can efficiently mitigate the noise effects and further improve the TOA estimation performance. However, such 2D super-resolution methods are limited by the requirement of an antenna array to obtain the channel space responses, and the assumption of quasi-linear phase variation property across the space dimension, which implies that the antenna array size should be small enough.
Without requiring an antenna array and suffering from the limitation of its relatively small size, a novel virtual antenna array and fractional Fourier transform (FRFT)-based 2D super-resolution TOA estimation method has been proposed in this paper. By using the built-in Kinect or inertial sensors in the mobile terminal (MT), the proposed method first obtains the CFRs at the equi-spaced positions on a line or quasi-line moving trajectory as the CFRs of a virtual antenna array. Meanwhile, considering that the size of virtual antenna array is not limited by the hardware and not precisely known, the channel multipaths’ phase variation across the space dimension is approximated by a chirp-like quadratic function, which is more reasonable than a traditional linear function. Based on the CFRs of a virtual antenna array, the proposed method employs the FRFT [13] to separate and extract multipath components in CFRs by exploiting the property of the chirp’s energy concentration in the FRFT domain [14], then it uses the existing time-delay estimation methods, such as the IFT method and the MP method [9], to estimate the TOAs of the separated multipaths. Therefore, by employing the virtual antenna array idea, modeling the multipath phase variation as a quadratic function, the proposed FRFT-based algorithm can make more use of the multipaths’ characteristics in the space domain and achieve more robust TOA estimation performance without requiring an antenna array.
The main contributions of this paper are as follows: (1) We present a novel virtual antenna array system model for positioning by regarding the CFRs at the equi-spaced positions on a line or quasi-line moving trajectory, which can be obtained by using built-in Kinect or inertial sensors in the MT, as CFRs of a virtual antenna array. (2) The channel multipaths’ phase variation across the space dimension is modeled as a chirp-like quadratic function instead of a traditional linear function. Due to being high-order approximation, the chirp-like quadratic function suffers from much less multipath phase model mismatch error for the unlimited and unspecified size of virtual antenna arrays. (3) An FRFT-based 2D super-resolution TOA estimation method has been proposed. By exploiting the chirp’s energy concentration property in the FRFT domain, the FRFT is employed to separate chirp-like multipath components in the space domain, then the existing TOA estimation methods, such as the IFT method and the MP method [9], are used to obtain multipaths’ TOA estimate. Thus, the FRFT-IFT algorithm and the FRFT-MP algorithm have been developed. (4) Simulation results demonstrate that the proposed method can achieve more robust TOA estimation performance than the 2D MP super-resolution TOA estimation method due to making more of the multipaths’ characteristics in the space domain.
The rest of the paper is organized as follows. In Section 2, the virtual antenna array system model is described. In Section 3, the properties of multipath components in the CFRs of a virtual antenna array will be analyzed first by FRFT, and a novel FRFT-based TOA estimation algorithm will be then developed based on these properties. Simulation results are presented in Section 4 to evaluate the performance of the proposed method. Finally, conclusions are drawn in Section 5.

2. Virtual Antenna Array System Model

For ease of exposition, we consider a basic OFDM WLAN system with single AP (Access Point) and single MT. The following system model can be extended directly to the WLAN system with multiple APs and multiple MTs. Suppose that the MT moves in a straight-line or quasi-straight-line during a short period, which is usually the case with pedestrians‘ walking way in indoor environments. If the MT’s acceleration, which can be measured by built-in Kinect or inertial sensors in the MT [15], is smaller than a chosen threshold, the MT can be assumed to make uniform linear motion with its speed imprecisely estimated due to the accumulated error. By simply taking positions at M equi-interval moments during such short period of relatively small acceleration or even employing the pedestrian dead reckoning method [16,17], we can obtain M equi-spaced positions on a straight-line or quasi-straight-line moving trajectory of the MT and regard the CFRs at such positions as the CFRs of a virtual antenna array. Suppose the central position of the trajectory as the position of interest and call the CFRs at the M positions on the trajectory as the CFRs of a virtual antenna array at this target position. Denote the distance between adjacent positions as ρ , i.e., the virtual antenna array element spacing, which is assumed constant but not precisely estimated due to the imprecisely estimated speed and the negligible incremental accumulated error in a short period. The channel impulse response (CIR) between the AP and the MT at the m-th position of the moving trajectory is given by [9]
h ( t , m ) = l = 1 L b l a ( τ m , l ) δ ( t τ m , l ) m = M 2 , , 0 , , M 2 1
where L is the number of distinct propagation paths, b l and τ m , l are the complex gain and the TOA of the l-th path, respectively, a ( τ m , l ) = e j 2 π f c τ m , l is the array response of the mth antenna to the l-th path.
Assume the direction of arrival (DOA) relative to the moving trajectory of the lth path at the central position, i.e., the 0-th position on the trajectory as θ l , which is shown in Figure 1, we have
τ m , l = 1 c ( m ρ ) 2 + ( c τ 0 , l ) 2 2 c o s θ l ( m ρ ) ( c τ 0 , l )    1 c ( c τ 0 , l ( m ρ ) c o s θ l + 1 2 ( m ρ s i n θ l ) 2 c τ 0 , l ) m = M 2 , , 0 , , M 2 1
where c is the speed of light, the approximation holds if the virtual antenna array size M ρ is deliberately chosen to satisfy M ρ 2 c τ 0 , l 1 . This approximation is obtained by using second order Taylor series expansion, thus it is more reasonable than the conventional linear function, i.e., the first order Taylor series expansion of τ m , l , considering the size of the virtual antenna array is unspecified and not limited by the hardware.
From (1) and (2), the estimated CFR for the nth subcarrier and the m-th position on the moving trajectory at the receiver side can be expressed as
H ^ m , n = H m , n + W m , n = l = 1 L b l e j ω c τ m , l e j 2 π Δ f τ m , l + W m , n l = 1 L { b l e j ω c ( τ 0 , l m ρ c o s θ l / c + ( m ρ s i n θ l ) 2 / ( 2 τ 0 , l c 2 ) ) · e j 2 π n Δ f ( τ 0 , l m ρ c o s θ l / c + ( m ρ s i n θ l ) 2 / ( 2 τ 0 , l c 2 ) ) } + W m , n m = M 2 , , M 2 1 ; n = 0 , 1 , , N
as in [10], where ω c = 2 π f c is the carrier angular frequency, Δ f , N and W m , n denote the OFDM subcarrier spacing, the subcarrier number and the additive white Gaussian noise (AWGN), respectively. For simplicity, the CFR can be also written as
H ^ m , n l = 1 L { b l e j φ l ( n ) + j 2 π f l ( n ) m ρ + j π μ l ( n ) ( m ρ ) 2 } + W m , n
where
φ l ( n ) = ( ω c + 2 π n Δ f ) τ 0 , l f l ( n ) = ( f c + n Δ f ) c o s θ l / c μ l ( n ) = ( f c + n Δ f ) ( s i n θ l ) 2 / ( τ 0 , l c 2 )

3. Fractional Fourier Transform Based TOA Estimation

It can be observed from (4) that the phases of multipath components in CFRs { H m , n } m = 1 , n = 1 M , N across the frequency dimension vary linearly, while the phase variations of multipath components in CFRs across the space dimension approximate quadratic functions, in other words, the CFRs over the space dimension can be equivalently viewed as L superimposed chirp-like multipath components. Consequently, the FRFT in the space domain and the frequency domain processing method, such as the IFT or MP method, can be jointly employed to extract multipath components’ parameters.
In this section, the multipath components of the CFRs will be analyzed first by employing FRFT across the space dimension, and a novel FRFT-based multipath TOA estimation algorithm will be then developed.

3.1. Properties of Multipath Components in FRFT Domain

To separate and extract the chirp-like components, a typical discrete fractional Fourier transform (DFRFT) [18] will be implemented on the CFRs at each subcarrier across the space dimension as
H ^ n ( α , u ) = F α ( { H m , n } m = M 2 M 2 1 ) + F α ( { W m , n } m = M 2 M 2 1 ) = A α M e j π u 2 c o t α m = M 2 M 2 1 { e j 2 π u c s c ( α ) / M e j π u 2 c o t ( α ) / M H m , n } + W n ( α , u ) n = 0 , 1 , , N
where α = p π / 2 and p is the FRFT order, A α = 1 j c o t ( α ) 2 π , F α ( ) denotes the DFRFT operation of order α = p π / 2 , W n ( α , u ) is the transformed AWGN by DFRFT.
By substituting (4) into the first part of (6), the noise-free transformed CFRs by DFRFT can be expressed as
H n ( α , u ) A α M e j π u 2 c o t α m = M 2 M 2 1 { e j 2 π u m c s c ( α ) / M e j π m 2 c o t ( α ) / M · l = 1 L { b l e j φ l ( n ) + j 2 π f l ( n ) m ρ + j π μ l ( n ) ( m ρ ) 2 } } = A α M e j π u 2 c o t α l = 1 L { b l e j φ l ( n ) m = M 2 M 2 1 { e j 2 π [ u c s c ( α ) / M + f l ( n ) ρ ] m · e j π [ c o t ( α ) / M + μ l ( n ) ρ 2 ] m 2 } } n = 0 , 1 , , N
Define the DFRFT signal in (7) contributed only by the l-th multipath component as S l , n ( α , u )
S l , n ( α , u ) = c l e j φ l ( n ) m = M 2 M 2 1 { e j 2 π [ u c s c ( α ) / M + f l ( n ) ρ ] m · e j π [ c o t ( α ) / M + μ l ( n ) ρ 2 ] m 2 } n = 0 , 1 , , N
where
c l = A α M e j π u 2 c o t α b l
Then it can be derived that S l , n ( α , u ) possesses the following properties.
Property 1.
Due to the chirp’s energy concentration property in the DFRFT domain [14], S l , n ( α , u ) has a distinguishable peak in the ( α , u ) plane, the peak, and its coordinate ( α ´ l ( n ) , u ´ l ( n ) ) correspond to
S l , n ( α ´ l ( n ) , u ´ l ( n ) ) A α ´ l ( n ) M e j π ( u ´ l ( n ) ) 2 c o t ( α ´ l ( n ) ) b l e j φ l ( n ) α ´ l ( n ) = a r c c o t ( μ l ( n ) M ρ 2 ) u ´ l ( n ) = f l ( n ) M ρ / c s c ( α ^ l ( n ) )
Property 2.
It can be proved that n = 1 N | S l , n ( α , u ) | 2 achieves its maximum value at ( α ´ l ( N 2 ) , u ´ l ( N 2 ) ) , i.e.,
( α ´ l ( N 2 ) , u ´ l ( N 2 ) ) arg max ( α , u ) n = 1 N | S l , n ( α , u ) | 2
and decreases rapidly with | α α ´ l ( N / 2 ) | or | u u ´ l ( N / 2 ) | increasing (See Appendix A for the derivation).
Property 3.
It can also be proved that if the Euclidean distance between any two of { ( α ´ l ( N 2 ) , u ´ l ( N 2 ) ) } l = 1 L is sufficiently large, n = 1 N | H n ( α , u ) | 2 = n = 1 N | l = 1 L S l , n ( α , u ) | 2 has L distinguishable peaks corresponding to the L multipath components in the ( α , u ) plane and the peaks’ coordinates correspond to { ( α ´ l ( N 2 ) , u ´ l ( N 2 ) ) } l = 1 L (See Appendix B for the derivation).

3.2. DFRFT Based Multipath TOA Estimation

Based on the above properties of multipath components in (8), the parameters of multipath components can be estimated through searching for the peaks of CFRs in the FRFT domain. However, the estimation error caused by interferences and/or AWGN cannot be neglected in cases of insufficiently large virtual antenna size M × ρ or at low SNRs.
Due to the Property 3 above, an averaging method over N subcarriers can be used to reduce the estimation error caused by noise. Furthermore, a successive interference cancellation (SIC) method can be employed to mitigate the interference effects as in [19]. Thus, the parameters of multipath components will be obtained by combining the SIC method and the averaging method.
Without loss of generality, we assume that the amplitudes { | b l | } l = 1 L of L paths in the channel satisfy | b 1 | > | b 2 | | b L | . Considering the SIC and averaging-based multipath parameters estimation is composed of L similar iterative steps, only the steps in the l -th iteration are exemplified as follows. Assume the initial filtered CFRs H ˜ m , n ( 0 ) = H ^ m , n .

3.2.1. Searching for the Peak

Given the filtered CFRs H ˜ m , n ( l 1 ) in the ( l 1 ) -th step, the peak of n = 1 N | H ˜ n ( l 1 ) ( α , u ) | 2 in the ( α , u ) plane at the l -th iteration can be expressed as
( α ^ ( l ) , u ^ ( l ) ) = arg max ( α , u ) n = 1 N | H ˜ n ( l 1 ) ( α , u ) | 2 = arg max ( α , u ) n = 1 N | F α ( { H ˜ m , n ( l 1 ) } m = M 2 M 2 1 ) | 2
The estimated peak corresponds to the strongest multipath component in the filtered CFRs H ˜ m , n ( l 1 ) , i.e., the l -th strongest multipath component in the CFRs { H ^ m , n } m = M 2 M 2 1 . According to the Property 3 above, it approximates the peak of the DFRFT of the l -th multipath component at the N / 2 -th subcarrier in the ( α , u ) plane.
( α ^ ( l ) , u ^ ( l ) ) ( α ´ l ( N 2 ) , u ´ l ( N 2 ) )
Since the DFRFT transformed signal for each subcarrier means averaging the noise over M positions on the trajectory and the summation averages the noise energy over N subcarriers, n = 1 N | H ˜ n ( l 1 ) ( α , u ) | 2 contains the averaged noise energy across both the frequency and space dimensions, thus the noise effects on the estimation ( α ^ ( l ) , u ^ ( l ) ) can be reduced greatly considering that the subcarrier number and the number of the virtual antenna array elements are usually large.
From (5) and (9), the peak of the DFRFT of the l -th multipath component at the n-th subcarrier in the ( α , u ) plane can be derived as
α ´ ^ l ( n ) = a r c c o t ( f c + n Δ f f c + N Δ f / 2 c o t ( α ^ ( l ) ) ) u ´ ^ l ( n ) = f c + n Δ f f c + N Δ f / 2 c s c ( α ^ ( l ) ) c s c ( α ´ ^ l ( n ) ) u ^ ( l )

3.2.2. Estimating the TOA of the Strongest Multipath

From the relationship between the peaks of the DFRFT of the CFRs and the parameters of L multipath-chirp components in (7), the peak of the DFRFT at the n-th subcarrier is compensated as
H ˜ ´ n ( l 1 ) ( α ´ ^ l ( n ) , u ´ ^ l ( n ) ) = H ˜ n ( l 1 ) ( α ´ ^ l ( n ) , u ´ ^ l ( n ) ) A α ´ ^ l ( n ) * e j π ( u ´ ^ l ( n ) ) 2 c o t α ´ ^ l ( n ) | A α ´ ^ l ( n ) | 2 M b l e j φ l ( n ) + W ˜ ´ n ( α ´ ^ l ( n ) , u ´ ^ l ( n ) ) + { l = l + 1 L S l , n ( α ´ ^ l ( n ) , u ´ ^ l ( n ) ) } · A α ´ ^ l ( n ) * e j π ( u ´ ^ l ( n ) ) 2 c o t α ´ ^ l ( n ) I n
where the approximation holds since the first l 1 strongest multipath component are almost removed in the previous l 1 iterations, the first term in the right hand represents the compensated DFRFT of the strongest multipath component, W ˜ ´ n ( α ´ ^ l ( n ) , u ´ ^ l ( n ) ) and I n are the compensated and filtered AWGN and interferences caused by the remaining multipath components, respectively.
By substituting (5) into (14), the compensated DFRFT at the peak ( α ´ ^ l ( n ) , u ´ ^ l ( n ) ) in (14) can be further expressed as
H ˜ ´ n ( l 1 ) ( α ´ ^ l ( n ) , u ´ ^ l ( n ) ) d l e j 2 π n Δ f τ 0 , l + W ˜ ´ n ( α ´ ^ l ( n ) , u ´ ^ l ( n ) ) + I n
with
d l = | A α ´ ^ l ( n ) | 2 M b l e j ω c τ 0 , l
From (15), H ˜ ´ n ( l 1 ) ( α ´ ^ l ( n ) , u ´ ^ l ( n ) ) can be viewed as a new channel frequency response at the n-th subcarrier. Due to the energy concentration property of DFRFT and the uniformly distributed noise in the DFRFT domain, the strongest multipath component dominates the compensated DFRFT of the n-th subcarrier at the peak ( α ´ ^ l ( n ) , u ´ ^ l ( n ) ) . Hence, this new channel has a dominated multipath of interest, which has the same TOA τ 0 , l as the l multipath.
Therefore, the existing TOA estimation methods, such as the IFT method or the MP-based method [9], can then be directly used on { H ˜ ´ n ( l 1 ) ( α ´ ^ l ( n ) , u ´ ^ l ( n ) ) } n = 1 N to obtain the estimate of the strongest multipath’s TOA in this iteration, which is denoted as τ ^ 0 , l . Moreover, the estimate of the strongest multipath’s TOA is irrelevant to the virtual antenna array element spacing ρ from (16). It is worth mentioning, although the MP-based TOA estimation is susceptible to the noise, it can remove most interference in such new equivalent channel, especially at high SNRs.

3.2.3. Cancelling the Strongest Component

Furthermore, the filtered CFRs H ˜ m , n ( l 1 ) at the n-th subcarrier in the ( l 1 ) -th step will be transformed by DFRFT of order α ´ ^ l ( n ) as
H ˜ n ( l 1 ) ( α ´ ^ l ( n ) , u ) = F α ´ ^ l ( n ) ( { H ˜ m , n ( l 1 ) } m = 1 M ) n = 0 , 1 , , N
For each u in a neighborhood of u ´ ^ l ( N / 2 ) , [ H ˜ 0 ( l 1 ) ( α ´ ^ l ( 0 ) , u ) , , H ˜ N 1 ( l 1 ) ( α ´ ^ l ( N 1 ) , u ) ] is formed as a new equivalent CFR vector over N subcarriers. Similarly, it can be derived from (7) that the l -th multipath component is dominant in such formed CFR vector. Again, the IFT or MP-based multipath parameters estimation method [9] can be employed on such a CFR vector to obtain the TOA and the complex gain of its strongest multipath component, which are denoted as τ ˜ 0 , l ( u ) and b ˜ 0 , l ( u ) . Consequently, the DFRFT at u can be filtered by cancelling this strongest multipath component from such a CFR vector as
H ˜ n ( l ) ( α ´ ^ l ( n ) , u ) = H ˜ n ( l 1 ) ( α ´ ^ l ( n ) , u ) b ˜ 0 , l ( u ) e j 2 π n Δ f τ ˜ 0 , l ( u ) n = 0 , 1 , , N
if the TOA estimate difference | τ ˜ 0 , l ( u ) τ ^ 0 , l | is smaller than a chosen threshold.
Then the filtered DFRFT H ˜ n ( t ) ( α ´ ^ l ( n ) , u ) is transformed back to the frequency domain by DFRFT of order α ´ ^ l ( n ) . Hence, the filtered CFRs in the l -th iteration can be expressed as
H ˜ m , n ( l ) = F α ´ ^ l ( n ) ( H ˜ n ( t ) ( α ´ ^ l ( n ) , u ) ) n = 0 , 1 , , N
The above procedures will repeat until the amplitude of the detected multipath component is smaller than a certain threshold. Finally, the TOA of the direct multipath can be estimated as the minimum one among { τ ^ 0 , l } l = 1 L .

4. Simulation Results and Discussions

4.1. Simulation Settings

To evaluate the performance of the proposed algorithm, we have conducted computer simulations in a relatively simple 9 m × 7 m × 4 m (length × width × height) three-dimension indoor WLAN environment, in which there are 1 AP located at ( 0.5 , 4.5 , 2.5 ) and four mixed cement and glass walls. The 2-dimension layout of the indoor WLAN environment is shown in Figure 2. The classical ray-tracing propagation model [20] is adopted in the simulations to generate the CIRs at indoor positions, and the dielectric constant and electrical conductivity of concrete and glass are respectively ϵ r 1 = 8 , σ 1 = 0.01 and ϵ r 2 = 2 , σ 2 = 0.001 , and the system measurement noise is Gaussian. For simplicity, only the indirect paths created by a single specular reflection from a side wall and the direct path are considered in the simulations. The carrier frequency, subcarrier number, and subcarrier spacing are assumed as 5.0 GHz, 64, and 625 KHz, respectively.
Two kinds of scenarios have been studied in the simulations: LOS scenarios and NLOS scenarios. In the LOS scenarios, all the multipaths’ gain are generated according to the ray-tracing propagation model, thus the LOS multipath is the strongest one among all the multipaths. In the NLOS scenarios, all multipaths except the LOS multipath have the same gain as in LOS scenarios, while the LOS path is additionally attenuated by 15 dB to simulate the presence of human body or some barrier between the AP and MT according to [21,22].
To average the LOS path TOA estimation performance over space, the target positions, i.e., the central positions of MT’s moving trajectories, are set as uniformly spaced positions with the distance of 1.5 m between neighboring central positions in this indoor environment, thus there are 20 uniformly spaced target positions also shown in Figure 2. Based on the CFRs at the positions equi-spaced on a line trajectory, the direct path’s TOA at the central position of the corresponding line trajectory is obtained as the target TOA estimate in the simulations. The mean square error (MSE) of the LOS path’s TOA estimate (unit: ns) is employed as the performance measurement and all results are obtained by averaging 400 independent simulations: 20 independent simulations at each one of 20 uniformly spaced target positions.
In the simulations, we compare the performance of the proposed algorithms, the FRFT-IFT algorithm and the FRFT-MP algorithm, with the 2D MP-based super-resolution TOA estimation algorithm [10]. In each iteration of the proposed algorithms, the DFRFT is implemented on the CFRs across the space domain at each subcarrier according to (6) with { α = p π / 2 , p = q / 800 } q = 800 800 , { u = m / M } m = M 2 M 2 1 and M = 60 , and the parameters of the strongest multipath components is extracted by searching the peak in the discretized FRFT domain according to (11), then the FRFT-IFT algorithm and the FRFT-MP algorithm respectively employ the IFT scheme and MP scheme [9] to estimate the strongest multipath’s TOA, the neighborhood of the peak in the u domain is set as the whole u domain and the TOA estimate difference threshold is set as 1 n s in the step of cancelling the strongest component. The iteration repeats until the energy of the extracted component is 20 times less than that in the first iteration.

4.2. Simulation Results

It is worth noting that the performance of the three algorithms are mainly dependent on the noise, the multipath components interference and the mismatch errors of the multipaths’ phase model. In other words, the performance of the LOS path TOA estimation depends on the ‘SINR’ of the LOS path component, i.e., the ratio of the LOS path energy to the sum of the energy of noise, the multipath components interference and the mismatch errors. Since the 2D MP method almost removes the multipath interference in both frequency and space domain, the performance of the 2D MP TOA estimation algorithm is mainly relevant to the noise and the multipaths’ phase model mismatch errors. Due to adopting the quadratic function as the multipath phase variation model across the space dimension, the two proposed algorithms: the FRFT-IFT algorithm and the FRFT-MP algorithm, suffer from much less multipaths’ phase model mismatch errors. Besides the mismatch errors and the noise, the FRFT-IFT algorithm is also subject to multipath components interference caused by the insufficient virtual antenna array size, while the FRFT-MP method can remove most of the multipath interference by employing MP-based TOA estimation method [9] in frequency domain, thus it is only subject to the residual multipath components interference. Moreover, the ‘SNR’ in the following simulations means the ratio of the sum of all the multipaths’ energy to the noise energy, whereas the ‘SNR of the LOS path component’ means the LOS path energy to the noise energy.
Figure 3 demonstrates the MSE of the TOA estimation versus SNR in LOS scenarios and NLOS scenarios, respectively. In the simulations, the virtual antenna array size are assumed as M ρ = 60 × 2 cm. From Figure 3, it is obvious that the proposed algorithms: the FRFT-MP algorithm and the FRFT-IFT algorithm generally outperform the 2D MP algorithm since the proposed algorithms can make more of the multipaths’ characteristics in the space domain. Figure 3a shows that the two proposed algorithms perform better than the 2D MP algorithm in NLOS scenarios, especially at low SNRs, because in such scenarios the noise dominates the MSE performance at low SNRs, and the proposed methods concentrate the LOS multipath energy by FRFT, thus can achieve higher ‘SINR’ of the LOS path component than the 2D MP method, while the 2D MP method suffers from noise threshold effects [23].
In LOS scenarios and high SNR cases in NLOS scenarios, the 2D MP algorithm and the FRFT-IFT algorithm are mainly subject to the mismatch error of the quasi-linear phase model and the multipath components’ interference, respectively, since the SNR of the LOS path component is sufficiently high in such cases. This explains that both these two methods have error floor effects at high SNRs in NLOS scenarios as shown in Figure 3a and their performance improves slightly with SNR increasing in LOS scenarios as shown in Figure 3b. The FRFT-MP method can mitigate the multipath interference by employing MP-based TOA estimation method [9], thus the FRFT-MP method has the lower error floor as shown in Figure 3a. On the other hand, since the MP scheme is susceptible to the noise, the FRFT-MP method performs worse than the FRFT-IFT method at low SNRs in both NLOS and LOS scenarios, as shown in Figure 3a,b.
Figure 4 demonstrates the MSE of the TOA estimation versus different sizes of virtual antenna array at SNR = 25 dB. Please note that the SNR of the LOS path component is not very high due to the relatively weak LOS path component in NLOS scenarios. As the virtual antenna array size increases, more multipaths’ characteristics in the space domain can be exploited to enhance the ‘SINR’ of the LOS path component until the multipaths’ phase model mismatch errors are not negligible at cases of relatively big virtual antenna array size in NLOS scenarios. This explains that the performance of the three algorithms improve with the virtual antenna array size increasing from 60 × 0.5 cm to 60 × 2.5 cm, but they are worsened at the case of virtual antenna array size 60 × 3 cm in NLOS scenarios, which are shown in Figure 4a. Moreover, Figure 4a shows that the two proposed methods outperform the 2D MP method, because the proposed methods concentrate the LOS path energy and suffer from less model mismatch errors, thus achieve higher ‘SINR’ of the LOS path component.
At high SNRs in LOS scenarios, the SNR of the LOS path component is much higher than that in NLOS scenarios due to the dominant LOS path, hence all the three methods perform much better in LOS scenarios than in NLOS scenarios by comparing Figure 4b with Figure 4a. It can be seen from Figure 4b that the performance of the 2D MP algorithm and the FRFT-MP method only improves with virtual antenna array size increasing from 60 × 0.5 cm to 60 × 1.0 cm and then it is worsened with virtual antenna array size increasing from 60 × 1.0 cm to 60 × 3.0 cm, and the FRFT-IFT method performs better with virtual antenna array size increasing from 60 × 0.5 cm to 60 × 2.0 cm and performs worse with virtual antenna array size increasing from 60 × 2.0 cm to 60 × 3.0 cm, because the multipaths’ phase model mismatch errors dominate the performance of the three methods in most virtual antenna array sizes considering the ‘SINR’ of the LOS path component is considerably high in LOS scenarios. Figure 4b also shows that the performance of the proposed algorithms are generally inferior to that of the 2D MP algorithm, especially at cases of relative small virtual sizes, since the performance of the 2D MP algorithm in LOS scenarios is dominated by only the LOS path’s phase model mismatch error, which is greatly less than the multipaths’ phase model mismatch errors in NLOS scenarios. Besides the LOS path’s phase model mismatch error, the proposed FRFT-IFT method and the FRFT-MP method are also affected by the multipath interference and the residual multipath interference, respectively.
Considering that the virtual antenna array size is imprecisely determined and is deliberately chosen as big as allowed to make more of the multipath characteristics, and whether the indoor environment is NLOS or LOS scenarios is unknown, the proposed algorithms, especially the FRFT-MP algorithm, can achieve more robust performance than the 2D MP algorithm by taking the above simulation results into account.

5. Conclusions

In this paper, a novel virtual antenna array and fractional Fourier transform-based 2-dimension super-resolution TOA estimation algorithm was presented without requiring a real antenna array. By employing the virtual antenna array idea and modeling the multipath phase variation as a quadratic function, the proposed algorithm can make more use of the multipaths’ characteristics in the space dimension, thus it can enhance the multipath resolution and achieve more robust TOA estimation performance, especially in NLOS scenarios. Simulation results confirm the effectiveness of the proposed algorithm.

Author Contributions

Z.C. proposed the TOA estimation algorithm, designed the simulation scheme and wrote the paper. L.W. gave suggestions on the proposed algorithm and helped to improve the quality of the final manuscript. M.Z. helped to perform the simulation tests and analyzed the results.

Funding

The work was supported in part by the National Natural Science Foundation of China under Grant No. 61871317 and the Natural Science Basic Research Plan in Shaanxi Province of China under Grant No. 2015JM6305.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
FRFTfractional Fourier transform
DFRFTdiscrete fractional Fourier transform
TOAtime-of-arrival
CFRchannel frequency responses
IFTinverse Fourier transform
MPmatrix pencil
MUSICmultiple signal classification
LOSLine-of-Sight
NLOSNon-Line-of-Sight
APAccess Point
MTmobile terminal
CIRchannel impulse response
AWGNadditive white Gaussian noise

Appendix A. Proof of the Property 2 of S l , n (α,u)

The part of the DFRFT signal in (7) contributed by the lth multipath component is given by
n = 1 N | S l , n ( α , u ) | 2 = n = 1 N { S l , n ( α , u ) } { S l , n ( α , u ) } * = n = 1 N { | c l | 2 m = 1 M k = 1 M { e j 2 π [ u c s c ( α ) / M + f l ( n ) ρ ] ( m k ) · e j π [ c o t ( α ) / M + μ l ( n ) ρ 2 ] ( m 2 k 2 ) } } = | c l | 2 m = 1 M k = 1 M { e j 2 π u c s c ( α ) ( m k ) M e j π c o t ( α ) ( m 2 k 2 ) M · n = 1 N { e j 2 π f l ( n ) ρ ( m k ) e j π μ l ( n ) ρ 2 ( m 2 k 2 ) } }
From (5), it can be observed that both f l ( n ) and μ l ( n ) are linear to n. Hence, we have
n = 1 N { e j 2 π f l ( n ) ρ ( m k ) e j π μ l ( n ) ρ 2 ( m 2 k 2 ) } d l , m , k · e j 2 π f l ( N / 2 ) ρ ( m k ) e j π μ l ( N / 2 ) ρ 2 ( m 2 k 2 )
with
d l , m , k = n = N 2 + 1 N 2 e ( j 2 π n Δ f ( c o s θ l ρ ( m k ) c + ( s i n θ l ) 2 ρ 2 ( m 2 k 2 ) 2 τ 1 , l c 2 ) ) = n = N 2 + 1 N 2 c o s ( 2 π n Δ f ( c o s θ l ρ ( m k ) c + ( s i n θ l ) 2 ρ 2 ( m 2 k 2 ) 2 τ 1 , l c 2 ) ) N c o s ( π N Δ f ( c o s θ l ρ ( m k ) c + ( s i n θ l ) 2 ρ 2 ( m 2 k 2 ) 2 τ 1 , l c 2 ) ) N
where the first approximation is derived by using trigonometric function series summation, the second approximation holds since M ρ is deliberately selected to satisfy M ρ τ 1 , l c and M ρ 2 c N Δ f generally holds in narrowband positioning systems.
By substituting (A2) into (A1), we have
n = 1 N | S l , n ( α , u ) | 2 | c l | 2 m = 1 M k = 1 M { d l , m , k e j 2 π [ u c s c ( α ) M + f l ( N 2 ) ρ ] ( m k ) · e j π [ c o t ( α ) M + μ l ( N 2 ) ρ 2 ] ( m 2 k 2 ) }
Hence, the coordinate ( α ^ , u ^ ) of the peak of n = 1 N | S l , n ( α , u ) | 2 is given by
α ^ = a r c c o t ( μ l ( N / 2 ) M ρ 2 ) = α ´ l ( N 2 ) u ^ = f l ( N / 2 ) M ρ / c s c ( α ^ l ( N / 2 ) ) = u ´ l ( N 2 )
from (A4).
By substituting (A2) into (A4), n = 1 N | S l , n ( α , u ) | 2 can be further approximated as
n = 1 N | S l , n ( α , u ) | 2 N | c l | 2 · m = 1 M e j π { [ 2 u c s c ( α ) M + f l ( N 2 ) ρ ] m e [ c o t ( α ) M + μ l ( N 2 ) ρ 2 ] m 2 } 2 = N | c l | 2 m = 1 M e j π { [ 2 u ´ l ( N 2 ) c s c ( α ´ l ( N 2 ) ) 2 u c s c ( α ) M ) ρ ] m e [ c o t ( α ) c o t ( α ´ l ( N 2 ) ) M ρ 2 ] m 2 } 2
According to the DFRFT property, it is not difficult to derive from (A6) that n = 1 N | S l , n ( α , u ) | 2 decreases rapidly with | α α ^ l ( N / 2 ) | or | u u ^ l ( N / 2 ) | increasing.
Therefore, the Property 2 of S l , n ( α , u ) is proved.

Appendix B. Proof of the Property 3 of S l , n (α,u)

The summation of DFRFT of CFRs over N subcarriers can be expressed as
n = 1 N | H n ( α , u ) | 2 = n = 1 N | l = 1 L S l , n ( α , u ) | 2 = n = 1 N { l = 1 L S l , n ( α , u ) l = 1 L { S l , n ( α , u ) } * } = n = 1 N { l = 1 L | S l , n ( α , u ) | 2 + l = 1 L l l L S l , n ( α , u ) { S l , n ( α , u ) } * } = n = 1 N l = 1 L | S l , n ( α , u ) | 2 + n = 1 N { l = 1 L l l L S l , n ( α , u ) { S l , n ( α , u ) } * } l = 1 L n = 1 N | S l , n ( α , u ) | 2
where the approximation holds because the cross term between any two of { S l , n ( α , u ) } l = 1 L is negligible due to the large Euclidean distance between any two of { ( α ´ l ( N 2 ) , u ´ l ( N 2 ) ) } l = 1 L . Furthermore, l = 1 L n = 1 N | S l , n ( α , u ) | 2 has L distinguishable peaks corresponding to the peaks of { n = 1 N | S l , n ( α , u ) | 2 } l = 1 L , i.e., { ( α ´ l ( N 2 ) , u ´ l ( N 2 ) ) } l = 1 L .
Therefore, the Property 3 is proved.

References

  1. Zou, H.; Huang, B.; Lu, X.; Jiang, H.; Xie, L. A Robust Indoor Positioning System based on the Procrustes Analysis and Weighted Extreme Learning Machine. IEEE Trans. Wirel. Commun. 2016, 15, 1252–1266. [Google Scholar] [CrossRef]
  2. Peral-Rosado, J.A.; Saloranta, J.; Destino, G.; López-Salcedo, J.A.; Seco-Granados, G. Methodology for Simulating 5G and GNSS High-Accuracy Positioning. Sensors 2018, 18, 3220. [Google Scholar] [CrossRef] [PubMed]
  3. Schrooyen, F.; Baert, I.; Truijen, S.; Pieters, L.; Denis, T.; Williame, K.; Weyn, M. Real time location system over Wi-Fi in a healthcare environment. J. Inf. Technol. Healthc. 2006, 4, 401–416. [Google Scholar]
  4. Li, X.; Pahlavan, K. Super-Resolution TOA Estimation with Diversity for Indoor Geolocation. IEEE Trans. Wirel. Commun. 2004, 3, 224–234. [Google Scholar] [CrossRef]
  5. Wen, C.; Chan, F. Adaptive AOA-Aided TOA Self-Positioning for Mobile Wireless Sensor Networks. Sensors 2010, 10, 9742–9770. [Google Scholar] [CrossRef] [PubMed]
  6. Cheon, J.; Hwang, H.; Kim, D.; Jung, Y. IEEE 802.15.4 ZigBee-Based Time-of-Arrival Estimation for Wireless Sensor Networks. Sensors 2016, 16, 203. [Google Scholar] [CrossRef] [PubMed]
  7. Saarnisaari, H. TLS-ESPRIT in a time delay estimation. In Proceedings of the IEEE 47th Vehicular Technology Conference, Phoenix, AZ, USA, 4–7 May 1997; Volume 3, pp. 1619–1623. [Google Scholar]
  8. Ding, R.; Qian, Z.; Jiang, H. TOA estimation for IR-UWB system using Matrix Pencil. In Proceedings of the 2009 World Congress on Computer Science and Information Engineering, Los Angeles, CA, USA, 31 March–2 April 2009; Volume 1, pp. 461–464. [Google Scholar]
  9. Vanderveen, M.; Van der Veen, A.J.; Paulraj, A. Estimation of multipath parameters in wireless communications. IEEE Trans. Signal Process. 1998, 46, 682–690. [Google Scholar] [CrossRef]
  10. Gaber, A.; Omar, A. A Study of Wireless Indoor Positioning Based on Joint TDOA and DOA Estimation Using 2-D Matrix Pencil Algorithms and IEEE 802.11ac. IEEE Trans. Wirel. Commun. 2015, 14, 2440–2454. [Google Scholar] [CrossRef]
  11. Liu, L.; Liu, H. Joint Estimation of DOA and TDOA of Multiple Reflections in Mobile Communications. IEEE Access 2016, 4, 3815–3823. [Google Scholar] [CrossRef]
  12. Han, S.; Li, Y.; Meng, W.; He, C. A New High Precise Indoor Localization Approach Using Single Access Point. In Proceedings of the 2017 IEEE Global Communications Conference, Singapore, 4–8 December 2017; pp. 1–5. [Google Scholar]
  13. Almeida, L.B. The Fractional Fourier Transform and Time-Frequency Representations. IEEE Trans. Signal Process. 1994, 42, 3084–3095. [Google Scholar] [CrossRef]
  14. Tao, R.; Li, B.Z.; Wang, Y. Spectral analysis and reconstruction for periodic nonuniformly sampled signals in fractional Fourier domain. IEEE Trans. Signal Process. 2007, 55, 3541–3547. [Google Scholar] [CrossRef]
  15. Feng, S.; Murray-Smith, R. Fusing Kinect Sensor and Inertial Sensors with Multi-rate Kalman Filter. In Proceedings of the IET Conference on Data Fusion Target Tracking: Algorithms and Applications, Liverpool, UK, 30 April 2014; Volume 10, pp. 1–8. [Google Scholar]
  16. Kang, W.; Han, Y. SmartPDR: Smartphone-Based Pedestrian Dead Reckoning for Indoor Localization. IEEE Sens. J. 2015, 15, 2906–2916. [Google Scholar] [CrossRef]
  17. Kuang, J.; Niu, X.; Chen, X. Robust Pedestrian Dead Reckoning Based on MEMS-IMU for Smartphones. Sensors 2018, 18, 1391. [Google Scholar] [CrossRef] [PubMed]
  18. Ozaktas, H.M.; Arikan, O.; Kutay, M.A.; Bozdagt, G. Digital computation of the fractional Fourier transform. IEEE Trans. Signal Process. 1996, 41, 2141–2150. [Google Scholar] [CrossRef]
  19. Qi, L.; Tao, R.; Zhou, S.; Wang, Y. Detection and Parameter Estimation of Multicomponent LFM Signal Based on the Fractional Fourier Transform. Sci. China (Ser. F) 2004, 47, 184–198. [Google Scholar] [CrossRef]
  20. Chen, S.; Jeng, S. An SBR/Image Approach for Radio Wave Propagation in Indoor Environments with Metallic Furniture. IEEE Trans. Antennas Propag. 1997, 45, 98–106. [Google Scholar] [CrossRef]
  21. Villanese, F.; Evans, N.E.; Scanlon, W.G. Pedestrian-Induced Fading for Indoor Channe-l s at 2.45,5.7 and 62 GHz. In Proceedings of the 2000 IEEE Vehicular Technology Conference Fall, Boston, MA, USA, 24–28 September 2000; Volume 1, pp. 43–48. [Google Scholar]
  22. Seesai, W.; Chamchoy, M.; Promwong, S. A Body-Shadowing Model for Indoor UWB Communication Environments. In Proceedings of the 2008 5th International Conference on Electrical Engineering/Electronics, Computer, Telecommunications and Information Technology, Krabi, Thailand, 14–17 May 2008; Volume 1, pp. 261–264. [Google Scholar]
  23. Hua, Y.; Sarkar, T.K. On SVD for estimating Generalized Eigenvalues of Singular Matrix Pencil in Noise. IEEE Trans. Signal Process. 1991, 39, 892–900. [Google Scholar] [CrossRef]
Figure 1. The TOAs of the l-th path at equi-spaced positions on a line moving trajectory.
Figure 1. The TOAs of the l-th path at equi-spaced positions on a line moving trajectory.
Sensors 19 00638 g001
Figure 2. The 2-dimension layout of the indoor WLAN environment.
Figure 2. The 2-dimension layout of the indoor WLAN environment.
Sensors 19 00638 g002
Figure 3. Mean square error versus SNR in NLOS scenarios and LOS scenarios, respectively.
Figure 3. Mean square error versus SNR in NLOS scenarios and LOS scenarios, respectively.
Sensors 19 00638 g003
Figure 4. Mean square error versus virtual antenna array size in NLOS scenarios and LOS scenarios, respectively.
Figure 4. Mean square error versus virtual antenna array size in NLOS scenarios and LOS scenarios, respectively.
Sensors 19 00638 g004

Share and Cite

MDPI and ACS Style

Chen, Z.; Wang, L.; Zhang, M. Virtual Antenna Array and Fractional Fourier Transform-Based TOA Estimation for Wireless Positioning. Sensors 2019, 19, 638. https://doi.org/10.3390/s19030638

AMA Style

Chen Z, Wang L, Zhang M. Virtual Antenna Array and Fractional Fourier Transform-Based TOA Estimation for Wireless Positioning. Sensors. 2019; 19(3):638. https://doi.org/10.3390/s19030638

Chicago/Turabian Style

Chen, Zhigang, Lei Wang, and Mengya Zhang. 2019. "Virtual Antenna Array and Fractional Fourier Transform-Based TOA Estimation for Wireless Positioning" Sensors 19, no. 3: 638. https://doi.org/10.3390/s19030638

APA Style

Chen, Z., Wang, L., & Zhang, M. (2019). Virtual Antenna Array and Fractional Fourier Transform-Based TOA Estimation for Wireless Positioning. Sensors, 19(3), 638. https://doi.org/10.3390/s19030638

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

Article Metrics

Back to TopTop