Abstract

We propose the use of an iterative receiver based on the Space Alternating Generalized Expectation maximization (SAGE) algorithm for crosstalk cancellation in upstream vectored VDSL. In the absence of alien crosstalk, we show that when initialized with the frequency-domain equalizer (FEQ) output, the far-end crosstalk (FEXT) can be cancelled with no more real-time complexity than the existing linear receivers. In addition, the suggested approach does not require offline computation of the channel inverse and thus reduces the receiver complexity. In the presence of alien crosstalk, there is a significant gap between the rate performance of the linear receivers as compared with the single-user bound (SUB). The proposed receiver is shown to successfully bridge this gap while requiring only a little extracomplexity. Computer simulations are presented to validate the analysis and confirm the performance of the proposed receiver.

1. Introduction

“Very-high-speed Digital Subscriber Lines (VDSLs)” is a broadband access technology that uses twisted pairs (TPs) as a medium for high-speed data transmission [1]. It is one of the key broadband technologies for solving the “last mile” problem. With the exploitation of a high bandwidth (in tens of a megahertz), it can provide a bidirectional data rate up to 200 Mbps over short loop lengths [2]. Several TPs corresponding to a number of users are contained in a binder, which ultimately connect the central office (CO) or optical network unit (ONU) to the customer premise equipment (CPE). Because of the electromagnetic coupling among closely packed TPs, a “crosstalk” termed as far-end crosstalk (FEXT) is introduced into the far-end signal at each TP. Such crosstalk that arises from subscribers enjoying similar type of services under VDSL systems is referred to as self-crosstalk and degrades the performance significantly, especially for shorter loop lengths [1]. Another major cause of rate degradation is crosstalk that originates from subscribers enjoying other services, and is referred to as alien crosstalk [3]. Such crosstalk exists in practical situations mainly due to the coexistence of broadband over power line (BPL) systems, radio frequency interference (RFI) ingress, and crosstalk from subscribers within the same binder enjoying other DSL services.

With the advent of vectored transmission in [4], which leverages the colocation of receiver modems at the CO, there has been a surge of research interest in receiver designs for crosstalk cancellation [5]. For the design of crosstalk cancelers, computational complexity is an important issue due to the presence of a large number of tones (typically 4096) as well as a number of users per vectored group. Recently, a near-optimal linear zero forcing (ZF) receiver has been proposed in [6] for self-crosstalk cancellation, which requires channel matrix inversion at each tone. Such matrix inversions are frequently required due to changes in user status or variations in crosstalk characteristics [7], and hence cause an increase in overhead on the computational cost of the ZF receiver. To avoid this, a low-order truncated series approximation of the inverse channel matrix was considered for downstream transmissions in [8]. But it was shown that such an approach does not provide performance as good as the ZF receiver when the loop lengths are short (which is often the case with VDSL systems). The authors in [9] have suggested a pilot-based least-mean-square (LMS) tracking algorithm, which requires a large training overhead for self-crosstalk cancellation.

All the above authors assume the absence of alien crosstalk. The presence of an alien noise, originating from an external source, introduces a high spatial correlation among noise at the receivers of different vectored users [1012]. It was shown in [3] that noise correlation between twisted pairs is more than the correlation between the tones. In theory, the receivers suggested for self-crosstalk can be applied after a whitening procedure. However, the prewhitening operation applied on the spatially correlated noise destroys the columnwise diagonal dominant (CWDD) characteristics of the channel (leading to poor performance of linear receivers for alien crosstalk cancellation). The authors in [10] suggested a nonlinear successive interference canceler to achieve higher data rates than ZF receiver. However, this receiver requires 𝐐𝐑 (here 𝐐 denotes a unitary matrix, and 𝐑 is an upper triangular matrix) decomposition of the channel matrix and is quite involved computationally because of the search required for 𝐐𝐑 ordering. A computationally expensive turbo receiver based on the minimum mean-square-error (MMSE) criterion for such crosstalk cancellation was suggested in [13]. An alien crosstalk canceler was considered in [14] by assuming perfect symbol estimation after self-crosstalk cancellation. In [3], a joint transmitter-receiver cooperation framework was shown to achieve capacity for alien crosstalk cancellation. However, the proposed algorithm is dependent on knowledge of the channel at the transmitter, and cooperation between the CPEs, which is not feasible in most situations. An algorithm to mitigate a single interference from home local area network using iterative soft cancellation was suggested in [15] for use in downstream DSL. The proposed scheme does not consider vectoring and, therefore, may not be suitable for use in the upstream transmissions.

In practice, there is a need for crosstalk canceler at each tone, which can support both conventional single-user as well as multiuser detection with a limited complexity. Considering a bit loading based on the crosstalk-free signal-to-noise ratio (SNR), the performance after frequency-domain equalization (FEQ) may be degraded in terms of bit error rate (BER). However, it is important to appreciate that the estimated symbols are located within a small radius (of the order of minimum distance between constellation points) of the true symbol value. When the crosstalk is present on the victim line, it can be mitigated by utilizing the estimates of the disturbers (after FEQ) iteratively to yield a relatively smaller variance of the residual crosstalk, ensuring the minimum BER level of 10−7 as per the DSL standard. The specific CWDD property of DSL channels facilitates achievement of crosstalk-free performance and hence motivates the deployment of an iterative canceler.

With the above motivation, we propose in this paper an iterative receiver based on a space-alternating generalized expectation maximization (SAGE) algorithm for cancellation of the crosstalk in VDSL systems. Basically, the SAGE algorithm is a variant of the expectation-maximization (EM) algorithm [16] that yields performance close to the maximum likelihood (ML) solution in situations where the ML solution is computationally intractable [17]. We initially consider a situation when alien crosstalk is not present. By employing an ordered SAGE (OSAGE), we derive simple bounds on the achievable signal-to-interference-noise ratio (SINR) to help facilitate the performance analysis. Based on the derived bounds, we show that our proposed canceler provides near crosstalk-free performance while eliminating the need for channel inversion. We also show that the proposed receiver cancels the self-crosstalk by initializing the receiver with the FEQ output and requires only a single iteration to come close to the optimal performance. We next consider the case when alien crosstalk is also present. By deriving a new bound on the CWDD parameter of equivalent channel after noise whitening, we show that the convergence conditions continue to be satisfied in most situations of practical interest and can be exploited to cancel the alien crosstalk. In the presence of alien crosstalk of high power and/or low correlation, the SAGE algorithm is still shown to require only one iteration for approximating the ML solution, though a few additional iterations may be required under the conditions of low power and/or high correlation of alien noise across the TPs. This offers an attractive trade-off between data rate and complexity while canceling the self-crosstalk and mitigating alien crosstalk. Computer simulations are conducted to demonstrate the relevance of the proposed method in practical VDSL deployments.

The organization of this paper is as follows. A description of the system model is given in Section 2. Section 3 presents an iterative receiver based on the SAGE algorithm for self-crosstalk cancellation while Section 4 describes an alien crosstalk cancellation algorithm with noise prewhitening. The performance of the proposed iterative receiver based on the derived bounds as well as computer simulations is assessed in Section 5. Finally, conclusions are drawn in Section 6.

Notation
Vectors (matrices) are denoted by boldface lower (upper) case letters. 𝐴𝑖𝑗 and [𝐀]𝑖𝑗 denote the 𝑖𝑗th element of the matrix 𝐀 while 𝑋𝑖 denotes the 𝑖th element of vector 𝐱. All vectors are column vectors. The variance of random variable 𝑋 is denoted by 𝜎2𝑥. The operators (), ()𝑇, and () denote the conjugate, transpose, and conjugate transpose, respectively. The operators 𝔼{}, tr(), ||, and det() represent expectation, trace, absolute value, and determinant, respectively. We use ̃𝑧 and ̂𝑧 to differentiate 𝑧 corresponding to self- and alien crosstalk cancellation, respectively.

