Crosstalk between subscriber loops used to transmit different broadband services through the same bundled telephone cable is a significant limitation to providing digital subscriber line services. A method for estimating the crosstalk, identifying the sources of the crosstalk and predicting additional sources of crosstalk are disclosed. The crosstalk sources are identified in the frequency domain by maximizing the correlation with a "basis set" of received crosstalk PSDs which consist of the cascade of a finite set of known transmit PSDs types times a representative set of crosstalk couplings. Multiple crosstalk types are identified with a technique of successive spectral subtraction. Once a type is identified the crosstalk disturber is compared against all other members of the complete set for that type. Additionally, Multiple Regression (MR) techniques and a Matching Pursuit (MP) algorithm are used to increase the ability of the system and method to identify various crosstalk disturbers. The identification of the crosstalk disturber is then used to perform crosstalk cancellation and spectrum management in a DSL system.

Claim:

We claim:

1. An apparatus, comprising: an input/output unit configured to receive data indicative of a power spectral density for a digital subscriber loop; and a processor configured to:correlate the received data indicative of the power spectral density for the digital subscriber loop with a predetermined set of power spectral densities corresponding to a set of possible crosstalk disturbers; select a crosstalk disturber from the setof possible crosstalk disturbers having a most closely correlated power spectral density; and refine the selected crosstalk disturber using a Best Basis Set technique.

2. The apparatus of claim 1, wherein the processor is further configured to: subtract the power spectral density for the refined crosstalk disturber from the data indicative of the power spectral density using spectral subtraction to generate aresidual power spectral density; correlate the residual power spectral density with the predetermined set of power spectral densities for the set of possible crosstalk disturbers; select an additional crosstalk disturber having the most closelycorrelated power spectral density to the residual power spectral density; and refine the selected additional crosstalk disturber using the Best Basis Set technique.

3. The apparatus of claim 2, wherein the processor is configured to subtract, correlate the residual power spectral density, select the additional crosstalk disturber, and refine the selection of the additional crosstalk disturber iterativelyuntil a stopping criterion is reached.

4. The apparatus of claim 3, wherein the stopping criterion is a predetermined correlation threshold.

5. The apparatus of claim 1, wherein the Best Basis Set technique is selected from the group consisting of the Matching Pursuit algorithm and the Focal Undetermined System Solver.

6. The apparatus of claim 1, wherein the processor is further configured to send the identity of the selected crosstalk disturber to a broadband spectrum management system.

7. The apparatus of claim 1, wherein the processor is further configured to send an estimated pair-to-pair coupling of the selected crosstalk disturber to a broadband spectrum management system.

8. An apparatus, comprising: means for receiving data indicative of a power spectral density for a digital subscriber loop; means for correlating the received data indicative of the power spectral density for the digital subscriber loop with apredetermined set of power spectral densities corresponding to a set of possible crosstalk disturbers; means for selecting a crosstalk disturber from the set of possible crosstalk disturbers having a most closely correlated power spectral density; andmeans for refining the selection of the crosstalk disturber using a Best Basis Set technique.

9. The apparatus of claim 8, further comprising: means for subtracting the power spectral density for the refined crosstalk disturber from the data indicative of a power spectral density using spectral subtraction to generate a residual powerspectral density; means for correlating the residual power spectral density with the predetermined set of power spectral densities for the set of possible crosstalk disturbers; means for selecting an additional crosstalk disturber having the mostclosely correlated power spectral density to the residual power spectral density; and means for refining the selection of the additional crosstalk disturber using the Best Basis Set technique.

10. The apparatus of claim 9, wherein the means for subtracting, the means for correlating the residual power spectral density, the means for selecting the additional crosstalk disturber, and the means for refining the selection of theadditional crosstalk disturber are configured to perform their respective functions iteratively until a stopping criterion is reached.

11. The apparatus of claim 10, wherein the stopping criterion is a predetermined correlation threshold.

12. The apparatus of claim 8, wherein the Best Basis Set technique is selected from the group consisting of the Matching Pursuit algorithm and the Focal Undetermined System Solver.

13. The apparatus of claim 8, further comprising: means for sending the identity of the selected crosstalk disturber to a broadband spectrum management system.

14. The apparatus of claim 8, further comprising: means for sending an estimated pair-to-pair coupling of the selected crosstalk disturber to a broadband spectrum management system.

15. A method comprising: measuring a power spectral density on a digital subscriber loop; using a processor to correlate the power spectral density for the digital subscriber loop with a predetermined set of power spectral densitiescorresponding to a set of possible crosstalk disturbers; selecting a crosstalk disturber, from the set of possible crosstalk disturbers, having the most closely correlated power spectral density; and refining the selected crosstalk disturber using aBest Basis Set technique.

16. The method of claim 15, further comprising: subtracting the power spectral density for the refined crosstalk disturber from the measured power spectral density of said digital subscriber loop using spectral subtraction to generate aresidual power spectral density; correlating the residual power spectral density with the predetermined set of power spectral densities for the set of possible crosstalk disturbers; selecting an additional crosstalk disturber having the most closelycorrelated power spectral density to the residual power spectral density; and refining the selected additional crosstalk disturber using the Best Basis Set technique.

17. The method of claim 16, wherein subtracting, correlating the residual power spectral density, selecting the additional crosstalk disturber, and refining the selection of the additional crosstalk disturber are iteratively performed until astopping criterion is reached, wherein the stopping criterion is a predetermined correlation threshold.

18. The method of claim 15, wherein the Best Basis Set technique is selected from the group consisting of the Matching Pursuit algorithm and the Focal Undetermined System Solver.

19. The method of claim 15, further comprising sending the identity of the selected crosstalk disturber to a broadband spectrum management system.

20. The method of claim 15, further comprising sending an estimated pair-to-pair coupling of the selected crosstalk disturber to a broadband spectrum management system.

Description:

FIELD OF THEINVENTION

