International Journal of Antennas and Propagation

Volume 2012 (2012), Article ID 803617, 6 pages

doi:10.1155/2012/803617

Research Article

Angle of Arrival Estimation Using Cholesky Decomposition

Saleh O. Al-Jazzar

Electrical Engineering Department, University of Ha'il, Ha'il 81451, Saudi Arabia

Received 23 November 2011; Revised 22 February 2012; Accepted 5 March 2012

Academic Editor: Stefano Selleri

Copyright © 2012 Saleh O. Al-Jazzar. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

An angle of arrival (AOA) estimator is presented. Many applications require accurate AOA estimates such as wireless positioning and signal enhancement using space-processing techniques. The proposed AOA estimator depends on the Cholesky decomposition of the received signal autocorrelation matrix. The resultant decomposed matrices are used to modify the crosscorrelation matrix of the received signals at the antenna array doublets. The proposed method is named the Cholesky-decomposition-based-AOA (CDBA) estimator. In comparison with the TLS-ESPRIT algorithm which utilizes the eigenvalue decomposition (EVD) of the received signal autocorrelation and crosscorrelation matrices, the CDBA method has better performance than the TLS-ESPRIT algorithm especially in low signal-to-noise-ratio (SNR) cases. Simulations for the proposed CDBA method are shown to assess its performance.

1. Introduction

Accurate estimation of the received signal angle of arrival (AOA) can be very beneficial for signal reception enhancement [1, 2] and wireless positioning [3, 4]. If the AOA estimates were not estimated accurately, then wireless devices would not be located accurately or the received signal bit error rate would be high. Thus, AOA estimators are developed and modified in such cases to increase their accuracy.

Many AOA estimators were proposed in the literature [5–9]. The authors in [5] present an iterative method to estimate the AOA. In [6], the authors propose estimating the AOA using the frequency domain of the received signals which requires taking their Fourier transform. The authors in [7] introduce an iterative AOA estimator which depends on the generalized expectation maximization algorithm. In [8], the authors present an AOA estimator for circular arrays. Another very popular AOA estimator is the TLS-ESPRIT algorithm proposed in [9] which utilizes subspace techniques to estimate the AOA.

In this paper, an AOA estimator is proposed which enhances the AOA estimation. The proposed method is named the Cholesky-decomposition-based-AOA (CDBA) estimator. The CDBA method utilizes the received signal at the two sides of the antenna array doublets. An autocorrelation matrix of the received signal at one side of the antenna array doublets is calculated. Then, the received signal autocorrelation matrix is decomposed using the Cholesky decomposition. A crosscorrelation matrix between the received signals at both sides of the antenna doublets is calculated. The decomposed matrices of the autocorrelation matrix together with the crosscorrelation matrix are used to form a new matrix from which the AOAs are estimated.

The CDBA method does not require taking the Fourier transform of the received signals nor it requires any iterative or searching procedure to estimate the AOA. Also, in comparison with the well-known TLS-ESPRIT algorithm of [9], the accuracy of the CDBA method is higher than that of the TLS-ESPRIT method especially in low signal-to-noise-ratio (SNR) cases.

The paper is organized as follows: Section 2 introduces the system model that forms the foundation for the proposed estimator. The proposed CDBA estimator is presented in Section 3. Section 4 presents the simulated performance of the CDBA estimator. Finally, conclusions are shown in Section 5.

2. System Model

This section presents the narrowband received signal model that will be utilized for the AOA estimation. The antenna array is formed from uniform antenna doublets (i.e., total antenna elements). Each antenna doublet is formed of two antenna elements spaced by a distance . We assume a BPSK sources signals, , impinging upon the antenna array. The signal is represented as , where is the carrier frequency , is an integer which represents the time index and is the complex low-pass equivalent of , is the symbol period, , where is discrete which represents the symbol parity and is a positive constant which represents the amplitude of , and is the (raised cosine) pulse-shaping function with

where is an integer so that .

