# Patent

## Title
Peak-to-average power ratio (PAR) reduction based on active-set tone reservation

## Abstract

In an embodiment, a transmitter includes first and second processing blocks, which may each include hardware, software, or a combination of hardware and software. The first processing block is operable to generate a first peak-reducing vector. And the a second first processing block is operable to receive a first data vector, the data vector comprising a plurality of samples, the first data vector having a first peak with a first index and a first magnitude, a second peak with a second index and a second magnitude that is less than the first magnitude, and a first peak-to-average power ratio, and to generate a second data vector having a second peak-to-average power ratio that is lower than the first peak-to-average power ratio by using the first peak-reducing vector.

## Background
PRIORITY CLAIM
The present application is a Continuation of copending U.S. patent application Ser. No. 12/860,142 filed Aug. 20, 2010, which application claims priority to Chinese Patent Application No. 200910265275.7, filed Dec. 28, 2009, all of the foregoing applications are incorporated herein by reference in their entireties.


BACKGROUND
Multicarrier transmission has been widely adopted in both wired and wireless communication systems such as asymmetric digital subscriber line (ADSL) systems, Digital Video Broadcast (DVB), and wireless local/metropolitan area networks (WLAN/WMAN). Exploiting Discrete Multitone Modulation (DMT) or Orthogonal Frequency Division Multiplexing (OFDM), these systems may achieve greater immunity to multipath fading and impulse noise with lower cost. However, they may also suffer from high peak-to-average power ratios (PAR). Without an appropriate process to counter this problem, the high PAR of a transmitted signal may cause a high-power amplifier (HPA) to operate in its nonlinear region (i.e., the peak-to-peak amplitude of the transmitted signal may be high enough to saturate the amplifier), leading to significant performance degradation.
OFDM effectively partitions overall system bandwidth into a number of orthogonal frequency subchannels. These subchannels are also interchangeably referred to throughout as “tones” or “subcarriers.” In an OFDM system, an input serial data symbol is separated into D groups. Each of the D groups may be mapped onto a quadrature amplitude modulated (QAM) constellation point, and then modulated onto a respective one of N subchannels (or tones) having approximately equal bandwidth and a frequency separation of approximately 1/T, where T is the time duration of an OFDM symbol during which all N groups are transmitted, and D≦N. Generally, the larger the value of D, the larger the system bandwidth, and, because of a resulting quasi-Gaussian distribution in the resulting time-domain signal, the higher the PAR (for example, peak amplitude to average amplitude).
Tone reservation, which modulates reserved or unused ones of the N tones within the signal space to produce data-block-dependent peak-canceling signals, is a technique for reducing the PAR for these systems. That is, where D<N, one may modulate one or more of the unused ones of the N tones to reduce the PAR of the transmitted signal.
SUMMARY
A challenge in tone reservation is how best to produce that peak-canceling signal. Unfortunately, known solutions have typically involved high computational overhead, inaccuracy, or both.
According to an embodiment, a transmitter includes first and second processing blocks, which may each include hardware, software, or a combination of hardware and software. The first processing block is operable to generate a first peak-reducing vector. And the second first processing block is operable to receive a first data vector, the data vector comprising a plurality of samples, the first data vector having a first peak with a first index and a first magnitude, a second peak with a second index and a second magnitude that is less than the first magnitude, and a first peak-to-average power ratio, and to generate a second data vector having a second peak-to-average power ratio that is lower than the first peak-to-average power ratio by using the first peak-reducing vector.



BRIEF DESCRIPTION OF THE DRAWINGS
Embodiments of the subject matter disclosed herein will become more readily appreciated as the same become better understood by reference to the following detailed description, when taken in conjunction with the accompanying drawings.
FIG. 1 illustrates an embodiment of an iterative method for PAR reduction using a peak-canceling vector.
FIG. 2 is a flowchart of the embodiment of the iterative method for PAR reduction depicted in FIG. 1.
FIG. 3 is a plot of an example time-domain waveform that may be generated for the purpose of PAR reduction using reserved OFDM tones.
FIG. 4 is a plot of an example time-domain waveform that may be appropriate for PAR reduction using reserved OFDM tones.
FIG. 5 is a plot of an example time-domain waveform that may result from using a single iteration of an embodiment of a PAR reduction method.
FIG. 6 is a plot of an example time-domain waveform that may result from using a second iteration of an embodiment of a PAR reduction method.



