INITIALIZATION/PREWINDOWING REMOVAL POSTPROCESSING FOR FAST RLS FILTER ADAPTATION

Abstract: The present invention, generally speaking, accelerates convergence of a fast RLS adaptation algorithm by, following processing of a burst of data, performing postprocessing to remove the effects of prewindowing, fictitious data initialization, or both. This postprocessing is part of a burst mode adaptation strategy in which data (signals) get processed in chunks (bursts). Such a burst mode processing approach is applicable whenever the continuous adaptation of the filter is not possible (algorithmic complexity too high to run in real time) or not required (optimal filter setting varies only slowly with time). Postprocessing consists of a series of "downdating" operations (as opposed to updating) that in effect advance the beginning point of the data window. The beginning point is advanced beyond fictitious data used for initialization and beyond a prewindowing region. In other variations, downdating is applied to data within a prewindowing region only. The forgetting factor of conventional algorithms can be eliminated entirely. Performance equivalent to that of GWC RLS algorithms is achieved at substantially lower computational cost. In particular, a postprocessing Fast Kalman Algorithm in effect transforms an initialized/prewindowed least squares estimate into a Covariance Window least squares estimate. Various further refinements are possible. Initialization may be cancelled completely or only partially. For example, in order to reduce the dynamic range of algorithmic quantities, it may be advantageous to, in a subsequent initialization, add an increment to a forward error energy quantity calculated during a previous burst. Postprocessing may then be performed to cancel only the added increment. Also, to reduce the usual large startup error transient, the desired response data can be modified in a way that dampens the error transient. The modified desired response data are saved for use in later postprocessing. Furthermore, to allow for more rapid adaptation without the use of an exponential forgetting factor, a weighting factor less than one may be applied to the forward error energy quantity during initialization from one burst to the next. This allows for the most efficient use of data but limited adaptation within a burst, but more rapid adaptation from one burst to the next.

