Field of the Invention The present invention relates generally to communication systems, and more particularly, the present invention relates to a method and apparatus for determining the angle of arrival of a transmitted signal in a communication system.

Background of the Invention It is well known that a remote unit's location within a communication system may be determined using an amplitude difference-based angle of arrival estimation (AD-AOA). Such a method of AD-AOA location using average energy values for the received prompt rays is described in commonly assigned U. S. Pat. No. 5,786,791, METHOD FOR DETERMINING AN ANGLE OF ARRIVAL OF A SIGNAL TRANSMITTED BY A REMOTE UNIT IN A COMMUNICATION SYSTEM, by Bruckert, the disclosure of which is hereby expressly incorporated herein by reference.

Other methods are described in U. S. Pat. No. 4,636,796, RADIO DIRECTION FINDING SYSTEM, by Imazeki and U. S. Pat. No. 4,833,478, AUTOMATIC DIRECTION FINDER ANTENNA ARRAY, by Nossen, the disclosures of which are hereby expressly incorporated herein by reference. Briefly, these later methods determine the angle of arrival of a signal transmitted from a remote unit by analyzing the amplitude differences of signals received at multiple antennas at a base site.

In accordance with known methods of locating a remote unit transmitter in code division multiple access (CDMA) based communication systems, a high precision estimation of the prompt ray time-delay and amplitude is required. High precision estimation of the received signal amplitude at the main sector can be achieved by the location searcher

with the main sector base station based on the method for prompt ray extraction with coherent integration, deconvolution and non-coherent integration procedures. However, the low received chip-to-noise ratio Ec/No at the adjacent sector (7 to 15 dB less than at the main sector) makes it impossible to use the location searcher at the adjacent sector to provide accurate AOA estimate.

Higher precision estimations may be obtained at the adjacent sector using the methods and apparatus described in the commonly assigned U. S. Pat. Application METHOD AND APPARATUS FOR ESTIMATING A CHANNEL METRIC by Henderson, et al. filed of even date herewith (attorney docket no. CE08056R-Henderson), the disclosure of which is hereby expressly incorporated by reference. Still, low signal- to-noise ratio (SNR) conditions can hamper the effectiveness of the location searcher in determining the AOA even given a higher precision estimation.

Therefore, there is a need for an improved method and apparatus for AD-AOA estimation in the presence of low signal-to-noise ratio (SNR) conditions at the adjacent sector location searcher.

Brief Description of the Drawings FIG. 1 is a schematic illustration of a communication system.

FIG. 2 is a block diagram illustration of a communication system incorporating a location searcher.

FIG. 3 is a block diagram illustration of a communication system incorporating a modified location searcher.

FIG. 4 is a block diagram illustration of a communication system incorporating a further modified location searcher.

FIG. 5 is a block diagram illustration of a communication system incorporating a location searcher adapted in accordance with a preferred embodiment of the present invention.

FIG. 6 illustrates several plots depicting convolution and non- convolution channel outputs.

FIG. 7 illustrates several plots depicting convolution and non- convolution channel outputs.

FIG. 8 illustrates the improvement in the amplitude estimation accuracy as a function of SNR.

Detailed Description of the Drawings Accurate estimations of the angle-of-arrival (AOA) of the prompt ray at an adjacent sector of a base station are provided by herein described methods and apparatus according to preferred embodiments of the invention. The invention encompasses a method of estimating an angle- of-arrival in which high precision parameter estimates from a main sector location searcher are provided to an adjacent sector location searcher in order to enhance the accuracy of the parameter estimates generated by the adjacent sector location searcher. In a preferred implementation of the method, the main sector time-delay estimates are used as time-delay offset estimates for the adjacent sector amplitude estimation process. In an alternate preferred implementation of the method, main sector time delay and channel phase (hence frequency) estimates are used to increase the coherent averaging window within the adjacent sector location searcher to in turn increase the accuracy of the amplitude estimation from the adjacent sector location searcher.

