# Patent

## Title
Adaptive predistortion for a transmit system

## Abstract

Systems, methods, and devices relating to the provision of deliberate predistortion to an input signal to compensate for distortions introduced by an amplifier subsystem. An input signal is received by a signal processing system which includes a predistortion subsystem. The input signal is decomposed and the fragments are then predistorted by the predistortion subsystem by applying a deliberate predistortion to the fragments. The predistorted fragments are then separately processed and recombined to arrive at the system output signal. The predistortion subsystem adaptively adjusts based on characteristics of the system output signal. Also, the predistortion subsystem is equipped with a control system that is state based—the state of the predistortion subsystem is dependent upon the prevailing conditions and, when required, the control system switches the state of the predistortion subsystem. A feedback signal, a replica of the system output signal, is used in updating lookup table entries used to determine the predistortion.

## Background
The present application is a division of and claims priority to U.S. patent application Ser. No. 11/874,878, filed Oct. 18, 2007, which is a continuation of and claims priority to U.S. patent application Ser. No. 11/401,603, filed Apr. 10, 2006, now U.S. Pat. No. 7,295,066, which is a division of and claims priority to U.S. patent application Ser. No. 10/641,372, filed Aug. 13, 2003, now U.S. Pat. No. 7,026,871, and is a continuation-in-part of U.S. patent application Ser. No. 10/613,372, filed Jul. 3, 2003, now U.S. Pat. No. 7,068,101, all of which are incorporated by reference.


FIELD OF THE INVENTION
The present invention relates generally to power amplification systems and is specifically applicable but not limited to power amplification systems using a Chireix architecture.
BACKGROUND TO THE INVENTION
The recent revolution in communications has caused a renewed focus on wireless technology based products. Mobile telephones, handheld computers, and other devices now seamlessly communicate using wireless technology. One component that forms the core of such technology is the amplifier. Wireless devices require high efficiency amplifiers to not only extend the range of their coverage but also to conserve the limited battery power that such devices carry.
One possible architecture which may be used for such a power amplifier is called a Chireix architecture. Named after Henry Chireix who first proposed such an architecture in the 1930s, the Chireix architecture has fallen out of favor due to its seemingly inherent limitations. However, it has recently been revisited as it provides some advantages that other architectures do not have.
While the Chireix architecture provides some advantages, the process which the input signal undergoes also introduces some drawbacks. Specifically, distortions are introduced into the signal by the components in the Chireix based amplifier/modulator system. These distortions may also change over time and may therefore lead to a time-varying “drift” or change in the signal. Such distortions, time-varying or not, have led to problems that are not only inconvenient but expensive as well.
Based on the above, there is therefore a need for an amplifier system which provides the benefits of a Chireix based amplifier but which also compensates for or avoids the distortions which a Chireix based amplifier introduces. Such an amplifier system should adjust to differing conditions, preferably with little or no user intervention. It is therefore an object of the present invention to provide alternatives which mitigate if not overcome the disadvantages of the prior art.
SUMMARY OF THE INVENTION
The present invention provides systems, methods, and devices relating to the provision of deliberate predistortion to an input signal to compensate for distortions introduced by an amplifier subsystem. An input signal is received by a signal processing system which includes a predistortion subsystem. The input signal is decomposed and the fragments are then predistorted by the predistortion subsystem by applying a deliberate predistortion to the fragments. The predistorted fragments are then separately processed and recombined to arrive at the system output signal. The predistortion subsystem adaptively adjusts based on characteristics of the system output signal. Also, the predistortion subsystem is equipped with a control system that is state based—the state of the predistortion subsystem is dependent upon the prevailing conditions and, when required, the control system switches the state of the predistortion subsystem. A feedback signal, a replica of the system output signal, is used in updating lookup table entries used to determine the predistortion.
In a first aspect, the present invention provides a system for processing an input signal, the system comprising:

    
    
        a signal processing subsystem receiving and processing said input signal and producing a system output signal,
        an adaptive predistortion subsystem receiving at least two internal input signals and producing at least two predistorted signals by applying a deliberate predistortion to said at least two internal input signals,

wherein

        said predistortion subsystem distorts said internal input signals to compensate for distortions in said system output signal;
        said signal processing subsystem decomposes said input signal into separate components to produce said at least two internal input signals, each of said separate components being processed separately;
        said processing subsystem combines said predistorted signals after processing to produce said system output signal; and
        said deliberate predistortion applied to said at least two internal input signals by said adaptive predistortion subsystem to produce said at least two predistorted signals is adjusted based on characteristics of said system output signal and said input signal.
    
    


In a second aspect the present invention provides a method of processing an input signal to produce a system output signal, the method comprising:

    
    
        a) receiving said input signal;
        b) decomposing said input signal into at least two component signals;
        c) applying a deliberate predistortion to each of said at least two component signals to produce predistorted signals;
        d) combining said predistorted signals to produce said system output signal;
        e) adjusting said deliberate predistortion based on said characteristics of said system output signal.
    
    


In a third aspect the present invention provides an adaptive predistortion subsystem for use with a signal processing system which produces a system output signal, the predistortion subsystem comprising:

    
    
        determining means for determining a deliberate predistortion to be applied to an input signal;
        adjustment means for applying said deliberate predistortion to said input signal;
        update means for periodically updating said determining means based on characteristics of said system output signal.
    
    


In a fourth aspect, the present invention provides a preprogrammed control device for use in controlling an adaptive predistortion subsystem, said device being programmed with a predefined set of states, each state having associated with it a predefined set of commands to be executed by said subsystem when said device is in said state, the device comprising:

    
    
        a processor means for receiving and processing data relating to a status of said subsystem;
        first memory means for storing said data; and
        second memory means for storing preprogrammed settings for said device,

wherein

        when said device detects one of a set of specific, predetermined conditions, said device switches from one state to another.
    
    





BRIEF DESCRIPTION OF THE DRAWINGS
A better understanding of the invention will be obtained by considering the detailed description below, with reference to the following drawings in which:
FIG. 1 is a block diagram of a Chireix architecture amplifier subsystem;
FIGS. 2A and 2B illustrate how a vector can be decomposed in two different but similar manners;
FIGS. 3A-3D illustrate how different decomposition methods used on the same vector can produce different characteristics;
FIG. 4 illustrates a block diagram of an amplifier system which incorporates a predistortion subsystem according to one aspect of the invention and which also uses a feedback signal used by the predistortion subsystem;
FIG. 5 is a detailed block diagram of the predistortion subsystem illustrated in FIG. 4;
FIG. 5A is a block diagram of one type of control block which may be used in the predistortion subsystem of FIG. 4;
FIG. 5B is a state diagram illustrating the different states through which the predistortion subsystem of FIG. 4 may pass through;
FIG. 5C is a block diagram of a leaky integrator which may be used in the distortion monitor illustrated in FIG. 4;
FIG. 6 is a block diagram of a delay line circuit which may be used in the predistortion subsystem of FIG. 4;
FIG. 6A is a block diagram of a magnitude adjustment circuit which may be used in the predistortion subsystem illustrated in FIG. 5; and
FIG. 7 is a block diagram of a signal processing system according to another embodiment of the invention,



DETAILED DESCRIPTION
For clarity, the following terms are to be used with the following definitions:

    
    
        AM (amplitude modulation) refers to the AM of an RF (radio frequency) signal and is equal to the magnitude of the RF signal's complex base band equivalent
        PM (phase modulation) refers to the PM of an RF signal and is equal to the phase of the RF signal's complex base band equivalent.
    
    


