Adaptive channel estimation in a wireless communication system

[0001] The present Application for Patent is a Continuation and claims priority to U.S. Pat. No. 6,493,329 entitled “ADAPTIVE CHANNEL ESTIMATION IN A WIRELESS COMMUNICATION SYSTEM,” issued Dec. 10, 2002, assigned to the assignee hereof and hereby expressly incorporated by reference herein.

[0002] I. Field of the Invention

[0003] The present invention relates to wireless communication systems. More particularly, the present invention relates to a novel and improved method and apparatus for adaptively estimating the channel conditions of a wireless communication channel.

[0004] II. Description of the Related Art

[0005] In a wireless radiotelephone communication system, many users communicate over a wireless channel. Communication over the wireless channel can be one of a variety of multiple access techniques that allow a large number of users in a limited frequency spectrum. These multiple access techniques include time division multiple access (TDMA), frequency division multiple access (FDMA), and code division multiple access (CDMA).

[0006] The CDMA technique has many advantages. An exemplary CDMA system is described in U.S. Pat. No. 4,901,307, entitled “SPREAD SPECTRUM MULTIPLE ACCESS COMMUNICATION SYSTEM USING SATELLITE OR TERRESTRIAL REPEATERS,” issued Feb. 13, 1990, assigned to the assignee of the present invention, and incorporated herein by reference. An exemplary CDMA system is further described in U.S. Pat. No. 5,103,459, entitled “SYSTEM AND METHOD FOR GENERATING SIGNAL WAVEFORMS IN A CDMA CELLULAR TELEPHONE SYSTEM,” issued Apr. 7, 1992, assigned to the assignee of the present invention, and incorporated herein by reference.

[0007] In each of the above patents, the use of a forward-link (base station to mobile station) pilot signal is disclosed. In a typical CDMA wireless communication system, such as that described in EIA/TIA IS-95, the pilot signal is a “beacon” transmitting a constant zero symbol and spread with the same pseudonoise (PN) sequences used by the traffic bearing signals. The pilot signal is typically covered with the all-zero Walsh sequence. During initial system acquisition, the mobile station searches through PN offsets to locate a base station's pilot signal. Once it has acquired the pilot signal, it can then derive a stable phase and magnitude reference for coherent demodulation, such as that described in U.S. Pat. No. 5,764,687 entitled “MOBILE DEMODULATOR ARCHITECTURE FOR A SPREAD SPECTRUM MULTIPLE ACCESS COMMUNICATION SYSTEM,” issued Jun. 9, 1998, assigned to the assignee of the present invention, and incorporated herein by reference.

[0008] A functional block diagram of a typical prior art forward link data formatter as used by a CDMA base station is shown in FIG. 1. Data source 102 may be, for example, a variable rate vocoder such as that described in U.S. Pat. No. 5,657,420, entitled “VARIABLE RATE VOCODER,” issued Aug. 8, 1997, assigned to the assignee of the present invention and incorporated herein by reference. Data source 102 generates traffic channel information in the form of frames of digital data. CRC and tail bit generator 104 calculates and appends cyclic redundancy check (CRC) bits and tail bits to the frames generated by data source 102. The frame is then provided to encoder 106, which provides forward error correction coding, such as convolutional encoding, upon the frame as is known in the art. The encoded symbols are provided to repetition generator 120, which repeats the reordered symbols to provide the appropriate modulation symbol rate. The repeated symbols are then provided to interleaver 108, which re-orders the symbols in accordance with a predetermined interleaver format. The repeated, interleaved symbol stream is then covered with one of a set of orthogonal Walsh sequences in traffic Walsh coverer 122, and gain adjusted in gain element 124. It should be understood that other forward link data formatters are also known in the art. For example, it is well known that the repetition generator 120 may be placed after the interleaver 108.

[0009] Pilot signal generator 128 generates a pilot signal, which may be a sequence of all ones. The pilot signal is then covered with the all-one Walsh sequence and combined with the output of gain element 124 in combiner 136. The combined pilot channel and traffic channel data (which may be plus or minus ones) is then spread in PN spreader 138 using a complex PN code generated by PN generator 140, and then transmitted by radio frequency transmitter 142 over antenna 144. A similar forward link data formatter is disclosed in U.S. Pat. No. 6,396,804, entitled “HIGH DATA RATE CDMA WIRELESS COMMUNICATION SYSTEM,” issued May 28, 2002, and assigned to the assignee of the present invention and incorporated by reference herein.

