Efficient Adaptive Algorithms for DOA Estimation in Wireless Communications ()
The problem of direction-of-arrival (DOA) estimation in mobile communication systems requires efficient algorithms with a high spatial resolution and a low computational complexity when moving sources are considered. MUSIC is known as a high resolution algorithm for DOA estimation but, this method demands high computational complexity to compute the signal subspace of the time-varying data of the correlation matrix. This paper focuses on MUSIC using subspace tracking methods, such as Bi-SVD and PAST, to carry out iterative DOA estimation. Accuracy and the processing time of both methods are evaluated and compared with the results of MUSIC. These results show the potential of Bi-SVD and PAST to reduce the processing time and to improve the accuracy when the number of snapshots and the source angular variation increases, assessing the source location for a dynamic mobile cellular environment.
1. Introduction
In mobile communication systems source localization with distributed sensor arrays can be performed by estimating the direction of arrival (DOA) of the signals coming from mobile terminals (sources). DOA is one of the most challenging problems which one has to solve for localizing and tracking multiple moving source as in radar, mobile communications and in other areas.
Within the class of high-resolution source localization methods, the MUSIC (MUltiple SIgnal Classification) method [1] has received the most attention and has been widely studied [2–4]. However, when the problem of moving sources is addressed, this method demands heavy computational load due to the decomposition of the subspace of the vectors of the correlation matrix, which is estimated based on an N-size sample [4–7]. For fast DOA tracking, the subspace tracking methods iteratively estimate the signal subspace which in this case is involved in the MUSIC spectrum. To address this problem, Bi-SVD (Bi-Iteration Singular-Value Decomposition) [4,5] and PAST (Projection Approximation Subspace Tracking) [6] have been proposed and investigated due to their capabilities of successively updating (tracking) eigenvectors in the signal subspace of a correlation matrix [7–9]. Compared with the rectangular window in MUSIC conventional, in this signal subspace algorithms an exponential window is used for the signal processing. This window considers a forgetting factor to reduce the effect of the past observations which are down-weighted.
In this work, we present and evaluate a scheme that incorporates those subspace tracking methods into MUSIC, and shows a comparison of DOA successive estimation performance. For a small sample size in [4] is found that all the methods show the same root mean square error (RMSE) after convergence, and that all the methods can satisfactorily estimate DOAs. In contrast, we use a larger sample size to estimate the DOA considering the source’s mobility. Compared with MUSIC conventional, our experimental analysis on simulated data demonstrates the accuracy improvement of both Bi-SVD and PAST algorithms to assess the sources DOA for a dynamic environment in terms of the root mean square error (RMSE).
2. System Model
Let x(t) be the Kx1 data vector at the output of an uniform linear array of K elements spaced a distance d between two consecutive elements. Assuming that the array of K sensors receive D narrowband signal waves from far-field sources with the same known center frequency, then the output vector x(t) can be described as:
(1)
where the Dx1 vector s(t)=[s0(t), …, sD-1(t)]T and the Kx1 vector n(t)=[n1(t), ..., nK(t)]T denote the complex amplitudes of the signals and the measurement noise at time t. The KxD steering matrix A is composed by D steering vectors (Kx1), with al(ql)=[1, exp(-j(2π/l)d senql), …, exp(-j(2π/l)d(K-1)senql)], ql is the direction of arrival and l is the carrier wavelength. The noise samples are assumed to be zero mean with variance s2.
From the eigen decomposition of the correlation matrix R=E[x(t)xH(t)] are obtained Es(t) and En(t) corresponding to the signal and noise subspaces respectively, [1]. The signal subspace is formed by the D eigenvectors corresponding to the D largest eigenvalues. Usually, the MUSIC spectrum is computed with the noise subspace to improve the spectral resolution. Besides that, the subspace tracking algorithms approaches the signal subspace Es(t), then, the identity En(t)En(t)H=I-Es(t)Es(t)H must be employed to track the DOA from the MUSIC spectrum as follows [4]:
(2)
In what follows, are explained Bi-SVD and PAST algorithms for the Es(t) signal subspace tracking.
2.1. Bi-SVD Algorithm
Let x(t) be a random complex process with the correlation matrix R=E[x(t)x(t)H] and the signal subspace Es(t) contains the D dominant eigenvectors. For subspace tracking, the data matrix X(t) can be updated in the time according [5]:
(3)
where x(t) is the current data vector, 0<β<1 is the forgetting factor and X(t-1) is the data matrix at time instant t-1. Consider the following bi-iteration applied on the (NxK) data matrix X(t) as given in [5]:
, (4)
where QA(0) is the initial value for QA(t). Here, QA(t) denotes an estimate of the signal subspace Es(t), RA(t) and RB(t) are D-dimensional upper-triangular matrices obtained from QR decomposition. Notice that the data matrix X(t) in (3) is updated recursively then it results in a growing matrix. A solution of this problem is to approximate X(t) as [5]:
(5)
With the suboptimal approximation in (4), a fast exponential window based Bi-SVD subspace tracking algorithm was developed in [5]. Now, projecting x(t) into the space spanned by QA(t-1) the complement of its orthogonal projection is
, (6)
where. The normalization of x^(t) results in
, (7)
with. From (5), the decomposition for x(t) becomes:
(8)
Using (5) and (8) into (3) and then substituting X(t) in (4), we obtain QA(t) after several computations as follows:
(9)
where
(10)
(11)
In Equation 9, the subspaces RA(t), RB(t) and QA(t) are updated successively with initial values RA(0)=0, RB(0)=ID and QA(0)=ID.
2.2. PAST Algorithm
The PAST algorithm is based on the idea that QA(t) is the signal subspace when is minimized (12).
(12)
An exponentially weighted sum is considered into J(QA(t)), where β is the forgetting factor. At time t, all sample vectors are involved in the estimation of the signal subspace. The algorithm PAST approximates QAH(t)x(i) in (12) with y(i)=QAH(i-1)x(i) by using the unknown projection of x(i) into the columns of QA(t) [6]. The approximated cost function is then given by:
. (13)
The choice of QA(t) that minimize J'(QA(t)) is [6]
(14)
(15)
(16)
Therefore, the minimization process yields that J(QA(t)) on only internal noise, resulting in QA(t) equal to the signal subspace matrix Es(t), obtained recursively from (14).
3. Simulation Results
Successive update algorithms using Bi-SVD and PAST into MUSIC are compared in estimation accuracy with the conventional MUSIC. The estimation accuracy is evaluated by RMSE (Root Mean Square Error), defined as follows:
(17)
wherelp defines the estimated DOAs, D denotes the
Figure 1. RMSE evaluation in DOA estimation (dynamic environment). SNR=20 dB spatial variation d=0.01°.
Figure 2. RMSE evaluation in DOA estimation (dynamic environment). SNR=20 dB spatial variation d=0.02°.
Figure 3. Normalized processing time. Signal subspace estimation considering variation of the number of array antenna elements.
number of sources and P the number of independent trials.
For the simulations, we consider a uniform linear array (ULA) of K=10 sensors separated a half wavelength of the actual source signals. We assume the case of two narrowband signals in the far-field, located at [-30°, 50°] and separated from the antenna array by 930m and 425m respectively, transmitting with the same power and the same carrier frequency of f=1900 MHz. For this case, the number of independent trials is set to P=1000, assuming a forgetting factor of β=0.9.
Simulation results are obtained by incorporating over time the sources’ spatial angular variation d as the dynamic component. Figure 1 compares the RMSE obtained by the PAST, Bi-SVD and MUSIC methods in a dynamic environment. For each snapshot, the angular variation d is increased a factor of 0.01 degrees, that is, the sources location at time t are given by [-30°+td, 50°+ td]. The results in Figure 1 show a small RMSE for PAST and Bi-SVD methods. After a 40 snapshots the accuracy of the MUSIC method is degraded with an accumulated RMSE value of 0.9° approximately. On the other hand, for Bi-SVD and PAST, their RMSE error increases slowly until a value close to 0.26°. In this case, comparing the subspace tracking algorithms into MUSIC with the conventional, they improve the accuracy as the number of snapshots increases.
In a second case, is considered an increase of angular variation (d=0.02°). From Figure 2, we found that after 20 snapshots MUSIC conventional rapidly increases its RMSE error from 0.25° to 2° on the maximum number of snapshots considered. In this case, the accuracy of the MUSIC method is remarkably worse compared with the maximum RMSE error of 0.3° obtained by the subspace tracking algorithms. The forgetting factor used by subspace tracking algorithms improved the accuracy into MUSIC to estimate the DOA, while MUSIC conventional fails as the number of snapshots and angular variation increases.
Finally, the processing time of the methods is evaluated and normalized taking into account the maximum processing time. In Figure 3 the processing time is normalized with the maximum time obtained during the simulations. We observe that the processing time of the PAST algorithm is smaller than both Bi-SVD and MUSIC. More interestingly, the results show that PAST and Bi-SVD processing times are independent of the number of antenna array elements, situation that does not hold by MUSIC which grows linearly with the number of antenna elements.
4. Conclusions
The PAST and Bi-SVD subspace-tracking algorithms were applied to MUSIC in order to face the accuracy and processing time problems in DOA estimation. The simulations show good results as far as the application of the subspace-tracking algorithms. Our experimental analysis on simulated data demonstrates the potential reduction in processing time and accuracy improvement to assess the sources DOA when the number of snapshots and angular variation increases. It is shown that as source mobility increases, it produces a higher error for MUSIC conventional. However, the PAST and Bi-SVD algorithms show a lower estimation error, varying quietly with the mobility of sources. On the other hand, the processing time of PAST algorithm is smaller than both Bi-SVD and MUSIC conventional. The processing time required by these subspace tracking algorithms is independent of the number of antenna array elements.
5. Acknowledgments
This work was supported by PROMEP and the Mexican National Science and Technology Council, CONACyT, under Grant 52374.
[1] R. O. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions Antennas and Propagation, Vol. 34, pp. 276–280, 1986
[2] D. H. Covarrubias and J. G. Arceo, “Improving resolution on DOA estimation in a multipath macrocell environment,” Proc. IEEE SympoTIC, Bratislava Slovak Republic, pp. 39–42, 2004.
[3] J. G. Arceo, D. H. Covarrubias, and J. M. Luna, “Efficiency evaluation of the unconditional maximum likelihood estimator for near-field source DOA estimation,” ETRI Journal, Vol. 28, pp. 761–769, 2006.
[4] N. Kikuma, “Iterative DOA estimation using subspace tracking methods and adaptative beamforming,” IEICE Transactions on Communication, Vol. E88-B, pp. 1818– 1828, 2005.
[5] P. Strobach, “Bi-Iteration SVD subspace tracking algorithms,” IEEE Trans. Signal Processing, Vol. 45, pp. 1222– 1240, 1997.
[6] B. Yang, “Projection approximation subspace tracking,” IEEE Transactions Signal Proceedings, Vol. 43, pp. 95– 107, 1995.
[7] P. Strobach, “Bi-Iteration multiple invariance subspace tracking and adaptive ESPRIT,” IEEE Transactions on Signal Processing, Vol. 48, pp. 442–456, 2000.
[8] A. Kuchar, M. Tangemann, and E. Bonek, “A real-time DOA-based smart antenna processor,” IEEE Transactions on Vehicular Technology, Vol. 51, pp. 1279–1293, 2002.
[9] S. Ouyang and Y. Hua, “Bi-Iterative least-square method for subspace tracking,” IEEE Transactions on Signal Processing, Vol. 53, pp. 2984–2996, 2005.