The AOA is taken between the antenna array axis, and the source signal arrival direction and will be given the notation . In this paper, the superscript notation is used to denote the estimated value of a variable (for example is the estimated value of ). So, to start developing the received signal model, let us first consider the sampled received signal at the antenna elements (after the matched filter stage), from to , that are located in the first side of the antenna doublets, and let us call it , and is represented as

where

where is the signal wavelength, and is the distance between consecutive antenna elements. In addition, is the noise vector added to the received signal at the antenna elements (from ) located along the first side of the antenna doublets, which is additive white Gaussian noise (AWGN) and has a covariance matrix of , where is the identity matrix. Finally, represents the transpose operation.

The received signal vector at the first side of the antenna doublets set can be written in matrix form as

where

From now on the term will be dropped from all terms for simplicity.

The received signal vector at the second side of the antenna doublets set will be given the notation . Thus,

where is an AWGN vector at the second side of the antenna array doublets, and

with

Next section will describe the proposed CDBA estimator.

3. Proposed Cholesky-Decomposition-Based AOA (CDBA) Estimator

To implement the CDBA method, the crosscorrelation matrix () between and is formulated as follows:

where

where is the power of the th source.

Another matrix considered by the CDBA algorithm is the autocorrelation matrix of , as follows:

Take the noise free version of , and let us call it , where , that is,

Then, the eigenvalue decomposition (EVD) of is

where is the eigenvector matrix of , and is an diagonal matrix with its diagonal that contains the eigenvalues of , that is,

where the nonzero values of the diagonal elements of correspond to the sources which will be called . Each element of corresponds to one of the sources.

Comparing (12) with (13), then we have

Comparing both sides of (15), then we have

where is the appropriate matrix to change the basis.

Substituting (16) into (15), then we have

The matrix is an matrix, and is a diagonal matrix. Also, the nonzero diagonal elements of are only the first diagonal elements. So, (17) can be written as

Looking at (18), it can be deduced that the matrix should be . Also, the lower part of should be all zero elements, that is,

where

But, the left hand side of (18) is a diagonal matrix, and the matrix is diagonal as well. So, from (18) it can deduced that the matrix is diagonal too, that is,

Thus, from (18), (19), and (21), each of the nonzero diagonal elements of is given by

which has a positive real value.

Consequently, the EVD of would be given by

where is an diagonal matrix and is given by , that is,

Substituting (16) into (23), then we have

where the notation is the pseudoinverse for a given matrix.

Now, taking the Cholesky decomposition of , then we have

where is a lower triangular matrix that results from the Cholesky decomposition.

Comparing (26) with (25), it can be deduced that

where is an orthogonal rotating matrix.

A new matrix is formed by multiplying the left and right hand sides of by and , respectively, as follows:

But, as shown in (9) and as shown in (27). Thus,

Now, let

Also, let . Then, (29) can be written as

Looking at (31), if it is proved that the matrix is a diagonal matrix, then (31) will be the EVD of , where and are the corresponding eigenvectors and eigenvalues matrices, respectively.

To prove that is a diagonal matrix, consider the matrix which is a diagonal matrix of size and the matrix which is an matrix where its upper part is diagonal, and its lower part contains zero elements, then

where is a diagonal matrix of size where its diagonal elements are the same as the first diagonal elements of , that is,

Thus, (30) can be written as

Since , and are diagonal matrices, then we have

But, the matrix is diagonal matrix where its elements are positive real. So, , and

Substituting for , , , and in (36), then we have

where

But, is defined in (22) as . Thus, (38) can be written as

Thus, it is clear that is a diagonal matrix and it is proved now that (31) is the EVD of .

Defining the term

then (39) can be written as

Now, because has a positive real value, then

From the definition of , it is clear that the AOAs are implicated in . Also, from (42) the phase of equals the phase of , and for are the diagonal elements of . Thus, to find , consider the th diagonal element of as follows:

where is the th diagonal element of the estimated .

Thus, the required AOAs are estimated.

3.1. Summary of the CDBA Method

