Abstract:
The present invention comprises a method of obtaining coefficients for a transmit filter such that the power spectral density of the output for the different frequencies comes close to but does not exceed a maximum power spectral density of a communication(s) standard. By first doing a convex optimization procedure to obtain the autocorrelation coefficients for the filter and then using the autocorrelation coefficients to determine the filter coefficients, a low-order filter that closely approximates the desired output power spectral densities can be produced.

Description:
BACKGROUND OF THE INVENTION 
     The present invention relates to a transmit filter for high-bit-rate digital subscriber line (HDSL2 systems). The HDSL2 standard (currently, Draft HDSL2 Standard, T1E1. 4/2000-006, ANSI) requires both the upstream and the downstream transmitted power spectrum densities (PSDs) to be below a certain level. Such transmitted power specifications are standard on present and future communication systems. It is desired to have a transmit power that comes close to but does not exceed the maximum allowable power at any frequency. Because the specification is relatively complex, it is difficult to come up with a low-order filter that allows the system to approach the maximum allowed power output. 
     Additionally, transmit filters which are produced using conventional methods to approach the maximum PSDs for the different frequencies of the HDSL2 specification tend to be quite complex having coefficients into the hundreds. This becomes quite difficult to implement in a realistic system. Typically, less accurate lower-order filters are used, giving up some of the allowed power spectrum density at certain frequencies. 
     It is desired to have an improved transmit filter and method for calculating transmit filter coefficients. 
     SUMMARY OF THE PRESENT INVENTION 
     This invention describes a two-step procedure for the design of transmit filter for the communication system with maximum allowable power spectral density. First, the relevant section of the transmit path is partitioned into two parts: (i) the transmit filter that needs to be designed, and ii) the remaining elements that affect the output but not part of the transmit filter, modeled as a fixed weighting function. In contrast to conventional approaches, in one embodiment of the present invention, the transmit filter is selected by optimizing the autocorrelation coefficients via convex optimization methods, and then extracting filter coefficients from autocorrelation coefficients. The structure of the convez optimization allows for a Linear Programming (LP) solution. In the system of the present invention, transmit filter with thirty-two coefficients provides a power spectral density close to the maximum allowed. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     FIG. 1A is a table illustrating the upstream and downstream maximum power spectral density for the HDSL2 system, as specified in the HDSL2 standard. 
     FIG. 1B is a diagram of the maximum power spectral density for the HDSL2 system. 
     FIG. 2 is a diagram of the transmit path for the HDSL2 system. 
     FIG. 3 is a simplified transmit path emphasizing the transmit filter. 
     FIG. 4 is a flowchart illustrating the method of the present invention. 
     FIG. 5 are equations used in the convex optimization procedure. 
     FIG. 6 is a diagram of one embodiment of the magnitude response of an interpolation filter (IF) which uses one embodiment of the present invention. 
     FIG. 7 is a diagram of the response of a sample-and-hold (SH) element used in one embodiment of the present invention. 
     FIG. 8 is a diagram of the frequency response of the analog filters (AF) used in one embodiment of the present invention 
     FIG. 9 is a diagram of the combined (IF, SH and AF) frequency response W(f) in one embodiment of the present invention. 
     FIG. 10 is a table illustrating the transmit filter coefficients for one embodiment of the present invention. 
     FIG. 11 is a diagram that illustrates the downstream PSD. 
     FIG. 12 is a diagram that illustrates the upstream PSD. 
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
     FIGS. 1A and 1B illustrate the PSD limitations of the HDSL2 system as specified in the HDSL2 standard. The present invention can also be used with other DSL specifications or with non-DSL communications specifications. FIG. 1A is a table illustrating the upstream and downstream maximal power spectral density. FIG. 1B is a graph illustrates both the upstream and downstream maximal PSDs. Note that the upstream PSD has a peak power at around 230 KHz and the downstream PSD has a peak power at around 350 KHz. The peaks of the upstream and downstream PSDs correspond to corresponding dips in the downstream and upstream PSDs, respectively, to avoid interference between the upstream and downstream signals. 
     FIG. 2 is a diagram of relevant elements of a DSL transmitter in one embodiment of the present invention. Up-sampler  22  up-samples the input signal by a factor of 2 using zero filling. The transmit filter  24  filters the signal such that the ultimate output is within the power spectral density of the specification. An up-sampler  26  then up-samples the system by a factor of 2 again, using zero filling. The digital low-pass filter  28  or interpolation filter (IF) then filters the output of the up-sampler  26 . Sample-and-hold unit  30  is placed after the digital interpolation filter (IF) and models the effects of the Digital-to-Analog Converter (DAC). Finally, the analog low-pass filter (AF)  32  filters the output of the sample-and-hold (SH) to produce the output of the HDSL2 unit. The effects of all filters except the transmit filter can be modeled as a fixed frequency weighting function, W(f). 
     FIG. 3 is a simplified diagram with fixed elements  34  before the transmit filter and fixed elements  36  after the transmit filter. As will be described below, the system will meet the PSD requirements as long as |W(f)H(f)| 2 ≲S max (f), where S max (f) is the maximum allowable power. 
     FIG. 4 is a flow chart that illustrates the system of the present invention. In step  40 , a communication system specification is obtained. This communication system specification will have a power output limitations for different frequencies. In step  41 , the transmit path is partitioned into the transmit filter and the other elements. In step  42 , in a convex optimization procedure described below, the autocorrelation coefficients of the transit filter are optimized. In step  44 , the filter coefficients are determined from the autocorrelation coefficients. A spectral factorization is preferably used to obtain the filter coefficients from the autocorrelation coefficients. Note that this for the HDSL2 systems the procedures can be done both for the upstream and downstream filters. The convex optimization procedure is described as follows: 
     Let S max (f) denote the specified maximum allowable power in Watts/Hz. Then 
     
       
         | W ( f ) H ( f )| 2   ≲S   max ( f ),  
       
     
     for all frequencies. W(f) can be thought of-as a frequency domain weighting function for the contributions of the other elements in the transmit path besides the transmit filter. 
     For the example described below              S     x   7            (   f   )       =       2       R   L          f   s         |         H   a          (   f   )              H     s                 h            (   f   )              H   1          (   f   )            H        (   f   )              |   2     |       X   1          (   f   )            |   2          ≤       S   max          (   f   )             ,                          
     for all frequencies. And            W        (   f   )            =   Δ              2       R   L          f   s                  H   a          (   f   )              H     s                 h            (   f   )              H   1          (   f   )              X   1          (   f   )           ,                          
     is the combined scaled frequency response of the transit path except the transmit filter. As will be described below, FIG. 9 shows the magnitude response of W(f) where we have assumed that |X(f)| 2 =1; i.e., X 1 (f) is white which is true at the output of a precoder. 
     Now, we describe the optimization procedure for selecting the filter coefficients. 
     Let, r(nT 2 ) be the autocorrelation coefficients associated with the filter impulse response h(nT 2 ), i.e.,            r        (     n                   T   2       )       =       ∑     m   =       -     N   tap       +   1           N   tap     -   1              h        (     m                   T   2       )            h        (       (     m   +   n     )          T   2       )             ,                          
     for  n=−N   tap +1, −N   tap +2 , . . . ,N   tap −1. Let          r        =   Δ            [       r        (   0   )                       r        (     T   2     )                     …                   r        (       (       N   tap     -   1     )          T   2       )         ]     T       ,                          
     be the autocorrelation vector. Then, we can write            |     H        (   f   )            |   2       =       r        (   0   )       +       ∑     n   =   1     N          2        r        (   n   )            cos        (     2      π                 f                 n                   T   2       )               ,                          
     where r(n) is the nth component of the vector r. Note that r(0) is the total power under the magnitude response |H(f)| 2 . Thus, the filter design problem can be recast as the following optimization problem 
     
       
         max r   r (0), such that | W ( f )| 2   |H ( f )| 2   ≲S   max ( f ).  
       
     
     The optimization problem can be written as finite dimensional linear programming (LP) if we use a discrete approximation of the upper bound constraint. We choose a set of N uniformly sampled frequencies            f   k     =       k                   f   s         2      N         ,     k   =   0     ,   1   ,   …              ,     N   -   1     ,                          
     and replace the upper bound constraint for all frequency by N inequality conditions as 
     
       
         | W ( f   k )| 2   |H ( f   k )| 2   ≲S   max ( f   k ),  k= 0,1, . . . , N− 1  
       
     
     For sufficiently large N this discretization yields a good approximation to the original upper bound condition. Using this discretization we rewrite the optimization problem as 
     
       
         max  r   r (0),  
       
     
     
       
         such that | W ( f   k )| 2   |H ( f   k )| 2   ≲S   max ( f   k ),  k= 0, 1, . . . , N − 1.  
       
     
     To write the above problem as an LP, we define an N×N θ&lt; matrix F, a size N diagonal matrix W and two vectors s (N×1) shown in the matrices of FIG.  5 . 
     Using these matrices, we write the optimization problem as an LP 
     
       
         min r    c   T   r,  subject to  Ar≲b    
       
     
     where          A        =   Δ          [           W                 F                 -   W                   F           ]       ,     b        =   Δ          [         s           0         ]       ,                          
     and 0 is a zero vector of compatible dimension (N×1). Note that we added a nonnegativity condition in the optimization problem. This is done so that the vector r have a spectral factor, i.e.,              r        (   0   )       +       ∑     n   =   1         N   tap     -   1            2        r        (   n   )            cos        (     2      π                 f                 n                   T   2       )             ≥   0     ,                          
     for all frequencies. However, due to discretization this is not guaranteed. 
     We solve the LP to obtain the optimal auto-correlation vector r. Once such a vector is obtained, the filter coefficients are then computed using spectral factorization. For problems with small Nap tap , spectral factors can be obtained via root-finding methods. For example, see X. Chen and T. Parks, “Design of optimal minimum phase FIR filters by direct factorization,” Signal Processing, 10:369-383, 1986. 
     FIGS. 6-12 illustrate one embodiment of the system of the present invention. 
     The transmit path for the HDSL2 system with all the relevant components is shown in FIG.  2 . The output of the precoder (not shown) x 1 (nT 1 ) is input to the up-sampler  22 . The sample rate of the input signal R=517 ⅓ ksamples/sec (T 1 =1/R sec). Let              X   1          (   f   )            =   Δ            ∑     n   =     -   ∞       ∞              x   1          (     n                   T   1       )            exp        (       -   j2                   π                 f                 n                   T   1       )             ,                          
     be the frequency response of the input signal x 1 (nT 1 ). This input signal if then up-sampled by a factor of 2 using zero-filling. Thus, the output of the up-sampler is given by            x   2          (     n                   T   2       )            =   Δ          {             x   1          (     m                   T   1       )                 if                 n     =     2      m       ,             0         otherwise   ,                                    
     where T 2 =T 1 /2. In the frequency domain we have                    X   2          (   f   )                       =   Δ            ∑     n   =     -   ∞       ∞              x   2          (     n                   T   2       )            exp        (       -   j2                   π                 f                 n                   T   2       )             ,                              =       ∑     n   =     -   ∞       ∞              x   2          (     n                     T   1     /   2       )            exp        (       -   j2                   π                 f                     T   1     /   2       )             ,                              =         X   1          (   f   )       .                                    
     The output of the up-sampler is then filtered by the transmit filter Txfil. Let the frequency response of the transmit filter be            H        (   f   )            =   Δ            ∑     n   =   0         N   tap     -   1              h        (     n                   T   2       )            exp        (       -   j2                   π                 f                 n                   T   2       )             ,                          
     where h(nT 2 ),n=0,1, . . . ,N tap −1, are the filter coefficients and N tap  is the maximum allowable number of taps. We need to design these coefficients such that the PSD of the output x 7 (t) satisfies the power requirements. 
     Using the filter response H(f), we write the output of the transmit filter in the frequency domain as                      X   3          (   f   )                       =   Δ            ∑     n   =     -   ∞       ∞              x   3          (     n                   T   2       )            exp        (       -   j2                   π                 f                 n                   T   2       )             ,                              =       H        (   f   )              X   2          (   f   )           ,                                =       H        (   f   )              X   1          (   f   )           ,                                       
     The output is then up-sampled again by a factor of 2 using zero-filling and the output of the second up-sampler is given by            x   4          (     n                   T   3       )            =   Δ          {             x   2          (     m                   T   2       )                 if                 n     =     2      m       ,             0           o                 t                 h                 e                 r                 w                 i                 s                 e     ,                                    
     where T 3 =T 2 /2=T,/4, and in the frequency domain we have                      X   4          (   f   )                       =   Δ            ∑     n   =     -   ∞       ∞              x   4          (     n                   T   3       )            exp        (       -   j2                   π                 f                 n                   T   3       )             ,                              =       ∑     n   =     -   ∞       ∞              x   3          (     n                     T   2     /   2       )            exp        (       -   j2                   π                 f                 n                     T   2     /   2       )             ,                                =       H        (   f   )              X   1          (   f   )           ,                              =       X   3          (   f   )                                           
     The output of the second up-sampler x 4 (nT 3 ) is filtered through an interpolating (low pass) filter H 1 (f) (which sits inside the FPGA). For the up-stream we have chosen an 8 th  order interpolating filter with the following coefficients 
     
       
           h   1 ( nT   3 )=[0.052277 0.26772 0.6602 1 1 0.6602 0.26772 0.052277],  
       
     
     and for the down-stream we have chosen a 5 th  order interpolating filter with the following coefficients 
     
       
           h   1 ( nT   3 )=[0.19811 0.6949 1 0.6949 0.19811],  
       
     
     The frequency response of these two filters are shown in FIG.  6 . The output of the interpolating filter in the frequency domain is given by                      X   5          (   f   )                       =   Δ            ∑     n   =     -   ∞       ∞              x   5          (     n                   T   3       )            exp        (       -   j2                   π                 f                 n                   T   3       )             ,                              =         H   1          (   f   )              X   4          (   f   )           ,                                =         H   1          (   f   )            H        (   f   )              X   1          (   f   )           ,                                       
     The output of the interpolating filter is then passed through a sample and hold (S/H) circuit followed by an analog filter H a (f). The output of the S/H is given by              x   6          (   t   )            =   Δ            x   5          (     nT   3     )         ,                  nT   3     ≤   t   &lt;       (     n   +   1     )            T   3     .                                
     Thus,                    X   6          (   f   )            =   Δ                       ∫     -   ∞     ∞              x   6          (   t   )            exp        (       -   j2                   π                 f                 t     )               t           ,                 =                  ∑     n   =     -   ∞       ∞            ∫     n                   T   3           (     n   +   1     )          T   3                  x   6          (   t   )            exp        (       -   j2                   π                 f                 t     )               t             ,                 =                  ∑     n   =     -   ∞       ∞            ∫     n                   T   3           (     n   +   1     )          T   3                  x   5          (     n                   T   3       )            exp        (       -   j2                   π                 f                 t     )               t             ,                                
     Where we have defined            H     s                 h            =   Δ              T   3          exp        (       -   j                   π                 f                   T   3       )              sin        (     π                 f                   T   3       )         π                 f                   T   3           =       T   3          exp        (       -   j                   π                 f                   T   3       )          sin                   c        (     f                   T   3       )             ,                          
     The magnitude plot the function H sh  is shown in FIG.  7 . Now, we write 
     
       
           X   6 ( f )= H   sh   H   1 ( f ) H ( f ) X   1 ( f ).  
       
       
         
           
             
               = 
               
                 
                   ∑ 
                   
                     n 
                     = 
                     
                       - 
                       ∞ 
                     
                   
                   ∞ 
                 
                  
                 
                   
                     
                       x 
                       5 
                     
                      
                     
                       ( 
                       
                         n 
                          
                         
                             
                         
                          
                         
                           T 
                           3 
                         
                       
                       ) 
                     
                   
                    
                   
                     
                       ∫ 
                       
                         n 
                          
                         
                             
                         
                          
                         
                           T 
                           3 
                         
                       
                       
                         
                           ( 
                           
                             n 
                             + 
                             1 
                           
                           ) 
                         
                          
                         
                           T 
                           3 
                         
                       
                     
                      
                     
                       
                         exp 
                          
                         
                           ( 
                           
                             
                               - 
                               j2 
                             
                              
                             
                                 
                             
                              
                             π 
                              
                             
                                 
                             
                              
                             f 
                              
                             
                                 
                             
                              
                             t 
                           
                           ) 
                         
                       
                        
                       
                          
                         t 
                       
                     
                   
                 
               
             
             , 
             
               
 
             
              
             
               = 
               
                 
                   
                     [ 
                     
                       
                         ∑ 
                         
                           N 
                           = 
                           
                             - 
                             ∞ 
                           
                         
                         ∞ 
                       
                        
                       
                         
                           
                             x 
                             5 
                           
                            
                           
                             ( 
                             
                               n 
                                
                               
                                   
                               
                                
                               
                                 T 
                                 3 
                               
                             
                             ) 
                           
                         
                          
                         
                           exp 
                            
                           
                             ( 
                             
                               
                                 - 
                                 j2 
                               
                                
                               
                                   
                               
                                
                               π 
                                
                               
                                   
                               
                                
                               f 
                                
                               
                                   
                               
                                
                               n 
                                
                               
                                   
                               
                                
                               
                                 T 
                                 3 
                               
                             
                             ) 
                           
                         
                       
                     
                     ] 
                   
                    
                   
                     
                       1 
                       - 
                       
                         exp 
                          
                         
                           ( 
                           
                             
                               - 
                               j2 
                             
                              
                             
                                 
                             
                              
                             π 
                              
                             
                                 
                             
                              
                             f 
                              
                             
                                 
                             
                              
                             
                               T 
                               3 
                             
                           
                           ) 
                         
                       
                     
                     
                       j2 
                        
                       
                           
                       
                        
                       π 
                        
                       
                           
                       
                        
                       f 
                     
                   
                 
                  
                 
                   
 
                 
                 = 
                 
                   
                     
                       X 
                       5 
                     
                      
                     
                       ( 
                       f 
                       ) 
                     
                   
                    
                   
                     T 
                     3 
                   
                    
                   
                     exp 
                      
                     
                       ( 
                       
                         
                           - 
                           j 
                         
                          
                         
                             
                         
                          
                         π 
                          
                         
                             
                         
                          
                         f 
                          
                         
                             
                         
                          
                         n 
                          
                         
                             
                         
                          
                         
                           T 
                           3 
                         
                       
                       ) 
                     
                   
                    
                   
                     
                       sin 
                        
                       
                         ( 
                         
                           π 
                            
                           
                               
                           
                            
                           f 
                            
                           
                               
                           
                            
                           
                             T 
                             3 
                           
                         
                         ) 
                       
                     
                     
                       π 
                        
                       
                           
                       
                        
                       f 
                        
                       
                           
                       
                        
                       
                         T 
                         3 
                       
                     
                   
                 
               
             
             , 
             
               
 
             
              
             
               = 
               
                 
                   
                     X 
                     5 
                   
                    
                   
                     ( 
                     f 
                     ) 
                   
                 
                  
                 
                   
                     H 
                     
                       s 
                        
                       
                           
                       
                        
                       h 
                     
                   
                    
                   
                     ( 
                     f 
                     ) 
                   
                 
               
             
             , 
           
         
                 
         
             
         
      
     
     For the analog filter H a (f). We have chosen a 4 th  order Butterworth filter with cut-off frequency f c =300 kHz for the up-stream and an 8 th  order Butterworth filter with cut-off frequency f c =440 kHz kHz for the down-stream. The squared magnitude response of an Nth order Butterworth filter with cut-off frequency f c  is given by          |       H   a          (   f   )            |   2       =       1     1   +       (     f   /     f   c       )       2      N           .                            
     The magnitude response of the up-stream and down-stream analog filters are shown in FIG.  8 . 
     Finally, the frequency response of x 7 (t), the output of the analog filter H a (f) is given by                    X   7          (   f   )            =   Δ                       ∫     -   ∞     ∞              x   7          (   t   )            exp        (       -   j2                   π                 f                 t     )                          t           ,                 =                    H   a          (   f   )              X   6          (   f   )           ,                 =                    H   a          (   f   )              H     s                 h            (   f   )              H   1          (   f   )            H        (   f   )              X   1          (   f   )           ,                                
     The HDSL2 standard specifies the PSD of x 7 (t) in terms of dBm/Hz and the power is measured with a load impedance of R L =135Ω. Now the one-sided PSD S x     7   (f) and the magnitude response X 7 (f) are related as follows              S     x   7            (   f   )       =       2       R   L          f   s                       X   7          (   f   )            2         ,                          
     where f s  is the sampling frequency in Hz. Thus,            S     x   7            (   f   )       =       2       R   L          f   s                         H   a          (   f   )              H   sh          (   f   )              H   1          (   f   )            H        (   f   )              2                     X   1          (   f   )            2     .                              
     Using the techniques described above, a thirty-two coefficient digital transmit filter for the upstream and downstream directions are shown in the Table of FIG.  10 . 
     FIG. 11 illustrates the overall downstream design PSD and specification PSD, assuming |X(f) 2 =1. 
     FIG. 12 illustrates both the designed and specified overall upstream power spectral densities. 
     Note that the transmit filter coefficients depend upon the other elements in the transmit path since the other elements affect the W(f) frequency domain weighting function. Thus, for example, if a different interpolating filter were used, the coefficients of the transmit filter would have to be modified. 
     It will be appreciated by those of ordinary skill in the art that the invention can be implemented in other specific forms without departing from the spirit or character thereof. The presently disclosed embodiments are therefore considered in all respects to be illustrative and not restrictive. The scope of the invention is illustrated by the appended claims rather than the foregoing description, and all changes that come within the meaning and range of equivalents thereof are intended to be embraced herein.