WO2014184775A1 - Baseband equivalent volterra series for behavioral modeling and digital predistortion of wideband transmitters - Google Patents
Baseband equivalent volterra series for behavioral modeling and digital predistortion of wideband transmitters Download PDFInfo
- Publication number
- WO2014184775A1 WO2014184775A1 PCT/IB2014/061477 IB2014061477W WO2014184775A1 WO 2014184775 A1 WO2014184775 A1 WO 2014184775A1 IB 2014061477 W IB2014061477 W IB 2014061477W WO 2014184775 A1 WO2014184775 A1 WO 2014184775A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- series
- continuous
- signal
- volterra
- kernels
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F1/00—Details of amplifiers with only discharge tubes, only semiconductor devices or only unspecified devices as amplifying elements
- H03F1/02—Modifications of amplifiers to raise the efficiency, e.g. gliding Class A stages, use of an auxiliary oscillation
- H03F1/0205—Modifications of amplifiers to raise the efficiency, e.g. gliding Class A stages, use of an auxiliary oscillation in transistor amplifiers
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F1/00—Details of amplifiers with only discharge tubes, only semiconductor devices or only unspecified devices as amplifying elements
- H03F1/32—Modifications of amplifiers to reduce non-linear distortion
- H03F1/3241—Modifications of amplifiers to reduce non-linear distortion using predistortion circuits
- H03F1/3247—Modifications of amplifiers to reduce non-linear distortion using predistortion circuits using feedback acting on predistortion circuits
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F3/00—Amplifiers with only discharge tubes or only semiconductor devices as amplifying elements
- H03F3/189—High frequency amplifiers, e.g. radio frequency amplifiers
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F3/00—Amplifiers with only discharge tubes or only semiconductor devices as amplifying elements
- H03F3/189—High frequency amplifiers, e.g. radio frequency amplifiers
- H03F3/19—High frequency amplifiers, e.g. radio frequency amplifiers with semiconductor devices only
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F3/00—Amplifiers with only discharge tubes or only semiconductor devices as amplifying elements
- H03F3/20—Power amplifiers, e.g. Class B amplifiers, Class C amplifiers
- H03F3/21—Power amplifiers, e.g. Class B amplifiers, Class C amplifiers with semiconductor devices only
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F3/00—Amplifiers with only discharge tubes or only semiconductor devices as amplifying elements
- H03F3/20—Power amplifiers, e.g. Class B amplifiers, Class C amplifiers
- H03F3/24—Power amplifiers, e.g. Class B amplifiers, Class C amplifiers of transmitter output stages
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/32—Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
- H04L27/34—Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
- H04L27/36—Modulator circuits; Transmitter circuits
- H04L27/366—Arrangements for compensating undesirable properties of the transmission path between the modulator and the demodulator
- H04L27/367—Arrangements for compensating undesirable properties of the transmission path between the modulator and the demodulator using predistortion
- H04L27/368—Arrangements for compensating undesirable properties of the transmission path between the modulator and the demodulator using predistortion adaptive predistortion
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F2200/00—Indexing scheme relating to amplifiers
- H03F2200/111—Indexing scheme relating to amplifiers the amplifier being a dual or triple band amplifier, e.g. 900 and 1800 MHz, e.g. switched or not switched, simultaneously or not
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F2200/00—Indexing scheme relating to amplifiers
- H03F2200/336—A I/Q, i.e. phase quadrature, modulator or demodulator being used in an amplifying circuit
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F2200/00—Indexing scheme relating to amplifiers
- H03F2200/451—Indexing scheme relating to amplifiers the amplifier being a radio frequency amplifier
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F2201/00—Indexing scheme relating to details of amplifiers with only discharge tubes, only semiconductor devices or only unspecified devices as amplifying elements covered by H03F1/00
- H03F2201/32—Indexing scheme relating to modifications of amplifiers to reduce non-linear distortion
- H03F2201/3206—Multiple channels are combined and amplified by only one amplifier
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F2201/00—Indexing scheme relating to details of amplifiers with only discharge tubes, only semiconductor devices or only unspecified devices as amplifying elements covered by H03F1/00
- H03F2201/32—Indexing scheme relating to modifications of amplifiers to reduce non-linear distortion
- H03F2201/3209—Indexing scheme relating to modifications of amplifiers to reduce non-linear distortion the amplifier comprising means for compensating memory effects
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03F—AMPLIFIERS
- H03F2201/00—Indexing scheme relating to details of amplifiers with only discharge tubes, only semiconductor devices or only unspecified devices as amplifying elements covered by H03F1/00
- H03F2201/32—Indexing scheme relating to modifications of amplifiers to reduce non-linear distortion
- H03F2201/3224—Predistortion being done for compensating memory effects
Definitions
- the present invention relates to wireless communication systems, and in particular to a method and system for use of Volterra series for behavioral modeling and digital pre-distortion of wideband transmitters.
- PAs PAs
- PAPR peak to average power ratio
- DDR dynamic deviation reduction
- the PA's memory effects gain intensity, and consequently a larger value for the dynamic parameter of the DDR technique is needed and eventually the complexity of the pruned LPE Volterra series undesirably approaches the un- pruned series.
- intra-band carrier aggregated signals proportional to the bandwidth of the individual component carriers and not to the frequency separation.
- the frequency separation is significantly lower.
- the spectrum regrowth (about five times the bandwidth of component carrier bandwidth) engendered at the different component carriers can overlap.
- This spectrum overlap hinders the applicability of multi-band models since these models require the acquisition of the individual envelopes of the component carriers at the output of the PA.
- intra-band carrier aggregated systems need to be handled as single-input single-output systems, and call for accurate and low- complex behavioral modelling schemes.
- the first approach approximates the behavior of the PA using its pass band input and output signals.
- the inherent complexity associated with this approach has limited the pass band model's adoption in the area of RF PA modeling and linearization.
- the LPE modeling approach is used. It is applied to the complex envelope of the RF input/output signals and requires much simpler measurement hardware and reduced computation as compared to its pass band counterpart.
- the LPE modeling approach capitalizes on the band limited characteristics of communication signals, thus limiting approximation efforts to the PA distortions which affect the signal envelope around the carrier frequency.
- the PA is treated as an envelope processing system.
- the LPE methodology as shown in FIG. 1, consists of modeling the equivalent envelope behavior of the PA via a low pass transformation, as shown by the PA 4. This is illustrated in equation (2)
- x(n) and y(n) denote the envelopes of the input and output signals, respectively, around the carrier and / designates the LPE model.
- x(t) and y (t) designate the input and output RF pass band signals respectively
- hp denotes the Volterra series' kernels.
- the LPE Volterra model of equation (4) has been used extensively in RF research.
- the model has been applied to develop solutions related to nonlinear communication system modeling and estimation, satellite communication, digital transmission channel equalization, multichannel nonlinear CDMA system equalization, analysis and cancellation of the inter-carrier interference in nonlinear OFDM systems, decision feedback equalization, nonlinear system and circuit analysis, data predistortion, PA modeling, and DPD.
- a method includes receiving the one of the wideband signal and the intra-band carrier aggregated signal and generating a discrete base band equivalent, BBE, Volte rra series based on the received signal, where the series has distortion products grouped according to determined shared kernels.
- the shared kernels are determined based on a transformation of a real-valued continuous-time pass band Volterra series without pruning of kernels.
- the shared kernels are determined based on the transformation of the real-valued continuous-time pass band Volterra series by steps that include transforming the real-valued continuous time pass band Volterra series to a multi-frequency complex-valued envelope series.
- the multi-frequency complex-valued envelope signal is transformed to a continuous-time pass band-only series, which is then transformed to a continuous-time baseband equivalent series.
- the continuous- time baseband equivalent signal is discretized to produce the discrete base band equivalent Volterra series.
- Shared kernels are identified, each shared kernel having distortion products in common with another shared kernel.
- transforming the continuous- time pass band-only signal to a continuous- time baseband equivalent signal includes expressing the continuous- time pass band-only series in convolution form.
- the Laplace transform is applied to the convolution form to produce a Laplace domain expression, which is frequency shifted to baseband to produce a baseband equivalent expression in the Laplace domain.
- the inverse Laplace transform is applied to the baseband equivalent expression to produce the continuous-time baseband equivalent series.
- a number of terms in the Laplace domain expression are reduced via symmetry.
- terms of the Laplace domain expression are grouped based on frequency intervals where distortion terms are not zero.
- discretizing the continuous- time baseband equivalent series to produce the discrete base band equivalent Volterra series includes truncating the continuous-time baseband equivalent series to a finite non-linearity order, and expressing the truncated series as summations of non-linear distortion terms, with upper limits of the summations being memory depths, different memory depths being assignable to different ones of the summations.
- a distortion term is a group of distortion products multiplied by a shared kernel.
- embodiments include a digital pre-distorter (DPD) system for pre-distortion of one of a wideband signal and an intra-band carrier aggregated signal.
- the system includes a Volterra series DPD modelling unit.
- the DPD modelling unit is configured to receive the one of the wideband signal and the intra-band carrier aggregated signal.
- the DPD modelling unit is also configured to calculate a discrete base band equivalent, BBE, Volterra series based on the received signal.
- the series has distortion products grouped according to determined shared kernels.
- the shared kernels are determined based on a transformation of a real-valued continuous-time pass band Volterra series without pruning of kernels.
- the DPD may further comprise a power amplifier configured to produce an output in response to the one of the wideband signal and the intra-band carrier aggregated received signal.
- the output of the power amplifier is provided to the Volterra series DPD modelling unit
- the Volterra series DPD modeling unit is configured to compute the shared kernels based on the output of the power amplifier.
- the DPD system further comprises a transmitter observation receiver configured to sample the output of the power amplifier and provide the sampled output to the Volterra series DPD modelling unit.
- the distortion products and their associated kernels are determined by transforming the real-valued continuous time pass band Volterra series to a discrete base band equivalent Volterra series according to a series of steps that include: transforming the real- valued continuous time pass band Volterra series to a multi-frequency complex- valued envelope series; transforming the multi-frequency complex-valued envelope signal to a continuous-time pass band-only series; transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent series; discretizing the continuous-time baseband equivalent signal to produce the discrete base band equivalent Volterra series.
- the shared kernels are identified using the envelopes of the input and output signals of the PA.
- transforming the continuous- time pass band-only signal to a continuous-time baseband equivalent signal includes the following steps:
- transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent signal further includes grouping terms of the Laplace domain expression based on frequency intervals where distortion terms are not zero.
- discretizing the continuous-time baseband equivalent series to produce the discrete base band equivalent Volterra series includes: truncating the continuous- time baseband equivalent series to a finite non-linearity order; and expressing the truncated series as summations of non-linear distortion terms, with upper limits of the summations being memory depths assigned to each order of the non-linear distortion terms.
- a distortion term is a group of distortion products multiplied by a shared kernel.
- embodiments include a Volterra series digital pre-distorter, DPD, modelling unit.
- the DPD modelling unit includes a memory module configured to store terms of a discrete base band equivalent, BBE, Volte rra series. The series is based on one of a received wideband signal and an intra-band carrier aggregated signal. Also, a grouping module is configured to group distortion products of the series according to determined shared kernels.
- the DPD modelling unit also includes a shared kernel determiner configured to determine the shared kernels based on a
- a series term calculator is configured to calculate the terms of the discrete base band equivalent Volterra series, the terms being the distortion products multiplied by their respective shared kernels.
- the shared kernel determiner is further configured to determine the shared kernels via a least squares estimate based on the one of the wideband signal and the intra-band carrier aggregated received signal, and an output of a power amplifier.
- the kernels and distortion products are derived from the real-valued continuous-time pass band Volterra series by: transforming the real-valued continuous time pass band Volterra series to a multi-frequency complex-valued envelope series; transforming the multi-frequency complex-valued envelope signal to a continuous-time pass band-only series; transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent series; discretizing the continuous-time baseband equivalent signal to produce the discrete base band equivalent Volterra series; and identifying the shared kernels, each shared kernel having distortion products in common.
- transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent signal includes: expressing the continuous-time pass band-only series in convolution form; applying a Laplace transform to the convolution form to produce a Laplace domain expression; frequency shifting the Laplace domain expression to baseband to produce a baseband equivalent expression in the Laplace domain; and applying an inverse Laplace transform to the baseband equivalent expression to produce the continuous-time baseband equivalent series.
- FIG. 1 is a diagram of low pass transformation for a low pass equivalent (LPE) power amplifier
- FIG. 2 is a column chart showing numbers of kernels for different exemplary Volte rra series implementations
- FIG. 3 is a spectrum of a PA output for several different linearization methods for a wideband input signal
- FIG. 4 is a spectrum a PA output for several different linearization methods for an intra-band carrier aggregated input signal, for a first PA;
- FIG. 5 is a spectrum a PA output for several different linearization methods for an intra-band carrier aggregated input signal, for a second PA;
- FIG. 6 is a block diagram of an exemplary power amplification system having a digital pre-distorter modelling unit implementing the BBE Volterra model presented herein;
- FIG. 7 is a block diagram of the pre-distorter modelling unit of FIG. 6;
- FIG. 8 is a flowchart of an exemplary process for modelling a power amplifier fed by a wideband or intra-band carrier aggregated input signal;
- FIG. 9 is a flowchart of an exemplary process for transforming a real- valued continuous time pass band Volterra series to a discrete base band equivalent Volterra series in which shared kernels are identified;
- FIG. 10 is a flowchart of an exemplary process of transforming the continuous- time pass band-only signal to a continuous-time baseband equivalent signal
- FIG. 11 is a flowchart of an exemplary process of discretizing the continuous- time baseband equivalent series to produce the discrete base band equivalent Volterra series.
- relational terms such as “first” and “second,” “top” and “bottom,” and the like, may be used solely to distinguish one entity or element from another entity or element without necessarily requiring or implying any physical or logical relationship or order between such entities or elements.
- a base band equivalent (BBE) expression of the Volterra series for linearization of a power amplifier driven by one of a wideband signal and an intra-band carrier aggregated signal is provided and the procedure to derive this new expression is described below in detail.
- a wideband signal and an intra- band carrier aggregated signal are distinguishable from a multi-band signal.
- a multi-band signal has component carriers in different communication bands. Communication bands may be defined by the third generation partnership project (3GPP) standards.
- 3GPP third generation partnership project
- a wideband signal has component carriers that are contiguous in the same communication band.
- An intra-band carrier aggregated signal has carriers that may be contiguous or non-contiguous, but are in the same communication band.
- Step 1 Continuous-time real-valued Volterra series modeling.
- Step 2 Real-valued to complex-valued envelope signal transformation.
- x(t) represents the complex baseband envelope signal that modulates the angular frequency ⁇ .
- y 0 (t) denotes the envelope at DC
- y 1 (t) denotes the envelopes of the first harmonic in-band signals
- y 2 (0 denotes the envelopes of the second order harmonic signals
- y WL (t) represents the NL th harmonic signal.
- Step 3 Multi-frequency to passband only transformation.
- a multi-harmonic model is derived which consists of several distinct equations relating the output envelopes y p (t) to (t) and ⁇ 3 ⁇ 4.. Since interest is in the relationship between the envelopes of the output and input signals around the carrier frequency, only the pass band component of the PA output is considered in the equation below:
- Step 4 Continuous time passband Volterra series.
- the following derivation is presented for the term around e J£Jct only; however, the same equation can be used to derive the conjugate term.
- the PA passband signal y ⁇ t) can be modeled as a summation of the Volterra series nonlinear order responses yi,2fc + i (t), or as the nonlinear order envelope responses y (t). around the angular frequency ⁇ ⁇ as follows:
- Equating the terms on the right sides of the expanded equations (8) and (10) yield a continuous BBE Volterra series that expresses ⁇ 2fe + i (i) as a function of x(t) and ⁇ u c .
- yi,i( J ⁇ ( ⁇ ) ⁇ x(.t - ⁇ ⁇ - ⁇ ,) .
- dTi (12) i,3( J J J l 3 (Ti,T 2 ,T 3 ) ⁇ ⁇
- a passband model In order to be computed by a digital processor, with manageable complexity, a passband model should be transformed to a BBE model.
- a BBE model includes the RF nonlinear dynamic distortion and computes the terms in baseband with a low sampling rate.
- To translate the passband model to baseband the continuous time passband Volterra series of equations (10) and (12) is first rewritten using convolution form.
- the convolution form for y 1 ⁇ 1 (t) is given by:
- the output function is re-assigned as follows:
- Step 6 Discrete-time baseband equivalent Volterra series model. Using the signal and system causality, the integral bounds of the model of equation (19) are adjusted to ⁇ t and - ⁇ o, respectively. The model is then given by:
- PA (t > T ⁇ ) (the transient-response, time-invariant Volterra series is defined as t ⁇ T ⁇ ), the corresponding model can be represented as:
- the model can be simplified to: ⁇ ( ( -J) ⁇ x ⁇ t— ⁇ - j ) ⁇ dr 1
- x 3 (n, i lt i 2 , j 3 ) x(n - i ⁇ xfji - i 2 )x*(n - i 3 )
- x 5 (n, i lt i 2 , i , i 4 , is) are linear combinations of 3 third, and 10 fifth, order distortion products.
- x(n)x(n— l) * (n— 1) and x(n— l)x(n— l) * (n) share the same kernel 3 ⁇ 4o,i,i m the BBE Volterra model of equation (26), however, they have different kernels (fio,i,i an d3 ⁇ 4i,i,o ⁇ respectively) in the LPE Volterra model of equation (5). More generally, several distortion products, which were assigned distinct kernels in equation (4), are linearly combined to form a novel term in the new model space proposed by equation (25), and hence, share the same kernel.
- the expressions for the number of coefficients were used to plot FIG. 2.
- the newly proposed, compact model involved a significantly lower number of coefficients than the full classical LPE Volterra series, especially as the memory depth increased. This results in lower complexity and cost of implementation.
- the pruned LPE Volterra using the DDR approach also required a smaller number of coefficients as it discarded distortion products with a dynamic order > three.
- Equation (25) The signal transformations used to derive the BBE expression of equation (25) preserve the linear property of the Volterra series with respect to its coefficients. Therefore, commonly used model identification algorithms such as the least square error (LSE) can be applied to identify the kernels in equation (25) for a given RF PA. Equation (27) suggests the expression used to deduce the LSE solution of equation (25)
- LSE least square error
- X, h and Y denote the input signal non-square matrix, the vector that contains the unknown kernels and the output signal vector as described below:
- N the total number of measurement data.
- a discrete BBE Volterra series is generated based on a received wideband or intra-band carrier aggregated signal.
- the series has distortion products grouped according to determined shared kernels.
- the shared kernels are determined based on a transformation of a real-valued continuous-time pass band Volterra series without pruning of kernels.
- the transformation includes transforming the real-valued continuous time pass band Volterra series to a multi-frequency complex-valued envelope series.
- the multi- frequency complex-valued envelope signal is then transformed to a continuous- time pass band-only series.
- the continuous-time pass band-only signal is transformed to a continuous-time baseband equivalent series.
- the continuous- time baseband equivalent signal is discretized to produce the discrete base band equivalent Volterra series.
- a 200 W Doherty PA (DPA) driven with a 2C 1001 WCDMA signal was linearized.
- the device under test (DUT) was excited with a 1 ms signal sampled at 100 MHz, where, out of 100k points, 10k were used to train the proposed model.
- the nonlinearity order and memory depth for both the proposed model and the LPE Volterra model were set to be equal.
- the DDR pruned Volterra model parameters were varied separately for optimum performance.
- an intra-band carrier aggregated signal was synthesized.
- This signal was mixed-standard, with a bandwidth of 40 MHz and a PAPR of 9.8 dB, containing a 15 MHz Long Term Evolution (LTE) signal and a 15 MHz 101 Wideband Code Division Multiple Access (WCDMA) signal separated by a 10 MHz guard band.
- LTE Long Term Evolution
- WCDMA Wideband Code Division Multiple Access
- Two different RF PA demonstrators were used as DUT. The first was a single ended Class J 45 W GaN amplifier and the second was a laterally diffused MOSFET (LDMOS) 200 W Doherty PA. The DUT was excited with a 1 ms signal sampled at 122.88 MHz, where, out of 122.88k points, 10k were used to train the proposed model.
- LDMOS laterally diffused MOSFET
- Tables III and IV summarize the linearization results found with each amplifier using both the proposed compact and the pruned LPE Volterra models. For all testing, the nonlinearity order was set to seven and the memory depth was set two and three.
- both the pruned LPE Volterra and the compact models allowed for similar linearization performance as the EVM was reduced from 5.5% to 2% and the ACLR was increased from 35 dB to about 50 dB (after application of the linearization).
- This linearization capacity was successfully achieved despite the significant reduction in the number of kernels required by the compact model (91 in the pruned LPE Volterra model versus 70 in the proposed model).
- a noticeable residual out of band spectrum regrowth in the signal obtained using the pruned DDR LPE can be observed.
- the DDR pruning approach may suffer from reduced linearization capability when the dynamic order is set to two as it eliminates a large number of distortion products.
- the compact BBE Volterra model described herein does not suffer from this limitation. All distortion products are maintained.
- FIG. 6 is a block diagram of a power amplification system 18 having a digital pre-distorter modelling unit 3 implementing the BBE Volterra model presented herein.
- the power amplification system 18 includes a digital pre- distorter 20.
- the DPD 20 receives the identified shared kernels from the pre- distorter modelling unit 34, and pre-distorts the input signal xto produce a pre- distorted signal.
- the pre-distorted signal is converted to analog by a digital to analog (DAC) converter 22, low pass filtered by a filter 24, and mixed to radio frequency (RF) by a mixer 26 to prepare the signal for amplification by an RF PA amplifier 28.
- DAC digital to analog
- RF radio frequency
- a directional coupler 30 couples the output of the RF PA 28 to a transmitter observation receiver (TOR) 32
- the TOR 32 samples the output of the power amplifier 28 in each band and produces a TOR output signal, y.
- the TOR output signal is used by the DPD modeling unit 34 to derive the kernel vector h according to equation (27).
- the DPD modelling unit 34 calculates a discrete baseband equivalent Volterra series having distortion products grouped according to determined shared kernels, where the shared kernels are based on a transformation of a real-valued continuous-time pass band Volterra series without pruning of kernels.
- FIG. 7 is a detailed view of the DPD modelling unit 34, which includes a processor 36 in communication with a memory module 38.
- the processor 36 may be implemented as a micro-processor operating according to computer instructions organized as a group of software modules, or each module of the processor 36 may be implemented by application specific integrated circuitry.
- the memory 38 may be implemented as random access memory (RAM), for example, and may contain non-volatile components, such as read only memory (ROM) that stores programmatic code to implement the functions described herein.
- the DPD modelling unit 34 receives the input signal xand the TOR output signal 'from the transmitter observation receiver 32 and derives the modelling vector h according to equation (27).
- the modelling vector h is input to the DPD 20 to pre-distortthe input signal xto produce a pre-distorted signal.
- the processor 36 includes a grouping module 40, a shared kernel determiner 42 and a series term computer 44.
- the grouping module 40 is configured to group distortion products of the series according to determined shared kernels, as shown in equation (25).
- the shared kernel determiner 42 is configured to determine the shared kernels based on a transformation of a real-valued continuous- time pass band Volterra series without pruning of kernels, as described above with reference to Step 6.
- the series term calculator 44 is configured to calculate the terms of the discrete base band equivalent Volterra series according to equation (27), the terms being the distortion products multiplied by their respective shared kernels.
- the memory module 38 is configured to store terms of the discrete BBE Volterra series 46, generated by the processor 36.
- FIG. 8 is a flowchart of an exemplary process for modelling a power amplifier 40 fed by a wideband or intra-band carrier aggregated input signal using the discrete BBE Volterra series described herein for improved
- the input signal is received by a digital pre-distorter 28 (block SlOO).
- a discrete BBE Volterra series is generated by the DPD modelling unit 34 based on the received input signal (block S102), according to equations (25) and (27).
- the series has distortion products that are grouped by the grouping module 40 according to determined shared kernels determined by shared kernel determiner 42, as shown in equation (25).
- the shared kernels are determined based on a transformation of a real-valued continuous- time pass band Volterra series without pruning of kernels, according to Steps 1 through 6 above.
- FIG. 9 is a flowchart of an exemplary process for transforming a real- valued continuous time pass band Volterra series to a discrete base band equivalent Volterra series according to Steps 1 through 6 above in which shared kernels are identified as set out in block S102 of FIG. 8.
- the real-valued continuous time pass band Volterra series is transformed to a multi-frequency complex-valued envelope series, as in Steps 1 and 2, described above (block S104).
- the multi-frequency complex-valued envelope signal is transformed to a continuous- time pass band-only series, as in Step 3 described above (block S106).
- the continuous-time pass band-only signal is transformed to a continuous- time baseband equivalent series, as in Steps 4 and 5 described above (block S108).
- the continuous-time baseband equivalent signal is discretized to produce the discrete base band equivalent Volterra series, as in Step 6 described above, (block S110).
- Shared kernels of the discrete base band equivalent Volterra series are identified, where a shared kernel has distortion products in common with another shared kernel, also described as in Step 6 described above (block SI 12).
- FIG. 10 is a flowchart of an exemplary process of transforming the continuous- time pass band-only signal to a continuous-time baseband equivalent signal as shown in block SI 08 of FIG. 9, as in Step 5 described above.
- the continuous-time pass band-only series is expressed in convolution form (block S114).
- the Laplace transform is applied to the convolution form to produce a Laplace domain expression (block S116).
- a number of terms in the Laplace domain expression may be reduced based on symmetry (block S118).
- the Laplace domain expression is frequency- shifted to baseband to produce a baseband equivalent expression in the Laplace domain (block S120).
- An inverse Laplace transform is applied to the baseband equivalent expression to produce the continuous- time baseband equivalent series (block S122).
- FIG. 11 is a flowchart of a process of discretizing the continuous-time baseband equivalent series to produce the discrete base band equivalent Volterra series, as shown in block SI 10 of FIG. 9, as in Step 6 described above.
- the process includes truncating the continuous-time baseband equivalent series to a finite non-linearity order (block S124).
- the process also includes expressing the truncated series as summations of non-linear distortion terms, with upper limits of the summations being memory depths assigned to each order of the non-linear distortion terms (block S126).
- the PA modelling method described herein using the discrete BBE Volterra series is inherently compact and calls for significantly fewer coefficients then its LPE counterpart and therefore avoids the accuracy-compromising pruning transformations widely applied to classical formulations in the literature.
- the validation described herein proved the excellent modeling and linearization performance of the method when compared to the classical non-pruned Volterra model.
- the new formulation also outperformed pruned Volterra models while using a lower number of coefficients.
Abstract
A method, system and apparatus for modelling a power amplifier and pre- distorter fed by one of a wideband signal and an intra-band carrier aggregated signal are disclosed. According to one aspect, a method includes receiving the one of the wideband signal and the intra-band carrier aggregated signal and generating a discrete base band equivalent, BBE, Volterra series based on the received signal, where the series has distortion products grouped according to determined shared kernels. The shared kernels are determined based on a transformation of a real-valued continuous- time pass band Volterra series without pruning of kernels.
Description
BASEBAND EQUIVALENT VOLTERRA SERIES FOR BEHAVIORAL MODELING AND DIGITAL PREDISTORTION OF WIDEBAND
TRANSMITTERS TECHNICAL FIELD
The present invention relates to wireless communication systems, and in particular to a method and system for use of Volterra series for behavioral modeling and digital pre-distortion of wideband transmitters.
BACKGROUND
Development of energy efficient radio frequency (RF) power amplifiers
(PAs) continues to be an area of focus for both industrial and academic researchers. This interest has been motivated in part by the ever increasing high peak to average power ratio (PAPR) and bandwidth of recent communication signal standards. Innovative PA architectures for use with extended back-off and bandwidth, such as Doherty PAs and envelope tracking, have already been devised and successfully demonstrated. However, the efficiency enhancement achieved by these architectures has been compromised by the accentuation of signal distortion problems. The search for a solution to address the distortion exhibited by these RF PAs has been a prolific area of research which has led to the development of several linearization techniques; namely feed-forward, feedback and predistortion techniques.
While the predistortion technique, particularly its baseband digital form, has been widely adopted to mitigate the distortions generated by RF PAs in 3rd generation (3G) wireless base stations, current techniques do not meet the challenges of 4th generation (4G) wireless networks using wideband and multi- standard signals. Investigators of digital predistortion (DPD) have struggled with linearization capacity and implementation complexity. Over time, DPD schemes have evolved from simple memoryless models to more sophisticated types which attempt to mitigate the memory effects that gain intensity as the signal bandwidth widens.
Several DPD schemes with memory have been proposed, mainly consisting of derivations of the low pass equivalent (LPE) Volterra series. In order to address the complexity burden of the LPE Volterra series, these derivations discard kernels that are considered negligible. One popular example is the memory polynomial model, where several polynomials are applied to the delayed input signals' samples separately and delayed cross-terms are discarded. Another approach, dynamic deviation reduction (DDR), has been employed to prune the LPE Volterra series. The approach begins by reformulating an LPE Volterra series expression to assemble terms with the same dynamic distortion order. By setting a value for the maximum allowed dynamic order, such a reformulation allows a simple elimination of distortion terms of higher dynamic order in the LPE Volterra series, consequently reducing its complexity. A DDR based LPE Volterra series with a preset value for the dynamic order equal to one or two has been used to linearize various RF PAs and significant linearization results were reported for different types of signals. Yet, as the signals' bandwidth increases, e.g., inter-band and inter-band carrier aggregated signals, the PA's memory effects gain intensity, and consequently a larger value for the dynamic parameter of the DDR technique is needed and eventually the complexity of the pruned LPE Volterra series undesirably approaches the un- pruned series.
Recently, several multi-band behavioral models suitable for inter-band carrier aggregated scenarios have been reported. The models are suited for carrier aggregation scenarios when the frequency separation between the two component carriers is significantly larger than the bandwidth of the individual component carriers. Consequently, these models operate at a speed
proportional to the bandwidth of the individual component carriers and not to the frequency separation. However, in the case of intra-band carrier aggregated signals, the frequency separation is significantly lower. Hence, the spectrum regrowth (about five times the bandwidth of component carrier bandwidth) engendered at the different component carriers can overlap. This spectrum overlap hinders the applicability of multi-band models since these models
require the acquisition of the individual envelopes of the component carriers at the output of the PA. Hence, intra-band carrier aggregated systems need to be handled as single-input single-output systems, and call for accurate and low- complex behavioral modelling schemes.
Some have distinguished two RF PA modeling strategies. The first approach approximates the behavior of the PA using its pass band input and output signals. The inherent complexity associated with this approach has limited the pass band model's adoption in the area of RF PA modeling and linearization. Instead, the LPE modeling approach is used. It is applied to the complex envelope of the RF input/output signals and requires much simpler measurement hardware and reduced computation as compared to its pass band counterpart. In fact, the LPE modeling approach capitalizes on the band limited characteristics of communication signals, thus limiting approximation efforts to the PA distortions which affect the signal envelope around the carrier frequency.
Essentially, the PA is treated as an envelope processing system.
Therefore, if the real RF PA behavior is expressed as equation (1): y(t) = f(x(f)), (1)
where x(t) is the real RF input signal, as shown by PA 2, y(t) is the real RF output signal and / is a describing function (modeling the RF PA behavior), then the LPE methodology, as shown in FIG. 1, consists of modeling the equivalent envelope behavior of the PA via a low pass transformation, as shown by the PA 4. This is illustrated in equation (2)
where x(n) and y(n) denote the envelopes of the input and output signals, respectively, around the carrier and / designates the LPE model.
Since a PA can be treated as a nonlinear dynamic system with fading memory, the Volte rra series outlined in equation (3) is a suitable modeling framework to approximate its behavior and/or synthesize the corresponding predistortion module.
NL V
y( = /_, I ··· I Λρ (τ!, . . , τρ) \ x(t - Tj) dTj . (3)
p = l 7 =1
In this model, x(t) and y (t) designate the input and output RF pass band signals respectively, and hp denotes the Volterra series' kernels. The direct application of the LPE strategy to the discrete input and output signals' envelopes yields the expression of equation (4) which is commonly used for RF PA behavioral modeling.
(p + D/2 p
[ x(n - ij - [ x*(n - ij) . (4)
7 = 1 (p+3)
' 2
In equation (4) (n) and ;y(«) designate the input and output signals' envelopes around the carrier sampled at fs =1/ Ts , with n, NL, Mand h
representing the sample index the nonlinearity order, the memory depth and the LPE complex Volterra kernels respectively.
The LPE Volterra model of equation (4) has been used extensively in RF research. The model has been applied to develop solutions related to nonlinear communication system modeling and estimation, satellite communication, digital transmission channel equalization, multichannel nonlinear CDMA system equalization, analysis and cancellation of the inter-carrier interference in nonlinear OFDM systems, decision feedback equalization, nonlinear system and circuit analysis, data predistortion, PA modeling, and DPD. Although
computationally more efficient than its pass band counterpart, the LPE Volterra series in its classical form (4) still suffers from a large number of kernels. As stated earlier, this impediment has been the key limitation to its widespread adoption for RF PA behavioral modeling and predistortion. Various attempts have been made to reduce the number of kernels required (e.g., Hammerstein,
Weiner, Parallel Hammerstein, Parallel Weiner, DDR based Pruned Volte rra series) but at the cost of reduced modeling accuracy.
SUMMARY
Methods, systems and apparatus for modelling a power amplifier and pre- distorter fed by one of a wideband signal and an intra-band carrier aggregated signal are disclosed. According to one aspect, a method includes receiving the one of the wideband signal and the intra-band carrier aggregated signal and generating a discrete base band equivalent, BBE, Volte rra series based on the received signal, where the series has distortion products grouped according to determined shared kernels. The shared kernels are determined based on a transformation of a real-valued continuous-time pass band Volterra series without pruning of kernels.
According to this aspect, in some embodiments, the shared kernels are determined based on the transformation of the real-valued continuous-time pass band Volterra series by steps that include transforming the real-valued continuous time pass band Volterra series to a multi-frequency complex-valued envelope series. The multi-frequency complex-valued envelope signal is transformed to a continuous-time pass band-only series, which is then transformed to a continuous-time baseband equivalent series. The continuous- time baseband equivalent signal is discretized to produce the discrete base band equivalent Volterra series. Shared kernels are identified, each shared kernel having distortion products in common with another shared kernel. In some embodiments, transforming the continuous- time pass band-only signal to a continuous- time baseband equivalent signal includes expressing the continuous- time pass band-only series in convolution form. The Laplace transform is applied to the convolution form to produce a Laplace domain expression, which is frequency shifted to baseband to produce a baseband equivalent expression in the Laplace domain. The inverse Laplace transform is applied to the baseband equivalent expression to produce the continuous-time baseband equivalent series. In some embodiments, a number of terms in the Laplace domain
expression are reduced via symmetry. In some embodiments, terms of the Laplace domain expression are grouped based on frequency intervals where distortion terms are not zero. In some embodiments, discretizing the continuous- time baseband equivalent series to produce the discrete base band equivalent Volterra series includes truncating the continuous-time baseband equivalent series to a finite non-linearity order, and expressing the truncated series as summations of non-linear distortion terms, with upper limits of the summations being memory depths, different memory depths being assignable to different ones of the summations. In some embodiments, a distortion term is a group of distortion products multiplied by a shared kernel.
According to another aspect, embodiments include a digital pre-distorter (DPD) system for pre-distortion of one of a wideband signal and an intra-band carrier aggregated signal. The system includes a Volterra series DPD modelling unit. The DPD modelling unit is configured to receive the one of the wideband signal and the intra-band carrier aggregated signal. The DPD modelling unit is also configured to calculate a discrete base band equivalent, BBE, Volterra series based on the received signal. The series has distortion products grouped according to determined shared kernels. The shared kernels are determined based on a transformation of a real-valued continuous-time pass band Volterra series without pruning of kernels.
According to this aspect, the DPD may further comprise a power amplifier configured to produce an output in response to the one of the wideband signal and the intra-band carrier aggregated received signal. The output of the power amplifier is provided to the Volterra series DPD modelling unit The Volterra series DPD modeling unit is configured to compute the shared kernels based on the output of the power amplifier. In some embodiments, the DPD system further comprises a transmitter observation receiver configured to sample the output of the power amplifier and provide the sampled output to the Volterra series DPD modelling unit. In some embodiments, the distortion products and their associated kernels are determined by transforming the real-valued continuous time pass band Volterra series to a discrete base band equivalent
Volterra series according to a series of steps that include: transforming the real- valued continuous time pass band Volterra series to a multi-frequency complex- valued envelope series; transforming the multi-frequency complex-valued envelope signal to a continuous-time pass band-only series; transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent series; discretizing the continuous-time baseband equivalent signal to produce the discrete base band equivalent Volterra series. The shared kernels are identified using the envelopes of the input and output signals of the PA. In some embodiments, transforming the continuous- time pass band-only signal to a continuous-time baseband equivalent signal includes the following steps:
expressing the continuous-time pass band-only series in convolution form;
applying a Laplace transform to the convolution form to produce a Laplace domain expression; frequency shifting the Laplace domain expression to baseband to produce a baseband equivalent expression in the Laplace domain; and applying an inverse Laplace transform to the baseband equivalent expression to produce the continuous-time baseband equivalent series. In some embodiments, a number of terms in the Laplace domain expression are reduced via symmetry. In some embodiments, transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent signal further includes grouping terms of the Laplace domain expression based on frequency intervals where distortion terms are not zero. In some embodiments, discretizing the continuous-time baseband equivalent series to produce the discrete base band equivalent Volterra series includes: truncating the continuous- time baseband equivalent series to a finite non-linearity order; and expressing the truncated series as summations of non-linear distortion terms, with upper limits of the summations being memory depths assigned to each order of the non-linear distortion terms. In some embodiments, a distortion term is a group of distortion products multiplied by a shared kernel.
According to another aspect, embodiments include a Volterra series digital pre-distorter, DPD, modelling unit. The DPD modelling unit includes a memory module configured to store terms of a discrete base band equivalent,
BBE, Volte rra series. The series is based on one of a received wideband signal and an intra-band carrier aggregated signal. Also, a grouping module is configured to group distortion products of the series according to determined shared kernels. The DPD modelling unit also includes a shared kernel determiner configured to determine the shared kernels based on a
transformation of a real-valued continuous-time pass band Volterra series without pruning of kernels. Also, a series term calculator is configured to calculate the terms of the discrete base band equivalent Volterra series, the terms being the distortion products multiplied by their respective shared kernels.
In some embodiments, the shared kernel determiner is further configured to determine the shared kernels via a least squares estimate based on the one of the wideband signal and the intra-band carrier aggregated received signal, and an output of a power amplifier. In some embodiments, the kernels and distortion products are derived from the real-valued continuous-time pass band Volterra series by: transforming the real-valued continuous time pass band Volterra series to a multi-frequency complex-valued envelope series; transforming the multi-frequency complex-valued envelope signal to a continuous-time pass band-only series; transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent series; discretizing the continuous-time baseband equivalent signal to produce the discrete base band equivalent Volterra series; and identifying the shared kernels, each shared kernel having distortion products in common. In some embodiments, transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent signal includes: expressing the continuous-time pass band-only series in convolution form; applying a Laplace transform to the convolution form to produce a Laplace domain expression; frequency shifting the Laplace domain expression to baseband to produce a baseband equivalent expression in the Laplace domain; and applying an inverse Laplace transform to the baseband equivalent expression to produce the continuous-time baseband equivalent series.
BRIEF DESCRIPTION OF THE DRAWINGS
A more complete understanding of the present invention, and the attendant advantages and features thereof, will be more readily understood by reference to the following detailed description when considered in conjunction with the accompanying drawings wherein:
FIG. 1 is a diagram of low pass transformation for a low pass equivalent (LPE) power amplifier;
FIG. 2 is a column chart showing numbers of kernels for different exemplary Volte rra series implementations;
FIG. 3 is a spectrum of a PA output for several different linearization methods for a wideband input signal;
FIG. 4 is a spectrum a PA output for several different linearization methods for an intra-band carrier aggregated input signal, for a first PA;
FIG. 5 is a spectrum a PA output for several different linearization methods for an intra-band carrier aggregated input signal, for a second PA;
FIG. 6 is a block diagram of an exemplary power amplification system having a digital pre-distorter modelling unit implementing the BBE Volterra model presented herein;
FIG. 7 is a block diagram of the pre-distorter modelling unit of FIG. 6; FIG. 8 is a flowchart of an exemplary process for modelling a power amplifier fed by a wideband or intra-band carrier aggregated input signal;
FIG. 9 is a flowchart of an exemplary process for transforming a real- valued continuous time pass band Volterra series to a discrete base band equivalent Volterra series in which shared kernels are identified;
FIG. 10 is a flowchart of an exemplary process of transforming the continuous- time pass band-only signal to a continuous-time baseband equivalent signal; and
FIG. 11 is a flowchart of an exemplary process of discretizing the continuous- time baseband equivalent series to produce the discrete base band equivalent Volterra series.
DETAILED DESCRIPTION
Before describing in detail exemplary embodiments that are in accordance with the present invention, it is noted that the embodiments reside primarily in combinations of apparatus components and processing steps related to use of baseband equivalent Volterra series for digital pre-distortion of power amplifiers driven by wideband or intra-band carrier aggregated signals.
Accordingly, the system and method components have been represented where appropriate by conventional symbols in the drawings, showing only those specific details that are pertinent to understanding the embodiments of the present invention so as not to obscure the disclosure with details that will be readily apparent to those of ordinary skill in the art having the benefit of the description herein.
As used herein, relational terms, such as "first" and "second," "top" and "bottom," and the like, may be used solely to distinguish one entity or element from another entity or element without necessarily requiring or implying any physical or logical relationship or order between such entities or elements.
A base band equivalent (BBE) expression of the Volterra series for linearization of a power amplifier driven by one of a wideband signal and an intra-band carrier aggregated signal is provided and the procedure to derive this new expression is described below in detail. A wideband signal and an intra- band carrier aggregated signal are distinguishable from a multi-band signal. A multi-band signal has component carriers in different communication bands. Communication bands may be defined by the third generation partnership project (3GPP) standards. In contrast to a multi-band signal, a wideband signal has component carriers that are contiguous in the same communication band. An intra-band carrier aggregated signal has carriers that may be contiguous or non-contiguous, but are in the same communication band.
The complexity and performance of the compact BBE Volterra series disclosed herein is compared to those of the widely used LPE Volterra series expression. For purposes of clarity, a simplified version of the LPE Volterra series given by equation (4), with nonlinearity order NL=3and memory effects order M = 1 ( i. e. x(t— rs)), will be used henceforward as per equation (5).
1
y(n) = ^ h^in-i^
ί! = 0
1 1 1
+ Σ Σ Σ ' (n - - )x*(n - i3)
^=0 i2 = H =0
= h0x(n)
+¾1x(n— 1)
+¾000x(n)x(n)x*(n)
+¾001x(n)x(n)x*(n— 1)
+h010x(n)x(n— l)x*(n)
+^o,i,i^(n)^(n— i)x*(n— 1)
+h110x(n— l)x(n— l)x*(n) +¾11 x(n- l)x(n- l)x*(n- 1). (5)
The implementation of the BBE Volterra series, suitable for the behavioral modeling and predistortion of RF PAs, uses a number of steps that are described in detail as follows:
Step 1: Continuous-time real-valued Volterra series modeling. The Volterra series framework is initially used to describe the relationship between the real passband signals located at the system input and output stages as follows: y( = hv (τΐ'■■ ' τν) \ (ί ~ τί)άτί (6)
p=l ' J 7=1
where x(t) and y (t) represent the PA input and output RF signals and NL is the nonlinearity order.
Step 2: Real-valued to complex-valued envelope signal transformation. The input signal x(t) is a band limited modulated signal that can be expressed as a function of its envelope around the carrier angular frequency ω as follows: x(t) = Re{x(t)e^^}
= i (x*(t)e ~ it)ct + X(t)e^ct), (7)
In the equation above, x(t) represents the complex baseband envelope signal that modulates the angular frequency ω . Substituting equation (7) into equation (6) yields an expression relating the output signal y(t) to x(t) and <¾. as follows: y(t) = f(m, e±^ct); V e {1, . - , NL], (8)
where the describing function / is used to represent the real valued Volterra series of equation (6). Since the output signal y(t) in equation (8) has been altered by the application of a nonlinear function to a band-limited RF signal, it now contains several spectrum components that involve multiple envelopes, represented by yp (t). The spectrum components modulate the mixing products of ω as follows:
p=NL p = -NL
+ ^ (y2 *(0 ■ e ~i2 )ct + 2( ^'2Wct)
(9)
where y0(t) denotes the envelope at DC, y1(t) denotes the envelopes of the first harmonic in-band signals, y2(0 denotes the envelopes of the second order harmonic signals and yWL(t) represents the NLth harmonic signal.
Step 3: Multi-frequency to passband only transformation. By equating the terms that share the same frequency range (fundamental, mixing product) from the right sides of equations (8) and (9), a multi-harmonic model is derived which consists of several distinct equations relating the output envelopes yp (t) to (t) and <¾.. Since interest is in the relationship between the envelopes of the output and input signals around the carrier frequency, only the pass band component of the PA output is considered in the equation below:
where y^t) = 1 (t)e,"'ct. The output signal envelope yVb(t) in equation (10) has the form of a band limited signal as given in equation (7).
Step 4: Continuous time passband Volterra series. The following derivation is presented for the term around eJ£Jct only; however, the same equation can be used to derive the conjugate term. The PA passband signal y^t) can be modeled as a summation of the Volterra series nonlinear order responses yi,2fc+i (t), or as the nonlinear order envelope responses y (t). around the angular frequency ωε as follows:
Note that only odd powered terms are retained; even terms are discarded since they do not appear around the carrier frequency. Equating the terms on the right sides of the expanded equations (8) and (10) yield a continuous BBE Volterra series that expresses ^ 2fe+i (i) as a function of x(t) and <uc. Below is the expression derived using y1 1(t) and y1 3(t):
yi,i( = J Λι(τι)■ x(.t - τι β^-τ,) . dTi (12) i,3( = J J J l3(Ti,T2,T3)■{
}·άτ3άτ2άτ1 (13)
te75; Continuous-time passband to baseband equivalent
transformation. In order to be computed by a digital processor, with manageable complexity, a passband model should be transformed to a BBE model. A BBE model includes the RF nonlinear dynamic distortion and computes the terms in baseband with a low sampling rate. To translate the passband model to baseband, the continuous time passband Volterra series of equations (10) and (12) is first rewritten using convolution form. The convolution form for y1}1(t) is given by:
As the kernel h3 is tri-variate, /ι3(τ1,τ2,τ3), and the outputy(t) is mono- variate, the output function is re-assigned as follows:
y( = y( , t2, t3)|tl=t2=t3=t y(t, t, t) .
The convolution form is given as:
y1,3(t1,t2,t3) = h3(t1,t2,t3)*{
(X(t1)e^^)(X(t2)e^^)(X(t3)e
+(χ(ί1)β'ωΛ)(^(ί2)^'ωΛ)*( (ί3)β^ω^3)
+(χ(ί1)β^ι)*(χ(ί2)β^ω^)(χ(ί3)β^ωΛ)}. (15)
Next, the Laplace transformation is applied to the convolution form of the Volterra series nonlinear order responses y12fe+i- The Laplace domain representation of equations (14) and (15) is given by:
= H1(s)-?(s-;'iuc). (16)
Yi,3 Oi, ¾ s3) = £ ( i,3 (fi, , hi)
= £(/3(ί1,ί2,ί3))■{ £(( (ί1)^'ωΛ)(^(ί2)^ωΑ)( (ί3)β^ωΛ)*) +£ (( (t1)e^^)( (i2)e^^)*( (i3)e^^))
+£((χ(ί1)^'ωΛ)*( (ί2)^'ωΛ)( (ί3)β^ωΛ))}
= £(/3(ί1,ί2,ί3))■ { £(χ(ί1)β^ωΛ)£(χ(ί2)β^ω^)£((χ(ί3)β^ωΛ)*) +£(x(t1)e^^)£((x(t2)e^^)*)£(x(t3)e +£((x(t1)e^ti)*)£(x(t2)e^^)£(x(t3)e = H3(s1,s2,s3)■ { -ja) c)X(s2 -ja) c)X*((s3— yac)*)
- ) *((s2 - )*)X(s3 -;<uc)
+-?*((s1 -ja>cY)X(s2 -ja>c-)X(s3 -j >c)} (17)
The above derivation uses the distributive and product transform properties of the multi-variate Laplace transform. The BBE form of equations (16) and (17) can be derived by applying a frequency translation of ja)c to the constituent terms Ύ1 and Y13 as follows:
¾,ιΟ) = >Ί,ιΟ+ )
= H1(s+joc)X(s
= H^Xds).
¾,3(si' ¾ ss) = ^I.SOI + /<¾., s2 + /<¾., s3 + /<¾.)
= H3 (Si + )ω0, s2 + )ω0, s3 + y wc)■ {
X(Sl)X(s2)X*((s3)*)
+X(Sl)X*((s2Y)X(s3)
+X*((s1Y)X(s3)X(s2)}
= H3(s1,s2,s3)-{
X(Sl)X(s2)X*((s3)*)
+X(Sl)^*((s2)*)X(s3)
+X*((s1)*)X(s3)X(s2)} (18)
where and ?13 designate the BBE outputs and ¾ and/?3 denote the BBE kernels, respectively. Applying the inverse Laplace to equation (18) yields the continuous-time BBE model as follows:
(x(t - τι))( (ί - τ2))( (ί - τ3))
+(x(t - Τ!))(χ(ί - τ2))*(χ(ί - τ3))
+(χ(ί-τ1))*(χ(ί-τ2))(χ(ί-τ3))
(19)
Step 6: Discrete-time baseband equivalent Volterra series model. Using the signal and system causality, the integral bounds of the model of equation (19) are adjusted to∞→ t and -∞→ o, respectively. The model is then given by:
(x(t - T!))( (i - τ2))( (ί - τ3))*
+(χ(ί - τι))(χ(ί - τ2))*(χ(ί - τ3))
+(χ(ί - τι))*(χ(ί - τ2))(χ(ί - τ3))
}■ άτ3άτ2άτ1 (20)
Using the fading memory assumption for the steady-state response of the
PA (t > T∞) (the transient-response, time-invariant Volterra series is defined as t < T∞), the corresponding model can be represented as:
Too Too Too
i,3 ( = J 0 J 0 I 0 ¾(¾,¾,¾)■{ (x(t - T )(x(t - T2))(x(t - τ3 Υ +(x(t - Ti))(X(t - τ2))*(χ(ί - τ3)) +(χ(ί - Ti))*(X(t - τ2))(χ(ί - τ3))
}■ άτ3άτ2άτ1 (21)
Using the symmetry of the integrated function (distortion components are symmetrical and Volterra kernels can be symmetrized), the model can be simplified to:
ι( ( -J)■ x^t— τ-j)■ dr1
(x(t - Tl))(X(t - τ2))(χ(ί - τ3)) +(x(t - Ti))(X(t - τ2))*(χ(ί - τ3))
+(χ(ί -τ1))*(χ(ί-τ2))(χ(ί -τ3))
}■ άτ3άτ2άτ1 (22)
Digitization of the model yields:
M3 M3 M3 ί1=012=1! i3=i2
Xi(n, = x(n— i-j
x3(n, ilt i2, j3) = x(n - i^xfji - i2)x*(n - i3)
+x(n— i^x^n— i2)x(n— i3)
+x*(n— i^xin— i2)x(n— i3) (23)
where MI,M3 denote the memory depth of the first and third, order Volterra series. Similarly to the above derivation, it can be shown that the fifth order Volterra kernel, li5(n), is given by:
Ks ) ■ Xs(n> > > > > is)
Xs(n> > i-2> i-3> i-4> is) =
x(n— i^)x(n— i2)x(n— i3)x*(n— i4)x*(n— i5)
+x(n - - i^xin - - i2)x*(n - - i3 x(n - - i4)x*(n -is)
+x(n - - i-i)x*(i - i2 x(n - - i3 x(n - - i4)x*(n -is)
+x*(n — i^xin - i2 x(n - - i3 x(n - - i4)x*(n -is)
+x(n - - i^^in - i2 x(n - - i3)x*(n — i4)x(n -is)
+x*(n — i^xin - i2)x(n - - i3)x*(n - j4)x(n -is)
+x(n - - i-i)x*(i - i2)x*(n - i3)x(n — i4)x(n -is)
+x*(n — i^xin - i2)x*(n - i3)x(n - j4)x(n -is)
where M denotes the memory depth of the fifth order Volterra series. The complex valued BBE Volterra series in equations (23) and (24) up to nonlinear distortion products 5 is then given by:
9i (n) = y1}1 (n) + y1}3 (n) + 1>5 (n) (25)
Higher order distortion products can be similarly derived. In addition, only odd powered terms are retained. The terms x3(n, ilt i2, i ) and
x5(n, ilt i2, i , i4, is) are linear combinations of 3 third, and 10 fifth, order distortion products. An example of the BBE Volterra model of equation (25) is given for NL = 3 and M = 1 as follows
1 1
y(n) = -h0x(n) + -h1x(n - 1)
3
+-ft000x(n)x(n)x*(n)
8
1
+ o^o,o,i(2 x(n)x(n— l)x*(n) + x(n)x(n)x*(n— 1))
8
1
+— h011(2 x(n)x(n— l)x*(n— 1) + x(n— l)x(n— l)x*(n))
3
+— hx ! ix( — l)x(n— l)x*(n— 1) (26)
8
Comparison of the classical expression of equation (5) with the newly proposed BBE Vol terra expression of equation (26) described herein results in several conclusions. For example, the proposed BBE Volte rra model effectively combines a number of distortion components to share a unique kernel, as per equations (23) and (24). For example, the two distortion components
x(n)x(n— l) * (n— 1) and x(n— l)x(n— l) *(n) share the same kernel ¾o,i,i m the BBE Volterra model of equation (26), however, they have different kernels (fio,i,i and¾i,i,o< respectively) in the LPE Volterra model of equation (5). More generally, several distortion products, which were assigned distinct kernels in equation (4), are linearly combined to form a novel term in the new model space proposed by equation (25), and hence, share the same kernel.
Table I summarizes the number of kernels needed for four different example modeling scenarios where the nonlinearily order has been set to seven and the memory depth M, equals Ml = M3 = MB = M. The expressions for the number of coefficients were used to plot FIG. 2. The plot shows the number of kernels for different memory depths for DDE LPE Volterra r=2 10, r=3 12, compact BBE Volterra 14 and Classical LPE Volterra 16.
As can be seen, the newly proposed, compact model involved a significantly lower number of coefficients than the full classical LPE Volterra series, especially as the memory depth increased. This results in lower complexity and cost of implementation. However, the pruned LPE Volterra using the DDR approach also required a smaller number of coefficients as it discarded distortion products with a dynamic order > three.
TABLE I
Thus, different orders of the nonlinear dynamic response, yi,i(n), y 1,3 (1) and y1 5(n) are results of different nonlinear integrals yi i(t), y1 3 (t) and yli5(t).
Depending on the dynamic property of the PA, and depending on the numerical integration scheme used, different integrals may call for different memory depths. Hence, in the proposed BBE Volterra model of equation (25), three memory depth parameters may be defined, Ml, 3 and MB, instead of the single memory depth parameter, M, used in the LPE Volterra model of equation (4). This allows a flexible control of the model's dynamics based on a theoretical framework. Note also that, due to the symmetry properties of the integrals in each order of the nonlinear dynamic response, the different summation operators in yi,i (w), ylj3 (n) and y1 5(n) use the same memory depth Ml, 3 and 5 respectively.
The signal transformations used to derive the BBE expression of equation (25) preserve the linear property of the Volterra series with respect to its coefficients. Therefore, commonly used model identification algorithms such as the least square error (LSE) can be applied to identify the kernels in equation
(25) for a given RF PA. Equation (27) suggests the expression used to deduce the LSE solution of equation (25)
X - h = Y (27)
where X, h and Y denote the input signal non-square matrix, the vector that contains the unknown kernels and the output signal vector as described below:
In the above expression, N represents the total number of measurement data. The solution to the problem defined in equation (27) is:
h = (X1 - ΧΥ1 - X1■ Y (29)
where h is the estimate of h.
To summarize the approach, a discrete BBE Volterra series is generated based on a received wideband or intra-band carrier aggregated signal. The series has distortion products grouped according to determined shared kernels. The shared kernels are determined based on a transformation of a real-valued continuous-time pass band Volterra series without pruning of kernels. The transformation includes transforming the real-valued continuous time pass band Volterra series to a multi-frequency complex-valued envelope series. The multi- frequency complex-valued envelope signal is then transformed to a continuous- time pass band-only series. The continuous-time pass band-only signal is transformed to a continuous-time baseband equivalent series. The continuous- time baseband equivalent signal is discretized to produce the discrete base band equivalent Volterra series. Shared kernels of the discrete base band equivalent Volterra series are identified, where each shared kernel has distortion products in common with another shared kernel.
To assess the performance of the newly proposed compact BBE Volterra model, the model was used to linearize different PAs driven with wideband and intra-band carrier aggregated signals. The linearization performance was compared with that of the classical LPE Volterra model and the DDR pruned Volterra model.
By way of example, for the wideband signal linearization test, a 200 W Doherty PA (DPA) driven with a 2C 1001 WCDMA signal was linearized. The device under test (DUT) was excited with a 1 ms signal sampled at 100 MHz, where, out of 100k points, 10k were used to train the proposed model. Being non-pruned Volterra models, the nonlinearity order and memory depth for both the proposed model and the LPE Volterra model were set to be equal. Hence, the memory depths of the different nonlinear dynamic responses of the BBE Volterra model were set at the same M1=M3=M5=M. However, the DDR pruned Volterra model parameters were varied separately for optimum performance. The different models' parameters and the resulting model complexities are summarized in Table II. It can be seen that both the DDR pruned Volterra and the BBE Volterra models have a significantly lower number of kernels than the classical LPE Volterra model. Table II also summarizes the linearization results of the DUT.
TABLE II
COMPLEXITY AND LINEARIZATION PERFORMANCE OF A DOHERTY POWER AMPLIFIER USING DIFFERENT VOLTERRA FORMULATIONS
The spectrum plot of the PA response without predistortion, and under the linearization of the different Volterra formulations, is presented in FIG. 3. It shows that all of the tested models successfully reduced the error vector magnitude (EVM) from 6 % to about 1.2 % and the adjacent channel leakage ratio (ACLR) from -38 dBc to about -55dBc. It can be concluded that both the DDR pruning of the LPE Volterra model and the kernel sharing of the BBE
Volterra model did not affect the linearization performance of the classical LPE Volterra formulation.
In order to validate the proposed model under recent 4G communication signals, an intra-band carrier aggregated signal was synthesized. This signal was mixed-standard, with a bandwidth of 40 MHz and a PAPR of 9.8 dB, containing a 15 MHz Long Term Evolution (LTE) signal and a 15 MHz 101 Wideband Code Division Multiple Access (WCDMA) signal separated by a 10 MHz guard band. Two different RF PA demonstrators were used as DUT. The first was a single ended Class J 45 W GaN amplifier and the second was a laterally diffused
MOSFET (LDMOS) 200 W Doherty PA. The DUT was excited with a 1 ms signal sampled at 122.88 MHz, where, out of 122.88k points, 10k were used to train the proposed model.
Tables III and IV summarize the linearization results found with each amplifier using both the proposed compact and the pruned LPE Volterra models. For all testing, the nonlinearity order was set to seven and the memory depth was set two and three.
TABLE III COMPLEXITY AND LINEARIZATION PERFORMANCE OF A CLASS AB GAN POWER AMPLIFIER USING DIFFERENT VOLTERRA FORMULATIONS
Number of
Model Parameters EVM ACLR (dBc)
Coefficients
Without
5.5% -35
DPD
Classical
LPE NL=1, M=2 231 2% -48
Volterra
DDR LPE NL=7, =3,
91 2.1 % -48
Volterra r=2
BBE
NL=1, M=2 70 2% -48
Volterra
TABLE IV
COMPLEXITY AND LINEARIZATION PERFORMANCE OF LDMOS DOHERTY POWER AMPLIFIER USING DIFFERENT VOLTERRA FORMULATIONS
According to Table III, both the pruned LPE Volterra and the compact models allowed for similar linearization performance as the EVM was reduced from 5.5% to 2% and the ACLR was increased from 35 dB to about 50 dB (after application of the linearization). This linearization capacity was successfully achieved despite the significant reduction in the number of kernels required by the compact model (91 in the pruned LPE Volterra model versus 70 in the proposed model).
It is worth mentioning that the full LPE Volterra model required 231 kernels. The linearization performance was further confirmed by examining the spectrum of the 45 W GaN PA output signals with and without predistortion. This is shown in FIG. 4, which is the spectrum of the single-ended GaN PA driven with the 40 MHz multi- standard intra-band carrier aggregated signal using the different Volterra series formulations and without linearization.
In the case of the second RF PA (see Table IV), the pruned LPE Volterra and the proposed model showed different linearization performance. While both the classical LPE Volterra model and the compact BBE Volterra model achieved similar linearization results in terms of EVM (1.9%) and ACLR (-50 dBc), the DDR LPE
Volterra led to a 0.4% higher EVM and 4 dB lower ACLR. The superiority of the linearization capacity of the new compact model is confirmed in FIG. 5, which displays the spectra of the LDMOS Doherty PA excited by the 40 MHz multi- standard intra-band carrier aggregated signal using different Volterra series formulations and without linearization.
A noticeable residual out of band spectrum regrowth in the signal obtained using the pruned DDR LPE can be observed. As previously mentioned, the DDR pruning approach may suffer from reduced linearization capability when the dynamic order is set to two as it eliminates a large number of distortion products. The compact BBE Volterra model described herein does not suffer from this limitation. All distortion products are maintained. The pruned LPE Volterra linearization performance could be further improved if the dynamic order R were increased further, however, this would result in a larger number of kernels, e.g., the number of kernels would be 162 for r=3 M=2.
FIG. 6 is a block diagram of a power amplification system 18 having a digital pre-distorter modelling unit 3 implementing the BBE Volterra model presented herein. The power amplification system 18 includes a digital pre- distorter 20. The DPD 20 receives the identified shared kernels from the pre- distorter modelling unit 34, and pre-distorts the input signal xto produce a pre- distorted signal. The pre-distorted signal is converted to analog by a digital to analog (DAC) converter 22, low pass filtered by a filter 24, and mixed to radio frequency (RF) by a mixer 26 to prepare the signal for amplification by an RF PA amplifier 28. A directional coupler 30 couples the output of the RF PA 28 to a transmitter observation receiver (TOR) 32 The TOR 32 samples the output of the power amplifier 28 in each band and produces a TOR output signal, y. The TOR output signal is used by the DPD modeling unit 34 to derive the kernel vector h according to equation (27). The DPD modelling unit 34 calculates a discrete baseband equivalent Volterra series having distortion products grouped according to determined shared kernels, where the shared kernels are based on a transformation of a real-valued continuous-time pass band Volterra series without pruning of kernels.
FIG. 7 is a detailed view of the DPD modelling unit 34, which includes a processor 36 in communication with a memory module 38. The processor 36 may be implemented as a micro-processor operating according to computer instructions organized as a group of software modules, or each module of the processor 36 may be implemented by application specific integrated circuitry. The memory 38 may be implemented as random access memory (RAM), for example, and may contain non-volatile components, such as read only memory (ROM) that stores programmatic code to implement the functions described herein. The DPD modelling unit 34 receives the input signal xand the TOR output signal 'from the transmitter observation receiver 32 and derives the modelling vector h according to equation (27). The modelling vector h is input to the DPD 20 to pre-distortthe input signal xto produce a pre-distorted signal. The processor 36 includes a grouping module 40, a shared kernel determiner 42 and a series term computer 44. The grouping module 40 is configured to group distortion products of the series according to determined shared kernels, as shown in equation (25). The shared kernel determiner 42 is configured to determine the shared kernels based on a transformation of a real-valued continuous- time pass band Volterra series without pruning of kernels, as described above with reference to Step 6. The series term calculator 44 is configured to calculate the terms of the discrete base band equivalent Volterra series according to equation (27), the terms being the distortion products multiplied by their respective shared kernels. The memory module 38 is configured to store terms of the discrete BBE Volterra series 46, generated by the processor 36.
FIG. 8 is a flowchart of an exemplary process for modelling a power amplifier 40 fed by a wideband or intra-band carrier aggregated input signal using the discrete BBE Volterra series described herein for improved
linearization with reduced complexity. The input signal is received by a digital pre-distorter 28 (block SlOO). A discrete BBE Volterra series is generated by the DPD modelling unit 34 based on the received input signal (block S102), according to equations (25) and (27). The series has distortion products that
are grouped by the grouping module 40 according to determined shared kernels determined by shared kernel determiner 42, as shown in equation (25). The shared kernels are determined based on a transformation of a real-valued continuous- time pass band Volterra series without pruning of kernels, according to Steps 1 through 6 above.
FIG. 9 is a flowchart of an exemplary process for transforming a real- valued continuous time pass band Volterra series to a discrete base band equivalent Volterra series according to Steps 1 through 6 above in which shared kernels are identified as set out in block S102 of FIG. 8. The real-valued continuous time pass band Volterra series is transformed to a multi-frequency complex-valued envelope series, as in Steps 1 and 2, described above (block S104). The multi-frequency complex-valued envelope signal is transformed to a continuous- time pass band-only series, as in Step 3 described above (block S106). The continuous-time pass band-only signal is transformed to a continuous- time baseband equivalent series, as in Steps 4 and 5 described above (block S108). The continuous-time baseband equivalent signal is discretized to produce the discrete base band equivalent Volterra series, as in Step 6 described above, (block S110). Shared kernels of the discrete base band equivalent Volterra series are identified, where a shared kernel has distortion products in common with another shared kernel, also described as in Step 6 described above (block SI 12).
FIG. 10 is a flowchart of an exemplary process of transforming the continuous- time pass band-only signal to a continuous-time baseband equivalent signal as shown in block SI 08 of FIG. 9, as in Step 5 described above. The continuous-time pass band-only series is expressed in convolution form (block S114). Then, the Laplace transform is applied to the convolution form to produce a Laplace domain expression (block S116). A number of terms in the Laplace domain expression may be reduced based on symmetry (block S118). The Laplace domain expression is frequency- shifted to baseband to produce a baseband equivalent expression in the Laplace domain (block S120). An inverse
Laplace transform is applied to the baseband equivalent expression to produce the continuous- time baseband equivalent series (block S122).
FIG. 11 is a flowchart of a process of discretizing the continuous-time baseband equivalent series to produce the discrete base band equivalent Volterra series, as shown in block SI 10 of FIG. 9, as in Step 6 described above. The process includes truncating the continuous-time baseband equivalent series to a finite non-linearity order (block S124). The process also includes expressing the truncated series as summations of non-linear distortion terms, with upper limits of the summations being memory depths assigned to each order of the non-linear distortion terms (block S126).
The PA modelling method described herein using the discrete BBE Volterra series is inherently compact and calls for significantly fewer coefficients then its LPE counterpart and therefore avoids the accuracy-compromising pruning transformations widely applied to classical formulations in the literature. The validation described herein proved the excellent modeling and linearization performance of the method when compared to the classical non-pruned Volterra model. The new formulation also outperformed pruned Volterra models while using a lower number of coefficients.
It will be appreciated by persons skilled in the art that the present invention is not limited to what has been particularly shown and described herein above. In addition, unless mention was made above to the contrary, it should be noted that all of the accompanying drawings are not to scale. A variety of modifications and variations are possible in light of the above teachings without departing from the scope of the following claims.
Claims
1. A method of modelling a power amplifier fed by one of a wideband signal and an intra-band carrier aggregated signal, the method comprising:
receiving the one of the wideband signal and the intra-band carrier aggregated signal (S 100);
generating a discrete base band equivalent, BBE, Volterra series based on the received signal, the series having distortion products grouped according to determined shared kernels (SI 02); and
the shared kernels being determined based on a transformation of a real-valued continuous-time pass band Volterra series without pruning of kernels (S102).
2. The method of Claim 1 , wherein the shared kernels are determined based on the transformation of the real- valued continuous-time pass band Volterra series by: transforming the real-valued continuous time pass band Volterra series to a multi- frequency complex- valued envelope series (SI 04);
transforming the multi-frequency complex-valued envelope signal to a continuous-time pass band-only series (SI 06);
transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent series (S108);
discretizing the continuous-time baseband equivalent signal to produce the discrete base band equivalent Volterra series (SI 10); and
identifying the shared kernels, each shared kernel having distortion products in common with another shared kernel (SI 12).
3. The method of Claim 2, wherein transforming the continuous-time pass band- only signal to a continuous-time baseband equivalent signal includes:
expressing the continuous-time pass band-only series in convolution form
(S114);
applying a Laplace transform to the convolution form to produce a Laplace domain expression (SI 16);
frequency shifting the Laplace domain expression to baseband to produce a baseband equivalent expression in the Laplace domain (SI 20); and
applying an inverse Laplace transform to the baseband equivalent expression to produce the continuous-time baseband equivalent series (S122).
4. The method of Claim 3, wherein a number of terms in the Laplace domain expression are reduced via symmetry (S 118).
5. The method of Claim 3, further comprising grouping terms of the Laplace domain expression based on frequency intervals where distortion terms are not zero.
6. The method of Claim 2, wherein discretizing the continuous-time baseband equivalent series to produce the discrete base band equivalent Volterra series includes:
truncating the continuous-time baseband equivalent series to a finite non- linearity order (SI 24); and
expressing the truncated series as summations of non-linear distortion terms, with upper limits of the summations being memory depths, different memory depths being assignable to different ones of the summations (S126).
7. The method of Claim 6, wherein a distortion term is a group of distortion products multiplied by a shared kernel.
8. A digital pre-distorter (DPD) system for pre-distortion of one of a wideband signal and an intra-band carrier aggregated signal, comprising:
a Volterra series DPD modelling unit (34), the DPD modelling unit (34) configured to:
receive the one of the wideband signal and the intra-band carrier aggregated signal (SI 00): and
calculate a discrete base band equivalent, BBE, Volterra series based on the received signal, the series having distortion products grouped according to determined shared kernels (SI 02); and
the shared kernels being determined based on a transformation of a real-valued continuous-time pass band Volterra series without pruning of kernels (S102).
9. The DPD system of Claim 8, further comprising:
a power amplifier (28), the power amplifier (28) configured to produce an output in response to the one of the wideband signal and the intra-band carrier aggregated received signal, the output of the power amplifier (28) provided to the Volterra series DPD modelling unit (34), the Volterra series DPD modeling unit (34) being configured to compute the shared kernels based on the output of the power amplifier (28).
10. The DPD system of Claim 9, further comprising a transmitter observation receiver (32) configured to sample the output of the power amplifier and provide the sampled output to the Volterra series DPD modelling unit (34).
11. The DPD system of Claim 8, wherein the distortion products and their associated kernels are determined by:
transforming the real-valued continuous time pass band Volterra series to a multi- frequency complex- valued envelope series (SI 04);
transforming the multi-frequency complex-valued envelope signal to a continuous-time pass band-only series (SI 06);
transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent series (S108);
discretizing the continuous-time baseband equivalent signal to produce the discrete base band equivalent Volterra series (SI 10); and
identifying the shared kernels, each shared kernel having distortion products in common with another shared kernel (SI 12).
12. The DPD system of Claim 11, wherein transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent signal includes:
expressing the continuous-time pass band-only series in convolution form
(S114);
applying a Laplace transform to the convolution form to produce a Laplace domain expression (SI 16);
frequency shifting the Laplace domain expression to baseband to produce a baseband equivalent expression in the Laplace domain (SI 20); and
applying an inverse Laplace transform to the baseband equivalent expression to produce the continuous-time baseband equivalent series (S122).
13. The DPD system of Claim 12, wherein a number of terms in the Laplace domain expression are reduced via symmetry (SI 18).
14. The DPD system of Claim 12, wherein transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent signal further comprises grouping terms of the Laplace domain expression based on frequency intervals where distortion terms are not zero.
15. The DPD system of Claim 11, wherein discretizing the continuous-time baseband equivalent series to produce the discrete base band equivalent Volterra series includes:
truncating the continuous-time baseband equivalent series to a finite non- linearity order (SI 24); and
expressing the truncated series as summations of non-linear distortion terms, with upper limits of the summations being memory depths assigned to each order of the non-linear distortion terms (S126).
16. The DPD system of Claim 15, wherein a distortion term is a group of distortion products multiplied by a shared kernel.
17. A Volterra series digital pre-distorter, DPD, modelling unit (34), comprising:
a memory module (38), the memory module configured to store terms of a discrete base band equivalent, BBE, Volterra series, the series being based on one of a received wideband signal and an intra-band carrier aggregated signal;
a grouping module (40), the grouping module configured to group distortion products of the series according to determined shared kernels;
a shared kernel determiner (42), the shared kernel determiner configured to determine the shared kernels based on a transformation of a real-valued continuous- time pass band Volterra series without pruning of kernels; and
a series term calculator (44), the series term calculator configured to calculate the terms of the discrete base band equivalent Volterra series, the terms being the distortion products multiplied by their respective shared kernels.
18. The Volterra series DPD modelling unit of Claim 17, wherein the shared kernel determiner (42) is further configured to determine the shared kernels via a least squares estimate based on the received one of the wideband signal and the intra-band carrier aggregated signal, and an output of a power amplifier.
19. The Volterra series DPD modelling unit (34) of Claim 17, wherein the kernels and distortion products are derived from the real-valued continuous-time pass band Volterra series by:
transforming the real-valued continuous time pass band Volterra series to a multi- frequency complex- valued envelope series (SI 04);
transforming the multi-frequency complex-valued envelope signal to a continuous-time pass band-only series (SI 06);
transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent series (S108);
discretizing the continuous-time baseband equivalent signal to produce the discrete base band equivalent Volterra series (SI 10); and
identifying the shared kernels, each shared kernel having distortion products in common (SI 12).
20. The Volterra series DPD modelling unit (34) of Claim 19, wherein transforming the continuous-time pass band-only signal to a continuous-time baseband equivalent signal includes:
expressing the continuous-time pass band-only series in convolution form (S114);
applying a Laplace transform to the convolution form to produce a Laplace domain expression (SI 16);
frequency shifting the Laplace domain expression to baseband to produce a baseband equivalent expression in the Laplace domain (SI 20); and
applying an inverse Laplace transform to the baseband equivalent expression to produce the continuous-time baseband equivalent series (S122).
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP14732395.0A EP2997656A1 (en) | 2013-05-16 | 2014-05-15 | Baseband equivalent volterra series for behavioral modeling and digital predistortion of wideband transmitters |
US14/888,876 US9727677B2 (en) | 2013-05-16 | 2014-05-15 | Baseband equivalent volterra series for behavioral modeling and digital predistortion of wideband transmitters |
Applications Claiming Priority (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201361824075P | 2013-05-16 | 2013-05-16 | |
US61/824,075 | 2013-05-16 | ||
US201361886907P | 2013-10-04 | 2013-10-04 | |
US201361887012P | 2013-10-04 | 2013-10-04 | |
US61/886,907 | 2013-10-04 | ||
US61/887,012 | 2013-10-04 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2014184775A1 true WO2014184775A1 (en) | 2014-11-20 |
Family
ID=50943355
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/IB2014/061476 WO2014184774A1 (en) | 2013-05-16 | 2014-05-15 | Baseband equivalent volterra series for digital predistortion in multi-band power amplifiers |
PCT/IB2014/061477 WO2014184775A1 (en) | 2013-05-16 | 2014-05-15 | Baseband equivalent volterra series for behavioral modeling and digital predistortion of wideband transmitters |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/IB2014/061476 WO2014184774A1 (en) | 2013-05-16 | 2014-05-15 | Baseband equivalent volterra series for digital predistortion in multi-band power amplifiers |
Country Status (4)
Country | Link |
---|---|
US (2) | US9727677B2 (en) |
EP (2) | EP2997655A1 (en) |
CN (1) | CN105393452B (en) |
WO (2) | WO2014184774A1 (en) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015045709A1 (en) * | 2013-09-26 | 2015-04-02 | 日本電気株式会社 | Signal transmission apparatus, distortion compensation apparatus, and signal transmission method |
US9853599B2 (en) * | 2015-12-04 | 2017-12-26 | The Boeing Company | Simultaneous linearization of multiple power amplifiers with independent power |
CA3059091C (en) * | 2017-02-15 | 2023-09-19 | Fadhel Ghannouchi | System and method for rf amplifiers |
US11531639B2 (en) * | 2018-07-26 | 2022-12-20 | The Boeing Company | Pipelined cognitive signal processor |
US10985951B2 (en) | 2019-03-15 | 2021-04-20 | The Research Foundation for the State University | Integrating Volterra series model and deep neural networks to equalize nonlinear power amplifiers |
US10715185B1 (en) * | 2019-03-29 | 2020-07-14 | Intel Corporation | Digital-to-analog conversion system |
CN110276100B (en) * | 2019-05-16 | 2022-11-22 | 中国工程物理研究院电子工程研究所 | Behavior model modeling and implementation method based on Volterra series |
CN111884602B (en) * | 2020-06-29 | 2022-05-20 | 北京邮电大学 | Power amplifier predistortion method based on single-output-node neural network |
FI20205889A1 (en) | 2020-09-15 | 2022-03-16 | Nokia Technologies Oy | Interference mitigation with Multi-Band Digital Pre-Distortion |
US20230291620A1 (en) * | 2020-10-02 | 2023-09-14 | Telefonaktiebolaget Lm Ericsson (Publ) | Method and device(s) for performing digital predistortion on signals for multiband transmission in a wireless communication network |
CN116015509B (en) * | 2022-12-27 | 2023-09-29 | 华中科技大学 | Volterra filter establishment method, device and system applied to DML/DD system |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130094610A1 (en) * | 2011-10-14 | 2013-04-18 | Fadhel Ghannouchi | Digital Multi-band Predistortion Linearizer with Nonlinear Subsampling Algorithm in the Feedback Loop |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6236837B1 (en) * | 1998-07-30 | 2001-05-22 | Motorola, Inc. | Polynomial Predistortion linearizing device, method, phone and base station |
JP4835241B2 (en) * | 2006-04-11 | 2011-12-14 | 株式会社日立製作所 | Digital predistortion transmitter |
US7729446B2 (en) * | 2006-12-01 | 2010-06-01 | Texas Instruments Incorporated | System and method for digitally correcting a non-linear element using a multiply partitioned architecture for predistortion |
CN102017553B (en) * | 2006-12-26 | 2014-10-15 | 大力系统有限公司 | Method and system for baseband predistortion linearization in multi-channel wideband communication systems |
US9705477B2 (en) | 2008-04-30 | 2017-07-11 | Innovation Digital, LLC | Compensator for removing nonlinear distortion |
US8030997B2 (en) * | 2008-11-11 | 2011-10-04 | Philip Brown | Resource efficient adaptive digital pre-distortion system |
US8532215B2 (en) | 2009-09-14 | 2013-09-10 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of Industry, Through The Communications Research Centre Canada | Multi-carrier amplifier linearization system and method |
CN101860326A (en) * | 2010-06-29 | 2010-10-13 | 奥维通信股份有限公司 | Indirect learning predistortion linearized system based on Volterra series |
CN102427336B (en) * | 2011-11-30 | 2015-07-08 | 钱骅 | Radio frequency power amplification system with function of adaptive digital predistortion linearization |
US8855175B2 (en) * | 2012-08-02 | 2014-10-07 | Telefonaktiebolaget L M Ericsson (Publ) | Low complexity all-digital PIM compensator |
-
2014
- 2014-05-15 WO PCT/IB2014/061476 patent/WO2014184774A1/en active Application Filing
- 2014-05-15 CN CN201480040494.4A patent/CN105393452B/en not_active Expired - Fee Related
- 2014-05-15 WO PCT/IB2014/061477 patent/WO2014184775A1/en active Application Filing
- 2014-05-15 US US14/888,876 patent/US9727677B2/en active Active
- 2014-05-15 US US14/889,057 patent/US9659120B2/en not_active Expired - Fee Related
- 2014-05-15 EP EP14730588.2A patent/EP2997655A1/en not_active Withdrawn
- 2014-05-15 EP EP14732395.0A patent/EP2997656A1/en not_active Withdrawn
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130094610A1 (en) * | 2011-10-14 | 2013-04-18 | Fadhel Ghannouchi | Digital Multi-band Predistortion Linearizer with Nonlinear Subsampling Algorithm in the Feedback Loop |
Non-Patent Citations (4)
Title |
---|
AHMED MOHAMED MAHMOUD MOHAMED ET AL: "Reconfigurable Doherty Power Amplifier for Multifrequency Wireless Radio Systems", IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, IEEE SERVICE CENTER, PISCATAWAY, NJ, US, vol. 61, no. 4, 1 April 2013 (2013-04-01), pages 1588 - 1598, XP011498931, ISSN: 0018-9480, DOI: 10.1109/TMTT.2013.2247617 * |
FAROUK MKADEM ET AL: "Wiener G-functionals for nonlinear power amplifier digital predistortion", MICROWAVE SYMPOSIUM DIGEST (MTT), 2012 IEEE MTT-S INTERNATIONAL, IEEE, 17 June 2012 (2012-06-17), pages 1 - 3, XP032217223, ISBN: 978-1-4673-1085-7, DOI: 10.1109/MWSYM.2012.6259707 * |
GHANNOUCHI FADHEL M ET AL: "Distortion and impairments mitigation and compensation of single- and multi-band wireless transmitters (invited)", IET MICROWAVES, ANTENNAS & PROPAGATION, THE INSTITUTION OF ENGINEERING AND TECHNOLOGY, UNITED KINGDOM, vol. 7, no. 7, 15 May 2013 (2013-05-15), pages 518 - 534, XP006045911, ISSN: 1751-8725, DOI: 10.1049/IET-MAP.2012.0663 * |
SEYED AIDIN BASSAM ET AL: "Linearization of Concurrent Dual-Band Power Amplifier Based on 2D-DPD Technique", IEEE MICROWAVE AND WIRELESS COMPONENTS LETTERS, IEEE SERVICE CENTER, NEW YORK, NY, US, vol. 21, no. 12, 1 December 2011 (2011-12-01), pages 685 - 687, XP011379862, ISSN: 1531-1309, DOI: 10.1109/LMWC.2011.2170669 * |
Also Published As
Publication number | Publication date |
---|---|
EP2997655A1 (en) | 2016-03-23 |
US9659120B2 (en) | 2017-05-23 |
US9727677B2 (en) | 2017-08-08 |
US20160117430A1 (en) | 2016-04-28 |
WO2014184774A1 (en) | 2014-11-20 |
CN105393452B (en) | 2019-04-02 |
US20160079933A1 (en) | 2016-03-17 |
EP2997656A1 (en) | 2016-03-23 |
CN105393452A (en) | 2016-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9727677B2 (en) | Baseband equivalent volterra series for behavioral modeling and digital predistortion of wideband transmitters | |
US20230370937A1 (en) | Method and system for baseband predistortion linearization in multi-channel wideband communication systems | |
EP3111554B1 (en) | Dual-band signal transmission | |
EP3005554B1 (en) | Digital predistortion of wideband power amplifiers with reduced observation bandwidth | |
US9432228B2 (en) | Digital pre-distortion filter system and method | |
US8989307B2 (en) | Power amplifier system including a composite digital predistorter | |
Fehri et al. | Baseband equivalent Volterra series for behavioral modeling and digital predistortion of power amplifiers driven with wideband carrier aggregated signals | |
KR101669173B1 (en) | Baseband-derived wideband rf digital predistortion system | |
US20110235748A1 (en) | Active antenna array having analogue transmitter linearisation and a method for predistortion of radio signals | |
Suryasarman et al. | A comparative analysis of adaptive digital predistortion algorithms for multiple antenna transmitters | |
US20110235749A1 (en) | Active antenna array having analogue transmitter linearisation and a method for predistortion of radio signals | |
CN106506417A (en) | A kind of digital pre-distortion system of narrowband feedback and method | |
Younes et al. | Three-dimensional digital predistorter for concurrent tri-band power amplifier linearization | |
CN105634415B (en) | Digital pre-distortion system and method for amplified signal | |
CN104639481B (en) | A kind of multi-band signal processing method and equipment | |
Jaraut et al. | Harmonically related concurrent tri-band behavioral modeling and digital predistortion | |
Tafuri et al. | Memory models for behavioral modeling and digital predistortion of envelope tracking power amplifiers | |
Bensmida et al. | Overlapped segment piece-wise polynomial pre-distortion for the linearisation of power amplifiers in the presence of high PAPR OFDM signals | |
US20140308908A1 (en) | Apparatus and method for signal predistortion | |
Wang et al. | Performance analysis of multi-stage cascaded digital predistortion | |
Naraharisetti et al. | 2D quasi exact inverse of PA model in digital predistorter for concurrent dual-band system | |
Changwei et al. | One-step model extraction method for joint polynomial/interpolating lookup-table two-box nonlinear model of power amplifier | |
Cabarkapa et al. | 2-D nonlinearity compensation technique for concurrent dual-band wireless transmitters | |
Xia et al. | Dual-band linear filter assisted envelope memory polynomial for linearizing multi-band power amplifiers | |
Reddy | Digital pre-distortion using Legendre polynomials for band 25 transmitters |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 14732395 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14888876 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2014732395 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |