The values of the 6 complex parameters of a large-signal S-parameter model of a high frequency device-under-test are determined by using a frequency-offset probe-tone method. A relatively large one tone signal is applied to the input port of the device and a relatively small one tone signal having a frequency offset relative to the frequency of this large one tone signal is applied to the output port of the device. The 6 large-signal S-parameters are found by measuring and processing the spectral components of the incident and the scattered voltage waves at the device signal ports. These spectral components appear at 3 frequencies: at the frequency of the large one tone signal, at the frequency of the small one tone signal and at the frequency of the large one tone signal minus the frequency offset of the small one tone signal.

Patent
   7038468
Priority
Jun 11 2003
Filed
Jun 01 2004
Issued
May 02 2006
Expiry
Jun 01 2024
Assg.orig
Entity
Large
11
14
all paid
1. A method for characterizing a device under test (DUT) comprising the steps of:
providing a first sinusoidal signal having a first frequency to a first signal port of said DUT as a power tone signal;
providing a second signal having a second frequency to a second signal port of said DUT as a probe tone signal, whereby the difference between said second frequency and said first frequency is equal to a third frequency and whereby the amplitude of said second signal is such that there exists a non-analytic substantially linear relationship between the spectral components of the traveling voltage waves that are incident to said DUT signal ports and that have said second frequency and the spectral components of the traveling voltage waves that are scattered by said DUT signal ports and that have a frequency equal to said first frequency minus said third frequency; and
determining the spectral components of the traveling voltage waves that are incident to said DUT signal ports and the spectral components of the traveling voltage waves that are scattered by said DUT signal ports, whereby said spectral components are determined for a frequency equal to said first frequency, equal to said second frequency and equal to said first frequency minus said third frequency.
2. The method of claim 1, further including the step of using said spectral components for calculating “large-signal S-parameters” of said DUT, whereby said “large-signal S-parameters” include coefficients that relate scattered voltage waves with a complex conjugate of a voltage wave that is incident on said second signal port of said DUT.
3. The method of claim 2, whereby said first signal and said second signal are provided to said DUT in a circuit simulator and whereby said determination of said spectral components is provided by performing a circuit simulation.

This application is entitled to the benefit of Provisional application Ser. No. 60/477,349, filed on Jun. 11, 2003.

Not Applicable

Not Applicable

Not Applicable

Not Applicable

1. Field of Invention

The present invention relates to the measurement of scattering parameters for microwave and radio-frequency (RF) devices-under-test (DUTs) under nonlinear operating conditions.

2. Description of the Related Art

It was shown in J. Verspecht, M. Vanden Bossche and Frans Verbeyst, “Characterizing Components Under Large Signal Excitation: Defining Sensible ‘Large Signal S-Parameters’?!,” 49th ARFTG Conference Digest, pp. 109–117, June 1997 that the nonlinear behaviour of a radio-frequency or microwave device-under-test with a periodic excitation signal at the input in a near matched environment can accurately be described by a ‘Large-Signal S-parameter’ model as defined in the above reference.

One practical method to measure the ‘Large-Signal S-parameter’ model for an actual DUT, restricted to fundamental frequency behaviour only, is described in Frans Verbeyst and Jan Verspecht, “Characterizing non-linear behavior,” European Patent Application “EP 1 296 149 A1”. With this method the ‘Large-Signal S-parameter’ model is identified by connecting to the output port of the DUT a matched load, an open, a short and a plurality of attenuators and delays and by each time measuring the fundamental spectral component of the incident and scattered voltage waves. The ‘Large-Signal S-parameters’ are then calculated by digitally processing the acquired data.

In Jon S. Martens, “Probe Tone S-parameter Measurements,” United States patent application Publication, Pub. No.: “US 2002/0196033 A1”, is described how ‘Hot S22’, which is one part of the complete ‘Large-Signal S-parameter’ model, can be measured by using a frequency-offset probe-tone method. With this probe-tone method a one-tone excitation signal is applied at the input signal port of the DUT (this is called the “power tone”, its frequency is called the fundamental frequency), while a small one-tone signal with a frequency equal to the “power tone” frequency plus a small frequency offset is applied at the output signal port of the DUT (this signal is called the “probe tone”). ‘Hot S22’ is calculated by taking the ratio between the measured reflected voltage wave at the output signal port having a frequency equal to the fundamental frequency plus the said frequency offset and the measured incident scattered voltage wave at the output signal port having a frequency equal to the fundamental frequency plus the said frequency offset.