2. System Model

We consider a system model for upstream transmissions as shown in Figure 1. As is well known, DMT modulation based on Inverse Fast Fourier Transform (IFFT) is used in VDSL. The available frequency band is divided into a number of parallel subcarriers or tones, and IFFT effectively loads symbols onto the multiple tones. The DMT receivers (at each TP in the CO) ignore the cyclic prefix portion of the received signal and use the Fast Fourier Transform (FFT) for demodulation. Since an adaptive bit loading is used, the QAM constellation size varies according to the receive SNR. We assume that all users are perfectly synchronized and that the impulse response length of channel is no longer than the cyclic prefix length. We now consider a VDSL system with 𝑁 vectored users as shown in Figure 1. At the receiver of the 𝑖th TP, the 𝑘th sample of the FFT 𝑌𝑖,𝑘 is given by𝑌𝑖,𝑘=𝐻𝑖𝑖,𝑘𝑋𝑖,𝑘+𝑁𝑗=1,𝑗𝑖𝐻𝑖𝑗,𝑘𝑋𝑗,𝑘+𝑉𝑖,𝑘,(1) where 𝑋𝑖,𝑘 is the data symbol of the 𝑖th user and 𝑉𝑖,𝑘 is a component of the additive noise at the 𝑘th tone that includes the alien crosstalk and thermal noise. The coefficient 𝐻𝑖𝑖,𝑘 basically arises due to attenuation of the 𝑘th tone in the TP, which is usually modeled by the transmission line theory. Crosstalk coupling coefficient 𝐻𝑖𝑗,𝑘 is the complex channel element of the 𝑘th tone from the 𝑗th interferer to the 𝑖th victim and is modeled as in [18]. The second term on the right-hand side of (1) represents the FEXT. As is well known, this crosstalk is the major factor limiting the performance of DSL systems. We consider signal-level coordination in such a way that the samples 𝑌𝑖,𝑘 at the 𝑘th tone for all TPs are processed together (referred to as vectoring in [4]). Therefore, the received vector on the 𝑘th tone can be expressed as 𝐲𝑘=𝐇𝑘𝐱𝑘+𝐯𝑘,(2) where 𝐲𝑘=[𝑌1,𝑘,𝑌2,𝑘,𝑌3,𝑘,,𝑌𝑁,𝑘]𝑇, 𝐱𝑘=[𝑋1,𝑘,𝑋2,𝑘,𝑋3,𝑘,,𝑋𝑁,𝑘]𝑇, and 𝐯𝑘=[𝑉1,𝑘,𝑉2,𝑘,𝑉3,𝑘,,𝑉𝑁,𝑘]𝑇 represent the vectors of received samples, transmitted symbols, and noise samples, respectively. 𝐇𝑘 is the channel matrix for the 𝑘th tone whose diagonal elements correspond to the direct paths between the CPEs and the CO while the off-diagonal elements represent the crosstalk. The maximum ratio of the nondiagonal element to the diagonal element is defined by a parameter 𝛼cwdd𝑘 as0||𝐻𝑗𝑖,𝑘||𝛼cwdd𝑘||𝐻𝑖𝑖,𝑘||,𝑖𝑗.(3) Equation (3) implies that the magnitude of the crosstalk channel coefficient |𝐻𝑗𝑖,𝑘| from the 𝑖th disturbing transmitter into the 𝑗th receiver is always weaker than the magnitude of the corresponding direct channel |𝐻𝑖𝑖,𝑘|, which indicates the CWDD characteristic of the channel, typically observed in this context [4, 6].

In what follows, we omit the index 𝑘 for notational simplicity since crosstalk cancellation is carried out tonewise. The analysis is similar for all tones. Henceforth, we refer to 𝑌𝑖,𝑘, 𝑌𝑘, 𝑋𝑖,𝑘, 𝑋𝑘, 𝐻𝑖𝑗,𝑘, and so forth, by 𝑌𝑖, 𝑌, 𝑋𝑖, 𝑋, and 𝐻𝑖𝑗, respectively. In the next section, we assume the vector of noise samples to be spatially white, while the case of spatially correlated noise (which arises due to the alien crosstalk) is discussed in Section 4.

3. SAGE Algorithm for Self-Crosstalk Cancellation

In this section, we investigate the performance of an iterative receiver based on the SAGE algorithm [16] for the cancellation of self-crosstalk in upstream VDSL. The SAGE algorithm with user ordering is considered in Section 3.1, while its special unitary subset case is dealt with in Section 3.2.

3.1. Iterative Receiver with Ordered SAGE

The magnitude of FEXT on any tone shows a statistical variation from TP to TP because of the variations in the characteristics of twisted pairs, the nature and line length of the disturbers, and so forth. Consequently, some selected TPs are more highly affected than others at a given frequency. By giving priority to these TPs we order the users accordingly in the SAGE algorithm to increase the convergence rate and refer to this as “Ordered SAGE (OSAGE)” algorithm. Here we design our iterative technique such that the crosstalk of an ordered subset of users (grouped according to their decreasing SINR) is cancelled sequentially. As such each step of the OSAGE algorithm updates only one component of each ordered subset at a time, while keeping the estimates of the other components fixed at their previous values. By ordering 𝑁 users of set 𝑆={1,2,3,,𝑁} into 𝑀 subsets 𝑆1,𝑆2,𝑆3,,𝑆𝑀 containing 𝑁𝑆1,𝑁𝑆2,𝑁𝑆3,,𝑁𝑆𝑀 users, respectively, the algorithm at each iteration is described as follows.(1)Definition step: 𝐿=1+𝑆mod𝑀 where 𝑆=0,1,2,3,,𝑀1. (2)Maximization step: 𝑋𝑖𝑆𝐿=1𝐻𝑖𝑖𝑌𝑖𝑁𝑗=1,𝑗𝑖𝐻𝑖𝑗𝑋𝑗𝑆,𝑋𝑗𝑆𝐿=𝑋𝑗𝑆,(4) where 𝑋𝑖𝑆𝐿 is the symbol update of the 𝑖th user of the 𝐿th subset (𝐿 is used for indexing subsets) and 𝑋𝑗𝑆 is the latest available estimate of 𝑗th disturber of the set 𝑆. Each iteration uses the prior iteration's estimates to generate new estimates of interference and subtracts these recent estimates from the received signal to produce new estimates with lower interference levels. In the rest of this section, we discuss initialization of OSAGE and its performance analysis (convergence behaviour).

3.1.1. Initialization and SINR after FEQ

