Automatic frequency correction method and apparatus for time division duplex modes of 3g wireless communications

CA 02461336 2006-07-31

[0001] AUTOMATIC FREQUENCY CORRECTION METHOD

AND APPARATUS FOR TIME DIVISION

DUPLEX MODES OF 3G WIRELESS COMMUNICATIONS

[0002] FIELD OF INVENTION

[0003] The present invention relates to the field of wireless communications. More specifically, the present invention relates to the field of third generation 3G wireless communications employing Time Division Duplex (TDD) and to frequency error detection and correction within a wireless communication system receiver.

[0004] BACKGROUND OF THE INVENTION

[0005] In typical wireless communications systems, a frequency difference between the transmitter and receiver local oscillators can prevent the transmission of data. Additionally, because many systems utilize the same Local Oscillator (LO) for both receiver and transmitter functions, a large frequency offset can cause significant out-of-band interference.

[0006] In order to overcome this problem, prior systems have utilized differential detection of phase or applied the Discrete Fourier Transform to estimate the frequency error and apply an update to the Local Oscillator. However, these prior systems either ignored the effects of multipath interference or combined the AFC with a RAKE receiver. Therefore, these prior techniques were not applicable for systems that use Multi-user Detection without a RAKE receiver.

[0007] SUMMARY OF THE INVENTION

[0008] The present invention enables detection and correction of the oscillator frequency error within a wireless communication system receiver. Moreover, the present invention provides robust performance in the presence of multipath interference. Furthermore, the present invention overcomes the interference problem as well as exploiting the diversity gains associated with large delay spread. Additionally the invention provides the capability of rejecting inter-cell and intra-cell interference sources, while operating effectively in the

-1CA 02461336 2004-03-23

WO 03/030400 PCT/US02/30802 presence of both RF carrier offset and sampling clock offset. Also, the present invention has adaptive tuning speed and functions with systems using Multi-user Detection algorithms without RAKE receivers and can operate on a discontinuous pilot (training) signal.

[0009] The present invention includes a frequency estimator having block correlators, conjugate product and sum block, accumulation block, multipath detection and a loop filter (adaptive bandwith). The multipath detection includes a search block, threshold detection block and a block for combining multipath components.

[0010] BRIEF DESCRIPTION OF THE DRAWINGS

[0011] The present invention will be understood from the following description and drawings in which like elements are designated by like numerals and, wherein:

[0012] Figure 1 is a block diagram of a automatic frequency control (AFC) algorithm employing the technique of the present invention.

[0013] Figure 2 is a block diagram illustrating the frequency estimation algorithm of the present invention.

[0014] Figures 3 and 4 are schematic diagrams which illustrate the structure included in each block correlator of Figure 2.

[0015] Figure 5 is a schematic diagram showing the conjugate product and sum block of Figure 2 in greater detail.

[0016] Figure 6 is a schematic diagram showing the details of the loop filter block of Figure 1.

[0017] Figure 7 is a flow diagram showing the algorithm performed by the apparatus of Figure 2.

[0018] DETAILED DESCRIPTION OF THE INVENTION

[0019] Figure 1 is a block diagram of the closed loop automatic frequency control (AFC) 10 wherein a received signal Rx is reduced to baseband at multiplier 12 by a voltage controlled oscillator (VCO) 14. The received baseband

CA 02461336 2004-03-23

WO 03/030400 PCT/US02/30802 signal Rx undergoes analog-to-digital conversion (ADC) at 16, automatic gain control (AGC) at 18 and then passes through a root-raised cosine (RRC) filter 20. [0020] After cell search at 22 and frequency estimation at 24, the frequency estimate is applied to loop filter 26. This digital output is converted by digital-toanalog converter (DAC) 28 to adjust the frequency of VCO 14, which is also used for transmission wherein baseband (BB) Tx data is converted by digital-to-analog converter (DAC) 30 which is used to modulate the carrier frequency provided by VCO 14 at multiplier 32.

[0021] Figure 2 is a block diagram showing the steps that are performed by frequency estimation block 24 and which shows the frequency estimation algorithm in greater detail.