(1)First step: calculate from (9).(2)Second step: calculate from (11).(3)Third step: calculate by taking the Cholesky decomposition of .(4)Fourth step: calculate the matrix by .(5)Fifth step: calculate the eigenvalues of matrix.(6)Sixth step: calculate the AOAs from the largest eigenvalues of matrix as shown in (43).

3.2. Analysis for the Effect of Noise Variance

The performance of the AOA estimators in high noise power (low SNR) cases is of interest. A measure of the noise power is the noise variance. A study of the noise variance effect on the performance of the CDBA algorithm is derived and compared with that of the TLS-ESPRIT.

As shown in (43), the AOAs are estimated from the eigenvalues (). Thus, to start let us consider the effect of noise variance (power) on .

Consider the magnitude of as follows:

But has a positive real value (recall that , and and have positive real values). Thus, . But, , so it is clear that . So, the magnitude of the eigenvalue () is given by

The noise variance effect on the can be measured by defining the variable which is defined as the absolute difference between and its noise-free version , where . Thus, is defined as follows:

Clearly, the minimum value for is when . Whereas its maximum value is as .

As for the TLS-ESPRIT algorithm, the AOA estimation starts by taking the EVD of an autocorrelation matrix of the received signals (call it ) [9]. Let be defined as the noise-free version of . Also, let for be one of the matrix eigenvalues. Then, each eigenvalue of (which is given the notation ) will be given as follows:

Defining as the absolute difference between and its noise free version (), then,

Thus the minimum value for is when . Whereas its maximum value approaches as .

Comparing and in (46) and (48), respectively, it is clear that in low noise variance cases, the deviation between the estimated eigenvalues and the true eigenvalues will be small for both algorithms (CDBA and TLS-ESPRIT); thus, it is expected that both algorithms will have close performance in low noise variance cases (i.e., high SNR).

In high noise variance cases, the deviation between the estimated eigenvalues and the true eigenvalues will have a maximum value of 1 in CDBA algorithm, whereas, in the TLS-ESPRIT the deviation will approach . Recall that the CDBA algorithm estimates the AOAs from the eigenvalues () for , and the TLS-ESPRIT estimates the AOAs from the eigenvectors of (see [9] for more insight), and any error in estimating the eigenvalues of will directly change their corresponding eigenvectors causing an error in estimating the AOAs. Also, it is clear that in high noise variance cases (i.e., low SNR) the EVD process in the CDBA algorithm is much less affected by the noise than the EVD process in the TLS-ESPRIT algorithm. Thus, it is expected that the CDBA algorithm will perform better than the TLS-ESPRIT in high noise variance cases (i.e., low SNR).

The performance of the proposed CDBA method will be shown in the simulation section (Section 4).

4. Simulation Results

Simulations for the proposed CDBA estimator were completed to assess its performance. The results in this section were averaged over 1000 ensemble runs. The elements of each antenna array were separated by a half-wavelength (i.e., ), and without any loss of generality the distance between the elements of each antenna doublet () was chosen to be as well. The number of sources was set to 2. was set to 1.2, and was set to 1. The proposed CDBA method was compared with the TLS-ESPRIT method of [9].

Figure 1 shows the root-mean-square error (RMSE) of the AOA estimation in degrees for the proposed CDBA method compared with the TLS-ESPRIT method for different SNRs. Since the power for both users was not equal, then the SNR was taken for the second user; that is, if the SNR was set to 5 dB, then the second user would have a 5 dB SNR and the first user would have SNR of dB. The angles for the two users were set to and . The number of antenna doublets was set to . The number of snap shots was set to and (where is the number of snap shots over which the correlation matrices were estimated). The results in Figure 1 show that the proposed CDBA method gave better performance than the TLS-ESPRIT method in low SNR cases. This result agrees with the discussion presented in Section 3.2.

Figure 1: Root mean square error (RMSE) of angle estimation in degrees versus signal-to-noise-ratio (SNR) in dB for .