This invention is related to the identification of crosstalk between subscriber loops transmitting high bandwidth telecommunication services through bundled telephone lines or cables, particularly to identification of crosstalk in digitalsubscriber line ("DSL") service loops. Identification of the crosstalk enables minimization of its effect through crosstalk cancellation, and as part of the spectrum management of the local plant in a telecommunications system. Identification of thecrosstalk disturbers is accomplished through a variety of methods including Multiple Regression (MR) analysis and an application of a Best Basis Set (BBS) problem with a sparse solution including use of a Matching Pursuit (MP) algorithm.

BACKGROUND

The mainstay of the telephone company local network is the local subscriber loop. The local subscriber loop is now being used to provide broadband digital telecommunication services such as DSL service. Such broadband DSL services includeintegrated services digital subscriber network ("ISDN"), high-rate digital subscriber line ("HDSL"), asymmetrical digital subscriber lines ("ADSL") and very high rate digital subscriber lines ("VDSL") technology. DSL services allow residential andbusiness customers to send and/or receive digital data at higher rates of speed than were previously possible using analog modem technology.

DSL technologies are engineered to operate over a class of subscriber loops, such as nonloaded loops (18 kft) or Carrier Serving Area (CSA) loops (9 to 12 kft). Digital Subscriber Line (DSL) technology exploits the existing, ubiquitous, coppertelephone loop plant to provide megabit per second (Mbps) high-speed Internet access and other services.

The great majority of residential customers and many business customers are served by such metallic (copper) twisted pair cables connected from a local switch in the central office ("CO") to the subscriber's land-line telephones. For eachsubscriber, telephone and DSL signals travel on a single twisted pair from a central office (CO) to the subscriber. Many, sometimes thousands, of twisted pairs are wrapped together in a single cable. The twisting of the pairs keeps the average amountof electromagnetic coupling between the balanced circuits on each pair to a low level, resulting in low crosstalk coupling between different circuits at voice frequencies. Twisted-pair cabling was designed and built to carry voice services, requiringonly a very low probability of intelligible crosstalk at voice frequencies, generally up to a few kHz.

Crosstalk generally increases, however, with increasing frequency, and because DSL frequencies extend into the MegaHertz ("MHz") range, crosstalk becomes the major limitation to high-speed DSL. Because a number of individual twisted pairs arewrapped together in a binder, and a number of binders are in each telephone cable (there are typically 12, 13, 25, 50 or 100 pairs in a binder) the crosstalk between pairs in the same binder is much higher than it is between pairs in different binders. Telephone cabling can be thought of as a multi-input, multi-output interference channel with crosstalk in-between from 12 to 100 channels, although there are typically a smaller number of high-power crosstalkers in each pair.

DSL technology is still fairly new. The general approach to date has been to treat each signal on each twisted pair entirely separately, and to assume that crosstalk is simply additive Gaussian noise of which nothing is known. With thisapproach, subscribers can enjoy relatively high-speed service, far faster than voice-band modems, up to multiple Megabits per second ("Mbps"). This is the current static approach to DSL spectrum management.

As time progresses it is expected that there will be many more DSL users each demanding higher speed service. This will result in more crosstalk and higher-bandwidth services that are more vulnerable to crosstalk. Without accurate systems andmethods for identifying crosstalk, canceling such crosstalk and thereby carefully managing the available frequency spectrum in local loop plants there may be no way to avoid problems that will ultimately result in customer disappointment, delay andhigher costs.

The American National Standards Institute (ANSI) Working Group T1E1.4, being responsible for DSL standards on ISDN Basic Access, HDSL and ADSL, initiated a project in 1998 to develop relevant standards for spectrum management. ANSI StandardT1.417-2001, Spectrum Management for Loop Transmission Systems, was approved and issued in January 2001. This ANSI standard applies relatively rigid rules uniformly across the entire loop plant. There have been a number of proposals for new DSLsystems, for example the class of symmetric DSL ("SDSL"). Unfortunately, this creativity complicates the problem of spectrum management, as is apparent from Table 1, which shows the nine currently agreed classes of DSL systems.

FIG. 8 shows currently adopted power spectral density ("PSD") templates for each of these spectrum management classes except Class 6 with line 110 representing Class 1 and line 180 representing Class 9 with the others represented respectively inbetween. To comply with the ANSI standard, the PSD of the signal transmitted by a piece of equipment must fall under the template at all frequencies. This will generally be verified during manufacture, but if problems arise in the field the localoperator will need to resolve incompatibilities and complaints, possibly including identification of systems that are transmitting in violation of the standard.

Current rules for DSL spectrum management assume static worst case crosstalk types and crosstalk couplings. These rules do not take into account the individual types of crosstalk sources or crosstalk couplings related to a particular pair in acable, which may be considerably different than the near worst-case couplings that are assumed in the ANSI standard. Dynamic spectrum management would take into account the individual types of crosstalk sources and crosstalk couplings of each particularcable, and could greatly increase the number of customers that can be provided DSL service and their bit rates. Thus, a system that can characterize crosstalk on a loop-by-loop basis has the potential to yield a better crosstalk characterization of theplant. Furthermore, the use of such crosstalk data in a loop spectrum management database, in turn has the potential to be mined, correlated and exploited to provide more optimal performances for individual subscriber loops.

It would be desirable to have a technique that could characterize the crosstalk environment on a loop-by-loop basis, in a mechanized and highly accurate manner without the need for special equipment or intervention at the subscriber's location.

Further it would be desirable to use the automated system for the identification of crosstalk for the spectrum management of DSL systems.

Therefore, it is desirable to have a means for identifying crosstalk between local loops that will enable a system operator to identify actual or potential crosstalk disturbers depending on the configuration of the DSL services being provided toa plurality of customers.

SUMMARY

In accordance with the present invention a method and system for the identification of crosstalk in high-bandwidth DSL systems are used to ensure spectral compatibility between a number of DSL lines/subscribers in a loop plant providing a meansfor spectrum management in a local loop plant.

A novel method for crosstalk identification in the frequency domain computes the correlation coefficient between the measured spectrum and a set of pre-defined PSDs that may be viewed as a crosstalk basis set. The size of the basis set may bemanaged through the use of a singular value decomposition of the complete basis set in accordance with the type of crosstalk disturber.

