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
and, on the other hand, the modulating frequency
given by the user-defined length of the modulating random symbol sequence
. 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]:
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).
and , 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 ( and ), respectively.
The total number of useful frequencies for the PRM-HB technique is then the following one:
From this last equation, the total number of frequencies of interest implemented in the PRM-HB technique can grow very quickly when the orders
and
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 (
) is generated thanks to the following:
The IQ mapping block (generation of a couple);
The two Delta shaping filters (δshaping(f) used to generate a );
The two half Nyquist filters (square root raised cosine filter usually used in telecommunication systems allowing the generation of a );
the quadrature IQ modulator;
The Delta shaping filter (δshaping(f)) is defined by the following equation:
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
) is generated thanks to a Lanczos sigma factor [
12] filtering function applied to the
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
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 and demonstrates a simulated 13.5 MHz modulation bandwidth 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,
is equal to
and the process is only demonstrated for
. As shown in
Section 2.3.6, when the AUT is connected, the EVM calculation process is applied to
and
.
2.3.1. Step 1: Implementation of the Raw Envelope Voltage/Current Calculation, around , as Quadrature Demodulation
It consists of calculating the raw envelopes, around , of the frequency-domain probe voltages at the input of the DUT and at its output . 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
(dot) overlaid with the RF modulated voltage
, subsampled to be better visualized, and called
.
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 and , 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 (: 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 and 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:
Figure 5 shows the magnitude of the raw voltage envelope
overlaid with the RC filtered voltage envelope called
.
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
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
when the linear voltage gain (
and the group delay (
) 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)
is estimated (estimation noted
) by the maximization of a cost function or by calculating and maximizing the cross-correlation function between
and
.
Figure 6 shows the fully demodulated input voltage before optimal sampling (
) and after optimal sampling (
).
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
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
and the reference sampled symbol sequence determined thanks to the
) 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 from the matched filtered and optimally sampled received voltage is performed. The calculation of the corrected trajectory from the received RC filtered trajectory is also achieved with the same coefficients.
The corrected time-domain complex envelope magnitude
is then plotted with
and
, 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
, 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 (
) is then calculated as follows:
The is calculated for and for . For the considered example with a zero-length transmission line (“through”) and a modulated 16-QAM, the 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 , 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 , with is connected between a node k of the circuit and the ground;
The non-linear circuit, driven by the two large-signal fundamental frequencies, (carrier) and (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 and , and one small-signal generator at . In both cases, the whole circuit works in the so-called mixer mode;
The (small signal) driving port admittance seen by the probing current generator is first simulated.
The admittance
seen by the small-signal current generator in the real frequency domain jω is identified in the complex frequency domain as
, 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 share the same numerator at all nodes, which is the determinant of the following whole circuit: .
The numerator
captures all the Floquet exponents [
24,
25,
26] (
) of the whole circuit. The software STAN
® determines all the Floquet exponents (
) of the LAPTV circuit:
If the real part of all zeros is negative, the circuit is locally stable;
If there is a positive real part of any zero, the circuit is locally unstable.
In order to study the local stability in the [0; ] 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
(highlighted in magenta in
Figure 11) will allow an appropriate fitting of the complex admittance
and an accurate identification of
.
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
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
). 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
as the following:
where
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 for the PRM-HB simulation is then equal to 4406. 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
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 ( 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 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
) 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
plane (
Figure 1) for the upper and lower sidebands around
. The last step is an absolute phase calibration based on a calibrated scope measurement standard [
25] for the upper and lower sidebands around
.
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
versus
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
versus
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
versus
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
versus
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.