DETAILED DESCRIPTION
The following discussion is presented to enable a person skilled in the art to make and use the subject matter disclosed herein. The general principles described herein may be applied to embodiments and applications other than those detailed above without departing from the spirit and scope of the subject matter disclosed herein. This disclosure is not intended to be limited to the embodiments shown, but is to be accorded the widest scope consistent with the principles and features disclosed or suggested herein.
In the following description, certain details are set forth in conjunction with the described embodiments to provide a sufficient understanding of the subject matter disclosed herein. One skilled in the art will appreciate, however, that the disclosed subject matter may be practiced without these particular details. Furthermore, one skilled in the art will appreciate that various modifications, equivalents, and combinations of the disclosed embodiments and components of such embodiments are within the scope of the disclosed subject matter. Embodiments including fewer than all the components or steps of any of the respective described embodiments may also be within the scope of the disclosed subject matter although not expressly described in detail below. Finally, the operation of well-known components and/or processes has not been shown or described in detail below to avoid unnecessarily obscuring the disclosed subject matter.
FIG. 1 illustrates an embodiment of an iterative method for PAR reduction using a peak-canceling vector. A time-domain data vector x(t), labeled x0 for notational convenience, is the result of performing an Inverse Discrete Fourier Transform (IDFT) with an oversampling factor L on a frequency-domain OFDM symbol X that includes N information blocks X0 . . . XN-1 respectively modulated on N subcarriers, plus D PAR reduction blocks p0 . . . pD-1 respectively modulated on D subcarriers, where D≦(L−1)N. A time-domain peak-canceling vector p0, discussed in further detail below, is generated by performing an IDFT on a series of frequency impulses (ones) at one or more of the D reserved tones. An example of a peak-canceling vector p0 is shown in FIG. 3.
If the PAR of the data vector x0 is above a target threshold, then a first peak reduction is performed on the data vector x0 by locating the sample n1 of x0 that has the greatest magnitude within the data vector; circularly shifting the peak-canceling vector p0 by that number of samples n1 to get a shifted peak-canceling vector pn1; scaling pn1 by a step size −μ11 and complex coefficient α11; and combining (e.g., summing) the original data vector x0 with the shifted and scaled peak-canceling vector pn1 to get a PAR-reduced data vector x1. In this manner, in the vector x1, the magnitude of the sample at n1 is reduced to the magnitude of the sample at n2 such that the magnitudes of n1 and n2 are substantially equal and are the greatest peak magnitudes of x1. The step size μ is a real quantity, used to modify the magnitude of the vector pn1. The complex coefficient α is a complex coefficient vector used to modify both the magnitude and phase of the data vector pn1, because the peak at n1 of x0 has a phase as well as a magnitude, and so both the phase and magnitude of pn1 are scaled so as to give the desired peak reduction at n1. Both μ and α are discussed in greater detail below.
In an embodiment, if the PAR of the data vector x1 is above the target threshold, then a second peak reduction is performed on x1 in a similar fashion. For this second iteration of peak reduction, the original peak-canceling vector p0 is shifted by that number of samples n2 to achieve a shifted peak-canceling vector pn2. Both μ and α are regenerated. The shifted vector pn2 is scaled by a step size −μ22 and a complex coefficient α22; the shifted vector pn1 is scaled by a step size −μ12 and a complex coefficient α12. The complex coefficients α cause the samples n1 and n2 of the combined shifted and scaled peak-canceling vectors pn1 and pn2 to each have an appropriate magnitude and the same phase (positive or negative) as the samples n1 and n2 of x1. The data vector x1 is combined with these shifted and scaled peak-canceling vectors pn1 and pn2 to achieve the PAR-reduced data vector x2. The real coefficients μ are used to reduce the magnitudes of x1's samples at n1 and n2 to substantially equal that of the sample at n3 of x2. These three peaks (located at samples n1, n2 and n3) will have approximately the same magnitude value, which is the greatest magnitude value in x2.
In an embodiment, individual peak cancellation is reiterated by repeating the above technique until the PAR of the resulting vector x equals or is less than the target threshold. In another embodiment, a maximum number of iterations number may be set even if the PAR of the resulting vector x is not less than the target threshold.
The time-domain waveform x(t) is the result of performing an IDFT on a frequency-domain OFDM symbol X, where the symbol X comprises N information blocks modulated on respective N subcarriers. For complex baseband OFDM systems, this time-domain waveform may be represented as:




  
    
      
        
          
            x
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              1
              
                N
              
            
            ⁢
            
              
                ∑
                
                  n
                  =
                  0
                
                
                  N
                  -
                  1
                
              
              ⁢
              
                
                  X
                  n
                
                ·
                
                  e
                  
                    j2π
                    ⁢
                    
                      
                    
                    ⁢
                    n
                    ⁢
                    
                      
                    
                    ⁢
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    f
                    ⁢
                    
                      
                    
                    ⁢
                    t
                  
                
              
            
          
        
        ,
        
          0
          ≤
          t
          ≤
          NT
        
      
    
    
      
        (
        1
        )
      
    
  