[0022] Initially, the frequency estimation algorithm performs four (4) block correlations of the received signal samples with a known reference (midamble) at 24-1. The output of the four (4) block correlators, at 24-2 are successively multiplied in the conjugate sense to produce three (3) complex numbers with angles representative of the phase shift, in time, from one correlator to the next. These three conjugate products are then summed together to produce a lower variance estimate of the phase change. The output of the accumulate block 24-3 is a function of window lag, i, the value being accumulated over N frames. After N frames of data have been processed, the accumulated D(i) (values) are searched for the three values that have the largest absolute value, DO (largest), DI and D2, at 24-4. Magnitudes of these values are computed at 24-5 in order to obtain the three largest D(i) values.

[0023] A detection threshold is then applied at 24-6 which threshold is based on the magnitude of the peak value (DO). If the magnitude of the second and third largest components exceed this threshold they are deemed sufficiently large to be included in the frequency estimation computation.

[0024] After the threshold detection is performed, the surviving multipath components are then coherently summed at 24-7 to provide a single complex number whose angle may be used as an estimate of the phase change between correlator blocks. The frequency estimate is computed, employing blocks 24-8

CA 02461336 2004-03-23

WO 03/030400 PCT/US02/30802 and 24-9, which utilize two (2) approximations, to be described hereinbelow, avoiding the need for explicit trigonometric calculation.

[0025] Figure 3 shows the sliding window block correlation operation. Due to possible corruption of the first portion of the midamble with multipath interference from the first data burst, the last 456 chips of the midamble are utilized in the frequency estimation. The window that is searched includes 49 leading, 49 lagging and the 0 lag alignments, the total number of samples executed by the sliding window block correlator being 1108. In a 3GPP TDD communications system, equal length frames of 10 ms are comprised of fifteen (15) equal length time slots each having 2560 chips.

[0026] At each lag the four B chip (2B sample) correlations are performed, as shown in Fig. 3.

[0027] Fig. 4 shows the details of the first block correlator that produces Ro,i. As shown in Fig. 4, each received sample is correlated with a known midamble and summed with the next successive correlation.

[0028] Figure 5 shows the conjugate product and summation operation 24-2 which is performed on the outputs of the sliding block window correlators of block 24-1. The correlator outputs R are complex vectors representing the centroid of the received samples with the midamble modulation removed. The next step is estimating the phase change from one correlator to the next which is accomplished by computing a conjugate product of successive correlator outputs. Each output from a conjugate product operation is a complex vector whose angle approximates the phase change from the center of one correlation to the next. The three conjugate products developed by product circuits Pl, P2 and P3 are summed together at SI and S2 to produce a lower variance estimate of the phase change from one correlator to the next.

[0029] The D(i) values of the conjugate product and sum block 24-2 are accumulated over N midambles before computing a frequency estimate.

[0030] The accumulation time constant N is initialized to be 10 and is subsequently determined based on the most recent estimate of the absolute value of the frequency error. The value of N is selected to minimize the variance of the

CA 02461336 2004-03-23

WO 03/030400 PCT/US02/30802 frequency estimate while preventing significant drift during the estimation interval.

[0031] After N midambles have been processed through the sliding window correlators 24-1, conjugate product, sum 24-2 and accumulator 24-3, a search is performed to find the lag, i, which maximizes the magnitude of D (i). Due to the fact that there may be multiple resolvable multipath components, the three (3) largest paths are sought, the number of paths sought being a compromise between additional signal-to-noise ratio (SNR) improvement and increased hardware complexity.

[0032] Since it is possible that there is only one resolvable multipath component available, the second (DI) and third (D2) largest components are tested for significance. DI and D2 are considered significant if they are greater than half DO in the magnitude square sense. Thus DI and D2 are accepted it if they are greater than DO divided by the square root of 2 (DO/V2) and rejected if otherwise.

[0033] The multipath components meeting the above requirements are then combined into a single complex vector at 24-7, whose angle is an estimate of the phase change of the carrier offset over one block time.

[0034] In order to extract angle information from the multipath combiner output, the complex variable is scaled to unit magnitude and an approximation of the complex absolute value function is utilized, the approximation being that the imaginary part of the complex vector is equal to the argument of the complex vector which is equal to 0, if Θ is much less than 1 (0 «1) and the absolute value of the complex vector is 1.