[0010] Other data formatting techniques also exist. For example, in the cdma2000 reverse link, the pilot signal is time-multiplexed with power control commands. Additionally, in W-CDMA, the forward link uses dedicated pilot signals that are time-multiplexed with other information.

[0011] FIG. 2 illustrates a functional block diagram of a typical prior art data demodulator for use in a CDMA mobile station. Receiver (CVR) 202 receives and downconverts the signals transmitted by transmitter 142 of FIG. 1. The digital baseband output of receiver 202 is despread in PN despreader 204 using the complex PN code generated PN generator 206, which is the same complex PN code as that generated by PN generator 140 of FIG. 1.

[0012] The despread signal is then Walsh uncovered in traffic channel Walsh uncoverer 208 using the same Walsh sequence as that of the traffic channel Walsh coverer 122 of FIG. 1. The Walsh-uncovered chips are then accumulated into Walsh symbols in Walsh chip summer 210 and provided as a traffic channel signal to dot product circuit 212. In some applications, an additional delay element (not shown) is introduced between Walsh chip summer 210 and dot product circuit 212 to account for delays introduced by pilot filter 216. However, if pilot filter 216 is a causal filter, such a delay element (not shown) is not necessary. The dot product circuit is also known as a “conjugate product” circuit. It performs the operation expressed mathematically by one of the following equivalent forms: <a,b>=a•b=ab* , where b* is the complex conjugate of b.

[0013] The despread signal is also provided to Walsh chip summer 214 where they are accumulated into Walsh symbols and provided to pilot filter 216 as pilot channel symbols. Note that since the pilot channel is covered with the all-one Walsh sequence in Walsh coverer 134 of FIG. 1, a vacuous operation, the corresponding uncoverer is also vacuous in operation. However, in the general case, the pilot signal may be uncovered using any same Walsh sequence as is used to cover it. The pilot filter 216 serves to reject the noise in the pilot symbols, providing a phase and scale reference for the dot product circuit 212.

[0014] Once per symbol, the dot product circuit 212 computes the component of the traffic channel signal in phase with the pilot channel signal generated by the pilot filter 216. As described in U.S. Pat. No. 5,506,865, entitled “PILOT CARRIER DOT PRODUCT CIRCUIT,” issued Apr. 9, 1996, assigned to the assignee of the present invention and incorporated herein by reference, the dot product adjusts both the received signal's phase and scale as needed for coherent demodulation.

[0015] The symbols output from dot product circuit 212 are de-interleaved in de-interleaver 218, using the same format used by interleaver 108 of FIG. 1. The de-interleaved symbols are then decoded in decoder 220 according to the error correcting codes employed by encoder 106 of FIG. 1. The resulting decoded symbols are analyzed on a frame-by-frame basis by quality indicator (CRC CHECK) 222 to ensure that the frame was properly decoded. If the frame was properly decoded, then that decoded frame is forwarded for further processing. Quality indicator (CRC CHECK) 222 typically would examine the CRC portion of the frame, but may also use other frame quality indications such as Yamamoto metrics.

[0016] A typical pilot filter 216 is implemented as an equal-weight finite impulse response (FIR) filter with all defining parameters (e.g., weighting, window width, window center) remaining constant regardless of the channel conditions. Alternately, an exponential decay infinite impulse response (IIR) filter having fixed parameters (e.g., decay time constant, scaling) may be used. In other words, the designer of a typical prior art pilot filter 216 will choose a static filter design that performs adequately for a given energy per bit to noise density ratio (E_b/N₀) under most channel conditions of interest, but not optimally over the entire range of conditions.

[0017] As a mobile station (e.g., a cellular telephone, PCS telephone or other wireless remote communication terminal) moves through the terrestrial environment, the signals it transmits and receives will experience various types of fading. The mobile environment is usually characterized by fading that can be either Rician or Rayleigh in nature. Other types of fading are also possible. The fading characteristic in the typical channel signal is caused by the signal being reflected from many different features of the physical environment, thus it is called multipath fading. At the UHF frequency bands usually employed for mobile radio communications, including those of cellular mobile telephone systems, significant phase differences in signals traveling on different paths may occur. The possibility for both constructive and destructive summation of the signals may result, with on occasion deep fades occurring.

