The present invention relates to methods and arrangements to minimize error contributions in a discrete multitone modulation transmission system where a received DMT-signal is affected by error contributions.
xDSL systems, e.g. Asymmetric Digital Subscriber Line ADSL and Very-High-Speed Digital Subscriber Line VDSL, use an ordinary phone line to transmit data at high speed. When the transmitted signal passes through the line this signal will be filtered by the channel transfer function. This will cause a frequency dependent attenuation and phase-shift. The signal will be smeared out in time. A long impulse response of the channel causes one symbol to disturb the next one.
The modulation scheme standardized for ADSL and VDSL is DMT modulation. The modulation scheme partions the available bandwidth into a number of subchannels (tones) with smaller bandwidth. Each of these subchannels is considered as an independent transmission channel. Each of these channels is modulated by an Quadrature Amplitude signal QAM. The size of the QAM constellation size may vary from channel to channel.
To simplify the demodulation process in the receiver, each DMT symbol is prefixed by the cyclic prefix in the transmitter prior digital filtering and D/A conversion. The demodulation process only works in case that the memory of the channel impulse response is shorter than or at most equal to the length of the cyclic prefix. If this assumption does not hold, the tail of the previous DMT symbol contributes to the block received for the current DMT symbol and undesired Intersymbol Interference ISI. Intercarrier Interference ICI also occurs due to disturbances caused by other tones within the current DMT symbol.
To summarize, an impulse response which does not fulfil the criterion of a short impulse response leads to undesired ISI and ICI. A received block representing one DMT symbol suffers therefore not only from additive noise, but also from Intersymbol Interference and Intercarrier Interference.
In the U.S. Pat. No. 6,097,763 is disclosed a method to provide a frequency-domain training algorithm, to obtain a minimum square error equalizer that accounts for Intersymbol Interference and Noise. The algorithm in the US patent works iteratively in frequency domain. It does not account for the particularities of a DMT receiver, i.e. the FFT used for demodulation. In the US patent, sudden change of transmission conditions will also lead to inconvenience.
The present invention solves problems related to minimization of error contributions in a Discrete Multitone Modulation transmission system. The problems with earlier known techniques are due to complex algorithms which have been necessary for the minimization. Another problem solved by the invention is difficulties in modifications of an already existing solution by a sudden change of transmission conditions.
The problems are solved by the invention by calculating a cost function represented by error contributions affecting the signal received in the modulation system. After introduction of a constraint, the cost function is minimized with respect to the coefficients of a time-domain equalizer.
More in detail, adding intersymbol interference average power, intercarrier interference average power and additive noise average power, affecting the received signal, calculate the cost function. A constraint is then introduced which avoids the trivial solution i.e. when a TEQ vector that represents the time-domain equalizer is equal to zero. The cost function is thereafter minimized by optimization of the TEQ vector.
The object of the present invention is to introduce an algorithm by which error contributions in a modulation system easily can be minimized.
One advantage with the invention is that complicated calculations involving inverse matrices are avoided.
Another advantage is the possible simple modification of an already existing solution in case of a sudden variation of for example noise.
Yet another advantage with the invention is that the trivial solution (q=0) is avoided.
In
Data transmission using DMT is degraded by disturbances from three sources, i.e. Intersymbol Interference ISI, Intercarrier Interference ICI and Noise N. When the transmitted signal passes through the line this signal will be filtered by the channel transfer function. This will cause a frequency dependent attenuation and phase-shift. The signal will be smeared out in time. If the channel has a long impulse response this will cause one symbol to disturb the next one. This is called the Intersymbol Interference ISI. Another effect by the impulse response being too long is that orthogonality between tones no longer is maintained and tones will interfere with each other. This is called the Intercarrier Interference ICI. Yet another error contribution to the transmitted signal is the noise N.
A goal is to minimize the overall contributions caused by these three error sources either in frequency domain or in time domain. At first Intersymbol Interference ISI and Intercarrier Interference ICI will be accurately shown in
In
The filtered DMT-symbol (l) starts to become sinusodial (Lh−1) samples after the transmission of DMT symbol (l) has been started, where Lh denotes the length of the channel impulse response. Since the impulse response is much longer than the cyclic prefix, the sequence up to sample 33 is not sinusoidal. The FFT applied to samples 24 to 39 sees in the beginning (samples 24 to 33) a non sinusoidal sequence. Samples 34 to 39 are a pure sinusoidal sequences. The non sinusoidal part causes errors in the demodulated DMT signal.
In the algorithm according to a non-limiting, first example embodiment, a cost function J for the overall error caused by ISI, ICI and additive noise N is formulated. This cost function is minimized with respect to the coefficients of the Time Domain Equalizer TEQ. The impulse response of the channel filtered by the calculated TEQ may not have the shortest impulse response, but the overall error is minimized.
The error contributions ISI, ICI and N in the received signal will now be calculated. The calculated power contributions will then be used in the cost function J that is to be minimized according to the invention.
As a first step, equations for the error due to Intersymbol Interference ISI are set up. The error is defined by:
YISI=RXISIHq
where
Now, the power of the ISI contributions accumulated over the used tones is defined by:
PISI=E{YISIHYISI}=qTHTE{XISITRHRXISI}Hq
where E{XISITRHRXISI}=EISI is an expectation matrix depending on the used tones, the average power applied to the used tones, and on a delay parameter Δ. For each received DMT symbol the samples belonging to its cyclic prefix are removed by the cyclic prefix remover NCP. In the receiver, Δ specifies now for DMT symbol 0, the beginning of that block, which is further processed by the receiver in order to demodulate DMT symbol 0. The block which must be processed in the receiver to demodulate DMT symbol l starts now with sample Δ plus l times the length of one DMT symbol uncluding its cyclic prefix. The use of the delay parameter Δ will be further explained later together with
The next step is to set up an equation for the error due to Intercarrier Interference ICI. The error on tone k due to Intercarrier Interference is defined by:
YICIk=rkTXICIkHq
where
Now, the power of the ICI contributions on tone k is defined by:
PICIk=E{YICIk*YICIk}=qTHTE{XICIk
where E{XICIkTrk*rkTXICIk}=EICIk is an expectation matrix that depends on the used tone k, the average power applied to the used tones, and on the delay parameter Δ.
The power of ICI contributions accumulated over the used tones is
where the set ku contains the indices of the used EICI tones.
The next step is to set up an equation for the error due to additive noise. The error is defined by:
YNOISE=RWq
where
The power of the noise accumulated over the used tones is
PNOISE=E{YNOISEHYNOISE}=qTE{WTRHRW}q
where E{WTRHRW}=ENOISE is an expectation matrix that depends on the used tones and the autocorrelation function (ACF) of the noise.
According to the invention, the Time Domain Equalizer TEQ should minimize the power of the three error contributions Intersymbol Interference ISI, Intercarrier Interference ICI and additive noise N. This is in a first step done by setting up the cost function J. The cost function is defined by:
J=PISI+PICI+PNOISE=qT(HT(EISI+EISI)H+ENOISE)q
The next step is to introduce a constraint. Without any constraint, minimizing the cost function J yields the trivial solution q=0. To avoid this solution, some constraints must be introduced. In the following, four different embodiments to realize this according to the invention will be presented.
In the first embodiment the constraint is realized by fixing one tap of the overall impulse response. The overall impulse response after the Time Domain Equalizer TEQ is given by the convolution hall[n]={h*q}[n], which can be written as hall=Hq where hall contains the samples of the overall impulse response hall[n], n=0, 1, . . . , Lh+Lq−2. The matrix H is the convolution matrix of the channel. Fixing tap v of the overall impulse response to unity, this constraint can be expressed by the equation hall[v]=hvTq=1 with hvT the v-th row of H. With this linear constraint the solution of the desired Time Domain Equalizer TEQ becomes:
The optimization problem is solved by the linear system of equations
λ(HT(EISI+EICI)H+ENOISE)qopt=hr
The scaling factor λ must be chosen to satisfy hall[v]=1.
The above described solution presents an algorithm where complicated calculations involving inverse matrices are avoided. The fact that the three different powers of error contributions have been separated from each other also makes it possible to easy modify an already existing solution in case of a sudden variation of transmission conditions.
The power is defined as the quadratic form qTHTEISIHq, where the expectation matrix EISI depends on the used tones, on the average power applied to the used tones and on the delay parameter Δ. This step is disclosed in
In the second example embodiment the constraint is applied to the Time Domain Equalizer TEQ. A tap v of the vector q, where the vector q contains the TEQ coefficients, is set to unity. The desired TEQ is obtained by the optimization problem:
which is solved by the linear system of equations:
λ(HT(EISI+EICI)H+ENOISE)qopt=ev
where the vector ev is equal to “one” at position v and equal to “zero” otherwise.
The scaling factor λ must be chosen to fulfil the constraint [qopt]v=1
In the third example embodiment the constraint is applied to the used tones, i.e. the accumulated power of the desired contributions received on the used tones is fixed to unity. The desired power of the data received over all used tones is equal to Pdes. Minimization of the cost function subject to the constraint Pdes=1 yields the following optimization problem:
To solve this problem, the eigenvector dmin corresponding to the smallest eigenvalue of the matrix
D=√{square root over (HTEdesH)}−1
must be calculated first. The desired Time Domain Equalizer TEQ is then obtained by the linear transformation
with
HTEdesH=√{square root over (HTEdesH)}T√{square root over (HTEdesH)}.
The matrix Edes is defined as
where Xdesk, is the convolution matrix of the sequence xdesk[n]. The sequence xdesk[n] is that sequence to be transmitted, which only contributes desired power to tone k at DMT symbol l in the receiver. xdesk[n] only transmits DMT symbol l , all; other DMT symbols are equal to zero. Within DMT symbol l only tone k is modulated.
In the fourth example embodiment the constraint is realized by fixing the energy of the Time Domain Equalizer TEQ.
This optimization problem is solved by the eigenvector dmin corresponding to the smallest eigenvalue λmin of the matrix D=qT(HT(EISI+EICI)H+ENOISE)q and scaling according to
The method used to find the constraints in the third and fourth embodiments, responding to the preferred minimized cost function, is similar to the method described together with
Different variations are of course possible within the scope of the invention. For example can the sample frequency in the receiver be increased so that the time domain equalizer is over-sampled which will lead to a possibility to discard samples after the TEQ and thereby obtain better performance.
Instead of minimizing errors observed on the used tones, it is also possible to minimize the error observed on all tones. This approach has the advantage that the matrices representing the estimates EISI, EICI and ENOISE will be less complex. In other words, the invention is not restricted to the above described and illustrated embodiments.
Number | Date | Country | Kind |
---|---|---|---|
01850061 | Mar 2001 | EP | regional |
This application claims the benefit of Provisional Application No. 60/279,718, filed Mar. 30, 2001, the entire content of which is hereby incorporated by reference in this application.
Number | Name | Date | Kind |
---|---|---|---|
5285474 | Chow et al. | Feb 1994 | A |
6882619 | Gerakoulis | Apr 2005 | B1 |
20020131537 | Storm et al. | Sep 2002 | A1 |
Number | Date | Country |
---|---|---|
0 838 928 | Apr 1998 | EP |
0 913 973 | May 1999 | EP |
2 330 499 | Apr 1999 | GB |
Number | Date | Country | |
---|---|---|---|
20030007573 A1 | Jan 2003 | US |
Number | Date | Country | |
---|---|---|---|
60279718 | Mar 2001 | US |