Method and apparatus for modeling and estimating the characteristics of a power amplifier by retaining even-order terms in estimating characteristics

Information

  • Patent Grant
  • 6504425
  • Patent Number
    6,504,425
  • Date Filed
    Thursday, June 7, 2001
    23 years ago
  • Date Issued
    Tuesday, January 7, 2003
    21 years ago
Abstract
Disclosed is an apparatus and method for modeling and estimating the characteristics of a power amplifier. A predistortion module generates a predistorted signal in response to a predistortion function and an input signal. A power amplifier receives the predistorted signal and generates an output signal. A polynomial module generates coefficients of a complex polynomial of order p (p is an integer greater than one) in response to the predistorted signal and the output signal. In particular, the complex polynomial is implemented with both even and odd terms, thereby improving the ability to accurately model the power amplifier.
Description




FIELD OF THE INVENTION




This invention relates generally to signal amplification and, in particular, to determining signal amplifier characteristics for intentionally induced distortion techniques utilized prior to and in conjunction with signal amplification.




BACKGROUND OF THE INVENTION




In the field of radio communication systems, it is a well-known problem that the power amplifiers present in transmission equipment operate in a non-linear fashion when the power amplifiers are operated near their peak output. As a result, the power amplifier introduces significant signal distortion that can appear in various forms. For example, if more than one signal is input into the power amplifier or power amplification stage, its non-linear characteristics can cause an undesirable multiplicative interaction of the signals being amplified, and the power amplifier's output can contain intermodulation products. These intermodulation products cause interference and crosstalk over the power amplifier's operational frequency range.




In power amplifier design, there is a trade off between distortion performance and efficiency. Linear amplifiers that operate under “Class A” conditions create little distortion but are inefficient, whereas non-linear amplifiers operated under “Class C” conditions are reasonably efficient but introduce significant distortions. While both efficiency and distortion are important considerations in amplifier design, efficiency becomes increasingly important at high power levels. Because of their efficiency, non-linear amplifiers are largely preferred, leaving the user with the problem of distortion.




In order to employ non-linear power amplifiers, techniques have been used to improve linearity and thereby reduce the effects of interference and crosstalk. Linearity can be achieved by application of various linearization techniques that reduce the distortion caused by non-linear amplification. Conventional amplifier linearization techniques can be broadly categorized as feedback, feedforward, or predistortion.




The last mentioned technique, predistortion, intentionally distorts the signal before the power amplifier so that the non-linearity of the power amplifier can be compensated. According to this technique, linearization is achieved by distorting an input signal according to a predistortion function in a manner that is inverse to the amplifier characteristic function. The predistortion technique can be applied at radio frequency (RF), intermediate frequency (IF), or at baseband.




In the baseband domain, the input signal information is at a much lower frequency, allowing digital methods to be employed. The predistortion function is applied to the input signal with the resulting predistorted signal being upconverted to IF and then finally to the RF carrier frequency. It is also possible to apply adaptive predistortion techniques where feedback from the output of the amplifier is used to update and correct the predistortion function.




The form of the predistortion function is dependent upon the model used to characterize the output of the amplifier. Predistortion functions in the baseband domain are typically implemented as a table of gain and phase weighting values within a digital signal processor. A Cartesian feedback method employs a quadrature representation of the signal being amplified. The incoming quadrature signals I and Q are compared to the feedback quadrature signals. Thus, there are two sets of coefficients, one for each quadrature channel, that are being updated to model the predistortion characteristics. In this manner, gain and phase non-linearities within the amplifier can be compensated. Performance is dependent on the size of the look up table and the number of bits used to represent the signal. Better performance and more adaptivity is achieved with larger look up tables and more bits albeit at the expense of longer processing times.




Predistortion functions are also modeled as polynomials. Ideal amplifiers have linear characteristics; consequently, amplifiers with slight non-linearities can be modeled as polynomials of only a few terms, with only odd terms being employed. Even terms are discarded because their use with negative-valued inputs can interfere with linearity. While processing demands are eased by excluding and limiting the number of terms in the polynomial modeling, performance is sacrificed.