In the method of the present invention the measured power spectral density of a subscriber loop is correlated with a predetermined set of power spectral densities for a group of possible crosstalk disturbers and the crosstalk disturber havingthe most closely correlated power spectral density is selected as the identified disturber. Once identified, the power spectral density for the selected crosstalk disturber is subtracted from the measured power spectral density of said subscriber loopusing spectral subtraction and the resulting power spectral density is again correlated with the predetermined set of power spectral densities for the group of possible crosstalk disturbers to identify additional disturbers.

The method may also employ Multiple Regression (MR) techniques or the Matching Pursuit (MP) algorithm for finding an optimal sparse representation of a vector from an overcomplete set of vectors. These techniques can be employed to increase theaccuracy of the crosstalk disturber identification method using correlation or can be used in place of the correlation method.

The method may be applied to a DMT receiver capable of receiving a signal, r.sub.i(k) and capable of generating a set of decisions {circumflex over (D)}.sub.i(n), from decision variables Q.sub.i(n), and channel estimates H(n). A data slicerextracts the set of decisions {circumflex over (D)}.sub.i(n), from decision variables Q.sub.i(n). A crosstalk PSD estimator estimates the power spectral density of the crosstalk disturber based on the set of decisions {circumflex over (D)}.sub.i(n), thedecision variables Q.sub.i(n) and the channel estimates H(n) and a crosstalk identifier identifies the type of the crosstalk disturber by correlating the estimated power spectral density of the crosstalk disturber with a set of possible disturbers andselecting the disturber exhibiting the maximum correlation. Furthermore, a crosstalk predictor may be added to predict the p-th order linear MMSE prediction of the crosstalk disturber. The prediction may be used to cancel the effect of the crosstalkdisturber in the DMT receiver.

The identity of crosstalk disturbers can also be sent to an Operations Support System (OSS) for the management of the broadband spectrum available for a plurality of subscriber loops. The information could also be used to cancel the effect ofthe crosstalk disturber in a DSL receiver.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 depicts a DSL subscriber loop system having two DSL subscriber loops connecting to a central office;

FIG. 2 depicts an illustrative embodiment of a crosstalk identification system according to the present invention;

FIG. 3 depicts a cable binder having N twisted pairs of copper wire;

FIG. 4 depicts the power sum of NEXT loss in dB for a binder having 25 twisted pairs of copper wire;

FIG. 5 depicts examples of measured pair-to-pair NEXT loss;

FIG. 6 depicts an example of a fourteen member canonical set of pair-to-pair NEXT couplings;

FIG. 7 depicts an example of a fourteen member basis set of T1 NEXT couplings;

FIG. 8 depicts the power spectral density ("PSD") templates for each of the spectrum management classes;

FIG. 9 depicts the NEXT PSD for a first unknown disturber;

FIG. 10 depicts the NEXT PSD for a second unknown disturber;

FIG. 11 is a flow diagram of the crosstalk identification and subtraction algorithms of the present invention.

FIG. 13 is a flow diagram of the crosstalk identification algorithm using Multiple Regression (MR) and the BBS sparse solution technique with both a reduced and original crosstalk basis set.

FIG. 14 is a flow diagram of the crosstalk identification algorithm incorporating the use of Multiple Regression (MR) or the BBS sparse solution techniques such as Matching Pursuit.

DETAILED DESCRIPTION

FIG. 1 depicts crosstalk in an example of a local loop environment. DSL loops 18 comprising twisted pair copper wire connect two subscribers 20 to a central office ("CO") 12 or remote terminal ("RT"). Near-end crosstalk ("NEXT") 30 between thetwo loops may be present at the end of loops where the disturber's transmitter is located (either at the CO or the RT) while far-end crosstalk ("FEXT") 40 may be present at the end of the loop 18 opposite to where the disturber's transmitter is located.

Twisted-pair wire 18 can be considered a multi-user channel by expressing the received signal vector as Y(f)=R(f)+N(f)+H(f)D(f) (1) where R(f) is the received message vector, N(f) is the vector of independent background noise, D(f) is the vectorof transmitted signals creating crosstalk on nearby pairs, and H(f) is the matrix of crosstalk couplings. In the case of near-end crosstalk (NEXT), R(f) is transmitted from one end of the cable and D(f) is transmitted from the other end. The receivedcrosstalk H(f)D(f) can be canceled, minimized or otherwise treated so as to be less interfering than white noise.

The twisted-pair multi-user channel is generally time invariant. Crosstalk couplings may vary as the temperature of the cable changes, but only very slightly. The crosstalk sources may turn on and off, but not very often because the mainapplication of broadband service is "always on" Internet access. Time invariance allows the possibility of very accurate estimates of crosstalk coupling using large sample sizes. Joint statistics can be obtained by a "third-party" operating a thirdparty system, which is connected to all lines through communications links to the modems, and the knowledge of all the transmitted data could be used to estimate the crosstalk couplings. A "third-party" or operations center having access to multi-usercrosstalk statistics would be far more capable of diagnosing crosstalk problems than an uninformed technician.

Received crosstalk can be measured by temporarily connecting measurement equipment to a single p'air while crosstalk couplings could be measured by connecting measurement equipment to multiple pairs. The accuracy of estimates depends stronglyon which elements of the multi-pair cable may be accessed, for example crosstalk couplings can be estimated accurately by accessing data from both transmitters and receivers. Crosstalk parameter estimation can be performed off-line using very manysamples of data. Knowledge of crosstalk parameters can be useful for spectrum management and maintenance purposes and can also be used at run time to enhance modem performance.

There are many vagaries in the telephone outside plant, some older cables have poor crosstalk performance. Unanticipated problems from crosstalk are likely to occur and may be extremely difficult to diagnose.