[0018] Multipath channel fading is a function very sensitive to the physical position of the mobile unit. A small change in position of the mobile unit changes the physical delays of all the signal propagation paths, which further results in a different phase for each path. Thus, the motion of the mobile unit through the environment can result in a rapid fading process. For example, in the 850 MHz cellular radio frequency band, this fading can typically be as fast as one fade per second for every mile per hour of vehicle speed. Fading this severe can be extremely disruptive to signals in the terrestrial channel resulting in poor communication quality, particularly as the speed of the mobile station increases beyond 150 km/hr.

[0019] As previously stated, the typical fixed-parameter pilot filter 216 is not optimized for such a broad range of channel conditions. It is typically designed to work adequately at speeds from stationary to about 120 km/hr, or about as fast as a mobile station might be expected to travel in a motor vehicle on the highway. However, since the fading characteristics of the channel are vastly different as between a slow-moving mobile station and a fast-moving mobile station, the typical fixed-parameter pilot filter 216 cannot be optimized for both extremes. Typically, this forces the designer to design a pilot filter 216 that works well only when the mobile station is stationary or moves slower than about 150 km/hr, and works poorly beyond 150 km/hr. As transportation such as bullet trains and airplanes exceeds this speed, it is unlikely that the user of a mobile station will be able to obtain reliable communications. Even when operational, the signal-to-noise ratio (or in other words, the E_b/N₀) of the communication link must be kept at a high enough level to be reliable in these severe fading conditions. Increasing the E_b/N₀ of the communication link decreases the total capacity of the wireless system, particularly in a CDMA system where one transmitter's transmissions comprise interference to all other transmitters in the same CDMA frequency band. As a result, the designer of a fixed-parameter pilot filter 216 generally adopts an unfavorable compromise in selecting the filter parameters when faced with such a broad range of channel conditions.

[0020] Thus, there is a need for a more optimal pilot filtering method and apparatus, particularly in the wireless communication environment, that avoids these shortcomings in the prior art.

[0021] The present invention is a novel and improved method and circuit for adaptively estimating channel conditions of a pilot channel in a wireless communication system. The method includes estimating channel statistics of the pilot channel, and adaptively filtering the pilot channel in response to the estimated channel statistics. The estimation is performed by filtering a channel signal derived from the pilot channel, or any channel bearing ambiguous data with the ambiguity removed after detection or decoding. This determines an estimated channel mean and an estimated channel covariance. In order to perform the adaptive filtering, the present invention partitions the pilot channel into one or more time slots and weights each time slot according to the estimated channel statistics. Thus, an advantage of the present invention is that it automatically and continually updates the pilot filter parameters in order to optimize the pilot filter performance over a broad range of channel conditions.

[0022] In one embodiment of the present invention, the channel signal is filtered in one or more infinite impulse response (IIR) filters to determine the channel statistics. In another embodiment, the channel signal is filtered in a combination of IIR and finite impulse response (FIR) filters to determine the channel statistics.

[0023] In an embodiment applicable to the cdma2000 system wherein the pilot channel comprises a pilot signal portion having a known sign and a power control bit portion having an unknown sign. The method of the present invention may also include determining the sign of the power control bit portion, correcting the sign of the power control bit portion, and then combining the sign-corrected power control bit portion with the pilot signal portion to generate the channel signal from which the channel statistics are estimated by the channel statistics estimator. In yet another embodiment, the method includes determining the sign of the traffic signal and generating the channel signal in response to the traffic signal.

[0024] In yet another embodiment applicable to a cdma2000 system, wherein the pilot channel comprises a sequence of power control groups, each power control group having a pilot signal portion having a known sign and a power control bit portion having an unknown sign, the time slots have a duration substantially equal to the duration of the pilot signal portion.

[0025] A circuit for performing the method of the present invention is also described herein.