Referring to FIG. 1, a block diagram of a Chireix architecture amplifier subsystem 10 is illustrated. A signal decomposer 20 receives an input complex baseband signal 30. Phase modulated RF signals 80A, 80B are produced after the decomposed output of the decomposer 20 are phase modulated by phase modulation circuitry 85A, 85B. These phase modulated signals 80A, 80B are received by power amplifiers 90A, 90B. The phase modulated signals are thus amplified by the power amplifiers 90A, 90B and are received by a signal combiner 100. The system output signal 110 (an RF signal corresponding to the input baseband signal 30) is output from the combiner 100 and is an amplified and modulated version of the input signal 30. Phase modulation of the phase modulated signals 80A, 80B is executed after the signal decomposer 20 separates input signal 30 into at least two components. These at least two components, after phase modulation, are the signals 80A, 80B.
As noted above, the Chireix architecture has been known to introduce distortions in the system output signal 110. Part of these distortions result from he decomposition and subsequent recombining of these components. The phase fragmentation circuitry in the decomposer 20 translates a complex signal sample x(k)=M(k)exp(jΘ(k)) to samples of a pair of phase signals α(k) and β(k). The individual phase signals are then translated into the complex signals a(k) and b(k):

a(k)=exp(jα(k))/2

b(k)=exp(jβ(k))/2

For each of the complex signals a(k) and b(k), the real part of the signal is mapped to the RF (radio frequency) in-phase channel and the imaginary part of the signal is mapped to the RF quadrature channel. For ease of reference and for use in later sections of this document, it should be noted that the expression exp(jα(k)) is termed a phasor and that the α(k) alone is termed a phase angle. When the complex variable exp(jα(k)) is represented by a pair of variables representing the real and imaginary parts, this will be termed a vector and will be denoted by ā(k).
Both signals a(k) and b(k), each having constant magnitude, are summed (by means of RF power amplification circuitry) by the combiner 100 to produce the system output signal 110. It should be noted that the factor 1⁄2 in both expressions a(k) and b(k) is a scaling factor which limits their sum to less than or equal to one (i.e. a(k)+b(k)≦1) on the assumption that x(k)≦1.
Two types of distinct decomposition and combination calculations may be used to relate the phase signals α(k) and β(k) to the complex signal x(k). For brevity, we denote x=M* exp(jΘ) and the sampled RF base band equivalent of the system output signal is denoted as c(k) which corresponds to x(k) such that c(k)=x(k) if there are no impairments in the system.
The first type of decomposition is termed Magnitude Linear (ML) Decomposition and it is defined by the following equations:



  
    
      
        Define
        ⁢
        
          
        
        ⁢
        φ
      
      =
      
        
          cos
          
            -
            1
          
        
        ⁡
        
          (
          M
          )
        
      
    
    ,
    
      
        then
        ⁢
        
          
        
        ⁢
        α
      
      =
      
        
          θ
          -
          
            φ
            ⁢
            
              
            
            ⁢
            and
            ⁢
            
              
            
            ⁢
            β
          
        
        =
        
          θ
          +
          φ
        
      
    
  




  
    
      ML
      ⁢
      
        
      
      ⁢
      
        Combination:
      
      ⁢
      
        
      
      ⁢
      c
    
    =
    
      
        
          exp
          ⁡
          
            (
            jα
            )
          
        
        +
        
          exp
          ⁡
          
            (
            jβ
            )
          
        
      
      2
    
  



Using trigonometric identities, it can be shown that |c|=M and ∠c=Θ as desired. Also, since M=cos(φ), the difference angle




  
    φ
    =
    
      
        β
        -
        α
      
      2
    
  
  ,




corresponds to the magnitude of the signal and the sum angle





  θ
  =
  
    
      α
      +
      β
    
    2
  




corresponds to the phase of the signal. For convenience, we define a phasor fragment opening angle Φ=2φ.

The second type of decomposition is termed Phase-Linear (PL) Decomposition. The Phase-Linear combination embeds a level of predistortion in the phasor fragmentation that is suitable for a combiner that provides an output magnitude that is substantially proportional to the opening angle, as opposed to a combiner that provides an output magnitude that is substantially proportional to the cosine of the opening angle. Such phasor fragmentation relieves the predistortion circuitry from linearizing the distortion effects due to the lack of the cosine characteristic in the combiner. It is believed that power amplifiers and combiners with such phase-linear characteristics have better power efficiency than those without such characteristics. Phase Linear Decomposition is defined by the following:



  
    
      
        PL
        ⁢
        
          -
        
        ⁢
        
          Decomposition:
        
        ⁢
        
          
        
        ⁢
        
          φ
          ′
        
      
      =
      
        
          π
          2
        
        ⁢
        
          (
          
            1
            -
            M
          
          )
        
      
    
    ,
    
      
        then
        ⁢
        
          
        
        ⁢
        α
      
      =
      
        θ
        -
        
          φ
          ′
        
      
    
    ,
    
      β
      =
      
        θ
        +
        
          φ
          ′
        
      
    
  




  
    
      PL
      ⁢
      
        -
      
      ⁢
      
        Combination:
      
      ⁢
      
        
      
      ⁢
      
        
        c
        
      
    
    =
    
      
        1
        -
        
          
            
              2
              ⁢
              
                φ
                ′
              
            
            π
          
          ⁢
          
            
          
          ⁢
          and
          ⁢
          
            
          
          ⁢
          ∠
          ⁢
          
            
          
          ⁢
          c
        
      
      =
      
        θ
        .
      
    
  



Similar to the ML Decomposition, the phasor fragment opening angle is defined as Φ=2φ′.
As a further refinement of the above decompositions, it should be noted that there are two possibilities for each decomposition. As is well-known, a complex number can be represented by a vector and, as shown in FIGS. 2A and 2B, this vector can be decomposed in two ways. Both FIGS. 2A and 2B illustrate how a vector x can be decomposed as the sum of two equal length vectors ā and b in two different, and initially equivalent ways —a summation triangle can be formed to the left (clockwise) of x (FIG. 2A) or to the right (counterclockwise) of x (FIG. 2B).
This realization allows for the flexibility of selecting between the left or the right decomposition for any sample of x(k). Such a flexibility provides for some reduction of the power spectral density of the phase signals α(k) and β(k). Such a reduction is desirable as it reduces out-of-band emissions and in-band distortions.
The distortions for which the predistortion subsystem is to compensate may come as a phase distortion, a magnitude distortion, or as a combination of both. It has been found that, without predistortion, the system output signal 110 has an amplitude modulation (AM) that is not equal to the expected and desired AM. Furthermore, the phase modulation (PM) of the system output signal 110, if predistortion is not present, deviates from the expected and desired PM. Experiments have found that the AM distortion or error (magnitude distortion) depends on the AM of the input signal. Also, it has been found that the PM distortion (or phase distortion) depends on the AM of the input signal.
It should be noted that the predistortion modification, defined as any deliberate distortion which has been applied or is to be applied to the input signal to change at least one original characteristic of the input signal, can take many forms. Two specific types of predistortion, phase predistortion and magnitude predistortion are currently envisioned although other types are possible. These two types, separately or together, can make up the predistortion modification. In some applications, only a magnitude type predistortion modification may be required while in others only a phase type predistortion is required.
One possible source of the AM/Am and AM/PM distortion is the gain and phase imbalance between the phasor fragments.
Referring to FIGS. 3A and 3B and as explained above, a vector x can be decomposed to the sum of vectors ā and b in two different, and initially equivalent manners. FIG. 3A illustrates the left or clockwise decomposition while FIG. 3B illustrates the right or counterclockwise decomposition. Both summations lead to ā+ b= x.
Phasor fragmentation operates under the assumption that at the point of summation the ratio between the lengths of the vectors (the magnitudes of the phasors) equals the ratio at the point of decomposition. Typically that ratio equals one, implying that both vectors are to be of equal magnitude at all times.
Any common rotation of both vectors results in a common modification of the angles and this implies a rotation of the sum x. If the common rotation is constant over time, then the resulting phase rotation of the decomposed signal is not a non-linear distortion.
Branch imbalance between vectors ā and b consists of a magnitude difference (gain imbalance) and rotation difference (phase imbalance) between the two. As an example, in FIGS. 3C and 3D, vector ā is rotated to ā′ and while vector b remains the same. The sum is now not only caused to rotate, but the length of vector x is altered, forming a different vector x′. Unfortunately, the effects of such phase imbalance is twofold: it implies magnitude distortion of the sum depending on the magnitude of the decomposed vector (AM/AM), and a constant (linear) phase distortion.
Similarly, a gain imbalance between the two vectors implies a non-linear magnitude (AM/AM) and non-linear phase (AM/PM) distortion depending on the magnitude of the decomposed vector.
Also, the effect of an imbalanced ā′ on the sum x′ depends on the geometrical orientation of the decomposition triangle. If the left decomposition triangle is elected, a different distortion results compared to when the right decomposition triangle is elected. As a consequence of this, predistortion of such non-linear effects requires knowledge of the orientation of the decomposition triangle, and for either left or right orientation, a different pair of predistorting functions is required. These predistorting functions are denoted as the Left AM/AM predistortion function FL(M), the Left AM/PM predistortion function GL(M), the Right AM/AM predistortion function FR(M), and the Right AM/PM predistortion function GR(M). Applying a predistortion with phase and magnitude predistortion that differentiates between left and right decomposition shall be termed dual-predistortion.
Since the selection of the predistortion functions FL(M) or FR(M), and GR(M) or GL(M), depends on the orientation of the decomposition triangle for each vector x corresponding to the complex data x(k), predistortion is ideally performed after decomposition, when the elected orientation is known.
AM/AM magnitude predistortion is accomplished through modification of the opening angle Φ, and AM/PM phase predistortion is accomplished through modification of the resultant phase Θ. Consequently, predistortion can be accomplished by adjusting the angles α(k) and β(k) of the phasors, in replacement of adjusting the magnitude and phase of x(k). Two pairs of functions AL(M), BL(M) and AR(M) and BR(M), which adjust the angles α(k) and β(k), are defined to additively conform

