Joint TOA and DOA Estimation for UWB Systems with Antenna Array via Doubled Frequency Sample Points and Extended Number of Clusters

(e ultra-wideband (UWB) system, which transmits information using nanosecond or even sub-nanosecond pulses, has been widely applied in precise positioning. In this paper, we investigate the problem of the time of arrival (TOA) estimation and the direction of arrival (DOA) estimation in the UWB systems with antenna array and propose a joint TOA and DOA estimation algorithm with doubled frequency sample points and extended number of clusters. Specifically, the proposed algorithm uses two antennas to receive impinging signals and utilizes the conjugate symmetry characteristic of the delaymatrices to extend the sample points as well as the number of clusters. Moreover, in order to obtain TOA estimates with low computational complexity, the proposed algorithm transforms the two-dimensional (2D) spectral search to one-dimensional (1D) searches. (e DOA estimates can then be achieved by using the TOA estimation results and the geometric information. Simulation results are given to testify the performance of the proposed algorithm.


Introduction
e ultra-wideband (UWB) technique is a kind of wireless communication technology which uses nanosecond or even sub-nanosecond pulses as carrier to transmit information. Benefiting from the high transmission rate, low power consumption, and the anti-multipath characteristics, the UWB system has been widely applied in various fields such as radar, imaging, and positioning [1][2][3]. Due to the extremely narrow pulse, the positioning precision of the UWB system can reach centimeter level or even millimeter level. Even under complex multipath conditions such as indoor environment, the UWB system can still achieve accurate positioning due to its strong anti-multipath and penetration ability [4][5][6]. At present, IEEE 802. 15.4a standard has taken UWB as the preferred technology for positioning application [7].
One of the basic problems in the UWB positioning system is the time of arrival (TOA) estimation. TOA estimation algorithms can be divided into two categories: one is the traditional algorithms based on the time domain, and the other is the high-resolution algorithms based on the frequency domain. e former mainly includes the coherent detection method using pulse template matched filter [8] and the incoherent TOA estimation algorithm based on threshold or energy detection [9][10][11]. e coherent algorithm based on matched filter can obtain TOA estimates with high accuracy but at the same time with high sampling rate, complex receiver structure, and expensive equipment cost. e incoherent TOA estimation algorithm has the advantages of low sampling rate, fast convergence speed, and low hardware resource occupation rate. However, the low sampling rate leads to low time resolution, which reduces the accuracy of the TOA estimation. e traditional time-domain-based algorithms obtain TOA estimates by estimating the arrival time of the direct path (DP) component in the received signal. However, due to the multipath effect and the nonline-of-sight (NLOS) condition, DP may not be the strongest path; thus the algorithm resolution declines. erefore, super-resolution estimation algorithms based on frequency-domain processing are proposed [12][13][14][15][16][17][18]. ese algorithms, such as the multiple signal classification (MUSIC) algorithm [12,13], the propagator method (PM) [14], the estimation of signal parameters via rotational invariance techniques (ESPRIT) algorithm [15,16], and the matrix pencil algorithm [17,18], model the channel impulse response in frequency domain and realize the TOA estimation using the orthogonality between the signal subspace and the noise subspace, which can achieve high estimation resolution.
Researches on DOA estimation in UWB systems were also done. Since the UWB system has high time resolution, some joint TOA and DOA estimation algorithms were proposed. e matrix pencil algorithm was applied to the joint estimation in [19]. In [20], the rough TOA estimates were obtained via energy estimation and the minimum distance criterion, the accurate estimates were then achieved by a low-complexity and high-resolution method based on the signal power delay spectrum, and the DOA estimates were finally obtained by the minimum-variance unbiased estimation using the TOA estimation results. Besides, the electromagnetic vector sensor (EMVS) with multiple-input multiple-output (MIMO) radar [21] also has the potential to be exploited in the UWB systems.
In this paper, we investigate the problem of joint TOA and DOA estimation in UWB systems and propose a computationally efficient algorithm. Specifically, the proposed algorithm uses two antennas to receive impinging signals and utilizes the conjugate symmetry characteristic of the delay matrices to double the equivalent frequency sample points and extend the number of clusters. Moreover, to obtain TOA estimates with low computational complexity, the proposed algorithm transforms the two-dimensional (2D) spectral search to one-dimensional (1D) searches. e DOA estimates can then be achieved by using the TOA estimation results and the geometric information. e main contributions of this paper are summarized as follows: (1) We propose an effective algorithm for UWB systems to jointly obtain the DOA and TOA estimation results. (2) We transform the time-consuming 2D spectral search to twice 1D searches in order to release the computational burden. (3) e proposed algorithm can obtain high estimation accuracy due to the doubled sample points and the extended number of clusters. e rest of this paper is organized as follows: Section 2 introduces the data model of the UWB system. e proposed joint TOA and DOA estimation algorithm is derived in Section 3. Section 4 then analyses the estimation performance of the proposed algorithm and Section 5 gives the simulation results to testify the algorithm effectiveness. Finally, Section 6 concludes this paper. Notation 1. In this paper, we use upper (lower) bold characters to represent matrices (vectors). (·) * , (·) T , and (·) H , respectively, denote the operation of conjugate, transpose, and conjugate transpose. * denotes the convolution operation, and ⊗ denotes the Kronecker product. A(:, m: n) generates a new matrix consisted by the m-th to the n-th columns of the matrix A. E · { } represents the expectation operation. I M stands for an M × M identity matrix.

Data Model
2.1. Transmit Signals of UWB Systems. In this paper, we utilize the second derivation of the Gaussian pulse as the UWB transmit signal and modulate it with the direct sequence binary phase shift keying (DS-BPSK). e transmit signal of the UWB system can then be expressed as [22] where b c ∈ −1, +1 { } is the modulated binary data symbol sequence, c n ∈ −1, +1 { } is the pseudorandom sequence used to realize the multiple access communication, T c denotes the repeat period of the pulse, T s represents the period of the binary data symbol, and N c is the pulse repetition times of a single binary data symbol. e term p(t) is the second derivation of the Gaussian pulse and can be further expressed as where Γ is the pulse forming factor related to the pulse width.

Channel Model of UWB Systems.
According to the Saleh-Valenzuela (SV) model [23], we consider that the transmit signal produces multiple multipath components after passing through the channel, and these multipath components arrive at the receiver in the form of clusters. Specifically, we assume there are K clusters and L multipath per cluster in the UWB channel, then the channel impulse response of the kth cluster is given by where α (k) l is the channel fading factor of the lth path in the kth cluster, which obeys the Rayleigh distribution. θ (k) l distributes uniformly in the range [0, 2π], δ(·) is the Dirac function, and τ (k) l is the channel delay of the lth path in the kth cluster. Normally, the change rate of the channel is slow compared with the pulse speed of the transmit signal; hence we have τ (k) to denote the random complex fading amplitude, then (3) can be rewritten as

Received Signals of UWB Systems.
According to the basic theory of the digital signal processing, the received signal of the kth cluster in the time domain can be expressed as [23] 2 Mathematical Problems in Engineering where w (k) (t) is the additive Gaussian white noise of the kth cluster. Transform (5) into frequency domain, i.e., where and w (k) (t). Subsequently, by sampling N (N > L) points with interval Δω � (2π/N) in the frequency domain, we can obtain the measurement data vector as [23] where vector is the delay matrix which contains the information of the signal multipath delay, where e( L ] T contains the complex channel fading factor of the kth cluster and In this paper, we use two array antennas to receive the impinging signals. As shown in Figure 1, we assume there are L far-field signals from direction θ 1 , θ 2 , . . . , θ L . e TOA of the two antennas is τ � [τ 1 , τ 2 , . . . , τ L ] and ς � [ς 1 , ς 2 , . . . , ς L ], respectively. e distance between the two antennas is d, and c denotes the light speed. According to (7), the received signal of the two antennas in the frequency domain can be expressed as [23] K . E 1 (τ) and E 2 (ς) are the delay matrices of the two antennas, respectively, which can be expressed as e estimation of the two channel impulse responses in the frequency domain can be achieved by where V 1 � (W 1 /S) and V 2 � (W 2 /S). In the next section, we propose a joint TOA and DOA estimation algorithm utilizing the channel impulse response H 1 and H 2 .

Sample Points and Clusters Extension.
e cross-correlation matrix of the channel impulse response can be constructed as where R B � E BB H ∈ C L×L is a diagonal matrix. We divide the cross-correlation matrix R H ∈ C N×N into two N × (N − 1) matrices, i.e., where X 1 and X 2 consist of the first and the last N − 1 columns of R H and E 2a (ς) and E 2b (ς) consist of the first and the last N − 1 rows of E 2 (ς). According to the form of the delay matrix, we have where Ψ(ς) is an L × L diagonal matrix given by According to (9) and (10)

Mathematical Problems in Engineering
where J N ∈ C N×N is the antidiagonal matrix, and the rotation matrices Φ(τ) and Φ(ς) are, respectively, given by By utilizing the above properties of the delay matrices, we can construct a new matrix using X 1 and X 2 as where Equation (19) can be seen as the equivalent channel impulse response with doubled sample points and increased number of clusters, which can improve the maximum number of detectable signals as well as the estimation performance.

(26)
According to (26), we can construct the cost function as where ρ is a constant value. In order to get the extremum, we can construct the partial derivation of L(τ, ς) with respect to q(ς), i.e., us, we have q(ς) � μF −1 (τ)u with μ � −0.5ρ. Considering u H q(ς) � 1, the constant μ can be further expressed as erefore, the vector q(ς) can be further transformed into Substituting (30) into (26), then the TOA estimation result of the first antenna is given by which means we can get the TOA estimation of the first antenna by a 1D spectral search with the spectral function Similarly, in order to obtain the TOA estimates of the second antenna, we can exchange the order of H 1 and H 2 when constructing the cross-correlation matrix, i.e., en, by following the same procedures described above, the 1D spectral function with respect to the TOA of the second antenna is given by After obtaining the TOA estimates of both antennas, we can estimate the DOA of the impinging signals by using the TOA estimates as well as the geometric information, i.e., e main steps of the proposed algorithm are summarized as follows: (1) Construct the cross-correlation matrix R H according to (12). (2) Divide R H into X 1 and X 2 according to (13) and (14), and construct the extended observation matrix X as (19). (3) Compute the correlation matrix of X, and perform the eigenvalue decomposition to obtain the noise subspace U v . (4) Construct the spectral function according to (32) and perform the 1D spectral search to obtain the TOA estimates τ.

Performance Analysis
4.1. Complexity. According to Section 3, when estimating the TOA of the first antenna τ, constructing R H needs O N 2 K , computing R X needs O 8N 2 (N − 1) , and performing eigenvalue decomposition on it requires O 8N 3 . Constructing the spectral function and conducting 1D spectral search needs O 4N 2 (2N − L) + 8n s (N 2 + N + 1) , where n s denotes the search times. erefore, the total complexity of the proposed algorithm is O 2N 2 K + 48N 3 − 8N 2 (L + 2) + 16n s (N 2 + N + 1) . If we perform 2D spectral search on equation (22), then the complexity is increased to e computational complexity of the proposed algorithm and some other algorithms are concluded in Table 1. Besides, the complexity comparison versus the number of frequency sample points of different algorithms is depicted in Figure 2, where K � 100, L � 3, and n s � 100. As shown in Figure 2, the complexity of the proposed algorithm is remarkably reduced compared with the algorithms with 2D spectral search and is close to the complexity of the 2D-PM.

Algorithm Advantages.
We conclude the advantages of the proposed algorithm as follows: (1) e proposed algorithm reduces the computational complexity by transforming the 2D spectral search to twice 1D spectral search. (2) e proposed algorithm doubles the equivalent frequency sample points and increases the equivalent number of clusters, which can improve the maximum number of detectable signals as well as the estimation accuracy. (3) e proposed algorithm can obtain higher estimation accuracy compared with some existed algorithms.

Simulation Results
In this section, we perform large number of Monte Carlo trials to examine the effectiveness as well as the superiority of the proposed algorithm. e root mean square error (RMSE) of the trials is used to measure the estimation accuracy, which is given by where Q is the number of Monte Carlo trials and τ l,q , ς l,q , and θ l,q are the estimation result of τ l , ς l , and θ l in the qth trial, respectively. In the following simulations, we set Q � 500, d � 0.09 m, and L � 3 with τ � [0.2 ns, 0.3 ns, 0.4 ns]. Figures 3 and 4 show the TOA estimation results of the proposed algorithm, where each point represents a trial, K � 100, N � 64, and the signal-to-noise ratio (SNR) is −10 dB and 10 dB, respectively. As shown in the two figures, the proposed algorithm can successfully obtain the paired TOA estimation results and the estimation accuracy is improved with higher SNR. Figures 5 and 6, respectively, compare the TOA and DOA estimation accuracy with different number of clusters, where N � 64. e simulation results show that the estimation accuracy of both TOA and DOA improves with increased number of clusters. Figures 7 and 8, respectively, depict the TOA and the DOA estimation performance comparison versus SNR of different algorithms, where K � 100 and N � 64. It is illustrated clearly that the proposed algorithm can achieve Mathematical Problems in Engineering

Conclusion
In this paper, we investigate the problem of joint TOA and DOA estimation in UWB systems and propose a computationally efficient algorithm with doubled frequency sample points and extended number of clusters. e proposed algorithm utilizes the conjugate symmetry characteristic of the delay matrices to extend the sample points as well as the number of clusters and then transforms the 2D spectral search into twice 1D search in order to reduce the computational complexity. Simulations testify that the proposed algorithm can obtain high estimation accuracy and large number of identifiable signals with low complexity.

Data Availability
e data used to support the findings of this study are included within the article.