[0031] The present invention will now be described in an application that is specific but not limited to the demodulation of a cdma2000 reverse link signal. As is described in the proposed IS-2000 standard, the cdma2000 reverse link signal comprises a traffic channal signal and a reverse link pilot signal. Punctured into the reverse link pilot channel (R-PICH) is a forward link power control subchannel that is a single bit that spans the last quarter of each 1.25 ms power control group (PCG). This forward link power control bit is a signal from the mobile station to the base station transmitter to either increase or decrease its transmit power based on the mobile station's reception quality on the forward link. Such a reverse link pilot is described in detail in U.S. Pat. No. 6,396,804, entitled “HIGH DATA RATE CDMA WIRELESS COMMUNICATION SYSTEM,” issued May 28, 2002, and assigned to the assignee of the present invention and incorporated by reference herein.

[0032] Because the R-PICH is used by the base station as a coherent reference for data demodulation, a frequency reference for frequency tracking, and a received power reference for power control measurements, the uncertainty introduced by the forward link power control subchannel in the otherwise determinate signal could degrade the reverse link performance. In other words, the puncturing of the forward link power control subchannel into the R-PICH could lead to poorer reverse link performance than a continuous pilot channel. For example, if we use only the non-punctured portion of the pilot channel to estimate the channel, the signal-to-noise ratio of the channel estimate calculated by the base station is decreased. A technique for resolving the sign ambiguity of the R-PICH and reconstructing a continuous pilot signal therefrom is given in U.S. Pat. No. 6,304,563, entitled “METHOD AND APPARATUS FOR PROCESSING A PUNCTURED PILOT CHANNEL,” issued Oct. 16, 2001, and assigned to the assignee of the present invention and incorporated here by reference.

[0033] It will be understood by a person of ordinary skill in the art that although the present invention is disclosed with reference to a cdma2000 system having a R-PICH with a punctured power control subchannel, the present invention is equally applicable to other wireless communication systems having a non-punctured pilot channel. Thus, the following figures are intended to be example applications of the present invention and are not intended to limit the present invention to the cdma2000 system.

[0034] Turning now to FIG. 3, a first embodiment of the present invention is illustrated as being implemented in a cdma2000 reverse link traffic channel demodulator. Channel statistic estimator 300 estimates channel statistics which are used by adaptive non-causal channel estimator 302 to adaptively filter the pilot channel as will be discussed in more detail below. It should be noted that in other embodiments, indeed in the embodiment of FIG. 4 described below, an adaptive causal channel estimator or filter may be used in place of adaptive non-causal channel estimator 302. However, in the embodiment of FIG. 3, a non-causal filter is used for channel estimation due to its improved performance over a causal filter when the delay of the non-causal filter is tolerable. In the exemplary embodiment of FIG. 3, adaptive non-causal channel estimator 302 is an FIR filter. However, in the general case, adaptive non-causal channel estimator 302 may be an IIR filter or a hybrid filter with both FIR and IIR characteristics. One simple example is a filter that outputs the difference of two exponential decay IIR filters of the same time constant but different weights, such that the effective impulse response of the filter is of finite extent. Another example would be a cascaded FIR and IIR filter.

[0035] The filtered pilot channel symbols are then input to conjugate product circuits 306 and 308. Conjugate product circuit 308 also receives as a second input, the despread, Walsh uncovered, and delayed traffic channel that has been delayed by delay circuit 310. Delay circuit 310 delays the traffic channel by an amount substantially equal to the delay introduced by adaptive non-causal channel estimator 302. Conjugate product circuit 308 produces traffic channel symbols (soft decisions) for further de-interleaving and decoding as is known in the art.

[0036] Conjugate product circuit 306 receives as a second input the delayed pilot channel that has been delayed by delay circuit 304. Delay circuit 304 delays the pilot channel by an amount substantially equal to the delay introduced by adaptive non-causal channel estimator 302. Conjugate product circuit 306 performs a conjugate multiplication of the estimator 302 output with the delayed pilot channel which contains a power control bit punctured into the reverse link pilot channel. Conjugate product circuit then forwards the resultant power control bit signal to past forward link (FL) power control bit (PCB) detector 314.

