Next Article in Journal
Cooperative Coverage Path Planning for Multi-Mobile Robots Based on Improved K-Means Clustering and Deep Reinforcement Learning
Previous Article in Journal
MgMViT: Multi-Granularity and Multi-Scale Vision Transformer for Efficient Action Recognition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Linear Simulation by Harmonic Balance Techniques of Load Modulated Power Amplifier Driven by Random Modulated Signals

XLIM Laboratory, UMR CNRS n°7252, University of Limoges, Faculté des Sciences et Techniques, Campus La Borie, 123 avenue Albert Thomas, 87060 Limoges CEDEX, France
*
Author to whom correspondence should be addressed.
Electronics 2024, 13(5), 947; https://doi.org/10.3390/electronics13050947
Submission received: 30 January 2024 / Revised: 23 February 2024 / Accepted: 25 February 2024 / Published: 29 February 2024
(This article belongs to the Section Microwave and Wireless Communications)

Abstract

:
The simulation of the steady state and the non-linear stability of a load modulated power amplifier (LMPA) driven by a random modulated generator, fully performed in the frequency domain by harmonic balance (HB) techniques, is presented. The non-linear microwave circuit and the driving pseudo-random modulated (PRM) generator are integrally defined in the frequency domain. The simulation is implemented and performed using commercially available circuit simulation software. The demodulation of the output signal of the LMPA is implemented with optimally matched filters as software-defined demodulation. The simulated dynamic results of a Quasi-MMIC GaN Doherty power amplifier (DPA) are shown and compared to the measured results with a 16-QAM driving signal at 10 MS/s. The time-domain measurement allows the validation of the new simulation technique through the comparison of both the measured and the simulated error vector magnitude (EVM), the left and right adjacent channel power ratios (ACPRs) versus the average output power. This new simulation is then called pseudo-random modulated harmonic balance (PRM-HB) simulation. The full PRM-HB simulation of an LMPA driven by a random modulated signal, performed in the frequency domain at the design circuit level, results in an advanced simulation tool in the frame of the design of RF circuits and subsystems for telecommunication applications.

1. Introduction

Today, with the new telecommunication standards (5G and beyond), the modulations with a large peak-to-average power ratio (PAPR) become conventional. Therefore, it is mandatory that circuit designers have an accurate and reliable software design tool allowing, in the same HB simulation frame, the simulation of RF components or subsystems with continuous wave (CW) or with modulated RF sources at the circuit level. Many efforts have already been made in the past decades to reach this objective. In 1996, J.F. Sevic et al. [1] presented an overview of the possibilities of simulating circuits driven by modulated generators in the frame of already published harmonic balance (HB) and behavioral techniques. This article indicates that, at that time, the HB solver could not give the solution using the classical Newton–Raphson technique and the direct inversion of the Jacobian matrix for HB simulations with a high number of frequencies. Later, to overcome the issue of the high number of equations to be solved generated by modulated signals, V. Rizzoli et al. [2] proposed an “inexact Newton Harmonic Balance” method. Most recently, refs. [3,4] give some examples of the multidimensional HB method. For instance, in the Journal of Neural Engineering, a team from Stanford University [3], after periodizing the input signal, managed to perform a neural transient stimulation by piecewise HB techniques. Today, thanks to the work of many mathematicians such as Galerkin [5,6], Krylov [7], or Saad [8], for instance, commercially available simulators (Keysight (Santa Rosa, CA, USA) PATHWAVE Advanced Design System (ADS)®, Cadence (San Jose, CA, USA) SpectreRF®, Mentor Graphics (Wilsonville, OR, USA) Eldo-RF®) allow the HB simulation of circuits with a number of nodes exceeding several hundred thousand and millions of equations to be solved. Therefore, these simulations can now be performed under periodic and almost-periodic regimes with reasonable computation time, with one to three non-harmonically related fundamental frequencies.
While transient and envelope simulations use time-domain integration and need time-domain models of the whole circuit, HB simulation leads to a system of equations solved in the frequency domain, with frequency-domain models, except for the non-linear elements modeled by time-domain equations coupled to Fourier transforms.
The main drawback of this design process is that the designer must continuously switch between simulations in time/frequency domains with different schematics and/or circuit models depending on the desired time/frequency domain simulation.
At microwaves and millimeter waves, passive elements include dispersive transmission lines, interconnections between circuits (wire bonding), or radiating elements, for instance. The models of these passive elements are generally extracted from electromagnetic (EM) or S-parameter ([S]) simulations. After verification of the Kramers–Kronig relationships to ascertain their feasibility conditions [9,10], they are well defined in the frequency domain. But, the results may not be easily transposable into the time domain.
With a unified frequency-domain simulation frame, designers will be able to switch easily from one kind of simulation (continuous wave (CW)) to another one (modulated carrier) without having to modify their workspace.
The application of a unified frequency-domain simulation method at the circuit level, leveraging the harmonic balance technique, is proposed here and constitutes a major advanced improvement for designers of microwave non-linear devices (amplifiers, mixers, and oscillators, for example). This method also saves simulation time.
In the hereby-proposed method, the pseudo-random modulated (PRM) signals are generated by periodizing the pseudo-random symbol sequences, which could have been generated with an external software, such as Matlab® (https://www.mathworks.com/products/matlab.html, accessed on 23 February 2024), for example.
The driving random modulated generator is first transformed in a pseudo-random modulated (PRM) generator with two fundamental frequencies: the carrier frequency f c a r and, on the other hand, the modulating frequency f m o d f r a m e given by the user-defined length of the modulating random symbol sequence T m o d f r a m e . The random modulated generator then becomes an almost-periodic signal generator. Figure 1 illustrates how an M-QAM modulated signal is generated from two initial pseudo-random sequences to the transmitted modulated signal in the plane πoutPA. The TX and the amplifier under test (AUT) are integrated in the same schematic. The equivalent Rx function is developed within the data display.

2. Simulation Principle of Non-Linear Devices with Pseudo-Random Modulated Signals in the Frame of Almost-Harmonic Balance

2.1. Generation of Pseudo-Random Modulated Signal in the Frame of Almost-Harmonic Balance

In the PRM-HB simulation of a non-linear circuit, all voltages or currents of the circuit are defined and represented by the following waveforms in the time domain [8]:
v m o d N , K t = R k = 0 k = K A k , 0 e j φ k , 0 e j 2 π k f m o d f r a m e t + R n = 1 n = N k = K k = K A k , n e j φ k , n e j 2 π k f m o d f r a m e t e j 2 π n f c a r t with n f c a r + k f m o d f r a m e 0
This equation represents the truncated 2D Fourier series expansion of a modulated voltage, in that case, with the following:
  • N: the number of harmonics of the carrier frequency fcar;
  • K: the number of harmonics of the periodic modulating signal with a fundamental frequency fmodrame (also considering a limited number of frequencies created by the non-linearities of the AUT).
N and K , chosen by the user, are the maximum values of order truncation of the harmonic frequencies of the two fundamental frequencies considered in the almost-periodic harmonic balance simulation ( f c a r and f m o d f r a m e ), respectively.
The total number of useful frequencies for the PRM-HB technique is then the following one:
N t o t _ f r e q = N 2 K + 1 + K + 1
From this last equation, the total number of frequencies of interest implemented in the PRM-HB technique can grow very quickly when the orders N and K increase. But, thanks to the GMRES technique [8] (with Krylov bases [11]) and an artificial frequency mapping, both implemented in current circuit solvers, the calculations with large numbers of frequencies involved can be easily achieved. As shown in Figure 1, the pseudo-random modulated electromotive force (EMF) in the πout_EMF plane ( E ~ e m f ( f ) ) is generated thanks to the following:
  • The IQ mapping block (generation of a I ~ f ; Q ~ f couple);
  • The two Delta shaping filters (δshaping(f) used to generate a I ~ δ f ; Q ~ δ f );
  • The two half Nyquist filters (square root raised cosine filter usually used in telecommunication systems allowing the generation of a I ~ δ _ S R R C f ; Q ~ δ _ S R R C f );
  • the quadrature IQ modulator;
  • The Delta shaping filter (δshaping(f)) is defined by the following equation:
δ s h a p i n g f = 1 / s i n c f S y m b R a t e
The demultiplexed part of the binary sequence is realized with Nbit generators of pseudo-random bit sequences clocked at the symbol rate value (called SymbRate).
A last I/Q reference couple called I ~ r e f f ; Q ~ r e f f ) is generated thanks to a Lanczos sigma factor [12] filtering function applied to the I ~ δ _ S R R C f ; Q ~ δ _ S R R C f data. This function drastically reduces the Gibbs phenomena so that the output signals can be used as reference signals to perform the later EVM calculation.

