Abstract:
Using a known or later developed time domain equalizer coefficient training algorithm, a least square solution for the time domain equalizer coefficients is taken at a starting point and iteratively improved on. In particular, the improvement is directed towards maximizing number of bits per frame loaded over the time domain equalizer coefficient choice. This can be accomplished by maximizing capacity directly rather than setting a goal to shorten the channel and hoping that the capacity will be maximized as a result.

Description:
RELATED APPLICATION DATA  
       [0001]    This application claims benefit of and priority to U.S. Application Serial No. 60/241,664 filed Oct. 19, 2000, entitled “Frequency Domain Supplemental Training of the Time Domain Equalizer for DMT,” which is incorporated herein by reference in its entirety. 
     
    
     
       BACKGROUND OF THE INVENTION  
         [0002]    1. Field of the Invention  
           [0003]    In general, the systems and methods of this invention relate to time domain equalizer training. In particular, this invention relates to systems and methods for supplemental frequency domain training of the time domain equalizer in a discrete multi-tone system.  
           [0004]    2. Description of Related Art  
           [0005]    In Discrete Multi-tone Modulation (DMT) systems, the time domain equalizer (TDQ) is a finite impulse response (FIR) filter located at the receiver side of a DSL modem. The TDQ is used to reduce the intersymbol interference (ISI). If the channel is shortened in time to have a length no greater than the length of the cyclic prefix, the intersymbol interference can be eliminated. Thus, a common method for training the TDQ in a DMT system is to jointly optimize the numerator and denominator of the autoregressive (AR) model for the channel where the order of the numerator is equal to the cyclic prefix length and the denominator is used as the TDQ setting. The training is based on transmission and reception of a known reference signal, such as the reverb signal in ADSL systems, using a least squares fit of the AR channel model.  
         SUMMARY OF THE INVENTION  
         [0006]    The supplemental training according to the exemplary systems and methods of this invention starts with the least squares solution of the time domain equalizer coefficients outlined above as its starting point, and iteratively improves on it.  
           [0007]    Specifically, the medley-based supplemental training which is the subject of this application takes as input the output of a reverb-based TDQ training algorithm. Examples of reverb-based training algorithms are described in Stuart Sandberg and Michael Tzannes, “Overlapped Discrete Multitone Modulation for High Speed Copper Wire Communications,” IEEE JSAC, vol 13, no. 9, December 1995, pg 1571-1585, incorporated herein by reference in its entirety, and include channel shortening schemes based on an AR fit to the transmission channel.  
           [0008]    The improvement is geared towards maximizing the number of bits per frame loaded over the TDQ choice. In particular, capacity is maximized directly rather than setting a goal to shorten the channel and hoping that the capacity would be maximized as a result. The supplemental training operates in medley transmission mode, and requires a number of pseudo-random data frames.  
           [0009]    Medley operation is selected in that the reverb data transmission, which is the repetitive transmission of the same reference frame, would not produce ISI in the received signal, and the SINR (Signal-to-Interference and Noise) determined in this way would not take into account the very component of error that the TDQ is intended to reduce. From the medley data, the SINR can be estimated over the bins used in the actual data transmission mode, and therefore, the number of bits per frame loaded.  
           [0010]    The systems and methods of this invention use a directed search method on the capacity function to obtain an improved TDQ. The function in question is highly nonlinear, and after linearization, e.g., using the first two terms in a Taylor series expansion around the starting TDQ point, a local extremum is sought. Since this does not guarantee the best solution, the same supplemental training can be repeated one or more times, each time starting with the TDQ solution resulting from the previous run.  
           [0011]    In accordance with an exemplary embodiment of the invention, an aspect of the invention relates to performing frequency domain supplemental training of a time domain equalizer.  
           [0012]    Additionally, aspects of the invention also relate to performing frequency domain supplemental training of a time domain equalizer in a discrete multi-tone environment.  
           [0013]    Additional aspects of the invention also relate to performing the supplemental training numerous times, with each instance of the supplemental training using the last determined time domain equalizer coefficients to improve the quality of the results.  
           [0014]    These and other features and advantages of this invention are described in, or are apparent from, the following detailed description of the embodiments. 
       
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0015]    The embodiments of the invention will be described in detail, with reference to the following figures wherein:  
         [0016]    [0016]FIG. 1 is a functional block diagram illustrating a DMT DSL modem according to this invention;  
         [0017]    [0017]FIG. 2 illustrates how supplemental training fits into a sequence of ADSL receiver initialization/training tasks; and  
         [0018]    [0018]FIG. 3 is a flowchart detailing the perform supplemental training step of FIG. 2. 
     
    
     DETAILED DESCRIPTION OF THE INVENTION  
       [0019]    The exemplary embodiments of this invention will be described in relation to the application and the invention to an ADSL transceiver environment. However, it should be appreciated that in general, the systems and methods of this invention will work equally well for any multi-carrier communication system including, but not limited to, DSL, VDSL, SDSL, HDSL, HDSL2, or any other discrete multi-tone, discrete wavelet multi-tone DSL or wireless OFDM system.  
         [0020]    As discussed above, the supplemental training according to an exemplary embodiment of this invention commences with a Least Squares solution for the TDQ and improves on the accuracy. Specifically, the number of bits per frame loaded is maximized over the TDQ choice. The function to maximize is the sum of the number of bits that can be loaded in the bins that are used for transmission, and the maximization is over the TDQ setting:  
         Max (a) (SUM (k) (log 10 (SINR k ))),  
         [0021]    where:  
         [0022]    a is the TDQ vector of size L×1,  
         [0023]    k is the bin index (out of N used bins, while N&lt;M where M is the size of the receiver Fourier Transform), and  
         [0024]    SINR k  is the signal to noise and interference ratio in bin k, expressed as a function of TDQ coefficients, a.  
         [0025]    The above is the equivalent of minimizing the following:  
         min (a) (SUM (k) (1n E[|e k | 2 ]−1n E[|s k | 2 ])),  
         [0026]    where:  
         [0027]    E[|e k | 2 ] is the mean square error in bin k,  
         [0028]    E[|s k | 2 ] is the mean square signal in bin k, and  
         [0029]    s k =u k  H k  A k  a,  
         [0030]    where:  
         [0031]    u k  is the medley 4-QAM reference symbol in bin k,  
         [0032]    H k  is an estimated (during reverb training) complex channel frequency response at bin k,  
         [0033]    A k  is a Fourier basis row vector of length L, having frequency 2πk/M, and e k =s k −F k  B a,  
         [0034]    where:  
         [0035]    F k  is a Fourier basis row vector of length M, having frequency 2πk/M,  
         [0036]    B is the received data matrix of size M×L, and  
         [0037]    each column of which is the received data frame (before the TDQ block), shifted by one sample as to represent the time passing operation. In the following, the dependence of B, u k  and s k  on frame will sometimes be shown explicitly as B(n), u k , (n), and s k (n), where n is the frame index.  
         [0038]    This function is highly nonlinear, and only the portion around the TDQ starting point is modeled by taking the first two terms of the Taylor series expansion. As a result, the function to minimize (over TDQ setting a) is:  
         SUM (k) (w k   e E[|e k | 2 ]−w k   s E[|s k | 2 ]),  
         [0039]    where the weights are:  
         [0040]    w k   e =1=E[|e k,0 | 2 ],  
         [0041]    w k   s =1/E[|s k,0 | 2 ], and  
         [0042]    e k,0  and s k,0  are e k  and s k , evaluated for the initial TDQ setting.  
         [0043]    E[.] is evaluated as an average, over medley frames.  
         [0044]    After some manipulation, the function to be optimized can be rewritten as:  
         min (a) (E[a′G e a−a′G s a])=min (a) {a′E[G e ]a−a′E[G s ]a} 
         [0045]    where:  
         [0046]    G e =G e     —     mat   + W e G e     —     mat , is a matrix of size L×L and where + is the conjugate transpose,  
         [0047]    G s =G s     —     mat   + W s G s     —     mat , is a matrix of size L×L,  
         [0048]    G s     —     mat =D u D H A, is a matrix of size N×L,  
         [0049]    D u =diagonal (u) and D H =diagonal(H), are both matrices of size N×N,  
         [0050]    A is a Fourier basis matrix of size N×L, consisting of previously described vectors A k ,  
         [0051]    G e     —     mat =G s     —     mat −FB, is a matrix of size N×L,  
         [0052]    F is a Fourier basis matrix of size N×M, consisting of previously described vectors F k ,  
         [0053]    W e =diagonal(w k   e ), W s =diagonal(w k   s ), are both matrices of size N×N, and  
         [0054]    B is a received data matrix of size M×L, as discussed above.  
         [0055]    The directed search for the minimum starts with the initial TDQ vector, a 0 , and for each iteration the TDQ is updated:  
         a i =min eigenvector {E[G e ]−E[G s ]} 
         [0056]    where E[G e ]−E[G s ] has been linearized/localized about a i−1  as described above.  
         [0057]    In practice, as discussed hereinafter, the iterations of the supplemental training are continued until arriving at a TDQ with satisfactory performance, or for some other predetermined number of iterations. Note that to obtain the TDQ for a new iteration, the TDQ from the previous iteration is used to estimate the signal, the error, and to obtain the updated matrix E[G e ]−E[G s ].  
         [0058]    [0058]FIG. 1 illustrates an exemplary DSL modem  5  according to this invention. In particular, the DSL modem  5  comprises a bit loading module  10 , an encoder  20 , an Inverse Fast Fourier Transform module  30 , a cyclic prefix module  40 , an echo canceller  50 , a digital-to-analog converter  60 , an analog-to-digital converter  70 , a time domain equalizer  80 , a training module  90 , a cyclic prefix module  100 , a Fast Fourier Transform module  110 , a frequency domain equalizer  120 , a decoder  130  and a bit loading module  140 . As will be appreciated by one of ordinary skill in the art, various other components may be present in a DSL modem, however have been omitted for the sake of clarity.  
         [0059]    While the exemplary embodiment illustrated in the FIG. 1 shows the modem  5  and various components collocated, it is to be appreciated that the various components of the modem can be combined or located at distant portions of a distributed network, such as a local area network, a wide area network, an intranet and/or the Internet, or within a modem. Thus, it should be appreciated, that the components of the modem  10  can be combined into one device or collocated on a particular node of a distributed network or combined into one or more of a CO or CPE modem. Thus, it will be appreciated from the following description, and for reasons of computational efficiency, that the components of the modem  10  can be arranged any location, such as in a general purpose computer or within a distributed network or dedicated modem without affecting the operation of the system. Furthermore, the term module as used herein is to be understood to include, but is not limited to, one or more of hardware components and/or associated software for performing a given function.  
         [0060]    In operation, the encoder, in cooperation with the bit loading module  10 , receives the input data bit stream and encodes it into M QAM constellation points. This encoding is accomplished in accordance with a bit loading table that is stored in the bit loading module  10 . The bit loading table defines the number of bits carried by each tone.  
         [0061]    The IFFT module  30  receives the encoded data and determines a sum of N carriers each modulated by a predetermined phase and amplitude. Specifically, the input to the IFFT module  30  is a vector of QAM constellation points—N complex numbers, defining the amplitude and phase of each carrier.  
         [0062]    The cyclic prefix module  40  receives the output of the IFFT module  30  and separates the received symbols in time in order to decrease the intersymbol interference (ISI). As is well known, the signal passing through the line is linearally convolved with the impulse response of the line. If the impulse response is shorter than the duration of the cyclic prefix as discussed above, each symbol can be processed separately, thereby eliminating the intersymbol interference.  
         [0063]    The echo canceller  50  generates a replica of the transmitted signal that leaks back into the receiver. Upon subtraction of the near-end echo-replica, the received far-end signal can be processed as if its only impairment has been the channel induced noise sources. In general, the echo cancellation in DSL systems considers the asymmetric upstream/downstream nature that results in different sampling rates for upstream and downstream communications. However, many variations and methods for reducing echo are well known to one of ordinary skill in the communications arts and will not be discussed herein.  
         [0064]    The time domain equalizer module  80  is a filter designed to minimize the intersymbol interference and interchannel interference (ICI). This is done by reducing the total impulse response of the line to the length of the cyclic prefix, as discussed above, such that one symbol does not interfere with the next symbol and accordingly intersymbol interference can be reduced or eliminated.  
         [0065]    The cyclic prefix module  100  complements the cyclic prefix module  40  and forwards its output to the FFT module  110 . The FFT module  110  complements the operation of the IFFT module  30  by transforming the received N carriers back into amplitude and phase information, which is then decoded back into bits in cooperation with the decoder  130  and the bit loading module  140 .  
         [0066]    The training module  90  manages a number of training features that are present in the ADSL modem system  5 . However, for the sake of clarity, only the training related to the application of this invention will be described. Clearly, one of ordinary skill in the art will appreciate that additional training will be present during the training and/or operating condition of a typical DSL modem.  
         [0067]    In particular, during a portion of initial training, the DSL modem  5  enters into reverb. During this reverb training, and in conjunction with the training module  90 , an estimate of the channel frequency response (H k ) is determined. For example, as discussed in the Sandberg article referenced above, the channel frequency response can be estimated.  
         [0068]    Next, a reverb based TDQ training algorithm, such as the one discussed the Sandberg article referenced above, is used to determine the initial TDQ coefficients. Upon determination of these coefficients, which are stored in a memory (not shown) in the training module  90 , medley is commenced. During medley, the training module  90 , operating on data received from the echo canceller  50 , performs one or more supplemental TDQ training sessions according to the systems and methods of this invention. Then, the updated time domain equalizer coefficients are provided from the training module  90  to the time domain equalizer  80 .  
         [0069]    Based on the determined TDQ, additional medley training is performed such as, but not limited to, FDQ training and SNR measurements for bit loading. At this point, the DSL modem  5  is ready to enter showtime.  
         [0070]    [0070]FIG. 2 outlines an exemplary method of performing supplemental training to determine updated time domain equalizer coefficients according to an exemplary embodiment of this invention. In particular, control begins in step S 100  and continues to step S 110 . In step S 110 , reverb is commenced. Next, in step S 120 , the channel frequency response (H k ) is estimated. Then, in step S 130 , a reverb based TDQ training algorithm is used to determine the initial TDQ coefficients. Control then continues to step S 140 .  
         [0071]    In step S 140 , medley is commenced. Next, in step S 150 , the supplemental TDQ training in accordance with this invention determines the improved time domain equalizer coefficients by maximizing the number of bits per frame. Then, in step S 160 , the updated time domain equalizer coefficients are provided to the time domain equalizer for use during showtime. Control then continues to step S 170 .  
         [0072]    In step S 170 , based on the determined TDQ, additional medley training such as frequency domain equalizer training and signal-to-noise ratio measurements for the determined time domain equalizer coefficients that are used for the bit loading is completed. Then, in step S 190 , the modem enters the showtime. Control then continues to step S 200 , where the control sequence ends.  
         [0073]    [0073]FIG. 3 outlines in greater detail the perform supplemental training block S 150  in FIG. 2. In particular, control begins in step S 500  and continues to S 510 . In step S 510 , the time domain equalizer coefficients are initialized. Next, in step S 520 , m is set to zero. Then, in step S 530 , the mean squared signal value is determined for each bin, for the given time domain equalizer coefficients. Note that the s k  is equal to the medley 4-QAM reference symbol in bin k multiplied by the estimated complex channel frequency response at bin k, obtained during reverb training, multiplied by a Fourier basis row vector of length L having frequency 2πk/m, multiplied by the time domain equalizer coefficients (a). Control then continues to step S 540 .  
         [0074]    In step S 540 , n is set equal to zero and the mean square error is also set to zero. Next, in steps S 550  through S 570 , the average error squared value is evaluated, over N 1  frames, for each bin. Thus, the signal s k  (n) and the received data B(n) are frame dependent. Then, in step S 580 , W s  and W e  will be established as diagonal matrices with elements W k   s  and W k   e . This allows localization of linearized metrics about the current time domain equalizer coefficients. Control then continues to step S 590 .  
         [0075]    In step S 590 , n is set equal to zero, and the matrices G e  and G s  are initialized to all-zeros matrices. Next, in steps S 590  through S 650 , G s  and G e , which are functions of the reference signal u k (n) and the received data B(n) for each frame, are averaged over N 2  frames. Note that a + {E[G e ]−E[G s ]} a is a linearized/localized approximation for  
         Σ log 10 (1/SNR k )  
         [0076]    which is the metric to be minimized. Then, for the next TDQ vector, the minimum Eigen vector solution is determined and the process is repeated using the updated determined localization. Control then continues to step S 690  where the control sequence ends.  
         [0077]    As illustrated in FIG. 1, the time domain equalizer coefficient determination system can be implemented either on a single program general purpose computer, or a separate program general purpose computer. However, the time domain equalizer coefficient determination system can also be implemented on a special purpose computer, a programmed microprocessor or microcontroller and peripheral integrated circuit element, an ASIC or other integrated circuit, a digital signal processor, a hard-wired electronic or logic circuit such as a discrete element circuit, a programmable logic device such as a PLD, PLA, FPGA, PAL, a modem, or the like. In general, any device capable of implementing a finite state machine that is in turn capable of implementing the flowcharts can be used to implement the time domain equalizer coefficient determination system according to this invention.  
         [0078]    Furthermore, the disclosed method may be readily implemented in software using object or object-oriented software development environments that provide source code that can be used on a variety of computer or workstation hardware platforms. Alternatively, the disclosed line time domain equalizer coefficient determination system may be implemented partially or fully in hardware using standard logic circuits or VLSI design. Whether software or hardware is used to implement the systems in accordance with this invention is dependent on the speed and/or efficiency requirements of the system, the particular function, and the particular software and/or hardware systems or microprocessor or microcomputer systems being utilized. The time domain equalizer coefficient determination system and methods illustrated herein, however, can be readily implemented in hardware and/or software using any known or later-developed systems or structures, devices and/or software by those of ordinary skill in the applicable art from the functional description provided herein and a general basic knowledge of the computer and communications arts.  
         [0079]    Moreover, the disclosed methods may be readily implemented as software executed on a programmed general purpose computer, a special purpose computer, a microprocessor, or the like. In these instances, the methods and systems of this invention can be implemented as a program embedded on a personal computer such as a Java® or CGI script, as a resource residing on a server or graphics workstation, as a routine embedded in a dedicated line characterization system, a modem, a dedicated time domain equalizer coefficient determination system, or the like. The time domain equalizer coefficient determination system can also be implemented by physically incorporating the system and method into a software and/or hardware system, such as the hardware and software systems of a time domain equalizer coefficient determination system or modem, such as a DSL modem.  
         [0080]    It is, therefore, apparent that there has been provided, in accordance with the present invention, systems and methods for determining time domain equalizer coefficients. While this invention has been described in conjunction with a number of exemplary embodiments, it is evident that many alternatives, modifications and variations would be or are apparent to those of ordinary skill in the applicable arts. Accordingly, the invention is intended to embrace all such alternatives, modifications, equivalents and variations that are within the spirit and scope of this invention.