The initialization process in the considered OSAGE algorithm plays an important role in the receiver performance. The algorithm is initiated by the use of FEQ (a one-tap equalizer) because it has the following advantages. Firstly, it can help to exploit knowledge of the approximate symbol estimates, and secondly, it utilizes easily the accurate knowledge of direct channel coefficients. We divide 𝑌𝑖 in (1) by 𝐻𝑖𝑖 to find the FEQ estimate 𝑋0𝑖,feq of the desired symbol 𝑋𝑖 corrupted by 𝑁1 disturbers as 𝑋0𝑖,feq=𝑋𝑖+̃𝜉0𝑖,feq, wherẽ𝜉0𝑖,feq=𝑁𝑗=1,𝑗𝑖𝐻𝑖𝑗𝐻𝑖𝑖𝑋𝑗+𝑉𝑖.(5) In (5) above, ̃𝜉0𝑖,feq and 𝑉𝑖=𝑉𝑖/𝐻𝑖𝑖 represent the FEXT plus noise and noise terms after FEQ, respectively. We define FEXT𝑖=𝑁𝑗=1,𝑗𝑖|𝐻𝑖𝑗|2𝜎2𝑥,𝑗 as the total crosstalk noise on the 𝑖th user which is the sum of individual crosstalks. It is noted that 𝑋𝑗's are independent equiprobable symbols taken typically from a QAM constellation. We assume that the central limit theorem holds for ̃𝜉0𝑖,feq [15], so that it can be modeled as a complex Gaussian random variable 𝑋0𝑖,feq𝒞𝒩(0,𝜓0𝑖,feq), where𝜓0𝑖,feq=𝑁𝑗=1,𝑗𝑖||𝐻𝑖𝑗||2||𝐻𝑖𝑖||2𝜎2𝑥,𝑗+𝜎2̃𝑣,𝑖=FEXT𝑖||𝐻𝑖𝑖||2+𝜎2̃𝑣,𝑖.(6) This is reasonable since the FEXT component consists of the weighted sum of independent symbols. The SINR after FEQ can be expressed as SINRfeq𝑖=𝜎2𝑥,𝑖/𝜓0𝑖,feq. We remark that the value of SINRfeq𝑖 is low when number of disturbers is large, and specially so at higher frequencies. This is due to the fact that crosstalk power depends on both crosstalk channel gains and the signal strength of the disturbers. The SINRfeq𝑖 needs to be maximized by employing FEXT cancellation techniques. In the ideal scenario, all the crosstalk gets removed to get a parameter SNRawgn𝑖=𝜎2𝑥,𝑖/𝜎2̃𝑣,𝑖 as the SNR with only additive white Gaussian noise (AWGN). Higher rate performance is of course indicated by the single-user bound (SUB). SUB is the capacity achieved when a single-user is assumed to be transmitting and whose signal can be detected by all the receivers, either from direct channel or coupling paths [6]. However, it has been established in [6] that the SUB is only marginally higher than the data rate obtained with SNRawgn𝑖. Therefore, since the proposed OSAGE algorithm is essentially a crosstalk cancellation algorithm, it aims to achieve the performance close to crosstalk-free value, that is, SNRawgn𝑖. To quantify the gain of our designed crosstalk canceler, we define SNRgain𝑖 asSNRgain𝑖=SNRawgn𝑖SINRfeq𝑖=𝜎2𝑣,𝑖+FEXT𝑖𝜎2𝑣,𝑖.(7) For insights into the SNR gain, we provide computer simulations for 8 users of different line lengths (300 m to 1000 m) within a binder, and at different tones (4 and 12 MHz) in Figure 2. It can be seen that the SNR gain per tone for most of the users is significantly high (about 21 and 19 dB at 4 and 12 MHz, resp.). To understand the impact on data rate, using a thumb rule of 3 dB per bit for every tone, we can see that a substantial improvement in data rate is effected by FEXT cancellation, for the typically large number of tones used in the upstream VDSL system scenario.

3.1.2. Performance of OSAGE Algorithm

We now consider the performance of OSAGE by obtaining an expression for the SINR after the 𝑞th iteration, for users in each of the 𝑀 subsets. This is done by first looking at the variance 𝜓1𝑖 of the estimation error ̃𝜉1𝑖 after the first iteration for the 𝑖th user and then generalizing these to the corresponding values 𝜓𝑞𝑖 (variance of ̃𝜉𝑞𝑖) after the 𝑞th iteration, for each of the subsets 𝑆𝐿, 𝐿=1,2,3,,𝑀.

We first observe that 𝔼{̃𝜉0𝑖,feq̃𝜉0𝑗,feq} is small as compared to 𝔼{|̃𝜉0𝑖,feq|2}. To see this, we can use (5) by assuming equal transmit signal power to write𝔼̃𝜉0𝑖,feq̃𝜉0𝑗,feq=𝑁𝑗𝑖𝑁𝑙𝑗𝐻𝑖𝑗𝐻𝑗𝑙𝐻𝑖𝑖𝐻𝑗𝑗𝜎2𝑥.(8) Since 𝐻𝑖𝑗/𝐻𝑖𝑖 and 𝐻𝑗𝑙/𝐻𝑗𝑗 have small amplitudes (in view of the CWDD nature of the channel) and arbitrary phase values, 𝔼{̃𝜉0𝑖,feq̃𝜉0𝑗,feq} for a large 𝑁 is small enough to be assumed to be negligible.