where Xn, {n=0, 1, . . . , N−1} is the data symbol modulated onto the nth subchannel, Δf=1/NT is the subcarrier spacing, T is the time duration of the OFDM symbol.

The PAR of this transmitted waveform x(t) is defined as:




  
    
      
        PAR
        =
        
          
            
              max
              
                0
                ≤
                t
                <
                NT
              
            
            ⁢
            
              
                
                
                  x
                  ⁡
                  
                    (
                    t
                    )
                  
                
                
              
              2
            
          
          
            
              1
              /
              NT
            
            ·
            
              
                ∫
                0
                NT
              
              ⁢
              
                
                  
                    
                    
                      x
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                    
                  
                  2
                
                ⁢
                
                  d
                  t
                
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




To use tone reservation, the LN subcarriers (wherein L is the oversampling factor) in the OFDM system are divided into two subsets: a subcarrier set U for useful data (the tones N are in the set U) and a subcarrier set Uc for D (L−1)N PAR reduction blocks selected to reduce the PAR. As an example, in an embodiment wherein the signal space comprises a total of LN=256 subcarriers, one might reserve for PAR reduction a set Uc of D=11 of those subcarriers with indexes k=[5, 25, 54, 102, 125, 131, 147, 200, 204, 209, 247]. Expressing the sampled x(t) waveform as a sampled data vector x, one may also express a desired peak-cancelling waveform as a sampled vector c. The time-domain output of the transmitter x is the sum of the data vector x and peak-canceling vector c, which is given by:

 x=x+c=IDFT(X+C)  (3)