In accordance with the present invention, a method and an apparatus are provided to determine the values of the parameters of a large-signal S-parameter model of a device-under-test by using a frequency-offset probe-tone method and by measuring and processing the spectral components of the scattered voltage waves with frequencies equal to the fundamental frequency, the fundamental frequency plus said frequency offset and the fundamental frequency minus said frequency offset. Said frequency offset can have a positive as well as a negative value. With said probe-tone method a one-tone excitation signal is applied at the input signal port of the DUT (this is called the “power tone”, its frequency is called the fundamental frequency), while a small one-tone signal with a frequency equal to the “power tone” frequency plus a small frequency offset is applied at the output signal port of the DUT (this signal is called the “probe tone”). The small probe tone in the incident voltage wave at the output signal port of the device gives rise to spectral components in the reflected voltage waves at both DUT signal ports which are proportional to the real and imaginary parts of the complex value of the probe tone phasor. These spectral components appear at 2 frequencies: the fundamental frequency plus the said frequency offset and the fundamental frequency minus the said frequency offset. An apparatus measures the complex values of the phasors of the incident and the scattered voltage waves at the three frequencies of interest: the fundamental frequency, the fundamental frequency plus the said frequency offset and the fundamental frequency minus the said frequency offset. The 6 ‘Large-Signal S-parameters’ are calculated using the measured complex values of the said phasor quantities. The 6 ‘Large-Signal S-parameters’ include the 4 classic coefficients that are related to the spectral components of the incident voltage waves and 2 new coefficients that are related to the complex conjugate of the spectral components of the incident voltage waves.

FIG. 1 shows the naming convention of the incident and the scattered voltage waves

FIG. 2 shows a schematic test setup in accordance with the present invention

FIG. 3 shows a typical spectrum of the a1 voltage wave during the operation of the test setup

FIG. 4 shows a typical spectrum of the a2 voltage wave during the operation of the test setup

FIG. 5 shows a typical spectrum of the b1 voltage wave during the operation of the test setup

FIG. 6 shows a typical spectrum of the b2 voltage wave during the operation of the test setup

FIG. 7 shows an implementation of the invention in an harmonic balance simulator

FIG. 8 shows the simulated amplitude of S22 and Sacc22 as a function of the amplitude of a1

FIG. 9 shows the simulated phase of S22 and Sacc22 as a function of the amplitude of a1

A. Test Method