2.2. Example of a Generated Pseudo-Random 16-QAM Modulated Signal in the Frame of Almost-Harmonic Balance

A 16-QAM modulated signal is generated with the parameters’ values of Table 1.
In order to explain the PRM-HB simulation of Figure 1, the AUT is firstly replaced by a through and a load equal to 50 Ω. This simulation is performed using an X64-based desktop PC with Intel® (Santa Clara, CA, USA) Core™ [email protected] GHz, 1.61 GHz, 6 Core(s), and 16 GB of RAM. The simulation runtime is equal to 2.39 s to deal with a total number of frequencies equal to 5506.
After this PRM-HB simulation, the magnitude of the modulated voltages V ~ m o d ( f ) defined in Figure 1 can be directly plotted in the data display window, as given in Figure 2.
This spectrum is centered at the carrier frequency f c a r and demonstrates a simulated 13.5 MHz modulation bandwidth ( B W m o d = 2 B W b b ) corresponding to the theoretical one.

2.3. Demodulation and EVM Calculation after PRM-HB Simulation with the Generated Pseudo-Random 16-QAM Modulated Signal in the Frame of Almost-Harmonic Balance

The demodulation process is not performed in the schematic. This process can be easily implemented, but it increases the number of nodes and the data memory size. A demodulation method is then proposed that drastically alleviates the computation time of the PRM-HB simulation.
The post-processing is carried out in a “data display” window, and it corresponds, in our example, to the digital reception of the 16-QAM modulated signal. It is then applied, for instance, to the calculations of the EVM.
Six post-processing steps are applied to the input and output voltages of the DUT, as shown in Figure 3, with their own outcomes.
In order to describe these six steps, they are summarized in the next subsections with the example of the through connection, and a load equal to 50 Ω. A carrier magnitude voltage equal to 5 V is used. In that case, V ~ m o d f is equal to V ~ r e c f , and the process is only demonstrated for V ~ m o d f . As shown in Section 2.3.6, when the AUT is connected, the EVM calculation process is applied to V ~ m o d f and V ~ r e c f .

2.3.1. Step 1: Implementation of the Raw Envelope Voltage/Current Calculation, around f c a r , as Quadrature Demodulation

It consists of calculating the raw envelopes, around f c a r , of the frequency-domain probe voltages V ~ m o d ( f ) at the input of the DUT and at its output V ~ r e c ( f ) . This calculation is a software-defined way to realize the IQ demodulation in the data display, and it does not need an additional implementation in the schematic, saving a lot of computation time.
Figure 4 shows the magnitude of the raw time-domain envelope v ~ m o d t (dot) overlaid with the RF modulated voltage v m o d N , K t , subsampled to be better visualized, and called v m o d K s u b s a m p t .

2.3.2. Step 2: Application of the Square Root Raised Cosine (SRRC) Filter

The second step of the post-processing consists of performing matched filtering. An SRRC filtering process (lowpass ½ Nyquist) is applied in the frequency domain to the real and imaginary parts, called V ~ m o d I ( f ) and V ~ m o d Q ( f ) , of the input voltage. This filtering process allows for avoiding the intersymbol interference at the optimal sample instant, and consequently, the signal-to-noise ratio is maximized for noisy AUT. The association of the voltage driving the AUT (Tx Part) and the SRRC filtering process of the receiver (lowpass ½ Nyquist) is equivalent to a matched filtering ( H R C f : raised cosine) of a system transmission applied at a circuit level.
This second SQRRC filtering process for the demodulation step is implemented in the frequency domain to V ~ m o d I ( f ) and V ~ m o d Q ( f ) through an in-house XLIM function using the Application Extension Language (AEL) programming language. This function performs a baseband process in the frequency domain applied to voltage variables resulting from the simulation based on the PRM-HB technique.
The equations defining the RC filter are the following ones:
f 1 = S y m b R a t e × 1 r o l l o f f / 2 f 2 = S y m b R a t e × 1 + r o l l o f f / 2 R C f = 1 if f < f 1 R C f = 0.5 × 1 + c o s π f f 1 r o l l o f f if f f 2 R C f = 0 if f > f 2
Figure 5 shows the magnitude of the raw voltage envelope v ~ m o d t overlaid with the RC filtered voltage envelope called v m o d _ R C t .
This matched filtering process is mandatory to ensure that the final demodulated voltages and currents meet the Nyquist criteria.