In a preferred apparatus, a main sector location searcher and an adjacent sector location searcher within a base station of a communication system are coupled to communicate parameter estimates. In a preferred implementation of the apparatus, an output of the main sector location searcher time-delay estimator is coupled to the adjacent sector location searcher to provide a dispreading time offset estimate. In an alternate preferred implementation of the apparatus a channel estimation result

from the main sector location searcher is multiplie with the dispreading signal in the adjacent location searcher.

Referring to the communication system shown in FIG. 1, a sectorized base site 101 includes a plurality of receive antennas 103-113 which receive signal 115 transmitted from remote unit 117. As is evident, antennas 103 and 105 in a main sector"a"will receive signal 117 at a different angle of arrival than antennas 107 and 109 in an adjacent sector " (3." Because of this, antennas in the main sector will also receive signal 117 at a different amplitude than antennas in the adjacent sector. It is the difference in received amplitudes of signal 117 at various antennas that is utilized in determining the angle of arrival of signal 117. Adjacent is not used herein to refer to the next physically adjacent sector or sectors to the main sector. Instead, adjacent is used to refer to sectors that are not otherwise aligned with the main sector from an azimuth perspective whether or not physically located next to the main sector.

As shown in FIG. 2, each of the main sector and the adjacent sector include an RF front end 201 and 203 coupled to antennas 103 and 105 and to antennas 107 and 109, respectively. Each RF front end 201 and 203 operate as is well known to operate on the received signals and to generate baseband 1, Q signals which in turn are communicated, respectively, to a main sector location searcher 205 and an adjacent sector location searcher 207. The output of each of the searchers 205 and 207 is an estimated amplitude of the signal 117 arriving at the respective antennas. These estimates are fed to an AOA computer 209 from which the AD-AOA estimate is computed.

The gain of the antennas 103-109 is directional and varies as a function of the azimuth and elevation angles. Thus, two adjacent sectors antennas with angular diversity can be represented as a converter of azimuth angle-of-arrival to the amplitude difference and can be used to determine the line of bearing between the remote unit 115 and the base station 101. Applying Gaussian curve approximation to the horizontal and

