ANTENNA ANGLE-OF-ARRIVAL ESTIMATION USING UPLINK WEIGHT VECTORS

Abstract: The angle ( theta A, theta B) at which a multipath finger (A, B) produced by a digital radio communication source arrives at an antenna array can be estimated using an uplink weight vector (w) associated with the multipath finger (A, B). Due to a characteristic of uplink weight vectors (w) used in antenna arrays having a shift invariance structure, the estimate of the angle ( theta A, theta B) can be calculated based on the uplink weight vector (w).

Description: ANTENNA ANGLE-OF-ARRIVAL ESTIMATION USING UPLINK WEIGHT VECTORS FIELD OF THE INVENTION The invention relates generally to the use of adaptive antenna arrays in digital radio communications and, more particularly, to angle-of-arrival (AOA) estimation for downlink use in such adaptive antenna arrays. BACKGROUND OF THE INVENTION In array signal processing of digital communications, an array of receive and transmit antennas is deployed in the base station (e.g. in a cellular telephone network) to combat interference and in that way increase capacity and performance of the system. There are two main concepts for exploiting the potential performance gain from an antenna array. One approach is to use the antenna to form a number of fixed beams in the uplink. The receiver algorithm then uses the output from one or more beams to receive a desired user. In the downlink, the data is then transmitted in the beam with the highest signal-to-noise ratio (SNR), or in several beams in the case of transmit diversity systems. Another possibility is to steer a narrow beam directly towards a desired mobile. This can be done in two fundamentally different ways. In a parametric approach, the angles-of-arrival of the multipaths (fingers) of the desired user are explicitly estimated, and the uplink beams are then steered in the determined direction. The AOA giving the highest SNR is then used for the downlink transmission. Again, in the case of transmit diversity, several downlink beams can be utilized. In a non-parametric approach, the AOAs of the desired user are not estimated explicitly, but only the spatial channels (or spatial signatures) are estimated, giving rise to an estimated array weight vector (obtained for example using a vector RAKE receiver) for each finger. The combination of these weights applied to the array output yields the signal used to receive the user data. The non-parametric methods do not yield any explicit information on the AOAs for the user, and these AOAs have to be estimated separately. There is much scientific literature on AOA estimation using adaptive antennas. The interest in AOA estimation was mainly spurred by the advent of the so-called super-resolution algorithms, such as MUSIC, ESPRIT, MODE and WSF, methods that can estimate the AOAs with higher accuracy than the classical beamforming algorithms. One of the major shortcomings of the popular super- resolution methods is their computational complexity. They often involve performing the singular value decomposition of a matrix the size of the array. Another drawback of these methods is their reported sensitivity to model errors; their ability to give good estimates of the AOAs relies heavily on a parametric model of the signal environment and an accurate calibration of the antenna array. The non-parametric methods, such as beamforming, basically work by scanning the area covered by the array with a narrow beam and selecting the AOA as the direction giving the highest output power from the array. As such, these methods are more robust, but their bad resolution can be a limiting factor in practice, where many users are present close to each other and where noise levels might be high. It is desirable in view of the foregoing to provide for AOA estimation without the aforementioned disadvantages of the prior art approaches. In an attempt to avoid the aforementioned disadvantages of the prior art approaches, the present invention exploits a characteristic of conventionally available uplink weight vectors used in antenna arrays having a shift invariance structure, and calculates AOA estimates directly from the uplink weight vectors. BRIEF DESCRIPTION OF THE DRAWINGS FIGURE 1 diagrammatically illustrates the angular relationship between an antenna normal vector of an antenna sensor and two multipath fingers resulting from a single digital radio communication source. FIGURE 2 diagrammatically illustrates how a multipath finger can be modeled as a large number of arrays impinging from different angles. FIGURE 3 graphically illustrates the angular dependency of antenna gain. FIGURE 4 diagrammatically illustrates the propagation delay between sensors of a uniform and linear antenna array. FIGURE 5 diagrammatically illustrates exemplary embodiments of an angle of arrival estimator according to the invention. FIGURE 6 diagrammatically illustrates an exemplary digital radio communications interface according to the invention. FIGURE 7 illustrates exemplary operations which can be performed by the angle of arrival estimator embodiments of FIGURE 5. FIGURES 8 and 9 graphically illustrate simulation results associated with a least-squares angle of arrival estimation according to the invention. FIGURES 10 and 11 graphically illustrate simulation results associated with a recursive least-squares angle of arrival estimation according to the invention. DETAILED DESCRIPTION The present invention draws upon the fact that, for some antenna array structures, the components of the estimated weight vector for one multipath finger are approximately (neglecting the effect of noise and errors in the estimation of the weight vector) related by a simple phase shift. This can be seen from the model described below. FIGURE 1 diagrammatically illustrates the angular relationship of two exemplary multipath fingers A and B to the antenna normal vector C of a sensor in an antenna array. The multipath fingers A and B are the well-known result of various physical obstructions (buildings, hills, etc.) which re-direct a radio signal produced by a radio source, for example a digital cellular telephone. The angle between multipath finger A and the antenna normal vector C is designated in FIGURE 1 as Θ A , and the angle between multipath finger B and the antenna normal vector C is designated as Θ B . Each multipath finger can be assumed to have some angular spread around its nominal direction of propagation. For example, FIGURE 2 illustrates that multipath finger B can be modeled as a large number of rays, four of which are illustrated in FIGURE 2. The illustrated rays R j -R 4 impinge upon the multipath finger B from respective angles θ, - θ 4 . Thus, each of a large number L of rays (four of which are shown in FIGURE 2) impinges from a separate angle θ,, for 1 = 1 to L. Each of these rays also has a complex amplitude designated herein as h,(t,θ,). Each of the angles θ, can be modeled as a random variable with mean Θ B and variance σ, 2 . Using these assumptions, the weight vector for a given multipath finger can be written as follows: Equation 1 > K )= Σ /=ι M ) β(θ ) + ( 0 The variable θ on the left side of Equation 1 represents the angle of arrival of the particular multipath finger. Thus, for multipath finger B of FIGURES 1 and 2, the variable θ on the left side of Equation 1 corresponds to Θ B of FIGURES 1 and 2. In Equation 1, h,(t,θ|) denotes the base band channel for ray 1 at time t. As is well known in the art, this base band channel includes the effect of propagation attenuation, receiver filtering and pulse shaping. Also, the antenna gain, which depends on the angle θ, (see FIGURE 3) is included in the base band channel h,(t,θ,). The factor a(θ,) in Equation 1 accounts for the effect of the time delay as the wave propagates through the array. As shown in FIGURE 4, the propagation time delay corresponds to the time necessary for the wave corresponding to ray R, to travel the distance Δsinθ,. The term E(t) in Equation 1 is due to the fact that the weight vector is estimated in a noisy environment, using a finite number of training symbols. For narrow-band signals, for example digital radio communication signals, a time delay corresponds to a phase shift. In Equation 1 , a(θ,) is a vector of these phase shifts called an array response vector. In the example of FIGURE 4, the antenna array is a conventional uniform and linear array (ULA) including M sensors 1, 2, 3 ... M. In this example, the M-dimensional array response vector a(θι) is given by Equation 2 1 ;Δsιnθ, α(θ,) = 2/Δsιnθ, (M-l)ιΔsmB, The M components of this vector respectively correspond to the M sensors in the ULA of FIGURE 4. If the angular spread between θ and each θ i (see also FIGURE 2) is small, for example 3 ° to 10°, a Taylor series expansion of a(θ,) around θ, = θ yields Equation 3 w(t E t) Λ(f, E{i) Because θ - θ, is small on average (much less than 1 in m Equation 3 that the part of n(t) other than E(t) (i.e., the summation portion) is small on average. Referring again to the example of FIGURE 4, the ULA illustrated therein is an example of an antenna array having a shift invariance structure. This means that each of the sensors 1, 2, ... M is spaced a distance Δ wavelengths apart from its neighboring sensors, and that the array includes two identical sub-arrays. For example, a first sub- array could be defined by the sensors, 1 , 2 and 3 of FIGURE 4, and a second sub-array could be defined by the sensors 2, 3 and 4. The conventional uplink weight vector associated with the antenna of FIGURE 4 includes M components which respectively correspond to the M sensors. Analogously, the weight vector for a given sub-array is a sub-vector of the M-dimensional weight vector, including a plurality of (less than M) components which respectively correspond to the sensors of the sub-array. Accordingly, sub-vectors w, and w 2 respectively associated with first and second sub- arrays can, from Equations 1 and 3 above, be expressed as: Equation 4 w,(t,θ) = h(t,θ)a,(θ) + n,(t) Equation 5 w 2 (t,θ) = h(t,θ)a 2 (θ) + n 2 (t) wherein a, (θ) represents an array response sub-vector of the array response vector a(θ) including those components of a(θ) corresponding to the sensors of the first sub-array, and a 2 (θ) is similarly an array response sub-vector including those components of a(θ) corresponding to the sensors of the second sub-array. Similarly, n,(t) includes the components of n(t) associated with the sensors of the first sub-array, and n 2 (t) represents the components of n(t) associated with the second sub-array. In antenna arrays having a shift invariance structure, for example the ULA of FIGURE 4, the array response sub-vectors of first and second sub-arrays are known to be related as follows: Equation 6 a 2 (θ) = ε ,2nΔs,nβ a,(θ) Substituting Equation 6 into Equation 5 yields the following: Equation 7 w 2 (t,θ) = h(t,θ)a,(θ)e l2πΔs,n0 + n 2 (t). Now, a comparison of Equations 7 and 4 yields: Equation 8 w 2 (t,θ) = w 1 (t,θ)e l2πΔsmθ + »(t) where »(t) is a modified noise term with mean zero. Equation 8 can be solved for θ in many different ways, given a sample of weight vectors (w,(t,θ), w 2 (t,θ)}, t=l to N. The most straightforward example is the least-squares (LS) method described below. Multiplying the left and right sides of Equation 8 with w, H (t) (where H denotes the complex conjugate transpose operator) and summing over the available samples (neglecting »(t)) yields: Equation 9 ∑ w χ H {t)w 2 {t)=∑ w λ "(t) χ {t)e' 2M t= \ (= 1 (the dependence on θ in the notation is dropped here for simplicity). Rearranging Equation 9 produces Equation 10 N , ι2πΔsιnθ . w χ n {t)w { {t) w x H {f)w 2 (t) t = l which can be readily solved for the angle of arrival θ of the multipath finger (for example Θ B of FIGURES 1 and 2). The solution of Equations 9 and 10 is conventionally referred to as the "least-squares" solution because it minimizes the following criterion: Equation 11 -w (t,θ)e' 2πAsmθ f In practical situations the angle-of-arrival θ of a given multipath component varies slowly with time. Then a recursive implementation of the algorithm is preferred. Below is described a recursive least-squares (RLS) method to estimate θ. Note however that the invention is in no way limited to the RLS method, and that other types of recursive solutions are possible, such as least-mean-squares (LMS), Kalman filter or Newton type implementations. To derive the RLS solution to the problem, let α(t) = e l2πΔs ' n °ω According to the above discussion, an estimate of α(t) (including a forgetting factor μ) is given by Equation 12 Introduce the covar ance matrix Equation 13 which clearly satisfies the relation Equation 14 From Equations 12 and 13, Equation 15 t- \ O( = Pit) ∑ μ' ~s w l ) w 2 (s) +w "(t) w 2 (t) s = l H, =α(f - l) + P(t)w x (ήw χ (t)a(t - \) + P(ήw " ( ,ήw 2 (ή H =α(t - l) + P(t)w x \t)(w t) + . {t) (t - \)) The complete recursive least-squares (RLS) algorithm, including a forgetting factor μ can now be summarized as follows: 1 ) Set μ, α(0), and P(0) (see below for suggestions on how to chose these parameters) 2) for t = 1 to N do: Equation 16 α(0 = α(f-l) + K(t)ε(ή Equation 17 where: ε(t) = w 2 {t) - w. (ή (t - l) Equation 18 and K(t) = P(t)w χ H (t) The scalar P(t) used in Equation 18, is obtained from Equation 14. Concerning the choices of the user parameters μ, α(0) and P(0), the following exemplary general rules can be given. The forgetting factor μ determines how long memory is used in the estimation, and should be chosen slightly less than 1 (e.g., between 0.95 and 0.99). It can be shown that a given μ corresponds to estimation using Equation 9 with N = 2/(1 - μ). Thus μ determines how many old samples are used in the estimation. The initial value α(0) is chosen according to any a priori information that might be available on the AOA of the desired user. If no such information is available, α(0) can be set to 1 (corresponding to θ = 0) or any other suitable value. The initial variance P(0) is chosen in accordance with the level of confidence in the initial estimate α(0). A large P(0) gives large initial variation in the estimate and should be used if it is believed that the initial estimate α(0) is bad (e.g. there is no a priori information on the AOA). An exemplary large value is P(0) = 10. On the other hand, if α(0) is believed to be accurate, a small P(0) should be chosen. An exemplary small value is P(0) = 0.01. Depending on the forgetting factor μ, P(t) will converge to different values. If μ is close to one, P(t) will coverage to a small value and tracking will be slow (P(∞) = 0 for μ = 1 which gives no tracking at all), while a small μ will give a higher value of P(t) which enables a faster tracking of changes in α(t). It can be seen from the foregoing discussion that the angle of arrival for a given multipath finger can, according to the invention, be readily estimated from the corresponding uplink weight vector for any antenna array having a shift invariance structure, for example a uniform and linear array. Once the angle of arrival has been estimated, this angle of arrival can be used to suitably control operation of the antenna array for downlink communication to the source of the signal whose angle of arrival has been estimated. FIGURE 5 diagrammatically illustrates exemplary embodiments of an angle of arrival estimator according to the present invention. The angle of arrival estimator 51 includes an input 53 for receiving information indicative of the uplink weight vector w that has been estimated for a given multipath finger, and an output 55 for providing information indicative of the estimated angle of arrival Θ EST . A sub-vector extractor 52 extracts from the weight vector first and second sub-vectors w, and w 2 , as discussed above with respect to Equations 4 and 5. These sub-vectors w, and w 2 are provided as input to an angle of arrival calculator 56, which calculates therefrom and outputs the angle of arrival estimate Θ EST . FIGURE 6 illustrates an exemplary digital radio communications interface according to the present invention. The interface of FIGURE 6 can be provided, for example, in a base station operating in a digital cellular telecommunications network. As shown in FIGURE 6, an antenna array 62 is coupled to duplex (Dx) filters, which, in this example, are in turn coupled to a beamforming network 63. The beamforming network is in turn coupled to receiver (Rx) circuitry, including downconversion circuitry and analog/digital converter circuitry. The receiver circuitry is in turn coupled to an equalization and combining section 65 which provides uplink data to, for example, a data processing section of a base station (not explicitly shown). The section 65 also serves as a weight vector producer that provides the weight vectors w associated with the respective multipath fingers. The angle of arrival estimator 51 receives the weight vectors, produces the angle of arrival estimates Θ EST , and provides these angle of arrival estimates to a downlink beamforming section 66. The downlink beamforming section 66 also receives downlink data from the data processing section. Transmit (Tx) circuitry including digital/analog conversion circuitry, power amplification circuitry and upcon version circuitry is coupled to receive the output of the downlink beamforming section 66, and the output of the transmit circuitry is coupled to the duplex filters. With the exception of the angle of arrival estimator 51 provided according to the present invention, the remaining components illustrated in FIGURE 6 are well known in the art. FIGURE 7 illustrates exemplary operations which can be performed by the exemplary angle of arrival estimator of FIGURE 5. At 71, the weight vector for a multipath finger is obtained. At 72, sub-vectors corresponding to sub-arrays of the antenna array are extracted from the weight vector. At 73 , the angle of arrival estimate is calculated based on the sub-vectors, using any of the above-described exemplary techniques. As indicated at 74 and 75, the procedures at 71-73 can be repeated until all multipath fingers have been processed. The following link simulations demonstrate the performance of the disclosed angle-of-arrival estimation method. Results are shown both for the least-squares solution of Equations 10-11 and for the recursive implementation of Equations 16-18. The results from the least-squares simulations are shown in FIGURES 8 and 9. The root mean-squared error (RMSE) is presented in degrees, defined by Equation 19 Λ as a function of angle-of-arrival θ (θ(n) denotes the nth estimate of θ). The results are averaged over N = 500 Monte Carlo simulations. The RMSE has two origins; first there will be a random error due to the noise. Furthermore, there will be a deterministic error due to the approximation in Equation 3. Indeed, according to the signal model, each spread source is modeled as a large number of point sources impinging on the array. These sources cannot be resolved by the array and an estimation error is the result. All the simulations use a uniform and linear array, ULA, with half wavelength spacing between the elements (Δ = 1/2). FIGURES 8 and 9 show the results for M = 4 and M = 8 antenna sensors respectively, for an angular spread of 3 ° and a signal-to-noise ratio (SNR) of 10 dB. It can be seen that the RMSE becomes larger as the end fire direction is approached, but it never becomes very large in these examples. Also, the RMSE is smaller for the larger antenna array of FIGURE 9. FIGURES 10 and 11 illustrate the ability of the RLS algorithm to track the AOA of a signal. The AOA is modeled to change every 15 samples (corresponding to the 15 slots in a WCDMA radio frame), according to the random walk model Equation 20 θ(t) = θ(t - 1) + v(t) where v(t) is a random variable taking values ±0.5 degrees with equal probability. In reality, the AOAs of the fingers are expected to vary more continuously. The results are shown in FIGURES 10 and 11 for M = 4 and M = 8 antenna elements respectively, for an angular spread of 3° and an SNR of 10 dB. Forgetting factors μ = 0.95 and μ= 0.99 are used. As expected, a lower value of μ allows for a better tracking of the AOA, at the price of a higher variance in the estimates. In any case, the RMSE is low; around one degree in the cases considered. It will be evident to workers in the art that an angle of arrival estimator according to the present invention can be implemented by a suitably programmed data processing device such as a digital signal processor, or by such data processing device in combination with additional external components. Although exemplary embodiments of the present invention have been described above in detail, this does not limit the scope of the invention, which can be practiced in a variety of embodiments.

Claims: WHAT IS CLAIMED IS: 1. A method of estimating an angle at which a multipath finger produced by a digital radio communication source arrives at an antenna array, comprising: receiving the multipath finger at the antenna array; obtaining an uplink weight vector associated with the multipath finger; and calculating an estimate of the angle in response to the weight vector. 2. The method of Claim 1 , including extracting from the weight vector a plurality of sub-vectors respectively corresponding to a plurality of sub-arrays within the antenna array, said calculating step including calculating the estimate from the sub- vectors. 3. The method of Claim 2, wherein said step of calculating the estimate from the sub-vectors includes using a least-squares approximation to calculate the estimate. 4. The method of Claim 2, wherein said step of calculating the estimate from the sub-vectors includes using a recursive algorithm to calculate the estimate. 5. The method of Claim 4, wherein said step of using a recursive algorithm includes using a recursive least-squares approximation to calculate the estimate. 6. The method of Claim 1 , wherein said receiving step includes receiving the multipath finger at an antenna array having a shift invariance structure. 7. The method of Claim 6, wherein said receiving step includes receiving the multipath finger at a uniform and linear antenna array. 8. The method of Claim 1, wherein the multipath finger is produced by a mobile telecommunication device operating in a cellular telecommunication network. 9. An apparatus for estimating an angle at which a multipath finger produced by a digital radio communication source arrives at an antenna array, comprising: an input for receiving an uplink weight vector associated with the multipath finger; and an angle-of-arrival calculator coupled to said input for calculating an estimate of the angle in response to the weight vector. 10. The apparatus of Claim 9, including an extractor coupled between said input and said angle-of-arrival calculator for extracting from the weight vector a plurality of sub-vectors respectively corresponding to a plurality of sub-arrays within the antenna array, said angle-of-arrival calculator operable for receiving said sub- vectors from said extractor and calculating therefrom the estimate. 11. The apparatus of Claim 10, wherein said angle-of-arrival calculator is operable for using a least-squares approximation to calculate the estimate. 12. The apparatus of Claim 10, wherein said angle-of-arrival calculator is operable for using a recursive algorithm to calculate the estimate. 13. The apparatus of Claim 12, wherein said recursive algorithm is a recursive least-squares algorithm. 14. The apparatus of Claim 9, wherein the antenna array has a shift invariance structure. 15. The apparatus of Claim 14, wherein the antenna array is a uniform and linear array. 16. The apparatus of Claim 9, wherein the digital radio communication source is a mobile telecommunication device operating in a cellular telecommunication network. 17. A digital radio communication apparatus, comprising: an antenna array for receiving a multipath finger produced by a digital radio communication source; a weight vector producer coupled to said antenna array for producing an uplink weight vector associated with the multipath finger; and an angle-of-arrival calculator coupled to said weight vector producer for calculating, in response to the weight vector, an estimate of an angle at which the multipath finger has arrived at said antenna array, said angle-of-arrival calculator including an output coupled to said antenna array for controlling downlink operation of said antenna array based on the estimate. 18. The apparatus of Claim 17, including an extractor coupled between said weight vector producer and said angle-of-arrival calculator for extracting from the weight vector a plurality of sub-vectors respectively corresponding to a plurality of sub-arrays within the antenna array, said angle-of-arrival calculator operable for receiving said sub-vectors from said extractor and calculating the estimate therefrom. 19. The apparatus of Claim 18, wherein said angle-of-arrival calculator is operable for using a least-squares approximation to calculate the estimate. 20. The apparatus of Claim 18, wherein said angle-of-arrival calculator is operable for using a recursive algorithm to calculate the estimate. 21. The apparatus of Claim 20, wherein said recursive algorithm is a recursive least-squares algorithm. 22. The apparatus of Claim 17, wherein said antenna array has a shift invariance structure. 23. The apparatus of Claim 22, wherein said antenna array is a uniform and linear array. 24. The apparatus of Claim 17, wherein said digital radio communication source is a mobile telecommunication device operating in a cellular telecommunication network.

Inventors: CEDERVALL MATS, GOERANSSON BO, SORELIUS JOAKIM, CEDERVALL, MATS, GOERANSSON, BO, SORELIUS, JOAKIM
Applicants: ERICSSON TELEFON AB L M, TELEFONAKTIEBOLAGET LM ERICSSON (PUBL)
Application year: 2001
Main classification: H01Q
IPC codes: H04B7, H04W64