As stated earlier, the OSAGE receiver performs the crosstalk cancellation on a victim user by initializing the iteration with FEQ output. We use (4) to get an estimation error ̃𝜉1𝑖𝑆1 for the 𝑖th user of the first subset (𝑆=0) after the first iteration (superscript denotes the iteration step) as̃𝜉1𝑖𝑆1=𝑋1𝑖𝑆1𝑋𝑖=𝑁𝑗=1,𝑗𝑖𝐻𝑖𝑗𝐻𝑖𝑖𝑋𝑗𝑋0𝑗,feq+𝑉𝑖.(9) The summation term in (9) corresponds to the residual crosstalk after cancellation while the second term represents AWGN after frequency equalization. It is assumed that residual crosstalk is Gaussian distributed, which is reasonable when the number of users in the binder is large due to the central limit theorem. With such an assumption on ̃𝜉1𝑖𝑆1, the residual crosstalk 𝜓1𝑖𝑆1 for the 𝑖th user can be expressed as follows:𝜓1𝑖𝑆1=𝑁𝑗=1,𝑗𝑖||𝐻𝑖𝑗||2||𝐻𝑖𝑖||2𝔼|||𝑋𝑗𝑋0𝑗,feq|||2+𝜎2̃𝑣,𝑖.(10) Substituting (6) for the 𝑗th disturber in the above expression, we get𝜓1𝑖𝑆1=𝑁𝑗=1,𝑗𝑖||𝐻𝑖𝑗||2||𝐻𝑖𝑖||2𝑁𝑝=1,𝑝𝑗||𝐻𝑗𝑝||2||𝐻𝑗𝑗||2𝜎2𝑥,𝑝+𝜎2̃𝑣,𝑗+𝜎2̃𝑣,𝑖.(11) Interchanging |𝐻𝑖𝑖|2 and |𝐻𝑗𝑗|2 in the denominator of (11) with 𝜎2𝑣,𝑖=𝜎2𝑣,𝑗 (note that we assume thermal noise variance of each TP to be equal but not so for the noise power after FEQ), the above equation can be expressed as𝜓1𝑖𝑆1=𝑁𝑗=1,𝑗𝑖||𝐻𝑖𝑗||2|𝐻𝑗𝑗|2FEXT𝑗||𝐻𝑖𝑖||2+𝜎2̃𝑣,𝑖+𝜎2̃𝑣,𝑖.(12) To simplify further, we define two channel parameters as follows:𝛼𝑖𝑗=||𝐻𝑖𝑗||||𝐻𝑗𝑗||,𝛼max𝑖,𝑗𝛼𝑖𝑗𝑖𝑗,(13) where 𝛼 is the CWDD parameter corresponding to the longest coupling length of the binder, which can be easily estimated without binder configuration. It is important to emphasise that the value of 𝛼 is independent of the binder provided the maximum reach of VDSL is fixed. By considering (6) and FEXT on the longest TP (maximum) together with (13) in (12), we get an upper bound on the residual crosstalk as𝜓1𝑖𝑆1(𝑁1)𝛼2𝜓0𝑖,feq+𝜎2̃𝑣,𝑖.(14) Writing SNRgain𝑖 of (7) as SNRgain𝑖=(𝜎2̃𝑣,𝑖+FEXT𝑖/|𝐻𝑖𝑖|2)/𝜎2̃𝑣,𝑖, the crosstalk power after FEQ of (6) can be represented in terms of SNR gain as𝜓0𝑖,feq=𝜎2̃𝑣,𝑖SNRgain𝑖.(15) On substituting (15) into (14), we have𝜓1𝑖𝑆1(𝑁1)𝛼2𝜎2̃𝑣,𝑖SNRgain𝑖+𝜎2̃𝑣,𝑖.(16) Consequently, a lower bound of SINR for the 𝑖th user of the subset 𝑆1 (𝑆=0) is given bySINR1𝑖𝑆1SNRawgn𝑖(𝑁1)𝛼2SNRgain𝑖+1.(17) Similarly, the residual crosstalk for the 𝑖th user of the second subset 𝑆2 can be expressed as𝜓1𝑖𝑆2=𝑁𝑆1𝑗=1||𝐻𝑖𝑗||2||𝐻𝑖𝑖||2𝔼|||𝑋𝑗𝑋1𝑗𝑆1|||2+𝑁𝑁𝑆1𝑗=1,𝑗𝑖||𝐻𝑖𝑗||2||𝐻𝑖𝑖||2𝔼|||𝑋𝑗𝑋0𝑗,feq|||2+𝜎2̃𝑣,𝑖.(18) By applying the similar approach as for the subset 𝑆1, the residual crosstalk for 𝑖𝑆2 (𝑆=1) is upper bounded as𝜓1𝑖𝑆2𝑁𝑆1𝛼2𝜓1𝑖𝑆1+𝑁1𝑁𝑆1𝛼2𝜎2̃𝑣,𝑖SNRgain𝑖+𝜎2̃𝑣,𝑖.(19) Doing so recursively for subsets, we can represent the residual crosstalk after the first iteration (𝑞=1) for 𝑆1 as𝜓1𝑖𝑆𝐿𝑆1𝑗=1𝑁𝑆𝑗𝛼2𝜓1𝑖𝑆𝑗+𝑁1𝑆1𝑗=1𝑁𝑆𝑗𝛼2𝜎2̃𝑣,𝑖SNRgain𝑖+𝜎2̃𝑣,𝑖(20) and the lower bound on SINR for the 𝑖th user of the subset 𝑆𝐿 asSINR1𝑖𝑆𝐿SNRawgn𝑖𝑆1𝑗=1𝑁𝑆𝑗𝛼2𝜎2̃𝑣,𝑖𝜓1𝑖𝑆𝐿+𝑁1𝑆1𝑗=1𝑁𝑆𝑗𝛼2SNRgain𝑖+1.(21)For a general iteration (𝑞2), we can find the residual crosstalk of the 𝑖th user of the first subset as𝜓𝑞𝑖𝑆1𝑀𝑗=2𝑁𝑆𝑗𝛼2𝜓𝑖𝑆𝑗𝑞1+𝑁𝑆11𝛼2𝜓𝑖𝑆1𝑞1+𝜎2̃𝑣,𝑖,(22) where 𝑀 denotes the total number of subsets. Similarly, the residual crosstalk of the 𝑖th user of the last subset (𝑀th) can be expressed as𝜓𝑞𝑖𝑆𝑀𝑀𝑗=1,𝑗𝑀𝑁𝑆𝑗𝛼2𝜓𝑖𝑆𝑗𝑞+𝑁𝑆𝑀1𝛼2𝜓𝑖𝑆𝑀𝑞1+𝜎2̃𝑣,𝑖.(23) It can be seen from (22) and (23) that users of the last subset enjoy the maximum benefit of convergence while those of the first subset the least. This, of course, is acceptable in practical scenarios since the first subset is associated with users with good SINR. For a general subset (1𝐿𝑀), residual crosstalk can be obtained as𝜓𝑞𝑖𝑆𝐿𝑀𝑗=1,𝑗𝐿,𝑗<𝐿𝑁𝑆𝑗𝛼2𝜓𝑖𝑆𝑗𝑞+𝑁𝑆𝐿1𝛼2𝜓𝑖𝑆𝐿𝑞1+𝑀𝑗=1,𝑗𝐿,𝑗>𝐿𝑁𝑆𝑗𝛼2𝜓𝑖𝑆𝑗𝑞1+𝜎2̃𝑣,𝑖.(24) Using (24), a lower bound on SINR after the 𝑞th iteration (𝑞2) can be obtained asSINR𝑞𝑖𝑆𝐿SNRawgn𝑖𝑀𝑗=1,𝑗𝐿,𝑗<𝐿𝑁𝑆𝑗𝛼2𝜎2̃𝑣,𝑖𝜓𝑞𝑖𝑆𝑗+𝑁𝑆𝐿1𝛼2𝜎2̃𝑣,𝑖𝜓𝑞1𝑖𝑆𝐿+𝑀𝑗=1,𝑗𝐿,𝑗>𝐿𝑁𝑆𝑗𝛼2𝜎2̃𝑣,𝑖𝜓𝑞1𝑖𝑆𝑗+1.(25)Combining (17), (21), and (25) together generalizes the analysis of the considered iterative receiver for each user of every subset after any given iteration step.

Some of the important implications and conclusions of this analysis are taken up in the following remarks.

Remark 1. The bound on SINR serves as a performance predictor because our focus is on achieving high data rate for a given quality of service (BER is usually fixed at 10−7). It can be seen that the data rate for a practical DSL system with an SNR gap of Γ and 𝐾 tones having spacing of Δ𝑓 is DataRate=𝐾𝑘=1Δ𝑓log21+Γ1SINR𝑞𝑖𝑆𝐿,𝑘.(26) It follows that the performance of OSAGE can be assessed by the manner in which SINR𝑞𝑖𝑆𝐿 approaches SNRawgn𝑖 at each tone. It is useful here to consider this behaviour through a practical example. Using typical values of 𝛼102, SNRgain𝑖20 dB, and 𝑁=25 in (17), it is found that the SNR loss (with respect to SNRawgn𝑖) for user of the first subset after the first iteration is (𝑁1)𝛼2SNRgain𝑖0.24 (approx. 0.93 dB) and eventually approaches zero for the users of subsequent subsets. The resulting effect on the data rate is very small, as computed from (26). Thus, our proposed algorithm effectively cancels the crosstalk with a single iteration, with an associated computational complexity of 𝒪(𝑁2) per tone.
In contrast to the ML receiver which requires a computational complexity of 𝒪(𝐾𝒞𝑁) (large constellation size 𝒞 is often used in VDSL), the complexity of the OSAGE algorithm of the 𝑞th iteration is 𝒪(𝑞𝐾𝑁2). With a single iteration, online complexity of the considered receiver is similar to that of the ZF receiver, which incurs a computational cost of 𝒪(𝐾𝑁2). Since the ZF receiver requires a computation of inverse of the channel matrix, the proposed iterative receiver promises to be computationally efficient due to its ability to avoid this offline computation.