The method of the present invention can be performed in any single receiver. The first step is identification of the type of DSL line that is generating the crosstalk. To accomplish this a receiver estimates which varieties of DSL aregenerating crosstalk. This can be accomplished because there is actually a small number of common crosstalker types at high frequencies: T1 lines, integrated services digital subscriber lines (ISDN), high bit-rate digital subscriber lines (HDSL),single-pair high bit-rate digital subscriber lines (HDSL2), asymmetric digital subscriber lines (ADSL), and symmetric digital subscriber lines (SDSL). Single-pair high-speed digital subscriber lines (G.shdsl) and very-high rate digital subscriber lines(VDSL) are also likely to become possible sources of crosstalk soon. Each of these crosstalker types has a unique crosstalk power spectral density (PSD), and this PSD is estimated with classic matched-filter correlation techniques. The basicimplementation of the method is described by steps 310, 320 and 330 of FIG. 11, which also shows the additional steps in the iterative process described below.

Referring to FIG. 2, a crosstalk identification system according to the present invention provides for a computer 200 having a processor 210 in communication with a memory 220 (which can comprise either primary memory, such as RAM, or secondarymemory, such as a magnetic disk drive or other storage media or a combination thereof) and input/output (I/O) unit 230. I/O unit 230 is adapted to receive data indicative of the power spectral density of one or more subscriber loops 18. Display 260 maybe used to display graphical and/or textual information related to the PSD measurements and crosstalk identification. Other input/output devices (not shown) may be used in conjunction with the computer 200 such as a keyboard, mouse, touchpad, trackballetc. in order to provide a user interface.

Processor 210 executes the method steps described herein as stored in memory 220 in order to identify the crosstalk on one or more subscriber lines 18. Additional software may also be executed to generate the PSD from the measurements inputthrough I/O unit 230. Alternatively, a separate test measurement device capable of measuring the PSD for the various subscriber loops may be attached to I/O unit 230.

The identification algorithm is extended to the case of multiple crosstalkers by a process of successive estimation and cancellation using spectral subtraction described in greater detail below.

Knowledge of the type of crosstalkers could be used for crosstalk cancellation in a DSL receiver. Once the type of the crosstalk source is determined, then the correlations of the source of the crosstalk are known, assuming correct decisions. These correlations are used to construct an optimal set of weightings for a type of predictive filtering that minimizes the received crosstalk power before detection in the receiver. The coefficients of the predictive filter are predetermined based onthe known properties of each of the types of crosstalk.

Sources of crosstalk are often called "disturbers." There are two types of crosstalk: near-end crosstalk (NEXT) and far-end crosstalk (FEXT). NEXT is more powerful than FEXT, particularly below about 1 MHz where many DSLs use overlappingspectra. If there is one crosstalker, then the received crosstalk PSD is the product of a PSD transmitted on a nearby pair times the crosstalk coupling transfer function. With multiple crosstalkers the received crosstalk PSD is the linear sum of eachcomponent. FIG. 3 depicts a cable binder 100 having N twisted pairs of copper wire 18. Lines 40 of FIG. 3 represent the NEXT between the first two twisted pairs of copper wire 18 while lines 41 represent the NEXT between the first and last pairs ofcopper wire 18.

Crosstalk can be characterized in terms of power-sums. For a given pair in a binder group, the power-sum NEXT is formed from the sum of the pair-to-pair NEXT coupling powers of the other pairs in the binder group into this given pair. For atypical 25-pair binder group, the 25 power sums are shown in FIG. 4. The power-sum is usually displayed as a power-sum "loss" and the lower the loss the higher the NEXT coupling. The NEXT power-sum loss is approximately linear with frequency on thelog-log scale. The NEXT model often used for studies in the industry is stated as expected 1% worst case power-sum crosstalk loss as a function of frequency. This means that on average, 1% of the pairs tested are subjected to power-sum crosstalk lossless than that assumed by the model at the given frequency. Such a model is a smooth curve versus frequency, in which the loss decreases at about 15 dB per decade of frequency. This statistical model ensures that the 1% worst case NEXT power sum levelfollows the curve over frequency but this noise level may appear on different pairs at different frequencies. The 24-disturber NEXT power level for one specific pair is not a smooth curve if plotted over frequency and the deviation from the statisticalmodel increases as the number of interferers used in the power sum is reduced. The power-sum loss for less than a full binder group depends on the distribution of the pairs on which the crosstalking signal appears.

For a 25-pair binder group there are 300 distinct pair-to-pair couplings (assuming symmetry). Examples of six measured pair-to-pair NEXT couplings from a binder group of twisted pair wires 18 are depicted in FIG. 5. It should be noted thatthere is wide variation of coupling strength for each pair-to-pair combination with frequency. The power summing process smoothes and averages some of this variation.

The received signal on pair number k can be expressed as

.function..function..function..noteq..times..function..times..function. ##EQU00001## where R.sub.k(f) is the received message signal, N.sub.k(f) is independent background noise, D.sub.i(f) is the signal transmitted on pair i, H.sub.ik(f) is thecrosstalk coupling between pair i and pair k, and the crosstalk received on pair k from a system transmitting on pair i is H.sub.ik(f)D.sub.i(f). Crosstalk has a number of elements which may be identified or estimated: the crosstalk coupling H.sub.ik(f)between each pair, the crosstalkers' transmitted spectrum D.sub.i(f), sequences of sampled received crosstalk, etc. DSL modems can measure crosstalk samples. Assuming correct decisions, the received sampled sum of crosstalk and noise is easily obtainedin each receiver by subtracting R.sub.k(f)={circumflex over (R)}.sub.k(f) from Y.sub.k(f). These samples and their statistics can be exploited to improve the performance of the receiver. Often a process of successive estimation and cancellation isperformed where the most powerful crosstalk disturber is first estimated and subtracted, then the next most powerful and so on.