Note that frequency-domain vector X=[X0, X1, . . . , XN-1] is only modulated over subcarriers N within the data-carrying subset U; likewise, frequency-domain vector C=[CO, C1, . . . , CN-1] is only modulated over subcarriers D within the reserved-tone subset Uc. The vectors X and C cannot both be nonzero on a given subcarrier; that is, a subcarrier may be a member of U or of Uc, but cannot be a member of both U and Uc:





  
    
      
        
          
            X
            k
          
          +
          
            C
            k
          
        
        =
        
          {
          
            
              
                
                  
                    X
                    k
                  
                  ,
                  
                    k
                    ∈
                    U
                  
                
              
            
            
              
                
                  
                    C
                    k
                  
                  ,
                  
                    k
                    ∈
                    
                      U
                      c
                    
                  
                
              
            
          
        
      
    
    
      
        (
        4
        )
      
    
  




C, therefore, is selected to reduce or minimize the PAR of the time-domain output signal x by way of its time-domain correlation c:





  
    
      
        PAR
        =
        
          
            
              max
              
                0
                ≤
                k
                ≤
                
                  NL
                  -
                  1
                
              
            
            ⁢
            
              
                
                
                  
                    x
                    ⁡
                    
                      [
                      k
                      ]
                    
                  
                  +
                  
                    c
                    ⁡
                    
                      [
                      k
                      ]
                    
                  
                
                
              
              2
            
          
          
            E
            ⁢
            
              {
              
                
                  
                  
                    x
                    ⁡
                    
                      [
                      k
                      ]
                    
                  
                  
                
                2
              
              }
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




The calculation of a suitable vector for the time-domain vector c, without high computational overhead, may therefore result in PAR reduction for an OFDM system.

FIG. 2 illustrates an embodiment of a method for calculating a peak-cancelling vector to achieve PAR reduction. As described above with respect to FIG. 1, PAR reduction is performed iteratively, with a separate stage for cancellation of each successive peak in the time-domain data vector x, until a PAR threshold (or iteration threshold) is reached. A first peak of greatest magnitude within data vector x is found and reduced using a shifted peak-canceling vector scaled by real and complex coefficients until it is of substantially equal magnitude with another peak, e.g., the peak with the next greatest magnitude. With a combined peak-canceling vector which has unit magnitude and the phase of the reduced data vector at the locations of the two peaks, those peaks are reduced by an equal amount such that they then have a magnitude substantially equal to that of a third peak having the third greatest magnitude, and so on. In an embodiment, individual peak cancellation is reiterated until a target value for PAR is achieved. In another embodiment, a maximum iteration number M may be designated.
Starting in step 205, an N-subcarrier frequency symbol X={X0, X1, . . . , XN-1} is converted into a first time-domain sampled data vector x0 using an Inverse Discrete Fourier Transform (IDFT) with L-times oversampling. A second L-oversampled vector p0, termed the peak-reduction kernel, may be computed in such a way as to provide a time-domain impulse at sample number n=0 with the D reserved tones within set Uc. That is, to form the vector p0, the D reserved tones (the D subcarriers within set Uc) are set to “1+j0, all other tones in the subcarrier space (the N subcarriers within set U) are set to “0+j0”, and an IDFT is performed to achieve the time-domain peak-reduction vector p0, which is then normalized/scaled so that the vector p0 has a maximum magnitude of 1 at the sample location n=0. The phases of the D reserved tones are adjusted such that ideally, p0 has a peak of highest magnitude at n=0, and the magnitudes of the peaks at all the other sample locations are significantly lower (e.g., zero) than the peak at n=0. But because the set Uc has a finite set of D tones, at least some of the peaks at the other sample locations n≠0 may have non-zero magnitudes.
In step 210, the first peak in data vector x0 is determined such that the sample at that peak has the greatest magnitude in the vector. Let an iteration index i=1, and let E0 be the maximum magnitude in x0. Further, let sample n1 be the location where E0 is found. This sample n1 is the location of the first peak within vector x0 to be partially cancelled (i.e., reduced) to achieve PAR reduction. Also at step 210, this first peak is used to establish an active set of peaks A. At this point, the active set A contains only the first peak at sample location n1.
In step 215, a first complex coefficient α11 is determined. Coefficient α11 is the complex ratio between the complex value of x0 at sample location n1 and the maximum magnitude E0, such that




  
    
      
        
          p
          1
        
        =
        
          
            
              α
              1
              1
            
            ⁢
            
              p
              
                n
                1
              
            
          
          =
          
            
              
                x
                
                  n
                  1
                
                0
              
              
                E
                0
              
            
            ⁢
            
              p
              
                n
                1
              
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




where pn1 is obtained by circularly shifting peak-reduction kernel vector p0 by n1 samples, and xn10 is the complex value of the sample n1 only α11 and p1 are determined with equation (8); successive iterations, wherein the more generalized pi is calculated, are discussed below.

At step 220, peak testing is performed to locate samples in xi−1 (which, for the first peak, is x0) which are not yet in the active set A that possess large magnitudes relative to the other samples. These samples are candidates for peak balancing in the next step, and are placed in a test set B. For the very first peak reduction, active set A contains only the one peak at n1. Because it is possible for multiple peaks to have the same magnitude, however, it is possible that multiple peaks may be added to active set A in a single iteration. In one embodiment, the magnitude of these samples may be approximated by simply taking the sum of the absolute values of a given sample's real and imaginary components. In still another embodiment, peak testing may be skipped entirely. However, if peak testing is skipped, then the test set B will include all samples not in active set A, which greatly increases the computational load necessary for the next step 225.
In step 225, the minimum step size μi is determined as the difference in magnitude between the maximum magnitude of current data vector xi−1 and that of the next iteration xi, such that Ei=Ei−1−μi, and μi=Ei−1−Ei. This minimum step size μi may also be found according to the following equation:




  
    
      
        
          μ
          i
        
        =
        
          
            min
            
              q
              ∈
              B
            
          
          ⁢
          
            (
            
              
                
                  
                    b
                    q
                  
                  -
                  
                    
                      
                        b
                        q
                        2
                      
                      -
                      
                        
                          a
                          q
                        
                        ⁢
                        
                          c
                          q
                        
                      
                    
                  
                
                
                  a
                  q
                
              
              ≥
              0
            
            )
          
        
      
    
    
      
        (
        9
        )
      
    
  




where aq=1−pqi(pqi)*, bq=Ei−1−xqi−1(pqi)*), cq=(Ei−1)2−xqi−1(xqi−1)*. The complex conjugate of any variable a is denoted by a*, and (•) denotes the real part of a complex number. In embodiments wherein step 220 (peak testing) was skipped, the determination of minimum step size μi may be more computationally intensive because test set B will be much larger. In either case, the peaks within xi−1 having the complex magnitude associated with μi are then added to active set A.