2.3.3. Step 3: Determination of the Optimal Sampling Instant to Extract the Symbol Sequence: Estimation of the AUT’s Group Delay Noted τ ^ A U T

The third step of post-processing consists of implementing the technique to determine the optimal sampling instant useful for recovering the symbol frame. This technique is based on the estimation of the delay and the correction of the static errors of all the symbols in v ~ r e c _ R C t when the linear voltage gain ( G ~ v l i n A U T ) and the group delay ( τ A U T ) of the AUT are not known beforehand. Note that the static errors are the gain and phase scaling imperfections. In this third step, the group delay between the planes πoutPA and πinPA (Figure 1) τ A U T is estimated (estimation noted τ ^ D U T ) by the maximization of a cost function or by calculating and maximizing the cross-correlation function between v ~ m o d _ R C t and v ~ r e c _ R C t .
Figure 6 shows the fully demodulated input voltage before optimal sampling ( v ~ m o d _ R C t ) and after optimal sampling ( v ~ m o d _ R C _ s a m p p ).

2.3.4. Step 4: Calculation of Static Error Coefficients

The fourth step consists of correcting the “static” errors of the matched filtered and optimally sampled received voltage v ~ r e c _ R C _ s a m p p in order to perform the EVM calculations based on the use of a similar scaled reference vector diagram [13,14]. These errors are considered “static” errors because they appear identically on all the different samples of the received constellations. The determination of six static error coefficients allows the correction of the gain and phase imbalances between the sampled received voltage v ~ m o d _ R C _ s a m p p and the reference sampled symbol sequence determined thanks to the I ~ r e f f ; Q ~ r e f f ) previously calculated, as described in Section 2.1.

2.3.5. Step 5: Correction of the Matched Filtered and Optimally Sampled Extracted Envelopes

With the previous knowledge of the six static error coefficients, the calculation of the corrected vector diagram v ~ m o d _ R C _ s a m p _ c o r r p from the matched filtered and optimally sampled received voltage v ~ m o d _ R C _ s a m p p is performed. The calculation of the corrected trajectory v ~ m o d _ R C _ c o r r t from the received RC filtered trajectory v ~ m o d _ R C t is also achieved with the same coefficients.
The corrected time-domain complex envelope magnitude v ~ m o d _ R C _ c o r r t is then plotted with v ~ m o d t and v ~ m o d _ R C t , as shown in Figure 7.

2.3.6. Step 6: Plot of the EVM Linearity Criterion

After applying these corrections, the EVM can finally be calculated by overlaying the reference 16-QAM vector diagram and the corrected 16-QAM vector diagram v ~ r e c _ R C _ s a m p _ c o r r p , as shown in Figure 8.
The two vector diagrams are perfectly matched and overlaid for this example realized with the AUT replaced by a through and a load equal to 50 Ω and for the previously defined modulated 16-QAM voltage.
The root mean square (RMS) normalized EVM ( E V M R M S ) is then calculated as follows:
E V M R M S ( % ) = p = 1 N s y m b v ~ m o d _ R C _ s a m p _ c o r r p v ~ r e f p 2 p = 1 N s y m b v ~ r e f p 2
The E V M R M S is calculated for v ~ m o d _ R C _ s a m p _ c o r r p and for v ~ r e c _ R C _ s a m p _ c o r r p . For the considered example with a zero-length transmission line (“through”) and a modulated 16-QAM, the E V M R M S is very low and equal to 0.009% at its input and output ports.
All the previous simulation results validate all the post-processes implemented in the final data display of the simulation by HB techniques of pseudo-random modulated signals (PRM-HB), i.e., SRRC filtering, determination of τ ^ D U T , and correction of static errors.
Finally, the PRM-HB simulation method and the implemented data display with the six steps are now the template required to accurately evaluate the EVM linearity criterion of a voltage or current extracted at the input and output of an AUT driven by a 16-QAM modulation. Obviously, other modulation schemes can be identically implemented with this template.

3. HEMT GAN Technology

3.1. 0.25 µm GaN HEMT Technology

Transistors based on 0.25 µm AlGaN/GaN technology (GH25-10) from UMS foundry [15] have been used to design the DPA. They are qualified on a 4-inch-diameter SiC substrate for up to 30 V drain biasing. They are characterized by a gate-to-drain breakdown voltage higher than 120 V. These transistors, when they are optimally matched, demonstrate, at 10 GHz, about 4.5 W/mm of saturated RF output power and a maximum CW peak power-added efficiency (PAE) above 50%.

3.2. Quasi-MMIC Technology in Overmold QFN Package

The designed DPA [16] combines GaN power bars and input/output matching circuits realized on passive GaAs MMIC technology specifically developed for high power density functions. These circuits have been then assembled in a 54 leads 8 × 8 mm overmold plastic package with a high thermal conductivity glue, allowing a lower thermal resistance as compared to standard solutions. A good trade-off between electrical (high frequency) and industrial constraints (assembly rules) is obtained with this low-cost solution. The interconnections between, on the one hand, GaAs passive MMICs and GaN power bars and, on the other hand, the ones between package leads and GaAs devices are realized with wire bonds. All these interconnection networks have been defined according to the industrial rules

3.3. Design Approach

The theory about DPA is well documented [17,18,19,20]. The design needs to consider the constraints of the quasi-MMIC technology thanks to a full 3D EM simulation of the interconnection network into the QFN environment (a key point of the design). The schematic used for the [3.2–3.6] GHz asymmetric packaged Quasi-MMIC SI-DPA is given in Figure 9 [1]. This power amplifier is based on the use of an 8 × 275 mm GH25-10 GaN HEMT for the main amplifier and a 2 × 8 × 275 mm GH25-10 GaN HEMT for the peak amplifier.

