The present invention relates to wide-band radio frequency (RF) receivers, and more particularly to a system for calibrating wide-band RF receivers for amplitude flatness and phase linearity.
Increased data bandwidth requirements have driven a need for RF receivers that not only receive signals with very wide-band modulation, but also that maintain a flat amplitude and linear phase response over a wide bandwidth. As shown in
The RF receiver 10 is usually calibrated as part of the manufacturing process, and a finite impulse response (FIR) filter 26 is placed at the output of the ADC 22 prior to input of the digital sample stream to the DSP 24. The FIR filter coefficients are set by the manufacturer's calibration process, and serve to correct for amplitude flatness and phase linearity under those ideal conditions. However, each of the mixers, filters, amplifiers and other circuitry in the signal path from the input to the DSP 24 contribute amplitude ripples and group delay ripples, i.e., deviations from linear phase. These active components often exhibit temperature dependence in gain as well as changes in their frequency response. Further, as components age the amplitude and phase response also may be affected. It is often required to improve amplitude accuracy, amplitude flatness and phase linearity of the RF receiver 10.
What is desired is a method of calibrating the RF receiver at the time of use, i.e., at run-time, to correct amplitude flatness and phase linearity for temperature and other environmental dependencies in order to improve the accuracy of the measurement or demodulation algorithms used in the DSP.
Accordingly the present invention provides for magnitude and phase response calibration of RF receivers that provides amplitude flatness and phase linearity over a wide frequency bandwidth using a simple square law diode detector and a frequency-stepped two-tone source. A two-tone source generator provides two sinusoidal signals separated by a specified frequency delta about a center frequency. The center frequency is stepped across the wide frequency bandwidth of the RF receiver. At each center frequency the two sinusoidal signals of equal amplitude and of the same phase are input to the RF receiver and to the square law diode detector in an integrated calibrator. The two sinusoidal signals are processed by both the receiver path and the calibrator path, and the magnitude and phase results from the receiver path are combined with the beat frequency results from the calibrator path to generate coefficients for a finite impulse response filter in the receiver path at each of the stepped center frequencies.
The advantages and novel features of the present invention are apparent from the following detailed description when read in conjunction with the appended claims and attached drawing figures.
Referring again to
The two-tone source 32 produces a signal centered at ωm with tone separation of ωΔ. The signal at the outputs of the power splitter 34 are described mathematically as
The diode 36 functions as a square law detector with a scaling factor GD(ω). The output of the diode 36 is
The output signal from the diode 36 is then filtered by the anti-aliasing filter 38 before being digitized by the ADC 40. The filter 38 removes the twice frequency terms, leaving
The gain of the diode 36 has both a magnitude, GD, and phase, φD, response. The diode 36 is a very broad band device whose magnitude response is flat over the spacing between the two tones.
GD(ω)=|GD(ω)|ejθ
|GD(ω1)|=|GD(ω2)|=|GD(ωm)|
Both the amplitude information and phase information are contained in the filtered signal. The DC level for each value of ωm is
The amplitude of the sinusoidal component is
|ybeat(ωm)|=|GD(ωm)|a1a2
The phase difference between the two tones also is computed from the phase of the sinusoidal beat, not at the output of the diode 36.
θbeat(ωm)=(Ø2−Ø1)+ØD(ω2)−ØD(ω1)
The signal in the receiver signal path is filtered and frequency converted, as described above. The mixer 14 translates the input frequency by ωLO. For simplicity, all frequency response terms are aggregated into an equivalent filter that encompasses all of the amplitude and phase responses of the entire signal path up to the ADC 22, Hr(ω).
Hr(ω)=|Hr(ω)|ejØ
The signal at the input to the ADC 22 is described by
xIF(t)=a1|Hr(ω)|cos [ωLO−ω1)t−Ø1+ØR(ω1)]+a2|Hr(ω2)|cos((ωLO−ω2)t−Ø2+ØR(ω2))
The signal contains two components separated by ωΔ.
Fourier processing is performed on the received signal to determine the amplitude and phase of the two components as measured by the receiver DSP 24. The amplitude and phase of the lower and upper of the two tones is
A1(ω1)ejθ
A2(ω2)ejθ
The product of the magnitude transfer function is calculated by taking the ratio of the products of the amplitudes of the two tones taken from the Fourier transform of the receiver DSP 24 to the level of the beat note measured by the square law detector, diode 36.
If the spacing between the two tones is small enough, then the magnitude response of the two tones is equal.
The magnitude response of the receiver 10 is computed for each value of ωm by stepping ωm across the receiver's bandwidth and solving the following equation:
The phase response also is computed from the phase readings taken from Fourier transform of the receiver DSP 24 and the phase of the beat note from the diode 36 as determined by the calibrator DSP 42.
θ1(ω1)−θ2(ω2)+θbeat(ωm)=[Ø1−ØR(ω1)]+[Ø2−ØR(ω2)]−[(Ø2−Ø1)+ØD(ω2)−ØD(ω1)]
θ1(ω1)−θ2(ω2)+θbeat(ωm)=ØR(ω1)−ØR(ω2)+ØD(ω1)−ØD(ω2)
ØR(ω1)−ØR(ω2)=θ1(ω1)−θ2(ω2)+θbeat(ωm)−ØD(ω1)+ØD(ω2)
The group delay of the receiver path is computed from
The phase response of the receiver path is computed from the group delay by performing an integration.
Ø(ωm)=−∫τg(ωm)dωm+Ø0
There are many techniques for determining FIR coefficients from a complex transfer function. One such method is to take the inverse Fourier transform of the complex transfer function. Another method is called the “window” method, as is well known to those in the art.
In summary as shown in
1. Divide the receiver channel BW into M discrete frequencies, ωm.
2. For each value of ωm.
3. Compute the magnitude response using
4. Compute the group delay response using
5. Compute the phase response by integrating the group delay
Ø(ωm)=∫τg(ωm)dωm+Ø0
The integration process can have an arbitrary fixed offset Ø0
6. Convert the phase response into FIR coefficients.
Thus the present invention provides an integrated calibrator for a receiver signal path that measures both magnitude and group delay (phase after integration), the calibrator having a two tone source with constant separations and selectable center frequency and including a calibrated diode detector and a method for digitizing the diode output, where the two tones are input to the receiver signal path and the calibrator path, measured in each path and combined by the calibrator to generate filter coefficients for a correction filter in the receiver signal path to produce amplitude flatness and phase linearity across the RF bandwidth of the receiver.
Number | Name | Date | Kind |
---|---|---|---|
5337014 | Najle et al. | Aug 1994 | A |
5656929 | Humpherys | Aug 1997 | A |
7477875 | Zhang et al. | Jan 2009 | B2 |
20030176174 | Seppinen et al. | Sep 2003 | A1 |
20090149137 | Rockwell et al. | Jun 2009 | A1 |
20100067634 | Furman et al. | Mar 2010 | A1 |
20120157005 | Bardin et al. | Jun 2012 | A1 |
Number | Date | Country |
---|---|---|
03021826 | Mar 2003 | WO |
Entry |
---|
European Search Report for Application No. 13190234.8 dated Nov. 28, 2013, 7 pages. |
Number | Date | Country | |
---|---|---|---|
20140120857 A1 | May 2014 | US |