[0037] Past FL PCB detector 314 compares the resultant power control bit signal to a threshold and thereby determines the sign of the punctured FL PCB. It should be noted that the term “past” is used in describing the operation performed by past FL PCB detector 314 because the power control bit being detected at any given time is the power control bit from a pas power control group, as delayed by delay 304 and adaptive non-causal channel estimator 302. The FL PCB decision determined by past FL PCB detector 314 is then forwarded to PCB sign corrector 316 where the sign of the PCB is flipped, if necessary, to match the sign of the remainder of the pilot channel (+1). The output of PCB sign corrector 316 is then provided, along with the further delayed pilot channel from delay 312, to time multiplexer (MUX) 318. Time MUX 318 provides, during the first ¾ of the power control group, an output equal to the further delayed pilot channel from delay 312, which compensates for the processing time of the past FL PCB detector 314. Time MUX 318 provides, during the remaining ¼ of the power control group, the output of PCB sign corrector 316 which is now equal in sign (+1) to the pilot channel. The resulting output signal from time MUX 318, is, therefore, a reconstructed pilot channnel signal of constant sign. A similar technique for continuous pilot channel reconstruction in given in the above-incorporated U.S. Pat. No. 6,304,563.

[0038] The reconstructed channel signal is provided to channel statistic estimator 300 which estimates channel statistics therefrom for use in setting the weighting factors used by adaptive non-causal channel estimator 302 as mentioned above. The generation of channel statistics by channel statistic estimator 300 and the application of weighting factors by adaptive non-causal channel estimator 302 will be discussed in further detail below.

[0039] It should be noted that the elements in the left-hand side of FIG. 3 that are involved soley in generating the reconstructed channel signal for use by channel statistic estimator 300 (i.e., delay 304, conjugate product circuit 306, delay 312, past FL PCB detector 314, PCB sign corrector 316 and time MUX 318) are necessary only for a punctured pilot channel such as the R-PICH of the cdma2000 reverse link. These elements may be omitted in applications having a non-punctured or otherwise continuous pilot channel. In such a case, the pilot channel signal itself may be sufficient for use by channel statistic estimator 300. Furthermore, the channel signal may also be generated from the additional energy contained in the traffic signal using the method described in U.S. Pat. No. 6,452,917, entitled “CHANNEL ESTIMATION IN A CDMA WIRELESS COMMUNICATION SYSTEM,” issued Sep. 17, 2002, assigned to the assignee of the present invention and incorporated herein by reference. In the just-mentioned patent application, the decoded traffic channel symbols are re-encoded, re-interleaved, weighted according to the frame quality indicator and combined with the pilot channel estimate, which in FIG. 3, would be the reconstructed channel signal output from time MUX 318. Thus, in the general case, the generation of the channel signal, i.e., a signal that represents the channel, may be a weighted combination of contributions from the pilot portion of the R-PICH, the sign-corrected FL PCB from the R-PICH, and the sign-corrected traffic from one or more traffic channels.

[0040] Turning now to FIG. 4, a second embodiment of the present invention illustrated as being implemented in a forward link power control bit demodulator. In FIG. 4, channel statistic estimator 400, which may be the same as channel statistic estimator 300 of FIG. 3, is illustrated as receiving the channel signal, which may be the reconstructed channel signal of FIG. 3 in the case of a the R-PICH demodulator of FIG. 3, or a reconstructed channel signal from a sign-corrected traffic channel after successful decoding, re-encoding, and re-interleavi g as detailed in the above-referenced U.S. Pat. No. 6.452,917, or a weighted combination of the two signals. In the general case, channel statistici estimator 400 may be different than channel statistics estimator 300. Both may use essentially the same algorithm, described in detail below, but they may operate on different data, or over different time slots, or even different channel signal inputs. The channel signal is operated on by channel statistic estimator 400 as will be described further below in order to estimate the channel statistics by which adaptive predictive causal channel estimator 402 will assign filter weighting coefficients. In the exemplary embodiment of FIG. 4, adaptive predictive causal channel estimator 402 is an FIR filter. However, in the general case, it may be an IIR filter, or a hybrid filter having both FIR and IIR characteristics.