α′(k)=α(k)+AL(M(k))

β′(k)=β(k)+BL(M(k)

or

α′(k)=α(k)+AR(M(k)

β′(k)=β(k)+BR(M(k)

for the left and right decomposition triangles respectively.

The common contribution of the left predistorting adjustment is




  
    
      
        A
        L
      
      ⁡
      
        (
        M
        )
      
    
    +
    
      
        B
        L
      
      ⁡
      
        (
        M
        )
      
    
  
  2




and it modifies the phase of the resultant vector based on the magnitude of x(k). The differential contribution of the left predistorting adjustment is





  
    
      
        A
        L
      
      ⁡
      
        (
        M
        )
      
    
    -
    
      
        B
        L
      
      ⁡
      
        (
        M
        )
      
    
  
  2




and it modifies the magnitude of the resultant vector based on the magnitude of x(k). The same holds for the right adjustments respectively.

The angle adjustment functions AL(M) and BL(M) and AR(M) and BR(M) may be implemented by look up tables (LUTs) with linear interpolation.
As noted above, one solution to the above issues is to predistort the input signal. Further details on this approach can be found in co-pending application entitled Predistortion Circuit For a Chireix Power Amplifier Transmit System and filed with the US Patent and Trademark Office, the whole of which is incorporated herein by reference. While the predistortion solution does work, it is not as robust and as fault tolerant as may be desirable. An adaptive predistortion subsystem 200, as illustrated in FIG. 4, would compensate for changing conditions and for other distortions which the system output signal may have.
Referring to FIG. 4, a block diagram of such an adaptive predistortion subsystem 200 is illustrated as being part of a signal processing system.
While an analog implementation of the subsystem 200 is possible, it has been found that a digital implementation was simpler to achieve. As such; the following description assumes that the input signal 30 is a digital signal having a digital representation of the desired AM and its PM of the output RF signal. Digital AM/AM predistortion modifies the magnitude of the complex digital input signal such that the RF output signal has the desired AM, despite the distortion. Digital AM/PM predistortion modifies the phase of the complex digital input signal such that the RF output has the desired PM, despite the distortion.
As can be seen in FIG. 4, multiple components are involved in the signal processing system which incorporates an adaptive predistortion subsystem. The input signal x(t) 30 is received by the decomposer 20. The decomposer 20, containing a phasor fragmentation engine 25 procudes the signal components α 45A and β 45B along with the data regarding whether left or right decomposition was used. These three signals are received by the adaptive predistortion subsystem 200 which contains the α lookup table (LUT) 230 and the β lookup table 250. The subsystem 200 also receives a feedforward 203 of the input signal x(t) along with the feedback signals 320, 310 from the feedback processing block 220. The feedback processing block 220 receives a replica of the system output signal 110 from the combiner 100. The feedback processing block 220 feeds the polar equivalent of the system output signal 110 to the predistortion block 200 as the magnitude 310 of the system output signal z(t) and the phase 320 of the signal z(t).
The predistortion subsystem predistorts the α component 45A and the β component 45B based on the decomposition data 55 and the input signal x(t) 30. The predistorted a component (denoted by α′ 47A) and the predistorted β component (denoted by β′ 47B) are received by the phase modulation blocks 85A, 85B respectively.
Each lookup table 230, 250 consists of two internal tables—one for left decomposition (201A, 201B) and one for right decomposition (202A, 202B) Table 201A details the values to be used for the received α values for a left decomposition while table 201B details the values to be used for the received β values for a left decomposition. Similarly, table 202A contains the values to be used for the received α values for a right decomposition while table 202B contains the values to be used for the received β values for a right decomposition. As an example, if table 201A in LUT block 200A has a phase adjustment value of 0.4 at an input magnitude of 0.5, then if the magnitude value received by LUT block 200A is 0.5 and if the α value received is 0.2 and left decomposition was used, the adjustment value of 0.4 is added to the α value of 0.2 to result in the adjusted α value of 0.6 as the output of the LUT block 200A. This value is found from table 201A due to the left decomposition. If the decomposition was a right decomposition, then table 202A would have been used to find the corrective value for the received α value.
Similar to the above, table 201B in LUT block 200B is used to obtain a predistortion adjustment for the received β value. Since the decomposition was a left decomposition, then table 201B was used. If the decomposition was a right decomposition, then table 202B would be used for the corrective value for the received β value.
The adaptive predistortion subsystem 200 receives the α and β components and predistorts these components based on the different inputs such as the feedforward signal 210, the feedback signals 310, 320, and the decomposition information 55. The adaptive predistortion subsystem 200 is illustrated in more detail in FIG. 5.
It should be noted that the predistortion modification, defined as any deliberate distortion which has been introduced to the input signal to change at least the phase or magnitude of the input signal, can take many forms, two of which are phase predistortion and magnitude predistortion. These two types, separately or together, can make up the predistortion modification. In some applications, only a magnitude type predistortion modification may be required while in others only a phase type predistortion modification is required. In the embodiment explained here, the two types of predistortion, together comprise the predistortion modification.
Referring to FIG. 5, a more detailed block diagram of the adaptive predistortion subsystem 200 is illustrated. The α LUT 220 and the β LUT 250 each receive the magnitude of the input signal 30 (|x(t)|). Each LUT (220 and 250) also receives the decomposition information 55. As noted above, it has been found that both AM and PM distortions depend upon the AM of the input signal. Based on these inputs, the α and β LUTs 220, 250 each produce a distortion modification value that additively modifies the α and β values (45A and 45B) through adders 255A, 255B. The additively modified α and β values are output as α′ (47A) and β′ (47B) signals. A control block 600 contains and executes the control software which handles exceptions and initialization of the subsystem. This control software will be discussed later in this document.
The adaptive nature of the subsystem 200 stems from the ability of the subsystem to adapt to changing conditions—the entries in the LUTs 220, 250 are adjustable based on prevailing conditions. This is accomplished by the use of update blocks 230, 260. The update blocks 230, 260 receive characteristics of the input signal 30 (by way of delay blocks 240, 270) along with characteristics of the feedback signal 345 and, based on these two inputs, updates the LUTs 220, 250 by way of the update block 256. The phase update block 260 receives the phases of both the feedback signal and the input signal and determines phase error and phase updates. The magnitude update block 230 receives the magnitudes of both the input signal and the feedback signal and determines magnitude error along with magnitude updates. These updates, both phase and magnitude, are jointly translated into joint α and β values by the update block 256 to be used by the respective LUTs (220, 250) as updates.
Part of the input to the update blocks 230, 250 are the characteristics of the input signal x(k). The magnitude 290 and the phase 300 of this input signal are received from the conversion unit 210. As noted above, this conversion unit converts the feedforward signal 203 from Cartesian coordinates to polar coordinates. The conversion unit 210 is convenient and makes the predistortion easier to accomplish. As is well known, signal representations using Cartesian coordinates take the form of z=xαj y where x and y are the real and imaginary components. Polar coordinates take the form of z=Aejφ where the magnitude of the signal is A and its phase is φ. Since both the magnitude and the phase of the signal is to be modified by the predistortion subsystem, it is clearly more convenient to accomplish this if the signal were in polar coordinates. Again as is well known, A=(x2+y2)1/2 while φ=tan−1 (y/x). Once the signal has been converted into polar coordinates, these characteristics of the input signal can be directly used by the LUT update blocks.
To account for changing conditions and to acquire appropriate LUT entries, a feedback mechanism is employed to adjust or adapt the lookup table entries in lookup table blocks 220, 250. Delay blocks 240, 270 ensure that the feedback sample is mated with the proper value of the input signal waveform when processing and updating the lookup table entries in lookup table blocks 220, 250. Delay block 70 also performs a delay line that delays the decomposition information 55.
The lookup tables 220, 250 internal to the predistortion subsystem 200 contain values to be added to the received α or β values based on the other inputs (left or right triangle decomposition and input signal feed forward information) with the output being the predistorted signal value. Each lookup table block 200A, 200B contains dual lookup tables a table for left decompositions (201A, 201B), and a table for right decompositions (202A, 202B). Table 201A details the values to be used for the received α values for a left decomposition while table 201B details the values to be used for the received β values for a left decomposition. Similarly, table 202A contains the values to be used for the received α values for a right decomposition while table 202B contains the values to be used for the received β values for a right decomposition. As an example, if table 201A in LUT block 200A has a phase adjustment value of 0.4 at an input magnitude of 0.5, then if the magnitude value received by LUT block 200A is 0.5 and if the α value received is 0.2 and left decomposition was used, the adjustment value of 0.4 is added to the α value of 0.2 to result in the adjusted α value of 0.6 as the output of the LUT block 200A. This value is found from table 201A due to the left decomposition. If the decomposition was a right decomposition, then table 202A would have been used to find the corrective value for the received αvalue.
Similar to the above, the table 201B in LUT block 200B is used to obtain a predistortion adjustment for the received β value. Since the decomposition was a left decomposition, then table 201B was used. If the decomposition was a right decomposition, then table 202B would be used for the corrective value for the received β value. The β lookup table block 250, based on the value of |x(t)| and on the other inputs, determines the proper amount of adjustment and adds the adjustment to the input β value 45B by way of adder 225B to result in the desired predistorted signal fragment.
While the above described LUT performs an additive adjustment to the received α or β values, other implementations are possible. Instead of an additive adjustment value, the LUT may also provide a multiplicative adjustment to the received value.
As noted above, the adaptive predistortion subsystem 200 in FIG. 5 is adaptive in that the values of the lookup table entries in the lookup table (LUT) blocks 220, 250 change over time to account for changing conditions or for acquiring appropriate LUT entries. This adaptability is implemented by way of a feedback signal 345 tapped from the system output signal 110. The magnitude 310 of the feedback signal 345 and the phase 320 of the feedback signal 345, are received by the predistortion subsystem 200. These two values result from processing of the system output signal 110 by the feedback processing block 220 which contains an analog/digital (A/D) converter 330, a Cartesian to polar conversion unit 340, a demodulation module 335, and a filtering module 337. As can be seen in FIG. 4, the system output signal (an analog signal) is tapped and this tapped signal 345 is received (after passing through the filtering module 337 and the demodulation module 335) by the A/D converter 330 for conversion from an analog to a digital signal.
After conversion to a digital signal, the feedback signal is converted from Cartesian to polar coordinates by the conversion unit 340. The tapped signal 345 is thus represented by the two values—the magnitude 310 of the feedback signal 345 and the phase 320 of feedback signal 345.
An automatic gain control (AGC) module 355 and a multiplier 357 combine to remove any gain in the magnitude feedback signal 310 as compared to the magnitude of the input signal 30. Similarly, a synchronization circuit 365 and an adder 367 adjust the incoming phase feedback signal 320 as compared to the phase 290 of the input signal 30.
Once the two digital feedback signals are received, they are then compared with the delayed feedforward signal coming from the delay blocks 240, 270. The updated values for the LUT entries are then calculated and entered into their respective lookup tables. It should be noted that the comparison may be done by subtracting the phase of the feedback signal from the phase of the input signal and by subtracting the magnitude of the feedback signal from the magnitude of the input signal.
To further elaborate on the above process, the update process is dependent on the difference between the tapped system output signal 345 and the input signal 30. This difference is, of course, taken after both signals are in polar coordinates. The magnitude and phase errors are defined as:

em(k)=|z(k)|−|x(k)|

eφ(k)=(∠z(k)−∠x(k))

where

    
    
        em=magnitude error
        eφ(k)=phase error
        z(k)=magnitude of feedback signal (signal 310)
        x(k)=magnitude of input signal (signal 300)
        ∠z(k)=phase angle of feedback signal (signal 320)
        ∠x(k)=phase angle of input signal (signal 290)
    
    


The lookup table entries in the LUTs are modified or updated based on the magnitude M=|x(k)|. Based on the magnitude value, one, two, or no entries in an LUT are updated as follows:

    
    
        1. Obtain the update quantity δ
        2. Identify the LUT entry addresses n=└λM┘ and n+1, calculate the interpolation distance s=λM−n.
        3. Update table entries as follows:

Tn(k+1)=Tn(k)+(1−s)·δT iff 0 ≦n≦(N−1)

Tn+1(k+1)=Tn+1(k)+(s)·δT iff 0 ≦n+1≦(N−1)

    
    


From the conditions in the above equations, it is clear that two, one, or no entries are updated depending on the value of λM. The product λM will be called the virtual address, since it decomposes into two physical addresses (indices) for the LUT and an interpolation weight. Specific to the predistortion adaptation process, the LUT update is based on a magnitude and phase error, defined as follows:

eM(k)=|z(k)|−|x(k)|

eφ(k)=(∠z(k)−∠x(k))

The LUTs for AL(M), BL(M) (for a left triangle decomposition) or AR(M), BR(M) (for a right triangle decomposition) are updated as follows:

    
    
        1. Define a magnitude update quantity δF=μF·eM(k) based on the magnitude error, which depends on the difference of magnitudes and an update speed parameter 0≦μF<1.
        2. Define a phase update quantity δG=μG·eφ(k) based on the angular error (radians) which depends on the difference of angles and an update speed parameter 0≦μG<1
        3. Determine whether z(k) originates from a ‘Left’ or a ‘Right’ decomposition triangle. This is accomplished by noting for each transmitted x(k) whether it is decomposed into a ‘Left’ or a ‘Right’ triangle.
        4. If the decomposition triangle is ‘Right’:
        
            Update AR(M) using δT=−δF−δG 
            Update BR(M) using δT=δF−δG 
        
        
        5. If the decomposition triangle is ‘Left’:
        
            Update AL(M) using δT=δF−δG 
            Update BL(M) using δT=−δF−δG 
        
        
    
    


In practice, AM/AM and AM/Pm distortions by the power amplifier/combiner are not the only contributors to e(k). However, gradual updating through proper selection of μF and μG allows the process to extract the distortions of interest by averaging out other distortion effects and noise.
As a general example, consider a LUT BR(M) that has been designed with N=6 entries for M≦1. It follows that α=5, allowing the full magnitude range to be mapped to the full address range of the LUT. For an update at x(k)=0.35·exp(j·2.7), we have M=0.35. The virtual LUT address is 1.75 and the LUT entries of interest are addressed at n=└1.75┘=1 and n+1=2. Interpolation, is performed with s=1.75−1=0.75. If the corresponding sensed sample is z(k)=0.2·exp(j·2.3), then the phase error is eφ(k)=−0.4 and the magnitude error is eM(k)=−0.15. Let μF=0.1 and μG=0.1. From the magnitude and phase errors, it follows that the updates equal δG=−0.04 and δF=−0.015, resulting in δT=0.015+0.04=0.055. A weighted correction of 25%·0.055≈0.01 is applied to T1 and a complementary weighted correction of 75%·0.055≈0.04 is applied to T2. This update is illustrated in the table below:
















Address
LUT content

LUT content



(n)
before update
Correction
after update





















0
−1
0
−1



1
2
0.01
2.01



2
−0.5
0.03
−0.47



3
−0.5
0
−0.5



4
0.5
0
0.5



5
2
0
2










To protect against instability under certain distortion conditions, some exceptions have to be made. The first of these is the LUT entry at address 0. The entry should be optionally frozen while allowing other LUT entries to be updated. This is done by setting (following the formula for Tn(k+1))

T0(k+1)=T0(k) if n=0

T1(k+1)=T1(k)+(s)·δT iff n+1=1

if a freeze is desired. Such a freeze has been found to allow for added subsystem robustness if residual branch imbalances cause LUT update instability for near-zero magnitude values.

The other exception is to limit the range of amplitude values for an LUT update. By default, the magnitude range is 0≦M≦1 but LUT updates should only occur for the magnitude range M0≦M≦M1. The definition of this operative range is not dependent on the magnitude values corresponding to the LUT entries. LUT entries representing magnitude values beyond the limited range are still affected by the updates due to the use of interpolating adaptation. Limiting the active magnitude range has been found to improve stability, especially in the case of branch imbalance where phase distortion at near zero magnitude is severe. It has also been found to protect against instability for near-maximum magnitude values caused by inadvertent flipping due to pre-distortion.
A flip guard may also be used to protect against inadvertent flipping. As noted before, a further opening of the angle between the phasor fragments increases the magnitude of the output signal while closing the angle decreases the magnitude. If the decomposition is a right triangle decomposition, then opening the angle is done by increasing α and decreasing β. Conversely, for a left triangle decomposition, opening the angle is done by decreasing α and increasing β. This method only holds true until the magnitude is at its maximum and, after this, the modifications to α and β accomplish the exact opposite—“opening” the angle is effectively closing the angle and thereby leads to instability of the adaptation process. In such an instance, the fragments are “flipped” and what was thought as a right decomposition is really a left decomposition and vice versa.
To guard against such an eventuality, a guard flag for a table entry can be set after the application of predistortion to a magnitude corresponding to the table entry if the fragments are flipped. When the guard flag for an individual table entry is set, that table entry is no longer updated. The flag is cleared once predistortion for a magnitude corresponding to that table entry no longer flips the fragments. This is possible for an individual table entry, even when that entry is frozen, because of the interpolating effects from surrounding table entries. Defining a left guard flag as GL(n) and a right guard flag as GR(n), the guards are set and cleared as follows:

    
    
        For right decomposition (α<β): if a′>β′ for the predistorted fragments, then set the guard GR(n)=1, otherwise, clear the guard GR(n)=0.
        For a left decomposition (α<β): if α′<β′ for the predistorted fragments, then set the guard GL(n)=1, otherwise, clear the guard GL(n)=0.
    
    


To examine the guards during adaptation at a particular magnitude 0≦M<1, the following process is used:

    
    
        For a right update: if δF<0 and GR(n)=1, force δF=0 (i.e. perform the LUT update only using the phase update δG)
        For a left update if δR<0 and GL(n)=1, force δF=0 (i.e. perform the LUT update only using the phase update δG)
    
    


From the above discussion, it should be fairly clear that two of the values underpinning most of the calculations are the values for em(k) and e100(k), the magnitude and phase errors. While the equations above for these two values will provide acceptable answers, better results may be had by taking into account the phase, magnitude, and delay effects of the feedback signal (i.e. the tapped system output signal 345). However, to be able to do this, the input signal must be properly delayed so that samples from the interpolated input waveform, as obtained from the delayed input signal samples, are mated with the relevant system output signal sample. Such proper delaying should therefore take into account most, if not all, the time delay involved in the processing production, and feedback of the system output signal (round trip delay). This round trip delay is denoted as t (seconds) and, before the LUT updating begins, the delay blocks 240, 270 acquires the delay and delays the input signal accordingly so as to time-align the input signal samples with the incoming feedback signal.
To implement this delay, a delay line is used with a depth of K, meaning K samples of the input signals may be stored in the delay line. As should be clear, each of the K samples were sampled at different instances in time. The value of K is predetermined and should be enough to allow for the maximum possible path delay between the input signal and the feedback signal. These delays are due to a combination of any of the following: digital pipelining, analog and digital filter group delays, analog propagation delays, and the system and implementation dependent delays.
Because of this delay, a time delayed version of the input signal, xδ(k) is defined and this is ideally

xδ(k)=x(k−δ)

    
    
        where
        
            δ=τ·FS 
            FS=signal sampling rate
            τ=delay (normal trip delay between input and system output signal feedback)
        
        
    
    


To obtain a better result for xδ(k), linear interpolation is used to allow for fractional values of δ. Thus, the delay is divided into two parts κ, the integer part of the sample (representing a discrete sample delay at the sample rate FS), and ν, a fractional sample of the delay.
Using this notation, the delayed portion of the input signal can be represented as:

xδ(k)=(1−ν)·x(k−κ)+ν·x(k−κ+1)

    
    
        where
        
            δ=τ·FS 
            κ=└δ┘
            ν=δ−κ

As can be seen, for an integer δ, xδ(k)=x(k−δ).

        
        
    
    


The above Cartesian equation can be applied to polar representations by having separate delay lines for magnitude (|xδ(k)|) and phase (∠xδ(k)) using the sequences |x(k)| and ∠x(k). These are given by

∠xδ(k)=(1−ν)·∠x(k−κ)+ν·∠x(k−κ+1)

|xδ(k)|=(1−ν)·|x(k−κ)|+ν·|x(k−κ+1)|

It should be fairly clear that xδ(k) is calculated from the samples x(k), x(k−1), x(k−2), . . . x(k−κ), samples of the input signal taken at time k, k−1, k−2, . . . k−κ.

The above equations for ∠xδ(k) has a peculiarity that is due to the way angle values work. Since ∠x(k−κ) and ∠x(k−κ+1) are represented by modulo 2π radians (360 degrees) and since −π≦∠x(k)≦π, then errors could easily occur.
Thus, if −π≦∠x(k)≦π, and if |∠x(k−κ+1)−∠x(k−κ)|≧π, then

∠xδ(k)=(1−ν)·∠x(k−κ)+ν·(∠x(k−κ+1)+2π) if ∠x(k−κ+1)≦∠x(k−κ)

∠xδ(k)=(1−ν)·∠x(k−κ)+ν·(∠x(k−κ+1)−2π) if ∠x(k−κ+1)≦∠x(k−κ)

The above described delay can be implemented by cascaded delay elements and associated sample storage: FIG. 6 illustrates such a delay subsystem which can be used as delay blocks 240, 270. As can be seen, delay elements 242A, 242B, 242C, 242D, 242E are cascaded and provide delays and storage for input signal samples 244A, 244B, 244C, 244D, 244E. Switches 245A, 245B, 245C, 245D, 245E allow any one of the signal samples 244A-244E to be switched so that it can be used. These samples 244A-244E can be weighted accordingly by programmable weighting blocks 246A, 246B, 246C, 246D, 246E. The weighted samples are then summed up by adder 248 to produce the delayed signal 249 to be used by the system. The switches 245A . . . 245E and the values in the weighting blocks 246A . . . 246E may be user/system controllable so that any combination of weighted samples may be produced.
As an example, if τ·FS=2.4 samples is required, then a value of 0.6 is used by the weighting block 246C and a value of 0.4 is used by the weighting block 245D. Then, by closing switches 245C and 245D then the sample x(k−2.4) is obtained.
The phase correction required for the feedback portion of the system output signal is also dependent on the delayed signal xδ(k). The complex phase difference between the delayed signal xδ(k) and the system output signal z(k), represented by γ, is due in part to the round trip delay and to other factors. The phase correction and delay can be found by executing the following process:

    
    
        1. Reset the running-maximum register: set Pmax=0
        2. Reset the phase correction: set γ=0
        3. Preset the variable interpolating delay line (discussed above as a delay subsystem): set δ=δ0 (i.e. κ=κ0 and ν=ν0)
        4. Perform a coarse search: select only interpolating delay line integer sample-delay increments of δ (maintain ν=ν0, increment κ only): set the coarse and fine delay increments to Δκ=1 and Δν=0.
        5. Start with the first correlation window n=0. Time instance k=0 signifies the commencement of the search.
        6. “Integrate”: calculate the inner product
    
    






  
    P
    δ
  
  =
  
    
      A
      MAC
    
    ·
    
      
        ∑
        
          k
          =
          nW
        
        
          nW
          +
          W
          -
          1
        
      
      ⁢
      
        exp
        ⁡
        
          (
          
            j
            ·
            
              (
              
                
                  ∠
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      x
                      δ
                    
                    ⁡
                    
                      (
                      k
                      )
                    
                  
                
                -
                
                  ∠
                  ⁢
                  
                    
                  
                  ⁢
                  
                    z
                    ⁡
                    
                      (
                      k
                      )
                    
                  
                
              
              )
            
          
          )
        
      
    
  




    
    
         This is a complex multiply-and-accumulate (MAC) covering a time window of W current samples. Note that the above notation strictly implies that the correlation windows n=0, 1, 2, . . . are adjoined, yet this is not an operational requirement. A spacing between windows translates to a lengthening of the search.
        7. “Dump”: compare the complex magnitude |Pδ| with the value Pmax, and if |Pδ|>Pmax then
        
            (a) update the maximum: set Pmax=|Pδ|
            (b) update the delay: retrieve the corresponding value of δ and store it as δmax and
            c) update the angle: if fine searching store the angle γ=∠Pδ
        
        
        8. Proceed to the next window: set n=n+1, and increment δ by Δκ+Δν
        9. Repeat from Step 4, until the full integer (coarse) delay range κ=κ0 . . . Kν of the interpolating delay line is exhausted. (Kν is the maximum variable delay, κ0 is the starting coarse delay value for the search)
        10. Store δmax as δ1.
        11. Read the programmed fine search step size Δν.
        12. Repeat Steps 1 and 4 to 8 for the fine search: employ fractional-sample increments of Δν, covering the delays δ1−1≦δ≦δ1+1, ensuring not to exceed 0≦δ≦Kν
        13. Store δmax as δ2.
        14. Set and freeze the interpolating delay line delay to δ2.
        15. Apply the phase correction γ
    
    


Once the full course range is exhausted, the stored δ value δmax is used as a starting point for a fine search. With Δκ=0, the Δν is incremented fractionally so as to search the delays δ1−1≦δ≦δ1+1. Once found (using an analogous loop to that used for the coarse search), the phase correction γ is found along with the proper delay value. It should be clear that γ=∠Pδ.
Updating of the phase LUT is only to be accomplished after the phase synchronization of the feedback signal. Freezing the phase LUT update can be done by setting μF=0 and μG=0.
After the phase correction is done, the magnitude of the feedback signal z(k) has to be adjusted as well. This adjustment is required to compensate for any gain in the system output signal relative to the input signal. Magnitude adjustment is accomplished by directly multiplying the feedback signal z(k) by a factor ASX(k). This is found using the following formula:

ASX(k+1)=ASX(k)+μA·(|xδ(k)|−ASX(k)·|z(k)|)

For kA≦k≦kA+WA−1

where

    
    
        kA=commencement time for the magnitude adjustments (preferably as early as possible)
        WA=number of samples for which the adjustments are performed
        μA=update step size and 0≦μA≦1

ASX(kA), μA, and WA can be programmable.

    
    


Referring to FIG. 6A, a circuit for implementing the magnitude adjustment is illustrated. As can be seen, the magnitude 400 of the portion of the system output signal is multiplied with the factor ASX(k) by way of multiplier 410. The result is to be used in calculating the desired predistortion modification. This result is also subtracted from the magnitude (|xδ(k)|) of the delayed signal by way of adder 420. The result of the subtraction is multiplied by the update step size μA by way of multiplier 430 and added to the factor ASX(k) by adder 440 to produce the next value in the sequence for the factor. The delay element 450 delays the resulting value until it is ready for use.
The feedback signal (the tapped system output signal 345 which is a replica of the system output signal 110) may be further processed to achieve better results. As an example, the gain and/or phase of this replica of the system output signal may be adjusted for better coordination with the delayed replica of the input signal.
While the above embodiment utilizes LUT entries that are adjusted based on the difference between the input signal and the system output signal, the system output signal alone may be used, if desired, in adjusting the L′UT entries. This alternative would require the modulation of the system output signal and the regeneration of a signal similar to the input signal. Successful recovery is possible if the amount of distortion is restricted such that demodulation can occur without errors. This approach is not desirable if the input signal is readily available to be used in determining the parameters of the adjustments. Another alternative involves measuring the distortion by measuring the amount of out-of-band energy. Adjustments can then be made to minimize the amount of out-of-band energy and, consequently, minimize the distortion.
Referring to FIG. 5, a distortion monitor 500 is illustrated as receiving input from both update blocks 230, 260. The distortion monitor 500 monitors the distortion by way of the error functions defined as:

eM(k)=ASX·|z(k)|−|xδ(k)|

eφ(k)=(γ+∠z(k)−∠xδ(k))

where γ is the complex phase difference between xδ(k) and z(k). By monitoring the distortion, the predistortion subsystem can detect LUT divergence. If such a condition is detected, a full or a partial re-initialization may be in order.

The error signals given above are weighted as:

η(k)=s·|eM(k)|+(1−s)|e100(k)|

with 0≦s≦1. The averaging of η(k) can be done using a leaky integrator

E(k)=λ·η(k)+(1−λ)(η(k−1)

where the update parameter is λ.

The leaky integrator may be implemented as in the block diagram in FIG. 5C where λ−2−n which reduces the two multiplication operations into two simpler shift functions and one add function. This leaky integrator is a first order infinite impulse response low pass filter.
The error functions given above in relation to the distortion monitor many be used in place of the previous error functions for better results.
The control software for the above described predistortion subsystem allows for different actions to be taken for different conditions. The software can be viewed as a state machine with each state allowing the subsystem to execute specified actions. Thus, when the subsystem enters a specific state, it is equivalent to running a subroutine as specific actions and tasks are executed by the subsystem.
The control software and its analogous state machine can therefore be viewed as a control system that determines the actions of the predistortion subsystem. The various states and the steps taken when the predistortion subsystem enters it are listed and explained below:

    
    
        State: Reset [RST]
        The reset state resets all variables in the subsystem to their initial values and commences the synchronization and adaptation of the LUTs without regard to previous results or actions.
    
    


This reset state is activated through a reset external to the subsystem or through a self reset condition. The actions taken for a reset are, in order:

    
    
        1) All LUT entries are reset to an all-zero state and the update function is disabled.
        2) The predistortion subsystem is disabled and the signals bypass the additive modification such that (referring to FIG. 5) α′=α and β′=β.
        3) The LUT entries are loaded with preset values or, alternatively, maintain all-zero entries for all LUT entries and disable the update function.
        4) Preset the automatic gain control module 355 to a specific, predetermined and preprogrammed value (such as 1) to prevent further adjustments to the feedback magnitude signal 310.
        5) Preset the synchronization circuit 365 to a specific, predetermined and preprogrammed value (such as zero) to prevent further adjustments to the feedback phase signal 320.
        6) Preset the distortion monitor 500 with a preprogrammed starting value by loading the leaky integrator with a predetermined value. Alternatively the distortion monitor may be set to zero.
        7) Preset the delay line to its starting value.
        8) Reset the initialization duty control (a control value to be discussed later).
        9) Reset the update duty control (a control value to be discussed later).
        10) Reset the distortion watch dog (to be discussed further below).
        11) Disable the feedback from the system output signal.
    
    


When in the reset state, the predistortion subsystem 200 executes all of the actions listed above. During this time period, the predistortion subsystem 200 is bypassed as noted in the second item above.
State: Pause [PSE]
During the Pause state, all the functions of the predistortion subsystem 200 are halted and all the settings, values, variables, and entries in the subsystem are frozen and feedback is off.
State: Hold [HLD]
The Hold state is analogous to the Pause state with the exception that the subsystem, except for feedback, is functioning. The predistortion subsystem 200 freezes all the variables and LUT entries and turns off the feedback from the system output signal which still predistorting the α and β signals. However, the α and β LUT entries are not being adjusted.
State: Warm Up [WRM]
The Warm Up state is the same as the Pause state with the exception that the feedback from the system output signal enabled.
State: Update [UPD]
For the update state, the settings for the synchronization circuit 365 and the gain control module 355 and the delay blocks 240, 270 are frozen while the reset of the subsystem 200 functions normally.
State: Full Initialization [FIN]
The full initialization state initializes the update and feedback circuitry. In this state, the LUT entries are frozen in that the update blocks 220, 260 are disabled while the feedback path is active along with the additive adjustments of the α and β values. However, the distortion monitor control value and update duty control are both reset along with the update speed control (to be discussed below).
State: Re-Initialization [RIN]
For re-initialization, the feedback gain and phase adjustments (performed by gain control circuitry 355 and the synchronization circuitry 365) are adjusted while the delay settings (for the delay blocks 240, 270) are not modified. For this state, the α and β update blocks are disabled but predistortion continues. The distortion monitor control value is reset while the gain control block 355 and synchronization blocks 365 are operative.
State: Cool Down [CLD]
The cool down state maintains the feedback function after an LUT update to avoid problems which may be caused by the feedback circuit. It initiates a cool down timer, the duration of which is user configurable to enable avoidance of the feedback related problems. During this state, the distortion monitor is disabled (frozen) and the other relevant settings, such as the LUT entries, the setting for the gain control 355, the delay blocks 240, 270, and the synchronization circuitry 365, are also frozen. More importantly, the α and β LUTs are also disabled such that α′=α and β′=β.
The state transitions which the control block 600 undergoes between the different states is detailed in the table below. To assist in the understanding of the transition table below, FIG. 5A is presented. FIG. 5B is a state transition diagram illustrating the different states and the sequence of state transitions possible. For clarity, a listing of the timers implemented by the block 600 is also presented below.