FIG. 1 depicts a microwave or radio-frequency (RF) device-under-test (DUT) (1) having an input signal port (2) and an output signal port (3). The electromagnetic stimulus and response are described by 4 travelling voltage waves a1, a2, b1 and b2. The voltage waves are defined in a characteristic impedance which is typically 50 Ohm. Voltage wave a1 is incident to the input signal port (2) of the DUT (1), voltage wave a2 is incident to the output signal port (3) of the DUT (1), voltage wave b1 is scattered by the input signal port (2) of the DUT (1) and voltage wave b2 is scattered by the output signal port (3) of the DUT (1). The aforementioned voltage wave definitions are in standard use with a vector network analyzer (VNA). For a linear DUT the relationship between the stimulus (a1, a2) and the response (b1, b2) is described by an S-parameter matrix which can be measured with a standard VNA. When the DUT is an active device, like e.g. a power amplifier, the aforementioned S-parameter description can typically only be used if the input signals are small compared to the operating range of the device. Under large signal stimulus conditions (large compared to the specified operating range of the DUT) the superposition theorem will no longer be valid and S-parameters can no longer be used. So called “Large-Signal S-parameters” are used to describe DUT (1) behaviour under large signal operating conditions, as explained in William H. Leighton, Roger J. Chaffin and John G. Webb, “RF Amplifier Design with Large-Signal S-parameters,” IEEE Transactions on Microwave Theory and Techniques, Vol. MTT-21, No. 12, December 1973. It is shown in Jan Verspecht, Marc Vanden Bossche and Frans Verbeyst, “Characterizing Components Under Large Signal Excitation: Defining Sensible ‘Large Signal S-Parameters’?!,” 49th ARFTG conference digest, June 1997 that a “Large-Signal S-parameter” description having 4 complex parameters (S11, S12, S21, S22), as it is used for a VNA and also in the aforementioned paper by W. H. Leighton et al., is incomplete and needs to be extended by adding two additional complex parameters, here called Sacc12 and Sacc22. These two additional parameters describe the relationship between the conjugate of the a2 incident voltage wave and the b1 and b2 scattered voltage waves. Sacc12 and Sacc22 are exactly equal to zero under linear operating conditions and are as such not measured by a VNA. The object of the present invention is a novel method and test setup to measure the 6 “Large Signal S-parameters” (S11, S12, S21, S22, Sacc12, Sacc22). Consider a one-tone excitation a1 with frequency fc at the input signal port (2), with a large amplitude, and a one-tone excitation a2 having the same frequency fc at the output signal port (3), with a relatively small amplitude. The signal a2 at the output port (3) should be small enough to ensure a linear response of the DUT (1) versus a2. Under such operating conditions, as can be derived from the aforementioned paper by Jan Verspecht et al., the relationship between the response (b1, b2) and the stimulus (a1, a2) is given by:
b1(fc)=S11(|a1(fc)|)a1(fc)+S12(|a1(fc)|)a2(fc)+Sacc12(|a1(fc)|)p2a2*a(fc)  (1)
b2(fc)=S21(|a1(fc)|)a1(fc)+S22(|a1(fc)|)a2(fc)+Sacc22(|a1(fc)|)p2a2*a(fc)  (2),
where
p=ejφ(a1(fc)),
φ(x) represents the phase of complex phasor x,
|x| represents the amplitude of complex phasor x and
x* stands for the conjugate of phasor x.

A method to measure the “Large-Signal S-parameters” S22 and Sacc22 is described in Frans Verbeyst and Jan Verspecht, “Characterizing non-linear behavior,” European Patent Application ‘EP 1 296 149 A1’. With the said method the said “Large-Signal S-parameters” S22 and Sacc22 are measured by connecting to the output signal port (3) of the DUT (1) a matched load, an open, a short and a plurality of attenuators and delays and by each time measuring the fundamental spectral component of the incident and scattered voltage waves. The “Large-Signal S-parameters” S22 and Sacc22 are then calculated by digitally processing the acquired data. In contrast to the present invention, the aforementioned method described in the said European patent application by Frans Verbeyst et al. requires many manipulations which are hard to automate and only determines three “Large-Signal S-parameters”, namely S21, S22 and Sacc22.

The test setup of the present invention is illustrated in FIG. 2. A signal generator (4), called the “Power Tone Generator” is used to generate the necessary one-tone a1 signal. In contrast with the description in the aforementioned European patent application by Frans Verbeyst et al. one will not with the present invention apply different loads to the output signal port (3) of the DUT (1). Instead one will use a signal generator (5), called the “Probe Tone Generator”, for the generation of an a2 one-tone signal with a frequency slightly offset from fc and equal to fc+fd. In order for the method to be accurate the frequency difference fd, called the offset frequency, should be that small that the S-parameters do not vary across a frequency difference of fd. The frequency difference fd can have a positive as well as a negative value. In practice one typically uses an offset frequency of several kHz. Since fc is a microwave or RF frequency it is typically much higher than 100 MHz, such that the S-parameters do not significantly change for a kHz frequency change. This method is known as an “offset frequency probe tone” method and is described in J. Martens, “Probe Tone S-parameter Measurements,” United States patent application Publication, Pub. No. US 2002/0196033 A1, Dec. 26 2002. The resulting spectra of the voltage waves a1, a2, b1 and b2 during one “offset frequency probe tone” experiment are illustrated in FIG. 3, FIG. 4, FIG. 5 and FIG. 6. For a1 (FIG. 3) there is only one spectral component present with a frequency fc and for a2 (FIG. 4) there is only one spectral component present with a frequency fc+fd. For b1(FIG. 5) and b2 (FIG. 6) there are spectral components present with frequencies fc−fd, fc and fc+fd. Only the presence of the spectral components with frequencies fc and fc+fd is recognized and used in prior art. The existence of b1 and b2 spectral components with a frequency fc−fd is not recognized nor used in prior art, in contradiction with the present invention.