4. Principle of Non-Linear Local Stability of the Amplifier under Test Driven by Pseudo-Random Modulated Carrier Generator

This paragraph is dedicated to another application of the here-developed PRM-HB simulation method. This can also be useful to evaluate the non-linear local stability of circuits driven by pseudo-random modulated carrier generator as the previously defined DPA.
For that purpose, the frequency domain procedure described for a large-signal CW generator case [21,22] has been extended to an almost-periodic driving generator and the use of the PRM-HB simulation method, as described in Figure 10.
The local stability of the DPA driven by the pseudo-random modulated carrier generator can then be determined in the frequency domain as follows:
  • A probing small-signal current generator, at frequency j Ω , with Ω q · ω c a r ± p · ω m o d f r a m e is connected between a node k of the circuit and the ground;
  • The non-linear circuit, driven by the two large-signal fundamental frequencies, j ω c a r (carrier) and j ω m o d f r a m e (length of the modulating random bit sequence), appears then as a linear almost periodically time-varying (LAPTV) circuit to the (small signal) probing current generator;
  • The circuit can be simulated in a large-signal small signal (LSSS) mode or in an almost-periodic HB three-tone mode with two large-signal (local oscillators) generators at j ω c a r and j ω m o d f r a m e , and one small-signal generator at j Ω . In both cases, the whole circuit works in the so-called mixer mode;
  • The (small signal) driving port admittance Y kk j Ω seen by the probing current generator is first simulated.