At this point, if the iteration index i=1 (i.e., we are in the process of reducing the magnitude of the first peak), then we have already calculated the complex coefficient α11 in step 215 and we can proceed to step 235 to calculate the next iteration of data vector xi.
In step 230, the complex coefficients α are determined. As stated above, these complex coefficients are used both to scale the magnitude of the relevant peaks and to adjust the phase of those peaks such that at the locations within active set A the combined peak-reduction vector pi has unit magnitude and the same phase as the corresponding samples in xi−1. To solve for the α coefficients, the following complex matrix equation is used:




  
    
      
        
          
            [
            
              
                
                  1
                
                
                  
                    p
                    
                      
                        n
                        1
                      
                      -
                      
                        n
                        2
                      
                    
                  
                
                
                  ...
                
                
                  
                    p
                    
                      
                        n
                        1
                      
                      -
                      
                        n
                        i
                      
                    
                  
                
              
              
                
                  
                    p
                    
                      
                        n
                        2
                      
                      -
                      
                        n
                        1
                      
                    
                  
                
                
                  1
                
                
                  ...
                
                
                  
                    p
                    
                      
                        n
                        2
                      
                      -
                      
                        n
                        i
                      
                    
                  
                
              
              
                
                  ⋮
                
                
                  
                    
                  
                
                
                  
                    
                  
                
                
                  ⋮
                
              
              
                
                  
                    p
                    
                      
                        n
                        i
                      
                      -
                      
                        n
                        1
                      
                    
                  
                
                
                  
                    p
                    
                      
                        n
                        i
                      
                      -
                      
                        n
                        1
                      
                    
                  
                
                
                  ...
                
                
                  1
                
              
            
            ]
          
          ⁡
          
            [
            
              
                
                  
                    α
                    1
                    i
                  
                
              
              
                
                  
                    α
                    2
                    i
                  
                
              
              
                
                  ⋮
                
              
              
                
                  
                    α
                    i
                    i
                  
                
              
            
            ]
          
        
        =
        
          [
          
            
              
                
                  S
                  
                    n
                    1
                  
                
              
            
            
              
                
                  S
                  
                    n
                    2
                  
                
              
            
            
              
                ⋮
              
            
            
              
                
                  S
                  
                    n
                    i
                  
                
              
            
          
          ]
        
      
    
    
      
        (
        10
        )
      
    
  




Where pn is the nth entry of p0, and

Snl=xnli−1/Ei−1  (11)

All of the coefficients α are recalculated during each iteration, with the leftmost matrix of equation (10) gaining another row and column each time. In an embodiment, this i×i system of complex matrix equations may be replaced with a 2i×2i system of real equations. This and other simplification techniques are discussed in “An Active-Set Approach for OFDM PAR Reduction via Tone Reservation”, Brian S. Krongold and Douglas L. Jones, IEEE Transactions on Signal Processing, Vol. 52, No. 2, pp. 495-509, February 2004, the contents of which article are incorporated herein by reference.

Once the coefficients [α1i, α2i, . . . , αii] are calculated in step 230, then:

pi=Σl=1iαlipnl  (12)

The iterative data vector xi is determined by adding pi, negatively scaled by the minimum step size μi, to the prior iterative vector xi−1, so that:

xi=xi−1−μipi  (13)

This has the result of reducing the magnitudes of a number, e.g., i, peaks in the active set by the magnitude of μi.