In the aforementioned patent application publication by J. Martens, the “offset frequency probe tone” method is used for the determination of “Hot S22” only. “Hot S22” corresponds to the “Large-Signal S-parameter” S22. The value of “Hot S22” is found by sensing and measuring the voltage waves b2(fc+fd) and a2(fc+fd) and by calculating the ratio of the two measured quantities:
S22=b2(fc+fd)/a2(fc+fd)  (3).

It is not uncommon to measure the values a1(fc), b1(fc) and b2(fc) with a similar setup in order to calculate the “Large-Signal S-parameters” S11 and S21 as illustrated in “Measurement of Large Signal Input Impedance During Load Pull,” Maury Microwave Corporation Application Note 5C-029, November 1998. The presence of the probe tone signal a2(fc+fd) is not necessary for the measurement of S11 and S21, which are calculated by using the following equations:
S11=b1(fc)/a1(fc)  (4)
and
S21=b2(fc)/a1(fc)  (5).

With the present invention one will use an “offset frequency probe tone” method and one will not only measure the components a1(fc), a2(fc+fd), b1(fc), b1(fc+fd), b2(fc) and b2(fc+fd), as is done in the described prior art, but one will include the measurement of the components b1(fc−fd) and b2(fc−fd). The mere existence of these frequency components is not recognized nor used in prior art. With the present invention one will measure the values of a1(fc), a2(fc+fd), b1(fc−fd), b1(fc), b1(fc+fd), b2(fc−fd), b2(fc) and b2(fc+fd) and one will calculate the values of the “Large-Signal S-parameters” S22, S11 and S21 as described in (3), (4) and (5). As will be explained in the following, the values of the “Large-Signal S-parameters” S12, Sacc12 and Sacc22 are given by:
S12=b1(fc+fd)/a2(fc+fd)  (6)
Sacc12=b1(fc−fd)/(p2a2*(fc+fd))  (7)
Sacc22=b2(fc−fd)/(p2a2*(fc+fd))  (8)
where
p=ejφ(a1(fc)).

The above results can be derived as follows.

First one rewrites equations (1) and (2) using a complex envelope representation for the signals a1, a2, b1 and b2 (using fc as the carrier frequency). Complex envelope representations are mixed time frequency domain representations of signals and are explained in equation (1) of the detailed description of the invention described in U.S. Pat. No. 5,588,142 “Method for simulating a circuit”. In short, X(t) is the complex envelope representation (with carrier frequency fc) of a signal x(t) if the following holds.
x(t)=Real(X(t)ej2πfct)  (9)

This results in the following equations.
B1(t)=S11A1(t)+S12A2(t)+Sacc12A2*(t)P2(t)  (10)
B2(t)=S21A1(t)+S22A2(t)+Sacc22A2*(t)P2(t)  (11)
where
P(t)=ejφ(A1(t)).

During a “frequency offset probe tone” measurements A1(t) and A2(t) are given by the following equations.
A1(t)=K  (12)
A2(t)=Lej2πfdt  (13)

In the above equations K and L are 2 complex constants having a real and an imaginary part.

Substituting (12) and (13) in (10) and (11) results in the following equations.
B1(t)=S11K+S12Lej2πfdt+Sacc12L*ej2φ(K)e−j2πfdt  (14)
B2(t)=S21K+S22Lej2πfdt+Sacc12L*ej2φ(K)e−j2πfdt  (15)