Remark 2. As convergence is an important concern for any iterative technique, we now discuss the condition for convergence for our proposed OSAGE-based iterative algorithm. For ease of presentation, we consider a case of single subset (𝑀=1) in (22) and apply 𝜓𝑞𝑖𝑆1<𝜓𝑞1𝑖𝑆1 to get (𝑁1)𝛼2<(SNRgain𝑖1)/SNRgain𝑖 as a necessary and sufficient condition for the convergence of the OSAGE algorithm. In a crosstalk-limited DSL system (SNRgain𝑖1), the condition reduces to (𝑁1)𝛼2<1 for the convergence of the considered iterative receiver. It is worth mentioning that this convergence criterion holds good even for the case of a large number of users since 𝛼 usually takes a small value (typically of the order of 10−2). In Figure 3, we plot (𝑁1)𝛼2 for different values of 𝑁 and tones for various loop lengths. It can be seen that the convergence condition is readily satisfied.

Remark 3. The rapid convergence shown in Remark 2 can be also explained intuitively. Considering an initial bit loading (assuming an AWGN channel) in the presence of crosstalk, the estimated symbols after FEQ may be in error. However, the true symbols lie within a small radius, of the order of minimum distance 𝑑min between the constellation points. If the crosstalk due to the 𝑗th TP is removed using the FEQ estimates (by subtracting 𝐻𝑖𝑗𝑋𝑗,feq), the variance of the residual crosstalk is a small multiple of |𝐻𝑖𝑗|2𝑑2min. This is much smaller than the variance of the original crosstalk (|𝐻𝑖𝑗|2𝜎2𝑥,𝑗) since 𝑑2min is a small fraction of 𝜎2𝑥,𝑗 due to the high bit loading (bit loading of 13-14 bits per tone is not uncommon in VDSL). This coupled with the fact that the crosstalk channel coefficients are small due to the CWDD property accounts for the surprisingly fast convergence for the OSAGE algorithm.

Remark 4. Since the bounds in (17), (21), and (25) are derived by considering the maximum possible values of the CWDD parameter and FEXT, these are tight in most of the situations. This can be verified from the fact that the CWDD parameter as well as the FEXT (hence SNR gain) is not very sensitive to the variation of line length (except at shorter line length as shown in Figure 2). The tightness of bound is obvious in the case of a typical VDSL equal length lines scenario as both 𝛼 and FEXT remain almost same for all users in this case. Furthermore, in the near-far problem (that occurs when the binder contains loops with widely varying lines), equal FEXT for different users can be achieved by employing upstream power back-off (UPBO) under which the upstream transmitters vary their power spectral densities (PSDs) in accordance with the line lengths.

Remark 5. From a system design viewpoint, our derived bounds are of considerable importance. The main design parameters considered here are SNR gain and 𝑞. The SNR gain can be changed at the level of designer discretion as per the requirement of user data rates. The iterating parameter 𝑞 indicates the complexity. The trade-off between data rate and complexity can be realized by varying these parameters. Thus, the bounds in (17), (21), and (25) are useful tools for theoretical analysis of the performance of the considered canceler. Moreover, the bounds depend on local parameters such as binder size, CWDD parameter, and SNRawgn𝑖, and not on the binder configuration. This together with independent single-user water filling algorithms on each TP greatly simplifies the adaptive bit loading for the SRA operation. Here, SRA stands for seamless rate adaptation whereby the receiver monitors the SNR of the channel, and according to channel conditions, sends a message to the transmitter to initiate a change in bit loading.

3.2. The USAGE Algorithm and Its Performance

From the discussion in the previous section, we observe that a near crosstalk-free performance can be achieved even by considering only the first subset in the OSAGE algorithm. Motivated by this observation, we consider a special case of OSAGE, termed as “Unitary-subset SAGE (USAGE),” algorithm in this section in which all users are grouped into a single subset. With only one subset, each user is updated at the same time to realize a simultaneous cancellation of crosstalk terms. It offers the advantage of parallel implementation with savings in computation time but at the price of some performance loss. For its performance analysis, we can express the symbol estimates for the 𝑖th user at the 𝑞th iteration using (4) as𝑋𝑞𝑖=1𝐻𝑖𝑖𝑌𝑖𝑁𝑗=1,𝑗𝑖𝐻𝑖𝑗𝑋𝑞1𝑗,(27) where 𝑋𝑞1𝑗 is the symbol estimate of a disturber after the (𝑞1)th iteration. It should be noted that 𝑋0𝑗=𝑋0𝑗,feq represents the post-FEQ estimate of the 𝑗th crosstalker. Using (27), we obtain the power 𝜓𝑞𝑖 of residual crosstalk as 𝜓𝑞𝑖=𝑁𝑗=1,𝑗𝑖||𝐻𝑖𝑗||2||𝐻𝑖𝑖||2𝔼|||𝑋𝑗𝑋𝑞1𝑗|||2+𝜎2̃𝑣,𝑖.(28) Since 𝔼{|𝑋𝑗𝑋𝑞1𝑗|2} is the residual crosstalk at the (𝑞1)th iteration, we can represent (28) in the form of a recursive relation via𝜓𝑞𝑖=𝑁𝑗=1,𝑗𝑖||𝐻𝑖𝑗||2||𝐻𝑖𝑖||2𝜓𝑞1𝑗+𝜎2̃𝑣,𝑖.(29) By representing all the terms of 𝜓𝑞1𝑗, we can reexpress (29) as𝜓𝑞𝑖=𝑁𝑗=1,𝑗𝑖𝑁𝑝=1,𝑝𝑗𝑁𝑚=1,𝑚𝑝||𝐻𝑖𝑗||2||𝐻𝑖𝑖||2||𝐻𝑗𝑝||2||𝐻𝑗𝑗||2||𝐻𝑝𝑚||2||𝐻𝑝𝑝||2𝑞summationterms×FEXT𝑚||𝐻𝑚𝑚||2+𝜎2̃𝑣,𝑚+𝒮𝑞𝑁+𝜎2̃𝑣,𝑖,(30)where𝒮𝑞𝑁=𝑁𝑗=1,𝑗𝑖||𝐻𝑖𝑗||2||𝐻𝑖𝑖||2𝜎2̃𝑣,𝑗+𝑁𝑗=1,𝑗𝑖𝑁𝑝=1,𝑝𝑗||𝐻𝑖𝑗||2||𝐻𝑖𝑖||2|𝐻𝑗𝑝|2||𝐻𝑗𝑗||2𝜎2̃𝑣,𝑝++𝑁𝑗=1,𝑗𝑖𝑁𝑝=1,𝑝𝑗𝑁𝑚=1,𝑚𝑝||𝐻𝑖𝑗||2||𝐻𝑖𝑖||2||𝐻𝑗𝑝||2||𝐻𝑗𝑗||2||𝐻𝑝𝑚||2||𝐻𝑝𝑝||2𝑞summationterms𝜎2̃𝑣,𝑚𝑞1summationterms.(31)Assuming equal thermal noise at each TP, we can represent (31) in terms of 𝜎2̃𝑣,𝑖 (the noise after FEQ for 𝑖th user) and 𝛼 (the CWDD parameter of the binder) by making use of 𝜎2̃𝑣,𝑗=𝜎2𝑣,𝑗/|𝐻𝑗𝑗|, 𝜎2̃𝑣,𝑝=𝜎2𝑣,𝑝/|𝐻𝑝𝑝|,, 𝜎2̃𝑣,𝑚=𝜎2𝑣,𝑚/|𝐻𝑚𝑚| and interchanging the diagonal terms to obtain𝒮𝑞𝑁=𝒮𝑞𝑁(𝛼)𝜎2̃𝑣,𝑖,(32) where𝒮𝑞𝑁(𝛼)𝑞1𝑟=0(𝑁1)𝛼2𝑟1=(𝑁1)𝛼2𝑞(𝑁1)𝛼2(𝑁1)𝛼21.(33) Using (30), we can obtain residual crosstalk of 𝑖th user after the 𝑞th iteration in terms of initial crosstalk (post-FEQ) iteration by interchanging |𝐻𝑖𝑖|2 and |𝐻𝑚𝑚|2 and replacing FEXT𝑚 with the FEXT on the longest wire (maximum) together with appropriate interchanges of other diagonal terms, as𝜓𝑞𝑖(𝑁1)𝛼2𝑞𝜓0𝑖+𝒮𝑞𝑁+𝜎2̃𝑣,𝑖.(34) By substituting 𝜓0𝑖=SNRgain𝑖𝜎2̃𝑣,𝑖 from (15) in (34) and using the result along with (32), we obtain a lower bound of SINR after the 𝑞th iteration asSINR𝑞𝑖SNRawgn𝑖(𝑁1)𝛼2𝑞SNRgain𝑖+𝒮𝑞𝑁(𝛼)+1.(35) Although the USAGE algorithm can achieve a near crosstalk-free performance with a single iteration in most situations, some performance loss for the case of shorter line lengths may occur, which is quantified as𝛿SNRdB10log10(𝑁1)𝛼2𝑞SNRgain𝑖+𝒮𝑞𝑁(𝛼)+1,(36) where𝛿SNRdB=SNRawgn𝑖dBSINR𝑞𝑖dB.(37) At shorter loop lengths, the loss can be as much as 3 dB for a typical number of users. Improvement in performance however can be realized by increasing the iterations since (𝑁1)𝛼2<1 ensures convergence. This convergence condition (obtained by applying 𝜓𝑞𝑖<𝜓𝑞1𝑖 from (34) and using SNRgain𝑖1) is similar to that obtained in Remark 2. Under the condition of convergence, USAGE algorithm achieves rate performance close to crosstalk-free rate since when 𝛼0 and 𝑞, 𝛿SNRdB0, andSINR𝑞𝑖SNRawgn𝑖.(38) The lower bound of (35) is simple and can be used to obtain the data rate without the explicit knowledge of the crosstalk coupling coefficients.