The iteration index i is incremented such that i=i+1. If neither the maximum number of iterations nor the desired PAR value has been achieved, then the next iteration of peak reduction is begun by returning to step 220 for peak testing (or, in embodiments where peak testing is skipped, step 225 for directly finding the minimum step size).
If the desired PAR or maximum iterations M have been achieved, then the output x is equal to current data vector xi.
Returning to FIG. 1, in the leftmost stage 101, a first Inverse Discrete Fourier Transform (IDFT) block 110 provides an LN-point oversampled time-domain vector x for each frequency block X={X0, X1, . . . , XN-1}. Also in stage 201, a second IDFT block 120 calculates an LN-point oversampled vector p0. This vector p0 may be computed in such a way as to provide a time-domain impulse at sample number n=0 with the D reserved tones within set Uc, and scaling p0 such that it possesses a value of 1 at the location n=0. Frequency-domain impulses occur for every reserved tone D in the set, so that P={P0, P1, . . . , PN-1}. In one embodiment, these reserved tones are unchanging, so that the same frequency domain vector P may be used for all blocks X. In other embodiments, the reserved tones may vary depending on the value of X, such that P (and, therefore, p) are recalculated for each block X.
Stage 102 illustrates the first iteration of peak reduction. A first peak within data vector x0, having the maximum magnitude E0, is located at sample n1 and added to an active peak set A (not shown). Reduction kernel vector p0 is circularly shifted by that sample number n1 so that the highest amplitude in p0, formerly located at n=0, is now aligned at n=n1. Test set B is achieved by peak testing to find those samples which are not yet in the active set A and possess large magnitudes relative to the other samples in x0. Using equation (9), for example, this allows the determination of the minimum step size μ1, which is the difference in magnitude between the first peak, with magnitude E0, and the peak with the second-greatest complex magnitude in x0.
The first complex coefficient α11 is then determined, p1 is found as α11 pn1, and c1 is found as pn1(−μ1α11). By scaling pn1 by the complex coefficient α11, p1 has the same phase as that of x0 at n1. By scaling p1 by the negative step size −μ1, the magnitude of the peak at n1 is reduced to substantially match that of the next peak or peaks to be reduced. This is reflected in the resulting vector x1, which possesses both the peak located at sample n1 as well as the secondary peak located at n=n2 when determining the minimum step size μ1—but in the data vector x1, each of those peaks possess substantially the same magnitude.
In stage 103, the process is repeated with iteration index i=2. Here, two peaks are being reduced: the first peak located at n1, and another secondary peak located at n2. The minimum step size μ2 may be found using equation (9) based on the peaks of data vector x1 within the test set B; complex coefficients α12 and α22 are determined using equation (10); and the peak-reduction vector c2 is found as −μ2p2=pn1(−μ2α12)+pn2(−μ2α22). The subsequent data vector is thus found as x2=x1+c2. The coefficients α12 and α22 are calculated so that even if pn1 has a nonzero magnitude at sample n2, the magnitude of the reduced peak at n2 still substantially equals the magnitude of the peak at n3, and even if pn2 has a nonzero magnitude at sample n1, the magnitude of the reduced peak at n1 still substantially equals the magnitude of the peak at n3.
Stage 104 illustrates the Mth iteration of an embodiment of the peak reduction method. A total of M peaks are being reduced in complex magnitude simultaneously, with the first M−1 peaks already having been reduced by previous iterations. The output data vector is x=x+c, where x is the original data vector given by IDFT(X) and c is the time-domain sum of c1+c2+ . . . +cM.
Referring again to FIG. 1, an embodiment of the above-described PAR reduction technique is described with less emphasis on the mathematics and more emphasis on the physical phenomena.
Referring to block 101, an IDFT of an OFDM data signal is generated, and an IDFT of an OFDM peak-reduction kernel is generated. The OFDM data signal includes subcarriers that are selected for data transmission, and the OFDM peak-reduction kernel includes subcarriers that are not used for data transmission (e.g., because of excessive channel interference, distortion at the frequencies of these non-data subcarriers, or a desire to reserve those subcarriers for purposes of PAR reduction). The first IDFT yields a time-domain data signal x0 having a duration equal to a symbol period, and the second IDFT yields a time-domain peak-reduction pulse p0 also having a duration of one symbol period. An example of p0 is plotted in FIG. 3. The first and second IDFTs may be performed by a processor or other circuitry in a transmitter.
Referring to block 102, the transmitter identifies the peaks of x0 having the greatest magnitude and the second greatest magnitude as being located at sample locations n1 and n2, respectively.
Then, the transmitter circularly shifts p0 to generate a peak-reduction pulse p0 having its main peak at sample location n1.
Next, the transmitter generates −μ1 having a value such that when added to the magnitude of x0 at sample location n1, −μ1α11pn1 causes the magnitude of the peak of x1 at sample n1 to equal, or approximately equal, the magnitude of the peak of x0 at sample n2, as shown in block 103. The transmitter then adds −μ1α11pn1 to x0 to obtain the time-domain waveform x1 in block 103.
Then, referring to block 103, the transmitter identifies the peaks of x1 at sample locations n1 and n2 as having the greatest magnitude, and identifies the peak of x1 at n3 as having the second greatest magnitude.
Next, the transmitter circularly shifts p0 to generate a peak-reduction pulse pn1 having its main peak at sample location n1, and also circularly shifts p0 to generate a peak-reduction pulse pn2 having its main peak at sample location n2.
Then, the transmitter generates −μ2 having a value such that when added to the magnitude of x1 at sample location n1, −μ2α12pn1 causes the magnitude of the peak of x2 at sample n1 to equal, or approximately equal, the magnitude of the peak of x1 at sample n3. Furthermore, the transmitter generates α12 and a22 such that the magnitude of x2 at sample location n2 is approximately equal to the magnitude of x1 at the sample location n3 even if the magnitude of p0 at sample location n2 is nonzero.
Because the magnitudes of x1 at sample locations n1 and n2 are substantially equal to one another, when added to the magnitude of x1 at sample location n2, −μ2α22pn2 causes the magnitude of the peak of x2 at sample n2 to equal, or approximately equal, the magnitude of the peak of at sample n3. Furthermore, the transmitter generates α12 and α22 such that the magnitude of x2 at sample location n1 is approximately equal to the magnitude of x1 at the sample location n3 even if the magnitude of pn2 at sample location n1 is nonzero.
The transmitter then adds −μ2α12pn1 and −μ2α22pn2 to x1 to obtain the time-domain waveform x2, which is output from the block 104.
Referring to block 104, the transmitter continues in this manner for M−3 additional peaks until it generates a time-domain data signal xM having no peaks larger than a selected peak threshold, or until M equals a selected iteration threshold.
FIG. 3 is a plot of an example of a time-domain peak-reduction waveform p0 that may be generated for the purpose of PAR reduction using reserved OFDM tones. In the embodiment depicted, the kernel p0 is generated using D=11 reserved tones in a signal space having a total of LN=256 subcarriers. As discussed above, the kernel may have, therefore, been generated using a 256-point IDFT. One can see that the kernel p0 has its greatest peak of approximately unity (normalized) magnitude at sample n0.
FIG. 4 is a time-domain magnitude plot of a complex-baseband OFDM signal using the same system parameters as in the embodiment depicted in FIG. 3. Again, the signal space comprises 256 subcarriers, 11 of which are reserved for the purpose of PAR reduction. In the waveform depicted by FIG. 4, the highest peak may be seen as occurring in the area of sample index n150. The waveform has a PAR of 10.45 dB.
FIG. 5 is a magnitude plot of the complex-baseband OFDM signal from FIG. 4 after a single iteration of the PAR reduction method described above using the peak-reduction kernel p0 depicted in FIG. 3. As a result of that single iteration, the PAR of the signal has been reduced to approximately 7.56 dB. In the depiction of FIG. 5, the peak located approximately at sample index n150 has been reduced in magnitude to approximate that of a second peak, located in the region around sample index n7.
FIG. 6 is a magnitude plot of the complex-baseband OFDM signal from FIG. 3 after two successive iterations of the PAR reduction method described above using the peak-reduction kernel p0 depicted in FIG. 3. As a result of these two successive iterations, the PAR of the signal has been reduced to 6.92 dB. In the depiction of FIG. 6, the peaks located approximately at sample indexes n150 and n7, respectively, have been reduced in magnitude to approximate that of a third peak, located in the area around sample index n238.
It is to be understood that even though various embodiments and advantages have been set forth in the foregoing description, the above disclosure is illustrative only, and changes may be made in detail, and yet remain within the broad principles of the disclosed subject matter. For example, the method and system may be implemented in either software or hardware embodiments, and may comprise one or more integrated circuit devices. In some embodiments, the methods or individual steps described may be performed in a hardware implementation. In other embodiments, a software implementation may be utilized. In still other embodiments, the methods or individual steps described may be performed by a combination of hardware and software modules. Furthermore, p0 may have its main pulse at any sample location other than n0. Moreover, there may be other techniques for generating p0.