[0041] Adaptive predictive causal channel estimator 402 provides a channel estimate through delay 404 to conjugate product circuit 406. Conjugate product circuit 406 performs a conjugate multiplication on the channel estimate and the pilot channel (in this case, the R-PICH) to demodulate the FL PCB punctured into the R-PICH. The resultant power control bit signal is then forwarded to current FL PCB detector 408 for determination of the sign of the current FL PCB. Current FL PCB detector 408 may be of similar construction to past FL PCB detector 314 of FIG. 3. Current FL PCB detector 408 estimates a current FL PCB decision regarding the sign of the FL PCB by comparing the power control bit signal from conjugate product circuit 406 to a threshold. The current FL PCB decision may then be used by the forward link transmitter (not shown) to increase or decrease its power as necessary.

[0042] As described above, in both the embodiments of FIG. 3 and FIG. 4, the channel statistic estimators 300 and 400 estimate channel statistics for use by adaptive non-causal channel estimator 302 and adaptive predictive causal channel estimator 402, respectively, in setting their respective filter coefficients. To explain how this process occurs, the following mathematical description will be introduced. First, consider the received channel signal of a particular Rake finger of a Rake receiver entering the mobile station channel statistic estimator 300 or 400 input. This received signal is of the form:
y[n]=a[n]+w[n]  (1)
where

• • [0043] y[n] is a column vector representing the actual received channel signal at time slot n;
  • [0044] a[n] is a column vector representing the actual information signal as a true representation of the channel at time slot n;
  • [0045] w[n] is a column vector representing the actual noise contained in the channel signal at time slot n; and
  • [0046] n is a vector index in the time domain, and therefore may represent an arbitrary time slot for sampling the received channel signal.

[0047] Therefore, the dimension of the column vectors a[n], w[n], y[n], or the number of integral values that n can take on, is the same as the number of time slots of received channel signal used as input for adaptive channel estimation. In the second exemplary embodiment applicable to the predictive estimation of a FL PCB in a cdma2000 R-PICH, each time slot may be ¾ of a single 1.25 ms power control group (i.e., the non-punctured portion of the pilot channel). Thus, if we use 2 PCG worth of channel signal as input for estimation, as n ranges from 1 and 2, y[1] represents the received channel signal of the current ¾ PCG of pilot whereas y[2] represents the received channel signal of the previous ¾ PCG of pilot. These y[1] and y[2] values can be obtained by simple averaging of pilot symbols over the said ¾ PCG of pilots.

[0048] Other embodiments may use larger or smaller time slots depending on the resolution and estimation accuracy desired as a tradeoff to the costs of computation. For example, in the cdma2000 embodiment, the punctured FL PCB may be one coherent time slot, with an ambiguous sign, or it can comprise several smaller coherent time slots if finer resolution is desired. In the general case, the time slots may be of unequal duration, though the received channel signal at each time slot could be weighted appropriately to normalize the information signal amplitude or, depending on the application, the signal-to-noise ratio. In other embodiments, the time slots may be a configurable parameter. If a partition of the channel signal into finer time slots is used, the index n can range from one to a higher number. If a non-causal estimator is desired, as in the embodiment of FIG. 3, some of the elements of the column vectors y[n], a[n], and w[n] may also represent values after the time of interest in the system where the channel estimate is desired.

[0049] We must estimate the actual channel value x at a given moment, as shown below in Equation (2). Equation (2) represents the operation performed by both adaptive non-causal channel estimator 302 and adaptive predictive causal channel estimator 402. One difference between the two applications is the different partitions of the channel signal into different time slots as discussed above. Note that all estimates herein are denoted by the ^ (“hat”) notation Let x be the desired channel value at a time slot of interest. For example, while the index n for the y, a, and w vectors takes on the values 1 and 2 representing the pilot portions of the pilot channel during the current and previous PCGs as in the exemplary embodiment of FIG. 4, x represents the channel value at the current PCB portion of the pilot channel,
{circumflex over (x)}={circumflex over (m)}_x+H(y[n]−{circumflex over (m)}_y[n])  (2)
where