4. Iterative Receiver for Alien Crosstalk Cancellation

So far we have studied the cancellation of self-crosstalk assuming the presence of only white thermal noise. However, as discussed earlier noise from other external sources (alien crosstalk) may significantly affect the functionality of the vectored VDSL systems. In practical scenarios, such a crosstalk may arise due to the use of TPs and power lines over the same spectrum for broadband access. Further, the TPs near the customer premises are generally unshielded so that RFI from nearby radio transmitter couples with these wires to induce alien crosstalk. There is also a crosstalk due to the presence of nonvectored DSL lines (usually unsynchronized) within the same binder that is carrying the vectored lines. This alien crosstalk manifests itself as spatially correlated (i.e., correlated across TPs) additive noise at each tone at the receiver of the considered VDSL system and hence degrades its performance even with the effective cancellation of self-crosstalk. Modeling of alien crosstalk is discussed in [3, 10]. Alien crosstalk is easily amenable in the vectoring framework (due to the spatial correlation between different noise samples across the TPs). In order to facilitate this, we study the use of iterative technique of the previous section for the self-crosstalk cancellation but preceded with a prewhitening filter. Our aim is to develop an effective crosstalk canceler for VDSL systems in order to achieve considerable improvement in performance in the presence of alien crosstalk. The issue here is whether the prewhitening operation affects the CWDD property of the effective channel to enable the SAGE algorithms to converge fast on the postwhitened signals. We, therefore, investigate analytically the behavior of the modified channel (channel together with the noise whitening filter) in the next section and the receiver performance based on the modified channel.

4.1. CWDD Characteristic of Equivalent Channel after Whitening

To proceed, we denote the magnitude and phase of the correlation coefficient between the noise samples on various TPs by 𝜌 and 𝜙, respectively and the alien noise power by 𝜎2𝑎. We define a noise covariance matrix 𝐑[𝑟𝑖𝑗] whose 𝑖𝑗th element is given as 𝑟𝑖𝑖=𝜎2𝑣,𝑖+𝜎2𝑎,𝑖,𝑟𝑖𝑗=𝜌𝑖𝑗𝜎𝑎,𝑖𝜎𝑎,𝑗𝑒𝚥𝜙𝑖𝑗𝑖𝑗,𝑟𝑗𝑖=𝑟𝑖𝑗𝑗𝑖,(39) where 𝚥=1. We use a noise whitening filter at the receiver which multiplies (2) with 𝐑1/2 to obtain̂𝐲=𝐑1/2𝐇𝐱+𝐑1/2𝐯=𝐇𝐱+̂𝐯,(40) where 𝐇=𝐑1/2𝐇 is the postwhitening, equivalent channel matrix and ̂𝐯=𝐑1/2𝐯 represents the whitened noise with correlation matrix 𝔼{̂𝐯̂𝐯𝐻}=𝐈𝑁. If the 𝑖𝑗th element of a matrix 𝐀 is 𝐴𝑖𝑗, then define |𝐀| as a matrix with the 𝑖𝑗th element |𝐴𝑖𝑗|. We also define 𝐀(𝑖,) and 𝐀(,𝑗) as the 𝑖th row and the 𝑗th column of the matrix 𝐀. By taking the absolute value and applying the triangular inequality on each element of the equivalent channel 𝐇=𝐑1/2𝐇, that is, |𝐻𝑖𝑗|=|𝐑1/2(𝑖,)𝐇(,𝑗)||𝐑1/2(𝑖,)||𝐇(,𝑗)| for all 𝑖,𝑗, we define the absolute value of the modified channel matrix 𝐇 as an inequality, given as|||𝐇|||||𝐑1/2||||𝐇||.(41) Using the definition of matrix inverse, we represent the elements of |𝐑1/2| as|||𝐑1/2𝑖𝑗|||=|||det𝐑1/2𝑗𝑖|||||det𝐑1/2||,(42) where 𝐑1/2𝑗𝑖 is the (𝑁1)×(𝑁1) submatrix by removing the 𝑗th row and the 𝑖th column from the square root of the covariance matrix 𝐑1/2. We define𝛾𝑟=max𝑖|||det𝐑1/2𝑖𝑖|||,𝛾𝑡=max𝑖𝑗|||det𝐑1/2𝑗𝑖,𝑗𝑖|||(43) as the maximum absolute value of determinant of principal and nonprincipal submatrices of 𝐑1/2, respectively. Now it can be shown using (39) that 𝛾𝑟>𝛾𝑡 for the given covariance matrix 𝐑. By using (42) in (41) along with the results of (43), we get a bound on the absolute value of elements of 𝐇 as||𝐻𝑗𝑖||𝐺𝛾𝑟+(𝑁1)𝛼𝛾𝑡if𝑗=𝑖,𝐺𝛼𝛾𝑟+𝛾𝑡+(𝑁2)𝛼𝛾𝑡if𝑗𝑖,(44) where 𝐺=|𝐻𝑖𝑖|/|det(𝐑1/2)| and 𝛼 as defined in (13). We use (44) to find an upper bound on ̂𝛼 for the equivalent channel aŝ𝛼𝛼𝛾𝑟+𝛾𝑡+(𝑁2)𝛼𝛾𝑡𝛾𝑟+(𝑁1)𝛼𝛾𝑡=𝛼+𝜂+(𝑁2)𝛼𝜂1+(𝑁1)𝛼𝜂,(45) where 𝜂=𝛾𝑡/𝛾𝑟. From (45), one can verify that ̂𝛼=𝛼 for the case of uncorrelated noise. Since 0𝜂<1, it is found that ̂𝛼𝛼 for all the values of 𝛼. It can be however shown that ̂𝛼<1 as long as 𝛼+𝜂(1𝛼)<1 holds for 𝛼<1 and 0𝜂<1. Although this does not clearly establish the CWDD property ̂𝛼1, it is seen that for typical practical values of 𝑁, alien noise powers and correlation parameters, the required convergence conditions for the convergence of OSAGE and USAGE are satisfied. This is confirmed later in the simulation results presented in Section 5.