State Transition Table











From
To
Condition







any
RST
Reset externally applied to





subsystem



RST
PSE
Completion of RST





operations



PSE
WRM
Transmit process enable and





update request




HLD
Transmit process enable





otherwise



WRM
SBY
Start of transmit packet





samples at predistortion





subsystem



SBY
FIN
Expiration of SBY timer and





1. First SBY since last RST





or





2. Full initialization





procedure incomplete since





last FIN or





3. Full initialization request




RIN
Expiration of SBY timer, and





SBY to FIN conditions do not





apply, and





1. Reinitialization procedure





incomplete since last RIN or





2. Reinitialization request




UPD
Expiration of SBY timer,





otherwise



UPD
CLD
End of transmit packet





samples at predistortion





subsystem



CLD
RST
Expiration of CLD timer and





Reset request




PSE
Expiration of CLD timer





otherwise



HLD
PSE
End of transmit packet





samples at predistortion





subsystem



FIN
CLD
End of transmit packet





samples at predistortion





subsystem




UPD
Otherwise



RIN
CLD
End of transmit sample





packets at predistortion





subsystem




UPD
Othewise




















Timers Implemented by Control Block











Configuration





Parameter
Description
Range







SBY-Timer
Timer duration
0, 2, 4, 8, 16, 32,