• • [0050] {circumflex over (x)} is the estimated channel value at the time slot corresponding to the PCB portion of the pilot channel,
  • [0051] {circumflex over (m)}_x is the estimated mean of x;
  • [0052] {circumflex over (m)}_y[n] is the estimated mean of the channel signal at time slot n; and
  • [0053] H is a filter vector described by H={circumflex over (K)}_xy({circumflex over (K)}_yy)⁻¹
    where:
  • [0054] {circumflex over (K)}_xy is the estimated covariance matrix of the x and y[n]; and
  • [0055] {circumflex over (K)}_yy is the estimated covariance matrix of the vector y[n], and the superscript −1 represents the matrix inversion operation.

[0056] In other words, H is a vector representation of a linear-filter that when the filter is applied to an arbitrary signal z[n], the output resulting from the operation can be described in the following in time domain notation: $\begin{array}{cc}H\left(z\left[k\right]\right)=\sum _{k}h\left[k\right]z\left[k\right]& \left(3\right)\end{array}$
where k is the vector index.

[0057] It is also possible to make x into a vector by packing more values into it when the channel estimate at multiple time instances are desired, especially if this provides extra implementational convenience in some embodiment of application due to the computational architecture used. Thus:
{circumflex over (x)}[n]={circumflex over (m)}_x[n]+H(y[n]−{circumflex over (m)}_y[n])  (2a)
$\begin{array}{cc}\hat{x}\left[n\right]={\hat{m}}_x\left[n\right]+H\left(y\left[n\right]-{\hat{m}}_y\left[n\right]\right)& \left(2a\right)\\ H\left(z\left[k\right]\right)=\sum _{k}h\left[j,k\right]z\left[k\right]& \left(3a\right)\end{array}$
where H now represents the filter matrix, and j and k are the vector indices for the row and column, respectively, with j being the index of a particular element {circumflex over (x)}[j] in the vector x to be estimated.

[0058] The values {circumflex over (m)}_y[n], {circumflex over (K)}_xy, and {circumflex over (K)}_yy are what are referred to herein as the estimated channel statistics. These are the values estimated by channel statistics estimator 300 and 400 for use by adaptive non-causal channel estimator 302 and adaptive predictive causal channel estimator 402, respectively, in determining the filter weighting used to estimate {circumflex over (x)}. Thus, these channel statistics are the needed components of Equation (2) above, given that y[n], the actual received channel signal, is already known.

[0059] Next, the generation of these channel statistics is described. Channel statistic estimators 300 and 400 generate the estimated channel signal mean, {circumflex over (m)}_y[n], according to the following relationship:
{circumflex over (m)}_y[n]=g₁*y[n]  (4)
where g₁ is a filter impulse response, which in the general case may be either an FIR or an IIR or hybrid filter that has both FIR and IIR characteristics; and the operator * represents convolution in the time domain.

[0060] Note that y[n] now represents a signal vector in time, and the time index over which the convolution operation is applied is not shown. If it were shown, both {circumflex over (m)}_y[n] and y[n] would have two indices, with the first one (shown) representing the time slots (e.g. “current PCG” and “previous PCG” as used for the second exemplary embodiment of the application above), and the second (not shown) representing the updating of the vector values as what is “current time” proceeds.

[0061] The same method is used to calculate the estimate {circumflex over (m)}_x[n] since the elements of {circumflex over (m)}_x[n] represent the estimated channel just as the elements of {circumflex over (m)}_y[n] do, though at perhaps different time slots, so that only a shuffling of time slot indices is required before the direct application of the same method. Note that the channel mean usually does not change often, and when the time slots used as estimation inputs are close together in time, we may use one single value to represent all the elements in {circumflex over (m)}_x[n], {circumflex over (m)}_y[n] because they are similar in value, thus possibly saving computation. In another embodiment of application of the present invention in which the channel is known to have zero-mean, {circumflex over (m)}_x[n], {circumflex over (m)}_y[n] can be taken to be zero, thus simplifying Equation (2).
{circumflex over (x)}[n]=H(y[n])  (2b)

[0062] Channel statistic estimators 300 and 400 generate {circumflex over (K)}_xy, the estimated covariance of x and the vector y[n] according to the following relationship:
{circumflex over (K)}_xy[m,n]=g₂*(x[m]y*[n])  (5)
where