Accordingly, there is a need for a device to more quickly and efficiently determine the characteristics of a non-linear amplifier.




SUMMARY OF THE INVENTION




The present invention teaches an apparatus and method for modeling and estimating the characteristics of a power amplifier. A predistortion module generates a predistorted signal in response to a predistortion function and an input signal. A power amplifier receives the predistorted signal and generates an output signal. A polynomial module generates coefficients of a complex polynomial of order p (p is an integer greater than one) in response to the predistorted signal and the output signal. The coefficients characterize the power amplifier. The complex polynomial is implemented with both even and odd terms. Even order terms are typically ignored for the polynomial modeling of amplifiers because the non-linear distortion signals caused by power amplifiers are strongest at odd order harmonic frequencies. Additionally, it is commonly thought that the use of even terms with negative valued inputs can interfere with the linearity, since negativity is lost at even powers. However, the inclusion of even order terms improves the ability to accurately model the power amplifier and thereby improves predistortion performance.




In another exemplary embodiment of the present device, the polynomial module employs a minimum mean squared error criteria to determine said polynomial coefficients, thereby allowing a very fast and efficient implementation.




By improving the ability to model power amplifiers, the present invention improves the ability to model the power amplifier predistortion function. The invention further enables power amplifiers to be operated in the non-linear region near saturation, yet suppresses undesirable intermodulation products. Resort to a larger amplifier, to keep operation within the linear region, is avoided. Power amplifier sizes are kept small with associated cost savings, particularly important in the field of wireless communications.




The above factors make the present invention essential for effective power amplifier predistortion.











BRIEF DESCRIPTION OF THE DRAWINGS




For a better understanding of the present invention, reference may be had to the following description of exemplary embodiments thereof, considered in conjunction with the accompanying drawings, in which:





FIG. 1

is a block diagram providing an overview of an exemplary system employing adaptive power amplifier predistortion;





FIG. 2

is a block diagram of the simplified baseband model for power amplifier predistortion; and





FIG. 3

is a block diagram illustrating the device of

FIG. 2

as used in a RF transmission system in accordance with the principles of the present invention.











DETAILED DESCRIPTION




The following description is presented to enable a person skilled in the art to make and use the invention, and is provided in the context of a particular application and its requirements. Various modifications to the disclosed embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be applied to other embodiments and applications without departing from the spirit and the scope of the invention. Thus, the present invention is not intended to be limited to the embodiments disclosed, but is to be accorded the widest scope consistent with the principles and features disclosed herein.




The specification initially discusses the general concept and principles of adaptive digital predistortion in view of the novel system for determining the characteristics of a power amplifier. Exemplary embodiments of the system for determining the characteristics of a power amplifier are then described.




Overview of Adaptive Power Amplifier Predistortion




The principal benefit of the present invention is the ability to more efficiently model the power amplifier characteristics in order to improve the ability to employ adaptive digital predistortion (ADPD). The structure of an exemplary ADPD system is seen in FIG.


1


. An initial baseband digital signal


10


is identified as u


n


, where n is the time index. The initial baseband digital signal


10


is fed into a predistortion system


20


that is described as a function A(·). The output of the predistortion system


20


is the baseband digital input signal


12


to the power amplifier


50


and is defined as x


n


. The baseband digital input signal


12


is processed by a digital to analog (D/A) converter


30


with the resulting baseband analog signal being processed by an up-conversion means


40


that is comprised of mixers and filters and operates in the intermediate frequency (IF) range. The up-conversion means


40


outputs a signal in the radio frequency (RF) range and feeds the signal to the power amplifier


50


. While there are many methods for ADPD, the approach with the exemplary invention can be divided into two steps. First, the characteristics of the power amplifier


50


are estimated. Then, the predistortion function based on the power amplifier's


50


characteristics is obtained. For proper characterization of the power amplifier


50


, the time domain inputs and outputs of the power amplifier


50


need to be compared. Thus, the output of the power amplifier


50


is tapped and fed back to a down-conversion means


70


. Like the up-conversion means


40


, the down-conversion means


70


requires mixers and filters in the IF range. The output of the down-conversion means


70


is fed into an analog to digital (A/D) converter