64, or 128 samples



CLD-Timer
Timer duration
0, 2, 4, 8, 16, 32,





64, or 128 samples



RTO-Timer
Recency time-out
0 . . . 1 s in steps of




timer
at least 1 ms and ‘inf’










To build in more flexibility and robustness in the predistortion subsystem 200, the update speed parameter μG and μF are, as noted above, configurable. Another lookup table, internal to the control block 600, provides 4 possible settings for each of the update speed parameters. These possible values and the time slot (expressed in terms of the number of symbols received) for which they would be used are listed in the table below.














Time slot (expressed in a




symbol count, during which


μF
μG
the LUTs are updated)







1/2
1/4
 0 . . . 15


1/4
1/8
16 . . . 47


1/8
1/8
48 . . . 79


 1/16
1/8
81 and beyond









As noted above, an update duty control (a control value) is also used by the control block 600. The update duty control is a value which determines how often an LUT update is to be performed. In one contemplated application of the invention, the data in the signal to be predistorted arrives in discrete packages or packets with each packet containing multiple OFDM (Orthogonal Frequency Division Multiplexing) symbols. The LUT update duty cycle is expressed in terms of a fraction of the number of packets during which updates are performed. For instance, at a duty cycle of 1 the system updates the LUT's during all transmit packets, at 1/9 the system updates during the first packet, and then freezes the LUT during the following 8 packets: transmission and pre-distortion then occurs by means of the Hold state. In this case the cycle is 9 packets, the duty is one packet, and the idle time is 8 packets. More precisely, the duty cycle is defined by the idle time, and after expiration of the idle time an update request is issued.
In this application of the invention, the duty cycle is defined by the idle time with a request for an LUT update being made after the idle time has expired. The duty cycle may be defined by other criteria such as the number of packets transmitted but, for this application, it was found that the idle time was a more convenient choice.
It should however be noted that the duty cycle is not operative at all times. It is only applied after a specified delay of a specified number of symbols has been received. Before the symbol count, an LUT update is requested for each packet. A reset of the update duty control value accomplishes a start of a new period starting with a duty and a start of the delay symbol count. The duty is therefore 100% for an initial number of symbols after a reset of the update duty control value and, subsequently, the duty is as configured.
The table below documents the different possible values for the delay symbol count configuration.