If crosstalk is measured while a DSL system is running, what is measured or estimated is the crosstalk in the bandwidth of the considered DSL system. For this reason, other DSL services running on adjacent pairs may not be detected if theirbandwidth is not significantly overlapping, with the bandwidth of the disturbed system. A typical example would be ISDN crosstalk into non-overlapped downstream ADSL. If this is the case, there is little performance degradation for the considered DSLservice but, from a spectrum management point of view, it may be important to have an accurate map of all the services that generate crosstalk into a given pair. Moreover, it may also be important to identify the services that are generating crosstalkon a pair that may not even be carrying DSL services. An accurate map of all the services that generate crosstalk into a given pair is an impossible goal to achieve via modem-based methods of system identification if the considered pair is not in use. In this case, the way to achieve such a complete description of the crosstalkers is to perform a preliminary PSD measurement, e.g. via a spectrum analyzer or a selective voltage meter.

An important aspect of the crosstalk identification algorithm is a calculation of crosstalk correlation in the frequency domain. The correlation receiver is known to be an optimum detector in the case of equally likely transmitted signals withknown waveforms transmitted over an AWGN channel. It is mathematically equivalent to the matched-filter receiver which can be shown to maximize the output signal to noise ratio by using Schwarz's inequality. The basic structure of a correlationreceiver is used here, except that it is adapted with known transmitted signals but unknown crosstalk couplings

.rho..sigma..sigma..times..times..times..function..times..times..mu..mu..- times..times..times..times..times..ltoreq..rho..ltoreq. ##EQU00002## The correlation coefficient is a measure of how well the two data sets move together (positivevalues) or move apart from one another (negative values).

The data array Y={y.sub.1, . . . , y.sub.N} is the measured crosstalk power spectrum density (PSD) of the unknown DSL disturber. The data array X={x.sub.1, . . . , x.sub.N} is a reference or basis crosstalk PSD profile caused by a known DSLdisturber. Each basis set of crosstalk PSDs is generated from a single canonical set of measured strong pair-to-pair crosstalk couplings and a specific type of transmitted DSL PSD. Each type of transmitted DSL PSD is multiplied by all of the canonicalcrosstalk couplings to generate a basis set for that DSL. The assumption is that a member of this basis crosstalk set is highly correlated with any crosstalk that is strong enough to be troublesome or measurable. An example of a canonical set ofpair-to-pair couplings is shown in FIG. 6. This canonical set comprises the fourteen out of 300 possible pair-to-pair NEXT couplings with the smallest crosstalk loss sum. The crosstalk loss sum is obtained by dB summing over 401 equally spacedfrequency samples in the 10 kHz-2 MHz band.

The basis set of T1 disturber NEXT PSDs for this example is shown in FIG. 7. Basis sets for ISDN, HDSL, ADSL, 2B1Q SDSL and HDSL2 can be generated similarly. Increasing the number of members belonging to the canonical set of crosstalk transferfunctions would increase the accuracy of the identification algorithm. Another example of a canonical set of pair-to-pair couplings could be all 300 measured independent pair-to-pair couplings. A particularly useful NEXT basis set is found to comprisethe ISDN, HDSL, ADSL, 2B1.OMEGA. SDSL and HDSL2 NEXT sets that have power exceeding or equaling the 49.sup.th disturber power for a given DSL type obtained from the ANSI 1% NEXT model as defined in the Spectrum Management Standard, T1.417.

There are several ways of determining the canonical crosstalk PSD basis set. The basis set .XI..varies..left brkt-bot.X.sup.(1) X.sup.(2) . . . X.sup.(P).right brkt-bot. is the full rank N.times.P matrix containing all the canonical crosstalkPSDs. Vector X.sup.(i).varies.[X.sub.1.sup.(i) X.sub.2.sup.i(i) . . . X.sub.2.sup.(i)].sup.T is the i-th crosstalk PSD belonging to the canonical basis and X.sub.l.sup.(i) is its sample at frequency f.sub.I. Vector a.varies.[a.sub.1 a.sub.2 . . .a_].sup.T contains weighting coefficients and E.varies.[e.sub.1 e.sub.2 . . . e.sub.N].sup.T is an error vector. The presence of the error vector E is due to the fact that the measured crosstalk Y is due to a pair-to-pair coupling function that may beclose, but not equal to, a coupling in the canonical basis.

Identification of the crosstalk disturbers is generally improved with increasing dictionary size. Large dictionaries, however, imply searches with high computational complexity. The chosen size for the canonical set is a trade-off betweenaccuracy and computational complexity. As described above, a canonical PSD basis set can be formed by using all the available measured pair-to-pair crosstalk couplings and combining them with all the DSL types under consideration. This approach,however, yields a very large basis set. To reduce the size of this large basis set, the crosstalk power within each DSL type can be ranked with only the first B vectors with the smallest crosstalk loss sum being chosen. The crosstalk loss sum isobtained by dB summing over the frequency samples in the band of interest.

The similarity of crosstalkers of the same type (but different pair-to-pair couplings) is much greater than the similarity of crosstalkers of different types. The present invention exploits this observation by applying the following method. For each type of crosstalkers, the dominant dimensions from the space spanned by all the PSD's for that type are extracted. Vectors corresponding to these dimensions are then used in lieu of the complete set of crosstalkers for that type. After thetype of crosstalker is identified, the search for pair-to-pair couplings in the basis that is closed to the actual one is done within the subset of the original (and large) basis pertaining to the identified crosstalker type.

To obtain the above-described vectors, the set or matrices of all possible PSD couplings .XI. is partitioned into A sub-matrices pertaining to each type of crosstalker or disturber, j. Performing a Singular value Decomposition (SVD) on each ofthese sub-matrices results in the following equation: .XI..sup.(j)=U.sup.(j)D.sup.(j)V.sup.(j)', 1.ltoreq.j.ltoreq.A, (4) where U.sup.(j) and V.sup.(j) are orthogonal and D.sup.(j) is a diagonal matrix containing the singular values. Given the number ofdimensions q.sub.j that are desired to be used to approximate .XI..sup.(j), the columns of U.sup.(j) are extracted which correspond to the q.sub.j largest singular values. These columns are then used as the basis for the crosstalkers of the given type. Note that this new basis has

.times. ##EQU00003## vectors and is much smaller than the original one (i.e., {tilde over (P)}<<P). Moreover, all vectors in each U.sup.(j) are orthogonal, thus mitigating any multicollinearity issues associated with the old basis.