## Claims

1. A transmitter comprising:
a first processing block configured to generate a first peak-reducing vector; and
at least one second processing block configured to
generate a second peak-reducing vector,
receive a first data vector, the first data vector comprising a plurality of samples, the first data vector having
a first peak with a first index and a first magnitude,
a second peak with a second index and a second magnitude that is less than the first magnitude, and
a first peak-to-average power ratio (PAPR), the second peak-reducing vector being generated from the first peak-reducing vector by shifting the first peak-reducing vector a number of samples equal to the second index, and

generate a second data vector having a second PAPR that is less than the first PAPR by adding the second peak-reducing vector to the first data vector.

2. The transmitter of claim 1, wherein the at least one second processing block is configured to use the first peak-reducing vector to reduce the first magnitude by a difference between the first and second magnitudes.

3. The transmitter of claim 1, wherein the first data vector also has a third peak with a third index and a third magnitude that is less than the second magnitude.

4. The transmitter of claim 3, wherein the first processing block is configured to use the first peak-reducing vector to generate a third data vector having a third PAPR that is less than the second PAPR.

5. The transmitter of claim 3, wherein the at least one second processing block is configured to reduce the first and second magnitudes each by a difference between the second and third magnitudes.

6. The transmitter of claim 1, wherein the first processing block is configured to reduce the first magnitude by combining the first data vector and the first peak-reducing vector.