The admittance Y kk j Ω seen by the small-signal current generator in the real frequency domain jω is identified in the complex frequency domain as Y kk p , thanks to an identification procedure performed with the STAN® software (https://www.amcad-engineering.com/stability-analysis/, accessed on 23 February 2024) [23].
The driving point admittances Y kk p share the same numerator at all nodes, which is the determinant of the following whole circuit: D e t M c p .
The numerator D e t M c p captures all the Floquet exponents [24,25,26] ( p z r = σ z r + j Ω z r ) of the whole circuit. The software STAN® determines all the Floquet exponents ( p z r = σ z r + j Ω z r ) of the LAPTV circuit:
  • If the real part of all σ z r zeros is negative, the circuit is locally stable;
  • If there is a positive real part of any σ z r zero, the circuit is locally unstable.
In order to study the local stability in the [0; ω c a r + 2 π · B W b b 2 ] frequency range, and since the circuit appears almost periodic with respect to the perturbation frequency (mixer mode), the small-signal perturbation generator should be theoretically swept in the same frequency range.
The chosen bandwidth 0 < Ω < ω c a r 2 (highlighted in magenta in Figure 11) will allow an appropriate fitting of the complex admittance Y kk j Ω and an accurate identification of Y kk p .
Please recall that the numerical values of the small-signal frequencies must exclude the large-signal frequencies present in the circuit.
The designed DPA is simulated with the configuration described in Table 2.
First, a spectral balance analysis of the large-signal state variables of the circuit is performed. Then, the circuit seen by the perturbation becomes a linear time-varying (LTV) circuit.
Figure 12 presents the Floquet exponents p z r obtained from the STAN® software for three power levels.
It can be easily seen in Figure 12 that the real parts of the calculated zeroes are negative ( σ z r < 0 ). It can be deduced then that the DPA is locally stable for the three input levels. The conclusion is exactly the same for all the swept carrier magnitudes of the PRM-HB local stability simulation of the DPA.
The proposed frequency-domain analysis of local stability extends the method proposed in [21,22] to circuits driven by random modulated RF generators.
The non-linear local stability analysis can be extended to OFDM modulated generators by writing ω m o d f r a m e as the following:
ω m o d f r a m e = 2 π T S y m b o l O F D M
where T S y m b o l O F D M is the length of an integer number of OFDM symbols.

5. Comparison of Simulated and Measured Dynamic Results

5.1. PRM-HB 16-QAM Simulation of the 20 W-S Band Asymmetric DPA

To perform the comparison between simulated and measured dynamic results, the designed DPA is first simulated with a 50 Ω and a 16-QAM PRM generator. The numerical values of the 50 Ω 16-QAM PRM generator can be found in Table 1. The additional or modified parameters of the PRM-HB simulation are detailed in Table 3.
The total number of useful frequencies N t o t _ f r e q for the PRM-HB simulation is then equal to 4406. B W b b is then equal to 15 MHz.
The PRM-HB simulation, performed on the same X64-based desktop PC as previously described, leads to a runtime equal to 9291 s to deal with a sweep of 61 levels of the EMF voltage source.
Figure 13 presents the CW and 16-QAM results for A c a r equal to 4.45 V (large signal).
When the DPA is driven with a large signal in the OBO region, the envelope of the received voltage and the associated constellation are more highly distorted, leading to lower ACPRs and the lowest EVM value. The DPA presents better efficiencies. The dynamic curve of the dynamic drain efficiency ( D E d y n ) , D E d y n t = f P o u t R F t follows the equivalent CW drain efficiency (DE) curve with more dispersion. The peak dynamic output power almost reaches the peak value of the CW curve. The dynamic curve D E d y n t = f P o u t R F t follows the equivalent CW curve with dispersion. The peak dynamic output power almost reaches the peak value of the CW curve.

5.2. Time-Domain Measurement System to Characterize the 20 W-S Band Asymmetric DPA Driven by a 10 MS/s 16-QAM Modulated Voltage

Thanks to a full spectrum time-domain measurement system (as opposed to a bandpass spectrum), the dynamic experimental results are compared with the outcomes obtained from the previously presented PRM-HB simulation results.
Figure 14 describes the proposed six-channel time-domain measurement system based on the use of four track-and-hold amplifiers (THA) [27].
The six-channel measurement test bench contains two 20 dB wideband bidirectional couplers, enabling the simultaneous measurements of the incident and reflected voltage waves at the input and the output of the DUT. Variable calibrated attenuators are required to avoid saturation of THA-based receivers. The signal generated with the vector signal generator is linearly amplified using a broadband high gain amplifier (linear PA in Figure 14) before feeding the DPA. It also employs a large bandwidth [DC: 1.6 GHz] digitizer to measure the fully calibrated [25] coherent RF and LF voltage and current responses.
Two different channels of the eight-channel digitizers are also used to measure the raw LF output voltage and current of the DPA simultaneously with the envelopes of the RF voltages and currents. Note that the raw main and the peak drain voltages and currents are measured simultaneously and not separately, thanks to a unique voltage probe and a unique current probe.
The RF calibration process performed for all the frequencies of the frequency grid (baseband, upper and lower sidebands around f c a r ) is based on three different steps: The first one is a classical SOLR [28,29] VNA calibration at all frequencies of interest. The second one is a power calibration based on the use of a power probe connected to the π i n _ D P A plane (Figure 1) for the upper and lower sidebands around f c a r . The last step is an absolute phase calibration based on a calibrated scope measurement standard [25] for the upper and lower sidebands around f c a r .
The six-channel measurement THA-based test bench presented in Figure 14 is then used to perform the measurement of the coherent RF and LF voltage and current responses of the DPA driven with the same 16-QAM 10 MS/s modulated signal as the one used in the simulation.
Before this measurement, the stability of the DPA driven by a CW and a 16-QAM PRM generator is verified by varying the input power. No precursors or spurious frequencies are found.

5.3. 16-QAM Microwave Measurements of the 20 W-S Band Asymmetric DPA

The DPA is measured with the test set-up described in Figure 14 with the configuration defined in Table 3, except for the modified parameters given in Table 4.
The measured RF metrics of this HPA driven with a modulated signal are based on the measured power waves defined in the frequency domain.
The measured and simulated dynamic P A E % versus P ¯ o u t f 0 dBm curves of the DPA driven with the 16-QAM 10 MS/s modulated voltage are compared in Figure 15.
The measured and simulated dynamic left A C P R versus P o u t R F d y n t ¯ curves of the DPA driven with the 16-QAM 10 MS/s modulated signal are compared in Figure 16.
The measured and simulated dynamic right A C P R versus P o u t R F d y n t ¯ curves of the DPA driven with the 16-QAM 10 MS/s modulated signal are compared in Figure 17.
The measured and simulated dynamic right E V M versus P o u t R F d y n t ¯ curves of the DPA driven with the 16-QAM 10 MS/s modulated signal are compared in Figure 18.
Figure 15, Figure 16, Figure 17 and Figure 18 show a rather good agreement between measured and simulated results obtained with the 16-QAM 10 MS/s modulated signal validating the developed new PRM-HB simulation.

6. Conclusions

The steady-state and the non-linear stability simulations of a load modulated power amplifier (LMPA) driven by a random modulated generator are presented. The simulation is fully performed in the frequency domain by harmonic balance techniques. The demodulation of the output signal of the DUT is implemented, with optimally matched filters, as a software-defined demodulation, saving a lot of computation time.
This article presents a general methodology to perform a whole frequency-domain simulation of non-linear circuits driven by PRM microwave signals in the frame of almost-periodic HB. In the hereby proposed methodology, the PRM microwave signals are generated by periodizing the pseudo-random bit sequences, allowing almost periodic HB.
The simulated dynamic results of a Quasi-MMIC GaN LMPA, a Doherty power amplifier (DPA), are shown and compared to the measured results with a 16-QAM driving signal at 10 MS/s.
The dynamic modulation criteria and power metrics (adjacent channel power ratio (ACPR), error vector magnitude (EVM)) simulated performances at the design circuit level in the frequency domain are calculated and compared to the ones obtained with a specific THA-based time-domain calibrated test bench developed at XLIM.
These comparisons between measured and simulated results present rather good agreements, especially in the OBO region where the EVM presents local minimum values.
To summarize, the PRM-HB simulation tool allows simulations of any high-power amplifier (HPA) driven with any random modulated signals at a circuit level in the frequency domain. It can be extended to multicarrier applications as, for example, satellite transmissions.
The full simulation in the frequency domain, including the steady state and the non-linear stability of an LMPA driven by a random modulated signal, is presented for the first time. This new pseudo-random modulated harmonic balance simulation of the entire DPA is validated by large-signal time-domain measurements of the realized DPA.
This PRM-HB is an advanced and original improvement of the simulation tool of RF circuits and subsystems at the circuit level for telecommunication applications. Effectively, it will help the circuit designers to have a deep insight and understanding of the operation of non-linear circuits when driven with modulated signals. The presented PRM-HB method will effectively complement the joint optimization of power-added efficiency, output power, RF bandwidth, and linearity, which is generally performed using continuous wave harmonic balance during the design process of power amplifiers. The PRM-HB simulation during the first step of the PA design process, at a circuit level, will then enable the electrical knowledge of the power amplifier in the presence of modulated signals as they are in real life with the extraction of criteria at system-level to better optimize their final power performances.

Author Contributions

Resources, C.H. and D.P.; Writing—review & editing, G.N. and D.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors are indebted to J. Obregon for the many enlightening discussions they shared with him. They also acknowledge M. Ayad and E. Richard from UMS Foundry, who provided them with the Gan DPA and its associated schematics.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sevic, J.F.; Steer, M.B.; Pavio, A.M. Non Linear Analysis Methods for the Simulation of Digital Wireless Communication Systems. Int. J. Microw. Millim.-Wave Comput.-Aided Eng. 1996, 6, 197–216. [Google Scholar] [CrossRef]
  2. Rizzoli, V.; Cecchetti, C.; Masotti, D.; Mastri, F. Nonlinear processing of digitally modulated carriers by the inexact-Newton harmonic balance technique. Electron. Lett. 1997, 33, 1760–1761. [Google Scholar] [CrossRef]
  3. Chen, Z.C.; Wang, B.; Palanker, D. Harmonic-Balance Circuit Analysis for Electro-Neural Interfaces. J. Neural Eng. 2020, 17, 3. [Google Scholar] [CrossRef] [PubMed]
  4. Shiri, D.; Nilsson, H.R.; Telluri, P.; Roudsari, A.F.; Shumeiko, V.; Fager, C.; Delsing, P. Modeling and Harmonic Balance Analysis of Parametric Amplifiers for Qubit Read-out. arXiv 2023, arXiv:2306.05177. [Google Scholar]
  5. Galerkin, B.G. Series Solutions of Some Problems of Elastic Equilibrium of Rods and Plates; Vestnik Ingenerov: Petrograd, Russia, 1915; p. 897. (In Russian) [Google Scholar]
  6. Bizzarri, F.; Brambilla, A.; Codecasa, L. Reduction of harmonic balance equations through Galerkin’s method. In Proceedings of the 2015 European Conference on Circuit Theory and Design (ECCTD), Trondheim, Norway, 24–26 August 2015; pp. 1–4. [Google Scholar] [CrossRef]
  7. Vaezi, A.; Abdipour, A.; Mohammadi, A. Nonlinear analysis of microwave amplifiers excited by multicarrier modulated signals using envelop transient technique. Analog Integr. Circ. Sig. Process 2012, 72, 313–323. [Google Scholar] [CrossRef]
  8. Saad, Y.; Schutz, M.H. GMRES: A generalized minimal residual method for solving non-symmetric linear systems. SIAM J. Sci. Stat. Comput. 1986, 7, 856–869. [Google Scholar] [CrossRef]
  9. Saunders, C.S.; Steer, M.B. Passivity Enforcement for Admittance Models of Distributed Networks Using an Inverse Eigenvalue Method. IEEE Trans. Microw. Theory Tech. 2012, 60, 8–20. [Google Scholar] [CrossRef]
  10. Charest, A.; Nakhla, M.; Achar, R.; Saraswat, D. Passivity Verification of Delayed Rational Function Based Macromodels of Tabulated Networks Characterized by Scattering Parameters. IEEE Trans. Compon. Packag. Manuf. Technol. 2011, 1, 386–398. [Google Scholar] [CrossRef]
  11. Krylov, A.N. On the numerical solution of the equation by which the frequency of small oscillations is determined in technical problems. Izv. Akad. Nauk SSSR, Ser.Fiz.-Mat., 4 (1931). In The Theory of Matrices; Gantmacher, F.R., Translator; Chelsea Publishing Co.: New York, NY, USA, 1959; Volume 1–2, pp. 491–539. (In Russian) [Google Scholar]
  12. Lanczos, C. Applied Analysis; Dover Books on Mathematics; Dover Publications: Mineola, NY, USA, 1988; Chapter IV; pp. 225–228. ISBN 13-978-0486656564. Reprint in Prentic-Hall, Inc.: Englewood Cliffs, NJ, USA, 1956. [Google Scholar]
  13. Mckinley, M.; Remley, K.; Myslinski, M.; Kenney, J.; Schreurs, D.; Nauwelaers, B. EVM calculation for broadband modulated signals. In Proceedings of the 64th ARFTG Conference Digest, Orlando, FL, USA, 30 November–3 December 2004; pp. 45–52. [Google Scholar]
  14. Sombrin, J.; Medrel, P. Cross-correlation method measurement of error vector magnitude and application to power amplifier non-linearity performances. In Proceedings of the 88th ARFTG Microwave Measurement Conference (ARFTG), Austin, TX, USA, 8–9 December 2016; pp. 1–4. [Google Scholar] [CrossRef]
  15. Floriot, D.; Brunel, V.; Camiade, M.; Chang, C.; Lambert, B.; Ouarch-Provost, Z.; Blanck, H.; Grünenpütt, J.; Hosch, M.; Jung, H.; et al. GH25-10: New qualified power GaN HEMT process from technology to product overview. In Proceedings of the 2014 9th European Microwave Integrated Circuit Conference, Rome, Italy, 6–7 October 2014; pp. 225–228. [Google Scholar] [CrossRef]
  16. Richard, E.; Ayad, M.; Camiade, M. Recent Development for Linear Amplification Based on GaN Technology. In Proceedings of the 12th European Microwave Integrated Conference (EUMIC) Conference of the EuMW, Recent Advancements in Wide-Band and Effcient GaN Power Amplifiers Workshop, WW-03, Nuremberg, Germany, 8–10 October 2017. [Google Scholar]
  17. Zhou, H.; Perez-Cisneros, J.-R.; Hesami, S.; Buisman, K.; Fager, C. A Generic Theory for Design of Efficient Three-Stage Doherty Power Amplifiers. IEEE Trans. Microw. Theory Tech. 2022, 70, 1242–1253. [Google Scholar] [CrossRef]
  18. Ramella, C.; Camarchia, V.; Piacibello, A.; Pirola, M.; Quaglia, R. Watt-Level 21–25-GHz Integrated Doherty Power Amplifier in GaAs Technology. IEEE Microw. Wirel. Compon. Lett. 2021, 31, 505–508. [Google Scholar] [CrossRef]
  19. Quaglia, R.; Greene, M.D.; Poulton, M.J.; Cripps, S.C. A 1.8–3.2-GHz Doherty Power Amplifier in Quasi-MMIC Technology. IEEE Microw. Wirel. Compon. Lett. 2019, 29, 345–347. [Google Scholar] [CrossRef]
  20. Grebennikov, B.A.; Bulja, S. High-efficiency Doherty power amplifiers: Historical aspect and modern trends. Proc. IEEE 2012, 100, 3190–3219. [Google Scholar] [CrossRef]
  21. Jugo, J.; Portilla, J.; Anakabe, A.; Suarez, A.; Collantes, J.M. Closed-loop stability analysis of microwave amplifiers. Electron. Lett. 2001, 37, 226–228. [Google Scholar] [CrossRef]
  22. Collantes, J.-M.; Mori, L.; Anakabe, A.; Otegi, N.; Lizarraga, I.; Ayllon, N.; Ramírez, F.; Armengaud, V.; Soubercaze-Pun, G. Pole-Zero Identification: Unveiling the Critical Dynamics of Microwave Circuits Beyond Stability Analysis. IEEE Microw. Mag. 2019, 20, 36–54. [Google Scholar] [CrossRef]
  23. Dellier, S.; Gourseyrol, R.; Soubercaze-Pun, G.; Collantes, J.-M.; Anakabe, A.; Narendra, K. Stability analysis of microwave circuits. In Proceedings of the WAMICON 2012 IEEE Wireless & Microwave Technology Conference, Cocoa Beach, FL, USA, 15–17 April 2012; pp. 1–5. [Google Scholar] [CrossRef]
  24. Floquet, G. Sur les équations différentielles linéaires à coefficients périodiques. Ann. Sci. Ec. Norm. Supér. 1883, 12, 47–88. [Google Scholar] [CrossRef]
  25. Cappelluti, F.; Traversa, L.F.; Bonani, F.; Guerrieri, S.D.; Ghione, G. Rigorous, HB-based nonlinear stability analysis of multi-device power amplifier. In Proceedings of the 5th European Microwave Integrated Circuits Conference, Paris, France, 27–28 September 2010; pp. 92–93. [Google Scholar]
  26. Bolcato, P.; Nallatamby, J.C.; Rumolo, C.; Larcheveque, R.; Prigent, M.; Obregon, J. Efficient algorithm for steady-state stability analysis of large analog/RF circuits. In Proceedings of the 2001 IEEE MTT-S International Microwave Sympsoium Digest (Cat. No.01CH37157), Phoenix, AZ, USA, 20–24 May 2001; Volume 1, pp. 451–454. [Google Scholar] [CrossRef]
  27. Ben-Sassi, M.; Neveux, G.; Barataud, D. Ultra-Fast (13ns) Low Frequency/Microwave Transient Measurements, Application to GaN Transistors Characterization of Pulse to Pulse Stability. In Proceedings of the 2019 IEEE MTT-S International Microwave Symposium (IMS), Boston, MA, USA, 2–7 June 2019; pp. 1383–1386. [Google Scholar] [CrossRef]
  28. Ferrero, A.; Pisani, U. Two-port network analyzer calibration using an unknown ’thru’. IEEE Microw. Guid. Wave Lett. 1992, 2, 505–507. [Google Scholar] [CrossRef]
  29. Basu, S.; Hayden, L. An SOLR calibration for accurate measurement of orthogonal on-wafer DUTs. In Proceedings of the 1997 IEEE MTT-S International Microwave Symposium Digest, Denver, CO, USA, 8–13 June 1997; Volume 3, pp. 1335–1338. [Google Scholar] [CrossRef]
Figure 1. Principle of the proposed PRM-HB frequency-domain simulation by HB techniques of pseudo-random modulated signals (this figure shows an example of an M-QAM with N different bit sequence generators).
Figure 1. Principle of the proposed PRM-HB frequency-domain simulation by HB techniques of pseudo-random modulated signals (this figure shows an example of an M-QAM with N different bit sequence generators).
Electronics 13 00947 g001
Figure 2. V ~ m o d ( f ) Spectra in the 50 Ω load after PRM-HB simulation.
Figure 2. V ~ m o d ( f ) Spectra in the 50 Ω load after PRM-HB simulation.
Electronics 13 00947 g002
Figure 3. Principle of the design flow graph of the demodulation and EVM calculation after PRM-HB simulation.
Figure 3. Principle of the design flow graph of the demodulation and EVM calculation after PRM-HB simulation.
Electronics 13 00947 g003
Figure 4. Magnitude of the raw complex envelope, around f c a r , v ~ m o d t (dot), overlaid with the subsampled RF voltage v m o d K s u b s a m p t (line).
Figure 4. Magnitude of the raw complex envelope, around f c a r , v ~ m o d t (dot), overlaid with the subsampled RF voltage v m o d K s u b s a m p t (line).
Electronics 13 00947 g004
Figure 5. Magnitude (zoom) of the raw voltage envelope v ~ m o d t (dot) overlaid with the RC filtered voltage envelope v m o d _ R C t (line).
Figure 5. Magnitude (zoom) of the raw voltage envelope v ~ m o d t (dot) overlaid with the RC filtered voltage envelope v m o d _ R C t (line).
Electronics 13 00947 g005
Figure 6. Trajectories of the fully demodulated input voltage before optimal sampling (lines) and after optimal sampling (circles).
Figure 6. Trajectories of the fully demodulated input voltage before optimal sampling (lines) and after optimal sampling (circles).
Electronics 13 00947 g006
Figure 7. (Zoom) of v ~ m o d _ R C _ c o r r t (green) with v ~ m o d _ R C t (line) and v ~ m o d t (dot).
Figure 7. (Zoom) of v ~ m o d _ R C _ c o r r t (green) with v ~ m o d _ R C t (line) and v ~ m o d t (dot).
Electronics 13 00947 g007
Figure 8. Corrected 16-QAM vector diagram v ~ m o d _ R C _ s a m p _ c o r r p and 16-QAM reference vector diagram.
Figure 8. Corrected 16-QAM vector diagram v ~ m o d _ R C _ s a m p _ c o r r p and 16-QAM reference vector diagram.
Electronics 13 00947 g008
Figure 9. Layout (a) and demonstration board (b) of the realized asymmetric packaged Quasi-MMIC SI-DPA [1].
Figure 9. Layout (a) and demonstration board (b) of the realized asymmetric packaged Quasi-MMIC SI-DPA [1].
Electronics 13 00947 g009
Figure 10. Principle of non-linear LAPTV local stability analysis.
Figure 10. Principle of non-linear LAPTV local stability analysis.
Electronics 13 00947 g010
Figure 11. Frequency sweep (magenta color) of the perturbating generator.
Figure 11. Frequency sweep (magenta color) of the perturbating generator.
Electronics 13 00947 g011
Figure 12. Floquet exponents p z r obtained from the STAN® software for the 3 chosen power levels in CW large-signal regime.
Figure 12. Floquet exponents p z r obtained from the STAN® software for the 3 chosen power levels in CW large-signal regime.
Electronics 13 00947 g012
Figure 13. Main results of the PRM-HB simulation of the DPA driven by CW and modulated 16-QAM EMF for A c a r = 4.45 V (Subfigure (a): V ~ r e c f spectrum—Subfigure (b): 16 QAM reference vector diagram v ~ r e f p overlaid with the corrected 16 QAM Vector Diagram v ~ r e c _ R C _ s a m p _ c o r r p and the received RC filtered Trajectory v ~ r e c _ R C _ c o r r t —Subfigure (c): Dynamic drain efficiency ( D E d y n t % ) versus dynamic output power ( P o u t R F t dBm ) from PRM-HB simulation overlaid with the CW drain efficiency ( D E % ) versus output power ( P ¯ o u t f 0 dBm ) from HB simulation).
Figure 13. Main results of the PRM-HB simulation of the DPA driven by CW and modulated 16-QAM EMF for A c a r = 4.45 V (Subfigure (a): V ~ r e c f spectrum—Subfigure (b): 16 QAM reference vector diagram v ~ r e f p overlaid with the corrected 16 QAM Vector Diagram v ~ r e c _ R C _ s a m p _ c o r r p and the received RC filtered Trajectory v ~ r e c _ R C _ c o r r t —Subfigure (c): Dynamic drain efficiency ( D E d y n t % ) versus dynamic output power ( P o u t R F t dBm ) from PRM-HB simulation overlaid with the CW drain efficiency ( D E % ) versus output power ( P ¯ o u t f 0 dBm ) from HB simulation).
Electronics 13 00947 g013
Figure 14. Six-channel time-domain measurement system for power measurements of non-linear devices, driven by the 10 MS/s 16-QAM modulated voltage and injected in the π i n D P A plane.
Figure 14. Six-channel time-domain measurement system for power measurements of non-linear devices, driven by the 10 MS/s 16-QAM modulated voltage and injected in the π i n D P A plane.
Electronics 13 00947 g014
Figure 15. Simulated and measured power performances of the DPA ( P A E d y n t ¯ v s . P o u t R F d y n t ¯ ) driven with the 16-QAM 10 MS/s modulated voltage.
Figure 15. Simulated and measured power performances of the DPA ( P A E d y n t ¯ v s . P o u t R F d y n t ¯ ) driven with the 16-QAM 10 MS/s modulated voltage.
Electronics 13 00947 g015
Figure 16. Simulated and measured power performances of the DPA ( l eft A C P R v s . P o u t R F d y n t ¯ ) driven with the 16-QAM 10 MS/s modulated voltage.
Figure 16. Simulated and measured power performances of the DPA ( l eft A C P R v s . P o u t R F d y n t ¯ ) driven with the 16-QAM 10 MS/s modulated voltage.
Electronics 13 00947 g016
Figure 17. Simulated and measured power performances of the DPA ( r ight A C P R v s . P o u t R F d y n t ¯ ) driven with the 16-QAM 10 MS/s modulated voltage.
Figure 17. Simulated and measured power performances of the DPA ( r ight A C P R v s . P o u t R F d y n t ¯ ) driven with the 16-QAM 10 MS/s modulated voltage.
Electronics 13 00947 g017
Figure 18. Simulated and measured power performances of the DPA ( EVM vs . P o u t R F d y n t ¯ ) driven with the 16-QAM 10 MS/s modulated voltage.
Figure 18. Simulated and measured power performances of the DPA ( EVM vs . P o u t R F d y n t ¯ ) driven with the 16-QAM 10 MS/s modulated voltage.
Electronics 13 00947 g018
Table 1. Parameters of the generated 16-QAM.
Table 1. Parameters of the generated 16-QAM.
ParameterValueUnit
Symbol Rate10MHz
Symbol Number100
RollOff0.35
Harmonics of the periodic modulating signal500
Carrier Magnitude 5V
Carrier Frequency3.5GHz
harmonics of the carrier frequency1
Table 2. Parameters of the CW local stability simulation of the DPA.
Table 2. Parameters of the CW local stability simulation of the DPA.
ParameterValue with Unit
Quiescent Main Voltage Drain V d s q m a i n = 30 V
Quiescent Main Voltage Gate V g s q m a i n = 3.62 V
Quiescent Peak Voltage Drain V d s q p e a k = 30 V
Quiescent Peak Voltage Gate V g s q p e a k = 5.5 V
CW Frequency f 0 = 3.5 GHz
CW Magnitude 0.1 V A 0 28 V (Step: 0.1 V)
f 0 HB Order N = 5
Table 3. Additional parameters for the PRM-HB simulation of the DPA.
Table 3. Additional parameters for the PRM-HB simulation of the DPA.
ParameterValue with Unit
Quiescent Main Voltage Drain V d s q m a i n = 30 V
Quiescent Main Voltage Gate V g s q m a i n = 3.62 V
Quiescent Peak Voltage Drain V d s q p e a k = 30 V
Quiescent Peak Voltage Gate V g s q p e a k = 5.5 V
CW Frequency f 0 = 3.5 GHz
CW Magnitude 0.1 V A 0 28 V (Step: 0.1 V)
f 0 HB Order N = 5
Table 4. Parameters for the 16-QAM measurement of the DPA.
Table 4. Parameters for the 16-QAM measurement of the DPA.
ParameterValue with Unit
Carrier Frequency f c a r = 3.5 GHz
Carrier Power Sweep 40 d B m P ¯ g e n e f c a r 8.5 d B m (Step: 0.5 dB)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Neveux, G.; Hallepee, C.; Passerieux, D.; Barataud, D. Non-Linear Simulation by Harmonic Balance Techniques of Load Modulated Power Amplifier Driven by Random Modulated Signals. Electronics 2024, 13, 947. https://doi.org/10.3390/electronics13050947

AMA Style

Neveux G, Hallepee C, Passerieux D, Barataud D. Non-Linear Simulation by Harmonic Balance Techniques of Load Modulated Power Amplifier Driven by Random Modulated Signals. Electronics. 2024; 13(5):947. https://doi.org/10.3390/electronics13050947

Chicago/Turabian Style

Neveux, Guillaume, Clément Hallepee, Damien Passerieux, and Denis Barataud. 2024. "Non-Linear Simulation by Harmonic Balance Techniques of Load Modulated Power Amplifier Driven by Random Modulated Signals" Electronics 13, no. 5: 947. https://doi.org/10.3390/electronics13050947

APA Style

Neveux, G., Hallepee, C., Passerieux, D., & Barataud, D. (2024). Non-Linear Simulation by Harmonic Balance Techniques of Load Modulated Power Amplifier Driven by Random Modulated Signals. Electronics, 13(5), 947. https://doi.org/10.3390/electronics13050947

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop