• Skip to main content
  • Skip to search
  • Skip to footer
Cadence Home
  • This search text may be transcribed, used, stored, or accessed by our third-party service providers per our Cookie Policy and Privacy Policy.

  1. Blogs
  2. RF Engineering
  3. Measuring 2-Tone Intermodulation Using Envelope-Following…
Tawna
Tawna

Community Member

Blog Activity
Options
  • Subscribe by email
  • More
  • Cancel
CDNS - RequestDemo

Try Cadence Software for your next design!

Free Trials
RF
RF Simulation
analog/RF
envelope
RFIC
HB
Spectre RF
spectre spectreRF
RF designer
Analog Simulation
MMSIM
RF Block Simulation
Virtuoso Spectre Simulator GXL
fast envelope
analog
ADE
RF spectre spectreRF
spectreRF
Spectre
RF design
harmonic balance
SpectreRF tutorials
simulation
wireless integrated circuit verification

Measuring 2-Tone Intermodulation Using Envelope-Following Analysis

16 May 2012 • 4 minute read

From time to time, SpectreRF users simulate very large, extracted-view circuits in 2+ tone QPSS. In many of those cases, memory requirements exceed the available resources. When that happens and small-signal approximations aren’t applicable, the user is typically stuck.

The solution below and attached database describe a technique that achieves an order-of-magnitude reduction in memory requirements by converting an N-tone QPSS problem to an equivalent (N-1)-tone envelope-following analysis.

Basic Idea

This is a 2-tone signal in the common HB representation:

imageSolution1.jpg

We can also look at this as a modulated carrier, with the modulation function:

im.jpg

We run an envelope analysis to compute the modulations at each node and harmonic and calculate the multi-tone representation using DFT.

Motivation:

* Numerically, envelope simulation is similar to HB simulation repeated at each time step.

* By looking at a 2-tone signal as a modulated carrier, we convert a 2-tone problem to many smaller 1-tone problems.

* In general, we can convert an N-tone HB analysis to a series of N-1 tone analyses.

* For very large circuits with huge memory requirements, this can reduce the required memory by an order of magnitude – effectively, the same problem is now computed using much smaller, N-1 tone systems.

Basic Approach in Virtuoso:

* Use a Verilog-A module and transient simulation to generate I/Q components of the 2-tone modulation function

* Apply the I/Q components using the PORT element’s Modulation parameters and run an envelope simulation

* Run envelope analysis over N periods (the period is just the inverse of a half of the frequency spacing); this lets all transients expire (we are looking to reach the envelope steady-state)

* Use adaptive step envelope for best simulation accuracy, but use equally spaced strobe output for best postprocessing accuracy using DFT.

Applicability:

* This currently works with 2-tone simulation.

* This should also work with 3+ tones, but currently does not. CCR 979403 has been filed for this issue. Once fixed, you will be able to apply this to 3+ tones.

* PAC
– just set the power of one of the 2 tones to 0 and the other to something suitably small. In theory, this technique uses less memory than PAC.

Example:

In this example, you wil see:

* Creating the IQ signal

* Running envelope simulation

* Postprocessing

* Results comparison

Creating the IQ signal

Use the attached database testcase.tar.gz.   This is available on Cadence Online Support in Solution 11774216.

Open library test_2t_env.

Simulate the cell generate_envelope from state spectre_state1.

 

im2.jpg

 

Generating the 2-tone Envelope

The ADE setup lets you specify tone spacing, power in each tone, the number of samples per period (N) and the number of periods (cycles) to store.

You get something like this:

im3.jpg

 

This now represents the envelope of our two-tone signal, and we’ll run and envelope analysis with it as the source next.

Verifying the Source Signal

First, we’ll run a very simple example, just to make sure that we get what we expect.

Open the test_port schematic view and observe the port’s relevant settings.

im4.jpg

 

Port Edit Properties form:

im5.jpg

 

 

Running the Simulation:

Open test_port->spectre_state1 and simulate.