vertical patterns of the two adjacent sectors antennas (antennas 103/105 and antennas 107/109) the amplitudes of the received signals Ai and A2 can be expressed by: A1=Am1exp{-[(#-#1)2/0/72B2hl] = [atan (h,/r) 2/0. 72B2v1} (1) A2=Am2exp{-[(#-#2)2/0/72B2h2]=[atan(h2/r)2/0.72B2v2} (2) where: (pi and (P2 are the central angle corresponding to the antennas' boresight; (p is a current angle of arrival relative to the boresight; Bh, and Bh2 are the 3db beamwidth of the horizontal antenna pattern; Bvr and Bv2 are the 3db beamwidth of the vertical antenna pattern; Ami and Am2 are the ultimate attenuations of the antennas with heights hX and h2, respectively ; and r is a distance between the base station and the remote unit.

The amplitude ratio calculated in the AOA computer 209 for identical antenna parameters (Bh1=Bh2, Bv1=Bv2, Am1=Am2) may be written as: An/A2=exp{-[(#-#1)2/0/72B2h]+ (#-#2)2/0/72B2h} (3) from which: sp = { [0. 36 B2h ln(A1/A2)]/ (#1-#2)}+0. 5 (T, +T2)] (4)

Equation (4) is the AD-AOA transfer function of the adjacent sectors antennas with the horizontal patterns beamwidth Bh, the angular diversity ((P-(P2) and the equal gain azimuth angle 0.5 (cp+cp2).

The AOA estimation error due to the errors in the amplitude A1 and A2 estimation can be derived as a differential operation, which results in: d#=S0d(A1/A2)/(A1/A2) (5) where So=0.36 B2h/((p-ep2). For small values of errors, differentiation can be replaced by a difference, which yields an upper bound on the AOA error as a function of the relative errors of the amplitude estimation AAI/Al and AA2/A2 in the location equipment of the two adjacent sectors, respectively: ##<={|#A1/A1|+|#A2/A2|}/0.36 B2h/(#1-#2).

(6) In general to evaluate AOA estimation accuracy it is necessary to treat A (p as a stochastic error with variance : Var (##) = So {var (AAi/Ai) + var (#A2/A2)-2R12 [var (#A1/A1) var (AA2/A2)] ° 5} (7) where R12 is the correlation coefficient of the random variables AAi/Ai and AA2/A2, which depends on the correlation signals in the two angular diversity antennas with small space separation. For R12 the phenomenon of mutual amplitude errors compensation occurs resulting in a decrease in the AOA error.

FIG. 3 illustrates a first modification to the location searcher shown in FIG. 2 and described above. The main sector location searcher 301

and the adjacent sector location searcher 302 work in two phases. The signals received from the main sector antennas 103/105 and the adjacent sector antennas 107/109 are despread 303,305 and coherently averaged 307,309 over N, and N2 Walsh symbols, respectively. An estimate of the maximum Doppler frequency is used to determine appropriate values for the coherent averaging window Nopti and Nopt2, respectively, which are used for coherent averaging. In a second stage, the signals are respectively non-coherently averaged 311,313 over M, and M2 averaging periods, respectively. These procedures are done for several time offsets of the dispreading PN sequences separated by one eighth of a chip period to build a power delay profile (PDP) with the results being retain in a memory (not depicted).

From the memory, the time-delay estimators 315 and 317 and amplitude estimators 319 and 321 select the maximum of the PDPs and assign the time offset and square root values of these offsets to the time of arrival (Ti and T2) and signal amplitude (Ai and A2) estimates, respectively.

The estimates are then fed to the AOA computer 323, which forms the ratio Ai/A2 and the AOA estimate (p=So (A1/A2).

This technique, however, is not able to resolve the prompt ray in the main or the adjacent sectors and requires very long total integration time for the adjacent sector (Nopt2*M2>>Nopt1*M1) due to the low input SNR.

Thus, the accuracy of the AOA estimation leaves much to be desired.

FIG. 4 illustrates a main sector location searcher 401 and an adjacent sector location searcher 402. The location searchers 401 and 402 are identical to location searchers 301 and 302, and like reference numerals are used to depict like elements, with the exception that the main sector location searcher 401 further includes a deconvolver 403.

The location searcher 401 produces a PDP using the technique described above with respect to FIG. 3. Using this PDP a new window is formed around its peak (i. e., Wc=2Tc). The received signal from the main sector antennas 103/105 is despread for the time offset in the window Wc

and then integrated coherently 307, deconvolved 403 with the known pulse shape and non-coherently integrated 311. The deconvolution 403 is carried out after the coherent integration process 307 where the SNR has been enhanced significantly (14-21 db for maximum Doppler frequency range from 80Hz-100Hz). Thus, location searcher 401 is able to resolve sub-chip spaced rays and achieve high precision estimation of the prompt ray time-delay T1 and amplitude Ai.

However, the location searcher 402 does not use a deconvolver to produce the amplitude estimate A2 of the peak ray (maximum PDP) due to the low input SNR. Thus, the errors of the amplitude and final AOA estimates are not decreased significantiy as compared to the technique described in connection with FIG. 3.

FIG. 5 illustrates a main sector location searcher 501 and an adjacent sector location searcher 502. It will be appreciated that the location searchers 501 and 502 may be implemented as software on magnetic or optical media, firmware, a programmable gate array and/or an application specific integrated circuit without departing from the fair scope of the invention. In accordance with a preferred embodiment of the invention, high precision parameters estimates obtained from results of the main sector location searcher 501 calculations are used to aid the adjacent sector location searcher 502. This can be accomplished because of the high signal correlation between the main sector antennas 103/105 and the adjacent sector antennas 107/109. This is especially true where the separation is small, for example, approximately 1 m for a typical cellular tower top antenna configuration. The location searcher 501 estimates the time-delay and amplitude using the techniques described in connection with FIG. 4, and like reference numeral are again used to identify like elements between FIGs 3,4 and 5. The location searcher 501 further inclues a max Doppler estimator 519 coupled to the despreader 303 for providing an estimate coherent averaging window for coherent averaging 307.

Assuming a two ray multipath channel, the PDP observed by the location searcher 501 can be expressed as PDP1=α211a(T-T11) + a 2, a (T-T2,) Where an, Tn, and a21 and T21 (T11<T21) are the amplitudes and time-delays of the prompt and peak rays in the main sector antennas 103/105, and a (r) is the impulse response of the chip pulse-shaping filters. Due to the deconvolution provided by the deconvolver 403, the location searcher 501 is able to resolve the prompt ray and achieve high precision estimation for the parameters in and an. The location searcher 502 observes the PDP PDP2=α212a(T-T12) + a222a (T-T22) where a12, T12, and a22, T22 are the amplitudes and time-delays of the prompt and peak rays in the adjacent sector antennas 107/109.

It is not possible to resolve and estimate the prompt ray parameters ai2 and T12 for IT22-Tr2l<Tc without the operation of a deconvolver, which can not be implemented in the adjacent sector location searcher because of the low input SNR. Thus, the location searcher 502 produces estimates a2 and T2 whose values are closer to the actual parameters of the peak ray a22 and T22 than the actual parameters of the prompt ray a, and T12. This is illustrated in the plots shown in FIG. 6 and FIG. 7, wherein the estimate a2 and T2 are more closely aligned with the actual parameters of the peak ray a22 and T22.

In order to reduce the amplitude estimation error, the location searcher 502 observes only one value of the PDP2 corresponding to the unknown prompt ray time-delay T12 using the location searcher 501 estimation in. Due to small separation between the main and adjacent

sector antennas, r » is close to T12 and therefore, in can be used as an estimate to assign the time offset for dispreading the received signal in the location searcher 502. In this regard, the estimate in is fed to a PN offset assignor 505 in the location searcher 502, which assigns a Walsh symbol duration to the despreader 305. After dispreading 305, averaging 309 and 313 and bias correction 507, location searcher 502 produces an estimate a, 2, which is close to the actual amplitude a, 2 of the prompt ray without its extraction from the total PDP.

The location searcher 501 includes a channel estimator 509 that generates a phase estimation result (pi, and this phase estimation is used to increase the coherent averaging window Nopt2 resulting in increased SNR processing gain and hence in increased accuracy of the amplitude estimation in the location searcher 502. The dispreading signals plus noise at the location searcher 501 and the location searcher 502, respectively, can be written xi (t) =ai (t) expj (t) +n, (t) (10) X2 (t) =a2 (t) expj (p2 (t) +n2 (t).

(11) The channel estimation results from the location searcher 501, after amplitude fluctuation suppression 511 and conjugation operation 513, is multiplied with the dispreading signal in the location searcher 502 by multiplier 515. The product can be written x3(t)=x2(t)x'1(t)= α2(t)expj(#2-#1)+n'2(t).

(12)

Due to typically small separation between the main sector antennas 103/105 and the adjacent sector antennas 107/109 and high correction between x1 (t) and x2 (t), (pipi and as result phase rotation may be removed from equation (12) yielding X3 (t) =a2 (t) +n'2 (t).

(13) Therefore, the coherent averaging window can be increased as compared to that used by location searcher 501 (Nopt2>>Nopt1). This is referred to as quasi-coherent averaging. After coherent averaging 309 and non- coherent averaging 313 with a small window (M2 « Mi) and bias correction 507, location searcher 502 produces an amplitude estimate a2 with high precision. The plots in FIG. 8 illustrate the improvement in the amplitude estimation accuracy as a function of SNR. For a given value of the mean- squared error, for example-28dB, and the correlation coefficient between xi (t) and X2 (t) R12=0.9, the location searcher 502, using only coherent averaging (N=2304), requires about 10dB less input SNR than a location searcher with dual coherent (N=36) and non-coherent (M=64) configuration. For given SNR=-30dB (a typical ratio for the adjacent sector) and R12=0.9, the location searcher 502, again using only coherent averaging, provides about 5dB less mean-squared error of the amplitude estimation than a location searcher with a dual averaging algorithm and the same total sample integration period.

While the invention has been particularly shown and described with reference to several preferred embodiments, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention, and it is intended that all such changes come within the scope of the following claims.