This new and reduced-size basis enables us to identify efficiently which type(s) of crosstalk exist in the sample, but not which particular pair-to-pair couplings are associated with the crosstalk. However, after identifying the crosstalktypes, we can rerun the analysis using the subset of the original crosstalk basis associated with that type of crosstalker. Practically, we have decomposed the original search problem among P=A.times.B vectors, into two search problems, the first oneamong {tilde over (P)}<<B vectors, the second one among B vectors

In the method of the present invention for the case of a single disturber, the measured crosstalk PSD is correlated as described in the last subsection with each of the basis PSDs, and the identification is simply the type of crosstalker whichhas the highest correlation. In the examples set forth herein, most common DSL types are simulated: ISDN Basic Rate Interface (BRI), HDSL, T1, ADSL, 400 kbps SDSL, 1040 kbps SDSL, 1552 kbps SDSL, and HDSL2. All transmit system PSDs are as defined inthe Spectrum Management Standard, T1.417. ADSL and HDSL2 have different PSDs upstream and downstream.

An unknown disturber's NEXT PSD is shown in FIG. 9. The calculated correlation with each member of each DSL crosstalk basis set (for the example of a fourteen member canonical crosstalk coupling set) is given in Table 2 set forth below. TheDSL with the highest correlation and thus the identity of the unknown disturber is correctly identified as 1552 kbps SDSL.

Another unknown disturber's NEXT PSD is shown in FIG. 10. The calculated correlation with each member of each DSL crosstalk basis set (for the example of a fourteen member canonical crosstalk coupling set) is given in Table 3. The DSL with thehighest correlation and thus the identity of this unknown disturber is correctly identified as downstream ADSL.

In realistic cases it has been found that the DSL type of a single disturber with significant power is identified correctly nearly every time.

The case of mixed crosstalk arises when at least two different services are present in a binder. If the measured crosstalk PSD consists of the contribution of different kinds of disturbers, the identification algorithm described above would notidentify the whole set of crosstalkers all at once but would identify the strongest candidate only.

In order to overcome this limitation in cases of mixed crosstalk, the method of the present invention uses a frequency-domain onion peeling technique based on spectral subtraction methods as set forth in FIG. 11. Spectral subtraction is amethod originally proposed in speech and music processing for the restoration of the power spectrum of a signal observed in additive noise through subtraction of an estimate of the average noise spectrum from the noisy signal spectrum.

At the start 300 of the method the iteration counter is set to zero 305 and at step 310 a measurement of the PSD of a subscriber loop 18 is taken. Next at step 320 the first estimated PSD is set to the measured PSD as set forth below: |{tildeover (C)}.sub.1(f)|.sup.2.varies.|C(f)|.sup.2 (5)

At step 325, the iteration counter is set at one and the estimated PSD of the received crosstalk is correlated at step 330 with the basis set to identify the first source of crosstalk, i.e., the one having the highest correlation in accordancewith the procedure described above.

The onion-peeling technique used to identify additional crosstalk disturbers can be explained by focusing on the crosstalk received on a single pair. The subscript k in equation (2) is dropped andy(m)=r(m)+n(m)+.SIGMA..sub.ih.sub.i(m)d.sub.i(m) is the noisy received signal in the sampled-time domain, where r(m) is the received message signal, n(m) is additive background noise, and y(m) is the noisy observation, respectively. To simplify further,define each crosstalk component as c.sub.i(m)=h.sub.i(m)d.sub.i(m) so that the composite received crosstalk is c(m)=.SIGMA..sub.ic.sub.i(m)=.SIGMA..sub.ih.sub.i(m)d.sub.i(m). Then ignoring the background noise and transforming to the frequency domain,Y(f)=R(f)+C(f). Power spectrum subtraction is defined as in the following: |{circumflex over (R)}(f)|.sup.2=|Y(f)|.sup.2-b |C(f)|.sup.2 (6) where |{circumflex over (R)}(f)|.sup.2 is an estimate of the signal PSD, |Y(f)|.sup.2 is the instantaneous PSD ofthe noisy observation, b is a coefficient chosen to minimize the residual error and |C(f)|.sup.2 is the time-averaged noise spectra. The basic assumption is that the noise is a slowly varying process so that averaging would not smear the noise spectrum. Let c(m)=c.sub.1(m)+c.sub.2(m) be a composite crosstalk constituted of two crosstalk terms pertaining to two different DSL systems. Assuming that the identification algorithm (described above) identifies the strongest disturber that generated crosstalkc.sub.1(m), an estimate of its PSD |C.sub.1(f)|.sup.2 can be obtained. At step 340 a decision is made as to whether the maximum correlation is greater than a set threshold. If it is not, then the algorithm has identified all of the crosstalkcomponents. If it is, the identified component must be subtracted for further processing. Basically, the composite crosstalk c(m) can be viewed as a noisy observation of the useful signal c.sub.2(m) embedded in the noise c.sub.1(m). Because the outputof the crosstalk identification algorithm is a true PSD obtained on the basis of a PSD mask and an estimated pair-to-pair coupling function, there is no need to perform an averaging operation on the estimated PSD of the strongest crosstalker (noise)|C.sub.1(f)|.sup.2. Moreover, the instantaneous PSD of the noisy observation |Y(f)|.sup.2 (necessary in the original formulation of the problem because of the non-stationary nature of speech and music signals) can be replaced in our case by the measuredPSD |C(f)|.sup.2. Therefore, the power spectrum subtraction for crosstalk identification can be defined as |{tilde over (C)}.sub.2(f)|.sup.2=|C(f)|.sup.2-|C.sub.1(f)|.sup.2, (7) yielding the residual PSD |{tilde over (C)}.sub.2(f)|.sup.2 at step 350 ofFIG. 11.

Due to the great variations in the pair-to-pair coupling functions, spectral subtraction can result in negative estimates of the power spectrum. Power spectra are non-negative functions of frequency, and any negative estimate of these variablesshould be mapped into a non-negative value. To avoid negative estimates, the power spectrum is post processed in steps 360 and 370 using a mapping function T[ ] of the form:

Empirically a good value of .beta. ranged from approximately 0.0 to approximately 0.1. The simplest choice for the function .PHI.[.cndot.] is .PHI.[|C(f)|.sup.2]=noise floor=-140 dBm/Hz. Another possible choice could be.PHI.[|C(f)|.sup.2]=.beta.|C(f)|.sup.2.

A flow chart for implementing the crosstalk identification algorithm is given in FIG. 11. The residual PSD |{tilde over (C)}.sub.2(f)|.sup.2 is fed to the crosstalk identification algorithm that is now able to unveil the nature of the seconddisturber and estimate its PSD |C.sub.2(f)|.sup.2. Steps 325 and 330 are reiterated until all disturbers are identified, i.e., the iteration ceases at step 380 when the maximum correlation is less than or equal to the threshold at step 340. It has beenfound by simulations that overall accuracy is enhanced by retaining and using identified crosstalkers only if the maximum crosstalk correlation is greater than a certain threshold that is empirically determined. The threshold is preferably betweenapproximately 0.7 and 0.99 and is more preferably approximately 0.9. An example of the technique follows.

Consider the measured mixed crosstalk given in FIG. 12 that happens to be a mixture of crosstalk from SDSL @1040 kbps and HDSL. The crosstalk identification algorithm yields the results given in Table 4.

The crosstalk identification method of the present invention as set forth in FIG. 11 first identifies the strongest disturber as SDSL @1040 kbps. On the next iteration it identifies HDSL as the second strongest disturber. Finally on the thirditeration, the maximum correlation of the remainder NEXT is less than 0.90, the set threshold.

The use of spectral subtraction may cause deterioration in the quality and the information content of the residual crosstalk PSD. In the literature, there are a number of variants of spectral subtraction that aim to provide consistentperformance improvement across a range of SNRs.

The spectral subtraction equation can be expressed as the product of the noisy signal spectrum and the frequency response of a spectral subtraction filter as: |{circumflex over (R)}(f)|.sup.2=|Y(f)|.sup.2- |C(f)|.sup.2=H(f)|Y(f)|.sup.2 (9)assuming for now that b=0, where the frequency response of the spectral subtraction filter is defined as:

.function..function..function. ##EQU00005##

The spectral subtraction filter H(f) is a zero-phase filter, with its magnitude response in the range 0.ltoreq.H(f).ltoreq.1. The filter acts as an SNR-dependent attenuator. The attenuation at each frequency increases with the decreasing SNR,and conversely decreases with increasing SNR.

The least mean squared error linear filter for noise removal is the Wiener filter. The implementation of the Wiener filter requires the power spectra (or equivalently the correlation functions) of the signal and the noise process. The Wienerfilter is based on ensemble average spectra of the signal and the noise, and the averaging operations are taken across the ensemble of different realizations of the signal and noise processes. On the other hand, the spectral subtraction filter, usesinstantaneous spectra of the noisy observation and the time averaged spectra of the noise. This is necessary since in spectral subtraction only one realization of the noise process is available.

Spectral subtraction is computationally efficient because it only requires an estimate of the noise power spectrum. Because it does not utilize the statistics of the signal process, however, it may use too little a priori information. However,for an ergodic process, the time-averaged spectrum approaches the ensemble averaged spectrum, so that the spectral subtraction filter asymptotically approaches the Wiener filter. This property does not hold when processing speech and music signalsbecause they are intrinsically non-stationary (and, therefore, non-ergodic). When spectral subtraction is applied to the present case of crosstalk identification, however, the quantities involved are true PSDs and not instantaneous or time-averagedones. Therefore, the implementation of spectral subtraction may be viewed as a sort of optimal Wiener filtering. This interpretation may explain why spectral subtraction seems to perform very well when applied to successive crosstalk identification. The major limitation encountered is the non-linear processing distortion due to the necessity of mapping negative values into non-negative values. This distortion increases at every iteration when equation (8) is implemented and limits the accuracy ofthe method to the identification of up to a small number of different crosstalker types. From a practical point of view, this does not look as a major limitation because it is very unlikely that more than a few interferers cause performance degradationin any one DSL system.

There are methods other than correlation that may be used to identify the source of crosstalk including the use of multiple regression techniques and Matching Pursuit (MP) algorithms. Such methods would replace the identification of crosstalkdisturber by correlation or could be used to refine the correlation method. Two other approaches to the problem of multiple disturber identification as described below, enhance the identification accuracy, thus by-passing the non-linear distortionproblems due to the spectral subtraction technique.

Linear Multiple Repression Approach

Let Y={y.sub.i}.sub.i=1.sup.i=N be the set of frequency samples of the measured crosstalk PSD caused by an unknown DSL disturber. Similarly, let X={x.sub.i}.sub.i=1.sup.i=N be the set of frequency samples of the basis crosstalk PSD profilecaused by a known disturber. Let P be the number of all the basis crosstalk PSD profiles, so that the set of frequency samples of the k-th basis crosstalk PSD profile is X.sup.(k)={x.sub.i.sup.(k)}.sub.i=1.sup.i=N, where 1.ltoreq.k.ltoreq.P. Theproblem is to find the single disturber that generates crosstalk Y given the set of all the crosstalk PSD profiles X.sup.(k)={x.sub.i.sup.(k)}.sub.i=1.sup.i=N, i.e. to determine the relationship between the two variables K.sup.(k) and Y. This is aclassical statistical problem that is solved with linear regression. In particular, the linear regression model is: y.sub.i=a.sup.(k)+b.sup.(k)x.sub.i.sup.(k)+e.sub.i.sup.(k), 1.ltoreq.i.ltoreq.N, 1.ltoreq.k.ltoreq.P (11) where the regressioncoefficients a.sup.(k) and b.sup.(k) are determined by the condition that the sum

.times..times..times. ##EQU00006## of the squared residuals is minimum. The fitting of the regression coefficients is done over the N frequency samples. In so doing, P residuals are obtained and the DSL with the lowest S.sup.(k) is identifiedas the disturber.