Next equations (14) and (15) are converted into the frequency domain by applying a Fourier transform. The function δ(f) will be introduced. It is defined as being equal to 0 for all f different from 0 and with δ(0) equal to 1. This results in the following.
b1(f)=S11Kδ(f−fc)+S12Lδ(f−fc−fd)+Sacc12L*ej2φ(K)δ(f−fc+fd)  (16)
b2(f)=S12Kδ(f−fc)+S22Lδ(f−fc−fd)+Sacc22L*ej2φ(K)δ(f−fc+fd)  (17)

For the described “frequency offset probe tone” method K and L are given by the following.
K=a1(fc)  (18)
L=a2(fc+fd)  (19)

Substitution of (18) and (19) in (16) and (17) results in the following.
b1(f)=S11a1(fc)δ(f−fc)+S12a2(fc+fd)δ(f−fc−fd)+Sacc12a2*(fc+fd)p2δ(f−fc+fd)  (20)
b2(f)=S21a1(fc)δ(f−fc)+S22a2(fc+fd)δ(f−fc−fd)+Sacc22a2*(fc+fd)p2δ(f−fc+fd)  (21)

Subsequent substitution of f by fc, fc+fd and fc−fd leads to the results described in equations (3), (4), (5), (6), (7) and (8).

B. Test Setup

A typical test setup which implements the aforementioned method in practice is described in FIG. 2. A test set (6) is connected to the input signal port (2) of the DUT (1), the output signal port (3) of the DUT (1), a “Power Tone Generator” (4), a “Probe Tone Generator” (5) and a microwave and RF receiver (11). The a1 signal generated by the “Power Tone Generator” (4) is guided by the test set (6) to the input signal port (2) of the DUT (1) and the a2 signal generated by the “Probe Tone Generator” (5) is guided by the test set (6) to the output signal port (3) of the DUT (1). The test set senses the voltage waves a1, a2, b1 and b2 and guides the sensed voltage waves to a microwave and RF receiver (11). Four couplers (7), (8), (9) and (10) are typically used for sensing the voltage waves. Other possibilities for sensing the voltage waves are electrical and magnetic field probes and resistive bridges.

The output of the receiver (11) are the measured values of a1(fc), a2(fc+fd), b1(fc−fd), b1(fc), b1(fc+fd), b2(fc−fd), b2(fc) and b2(fc+fd). These values are transfered to a signal processor (12) for calibration purposes and for the calculation of the 6 “Large-Signal S-parameters” S11, S12, Sacc12, S21, S22 and Sacc22 according to equations (3), (4), (5), (6), (7) and (8).

In a typical experiment the amplitude of a1(fc) and the frequency fc are swept across a range specified by the user.

In a practical test setup a1(fc−fd), a1(fc+fd), a2(fc−fd) and a2(fc) may be slightly different from zero due to parasitical reflections in the test set. In that case the use of equations (3), (4), (5), (6), (7) and (8) will result in a small error. This error can be corrected by using calibration techniques which make use of calculations based upon the 6 measured quantities a1(fc), a2(fc+fd), a1(fc−fd), a1(fc+fd), a2(fc−fd) and a2(fc).

To those skilled in the art it is clear that the presented novel method and test setup can easily be extended to the measurement of “Large-Signal S-parameters” for DUTs having more than two signal ports.

C. Simulation

The present invention for determining the 6 “Large-Signal S-parameters” of a DUT (1) can also be used in circuit and system simulators. The described “Large-Signal S-parameter” model can be used as a fast and accurate black-box behavioral model which can easily be derived from complex transistor circuits by simulating the test setup as described above. An example of the implementation of the present invention in a commercial harmonic balance simulator is shown in FIG. 6. FIG. 7 shows the simulated amplitude of S22 and Sacc22 for a swept a1(fc) amplitude. FIG. 8 shows the simulated phase of S22 and Sacc22 for the same swept a1(fc) amplitude. Note that the amplitude of Sacc22 tends to zero for a small amplitude of a1(fc). This is what one expects since the linear small-signal S-parameter model is valid under those conditions, which implies a non-existent or in other words zero-valued Sacc22.