Configuration





Parameter
Description
Range









UDC-IDL
Idling time between
0, 21 . . . 26, and




update requests
‘inf’ packets





(8 possible values)



UDC-DLY
Delay in terms of a
0, 21 . . . 26, and




transmit symbol count
‘inf’ symbols




after which the
(8 possible values)




duty cycle applies










For clarity, a value of “inf” for the idle time denotes no duty (no LUT updates requested as ever) and a value of 0 for the idle time denotes LUT updates for each packet. A delay value of 0 denotes that the duty cycle is effective immediately.
Another control value used by the control block 600 is the initialization duty control. This initialization duty control value control how many initializations are performed for every update of the LUTs. The initialization duty cycle is expressed in terms of a fraction of the number of update packets at the start of which an initialization is performed. For instance, at a duty cycle of 1 the system performs an initialization at the start of all update transmit packets, at 1/9 the system initializes at the first update packet, and performs an update without initialization during the following 8 update packets. In this case the cycle is 9 update, the duty is one update, and the idle time is 8 updates.
The distortion monitor 500 may, depending on the implementation, be a simple component that issues an alarm if the distortion exceeds a certain preset distortion value. However, a more complex yet more useful and robust implementation is possible. Such an implementation may use multiple states and state transitions for the distortion monitor. These states can be as follows:

    
    
        Reset (resets the monitor)
        Measure (monitors the distortion metric and requests a predistortion subsystem recovery of the distortion metic exceeds the preprogrammed threshold)
        Freeze (holds the state of the distortion metric between packets or during initialization or re-initialization of the predistortion subsystem)
        No Monitor (disables the alarms and does not monitor the distortion metric. This allows the predistortion subsystem to converge and the distortion to settle)
        Monitor (compares the distortion metric with the programmed threshold and enables the triggers to be tripped if the threshold is exceeded)
    
    