Plot the ‘right’ and ‘left’ outputs: (Note that 'right' and 'left' are defined below).

im6.jpg

 

What are 'right' and 'left'??

* ‘Right’/’Left’ is what you would see to the right/left of the carrier when you look at the spectrum analyzer.

* This is just an artifact of how we post-process the envelope (to be explained later).

* Since the ViVA DFT works on purely real waveforms, ‘left’ is really the image of the negative frequency spectrum.

* So, we have a -20 dBm tone, offset 1 MHz from the carrier, and another of equal power, offset -1 MHz from the carrier. This is exactly what you’d expect.

 

Running a 2-tone HB Simulation for Reference

Open test_hb_2tone, which is just a behavioral LNA simulation.

Run test_hb_2tone, which is the excitation we saw in the above paragraphs, applied in an HB simulation.

 

im7.jpg

Running an Envelope Simulation on the LNA

Run test_env_2tone->spectre_state1.

This is the same as the previous example, except that uses envelope simulation with the PORT element as set up earlier.

Since the LNA has no memory and the excitation is symmetric, ‘left’ and ‘right’ are the same.

 

im8.jpg

Comparison Summary

HB and envelope give nearly identical results!

 

im9.jpg

 

Post Processing Background

* To calculate the spectrum, we would ideally take the complex DFT of the waveform and convert to dBm. The formula goes like:

  Spectrum=dbm(DFT(X(t)))

* Since there is no complex DFT in ViVA, we do something like:

  Right Spectrum = dbm(DFT(real(X(t)/2)+j*DFT(imag(X(t)/2))

  Left Spectrum = dbm(conj(DFT(real(X(t)/2)))-j*conj(DFT(imag(X(t)/2))))

* Note that, in all cases, we want to process only the last cycle of the waveform (when the steady-state is reached).

* The expressions are a bit messy. They are shown only for reference below. Use the ADE state as a bench rather than attempting to enter the expression from scratch.

Spectrum Post Processing Expressions 'Right'

(10 * log10((pow(abs(((dft(real(harmonic(v("/out" ?result "envlp_fd") '1)) ((pv("/cycles" "value" ?result "variables") - 1) / pv("/deltaf2" "value" ?result "variables")) (pv("/cycles" "value" ?result "variables") / pv("/deltaf2" "value" ?result "variables")) pv("/N" "value" ?result "variables") "Rectangular" 1 "default") + (sqrt(-1) * dft(imag(harmonic(v("/out" ?result "envlp_fd") '1)) ((pv("/cycles" "value" ?result "variables") - 1) / pv("/deltaf2" "value" ?result "variables")) (pv("/cycles" "value" ?result "variables") / pv("/deltaf2" "value" ?result "variables")) pv("/N" "value" ?result "variables") "Rectangular" 1 "default"))) / 2)) 2) * 10)))

Spectrum Post Processing Expressions 'Left'

(10 * log10((pow(abs(((dft(real(harmonic(v("/out" ?result "envlp_fd") '1)) ((pv("/cycles" "value" ?result "variables") - 1) / pv("/deltaf2" "value" ?result "variables")) (pv("/cycles" "value" ?result "variables") / pv("/deltaf2" "value" ?result "variables")) pv("/N" "value" ?result "variables") "Rectangular" 1 "default") - (sqrt(-1) * dft(imag(harmonic(v("/out" ?result "envlp_fd") '1)) ((pv("/cycles" "value" ?result "variables") - 1) / pv("/deltaf2" "value" ?result "variables")) (pv("/cycles" "value" ?result "variables") / pv("/deltaf2" "value" ?result "variables")) pv("/N" "value" ?result "variables") "Rectangular" 1 "default"))) / 2)) 2) * 10)))

Have fun simulating!

Best Regards,

Tawna

 


CDNS - RequestDemo

Have a question? Need more information?

Contact Us

© 2025 Cadence Design Systems, Inc. All Rights Reserved.

  • Terms of Use
  • Privacy
  • Cookie Policy
  • US Trademarks
  • Do Not Sell or Share My Personal Information