Verspecht, Jan

Patent Priority Assignee Title
10003453, Nov 13 2015 Anritsu Company Phase synchronization of measuring instruments using free space transmission
10116432, Nov 13 2015 Anritsu Company Real-time phase synchronization of a remote receiver with a measuring instrument
11199568, Sep 28 2018 Maury Microwave, Inc. Millimeter wave active load pull using low frequency phase and amplitude tuning
7437693, Mar 31 2005 Cadence Design Systems, Inc. Method and system for s-parameter generation
7486067, Jun 07 2004 National Instruments Ireland Resources Limited Real-time device characterization and analysis
7545152, Jan 11 2006 Thermo Keytek LLC Transmission line pulse measurement system for measuring the response of a device under test
7548069, Jun 10 2005 MAURY MICROWAVE, INC Signal measurement systems and methods
7671605, Jan 17 2008 Keysight Technologies, Inc Large signal scattering functions from orthogonal phase measurements
7924026, Mar 24 2008 Keysight Technologies, Inc Method and apparatus for determining a response of a DUT to a desired large signal, and for determining input tones required to produce a desired output
9214718, Mar 08 2012 Apple Inc. Methods for characterizing tunable radio-frequency elements
9541592, Sep 03 2013 TSIRONIS, CHRISTOS Noise parameter measurement system
Patent Priority Assignee Title
4816767, Jan 09 1984 Agilent Technologies Inc Vector network analyzer with integral processor
5036479, Apr 20 1989 Northrop Grumman Corporation Modular automated avionics test system
5336988, Feb 10 1993 Picker International, Inc. Scalar S-parameter test set for NMR instrumentation measurements
5588142, May 12 1995 Keysight Technologies, Inc Method for simulating a circuit
6066953, Dec 01 1994 Teradyne, Inc. Architecture for RF signal automatic test equipment
6292000, Sep 02 1998 Anritsu Company Process for harmonic measurement with enhanced phase accuracy
6300775, Feb 02 1999 Keysight Technologies, Inc Scattering parameter calibration system and method
6348804, Jun 10 1999 Rohde & Schwarz GmbH & Co. KG Vector network analyzer
6529844, Sep 02 1998 Anritsu Company Vector network measurement system
6594604, Oct 13 1999 SILICON VALLEY BANK, AS ADMINISTRATIVE AGENT S-parameter measurement system for wideband non-linear networks
6766262, May 29 2002 Anritsu Company Methods for determining corrected intermodulation distortion (IMD) product measurements for a device under test (DUT)
20020196033,
20030027197,
EP1296149,
//
Executed onAssignorAssigneeConveyanceFrameReelDoc
Dec 21 2011VERSPECHT, JANAgilent Technologies, IncASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0298370988 pdf
Aug 01 2014Agilent Technologies, IncKeysight Technologies, IncASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0337460714 pdf
Date Maintenance Fee Events
May 08 2009M2551: Payment of Maintenance Fee, 4th Yr, Small Entity.
Oct 02 2013M1552: Payment of Maintenance Fee, 8th Year, Large Entity.
Oct 31 2013R2552: Refund - Payment of Maintenance Fee, 8th Yr, Small Entity.
Oct 31 2013STOL: Pat Hldr no Longer Claims Small Ent Stat
Oct 19 2017M1553: Payment of Maintenance Fee, 12th Year, Large Entity.


Date Maintenance Schedule
May 02 20094 years fee payment window open
Nov 02 20096 months grace period start (w surcharge)
May 02 2010patent expiry (for year 4)
May 02 20122 years to revive unintentionally abandoned end. (for year 4)
May 02 20138 years fee payment window open
Nov 02 20136 months grace period start (w surcharge)
May 02 2014patent expiry (for year 8)
May 02 20162 years to revive unintentionally abandoned end. (for year 8)
May 02 201712 years fee payment window open
Nov 02 20176 months grace period start (w surcharge)
May 02 2018patent expiry (for year 12)
May 02 20202 years to revive unintentionally abandoned end. (for year 12)