Different alarm levels and states may also be implemented such that, depending on the alarm level and state, different actions may be taken to address the distortion (e.g. full initialization of the subsystem, a reset of the subsystem, a re-initialization etc.).
To execute and implement the states noted above and to execute the software, the control block 600 may be implemented as a general purpose central processing unit with attendant circuitry. Referring to FIG. 5B, a block diagram of such a control block 600 is illustrated. The CPU 610 is coupled to some random access memory (RAM) 620, a register file 630, and read only memory (ROM) 640. The CPU 610 performs the execution of the software and interfaces with the predistortion subsystem 200 while the registers 630 and the RAM 620 hold the temporary values (such as the various control values) required by the software in the CPU 610. The ROM 640 holds the software being executed by the CPU 610.
While the control block 600 above is illustrated as a general purpose data processing system, other implementations, such as a application specific integrated circuit (ASIC) or combinational logic circuit in which the logic of the software is hardwired in the hardware, are possible.
Regarding the amplifier subsystem 10, FIG. 4 illustrates one embodiment of the subsystem 10. In FIG. 4, the signal decomposer 20 of FIG. 1 comprises a phasor fragmentation engine 20A along with phase modulation units 60A, 60B. The fragmentation engine 20A receives the magnitude (M) and phase (φ) representing the predistorted signal. The phasor fragmentation engine 20A deconstructs a predetermined modulation waveform (the predistorted signal) into signal components which are of equal and constant magnitude. Further information regarding the phasor fragmentation engines may be found in the applicant's co-pending application U.S. application Ser. No. 10/205,743 entitled COMPUTATIONAL CIRCUITS AND METHODS FOR PROCESSING MODULATED SIGNALS HAVING NON-CONSTANT ENVELOPES, which is hereby incorporated by reference. In FIG. 4, these signal components are denoted by angles α and β. These components are each received by the predistortion circuit blocks 120A, 120B which, respectively, contain LUT blocks 200A, 200B. The predistortion circuit blocks 120A, 120B also receive the input signal 30 along with the decomposition information 45A, 45B from the decomposer 20. The output of these predistortion circuit blocks 120A, 120B are received by phase modulation and filtering blocks 60A, 60B which process the predistorted components to produce phase modulated and filtered versions of the components. The signal component 70A is an RF signal with predistorted phase α while signal component 70B is an RF signal with predistorted phase β. These components 70A, 70B are then amplified by amplifiers 90A, 90B. The amplified components are then recombined using combiner 100. Signal decomposition methods other than the phasor fragmentation referred to above may also be used by the signal decomposer 20.
Regarding the Chireix architecture amplifier subsystem 10, it has been found that, for higher amplification efficiencies, switch mode amplifiers are preferred for the amplifiers 90A, 90B. Such switch mode amplifiers, specifically Class D and Class F power amplifiers, provide low output impedances that allow higher amplification efficiencies. A co-pending application filed on Oct. 16, 2002 and having U.S. Ser. No. 10/272,725 entitled CHIREIX ARCHITECTURE USING LOW IMPEDANCE AMPLIFIERS provides further information on the desirable components and is hereby incorporated by reference. Such types of amplifiers are not required for the invention to function but they have been found to provide performance at a desirable level.
It should further be noted that while there are only two parallel amplifiers 90A, 90B in FIG. 1 and FIG. 4, multiple parallel amplifiers may be used as long as the decomposer 20 decomposes the input signal 30 into enough components so that each component is separately amplified and phase modulated in parallel with the other components and as long as each component is also predistorted in parallel by multiple predistortion circuit blocks.
It should also be noted that the predistortion subsystem 10 explained above does not linearize a power amplifier as is well-known in the field. Instead, the predistortion subsystem linearizes a whole power amplifier system—the output of the whole amplifier system is linearized and not simply the output of a single amplifier. Also, unlike the linearizing systems for power amplifiers that are currently known, the amplifier system discussed in this document compensates for distortions that mostly occur at mid signal amplitudes. Current single amplifier linearization systems linearize distortions that occur at large signal amplitudes.
It should further be noted that the invention may be applied to any signal processing system which decomposes a signal into components and recombines them. It has been found that signal combiners (block 100 in FIG. 1) invariably cause distortions. These combiners use addition to recombine the components and improper signal addition, such as when recombining sinusoidal components, has been found to be one cause of the distortions in the system output signal. In the above embodiment, the phasor fragmentation engine decomposes the incoming signal into vectors and the improper addition of these vectors by the combiner 100 lead to distortions in the output signal.
While the above embodiment amplifies the input signal, albeit separately for each component, this need not be the only signal processing accomplished after the input signal is decomposed. Referring to FIG. 7, such a generalized system 10A (which could be part of a larger signal transmission system) is illustrated. The signal decomposer receives an input signal 30 and decomposes it into multiple components. The predistortion subsystem 200 predistorts these components and compensates for distortions introduced in the system output signal 110 by the improper or imperfect recombining of the input signal's components. These components are produced by the signal decomposer 20 and are separately processed by signal component processor blocks 75A, 75B. The processing executed by the blocks 75A, 75B may take the form of amplification (as in the embodiment above), phase modulation, a combination of the two, or any other signal processing which may be desired. As an example, each of the signal components illustrated in FIG. 4 may be separately phase modulated in addition to being amplified by amplifiers 90A-90B. The phase modulation may be accomplished separately or be incorporated in the signal decomposer or, as contemplated for the implementation illustrated in FIG. 4, incorporated into the modulation and filtering blocks 60A, 60B.
As can be seen in FIG. 7, the signal processing subsystem 10A receives the input signal and decomposes it by way of the signal decomposer 20 into components. The predistortion subsystem then predistorts the components. These predistorted components are then separately processed by the signal component processor blocks 75A, 75B and are then recombined by the recombiner 100.
A feedback signal processing block 400 receives a portion of the system output signal 110 and processes this portion so it may be used by the adaptive predistortion subsystem 120. As an example, the feedback signal processing block 400 may contain the AO converter 330, the conversion unit 340, the filtering module 337, and the demodulation module 335 illustrated in FIG. 4. The same block 400 may also contain processing blocks 410, 420 for adjusting the magnitude and/or phase of the feedback signal.
One advantage using the above invention is that it allows less stringent tolerances to be used for the system components. Previously, components had to be substantially matched so that signal processing could produce acceptable results. By using the above invention, less than substantially matched components may be used together. Errors due to a mismatch may be measured and compensated for by the predistortion subsystem.
A person understanding this invention may now conceive of alternative structures and embodiments or variations of the above all of which are intended to fall within the scope of the invention as defined in the claims that follow.

## Claims

1. A method, comprising:
decomposing an input signal into at least two separate components, the at least two separate components being processed separately;
deliberately predistorting the at least two separate components, responsive to an entry in a lookup table, to produce at least two predistorted signals, the deliberate predistortion being configured to compensate for distortion in an output signal; and
combining the predistorted signals to produce the output signal.

2. The method of claim 1, further comprising amplifying the at least two predistorted signals.

3. The method of claim 2, further comprising non-linearly amplifying the at least two predistorted signals.

4. The method of claim 2, further comprising switch mode amplifying the at least two predistorted signals.

5. The method of claim 1, further comprising:
retrieving the entry from the lookup table according to a magnitude of the input signal; and
updating the entry in the lookup table according to the output signal.

6. The method of claim 1, further comprising introducing at least some distortion responsive to the combining of the predistorted signals.

7. The method of claim 1, further comprising processing each one of the separate components using a different amplifier.

8. A method, comprising:
decomposing an input signal into at least two component signals;
applying a deliberate predistortion to each of the at least two component signals to produce at least two predistorted signals according to an entry retrieved from a lookup table;
combining the at least two predistorted signals to produce an output signal; and
adjusting the deliberate predistortion based on at least one characteristic of the output signal.

9. The method of claim 8, further comprising:
separately processing each of the at least two predistorted signals; and
combining the at least two separately processed predistorted signals to produce the output signal.

10. The method of claim 9, further comprising amplifying the at least two predistorted signals.

11. The method of claim 9, further comprising non-linearly amplifying the at least two predistorted signals.

12. The method of claim 9, further comprising switch mode amplifying the at least two predistorted signals.

13. The method of claim 9, further comprising:
retrieving the entry from the lookup table according to a magnitude of the input signal; and
updating the entry in the lookup table according to the output signal.

14. The method of claim 9, further comprising introducing at least some distortion responsive to the combining of the predistorted signals.

15. The method of claim 9, further comprising processing each one of the component signals using a different amplifier.