Description: INITIALIZATION/PREWINDOWING REMOVAL POSTPROCESSING FOR FAST RLS FILTER ADAPTATION BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to adaptive filters and more particularly to fast recursive least squares (fast RLS) adaptive filters. 2. State of the Art Adaptive filters find widespread use in communications. A known class of adaptive filter is the fast RLS adaptive filter. This class of adaptive filter, or more accurately, this class of filter adaptation algorithm, is recognized for its properties of fast convergence and reduced computational complexity compared to earlier RLS algorithms (although the computational complexity is still very considerable). In one aspect, the attractiveness of RLS algorithms lies in the ability to compute updated filter settings using a new data input together with the old filter settings. This "recursion" is performed differently in different variations of the algorithm. In one variation, designated by the term "growing window covariance" (GWC), the adaptation algorithm takes into account a growing window of data that begins at time zero, for example, and is incrementally extended at each sample time until the adaptation has concluded. In another variation, designated by the term "sliding window covariance" (SWC), the adaptation algorithm takes into account a window of data that at a particular point in the algorithm becomes fixed in size. In the SWC algorithm, once the size of the data window has been fixed, as updates are computed, they incorporate knowledge from a new sample and "disincorporate" knowledge from an oldest sample. In general, various different methods of windowing the input data of an adaptive filter are known. In order to achieve a favorable mathematical structure for computational efficiency, the most common implementations of the fast RLS algorithm use prewindowing. The prewindowing method makes the assumption that the input data prior to time zero are zero. This assumption results in a "prewindowing transient" during which the error quantity is unusually large. In other words, prewindowing perturbs the algorithm, delaying convergence. The GWC version of RLS does not require prewindowing. Without prewindowing, however, the RLS algorithm becomes computationally more burdensome and more implementation-sensitive. Initialization also perturbs the algorithm. As described in J.M. Cioffi and T. Kailath, "Fast, recursive least squares transversal filters for adaptive filtering", IEEE Trans, on ASSP, ASSP-32(2):304-337, April 1984, more rapid convergence may be obtained by initializing the input data matrix with a sparse "fictitious" data submatrix located within a region corresponding to negative time. A sample covariance matrix is then initialized so as to be in agreement with the fictitious data. A desired response data vector (or matrix) may also be initialized with some fictitious data to reflect prior information on the filter settings. These initialization techniques may be referred to as "soft constraint initialization". As compared to initializing the same quantities with zeros, the foregoing technique may significantly reduce the amount of perturbation experienced by the algorithm if the prior information is correct. Substantial perturbation remains, however, if the assumed prior information is incorrect, as is usually the case. In order to overcome the deleterious effects of prewindowing and initialization, an exponential "forgetting factor" is commonly used. As more and more data is processed, the influence of old data becomes exponentially more attenuated. Besides forgetting "incorrect" data, however, the same forgetting also gets applied to actual data, with the result that the total amount of data available is used less than efficiently. If this forgetting factor could be eliminated entirely, the result would be better estimation performance and fewer numerical problems. Further background information concerning RLS adaptation algorithms may be found in the following references, incorporated herein by reference: T. Kailath, Lectures on Wiener and Kalman Filtering, Springer- Verlag, Wien -- New York, 1981. S. Haykin, Adaptive Filter Theory, Prentice-Hall, Englewood Cliffs, NJ, 1995, third edition. B. Widrow and S.D. Stearns, Adaptive Signal Processing, Prentice-Hall, Englewood Cliffs, NJ, 1985. B. Widrow and E. Walach, "On the Statistical Efficiency of the LMS Algorithm with Nonstationary Inputs", IEEE Trans, on Information Theory, IT-30(2):211-221, March 1984, Special Issue on Adaptive Filtering. D.T.M. Slock, "On the Convergence Behavior of the LMS and the Normalized LMS Algorithms", IEEE Trans, on Signal Processing, 41(9):2811-2825, Sept. 1993. S.M. Kay, Fundamentals of Statistical Signal Processing: Estimation Theory, Prentice Hall 1993. E. Eleftheriou and D. Falconer, "Tracking Properties and Steady-State Performance of RLS Adaptive Filter Algorithms", IEEE Trans. ASSP, ASSP-34(5): 1097- 1110, Oct. 1986. D.D. Falconer and L. Ljung, "Application of Fast Kalman Estimation to Adaptive Equalization", IEEE Trans. Com., COM-26(10): 1439-1446, Oct. 1978. J.M. Cioffi and T. Kailath, "Fast, recursive least squares transversal filters for adaptive filtering", IEEE Trans, on ASSP, ASSP-32(2):304-337, April 1984. D.T.M. Slock and T. Kailath, "Numerically Stable Fast Transversal Filters for Recursive Least Squares Adaptive Filtering", IEEE Trans. Signal Proc, ASSP-39(1):92-114, Jan. 1991. D.T.M. Slock, "Backward Consistency Concept and Round-Off Error Propagation Dynamics in Recursive Least Squares Algorithms", Optical Engineering, 31(6): 1153-1169, June 1992. J.M. Cioffi and T. Kailath, "Windowed Fast Transversal Filters Adaptive Algorithms with Normalization", IEEE Trans, on ASSP, ASSP-33(3):607-625, June 1985. J.M. Cioffi, "The Block-Processing FTF Adaptive Algorithm", IEEE Trans, on ASSP, ASSP 34(l):77-90, Feb. 1986. SUMMARY OF THE INVENTION The present invention, generally speaking, accelerates convergence of a fast RLS adaptation algorithm by, following processing of a burst of data, performing postprocessing to remove the effects of prewindowing, fictitious data initialization, or both. This postprocessing is part of a burst mode adaptation strategy in which data (signals) get processed in chunks (bursts). Such a burst mode processing approach is applicable whenever the continuous adaptation of the filter is not possible (algorithmic complexity too high to run in real time) or not required (optimal filter setting varies only slowly with time). Postprocessing consists of a series of "downdating" operations (as opposed to updating) that in effect advance the beginning point of the data window. The beginning point is advanced beyond fictitious data used for initialization and beyond a prewindowing region. In other variations, downdating is applied to data within a prewindowing region only. The forgetting factor of conventional algorithms can be eliminated entirely. Performance equivalent to that of GWC RLS algorithms is achieved at substantially lower computational cost. In particular, a postprocessing Fast Kalman Algorithm in effect transforms an initialized/prewindowed least squares estimate into a Covariance Window least squares estimate. Various further refinements are possible. Initialization may be cancelled completely or only partially. For example, in order to reduce the dynamic range of algorithmic quantities, it may be advantageous to, in a subsequent initialization, add an increment to a forward error energy quantity calculated during a previous burst. Postprocessing may then be performed to cancel only the added increment. Also, to reduce the usual large startup error transient, the desired response data can be modified in a way that dampens the error transient. The modified desired response data are saved for use in later postprocessing. Furthermore, to allow for more rapid adaptation without the use of an exponential forgetting factor, a weighting factor less than one may be applied to the forward error energy quantity during initialization from one burst to the next. This allows for the most efficient use of data but limited adaptation within a burst, but more rapid adaptation from one burst to the next. BRIEF DESCRIPTION OF THE DRAWING The present invention may be further understood from the following description in conjunction with the appended drawing. In the drawing: Figure 1 is a block diagram of an equalizer with which the present invention may be used; and Figure 2 is a diagram illustrating initialization of the adaptive filters of the equalizer of Figure 1. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS In the following description, details are described of an advantageous postprocessing algorithm that may be used to remove the effects of prewindowing and/or initialization and hence accelerate convergence of a fast RLS adaptation algorithm. For purposes of illustration only, the example of a fractionally-spaced (FS) decision feedback equalizer (DFE) is considered. Such an FS DFE finds particular application in xDSL modems; however, the technique described is broadly applicable to filter adaptation in all fields of application. In a concrete embodiment, an example is given of a fast RLS algorithm that can be used in conjunction with the postprocessing procedure. Many of the algorithmic details of the fast RLS algorithm and the postprocessing procedure are themselves familiar from the published literature. Different internal variables and recursive relations can be used to implement the same strategies. To explain the initialization removal operation, the interpretation of initialization as the result of processing fictitious data is explained in detail. Fractionally-Spaced Decision Feedback Equalization Consider the channel equahzation application depicted in Fig 1 Bits get mapped mto (scalar valued) symbols a k belonging to some finite alphabet (e g {-1, +1}), which then modulate (linearly) a certain pulse shape The resulting signal gets transmitted over a channel that is assumed to be linear and time-invariant and hence can be represented by a certain channel impulse response Noise and interference get added dunng the transmission At the receivei the signal gets lowpass filtered and sampled Assume the sampling rate to be twice the symbol rate Then the even and odd received samples can be considered as two symbol rate sequences which together constitute one vector valued symbol rate sequence y(hX kT) Vk y(t 0 + {k + _)T) where T is the symbol period and the tune offset to determines the phase of the sampling operation The whole svstem taking the symbols α*. as input and producing the 2 x 1 vector samples Ϊ/J. as output is a discrete-time system operating at the symbol rate The 2 x 1 transfer function ll(z) indicated m Figure 1 represents the even and odd samples of the convolution of a transmission filter (pulse shape), channel transfer function and receiver filter Consider now a Decision- Feedback Equalizer (DFE) with feedforward filter F(z) (1 x 2) and feedback filter l-B(z) In Figure 1, F(z) is assumed to be potentially non-causal and hence leads the DFE to produce undelayed symbol estimates a k However in practice and m the equations below, the filter T(z) is assumed to be causal so that the DFL produces symbol estimates k - d with a certain delay d Then the DFE output becomes ojfc_ d = (1 - B(z)) a k _ d + F{z) y k f(ι) y k -ι where the symbol estimate is delayed with respect to the received signal An equahzation error a _ d can be formed by subtracting the equalizer output a k - d (symbol estimate) from the true symbol a k _ d I he equalization error can be written as a k _ d = a k _ d - a k - d = \i(z) a k _i - F{z) y k = ∑ b(ι) a k _ d→ - ∑ f(ι) y k _ t = B A k _ d - F Y k , 5(0) = 1 i-l) i 0 where the b(τ) are scalars and the /(.) are 1 X 2 vectors The true symbols <_ _ < _ aie available during the so called training period of the equalizer After the training period, when processing actual data the true symbol a k _ d can be replaced by a decision a k _ d based on the equalizer output a k -. d '1 he decision is obtained (in the absence of channel coding) by passing the equalizer output through a sheer, as indicated in Fig 1 The sheer takes the element from the symbol constellation that is closest to the equalizer output Unless a decision error occurs, ά k _ d — a k - The equalization error is used by adaptive filtering algorithms to adapt the equa zei filter settings The adaptation of the equalizer with actual data through the use of a k→i is called the decision directed mode Define the following d k - a - 6jc = a k - -\ where d denotes the desired-response signal for the equalizer and s k is the input signal to the feed back filter Note that s^ - d _ \ th input to the feedback filter consists of past svmbols/decisions With the r/2 feedforward filter spacing considered here there results, m fast RLS terminologv, 3 channels (symbol rate filter input streams) past symbols (s k ) and even and odd received samples (the two components of y ) The filtering (equalization) error can be rewritten as N e k (W) = d k - W U k = d k - ∑ w(ι) U k (t) 2= 1 where W = [-6(1) -b(n B ) /(0) and super script τ denotes matrix transposition The input data vector U k is a column vector, whereas filter coefficient vectors such as W will be denoted by row vectors This allows filter outputs to be written as )/\ U k without tranposes N = n # + 2τip = N/ι > f- NJJ is the total filter order (total number of filter coefficients) Prewindowed Least-Squares Data Structure To understand the origins of fast RLS algorithms, it is essential to pay particular attention to the structure of the data in the least-squares criterion used to determine the filter settings The structure of the data for the prewindowed least-squares criterion applicable to the DFE filter coefficients is shown m Fig 2 The top line indicates (discrete) time, progressing to the left, with variable tun. index i Assume that the LS problem is formulated with the data available at time i — k In the prewindowing formulation, it is assumed that the signals are zero before tune ? = 0 So the signals are assumed to start at time % - 0 For the equalizer application, if the filter adaptation is to start as soon as training data are available then the time origin should be chosen such that υ = a_ d _ 1 is the first available training sequence symbol The second line in Fig 2 contains the desired response signal in the row vector d k The structure of d k reflects the prewindowing Also indicated in d is an exponential forgetting/weighting factor λ € (0, 1] Actually, the data are weighted by so that the squared errors m the least squares cπteπon to be considered below will be weighted by = λ The exponential weighting reduces the influence of past data exponentially fast so that the LS filter solution will be reflecting the characteristics of the most recent data, in case these characteristics vary with time A corresponding time constant ^ is implied by the exponential weighting The matrix Ui^ -i k m Fig 2 is the extended input data matrix f/jv+3 k is a structured matrix For the DFE problem at hand, it consists of two structured submatπces The top submatπx is filled with the signal s, and has a Hankel structure (if λ = 1) elements along any anti-diagonal are constant This Hankel structure is perturbed by the exponential weighting, which weighs consecutive columns more and more The bottom submatπx of (/ JV - J . ls block Hankel (if A = 1) and is filled with the 2x 1 vector signal (blocks) y l The data matrix U N < 3 k is called extended because it contains nβ + 1 rows of signal s, and rip+l (block) rows of signal y x so it has m fact n β +l-|-2(n / . +1) = N+3 rows (one extra row per input channel) The extended data matrix will be required for the fast RLS algorithms If the bottom row with signal s k _ n _ and the bottom (block) row with signal y k _ nF are removed from U^+s k , then there results the regular input data matrix U^ k — U k with N rows The columns of the matrix U are the weighted input data vectors \ ~ U t which upon multiplication with the filter coefficient vector W gives the weighted DFE output \X W U % at time 1 The n that is indicated in Fig 2 denotes n — max{nβ, nj } (though the figure is drawn assuming that TI B = rip = n) The extended data vectors f/jv+. 1 Wltn tιme m tne ran e 1 6 [0,n-l] contain some samples of signal 6, and/or y t situated at negative time Those samples are not available if sampling occurs beginning from time . = 0 onwards In the prewindowing assumption those samples are set equal to zero This means that the data vectors [/ / v+ 3 , contain zeros over the indicated time range, which can be called the prewindowing transient period Hence within this peπod of n discrete time instants, the extended data vectors /^ +3 , do not fully contain genuine data In order to make an RLS algorithm start from well-defined quantities some initialization needs to be provided Such an initialization can be provided by introducing fictitious data at negative time as long as the prewindowing (data is zero before the fictitious data) and the particular structure of the (extended) input data matrix remain preserved The structure of the fictitious data should preferably be such that the filter coefficients and all other quantities to be initialized at time k - -1 can be determined without computation The simplest structure for such fictitious input data consists of one non-zero sample per input channel, positioned in such a way that there is no overlap between channels in the extended data matrix The only non- zero input samples at negative time % < 0 are -n B -1 - This results in the following diagonal initial sample covariance matrix 7 -nT Λ_ι = U- \ U_ λ = blockdιag{μ B \ nB ) \ nB+l μ F A nF ® I_}, Λ„ = dιag{\ n , λ} where ® denotes the Kronecker product (A ® B = [z4 y S]) Figure 2 is simplified by defining I ____±_ - in μ FF = λ μ p Also whenever λ 2 with variable % appears, the value for % is actually the value indicated the top row of the figure, in the same column If the DFE filter coefficients are to be initialized to a nonzero value, W_ ≠ 0, then the portion of the desired response signal at negative time should reflect this d_ι = W- \ J- \ The Least-Squares Filtering Problem An optimal (Wiener filtering) design of the filter W corresponds to minimizing the expected value (I e , mean value or statistical average) of the squared filtering error (Minimum Mean Squared Error (MMSE) design) W° - argmm E The optimization of the MSE criterion however requires knowledge of the joint second older statis tics of the desired-response and input signals In practice, these statistics are usually not available However (portions of) the desired-response and input signals are usually available So in practice the statistical averaging operation in the MSE criterion must be dropped or replaced by another averaging operation In the Least-Mean-Square (LMS) algorithm, the instantaneous squared error is taken as the criterion, and a simple steepest-descent optimization strategy with one iteration per sample period is applied The LMS algorithm uses data quite inefficiently and hence requires quite a bit of data to converge to (approximate ly) Hit optimal filter setting W° In the Least-Squares (LS) approach, the statistical averaging is replaced bv a tempoial aveiagmg So m the LS approach the filter settings are determined by minimizing the temporally averaged squared filtering error or equivalently, the sum of squared errors If signals have been accumulated from time % = 0 to time ι = k then genuine data vectors are available m the time window [rt-l, ..] (also called covariance window) and hence the filter coefficients W can be determined from the LS cost function Vv k = aig mm ∑ (W) = argmm ∑ (d. - W lh) 2 w w i— — 1 Note that for the actual data vectors (/,, the prewindowing transient period lasts one instant less than for the extended data vectors U f i+3 % The LS criterion is often solved recursively, meaning that as a now signal sample becomes available, the new filter solution is determined from the old one instead of solving the new LS criterion all over again The resulting algorithm is called the Recursive LS (RLS) algorithm In the context of the RLS algorithm it is important to note that the solution of the LS criterion only becomes well-defined when the number of terms in the sum of squared errors becomes at least equal to the number of coefficients N to be determined To make the criterion well-defined even with only one or even zero terms in the sum of squared errors, an initialization term gets added to the sum of squared errors to yield the augmented LS criterion W k - arg mm I ∑ e*(W) + (W - W- )l (W - W- X ) τ J where R_ι is some positive definite symmetric matrix Note that with zero terms in the sum of squares, the solution of this criterion yields W = VE_ι The augmented LS cost function has an optimally property from the point of view of estimation theory, applied to the estimation of the filter coefficients W Indeed this LS criterion leads to the so-called Bayosian estimate of W if the optimal error signal e,(VE°) is white Gaussian noise (with variance ξ°) that is independent of the input sιgnal(s) and prior information exists on the filter settings H 7 m the form of a Gaussian distribution with mean W_j and covariance matrix ξ°RZ) So, the larger R_ι . the closer 14 _] should be to the optimal filter W° In practice however often little information is available on W and hence W_ \ is often chosen arbitrarily (e g W- \ = 0) and R_ι is taken as small as possible (typically as small as the numerical implementation of the algorithm safely permits) In general regardless of the initialization, the solution W k of the augmented LS cost function converges to the optimal solution W° as the amount of data k tends to infinity. This is because time averages converge to statistical averages for most signals. In order to dampen even further the influence of the initialization (or in general, to permit the RLS algorithm to track possible changes with time in the joint statistics of input and desired- response signals), exponential weighting can be introduced. The effect of this is that the past gets forgotten exponentially fast. With an exponential weighting factor λ e (0, 1). the LS criterion becomes W k = arg min - W._f With λ < 1 however, the filter solution W k will never exactly converge to the optimal filter W° because even as k tends to infinity, the total effective amount of data in the cost function becomes constant, commensurate with the time constant _-^ of the exponential window. This is the price to be paid to be able to forget the past quickly. So if the goal is to optimize the filter setting W k at a certain finite time instant k, the forgetting factor λ should be carefully designed, taking both estimation quality and numerical effects into account. Finally, in order to more easily allow the derivation of a so-called fast RLS algorithm, the prewindowing transient period may be incorporated into the LS cost function, leading to: W k = arg min - W-JR- W - W__f k arg mm ∑ ^Ϋ l (W) = arg mm Σ λ'-cϊon (1) ι=-N-i = arg min 11^(^)1 = arg min \\d k - W U k w The optimal error signal m this case is d k e k ^l k {W k ) = [l -W k ] t k x t k -ι(W k ) ϋ k -- k (W k ) - ( ■k \-f-k ι(Wk-ι) -iT The orthogonality principle of LS leads to e k U k = 0 =- e e k = e k d k — e U k W k - e k d k The =o normal equations are hence where R k is an extended sample covariance matrix (as compared to the regular sample covariance matrix R k - U k U k ), and ξ k is the minimum value of the LS criterion The solution to the normal equations is k = d k d τ k - w k v k ~ i - d k d τ k - d Jliϋji y k ξ Recursive Least-Squares (RLS) Algorithm From the recursive relations the following recursion is obtained for W k W k XV fc _ι + = W k _ λ + el C k where is the predicted ' ' (or a priori) error signal (predicted using the previous filter setting W f c-i) and C is the Kalman gam Note that C k is like W k with d k replaced by σ = [1 ()• ] C k = σU k R k λ = U k R k ] The quantity a is called the pmmng vector- it pms down the most recent sample, e g , U k σ τ - U k The Kalman gain C k is the filter for estimating σ from U k with error signal given by σ-C k U k and error energy given by (σ-C k U k )σ τ = 1 - C k U k = 1 - U k R k X k = % which equals the last component of the corresponding error signal The quantity G (0, 1] is often called the hkehhood variable It can be considered to be the RLS adaptation gam Indeed, we get for the a posteriori (after updating W k ) error signal ( k = c k (W k ) = d k -W k U k = d k -W k . } U k -( p k C k U k To find a recursion for R k the Matrix Inversion Lemma is applied on R k = XR k - \ + U k U k , which leads naturally to the introduction of C k = U k λ _1 R^ j = λ C k = overnormahzed Kalman gain The RLS algorithm can then be formulated, which is a recursion for W k -ι,R k ^, -. W k ,R k l with initialization W_ \ ,Rz Ck Riccati equation The algorithm described above leads to the lowest complexity but some algorithmic variations involving C k instead of C k are possible by exploiting the relations e p k η k C k = ι\C k = e k C k c k τ lk c k = c k τ χ k c k An update equation for the LS error energy can also be obtained (if desired) where the latter transposition operation is required only if e F k would have been a vector signal (not the case for a DFE) Fast RLS Algorithms Study of the input data structure reveals that permu — V can bo introduced to isolate the latest data u k = in the extended data matrix LX+'i k with a regular data matrix U k or U k - \ being the complementary part P * A U V N+i k This leads to the following 'shift time update" structure Uk V D U N+ , k = V n V^ lh U fc- l which means that going from the old data U k -ι to the new data U k can be accomplished by shifting in new data u k (temporarily leading to an extended data matrix [/ Λ r+ 3, j t ) and shifting out old data u k . This shift time update structure allows for the derivation of fast RLS algorithms. These {ilgorithms require the introduction of forward and backward prediction quantities to account for the introduction of new and the discarding of old data, as follows: Forward Prediction Backward Prediction Error Signal e k = u k - A_ U k - \ f k = u k - Dj. U k e k = e k σ τ r k = r k σ τ Filter Error Energy k = e k e τ k - e k u k β k = r k r τ k = r k u k l The introduction of these prediction quantities allows for an alternative update of the Kalman gain and the likelihood variable that does not involve the Riccati equation. The update of the prediction quantities themselves can then be done in the same fashion as the update for the filter W and the error energy ξ k . The update of the Kalman gain and likelihood variable requires the mtroduction of an extended Kalman gam and likelihood variable C +3,fc = N+3 k R N+2)k = v+3 k O +3,λ Cw+3 i = U N+3 k X ~ R Λ +3 J f c- l , where / / v+ 3 j t = U N- t s X 7 is the extended input data vector at time k Fast Kalman algorithm The first Fast RLS algorithm (introduced in 1976) was the Fast Kalman algorithm, shown m Table 1 The Fast Kalman algorithm works with C k , the entries of which have a lower dynamic range than the entries of C k , a desirable property for fixed-point implementation In the table, the quantity C k D denotes the entries m O / v+ . f c that are m the same (column) positions as the entries of 7 3 in D k = [V 3 -D k \ D Similarly, the quantity C k A will be used, which are the entries in Cjv t 3 , f c that are in the same (column) positions as the entries of I-> in A k = \ —A^ VA The table shows the computational complexity of the algorithm, measured by the number of multiplications or additions (or approximately the number of multiply-accumulates) that need to be perfoimed These complexities are indicated for the general case of q channels In the FS DFE example, q = 3 Fast Transversal Filter (FTF) algorithm The Fast Tiansversal Filter algorithm (see Table 2) works with C k (just like the RLS algorithm), which results in the lowest computational complexity Whereas the FKA computes the backward prediction error r as the output of the backward prediction filter D f c__ the FTF algorithm avoids this inner product and finds by the manipulation of scalai (in the single channel case) quantities Table 1 The Fast Kalman Algorithm (FKA) As a result, however, the FTF algorithm is numerically less stable than the FKA in the sense that lound-off errors grow faster It is possible to combine the two ways of computing the signal r p k to create a measurement of the round-off error subsystem associateel with the algorithm, and to feed back this error measurement to stabilize the unstable oπoi propagation subsystem of the algonthm Initialization of Fast RLS In the foregoing Fast RLS algorithm, the following initialization is performed Filtering pait 1 (implicitly) 2 Usually, d_< = 0 and W-ι = 0 Pi edict ion part Table 2: The Fast Transversal Filter (FTF) algorithm. σ 1. U N+q -ι or are permutations of diagonal matrices. f -ι 2. Hence the LS estimate of any row from any other row(s) is 0. 3. For the filters, A_ x = D__ = C_, = C_ι = 0 4 For the error energies (the energy of the signal that is being estimated) μ B 0 β- — r_γr_ \ — u_ l u_ — 0 X nB + l μ F I 2 tv_ι = e_ t > i T U;_ ! TUΓ1 J - λ np /, 7-ι σσ Lowest Complexity Fast Kalman algorithm The FTF algorithm works with C k wheieas the FKA works with C k For C k — O. X ~ R k _ λ in the piewmdowmg transient phase the last non-zero coefficients for the different channels put together give which shows that for small initialization constants μ B , μp, the ovornormahzed Kalman ga C k will take on large values during the transi nt phase The uimoπnah/ed Kalman gam Cy — R k l suffeTS much less from this problem since the 'deiiommatoi" R contains the latest data U k I his problem of larger dynamic range for C k compared to C manifests itself more generally whenevei a sudden increase in the input signal level occurs This motivates (with an eye toward fixed point implementation) the use of C k over C k The algorithm shown in Table 3 uses C k as m the FKA It also computes the backwaid prediction en or by straight lorward filteiing as m the FKA (as opposed to avoiding this inner product after the manner of the FTF algonthm) This leads to a less unstable algorithm Finally, note that as a result the algorithm does not require β k (required in the FTr algoπthm) These are the most important characteristics of the FKA, and hence the resulting fast RLS algorithm can appropriately be called an FKA, albeit the fastest version of the FKA Table 3 Lowest complexity Fast Kalman Algorithm In Matlab notation (for portions of vectors) C k D = C / . +, j. [N/_ + 1, N I 2 /V I 3| whereby line 10 in Table 3 can be rewritten as C = [(C Λ + , fc PJ) [l 4 N fb \ 2 2 3, N fb +4 N-2] . C k v " D k _ Since > 1 the implementation of t he multiplication with this number on a fixed point processor may take more than one cycle The least complex computation of the LS filter estimate would be to solve the normal equations with the generalized Levmson algorithm this may be more sensitive m a hxed-pomt Prewindowing/Initialization Removal Postprocessing It is clear that avoiding the prewmdowmg assumption and the initialization (when W-ι is not particularly close to W°) would lead to much faster and better convergence (The influence of the initialization can be made small by choosing μg μp small but a fixed-point implementation especially, will pose constraints on how small these quantities can be made ) With initialization and prewmdowmg, exponential weighting needs to be applied, not so much to adapt to changing statistics but simply to forget the initialization and the prewmdowmg Without the perturbation of initialization or prewmdowmg, no exponential weighting would need to be applied This would lead to better estimation perfoimance and fewer numerical problems To be robust against changing statistics later on in time tracking performance could still be obtained without exponential forgetting by processing the data in bursts and at each burst m corporatmg as initialization the filter setting resulting from processing the previous burst with appropriate weighting (In this case the initialization would not have to be gotten rid off, only the prewmdowm ) A straightforward way to avoid prewmdowmg is to use the Growing Window Covariance (GWC) version of the FKA In this last RLS algorithm the LS problem corresponding to the covariance window m Fig 2 is solved at each time instant k This algorithm is more complex than the prewindowed version Indeed, it requires the introduction of an additional filter G the dual Kalman gam which needs to be updated also at every time instant The GWC FKA may also be more sensitive for a fixed-point implementation (due to larger dynamics of quantities) As mentioned earlier, in many applications it may not be possible to run a fast RLS algorithm in real time to get initial convergence, due to complexity considerations On the other hand, it may also be overkill to keep a fast RLS algorithm running to adapt to slow variations in the optimal filter settings In this case data should initially be acquired and stored Then a fast RLS algorithm can be run off line to converge on these data The amount of data to be stored and processed should be such that it allows the LS filter estimate to have close to optimal performance (small excess MSE) For instance, an amount of data corresponding to k - ION should be largely sufficient for equalization applications If a covanance window is used less data can be used if that leads to a more favorable memory /complexity /performance tiade off If an FKA with initialization and prewmdowmg is used, however, much more data would normally be required to get comparable performance Now, given (as described above) that fast RLS algorithms with prewmdowmg have the lowest complexity, an effective strategy is to simply run the prewindowed FKA, possibly with relatively large initialization constants μ to dampen the dynamics of the (fixed-point) quantities This will be done over a burst of data such that the corresponding covariance window LS problem gives accept able performance (so the data length is relatively short) \t the end of the burst, a postprocessing FKA algorithm is run which will now remove the initial data vectors (hence the initialization con stants m R_ι and the prewmdowmg) This will have a lower overall complexity than running a GWC version of the FKA Downdating RLS Algorithm To formulate the postprocessing algorithm data matrices are needed which not only have an end time k but also a starting time J (j = -oo m principle m the prewindowed case) With the initialization of Fig 2 a starting time _ — -oo is in fact equivalent to j — -N-q since all data prior to that instant are /ero So the filter becomes the solution of the following LS problem W jk _ e.ιgm ∑X k' i (W) ι=l = argmm || Jjt (lF)|| 2 = argmin \\d ]k - WU ]k \\ 2 w w To obtain an RLS algorithm, the problem then is to find recursions to move j forward in time until the prewmdowmg period is passed so from about_ ■ — -N-q to j - τt-1 where n = max{r.β rij } if both initialization and prewmdowmg need to be removed, or from 0 to r. — 1 if only prewmdowmg needs to be removed Since the postprocessing consists of making the data window shorter by removing data at the beginning of the window, the resulting processing can be called downdating (as opposed to the usual updating) From the recursive relations d.+ifct .- it d ] u; R J+lk = R jk -X k~J U_Uj the following recursion is obtained for W } k W ]+1 k = V, jk - X k~ > = W }k - X k - ] e ]k G ]k where G 3 k is a (eweiiioimahzed) dual Kalman gam this time not to reflect the effect of adding a data point at the end but of removing a data point at the beginning furthermore the (unnormah/ed) dual Kalman gam is G }k = U R ~ \ Note that X k~J G ]k is like W }k with d ]k replaced by τ = _ τT [0---0 λ-^] : X k~3 G k = πU k R ~ l = X k~ W R ~ j.. The quantity π is called the dualpinnmg vector it pins down (and weighs) the oldest sample. The dual Kalman gain G k is the filter for estimating π from U J<k with error signal = π- X k~3 G k U_ k and error energy = (π-X k~3 G k U k T ~ X k~~J ^ it where S k = χx k - > c J , k u, = e(o.i] Again, a relation of the form G k - S ]ιk G } , k is obtained. To find a recursion for R ~ l, the Matrix Inversion Lemma is again applied on R J+ ι. k - R k - X k~3 U ] U j l . The downdating RLS algorithm, which is a recursion for W k ,R ~ l → = RI "1 can then be formulated as follows: Riccati equation. The corresponding downdate equation for the LS error energy (if desired)is given by: where the latter transposition operation is required only if t jk is a vector signal (not the case for a DFE). Downdating Fast Kalman Algorithm To find a downdating fast RLS algorithm, the shift structure in the data can again be exploited, of which the details are now: .,* f = ? D U N+ ,, k = V D V A T Uμ u,- k-l which means that going from the old data t/ j _ ι, f c_ ι to the new data U jιk can be accomplished by shifting in new data u k (temporarily leading to an extended data matrix c / v+ 3) ) and shifting out old data U j k . To get a complete downdating strategy which will be used to downdate G j^ i k to G j . k , the regular time updating strategy (involving the regular Kalman gain j ,jt) is first used to downdate G j _i.i t to G -i t-i , and then the shift structure is used (involving foward and backward prediction filters) to update G j - \t - to G j . k . The derivations are along the same lines as for the FKA/FTF algorithms. The resulting algorithm is shown in Table 4. Just as the updating algorithm requires initialization, so also the downdating algorithm requires initialization. The following initialization (at - - 1. if the desired is to get rid of prewindowing only) may be used: W ,A,D,C,a ~1 ,^ are initialized by their values obtained from the prewindowed algorithm, and G_ = 0, 1. The initialization can be placed earlier to get rid of the soft constraint initialization (W- \ , R- \ ) also. (This is done using the same initialization values, only the signal values to be processed during the downdating become more numerous). When downdating the soft constraint initialization, the complexity can be reduced by 3N-3 if the fact that the U j contain only one non-zero entry gets exploited. Prewindowing Suppression for the Filtering Section The filtering (equalization) error signal shows large values during the initial convergence of the algorithms One would expect however, that this would not bo the case when the filter W is initialized to some setting W- \ that is relatively close to the optimal setting (as a result of previous adaptation) However things are pertuibed due to the prewmdowmg, which will load good filter settings to produce a large error signal during the prewmdowmg transient The prior art (Eleftheiiou and Falconer) describes a solution to this problem I he solution consists of still using the piewindowed input data vector U k for the prediction part of the fast RLS algorithm, so that a last algorithm is possible, but to introduce a different data vector, U k , for the filtering section U k is the same as U k except that it's not prewindowed So U k differs from U k only during the prewmdowmg transient, m which period U k contains some actual data from before time zero (assuming this data is available but was simply not used due to the algorithmic constraints of prewindowed last RLS algorithms) So in this case, the filtering error signal is actually computed as = {d k - Wk-i {U k - Uk)) - Wk-, Uk = d k - W k ^ U k where U k — U k is non- ero, and d k is diffeiont from d k , only during the prewmdowmg transient So things happen as if the prewmdowmg assumption is still used, but with a (temporarily) modified desired response signal d k which compensates as well as it can for the prewmdowmg assumption (since before updating at time k, W k - \ is our best estimate for W) I he filtering error thus computed does not get perturbed by prewmdowmg and hence does ηpt show particularly largo value s duπng the prewmdowmg (not larger than what W- \ can dehvei ) To incorporate this modification m the postprocessing it suffices to use d k ' instead of d k (U k does not need to be used) Indeed, the filter estimate would have been exactly the same if prewmdowmg (U k ) in the filtering section was used, with d k as desired response signal Partial Cancellation of the Initialization It is trivially possible to cancel the initialization only paitially This may be desirable if the initialization constants wore chosen on the basis of numeric al considerations whereas after cancellation the objective is lor the constants to take on (smalki) values that conespond to certain statistical estimation considerations (Bayesian point of view) Indeed, if for the postprocessing purposes μ_ , μ B aie set to values μ F < μ F , μ' B < μ B , then after postprocessing the initialization will still influence the filtei estimate but new according to the constants μp - μ_ , μ - μ B Updating and Tracking by Processing Data in Bursts To make up lor the fact that the available processing power may be insufficient to run last RLS algorithms in real time, data can be stored and processed in bursts of length and frequency such that the processing of a burst is finished befoie the next burst of data arrives To incorporate the results from processing previous data bursts, when processing the next burst, the (soft constraint) initialization mechanism can be used To have an appropriate weighting of the previous results the following initialization values can be used put 4 . , D C, U to 0, to 1, and leave α "1 and H unchanged This coπesponds to a pie- and postv melowed data matrix [/_] with block-diagonal (3x 3 blocks) ii_ι The advantage of such a simple initialization is that the numerical errois associated with critical (unstable) quantities are flushed out. The appropriate weighting is provided by leaving ~l unchanged. If λ = 1 is used, then the filter estimate is the result of all the data bursts processed (not in the usual LS sense though, but in some weighted version of it). Forgetting, to track possible changes in the environment, can be introduced by taking λ < 1. Alternatively, one can take λ = 1 but weigh ~] : a '1 +- λ^α -1 at initialization, where λ[ is some burst-level forgetting factor. Or one can use a combination of both weighting mechanisms. Obviously, when incorporating information from previous bursts, the initialization should not be removed by postprocessing (postprocessing should then only be applied to the prewindowing transient). Table 4: Postprocessing Fast Kalman Algorithm to transform the initialized/prewindowed LS estimate into a Covariance Window LS estimate. Use of the described postprocessing fast Kalman algorithm shows that equivalent results as conventional algorithms can be obtained using from six to sixteen times less data than conventional algorithms. It will be appreciated by those of ordinary skill in the art that the invention can be embodied in other specific forms without departing from the spirit or essential character thereof. The presently disclosed embodiments are therefore considered in all respects to be illustrative and not restrictive. The scope of the invention is indicated by the appended claims rather than the foregoing description, and all changes which come within the meaning and range of equivalents thereof are intended to be embraced therein.

Claims: WHAT IS CLAIMED IS: 1. A method of performing recursive adaptation of a digital filter computed using at least a filter coefficient vector or matrix, a desired response data vector or matrix and an input data matrix, the method comprising: storing a burst of input data; storing or deriving corresponding desired response data; using prewindowing to initialize the desired response data vector and the input data matrix; running a prewindowed fast RLS-type adaptation algorithm on the burst of data to compute filter coefficients at the end of the burst of data; and performing postprocessing to adjust the filter coefficients to remove the effects of prewindowing. 2. The method of Claim 1, wherein the filter coefficients are adjusted to be more nearly optimal in accordance with a covariance window Least Squares error criterion. 3. The method of Claim 1, wherein the digital filter is computed using a sample covariance matrix, the method further comprising: initializing the input data matrix with at least some fictitious non-zero data at positions in the input data matrix corresponding to instants i in time prior to i = 0; initializing at least one of the desired response data vector and the sample covariance matrix to agree with the input data matrix, including the fictitious data; and after running the filter, performing postprocessing to adjust the filter coefficients to remove the effects of the fictitious data. 4. The method of Claim 3, further comprising: storing a different burst of data; initializing the filter coefficient matrix to values derived from values obtained from running the filter on a previous burst of data; using prewindowing to initialize the input data matrix; running the filter on the different burst of data to compute filter coefficients at the end of the burst of data; and performing postprocessing to adjust the filter coefficients to remove the effects of prewindowing. 5. The method of Claim 4, comprising initializing at least one algorithmic quantity by adjusting a value obtained from running the filter on the previous burst of data. 6. The method of Claim 5, wherein the algorithmic quantity is the forward prediction error energy. 7. The method of Claim 6, wherein the forward prediction error energy is adjusted upward to reduce the dynamic range of at least one quantity. 8. The method of Claim 7, further comprising, during postprocessing, removing an effect of some portion of the forward prediction error energy while preserving the effect of another portion of the forward prediction error energy. 9. The method of Claim 8, comprising removing the effect of an excess portion of the forward prediction error energy equal to an amount by which the forward prediction error energy was adjusted upward. 10. The method of Claim 6, wherein the forward prediction error energy is adjusted downward to reduce, from one burst to another, the influence of prior data, while using an exponential forgetting factor that is zero or substantially less than one. 11. The method of Claim 4, further comprising modifying desired response data vector values to reduce a pronounced error transient experienced during initial stages of processing a burst of data. 12. The method of Claim 11 , further comprising saving modified desired response data vector values for use during postprocessing. 13. The method of Claim 12, further comprising using the saved modified desired response data vector values to remove the effect of modification of the values.

Inventors: PAL DEBAJYOTI, SLOCK DIRK, PAL, DEBAJYOTI, SLOCK, DIRK
Applicants: PAL DEBAJYOTI, SLOCK DIRK, VIRATA CORP, PAL, DEBAJYOTI, SLOCK, DIRK, VIRATA CORPORATION
Application year: 2001
Main classification: H03H
IPC codes: H03H21