I. INTRODUCTION
Modern telecommunication standards make use of multilevel and/or multicarrier modulation schemes in order to accommodate high information rates in a limited bandwidth and to provide immunity to channel distortion and multipath effects. Spectral efficiency and immunity to channel conditions, however, are achieved at the expense of a high peak-to-average ratio (PAPR) transmitted signal and of a high sensitivity to intermodulation distortion (IMD). IMD is normally generated by nonlinear dynamic effects in the power amplifier (PA) at the end of the transmission chain. Moreover, the modulation process and PA nonlinearity generate out-of-band distortion, which must remain inside the mask specified by the communication standard. If no linearization scheme is adopted at the transmitter, the only way of meeting the in-band and the out-of-band distortion requirements for signals with high PAPRs is to operate the PA at a high level of backoff, which normally results in poor power efficiency. On the contrary, if the proper linearization scheme is used, significant improvements in system performance and costs can be obtained. In this paper, we address the linearization of a millimeter (mm)-wave PA using digital baseband predistortion [Reference Liu, Boumaiza and Ghannouchi1–Reference Melczarsky, Gilabert, Montoro, Bertran, Di Giacomo and Filicori3]. Possible applications include current mm-wave high-bandwidth applications such as the bidirectional satellite service at Ka band, IEEE 802.16-SC, LMDS, Point-to-Point Microwave Service and Private Operational Fixed Service (i.e., TV Broadcast, 38.6–40.0 GHz), HAPS, Amateur Radio (licensed), and defense applications. To this aim, and to preliminarily test the validity of the approach, we employed a circuit-level model of a PA operating at 38 GHz, which is based on a simple yet accurate non-quasi-static model of a GaAs pHEMT device. Computer simulations were then performed in order to identify a nonlinear auto-regressive moving average (NARMA) behavioral model for the PA from which the predistorter was derived. The predistortion function, also of NARMA type, was expressed as a combination of basic predistortion cells (BPCs) in order to be efficiently implemented in an FPGA, which is normally already present in most digital transmitters. The performance of the proposed solution was then evaluated for different single- and multi-carrier modulation schemes using computer simulations.
The rest of the paper is organized as follows: section II presents a review of linearization techniques, section III describes the device model and the main considerations made in the design of the mm-wave PA chosen for linearization, section IV presents the NARMA behavioral model used to represent the PA in the predistortion solution, section V describes the predistortion solution, and section VI presents validation results obtained from computer simulations with different modulated test signals.
II. PA LINEARIZATION
A number of linearization solutions have been proposed in the literature [Reference Kenington4]. However, an optimal linearization technique does not exist and so the choice of the linearization scheme is strongly dependent on the application and on a number of factors such as the linearity requirements, the complexity added by the linearizer (and the trade-off between accuracy and complexity), the signal bandwidth, the technology used (e.g., analog or digital), the power efficiency requirements, the robustness against system changes, etc. Feedback-type linearizers (Envelope, Cartesian, Polar) have been extensively used in the past and are still used in specific applications (e.g., circuit level linearization). However, stability problems make them unsuitable for current high-bandwidth applications. In addition, feedback linearizers reduce overall system gain, which can affect system efficiency. On the other hand, they are less sensitive to disturbances or changes in the transmitter subsystem and therefore good IMD cancellation can be achieved for low or moderate signal bandwidths. Different from feedback linearizers, Envelope Elimination and Restoration and the outphasing techniques (LINC, CALLUM) have revived thanks to the use of current high-speed digital signal processors (DSPs), which have permitted one to overcome certain drawbacks associated with analog circuitry. One of the biggest challenges of these techniques is with respect to power efficiency, since power losses associated with the combination of the signal components can in some cases outweigh the benefits of employing highly efficient PAs. On the other hand, the processing speed of current DSPs limits the linearization bandwidth, although this may change in the future as processing speed improves. Feedforward together with digital predistortion linearization are the most suitable linearization techniques to cope with wideband signals. Feedforward linearization is unconditionally stable, can achieve effective cancelation of IMD, and is immune to PA memory effects. However, its open-loop nature makes it sensitive to possible unbalances in cancelation loops, effectively degrading the overall linearity performance whenever imperfect cancelation occurs. This can be compensated by means of a proper loop control, but it requires additional circuitry. Moreover, the error amplifier and the losses associated with the output combiner strongly impact on the IMD cancelation performance and on overall power efficiency. Analog predistortion is a simpler solution with respect to digital predistortion, but although it can cope with high bandwidths, it only gives moderate IMD suppression. Moreover, analog solutions operating at RF are very sensitive to the electrical lengths in the circuit, which become critical for mm-wave applications. A similar dependence on implementation maladjustments appears in LINC transmitters, due to the need for a good equilibrium between the PA in the two branches. Digital predistortion performed at baseband can achieve better correction capabilities (it is independent of the final carrier frequency of operation) than analog solutions and can be adaptive, thus allowing for corrections of possible unexpected unbalances. The main drawbacks are related to power consumption associated to the DSP device and bandwidth limitations due to the speed of the DSP. Since DSPs are currently present in most communications transmitters for codification and digital modulation purposes, digital baseband predistortion can be a good linearization solution for a high-bandwidth mm-wave application such as the one we are considering, especially when memory effects in the PA are of significant importance in IMD cancelation performance.
III. PA DESIGN
To avoid the uncertainties associated with high-frequency prototyping and measurements, and thus to be able to concentrate on the performance of the PA behavioral model and predistortion scheme, it was decided to preliminarily validate the linearization approach using an accurate circuit-level description of the PA. Accordingly, the main goal in the design of the PA was to obtain a realistic mm-wave PA which was, at the same time, computationally efficient even with complex modulated signals.
This section describes the main characteristics of the PA and the considerations made during the design. The first part presents the non-quasi-static device model that was employed, whereas the second part illustrates the main features of the PA.
A) Non-quasi-static device modeling
To obtain a realistic circuit-level model of a high-frequency PA, an accurate device model is needed that is capable of accounting for short-memory (i.e., non-quasi-static) effects occurring at high frequencies [Reference Santarelli5]. Moreover, for non-constant envelope applications, low-frequency dispersive effects need to be accounted for. Besides accuracy in predicting the device behavior, the numerical efficiency of the model is also important, as it will be used with complex modulated signals and simulated in the mixed, time-frequency, envelope domain. This section briefly describes the nonquasi-static equivalent-voltage model (NQS-EVM) [Reference Santarelli5] of a FET device, which was used to account for nonlinear short memory effects at the intrinsic active device of the PA. Long-term memory effects were accounted for by the empirical model described in [Reference Filicori, Vannini, Santarelli, Mediavilla, Tazon and Newport6].
The NQS-EVM model is an empirical electron device model featuring a reliable identification procedure and a simple CAD implementation. It also provides good prediction capabilities up to very high frequencies under mildly nonlinear operation. Although the model is technology independent and can be implemented using either look-up tables (LUTs) or equivalent circuits, only the LUT-based version for FET devices will be considered in the following.
The basic idea of the model is to use an associated QS device controlled by suitably defined equivalent voltages: in such a way, NQS effects arising at high frequencies are modeled by means of voltage perturbations imposed by the equivalent voltage sources.
In QS conditions, the behavior of an intrinsic FET device can be simply described by the well-known voltage-controlled equations:
![{\bf i} \lpar t\rpar = {\bf F} \lsqb {\bf v} \lpar t\rpar \rsqb + {\bf C} \lsqb {\bf v} \lpar t\rpar \rsqb {\hbox{d} {\bf v} \lpar t\rpar \over \hbox{d}t}\comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn1.gif?pub-status=live)
where i and v are the vectors of the port currents and voltages, respectively, and F and C are purely algebraic functions. At high frequencies, NQS phenomena, which are ignored in (1), become important. We define suitable equivalent voltages c and
d:
![{\hat {\bf v}}_c \lpar t\rpar \ {\dot =} \ {\bf v} \lpar t\rpar + {\Delta {\bf v}}_c \lpar t\rpar,](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn2.gif?pub-status=live)
![{\bf \hat v}_d \lpar t\rpar \ {\dot =} \ {\bf v} \lpar t\rpar + { \Delta {\bf v}}_d \lpar t\rpar \comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn3.gif?pub-status=live)
in order to provide the proper NQS currents in a virtually QS device described by (1) when the actual controlling voltages are replaced with the equivalent voltages:
![{\bf i} \lpar t\rpar = {\bf F} \lcub {\bf \hat v}_c \lpar t\rpar \rcub + {\bf C} \lcub {\bf \hat v}_d \lpar t\rpar \rcub {\hbox{d} {\bf \hat v}_d \lpar t\rpar \over \hbox{d}t}.](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn4.gif?pub-status=live)
In particular, the equivalent voltage domain can be obtained according to (1) by means of the voltage perturbations Δvd and Δvc, which take into account all the fast dynamics of the NQS phenomena.
General considerations and empirical evidence described in [Reference Santarelli5] suggest that a linear dependence can be assumed for Δvd and Δvc, that the two voltage perturbations are not mutually independent, and that in the case of a strictly intrinsic GaAs pHEMT the voltage perturbations appear only in series to the gate. In such a way, just two voltage corrections, linearly voltage dependent, account for all NQS effects, according to Fig. 1.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629050858-35599-mediumThumb-S1759078709000075_fig1g.jpg?pub-status=live)
Fig. 1. Topology of the NQS-EVM model.
Provided that the duration of memory effects is short with respect to the period of the operating signals, a “modified” dynamic Volterra-series approach [Reference Mirri, Iuculano, Filicori, Pasini, Vannini and Pellegrini7] can be adopted to model the voltage perturbations, which, by discretizing the memory time in N D time slots Δτ, can be written by means of linear discrete convolutions:
![\eqalignno{\Delta v_{GS\comma d} \lpar t\rpar &= \sum_{p = 1}^{N_D} a_{1p} \lsqb v_{GS} \lpar t - p \Delta \tau\rpar - v_{GS} \lpar t\rpar \rsqb \cr &\quad + \sum_{p = 1}^{N_D } a_{2p} \lsqb v_{DS} \lpar t - p \Delta \tau\rpar - v_{DS} \lpar t\rpar \rsqb \comma}](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn5.gif?pub-status=live)
![\hskip -8pt \Delta v_{GS\comma c} \lpar t\rpar = \rho \Delta v_{GS\comma d} \lpar t\rpar .](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn6.gif?pub-status=live)
The dependence on the dynamic voltage deviations v(t −τ) − v(t) allows for a vanishing contribution of (5) and (6) when the signal frequencies involved are not high enough for NQS effects to manifest.
B) Model identification and validation
The first part of model identification involved the derivation of the “strictly intrinsic” device measurements. These were obtained by de-embedding the measured S parameters at the extrinsic ports both from the extrinsic parasitic elements and from all the linear capacitances in parallel to the intrinsic device. Parasitic elements were, in this case, identified through standard cold-FET procedures in pinched-off conditions, and also was the purely capacitive parallel network. The strictly intrinsic measurements were then used in order to identify the quasi-static device and the non-quasi-static voltage corrections.
The chosen device was a 0.25 µm GaAs pHEMT having 10 fingers and 600 µm of total gate periphery. In this case, a LUT-based definition was chosen for the associated quasi-static device both for the conductive currents and for the displacement currents. The conductive currents include a low-frequency dispersive term that accounts for self-heating and surface trap level states, and is modeled using an empirical approach proposed in [Reference Filicori, Vannini, Santarelli, Mediavilla, Tazon and Newport6]. The quasi-static capacitive currents are modeled through nonlinear LUT-based capacitors that are identified from the imaginary part of the strictly intrinsic admittance parameters, in a frequency range where reactive effects are evident, but non-quasi-static phenomena are still negligible.
The second part of model identification concerns the identification of parameters a 1p, a 2p, and ρ, defining the non-quasi-static voltage corrections. This was achieved by resolving a simple overdetermined system of equations based on
![\eqalignno{{\bf Y}\lpar {\bf V}\comma \; \omega \rpar &= {\bf G}\lpar \hat{\bf V}_c \rpar \left[{\bf 1} + \rho \sum_{p = 1}^{N_D} {\bf a}_p \lpar \hbox{e}^{ - j\omega p\Delta \tau } - 1\rpar \right]\cr &\quad + j\omega {\bf C}\lpar \hat{\bf V}_d\rpar \left[{\bf 1} + \sum_{p = 1}^{N_D} {\bf a}_p \lpar \hbox{e}^{ - j\omega p \Delta \tau} - 1\rpar \right]\comma}](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn7.gif?pub-status=live)
which gives the strictly intrinsic bias-dependent admittance parameters in the frequency domain and which is derived by inspection of Fig. 1 and simple algebraic manipulations. For the device considered, a total memory length of 1.8 ps divided into three elementary delays gave good results. Equation 7 was evaluated over a reduced bias grid (V GS: −1.2 to −0.2 V; V DS: 2–10 V), so as to have maximum accuracy over the region defined by the dynamic load curve of the device in large-signal operation of the PA. Moreover, in order to properly identify the non-quasi-static corrections, (7) was evaluated in a frequency range covering both quasi-static and non-quasi-static phenomena (3.5 GHz < f < 76 GHz).
Figure 2 shows some small-signal validation results that refer to the intrinsic admittance parameters Y 11 and Y 21 up to twice the fundamental frequency of the PA (76 GHz) with the device biased at the nominal quiescent point of the PA (V GS: −0.55 V; V DS: 6.5 V). The figures show the measured and modeled parameters, as well as the results that would have been obtained if a conventional quasi-static model had been employed. As can be seen from the figures, for this bias condition where the device is very conductive, the behavior at high frequencies is far from quasi-static and is accurately accounted for by the model thanks to the voltage-controlled dynamic deviation generators. Similar results are obtained for Y 12 and Y 22, which however are less affected by non-quasi-static phenomena. On the contrary, for bias conditions where the device channel is pinched-off, an almost quasi-static behavior is observed (not shown here due to lack of space) and the model is still able to predict it owing to the fact that it was identified at the strictly intrinsic ports.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629051044-11154-mediumThumb-S1759078709000075_fig2g.jpg?pub-status=live)
Fig. 2. Measured and modeled intrinsic Y parameters (V GS, V DS) = (−0.55 V, 6.5 V) exhibiting non-quasi-static behavior. Results obtained with a quasi-static model are included for comparison.
Figure 3, instead, shows, for the fundamental frequency of the PA (38 GHz), the measured and modeled admittance parameters as a function of the bias voltages in the region where the model was identified. As can be seen from the figures, even if bias-independent coefficients were assumed for the dynamic deviation generators, the model is able to predict the measured parameters inside a relatively large area of operation for the device including off- as well as on-state conditions.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629051329-61278-mediumThumb-S1759078709000075_fig3g.jpg?pub-status=live)
Fig. 3. Measured and modeled intrinsic Y parameters at 38 GHz for V DS = 2, 6.5, and 10 V.
C) Model implementation
After the identification procedure, the model was implemented in an Agilent ADS CAD environment, as shown in Fig. 4 for the intrinsic device. The conductive currents, including low-frequency dispersive effects due to traps and self-heating, are modeled using LUTs and the empirical low-frequency model described in [Reference Filicori, Vannini, Santarelli, Mediavilla, Tazon and Newport6].
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629051657-38347-mediumThumb-S1759078709000075_fig4g.jpg?pub-status=live)
Fig. 4. CAD implementation of the non-quasi-static device model at the intrinsic ports including displacement (top) and conductive low-frequency-dispersive currents (bottom). Non-quasi-static effects are accounted for by dynamic-deviation-based VCVSs. A purely capacitive linear network separates the “strictly intrinsic” ports from the intrinsic ones.
The displacement currents of the quasi-static device are modeled with nonlinear LUT-based capacitors and implemented using a symbolically defined device. The purely dynamic voltage corrections applied to both conduction and displacement parts of the quasi-static device are implemented through standard linear VCVSs, which depend on intrinsic port voltages and on model parameters a 1p, a 2p, ρ, Δτ, and N D, which were previously identified. A linear parallel capacitive network separates the strictly intrinsic from the intrinsic device ports, and the standard lumped parasitic network depicted in Fig. 5 is embedded in order to obtain the extrinsic device ports.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629051819-01502-mediumThumb-S1759078709000075_fig5g.jpg?pub-status=live)
Fig. 5. Parasitic linear network embedded onto the intrinsic device to obtain the extrinsic ports.
D) Power Amplifier
This section describes the design of an integrated mm-wave PA using the non-quasi-static device model identified and implemented in the previous sections. In keeping with the objective of obtaining an efficient yet realistic circuit-level model of the amplifier, a single-stage single-cell PA was considered, which can represent the typical output stage of a real mm-wave PA. Due to linearity requirements, class-AB operation was selected, and the device was biased at 35% of IDSS. The load and source impedance were optimized taking into account IMD and gain, using a two-tone source- and load-pull measurement setup and were also validated with simulations using the non-quasi-static device model [Reference Santarelli5]. Matching networks were synthesized using distributed networks, as would be the case in any integrated PA operating at such a frequency. Figure 6 shows a schematic of the designed PA and Table 1 reports its main performances.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629052036-56815-mediumThumb-S1759078709000075_fig6g.jpg?pub-status=live)
Fig. 6. Schematic of the designed mm-wave PA.
Table 1. Main performances of the power amplifier.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_tab1.gif?pub-status=live)
High-frequency PAs typically exhibit both short- and long-term memory effects due to electrical and thermal phenomena [Reference Vuolevi8], which generally make their linearization more difficult. Any realistic PA model, therefore, needs to account for these effects as they will have an impact in the design and performance of the predistorter. In the application under consideration, however, and taking into account that the signal bandwidth even for very high bit rate communication standards is in the order of tens of MHz, and that therefore the relative bandwidth is small with respect to an mm-wave carrier, only long-termFootnote 1 memory effects become important.Footnote 2 This is so because the response of the matching networks and the device are substantially flat over such a small relative bandwidth and therefore they do not exhibit any dynamic effects. Long-term memory effects, instead, are generally associated with thermal phenomena, surface charge traps in the active devices, and low-frequency dispersive bias circuits. The first two were modeled by the low-frequency part of the NQS model whereas bias networks effects were accounted for by careful design of the bias circuit. In fact, if ideal bias networks are considered, long-term memory effects are hardly noticeable, and, as shown in Fig. 7, the two-tone IM3 term is almost independent of tone spacing (for tone spacings up to 20 MHz). The same figure also shows that when a more realistic bias system is considered, not only the IM3 terms depend on tone spacing but also some asymmetry between the left and right IM3 terms appear, as is typically observed in real PAs. This variation dependence over such a narrow frequency range is a symptom of long-term memory effects.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629052158-18081-mediumThumb-S1759078709000075_fig7g.jpg?pub-status=live)
Fig. 7. PA two-tone IMD tests for different input powers and bias systems.
IV. BEHAVIORAL MODELING
The PA is a nonlinear dynamic system whose memory effects can be of significant importance when wideband signals are considered. The digital predistortion (DPD) linearization technique can compensate for both nonlinear and dynamic PA unwanted effects. The DPD technique requires nonlinear dynamic behavioral models capable of reproducing the PA nonlinear behavior as well as its memory effects. In the last few years, a significant number of behavioral models for PAs have been proposed, as discussed in [Reference Pedro and Maas9]. However, to be suitable for DPD purposes behavioral models not only need to be accurate enough to mimic PA behavior but also have to be invertible and computationally efficient.
For our particular application, a NARMA model, described in [Reference Gilabert, Montoro and Cesari10], was chosen as a good compromise between accuracy and computational efficiency. The NARMA model can be seen as an extension of a simpler nonlinear moving average (NMA) model that incorporates a nonlinear auto-regressive block. The advantage of using a NARMA model is the introduction of a nonlinear feedback path (infinite impulse response – IIR – terms) that may allow one to relax the number of delayed samples considered to model the PA, in comparison with a model using only finite impulse response (FIR) terms.
One of the main weaknesses of the NARMA model regards its stability. The use of nonlinear feedback paths can derive in the overall system instability. Fortunately, the stability of a NARMA model can be studied and guaranteed using a technique based on the small-signal gain theory, which is advanced in [Reference Gilabert, Montoro and Cesari10]. Figure 8 depicts the block diagram that describes the NARMA model.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629052416-99983-mediumThumb-S1759078709000075_fig8g.jpg?pub-status=live)
Fig. 8. NARMA model structure.
The general input–output mathematical expression in a NARMA model can be described as
![\hat y_{\rm NARMA} \left(k \right)=\sum_{i=0}^N {f_i \left({x\left({k - \tau _i } \right)} \right) - \sum_{j=1}^D {g_j \left({\hat y\left({k - \tau _j } \right)} \right)} }\comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn8.gif?pub-status=live)
with f i (·) and g j (·) being nonlinear memoryless functions, and where τ i (τ 0 = 0) and τ j (τ ⊂ ℕ) are the most significant sparse delays of the input and output, respectively, contributing to the description of the PA memory effects. Further details on how to obtain the best sparse delays to characterize the PA dynamics can be found in [Reference Gilabert, Silveira, Montoro, Gadringer and Bertran11, Reference Gilabert, Montoro and Cesari12].
By considering memoryless nonlinear functions f i (·) and g j (·) characterized by polynomials of order P,
![f_i \lpar x \lpar k - \tau_i\rpar \rpar = \sum_{p = 0}^P \alpha_{pi} x\lpar k - \tau_i\rpar \vert x \lpar k - \tau_i\rpar \vert^p\comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn9.gif?pub-status=live)
![g_j \lpar \hat{y} \lpar k - \tau_j\rpar \rpar = \sum_{p = 0}^P \beta_{pj} \hat{y} \lpar k - \tau_j\rpar \vert \hat{y} \lpar k - \tau_j\rpar \vert^p\comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn10.gif?pub-status=live)
it is possible to rewrite (7) in a more compact matrix notation:
![\hat{y}_{\rm NARMA}\lpar k\rpar = \hat{y} \lpar k\rpar = {\bf \hat \delta}^H {\bf \Phi}\comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn11.gif?pub-status=live)
where H denotes Hermitian (conjugate transpose) and where the vector of parameters and the data matrix Φ are defined in (11).
![{\bf \hat \delta} = \left[\matrix{\alpha_{00} \hfill \cr \vdots \hfill \cr \alpha _{P0} \hfill \cr \vdots \hfill \cr \alpha _{0\, N} \hfill \cr \vdots \hfill \cr \alpha _{PN} \hfill \cr - \beta _{01} \hfill \cr \vdots \hfill \cr - \beta _{P1} \hfill \cr \vdots \hfill \cr - \beta _{0D} \hfill \cr \vdots \hfill \cr - \beta _{PD}} \right]\comma \; \quad {\bf \Phi} = \left[\matrix{x\lpar k\rpar \hfill \cr \vdots \hfill \cr x\lpar k\rpar \vert x\lpar k\rpar \vert^P \hfill \cr \vdots \hfill \cr x\lpar k - \tau_N\rpar \hfill \cr \vdots \hfill \cr x\lpar k - \tau _N \rpar \vert x\lpar k - \tau_N\rpar \vert^P \hfill \cr \hat{y} \lpar k - \tau_1\rpar \hfill \cr \vdots \hfill \cr \hat{y} \lpar k - \tau_1 \rpar \vert \hat{y} \lpar k - \tau_1\rpar \vert^P \hfill \cr \vdots \hfill \cr \hat{y} \lpar k - \tau_N\rpar \hfill \cr \vdots \hfill \cr \hat{y} \lpar k - \tau_N \rpar \vert \hat{y} \lpar k - \tau_N\rpar \vert^P} \right].](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn12.gif?pub-status=live)
In order to extract the vector of parameters , we performed a least-squares (LS) identification from input–output (x PA, y PA) observations of the PA. The identification error is defined as
![e\lpar k\rpar = y_{PA} \lpar k\rpar - \hat{y} \lpar k\rpar \comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn13.gif?pub-status=live)
and the LS solution for is therefore
![{\bf \hat \delta} = \lpar {\bf \Phi}^H {\bf \Phi}\rpar ^{ - 1} {\bf \Phi}^H {\bf y}_{{\bf PA}}\comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn14.gif?pub-status=live)
with yPA = [y PA(0), y PA(1), … , y PA(L) ]T being the vector of output observations of length L.
The normalized mean square error achieved using the NARMA model, just considering one delayed sample of the input and another delayed sample of the output in the NARMA structure (small influence of memory effects), was about −43 dB.
V. DIGITAL PREDISTORTION LINEARIZATION
Figure 9 presents a block diagram of the digital baseband predistortion scheme that was adopted. To derive the predistortion function from the original NARMA behavioral model, it is necessary to rewrite the equations to obtain x A(k) in terms of the desired linear output y D(k). A detailed explanation on this procedure can be found in [Reference Montoro, Gilabert, Bertran, Cesari and Silveira13].
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629052657-10235-mediumThumb-S1759078709000075_fig9g.jpg?pub-status=live)
Fig. 9. Predistortion block diagram.
To implement the DPD in an FPGA, it is necessary to express the predistortion function as a combination of fundamental building blocks called BPCs. A BPC requires simple hardware blocks: a complex multiplier, a dual port RAM memory block acting as a LUT, an address calculator, and two ports of control. Further details on the implementation in an FPGA device can be found in [Reference Gilabert, Cesari, Montoro, Bertran and Dilhac14]. The structure of the predistortion function, in terms of BPCs, is depicted in Fig. 10.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629045815-61396-mediumThumb-S1759078709000075_fig10g.jpg?pub-status=live)
Fig. 10. Predistortion function.
The predistortion function can be expressed in terms of combinations of BPCs as
![x_A \lpar k\rpar = z\lpar k\rpar G_{{\rm LUT} \hat{f}_0^{ - 1}} \lpar \vert z\lpar k\rpar \vert\rpar \comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn15.gif?pub-status=live)
where z(k) is defined as
![\eqalignno{z\lpar k\rpar &= y_D \lpar k\rpar + \sum_{j = 1}^D {\underbrace{y_D \lpar k - \tau_j\rpar G_{\rm LUT} \hat{g}_j \lpar \vert y_D \lpar k - \tau_j\rpar \vert\rpar }_{\hat{g}_j \lpar y_D \lpar k - \tau_j\rpar \rpar }} \cr &\quad - \sum_{i = 1}^N {\underbrace{x_A \lpar k - \tau_i\rpar G_{{\rm LUT} \hat{f}_i} \lpar \vert x_A \lpar k - \tau_i\rpar \vert\rpar }_{\hat{f}_i \lpar x_A \lpar k - \tau_i\rpar \rpar}},}](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn16.gif?pub-status=live)
where G LUT (for both ĝ j and i) are complex gains stored in their corresponding LUT, x A(k) is the output of the DPD, and y D(k) is the desired output defined as the signal to be transmitted (x T(k)) multiplied by a linear amplification (G linear),
![y_D \lpar k\rpar = x_T \lpar k\rpar G_{\rm linear}.](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151022065804067-0736:S1759078709000075_eqn17.gif?pub-status=live)
Then, the adaptive process followed by the digital adaptive DPD in order to perform the PA linearization consists in the following four steps: (1) to identify new nonlinear functions i and ĝ j by monitoring current PA input and output samples, (2) to test the stability of both PA and DPD NARMA models using the method proposed in [Reference Gilabert, Montoro and Cesari10], (3) to invert the nonlinear function
0 to find the DPD output, and (4) to generate all necessary LUT contents to implement the predistorter function; and then back to step 1. The identification of the PA behavioral model is performed obtaining the LS solution presented in (Reference Gilabert, Cesari, Montoro, Bertran and Dilhac14). For this particular implementation, the LS takes complex data buffers of 2048 samples and a forgetting factor is adopted in order to prevent abrupt changes in the parameters that define the PA model.
In the case of the designed PA, good results were obtained with a predistorter function having only three BPCs to compensate for the nonlinear distortion and memory effects: in particular: the static BPC (associated to f0−1), one BPC-FIR (associated to g1), and one BPC-IIR (associated to f1).
VI. VALIDATION RESULTS
In a first approach toward a final implementation, it is necessary to analyze the linearization capabilities shown by this NARMA-based DPD. The DPD linearizer has to be capable of compensating both in-band and out-of-band nonlinear distortion. To validate the functioning of the adaptive predistorter, an FPGA DPD simulator has been implemented in Matlab. This simulator emulates the DPD function as it is performed in an FPGA, that is, by means of a set of BPCs. Moreover, a NARMA behavioral model has been used to reproduce the nonlinear dynamics of the designed PA operating at 38 GHz. Both single carrier (SC) and OFDM 16-QAM test signals (roll-off = 0.22) presenting PAPRs around 9–10 dB have been considered.
To ensure a fair comparison, it is necessary to take into account unlinearized backed-off signals presenting similar mean output power to the linearized output signal. Figure 11 shows the AM–AM characteristics for: unlinearized PA, unlinearized PA with backoff operation, inverse PA characteristic, and the linearized PA. It is possible to observe that the DPD is capable of counteracting nonlinear distortion and get into compression in a linear way. Moreover, Fig. 12 and Fig. 14 shows the in-band distortion compensation achieved by the DPD, measured in terms of EVM. The 16-QAM constellation of the backed-off signal presents an EVM of around 11% (SC) and 37% (OFDM), respectively, while thanks to predistortion linearization it can be reduced up to less than 1%.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629050015-46267-mediumThumb-S1759078709000075_fig11g.jpg?pub-status=live)
Fig. 11. PA AM–AM characteristics: unlinearized PA, unlinearized PA with back-off operation, inverse PA characteristic, and linearized PA.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629050107-78805-mediumThumb-S1759078709000075_fig12g.jpg?pub-status=live)
Fig. 12. SC 16-QAM constellation for both linearized (black) and backed-off output signals (gray).
Finally, Figs 13 and 15 show the spectral regrowth reduction (out-of-band distortion compensation) of ~6 dB achieved by the DPD in comparison to an unlinearized backed-off signal for both SC and OFDM configurations (Figs 13, 15).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629050344-42891-mediumThumb-S1759078709000075_fig13g.jpg?pub-status=live)
Fig. 13. Output power spectra (SC 16-QAM) for linearized PA, unlinearized PA, and unlinearized PA with back-off operation.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629050524-69615-mediumThumb-S1759078709000075_fig14g.jpg?pub-status=live)
Fig. 14. OFDM 16-QAM constellation for both linearized (black) and backed-off output signals (gray).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629050625-54306-mediumThumb-S1759078709000075_fig15g.jpg?pub-status=live)
Fig. 15. Output power spectra (OFDM 16-QAM) for linearized PA, unlinearized PA, and unlinearized PA with back-off operation.
VII. CONCLUSION
This paper has addressed the linearization of an mm-wave PA by means of an adaptive digital predistorter operating at baseband. This DPD is based in a NARMA structure that has been used for both characterizing and compensating PA nonlinearities and memory effects. Simulation results have evidenced that by means of the adaptive DPD it is possible to counteract the PA nonlinear behavior and thus to compensate both in-band and out-of-band distortion.
ACKNOWLEDGEMENT
This work was partly supported by the former EU network TARGET (IST-1–507893-NOE), by the Italian Ministry of University and Research (MIUR), and by the Spanish Government under project TEC2005-07985-C03–02.
Ilan Melczarsky received an electronics engineering degree (cum laude) from the University of Mar del Plata, Buenos Aires, Argentina, in 2000, and a Ph.D. degree in electronics and computer science from the University of Bologna, Italy, in 2005.
From 2000 to 2001, he has been a Senior Product Engineer with Ericsson Argentina. From 2002 to 2005, he was a Research Fellow with the Italian Research Council (CNR). Since 2005, he has been with the Dipartimento di Ingegneria Elettronica (DEIS), University of Bologna, as a Post-Doctoral Researcher and Teaching Assistant. His research interests include nonlinear characterization and modeling of high-frequency electron devices and circuits, electronic measurements, and high-frequency circuit design.
Pere Ll. Gilabert Pinal was born in Vilanova i la Geltrú, Barcelona (Spain) in 1978. He received the degree in telecommunication engineering from the Technical University of Catalonia (UPC) in 2002, and he developed his master's thesis at the University of Rome “La Sapienza” with an exchange grant. He joined the Department of Signal Theory and Communications (TSC) in 2003 and received his Ph.D. from the UPC in February 2008. His research is focused on system level linearization techniques for RF power amplifiers. Other research interests include digital signal processing for communications, communication electronics, and circuit theory. He is currently lecturing on undergraduate courses at the Castelldefels School of Technology (EPSC) as Assistant Professor.
Valeria Di Giacomo was born in Salerno, Italy. She received the Laurea degree in electronic engineering from the University of Bologna, Italy, in 2005.
Since then, she has collaborated with the Department of Electronics of the same University, and in 2006 she joined the Engineering Department of the University of Ferrara as a Ph.D. student. Her research activity is mainly oriented to nonlinear electron device modeling for microwave applications.
Eduard Bertran received the Engineer (option communication electronics) and Doctor Engineer degrees in telecommunication engineering, both from the Universitat Politècnica de Catalunya (UPC), Barcelona, Spain, in 1979 and 1985, respectively. He joined the Department of Signal Theory and Communications (TSC) in 1987, where he is currently a full professor. He has been the head of studies of the TSC department and an Associate Dean in different telecommunications schools. His research interests include control, communication electronics, signal processing, and circuit theory, topics where he has produced different books and peer-reviewed journal and conference papers. He is a member of IFAC and a senior member of IEEE. He has collaborated in different national and European research projects and teams.
Fabio Filicori received the Dr. Ing. degree in electronic engineering in 1974 from the University of Bologna, where he joined the Faculty of Engineering as Assistant Researcher and, later, as Associate Professor. In 1990 he became Full Professor of Electronics at the University of Perugia. In 1991 he moved to the University of Ferrara, where he was coordinator of the degree course in electronic engineering. Currently, he is Full Professor at the University of Bologna, where he has been coordinator of the Ph.D. course in electronics, computer science, and telecommunications. He has been the coordinator of research projects in electronic engineering promoted by the Ministry of University and Research. In 2007 he was appointed as member of the Technology Commission of the Italian Space Agency. He has been workpackage leader for the European NoE TARGET. He has been TPC chairman for the EUMIC Conference and editorial board member for IEEE MTT Transactions. He is the author of about 200 papers concerning his research in nonlinear microwave circuit design, electron device modeling, electronic measurements, and industrial electronics.