80


. The output of the A/D converter


80


is input into a means for delay adjustment


82


with its output representing the baseband digital output signal


14


identified as y


n


. The baseband digital output signal


14


and the baseband digital input signal


12


are input to the polynomial module


25


in order to determine the coefficients that characterize the power amplifier


50


. The output of the polynomial module


25


is coupled to the predistortion polynomial module


15


that determines the predistortion coefficients that are fed into the predistortion module


20


. The polynomial module, predistortion polynomial module and predistortion module may be implemented in hardware, or in other forms such as software or firmware.




As implemented in

FIG. 1

, the baseband digital input signal


12


to the power amplifier


50


as well as the baseband digital output signal


14


of the power amplifier


50


are easily accessible. However, the up-conversion means


40


and the down-conversion means


70


distort the signals. Mixers are non-linear devices and will add non-linear distortions. Furthermore, the phase responses of analog filters are not linear, thereby causing different time delays for different frequency components. Generally, these distortions can be considered negligible or can be corrected by using linear filters, and considered to be part of the baseband model for the power amplifier.




By neglecting the effect of up-conversion and down-conversion process, the whole predistortion process can be considered in the baseband domain. In

FIG. 2

, an exemplary baseband model for predistortion processing is illustrated.




The power amplifier


150


is defined as a baseband function B(·) with complex inputs and complex outputs.




Polynomial Modeling of the Power Amplifier




Predistortion requires the information on the characteristics of the power amplifier


150


. Proper baseband modeling of the power amplifier


150


is described herein. Since the power amplifier


150


is operating in radio frequency (RF) domain, the baseband model of the power amplifier


150


must be considered in complex numbers. Letting x and y, as previously defined, be the input and output of the power amplifier


150


, the following relationships can be defined,












y
=

B


(


b


,
x

)








=



b
1


x

+


b
2



&LeftBracketingBar;
x
&RightBracketingBar;


x

+


b
3




&LeftBracketingBar;
x
&RightBracketingBar;

2


x

+

+


b
p




&LeftBracketingBar;
x
&RightBracketingBar;


p
-
1



x









(
1
)













where p is the order of the polynomial, b


k


=b


kr


+jb


ki


, kε{1, 2, . . . ,p} and {right arrow over (b)} is a size 2p vector of complex polynomial coefficients defined as [b


1r


, b


2r


, . . . ,b


pr


, b


1i


, b


2i


, . . . b


pi


]. In general, p=5 is sufficient to model the power amplifier


150


.




As seen in equation (1), the complex polynomial is implemented with both even and odd terms, thereby improving the ability to accurately model the power amplifier. Typically, even order terms are ignored for polynomial modeling because the non-linear distortion signals caused by power amplifier are located at odd order harmonic frequencies. Moreover, the use of even order terms is generally bypassed because it is thought that negative-valued inputs can interfere with linearity. However, inclusion of even order terms as seen in the present invention retains the negativity of those terms and allows for better modeling of the power amplifier, thereby improving the predistortion performance.




Estimation of Power Amplifier Characteristics




An estimation is performed in order to obtain an optimum {right arrow over (b)} that describes the characteristics of the power amplifier


150


. A minimum mean squared error criteria is employed based upon the complex input and output samples of the power amplifier


150


. x


n


is the input sample and y


n


is the corresponding output sample, where nε{1, 2, . . . , N}. The error function f({right arrow over (b)}) is defined as













f


(

b


)


=

E


[


&LeftBracketingBar;


y
n

-

B


(


b


,

x
n


)



&RightBracketingBar;

2

]








=


1
N






n
=
1

N




&LeftBracketingBar;


y
n

-

B


(


b


,

x
n


)



&RightBracketingBar;

2









=


1
N






n
=
1

N



(



&LeftBracketingBar;

y
n

&RightBracketingBar;

2

-


y
n
*



B


(


b


,

x
n


)



-


y
n

·


B


(


b


,

x
n


)


*


+


&LeftBracketingBar;

B


(


b


,

x
n


)


&RightBracketingBar;

2


)










(
2
)













where E[x] is the mean of x, and x* is the complex conjugate of x. Minimizing the error or merit function results in the most accurate modeling of the power amplifier


150


and thereby the optimal polynomial coefficients. From equations (1) and (2), it is known that f({right arrow over (b)}) is a quadratic function of {right arrow over (b)}. Thus, f({right arrow over (b)}) can be expressed by a Taylor series as








f


(


{right arrow over (b)}+{right arrow over (d)}


)=


f


(


{right arrow over (b)}


)+∇


f


(


{right arrow over (b)}





{right arrow over (d)}




t{right arrow over (b)}


+{fraction (


1


/


2


)}


{right arrow over (b)}·H


(


{right arrow over (b)}





{right arrow over (d


t


)}


  (3)






where “t” is the transpose of the matrix, and ∇f({right arrow over (b)}) is the gradient of f({right arrow over (b)}), defined as












f


(

b


)





[




f




b

1

r




,



f




b

2

r




,





,



f




b
pr



,



f




b

1

i




,



f




b

2

i




,





,



f




b
pi




]





(
4
)













H({right arrow over (b)}) is the Hessian or second order derivative of f({right arrow over (b)}), and is defined as










H


(

b


)




[







2


f




b

1

r

2









2


f





b

1

r







b

2

r














2


f





b

1

r







b
pr










2


f





b

1

r







b

1

i











2


f





b

1

r







b

2

i














2


f





b

1

r







b
pi












2


f




b

2

r

2









2


f





b

2

r







b

2

r














2


f





b

2

r







b
pr










2


f





b

2

r







b

1

i











2


f





b

2

r







b

2

i














2


f





b

2

r







b
pi






































2


f




b
pr
2









2


f





b
pr






b

2

r














2


f





b
pr






b
pr










2


f





b
pr






b

1

i











2


f





b
pr






b

2

i














2


f





b
pr






b
pi












2


f




b

1

i

2









2


f





b

1

i







b

2

r














2


f





b

1

i







b
pr










2


f





b

1

i







b

1

i











2


f





b

1

i







b

2

i














2


f





b

1

i







b
pi












2


f




b

2

i

2









2


f





b

2

i







b

2

r














2


f





b

2

i







b
pr










2


f





b

2

i







b

1

i











2


f





b

2

i







b

2

i














2


f





b

2

i







b
pi






































2


f




b
pi
2









2


f





b
pi






b

2

r














2


f





b
pi






b
pr










2


f





b
pi






b

1

i











2


f





b
pi






b

2

i














2


f





b
pi






b
pi







]





(
5
)













Since f({right arrow over (b)}) is a quadratic function, it has a global minimum if H({right arrow over (b)}) is a positive definite matrix. By the same reason, the Hessian is not the function of {right arrow over (b)} and the notation of H can be used instead of H({right arrow over (b)}). The estimation obtains the optimum {right arrow over (b


opt


)} that makes f({right arrow over (b


opt


)}) on its minimum. Using the Newton Method, the optimum {right arrow over (b


opt


)} can be obtained. Letting {right arrow over (b)}


0


be an initial value, then {right arrow over (b


opt


)} can be expressed as






ƒ


{right arrow over (b


opt


)}=




{right arrow over (b


0


)}+




{right arrow over (d)}


  (6)






Equation (3) can then be re-written as








f


(


{right arrow over (b


opt


)})=




f


(


{right arrow over (b


0


)}+




{right arrow over (d)}) =




f


(


{right arrow over (b


0


)})+∇




f


(


{right arrow over (b


0


)}) ·




{right arrow over (d


t


)}+


½


{right arrow over (d)}·H·{right arrow over (d





t


)}  (7)






If H is a positive definite matrix, the optimum {right arrow over (d)} is found by making the gradient of f({right arrow over (b


0


)}+{right arrow over (d)}) based on {right arrow over (d)} equal zero.








{right arrow over (d)}=−H




−1




·∇f


(


{right arrow over (b


0


)})


  (8)






Thus, the optimum polynomial coefficients {right arrow over (b


opt


)} is








{right arrow over (b


opt


)}=




{right arrow over (b


0


)}−




H




−1




·∇f


(


{right arrow over (b


0


)})


  (9)






The gradient of f({right arrow over (b


0


)}) is calculated from equations (1), (2), and (4).












f


(


b
o



)



=



(


1
N






n
=
1

N




&LeftBracketingBar;


y
n

-

B


(



b
o



,

x
n


)



&RightBracketingBar;

2



)






(
10
)













Where R[x] and I[x] are the real components of x and the imaginary component of x, respectively. From the results,









B


(


{right arrow over (b


01


)}




x




n


)=[


x




n




,|x




n




|x




n




, . . . ,|x




n


|


p−1




x




n




, jx




n




, j|x




n




|x




n




, . . . j|x




n


|


p−1




x




n


]  (11)






then,












f


(


b
o



)



=


[






b
o1r



M
2


+


b
o2r



M
3


+

+


b
opr



M

p
+
1



-




[

C
1

]










b
o1r



M
3


+


b
o2r



M
4


+

+


b
opr



M

p
+
2



-




[

C
2

]















b
o1r



M

p
+
1



+


b
o2r



M

p
+
2



+

+


b
opr



M

2

p



+




[

C
p

]










b
o1i



M
2


+


b
o2i



M
3


+

+


b
opi



M

p
+
1



+

𝔍


[

C
1

]










b
o1i



M
3


+


b
o2i



M
4


+

+


b
opi



M

p
+
2



+

𝔍


[

C
2

]















b
o1i



M

p
+
1



+


b
o2i



M

p
+
2



+

+


b
opi



M

2

p



+

𝔍


[

C
p

]






]

t





(
12
)













where










M
k




2
N






n
=
1

N




&LeftBracketingBar;

x
n

&RightBracketingBar;

k







(
13
)







C
k




2
N






n
=
1

N





&LeftBracketingBar;

x
n

&RightBracketingBar;


k
-
1


·

x
n

·

y
n
*








(
14
)













Thus, the Hessian of f({right arrow over (b


0


)})









H
=


2
N






n
=
1

N



[





&LeftBracketingBar;

x
n

&RightBracketingBar;

2





&LeftBracketingBar;

x
n

&RightBracketingBar;

3








&LeftBracketingBar;

x
n

&RightBracketingBar;


p
+
1




0


0





0






&LeftBracketingBar;

x
n

&RightBracketingBar;

3





&LeftBracketingBar;

x
n

&RightBracketingBar;

4








&LeftBracketingBar;

x
n

&RightBracketingBar;


p
+
2




0


0





0
































&LeftBracketingBar;

x
n

&RightBracketingBar;


p
+
1






&LeftBracketingBar;

x
n

&RightBracketingBar;


p
+
2









&LeftBracketingBar;

x
n

&RightBracketingBar;


2

p




0


0





0




0


0





0




&LeftBracketingBar;

x
n

&RightBracketingBar;

2





&LeftBracketingBar;

x
n

&RightBracketingBar;

3








&LeftBracketingBar;

x
n

&RightBracketingBar;


p
+
1






0


0





0




&LeftBracketingBar;

x
n

&RightBracketingBar;

3





&LeftBracketingBar;

x
n

&RightBracketingBar;

4








&LeftBracketingBar;

x
n

&RightBracketingBar;


p
+
2
































0


0





0




&LeftBracketingBar;

x
n

&RightBracketingBar;


p
+
1






&LeftBracketingBar;

x
n

&RightBracketingBar;


p
+
2









&LeftBracketingBar;

x
n

&RightBracketingBar;


2

p





]







(
15
)













Defining a p×p matrix K as









K


[




M
2




M
3







M

p
+
1







M
3




M
4







M

p
+
2





















M

p
+
1





M

p
+
2








M

2

p





]





(
16
)













The Hessian can now be expressed in terms of K.









H
=

[



K


0




0


K



]





(
17
)













Where 0 is a zero matrix of size p×p. Finally, the Inverse Hessian can be obtained as










H

-
1


=

[




K

-
1




0




0



K

-
1





]





(
18
)














FIG. 3

illustrates a base station


310


with power amplifiers employing a predistortion linearization technique in accordance with the principles of the present invention.




As shown in

FIG. 3

, base station


310


comprises a pair of transmitters


320


each having a power amplifier as shown in FIG.


2


. Base station


310


can comprise, if necessary, a single transmitter or additional transmitters. In addition, base station


310


includes a suitable transmit antenna


315


for transmission in a RF transmission system that comprises both wireless and wired equipment. Base station


310


can utilize any equipment suitable for sending and receiving RF transmissions, such as those employing Code Division Multiple Access (CDMA) communications. In

FIG. 3

, a mobile radio


350


is shown as well as the base station


310


including a receive antenna


325


, a pair of receivers


330


and a multiplexer


340


. Additional mobile radios may be serviced by the base station


310


, and it will be apparent to one of ordinary skill that base station


310


can be used for providing wireless communications in any desired manner and for any type of wireless communications protocol or standard.




Numerous modifications and alternative embodiments of the invention will be apparent to those skilled in the art in view of the foregoing description. Accordingly, this description is to be construed as illustrative only and is for the purpose of teaching those skilled in the art the best mode of carrying out the invention. Details of the structure may be varied substantially without departing from the spirit of the invention and the exclusive use of all modifications which come within the scope of the appended claim is reserved.



Claims
  • 1. An apparatus for modeling and estimating the characteristics of a power amplifier, comprising:a predistortion module responsive to a predistortion function and an input signal by generating a predistorted signal; the power amplifier responsive to said predistorted signal by generating an output signal; and a polynomial module responsive to said predistorted signal and said output signal by generating coefficients of a complex polynomial of order p (p is an integer greater than one), said complex polynomial having both even order terms and odd order terms, both of said even and odd order terms being used by said polynomial module to model said power amplifier.
  • 2. The apparatus according to claim 1, wherein said polynomial module employs a minimum mean squared error criteria to determine said polynomial coefficients.
  • 3. A wireless radio frequency communications system including an apparatus for modeling and estimating the characteristics of a power amplifier, said system comprising:a predistortion module responsive to a predistortion function and an input signal by generating a predistorted signal; the power amplifier responsive to said predistorted signal by generating an output signal; and a polynomial module responsive to said predistorted signal and said output signal by generating coefficients of a complex polynomial of order p (p is an integer greater than one), said complex polynomial having both even order terms and odd order terms, both of said even and odd order terms being used by said polynomial module to model said power amplifier.
  • 4. The system of claim 3, wherein said polynomial module employs a minimum mean squared error criteria to determine said polynomial coefficients.
  • 5. A method for modeling and estimating the characteristics of a power amplifier, comprising:generating a predistorted signal in response to a predistortion function and an input signal; amplifying said predistorted signal to generate an output signal; and generating coefficients of a complex polynomial of order p (p is an integer greater than one), in response to said predistorted signal and said output signal, said complex polynomial having both even order terms and odd order terms, both of said even and odd order terms being used to model said power amplifier.
  • 6. The method according to claim 5, wherein said coefficients generating step employs a minimum mean squared error criteria to determine said polynomial coefficients.
RELATED APPLICATIONS

The present patent application is related to U.S. patent application Ser. No. 09/876,772, entitled, “METHOD AND APPARATUS FOR MODELING AND ESTIMATING THE CHARACTERISTICS OF A POWER AMPLIFIER”, being concurrently filed herewith and having a filing date of Jun. 7, 2001; to U.S. patent application Ser. No. 09/878,037, entitled, “METHOD AND APPARATUS FOR CALCULATING THE PREDISTORTION FUNCTION FROM A POWER AMPLIFIER”, being concurrently filed herewith and having a filing date of Jun. 8, 2001; and to U.S. patent application Ser. No. 09/878,038, entitled, “METHOD AND APPARATUS FOR CALCULATING THE PREDISTORTION FUNCTION FROM A POWER AMPLIFIER”, being concurrently filed herewith and having a filing date of Jun. 8, 2001; all of which have a common inventor and assignee and being incorporated herein by reference.

US Referenced Citations (3)
Number Name Date Kind
5900778 Stonick et al. May 1999 A
6236837 Midya May 2001 B1
6313703 Wright et al. Nov 2001 B1