• • [0063] g₂ is a filter impulse response, which in the general case may be either an FIR or an IIR or hybrid filter that has both FIR and IIR characteristics;
    • [0064] m and n are the row and column vector indices, respectively; and
      y*[n] is the conjugate of y[n].
      Again, the time index over which the convolution operation is applied is not shown.

[0065] Channel statistic estimators 300 and 400 generate {circumflex over (K)}_yy, the estimated covariance of the y[n] vector according to the following relationship:
{circumflex over (K)}_yy[n,m]=g₃*(y[n]y*[m])  (6)
where g₃ is a filter impulse response, which in the general case may be either an FIR or an IIR or hybrid filter that has both FIR and IIR characteristics; and y*[m] is the conjugate of y[m].

[0066] Again, the time index over which the convolution operation is applied is not shown. The present invention also contemplates taking advantage of the fact that {circumflex over (K)}_yy may be conjugate symmetric, that is,
{circumflex over (K)}_yy[m,n]={circumflex over (K)}_yy*[n,m]  (7)
such that only half of the values in the matrix {circumflex over (K)}_yy need to be computed.

[0067] In addition, if part of the received channel signal vector y[n] has time slots of equal durations representing values equally spaced in time, we may also take advantage of the shift-invariant property of the covariance function and save computation further,
{circumflex over (K)}_yy[m,n]={circumflex over (K)}_yy[m+k,n+k]  (8)
where k is any integer value as long as m+k and n+k are within the scope of valid indices of the matrix {circumflex over (K)}_yy and where the above assumption about the partition of time slots is true.

[0068] The time constants of these three filters, g₁, g₂, g₃, are usually picked for the particular application at hand depending on the system parameters. For example, if the channel statistics do not change over a duration of time on the order of one second, the time constants of the filters g₁, g₂, g₃, may be chosen to be one second or shorter.

[0069] When multiple Rake fingers are used, a more sophisticated application of the present invention can make sure that the channel statistics from the different fingers are the same or similar, and use all the channel signals y[n] from different fingers for the statistics estimation, giving more accurate estimates due to the increased amount of available input for estimating the same or similar values of channel statistics.

[0070] So, as can be seen from the above mathematical relationships, channel statistic estimators 300 and 400 generate the channel statistics {circumflex over (m)}_x, {circumflex over (m)}_y, {circumflex over (K)}_xy, and {circumflex over (K)}_yy. These are in turn provided to the respective channel estimators (both adaptive non-causal channel estimator 302 and adaptive predictive causal channel estimator 402, which may also be referred to generically as “pilot filters”). The channel estimators of FIG. 3 and FIG. 4 then use these channel statistics in carrying out the operation described by Equation (2) above, i.e., adaptively estimating the channel conditions.

[0071] Central to the present invention is the concept that the channel statistics are slowly but constantly changing as the mobile station moves through the fading environment. These channel statistics are continually updated by channel statistic estimators 300 and 400 for use by adaptive non-causal channel estimator 302 and adaptive predictive causal channel estimator 402, respectively. Because these channel statistics form the basis for Equation (2), the adaptive filtering operation, both adaptive non-causal channel estimator 302 and adaptive predictive causal channel estimator 402 are changing their filter parameters in response to changes in channel conditions.

[0072] Advantageous to the present invention, Equation (2) gives an estimator that minimizes the mean squared error of the estimate output. Moreover, the present invention provides convenient and efficient methods to automatically calculate all the ingredients for employing this equation. Thus, both adaptive non-causal channel estimator 302 and adaptive predictive causal channel estimator 402 are continually updated by channel statistic estimators 300 and 400, respectively, such that they are optimized under the then-prevailing channel conditions. This allows both adaptive non-causal channel estimator 302 and adaptive predictive causal channel estimator 402 to generate the best channel estimate for a given channel condition under a broad range of channel conditions. In sharp contrast to the prior art, the present invention provides a method and apparatus for adaptively estimating the channel conditions in a wireless communication system based on the use of estimated channel statistics to determine the channel estimator's filter parameters. Thus, the present invention is always optimized for the current channel conditions, resulting in a significantly lower communication link E_b/N₀ requirement for a given bit error rate.

[0073] The previous description of the preferred embodiments is provided to enable any person skilled in the art to make or use the present invention. The various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without the use of the inventive faculty. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.