The present invention relates generally to RAKE receivers for code division multiple access (CDMA) systems, and more particularly, to a method and apparatus for determining the finger placement of RAKE fingers in a RAKE receiver.
In a wireless communication system, a transmitted signal may travel through multiple propagation paths so that the received signal is a composite of multiple time-shifted versions of the signal. The different time-shifted versions of the received signal, referred to herein as signal images, suffer from different phase and attenuation effects. The multiple time-shifted signal images combine in an unpredictable manner resulting in signal fading.
CDMA receivers typically employ a RAKE receiver to combat signal fading due to multi-path propagation. The goal of the RAKE receiver is to detect the individual signal images and combine them coherently. A RAKE receiver typically includes a plurality of correlators, sometimes referred to as fingers, to separately despread different time-shifted signal images, and a combiner to combine the correlator outputs. For example, a RAKE receiver may detect and combine the M strongest signal images. A path searcher processes the received signal to identify the delays corresponding to the strongest signal images, and a finger placement processor determines the finger placement based on those delays. The process of finger placement comprises the assignment of a delay to each RAKE finger to align the RAKE finger in time with a signal image. A simple finger placement strategy is to assign the delays of the J strongest signal images found by the path searcher to respective RAKE fingers.
Typically, finger placement starts by generating an estimated power delay profile (PDP) over a defined search window that gives the signal power as a function of delay. The path searcher measures the signal power of the received signal samples. The spacing between samples defines a search grid and the signal power measurements define the PDP. One approach to finger placement, referred to herein as the “peak” approach, is to place fingers at or near the peaks or local maximas in the PDP. Ideally, the RAKE fingers would be placed at the exact delays corresponding to peaks in the PDP. Exact placement of the RAKE fingers at the peaks of the PDP is not always possible because the search grid does not always align with the peaks in the PDP. The actual peak in the PDP may fall between the grid points of the search grid.
In some circumstances, the path searcher may detect a dominant signal image in the received signal. This situation may occur, for example, when there is a direct line-of-sight path between the transmitter and the receiver, when there is a non-dispersive channel, or when one signal image is much stronger than other signal images. When a dominant signal image exists, the RAKE fingers are typically placed on a finger grid centered on the delay of the dominant signal image reported by the path searcher. Delay hypothesis testing can be used to determine the relative timing delay between the current finger placement and the dirac-like channel impulse response in order to track the path of the dominant signal image. Knowing this relative delay allows for interpolation of despread values at non-sample positions thus increasing the SINR at the output of the combiner compared to non-interpolating schemes.
Conventional methods for delay hypotheses testing suffer from a high computational complexity to achieve a given accuracy with respect to the relative delay estimate. The reason for the computational complexity is that the number of hypotheses to be tested depends on the uncertainty interval of the path searcher result, which is at least based on the sampling period distance and on the required delay resolution. Assuming typical values, i.e., a sampling period distance of ¼ chip duration and a required delay resolution of 1/32 chip to support 64-QAM, 9 different hypotheses needs to be tested.
One exemplary embodiment of the present invention provides a low complexity method for adapting a current path delay estimate for path delay tracking. The starting path delay estimate may be a previously determined path delay estimate or may be set to the most recent path searcher output to initialize or re-initialize path delay tracking. The method begins by selecting an operating mode. Next, a metric is computed for each hypothesized path delay. Finally, the current estimated path delay is adapted based on the computed metrics for the hypothesized path delays.
In one exemplary embodiment, two modes are contemplated: fine mode and coarse mode. In fine mode, a hypothesized path delay selected based on the computed metrics may be filtered to obtain a filtered delay estimate. In this case, the new path delay estimate is set equal to the filtered delay estimate. In coarse mode, a hypothesized path delay selected based on said computed metrics and an adjacent hypothesized path delay may be interpolated to obtain an interpolated delay estimate. In this case, the new path delay estimate is set equal to said interpolated delay estimate. Interpolation may be conditionally performed depending on whether the selected hypothesized path delay is on an edge of a test window. Also, those skilled in the art will readily appreciate that interpolation may also be used in fine mode to improve accuracy of path delay estimation.
Another exemplary embodiment of the invention provide a receiver having a processor configured to perform delay hypothesis testing for delay tracking. More specifically, the processor is configured to select an operating mode (e.g., fine or coarse), determine a set of hypothesized path delays based on a current path delay estimate and the operating mode, compute a metric for each hypothesized path delay, and adapt said path delay estimate based on said computed metrics for said hypothesized path delays. In some embodiments, the receiver may comprise a RAKE receiver. In other embodiments, the receiver may comprise a chip equalization receiver.
The present invention relates to a method and apparatus for determining finger placement in a RAKE receiver or tap delays in a chip equalizer receiver. As used herein, the term RAKE receiver includes a generalized RAKE (G-RAKE) receiver as described in U.S. Pat. No. 6,363,104, which is incorporated herein by reference. The invention has application in single-input single-output (SISO) receivers, multiple-input, single-output (MISO) receivers, single-input, multiple-output (SIMO) receivers, and multiple-input, multiple-output (MIMO) receivers.
RAKE processor 30 determines the number and placement of the RAKE fingers 16 and the combining weights applied to the individual RAKE finger outputs by weighting elements 24.
Finger placement processor 34 includes a finger locator 35, channel estimator 36, combining weight generator 38, and metric calculator 40. The finger locator 35 determines the locations of the RAKE fingers 16 based on the estimated path delays reported by the path searcher 32. In a conventional RAKE receiver, a RAKE finger 16 is typically assigned to the L strongest delays. In a G-RAKE receiver, additional RAKE fingers 16 may be used to detect the received signals that do not correspond to any path delay. The channel estimator 36 generates estimates of the propagation channel from the transmitter to the receiver 5 for each assigned RAKE finger 16. The channel estimates are provided to the combining weight generator 38. Combining weight generator 38 computes the combining weights applied to the RAKE finger outputs. The combining weight generator 38 may, for example, compute combining weights based on a maximal ratio combining (MRC) criteria. For MRC combining, the combining weights are based upon the signal power or signal power to interference power ratio (SIR) at the output of the correlator 20 for each RAKE finger 16. If the SIR for a particular RAKE finger 16 is low, it will be assigned a low weighting factor. Conversely, if the SIR for a particular RAKE finger 16 is high, it will be assigned a large weighting factor. With full G-RAKE combining, the combining weight generator 38 computes impairment correlations across the RAKE fingers 16 and generates an impairment covariance matrix R. The combining weight generator 38 multiplies the vector of channel estimates ĥ from the channel estimator 36 by the inverse of the impairment covariance matrix R to generate a weight vector w, whose elements are the weighting factors for the outputs of RAKE fingers 16.
In some circumstances, path searcher 32 may detect a dominant signal image in the received signal. This situation may occur, for example, when there is a direct line-of-sight path between the transmitter and the receiver. A dominant signal image will exist if there is only one detected signal image indicating a non-dispersive channel. If the channel is dispersive and more than one signal image is detected, the presence of a dominant signal image may be determined based on the relative signal power or SIR of the detected signal images. If the signal power or SIR of the strongest signal image exceeds all others by a predetermined amount, then a dominant signal image may be deemed to exist.
When a dominant signal image exists, RAKE fingers 16 are typically placed on a finger grid centered on the delay of the dominant signal image reported by the path searcher 32. Delay hypothesis testing may be used to determine the relative timing delay between the current finger placement and the dirac-like channel impulse response in order to track the path of the dominant signal image. An exemplary method of delay hypothesis testing is described in U.S. Published Patent Application No. 2006/0268962, which is incorporated herein by reference. The present invention provides a low complexity method of performing delay hypothesis testing.
In fine mode, a hypothesized path delay selected based on the computed metrics may be filtered to obtain a filtered delay estimate. In this case, the new path delay estimate is set equal to the filtered and possibly quantized delay estimate. If quantization is used, the quantization levels may have a higher resolution than the sampling period. In coarse mode, a hypothesized path delay selected based on said computed metrics and adjacent hypothesized path delay may be interpolated to obtain an interpolated delay estimate. In this case, the new path delay estimate is set equal to said interpolated delay estimate. Interpolation may be conditionally performed depending on whether the selected hypothesized path delay is on an edge of a test window. Also, those skilled in the art will readily appreciate that interpolation may also be used in fine mode to improve accuracy of path delay estimation.
For illustrative reasons and without loss of generality it is assumed that six fingers (F=6) spaced at a distance of
are used for CPICH cross correlation of the WCDMA Common Pilot Channel (CPICH), where TC is the chip period. The sampling interval is
The Rake fingers 16 are initially placed around the path searcher result as shown in
The processing is then started by setting the operating mode to coarse (block 104) and initializing the current path delay estimate τmin(s) (block 106).
Before delay hypotheses testing begins, net channel estimates ĥ are computed by removing the pilot sequence s(m) from the CPICH cross-correlation results x(m) (block 108). Averaging over e.g., ten subsequent estimates yields an estimated channel response vector ĥ(s) as follows:
Like x(m), the vector ĥ(s) is a complex-valued column vector having six elements.
The number K and a delay resolution or distance D (normalized with respect to the chip period TC) of delay hypotheses tests is initialized depending on the operating mode (block 110). In coarse mode, which is the initial mode, the distance D is generally larger to cover a wider uncertainty interval with a comparable (or equal) number of tests as in fine mode.
in coarse mode and KF=5 and
in fine mode demonstrates realistic figures. The uncertainty interval covered amounts to three-fourths of a chip period
in coarse and one eighth of a chip period
in fine mode. A need for that may either be caused by bad path searcher accuracy, especially in case of MIMO. Since the path searcher 32 generally determines the delay between a specific transmitter/receiver (Tx/Rx) antenna pair, radio frequency (RF), or propagation delay differences in other Tx/Rx pairs have to be taken into account as additional uncertainty. Note that due to these differences in case of MIMO, delay hypotheses testing needs to be done for each Tx/Rx antenna combination independently. In fine mode, the test window is reduced compared to coarse mode since the channel path delay is generally slowly time varying.
A variable C is also initialized independently for both the coarse and fine operating modes (block 112). The variable C is a threshold that is used to determine when delay hypothesis testing needs to be re-initialized as hereinafter described.
Delay hypothesis testing begins by setting the loop counter k=0 (block 114). Hypothetical delays τhyp,k(s) for current slots (slot number s) are set depending on the initialization of K and D (block 116). Hypotheses testing is centered around the previously estimated, quantized path delay τmin(s−1) at slot number s−1:
Quantization is performed with an even finer resolution than fine delay hypotheses testing. Thus, the maximum error may be further decreased from 1/64 to 1/128 without an increase in complexity.
The set of hypothetic delays τhyp,k(s) for all kε{0, 1, . . . , K−1} is depicted in
A medium channel coefficient g(k) is computed for each hypothetical delay τhyp,k(s) according to:
gk=rk,invĥ(s), (4)
where rk,inv is pseudo inverse of the real valued vector of overall transmit and receiver filter impulse response rk at each finger position (block 118). The pseudo inverse rk,inv of the filter impulse response may be computed according to:
rk,inv=(rkHrk)−1rkH. (5)
The vector rk is defined by:
is the autocorrelation function for the combined transmit and receive filters. Due to quantization of τhyp,k(s) to multiples of 1/64, the number of different vectors rk,inv is limited and may be read from a look-up table.
The hypothetical net channel coefficients {tilde over (h)}k for the current hypothetical delay τhyp,k are computed using the previously estimated medium channel coefficient g(k) and the overall filter response rk (block 120). The hypothetical net channel coefficients {tilde over (h)}k may be computed according to:
{tilde over (h)}k=gkrk. (8)
A distance metric ρk is computed for each delay hypothesis (block 122). The distance metric may be computed according to:
ρk=(ĥ(s)−{tilde over (h)}k)HRN−1({tilde over (h)}(s)−{tilde over (h)}k) (9)
using the inverse of the autocorrelation matrix of thermal noise RN−1. Those skilled in the art will appreciate that other metrics may be used in place of or in addition to the distance metric ρk to evaluate delay hypotheses. After the metrics are computed, the computed metrics are used to adapt the current path delay estimate τmin(s).
Adaptation of the current path delay estimate τmin(s) begins by selecting a best path delay estimate {circumflex over (τ)}min(s) based on the computed metric ρk (block 124). The best path delay estimate {circumflex over (τ)}min(s) is defined by:
Applying equation (10) in coarse mode provides more inaccurate estimates due to the larger hypotheses test distance. However, this initial inaccuracy may be reduced considerably by means of interpolation. Provided that ρk−1,ρk+1 exists, e.g., {circumflex over (k)}≠0 and {circumflex over (k)}≠K−1, an interpolated path delay estimate {circumflex over (τ)}int(s) may be computed by adding an offset τoffs(s) to the best path delay estimate {circumflex over (τ)}min(s). The offset τoffs(s) is defined by:
and the interpolated path delay estimate becomes:
{circumflex over (τ)}int(s)={circumflex over (τ)}min(s)+τoffs(s). (12)
Equation (10) is only used as fallback in the coarse operating mode in case the best path delay estimate {circumflex over (τ)}min(s) falls on the edge of the test window, e.g., {circumflex over (k)}=0 or {circumflex over (k)}=K−1. If desired, interpolation may also be performed in fine mode to improve the accuracy of the estimate.
After the best path delay estimate {circumflex over (τ)}min(s) or interpolated delay estimate {circumflex over (τ)}int(s) is determined, the internal counter variable c is incremented if the instantaneous (best or interpolated) path delay estimate {circumflex over (τ)}min(s) or {circumflex over (τ)}min(s) from the delay hypotheses testing lies at the edge of the test window, or it is reset to zero otherwise (block 128). If C subsequent delay estimates are located at the edge of the test window, the current assumption for the dominant path position may be wrong and the entire delay hypothesis testing procedure is reinitialized. This means that the finger position delivered from the path searcher 32 will be the new starting point for next round of delay hypotheses testing. Note that the parameter C may be configured independently for each mode. For coarse mode, the variable C=CC and for fine mode the variable C=CF.
In the fine operating mode, the best path delay estimate {circumflex over (τ)}min(s) or interpolated estimate {circumflex over (τ)}int(s) obtained from delay hypothesis testing is filtered by an infinite impulse response (IIR) filter or low pass filter to obtain a filtered path delay estimate {circumflex over (τ)}fil(s) (block 130). The filter process may be given by:
{circumflex over (τ)}fil(s)={circumflex over (τ)}fil(s−1)+α[{circumflex over (τ)}min(s)−{circumflex over (τ)}fil(s−1)]. (13)
or alternatively by:
{circumflex over (τ)}fil(s)={circumflex over (τ)}fil(s−1)+α[{circumflex over (τ)}int(s)−{circumflex over (τ)}fil(s−1)] (14)
The new path delay estimate τmin(s) is set equal to the filtered path delay estimate {circumflex over (τ)}fil(s) to complete the adaptation of the current path delay estimate (block 131).
In coarse mode, the new path delay estimate τmin(s) is set equal to the interpolated path delay estimate {circumflex over (τ)}int(s) if available or the best path delay estimate {circumflex over (τ)}min(s) if an interpolated path delay estimate {circumflex over (τ)}min(s) is not available (block 132). Note that setting the new path delay also initializes the filter for the next iteration of delay hypothesis testing. Further note that as long as delay hypotheses testing remains in coarse mode (when {circumflex over (τ)}min(s) is at the edge of the test window), the IIR filter is reinitialized over and over again until {circumflex over (τ)}min(s) lies inside the test window. If in coarse operating mode and the delay estimate {circumflex over (τ)}min(s) is inside the test window, e.g., {circumflex over (k)}≠0 or {circumflex over (k)}≠K−1, the operating mode is changed to fine (block 134).
Before the next iteration of delay hypothesis testing starts, the new path delay estimate τmin(s) is compared to the path searcher output and adaptation of the finger positions is made if needed. Note that the rate of adaptation may be greater than the rate at which the path searcher 32 outputs new path delay estimates. For example, adaptation may take place once per slot and the path searcher 32 may output a new estimate every 240 slots. Finger adaptation begins by comparing the distance TD measured in time between the most recent path delay output by the path searcher 32 and the new path delay estimate τmin(s) (block 136), and comparing the computed distance to the chip period TC (block 138). If the delay is greater than one half the chip period TC, the current assumption for the dominant path position may be wrong and the next round of delay hypothesis testing begins by re-centering the finger grid on the output of the path searcher (block 102). Next, it is checked whether the absolute value of the new path delay estimate τmin(s) exceeds ¼ of a chip period TC (block 140). When the absolute value of the new path delay estimate τmin(s) exceeds ¼ of a chip period TC, the true, e.g., denormalized, value exceeds the sampling distance Ts. In this case, a finger shift is invoked externally and the new path delay estimate τmin(s) is adapted accordingly (blocks 142, 144). To be more precise, if τmin(s)>¼, the fingers shall be shifted by +Ts, and if τmin(s)<−¼, the fingers shall be shifted by −Ts (block 144) and τmin(s) is recomputed by:
Following the finger adaptation, the next round of delay hypothesis testing begins with the computation of net channel coefficients (block 108).
The present invention provides a reduction in computational complexity and a reduction of the maximum and mean delay estimation error. The reduction in complexity is achieved using two operation modes (coarse and fine) with different resolutions for delay hypotheses testing. Initially, when there is a large uncertainty where the true channel path is located, the coarse mode allows for evaluating this large uncertainty interval with relatively few hypotheses tests. For example having an uncertainty interval of ± a quarter of a chip period would require seventeen hypothesis tests to be performed if the high resolution testing grid with a distance of 1/32 chip period is used from the beginning. Starting with a grid resolution of ⅛ chip period reduces the number of required tests. The inherent loss in resolution is overcome by interpolation between distance metrics. Positioning the delay hypotheses test window around the current (filtered) delay estimate also allows for minimizing the number of hypotheses to be considered.
A reduction of the maximum and mean delay estimation error is achieved by using a higher resolution for the position of the testing grid than for the testing grid itself. For example, using a testing grid with a resolution of 1/32 chip period and positioning this grid only at delays being an integer multiple of this period yields a maximum potential error of 1/64 chip period (half the resolution). Letting the grid be positioned at multiples of 1/64 chip period, for example, reduces the maximum potential error to 1/128 of a chip period without needing more hypothesis tests.
Those skilled in the art will appreciate that G-RAKE processing is equivalent to linear minimum mean square error (LMMSE) based chip equalization (CE).
The present invention may, of course, be carried out in other specific ways than those herein set forth without departing from the scope and essential characteristics of the invention. The present embodiments are, therefore, to be considered in all respects as illustrative and not restrictive, and all changes coming within the meaning and equivalency range of the appended claims are intended to be embraced therein.
This application claims the benefit of U.S. Provisional Patent Application 61/023,939, filed Jan. 28, 2008, which is incorporated herein by reference.
Number | Name | Date | Kind |
---|---|---|---|
6363104 | Bottomley | Mar 2002 | B1 |
20020136234 | Eriksson et al. | Sep 2002 | A1 |
20060268962 | Cairns et al. | Nov 2006 | A1 |
20090268787 | Cairns et al. | Oct 2009 | A1 |
Number | Date | Country | |
---|---|---|---|
20090190642 A1 | Jul 2009 | US |
Number | Date | Country | |
---|---|---|---|
61023939 | Jan 2008 | US |