[0035] This approximation simplifies the implementation of the algorithm, alleviating the need to perform trigonometric computations and it has been found that the error introduced by the approximations tends to zero as the AFC algorithm converges (0->l).

-5CA 02461336 2004-03-23

WO 03/030400 PCT/US02/30802

[0036] Loop filter 26 takes the estimated frequency error ε and performs an integration operation in order to obtain v(t) which is represented as: v(t)=v(t1)+λε(ί).

[0037] This is also depicted in Figure 6 wherein input ε is applied to an amplifier having a gain of -1 and summed with the previous value v(t-l) obtained at Dn, at summer S.

[0038] It should be noted that the integration is performed only when the errors is dumped from the previous block. Therefore, the value of v changes after N midambles are processed. A convergence detection algorithm (CDA) may be employed to determine convergence.

[0039] One technique is to compare the frequency estimate generated at the output of 24-9 against a threshold and if the estimated frequency error is smaller than |a|, convergence has been reached. The algorithm is considered memoryless because convergence is based only on the current estimate of frequency error.

[0040] An alternative arrangement is to declare convergence when two (2) successive frequency estimates are below a detection threshold a. Alternatively, the two frequency estimates need not be successive.

[0041] In still another alternative, convergence is detected based on a twopoint moving average of frequency estimates powering below a detection threshold a by successively averaging the last two frequency estimates obtained at 24-9 and comparing them against a threshold.

[0042] Regarding the detection threshold optimization employed by block

24-6, based on tests performed, an optimal choice of relative detection threshold is 0.56 (i.e. 0.56D(o)) which provides an improvement in the probabihty of p=0.99 in a convergence time of 0.65 seconds.

[0043] The optimum choice of loop gain λ is dependent upon the SNR and channel conditions. The optimal choice for loop gain is 0.26 which provides a significant improvement and success probability for the AWGN channel with an SNR of -3dB and two (2) active midambles.

-6CA 02461336 2004-03-23

WO 03/030400 PCT/US02/30802

[0044] In order to prevent the loss of coherency during the accumulation interval the relationship between N and estimated frequency error has been adjusted. The enhanced values prevent the drift of the clock from exceeding 0.25 chips over the accumulation periods, the value of N varies from 1 to 30 as a function of absolute frequency error of 6,000 to 0, the lower the absolute frequency error, the higher the number N of midambles accumulated.

[0045] Based on a comparison of the use of 456 versus 512 chips of the midamble in the correlation stage wherein the elimination of the first 56 chips of the received midamble which may be corrupted by multipath interference of the first data burst is offset by a reduction SNR of about 0.5dB, it has been determined that for all 3 WG4 test channels, the use of all 512 chips of the midamble is desirable. In a burst type 1, for example, each time slot has two (2) groups of data symbols, each having 976 chips separated by a 512 chip midamble and a 96 chip guard period following the last group of the two (2) groups of data symbols.

[0046] Previously, the window searched is described as including 49 leading, 49 lagging and no (i.e. 0) lag alignments. A more reasonable leading path search was determined to be 10 chips making the total number of samples required for the sliding window block correlator to execute as being 1142 samples, which reduction in window size is still acceptable for the worst case multipath WG4 channel model (case 2) in which the largest resolvable path is 46Tc delayed relative to the direct path.

[0047] An approach to multipath combining employed in blocks 24-6 and 24-7 in which the largest path DO is twice the weight of the second largest when only two paths survive, has been compared with a multipath combiner that treats the two surviving paths with equal gain and it has been found that equal gain combining resulted in slightly better performance in the WG4 case 1 and essentially the same performance for the other cases and hence is the preferred approach where only D(0) and D(l) are to be combined.

[0048] The present invention can also be implemented using an alternative method of estimating the phase difference (based on a plurality of multipath

-7CA 02461336 2004-03-23

WO 03/030400 PCT/US02/30802 components). In this case the phase estimation will still include a quality measure similar to the correlation magnitude used in the invention.

[0049] The same approach used to adjust the accumulation period (adaptive tuning rate) could be applied to the correlation block size. For large frequency offsets, smaller correlation block sizes are preferable because of the possibility of abasing and the loss of coherency in the estimate. As the frequency error diminishes, the correlation block size could be increased to improve the processing gain of the correlation and obtain more refined estimates of frequency error.

* * *