7. A method of reducing peak-to-average power ratio (PAPR) in a signal, the method comprising:
operating a first processing block to generate a first peak-reducing vector from the signal; and
operating at least one second processing block to
generate a second peak-reducing vector,
receive a first data vector, the first data vector comprising a plurality of samples, the first data vector having
a first peak with a first index and a first magnitude,
a second peak with a second index and a second magnitude that is less than the first magnitude, and
a first peak-to-average power ratio (PAPR), the second peak-reducing vector being generated from the first peak-reducing vector by shifting the first peak-reducing vector a number of samples equal to the second index, and

generate a second data vector having a second PAPR that is less than the first PAPR by adding the second peak-reducing vector to the first data vector.

8. The method of claim 7, further comprising operating the at least one second processing block to use the first peak-reducing vector to reduce the first magnitude by a difference between the first and second magnitudes.

9. The method of claim 7, wherein the first data vector also has a third peak with a third index and a third magnitude that is less than the second magnitude.

10. The method of claim 9, further comprising operating the first processing block to use the first peak-reducing vector to generate a third data vector having a third PAPR that is less than the second PAPR.

11. The method of claim 9, further comprising operating the at least one second processing block to reduce the first and second magnitudes each by a difference between the second and third magnitudes.

12. The method of claim 7, further comprising operating the first processing block to generate the first peak-reducing vector using an Inverse Discrete Fourier Transform (IDFT) on a set of subcarriers that all have a magnitude of one.

13. A non-transitory computer-readable medium storing instructions for performing a method of reducing peak-to-average power ratio (PAPR) in a signal, the method comprising:
operating a first processing block to generate a first peak-reducing vector from the signal; and
operating at least one second processing block to
generate a second peak-reducing vector,
receive a first data vector, the first data vector comprising a plurality of samples, the first data vector having
a first peak with a first index and a first magnitude,
a second peak with a second index and a second magnitude that is less than the first magnitude, and
a first peak-to-average power ratio (PAPR), the second peak-reducing vector being generated from the first peak-reducing vector by shifting the first peak-reducing vector a number of samples equal to the second index, and

generate a second data vector having a second PAPR that is less than the first PAPR by adding the second peak-reducing vector to the first data vector.

14. The non-transitory computer-readable medium of claim 13, further comprising operating the at least one second processing block to use the first peak-reducing vector to reduce the first magnitude by a difference between the first and second magnitudes.

15. The non-transitory computer-readable medium of claim 13, wherein the first data vector also has a third peak with a third index and a third magnitude that is less than the second magnitude.

16. A system comprising:
a processor; and
a transmitter coupled to said processor and comprising
a first processing block configured to generate a first peak-reducing vector, and
at least one second processing block configured to
generate a second peak-reducing vector,
receive a first data vector, the first data vector comprising a plurality of samples, the first data vector having
a first peak with a first index and a first magnitude,
a second peak with a second index and a second magnitude that is less than the first magnitude, and
a first peak-to-average power ratio (PAPR), the second peak-reducing vector being generated from the first peak-reducing vector by shifting the first peak-reducing vector a number of samples equal to the second index, and

generate a second data vector having a second PAPR that is less than the first PAPR by adding the second peak-reducing vector to the first data vector.

17. The system of claim 16, wherein the at least one second processing block is configured to use the first peak-reducing vector to reduce the first magnitude by a difference between the first and second magnitudes.

18. The system of claim 16, wherein the first data vector also has a third peak with a third index and a third magnitude that is less than the second magnitude.

19. The system of claim 18, wherein the first processing block is configured to use the first peak-reducing vector to generate a third data vector having a third PAPR that is less than the second PAPR.

20. The system of claim 18, wherein the at least one second processing block is configured to reduce the first and second magnitudes each by a difference between the second and third magnitudes.