4.2. Iterative Cancellation of Alien Crosstalk

Alien crosstalk cancellation after noise whitening follows analysis similar to that for self-crosstalk. To highlight the difference, we provide the distinct parameters in terms of the characteristics of alien noise. The CWDD characteristic of the equivalent channel has been already considered in the previous section. The noise variance is given by ̂𝜎2𝑣,𝑖=1 after whitening, while after FEQ it is given as ̂𝜎2̃𝑣,𝑖=|𝐻𝑖𝑖|2. Making use of (44), we find an upper bound on SNR after noise whitening (denoted by SNRwhite𝑖) asSNRwhite𝑖=||𝐻𝑖𝑖||2𝜎2𝑥,𝑖||𝐻𝑖𝑖||2𝜎2𝑥,𝑖||det(𝐑)||𝛾𝑟+(𝑁1)𝛼𝛾𝑡2.(46) We provide the bounds on |det(𝐑)|, 𝛾𝑟, and 𝛾𝑡 in (A.1) and (A.6), (A.2), and (A.10), respectively, with proofs given in the appendix. The FEXT on the 𝑖th TP after whitening filter can be represented asFEXT𝑖=𝑁𝑗=1,𝑗𝑖||𝐻𝑖𝑗||2𝜎2𝑥,𝑗=𝑁𝑗=1,𝑗𝑖̂𝛼2𝑗SNRwhite𝑗.(47) Substituting ̂𝜎2𝑣,𝑖=1 in (7) along with the result of (47), we can express SNR gain of the 𝑖th user after noise whitening asSNRgain𝑖=1+𝑁𝑗=1,𝑗𝑖FEXT𝑖=1+𝑁𝑗=1,𝑗𝑖̂𝛼2𝑗SNRwhite𝑗.(48) The SINR bounds for the present case can be obtained by simply invoking the above terms in (17), (21), (25), and (35). To show the effectiveness of the proposed iterative receiver for alien crosstalk cancellation, we study the case of equidistant (equal line length) VDSL system. By representing (48) for equidistant TPs as SNRgain𝑖1+(𝑁1)̂𝛼2SNRwhite𝑖 and substituting this in (17), we find the SINR of the OSAGE algorithm for the 𝑖th user of the first subset after first iteration asSINR1𝑖𝑆1SNRwhite𝑖(𝑁1)2̂𝛼4SNRwhite𝑖+(𝑁1)̂𝛼2+1.(49) Similarly, using (48) along with the result of (32) in (35), a simple lower bound on the performance of the USAGE algorithm for the case of alien crosstalk can be obtained asSINR𝑞𝑖SNRwhite𝑖(𝑁1)̂𝛼2𝑞+1SNRwhite𝑖+𝑞1𝑟=0(𝑁1)̂𝛼2𝑟+1.(50) It can be seen from the above expressions that performance depends on noise correlation and power of alien disturbers, since SNRwhite𝑖 and ̂𝛼 are functions of these quantities. Specifically, SNRwhite𝑖 and ̂𝛼 increase with correlation while the former decreases significantly and the latter increases marginally with alien power. With the help of these facts, the proposed iterative receiver can be designed to achieve the near crosstalk-free performance. It may be expected from (49) and (50) that fast convergence (i.e., within one iteration) may be exhibited under low correlation and/or high alien power. However, additional iterations may be required under conditions of high correlation and/or low alien noise power scenarios.

5. Numerical and Simulation Results

In this section, we present numerical results based on the analytical expressions derived here, along with MATLAB simulation experiments to investigate the performance of the proposed iterative receiver. For performing the simulations, we adopt the stochastic channel model [19] and consider the binder to consist of eight lines such that there are seven disturbers per line. We have simulated the proposed algorithm for various VDSL scenarios in the context of real VDSL deployments. Scenario 1 deals with the distributed case consisting of eight VDSL users with lines varying from 300 m to 1000 m in 100 m increments. Scenario 2 includes the case of equidistant lines. In crosstalk limited DSL systems, a near-far problem may occur with the coexistence of long lines with short lines in the same binder. This case is considered under scenario 3 which includes 4 near-end users (at 400 m) and 4 far-end users (located at 800 m) from the CO. For the analysis of OSAGE algorithm, we assume four subsets with sizes 2, 2, 3, 1 with appropriate ordering, for equidistant scenario and with 2 elements per subset for other scenarios at each tone by giving priority to the users with shorter line lengths. A bandplan 998 ADE is incorporated with 3 upstream bands US0 (25–138 kHz), US1 (3.75–5.2 MHz), and US2 (8.5–12 MHz). Other simulation parameters are used in accordance with DSL standards [2] as listed in Table 1.

5.1. Self-Crosstalk Cancellation

In this section, the performance of the proposed (OSAGE and USAGE) receiver is investigated for self-crosstalk cancellation. We consider mainly the performance after a single iteration, in order to keep the real-time complexity of 𝒪(𝑁2) per tone. We consider SINR and the corresponding data rate characteristics for the numerical studies. The analytical lower bounds, derived in (17) and (21) for OSAGE algorithm and in (35) for USAGE algorithm (with 𝑞=1), are also validated through the simulations.

The main observations that can be drawn from the numerical investigations and simulations are summarized below. (i)First, we consider a practical VDSL scenario such that the binder consists of distributed users with varying line lengths from the CO. For this case, we have simulated the SINR and the corresponding data rate characteristics for the considered algorithms, and the results are illustrated in Figures 4 and 5, respectively. From these figures, we see that both algorithms achieve performance close to crosstalk-free performance for a broad range of line lengths and bandwidths. It is further observed that the proposed algorithms have excellent performance for lower- and medium-frequency tones. Although there is a small loss in SNR with the USAGE algorithm (Figure 4) at the higher-frequency tones, the overall effect on data rate is negligible over the entire bandwidth as shown in Figure 5. (ii)Next, we consider a simplified scenario of equidistant lines, which can be feasible in some scenarios. Figures 6 and 7 incorporate this case and show that the USAGE algorithm works well and achieves close to crosstalk-free performance for longer loop lengths. Even for shorter loop lengths, it gives satisfactory performance for lower-frequency tones. However, the OSAGE algorithm always outperforms USAGE and achieves date rate close to crosstalk-free performance as shown in Figure 7. (iii)Finally, we address the performance of our proposed algorithms in a mixed scenario, to consider the impact of the near-far problem. This problem occurs when the binder contains loops with widely varying line lengths, which can be alleviated in practice by effective transmit power controls of nearby users. However, as shown in Figure 8, our algorithms offer far-end users to achieve crosstalk-free performance without compromising the data rates for the near-end users, thus potentially avoiding the need for such power controls. (iv)The ZF receiver exploits the CWDD property of DSL channels and achieves close to self-crosstalk-free performance. As stated earlier, it requires additional computation of channel matrix inversions. Further, it has been shown in many publications [6, 20] that linear receivers (ZF and MMSE) experience significant performance degradation in the presence of alien crosstalk. The SAGE algorithm can effectively mitigate the effect of alien crosstalk as discussed in next section.