Figure 2 shows the RMSE of the AOA estimation in degrees for the proposed CDBA method compared with the TLS-ESPRIT method for different number of antenna doublets. The angles for the two users were set to and . The number of snap shots was set to . The results in Figure 2 show that the CDBA method outperformed the TLS-ESPRIT algorithm for different number of antenna doublets.

Figure 2: Root mean square error (RMSE) of angle estimation in degrees versus number of antenna doublets for .

Another important issue in comparing AOA estimators is their performance when the AOAs of the received signals are close to each other. Table 1 shows the RMSE of the AOA estimation in degrees for the proposed CDBA method compared with the TLS-ESPRIT method for different AOA separations for two users. The number of antenna doublets was set to . The number of snap shots was set to . The angular deviation was performed by setting to be taken from the following equation: where is the angular deviation with . Table 1 shows that the CDBA method outperformed the TLS-ESPRIT algorithm for different AOA separations.

Table 1: Root mean square error (RMSE) of angle estimation in degrees for different AOA separations.

AOA separation
TLS-ESPRIT
CDBA

1°
11.5223°
3.1770°

1.5°
6.8057°
1.6336°

2°
4.2739°
0.7330°

Also, to compare both methods computational complexity, two functions in MATLAB (tic.m and toc.m) were used to measure the time it takes each method to estimate the AOA. The CDBA and the TLS-ESPRIT methods took 0.281 and 0.297 seconds, respectively, to perform one run of AOA estimation for two users. Clearly, this indicates that the CDBA method has less computational complexity than the TLS-ESPRIT method.

5. Conclusion

In this paper, an AOA estimator is proposed which is named the CDBA method. The CDBA method is applied by taking the Cholesky decomposition of the received signal autocorrelation matrix. The resultant decomposed matrices are used to modify the crosscorrelation matrix of the received signals at the antenna array doublets. The proposed CDBA method has better performance than the TLS-ESPRIT method in estimating the AOAs especially in low SNR cases. The performance of the proposed CDBA method was assessed and compared to the TLS-ESPRIT method.

References

N. Hew and N. Zein, “Space-time estimation techniques for UTRA system,” in Proceedings of the IEE Colloquium on Capacity and Range Enhancement Techniques for the 3rd Generation Mobile Communications and Beyond (Ref. No. 2000/003), pp. 6/1–6/7, London, UK, 2000.

Y. F. Chen and M. D. Zoltowski, “Joint angle and delay estimation for DS-CDMA with application to reduced dimension space-time rake receivers,” in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '99), vol. 5, pp. 2933–2936, March 1999.

S. Al-Jazzar, M. Ghogho, and D. McLernon, “A joint TOA/AOA constrained minimization method for locating wireless devices in non-line-of-sight environment,” IEEE Transactions on Vehicular Technology, vol. 58, no. 1, pp. 468–472, 2009.

P. Deng and P. Fan, “An AOA assisted TOA positioning system,” in Proceedings of the 2000 International Conference on Communication Technology Proceedings (WCC-ICCT '00), vol. 2, pp. 1501–1504, Beijing, China, August 2000.

Y. H. Choi, “Maximum likelihood estimation for angles of arrival of coherent signals using a coherency profile,” IEEE Transactions on Signal Processing, vol. 48, no. 9, pp. 2679–2682, 2000.

X. Huang and Y. J. Guo, “Frequency-domain AoA estimation and beamforming with wideband hybrid arrays,” IEEE Transactions on Wireless Communications, vol. 10, no. 8, pp. 2543–2553, 2011.

J. Verhaevert, E. V. Lil, and A. V. de Capelle, “Direction of arrival (DOA) parameter estimation with the SAGE algorithm,” Signal Processing, vol. 84, no. 3, pp. 619–629, 2004.

Q. Cheng, R. Yang, and H. Zhang, “Optimally weighted ESPRIT using uniform circular arrays,” Computers and Electrical Engineering, vol. 31, no. 4-5, pp. 272–281, 2005.

R. Roy and T. Kailath, “ESPRIT-estimation of signal parameters via rotational invariance techniques,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, no. 7, pp. 984–995, 1989.