It is easy to show that the sum of squared residuals S.sup.(k) can be expressed in terms of the correlation coefficient in (3). In fact, plugging into S.sup.(k) the expression of the optimal regression coefficients, after some manipulation itcan be shown that:

.rho..times..mu. ##EQU00007##

Now, since the term

.times..mu. ##EQU00008## is constant, the search for the minimum sum of squared residuals S.sup.(k) is equivalent to the search of the highest correlation coefficient .rho..sub.X.sub.(k).sub.,Y. Therefore, the correlation method described forthe single disturber case is.sub.y equivalent to solving a single linear regression problem.

While standard techniques are available to solve such multiple regression problems the application of multiple regression techniques requires the use of several non-standard limitations. All non-zero coefficients must be positive, i.e. thepower contribution to the spectral density by each disturber is positive. Moreover, there are many more potential disturber spectra than frequencies at which the power is evaluated, whereas standard regression techniques require that the number ofindependent variables be strictly less than the number of observations. Independent variables (the disturbers) may be disqualified from the model in pre-specified groups based upon a group's total power contribution to the modeled spectra. The problemof having more independent variables than observations can be solved by using the Singular Value Decomposition method described in the present invention that allows a dramatic reduction of the size of the crosstalk basis set. A way of checking that thecoefficients remain positive may be efficiently addressed by computing the so called t- or f-statistics for every coefficient, so that we can eliminate a certain number of columns at every iteration on the basis of the variance of their estimation. These quantities cal also be used for the stopping criterion. One of the main advantages of using multiple regression techniques is the fact that no intermediate residual needs to be computed, thus alleviating the non-linear distortion effects due tothe calculation of the residuals.

Matching Pursuit Approach

Rewriting equation (13) with M=P, i.e. considering all the PSD profiles of the basis, and by posing a.sub.0.sup.(k.sup.1.sup., . . . , k.sup.M.sup.)=0 it can be seen that: y.sub.i=a.sub.1x.sub.i.sup.(1)+a.sub.2x.sub.i.sup.(2)+ . . .+a.sub.Px.sub.i.sup.(P)+e.sub.i, (14)

In the majority of the cases of practical interest for crosstalk identification, it is the case that N<P. For example, the inventors have previously considered cases for which N=401 and P ranges between 80 and 2,560 (obviously, higher valuesof P yield lower misidentification rates). Since N<P, matrix .XI. constitutes an overcomplete set of vectors for .sup.N. Moreover, in practice, only combinations of at most 10 crosstalkers can be considered measurable (or harmful) within a binder. Therefore, among the P weighting coefficients a.sub.j, 1.ltoreq.j.ltoreq.P, only few should be expected to be non-zero.

On the basis of the previous considerations, the problem of identification of multiple crosstalkers can be considered equivalent to the problem of finding an optimal sparse representation of a vector from an overcomplete set of vectors. Unfortunately, it has been recently shown that the optimal solution to this problem is an NP-complete problem and requires a combinatorial search of prohibitive cost. Thus, more practical but suboptimal vector selection algorithms have been developed. Among these algorithms, those suitable for the problem of multiple crosstalker identification are the Matching Pursuit (MP) algorithm and its variations, and the FOCal Undetermined System Solver (FOCUSS). Given a signal vector Y and an overcompletebasis .XI., these algorithms solve the problem of finding the most compact representation of Y within a given tolerance using the basis vectors in the dictionary .XI..

In order to apply these techniques to our problem, several modifications to the original algorithms must be implemented. In particular, in addition to the non-standard limitations mentioned above, some or all of the following modificationsshould be considered. At each iteration, the difference between the measured data and its projection along the direction of the chosen columns must be carried out. This difference must always be positive. When a column is chosen, it is added to thedictionary .XI. and no residual is computed. When a column is chosen, the whole dictionary is projected along that direction and no residual is computed.

In FIG. 13 the use of the MR or BBS techniques in a method of identifying an crosstalk disturber is depicted. In FIG. 13 at step 410 the power spectral density (PSD) for a given line is measured. At step 420 the basis set of possible crosstalkdisturbers is reduced using single value decomposition. The reduced crosstalk basis set is then used to identify, at step 430, the of the identity of the crosstalk disturbers using either the MR or BBS technique described above. At step 440 the actualpair-to-pair coupling of the crosstalk disturber is then estimated using the MR or BBS techniques but utilizing the original complete basis set pertaining to the identified disturber.

It is also possible to combine the correlation method with the two methods of crosstalk identification previously described. This can be done by using the correlation method for an initial estimate of the first crosstalker, and then refiningsuch estimate using the MR and BBS techniques. The correlation method would identify the type of crosstalker, and the MR or BBS algorithms would exploit this information to refine the estimate of the crosstalker. This combination of methods has provento be very effective in the construction of good residuals so that the overall accuracy of the identification is enhanced. An example on how to combine these methods is given in FIG. 14. In FIG. 14 at step 510 the power spectral density (PSD) for agiven line is measured. At step 520 an iteration counter is set and at step 530 the identity of the crosstalk disturber is first estimated using the correlation technique described above. At step 535 the estimated crosstalk is refined using the MR orBBS techniques. At step 540 the stopping criterion is applied and is the identified crosstalk disturber falls within the stopping criterion then the method ends at step 580. If the stopping criterion has not been met the residual is computed usingspectral subtraction at step 550 as described above and the iteration counter is increased by a count of one at step 525. The next most likely crosstalk disturber is then estimated in the next iteration which process continues until the stoppingcriterion is met.

The above description has been presented only to illustrate and describe the invention. It is not intended to be exhaustive or to limit the invention to any precise form disclosed. Many modifications and variations are possible in light of theabove teaching. The applications described were chosen and described in order to best explain the principles of the invention and its practical application to enable others skilled in the art to best utilize the invention on various applications andwith various modifications as are suited to the particular use contemplated. For example, it is known that the measurement of the power spectral density on a subscriber line could either be a direct measurement taken through a network analyzer or otherequivalent device but could also include estimated measurements based on information provided to the system from a modem.