5.2. Alien Crosstalk Cancellation

For this scenario, we investigate the performance of the proposed OSAGE receiver for various alien crosstalk powers and noise correlations by assuming equal length TPs. We consider correlation coefficient values for the alien noise components in the various TPs, in the range of 0.50–0.99. We have also taken an alien power of −100 dBm/Hz as it lies between the levels of −120 dBm/Hz (alien-free environment) and −80 dBm/Hz (indicating a very strong alien crosstalk) [21]. We assume known covariance matrix 𝐑 and obtain numerically the parameters 𝛾𝑟 and 𝛾𝑡 in (43) and |det(𝐑)|.

We demonstrate the effect of whitening on the original channel through the CWDD behaviour of modified channel (after noise whitening) in Figures 9 and 10. Simulation results and upper bound (derived in (45)) of the CWDD parameter are also shown in Figure 9. It can be observed that the CWDD parameter increases with the increase in both alien power and correlation, but the effect of the former is insignificant. It can be seen that the CWDD characteristic of the modified channel is weaker as compared to that of the original channel. However, even the weaker CWDD property is seen to satisfy the convergence condition (for typical values of 𝑁) and hence can be utilized for alien crosstalk cancellation, via the OSAGE receiver.

Figures 11 and 12 highlight the receiver performance in the presence of alien crosstalk. It is observed from the figures that high spatial correlation can be efficiently utilized in alien crosstalk cancellation. The number of iterations required for such a cancellation follows the similar features as for self-crosstalk cancellation, depending on noise correlation and alien power. For achieving a higher data rate improvement at shorter loop length, the proposed iterative receiver can take 2-3 iterations when the alien crosstalk has a low power and high noise correlation. However, a single iteration is sufficient to mitigate crosstalk of high alien power and/or low correlation (which corresponds to a lower data rate improvement after noise whitening).

6. Conclusions

We have investigated the use of an iterative receiver based on the SAGE algorithm for crosstalk cancellation in upstream VDSL. The proposed receiver was shown to be practically feasible and computationally efficient. By employing user ordering, a lower bound (after each iteration) on the performance of the SAGE receiver was derived. Since the derived bound depends only on the binder size, the longest line, and crosstalk-free SNR, it is easier for designer to predict the data rate and thus achieve a specified quality of service. Analytical and simulation results confirm that the SAGE-based receiver operates close to self-crosstalk-free performance with a single iteration (and complexity comparable to that of the linear ZF receiver) while eliminating the need for channel inversion. An upper bound on the CWDD parameter of the modified channel (after noise whitening) was obtained. This CWDD property after noise whitening was shown to be retained, a fact that was exploited for the cancellation of alien crosstalk. Performance of the proposed receiver was shown to be dependent on noise correlation and alien crosstalk power. As such mitigation of alien crosstalk requires only one iteration for high alien power and/or low noise correlation while a few more iterations may sometimes be required for the case of low alien noise power and/or high noise correlation.

Appendix

In this appendix, we present the bounds on modulus of the determinant of the covariance matrix and its principal and nonprincipal submatrices. For simplicity, we assume equal alien powers and thermal noises at each TP. Applying Hadamard inequality (|det(𝐑)|(𝑁𝑖=1𝑁𝑗=1|𝑟𝑖𝑗|2)1/2) yields an upper bound||det(𝐑)||𝜎2𝑣+𝜎2𝑎2+(𝑁1)𝜌2𝜎4𝑎𝑁/2.(A.1) Similarly, an upper bound on the modulus of determinant of the principal submatrix (size (𝑁1)×(𝑁1)) can be obtained as𝛾𝑟=|||det𝐑1/2𝑖𝑖||||||det𝐑𝑖𝑖|||𝜎2𝑣+𝜎2𝑎2+(𝑁2)𝜌2𝜎4𝑎(𝑁1)/4.(A.2) For a lower bound on |det(𝐑)|, we first represent the bounds on minimum and maximum eigenvalues (|𝜆𝑅,min|, |𝜆𝑅,max|) in terms of mean (𝜇𝜆𝑅) and variance (𝜎2𝜆𝑅) of eigenvalues along with the trace of the covariance matrix 𝐑 [22] as𝜇𝜆𝑅𝜎𝜆𝑅(𝑁1)1/2||𝜆𝑅,min||tr(𝐑𝐑)𝑁1/2,𝜇𝜆𝑅||𝜆𝑅,max||𝜇𝜆𝑅+𝜎𝜆𝑅(𝑁1)1/2,(A.3) where mean and variance are given by𝜇𝜆𝑅=tr(𝐑)𝑁=𝜎2𝑣+𝜎2𝑎,𝜎2𝜆𝑅=tr𝐑𝐑𝑁||𝜇𝜆𝑅||2=𝜎2𝑣+𝜎2𝑎2+(𝑁1)𝜌2𝜎4𝑎||𝜇𝜆𝑅||2=(𝑁1)𝜌2𝜎4𝑎.(A.4) Using (A.4) in (A.3), lower (𝑙) and upper (𝑢) bounds on minimum and maximum of eigenvalues are obtained as follows:||𝜆𝑙𝑅,min||=𝜎2𝑣+𝜎2𝑎(𝑁1)𝜌𝜎2𝑎||𝜆𝑅,min||𝜎2𝑣+𝜎2𝑎+(𝑁1)𝜌2𝜎4𝑎,𝑖1/2=||𝜆𝑢𝑅,min||,||𝜆𝑙𝑅,max||=𝜎2𝑣+𝜎2𝑎||𝜆𝑅,max||𝜎2𝑣+𝜎2𝑎+(𝑁1)𝜌𝜎2𝑎=||𝜆𝑢𝑅,max||.(A.5) Now, we use [23, Theorem 1] to find a lower bound on |det(𝐑)| as||det(𝐑)||||𝜆𝑙𝑅,min||Ω𝑅||𝜆𝑢𝑅,max||𝑁Ω𝑅,(A.6) where a parameter Ω𝑅 is given byΩ𝑅=𝑁||𝜆𝑢𝑅,max||||𝜇𝜆𝑅||||𝜆𝑢𝑅,max||||𝜆𝑙𝑅,min||.(A.7) An upper bound on the modulus of determinant of the nonprincipal submatrix can be represented in terms of eigenvalues (𝜆𝑅) of 𝐑1/2 [24] as𝛾𝑡=|||det𝐑1/2𝑗𝑖,𝑗𝑖|||12|||𝜆𝑢𝑅||||||𝜆𝑙𝑅||||||𝜆𝑢𝑅|||𝑁2.(A.8) Here, eigenvalues (𝜆𝑅) of 𝐑1/2 can be represented in terms of eigenvalues (𝜆𝑅) of 𝐑 using the eigenvalue decomposition as𝐑1/2=𝐐Λ1/2𝐐,(A.9) where 𝐐 is the matrix of eigenvectors and Λ is the diagonal matrix of eigenvalues of 𝐑. Using (A.8) and (A.9) and eigenvalue bounds of (A.5), an upper bound on 𝛾𝑡 can be obtained as𝛾𝑡=|||det𝐑1/2𝑗𝑖,𝑗𝑖|||12||𝜆𝑢𝑅,max||||𝜆𝑙𝑅,min||||𝜆𝑢𝑅,max||𝑁2.(A.10)

Acknowledgments

This work was carried out at IIT Delhi as part of the research project sponsored by Conexant Systems Inc. A few preliminary results of this paper were presented at the Annual IEEE Indicon Conference, INDICON 2009, Ahmedabad, India, in December 2009.