Note: Descriptions are shown in the official language in which they were submitted.
205~9~5
SYSTEN AND METHOD FOR COMPENSATION OF
IN-PHASE AND QUADRATURE PHASE AND GAIN TMR~T.~NCE
BACRGROUND OF THE INVENTION
Field of the Invention:
The present invention relates to radar signal
processing. More specifically, the present invention
relates to the calibration and correction of phase and
gain imbalances between the in-phase (I) and quadrature
(Q) outputs of a radar receiver channel.
While the present invention is described herein with
reference to illustrative embodiments for particular
applications, it should be understood that the invention
is not limited thereto. Those having ordinary skill in
the art and access to the teachings provided herein will
recognize additional modifications, applications, and
embodiments within the scope thereof and additional
fields in which the present invention would be of
significant utility.
DescriPtion of the Related Art:
Coherent radars measure target velocity along the
line-of-sight by determining the Doppler frequency shift
of the received radar return. Typically, an I/Q
synchronous detector is used to form I and Q components
of the received signal. These components are then
digitized and processed by a digital FFT (Fast Fourier
Transform) to form narrowband Doppler filters. If a
~'
2054995
single frequency signal is input into an ideal receiver
with no phase or gain errors, the signal will appear in
one filter of the Doppler filter output. If the I and Q
channels do not have exactly the same gain or if the
phase shift between the two channels is not exactly 90
degrees, a spurious signal, usually referred to as the
"image", will be generated at the output of the Doppler
signal processor. The magnitude of the spurious signal
is directly proportional to the gain and phase errors
between the two channels. The frequency of the spurious
signal is equal to the frequency of the true signal, but
with the opposite phase rotation.
In a digital Doppler filter bank of N filters,
numbered 0 through N-1, with a signal in filter k, the
image signal will appear in filter N-k-1. This spurious
signal may be classified in error as a real target or may
obscure real target returns. There is therefore a need
in the art for a technique for reducing the magnitude of
the spurious signal.
One technique for controlling the phase and gain
imbalance involves the building of the radar receiver
hardware in accordance with very tight specifications.
For example, those skilled in the art will appreciated
that to ensure that the image is 25 dB below the return
signal, the gain and phase errors must be held to about
0.5 dB and 5 degrees, respectively. In applications
where these tolerances must be maintained over a wide
temperature range and/or where extensive measurement and
trimming of circuit component values is not cost
effective, additional design complexity must be added to
the receiver design to ensure that the tolerances are
met. The smaller the requirement for the image, the
more receiver design complexity is required.
Thus, there is a need in the art for a simple
technique that would allow the phase and gain imbalance
20549!~S
to be measured and the error removed in real time. This
would allow the image requirement to be met using a
simplified and less expensive analog receiver design
built to looser tolerances.
8UMMARY OF THE INVENTION
The present invention is a system and method for
correcting gain and phase imbalance errors between I and
Q channels of a synchronous detector. In a most general
sense, the system includes a first circuit for detecting
and measuring an imbalance error, a method for computing
numerical coefficients to correct the error, and a
second circuit for applying the correction factors to
received signal data to eliminate the effects of the
imbalance error. More specifically, the invention
provides a method for correcting phase and gain imbalance
errors in a synchronous detector which includes the steps
of:
a) inputting a signal into the detector and
extracting therefrom a received signal and a
corresponding image signal;
b) inputting the received signal in a first Doppler
filter;
c) inputting the image signal in a second Doppler
filter;
d) forming a first discriminant g=Re~VTxVI/¦VT¦2};
e) forming a second discriminant h=Im~VTxVI/¦VT¦2~;
where:
VT ~ the complex (real and imaginary parts)
voltage output of a target at the output of the first
filter and
VI = the complex (real and imaginary parts)
20549~5
voltage output of an image of the target at the output
of the second filter;
f) computing phase (0) and gain (p) errors from the
first and second discriminants;
g) calculating first and second correction factors
(m and n) based on the phase and gain errors; and
h) using the correction factors to correct the
phase and gain imbalance errors.
Another aspect of this invention is as follows:
A method for correcting phase and gain imbalance
errors in synchronous detector including the steps of:
a) inputting a signal into said detector and
extracting therefrom a received signal and a
corresponding image signal;
b) inputting said received signal in a first
Doppler filter;
c) inputting said image signal in a second Doppler
filter
d) forming a first discriminate g=Re{VTgVI/ ¦VT¦2};
e) forming a second discriminant h=Im {VTXVI/ ¦VT ¦2};
wherein:
VT = the complex voltage output of a target at the
output of the first filter and
Vl = the complex voltage output of an image of the
target at the output of the second filter;
computing phase (0) and gain (p) errors from the
first and second discriminants according to the
relationships:
p = (l+q2+h2) - ~(1+q2+h2)2_4g2~l/2
2g
and
0 = sin-l[(g-p)/g (l_p2) ]112;
9 ~
4a
g) calculating first and second correction factors
based on said phase and gain errors according to the
relationships
m - (l+p)/[(l-p)cOs(20)]
h) using the correction factors to correct the phase and
gain imbalance errors.
BRIEF DE8CRIPTION OF T~E DRAWING8
Fig. 1 illustrates the receiver of the present
invention.
Fig. 2 shows an illustrative implementation of the
compensation circuit of the present invention.
DE8CRIPTION OF THE INVENTION
Illustrative embodiments and exemplary applications
will now be described with reference to the accompanying
drawings to disclose the advantageous teachings of the
present invention.
Fig. 1 illustrates the receiver lO of the present
invention including a synchronous detector 12 with a gain
imbalance of +p and -p and a phase imbalance of +~ and -
~in the I and Q channels respectively. The received
signal is of the form cos~(~+~)t+e] where w represents
the IF (intermediate frequency) and the LO (local
oscillator) frequency, ~ represents the Doppler
frequency, and e represents the initial phase of the IF
input signal. Quadrature mixers 14 and 16 serve to
2054995
create the in-phase and quadrature signal components of
the received signal by mixing the received signal with
cosine and minus sine functions of the local oscillator
frequency ~. A phase imbalance is represented by the ~
terms in the cosine and sine terms input into the
quadrature mixers 14 and 16. Amplifiers 18 and 20 serve
to amplify the quadrature signals. A gain imbalance is
represented~Ythe p gain terms shown in the amplifiers 18
and 20. The synchronous detector 12 is completed with I
and Q channel low pass filters 22 and 24.
In accordance with conventional teachings, the
outputs of the synchronous detector 12 are digitized by
analog-to-digital converters 26 and 28 and input to a
digital signal processor 32. The signal processor
performs complex Fast Fourier Transforms to convert the
time domain digital signal into frequency domain data.
In accordance with the present teachings, the
complex outputs of the FFT are sent to a general purpose
computer (34) and are used to generate two correction
factors. In the illustrative embodiment, the two
correction factors are used in a compensation circuit 30
to correct the digitized output of the receiver to remove
the errors induced by the I/Q phase and gain imbalance.
Fig. 2 shows an illustrative implementation of the
compensation circuit 30 of the present invention. The
compensation circuit 30 includes two digital multipliers
36 and 38 two data latches 40 and 42 and a subtractor 44.
The data latches 40 and 42 are loaded with the correction
factor coefficients received over a data bus and input
into the compensation circuit 30 from the general purpose
computer 34 through a computer interface 46. Hence, the
digitized I and Q signals received from the analog-to-
digital converters 26 and 28 are multiplied by the m and
n correction factors and subtracted from each other to
provide the compensated quadrature output Q'.
20S4995
Returning to Fig. 1, with no compensation applied
(i.e., m = 1 and n = 0) a signal is inserted into the
receiver that places the receive signal and the image
signals in separate Doppler filters. Two discriminants
are then formed:
g = Re{VTXVI/lvTl }
h = Im{VTxVI/IVTl2}
where:
VT = the complex output of the target at the FFT
output and
VI = the complex output of the image at the FFT
output.
From g and h, the phase (~) and gain (p) errors can be
computed. Knowing ~ and p, two correction factors, m and
n, can be determined. Finally, the value of Q is
corrected using the formula:
Q' = mQ - nI
where Q' is the corrected quadrature signal. The in-
phase (I) signal input into the compensation circuit is
unchanged and is passed directly to the output as I'.
The equations to perform the correction represent a
closed form solution and will effectively eliminate the
spurious image.
In operation, the I and Q outputs generated by the
synchronous detector 12 with gain and phase errors are:
20s4g9~
I = (l+p)cos[(w+~)t+e]cos(wt+~) [1]
= (l+p){cos[(2w+~)t+e+~]+cos[at+e-~]}/2 [2]
Q = -(l-p)cos[(~+~)t+e]sin(wt-~) [3]
5= -(l-p)(sin[(2w+~)t+e-~]-sin[~t+e+~]}/2 [4]
The 2w terms are filtered out by the low pass
filters 22 and 24 leaving:
10I = (l+p)cos[~t+e-~]/2 [5]
Q = (l-p)sin[~t+e+~]/2 [6]
Let the I channel be viewed as a reference, and passed
directly through the compensation circuit 30. Then:
I' = (l+p)cos[~t+e-~]/2
If:
Q' = (l+p)sin[~t+e-~]/2 [7]
then, the amplitudes of I' and Q' would be identical,
viz., equal gain of (1 + p) and the phase difference
would be restored to 90 degrees (quadrature phase of
cos[~t+e-~] and sin[at+e-~]). With the proper gain and
phase relationships restored, the unwanted spurious image
signal would be eliminated.
Derivinq 0' from I and Q:
From equation [6], regrouping Q, multiplying each
side by (l+p)/(l-p) and expanding:
[(l+p)/(l-p)]Q = [(l+p)(l-p)sin[(~t+e-~)+(2~)]]/[2(1-p)] [8]
205~995
= (l+p)[sin(at+e-~)cos2~+cos(at+e-~)sin2~]/2 [9]
= [(l+p)sin(at+e-~)]cos2~/2
+ [(l+p)cos(at+e-~)]sin2~/2 [10]
However,
[(l+p)[sin(~t+e-~)/2] = Q-
in accordance with equation [7] and
[(l+p)[cos(~t+e-~)/2] = I
in accordance with equation [5], so:
[(l+p)/(l-p)]Q = Q'cos(2~) + I sin(2~) [11]
Q' = (l+p)Q/[(l-p)cos(2~)] - tan(2~)I
or:
Q' = mQ - nI [13]
where:
m = (l+p)/[(l-P)cos(2~)] [14]
and
n = tan(2~) [15]
This illustrates that if it is possible to process
the uncorrected I and Q signals to determine the gain
20S~99S
g
and phase errors, p and ~, the two correction factors, m
and n, can be computed and applied to I and Q to form Q l,
a Q channel output rescaled and adjusted in phase so that
the gains of the two channels are equal (l+p) and the
phases of the two channels are separated by 90 degrees
(cos[~t+e-~] and sin[~t+e-~]).
Derivation of Correction Factor~:
In brief, by processing the I and Q signals through
a complex FFT, the target and image signals may be
separated into separate Doppler filters. A discriminant
D is formed by manipulating the outputs of the two
Doppler filters as follows:
D = VTxVI/¦VT12 = (IT+iQT)x(II+jQI)/lIT+jQTl [16]
where:
VT = the complex output of the Doppler filter
containing the target;
VI = the complex output of the Doppler filter
containing the image;
IT = the real component of VT;
QT = the imaginary component of VT;
II = the real component of VI; and
QI = the imaginary component of VI.
Hence,
D = (ITII-QTQI)+i(ITQI+QTII)/(IT+QT)
= g + jh [17]
205~995
- 10
where:
g = Re{D}
ITII QTQI
_ [18]
IT + Q2
and
lo h = Im~D~
ITQI+QTII
[19]
IT + QT
It will be shown below that:
P = ~[(l+g2+h2)]-[(1+g2+h2~-4g2]1/2}/2g [20]
and
~ = sin~l[(g-p)/g(l-p2)] [21]
Knowing the phase and gain errors, the correction
25 factors, m and n, can be found by substituting [20] and
[21] into [14] and tl5]. Then, the value of Q~ can be
computed from [13].
So by properly manipulating the outputs of the
Doppler filters containing the target and image signals,
correction factors may be computed which allow the target
image to be removed.
Evaluation of p and ~:
3 5 Expanding equations [5] and [6] as sums of functions
20~4995
11
of (~t+e) and ~:
I = [(l+p)/2]cos(~t+e-~) [5]
= [(l+p)/2][cos(~t+e)cos(-~)-sin(~t+e)sin(-~)]
= [(l+p)/2][cos(~t+e)cos~+sin(at+e)sin~] [22]
and
Q = [(1-p)/2]sin(~t+e+~) [6]
= [(l-p)/2][cos(~t+e)sin~+sin(~t+e)cos~] [23]
Therefore:
I+jQ = (1/2)[(1+p)cos(~t+e)cos~+(1+p)sin(~t+e)sin~]
+ (j/2)[(1-p)cos(~t+e)sin~+(l-p)sin(at+e)cos~] [24]
= (1/2)[cos(~t+e)cos~+pcos(~t+e)cos~+sin(~t+e)sin~
+ psin(~t+e)sin~+jcos(~t+e)sin~-jpcos(~t+e)sin~
+ jsin(~t+e)cos~-jpsin(~t+e)cos~] [25]
Regrouping the terms with a common factor of cos~ or
jsin~:
I + jQ = (1/2)cos~{[cos(~t+e)+jsin(~t+e)]
+p[cos(~t+e)-jsin(~t+e)]}
- + (j/2)sin~[cos(~t+e)-jsin(~t+e)]
-p[cos(~t+e)+jsin(~t+e)]} [26]
And finally factoring out the [cos(~t+e)+jsin(~t+e)] and
[cos(~t+e)-jsin(~t+e)] terms:
I + jQ = (1/2)[cos(~t+e)+jsin(~t+e)][cS~-jpsin~]
+ (l/2)[cos(~t+e)-jsin(~t+e)][pcos~+jsin~] [27]
Those familiar with the art will recognize the term of
2054995
12
the form A[cos(~t+e)+jsin(at+e)] as the complex form of a
sinusoidal wave of amplitude A, frequency ~/2~ and
initial phase e. Similarly, the term of the form
A[cos(~t+e)-jsin(~t+e)] is the complex form of a
sinusoidal wave of amplitude A, frequency -~/2~ (opposite
sense of rotation), and initial phase -e. These two
frequencies will be extracted by two distinct filters of
the FFT 32. The output response of each filter will be
proportional to the signal amplitude, A, and the phase
will equal the initial phase, e.
For the signal I+jQ of equation [27], the FFT 32
extracts the +~ and -~ frequencies, so the complex FFT
outputs are the target voltage, VT:
V(~) = VT = (1/2)eie{cos~-jpsin ~} [28]
and the image voltage, VI:
V(-~) = VI = (1/2)e~je{pcos~+jsin ~} [29]
where the exponential form, e+je is used to express the
initial phase, cos(e)+jsin(e).
Now, using equation [16] to form the discriminant D
yields:
D = VTXVI/IVTI [16]
= (eje/2) (e~je/2) ~ [pcos2~+psin2~]
+j[cos~sin~-p2cos~sin~]}/[(1/2)(cos2~+p2sin2~)1/2]2
[3o]
= [p+j(l-p2)cos~sin~]/(cos2~+p2sin2~) [31]
Let:
g=Re~D~=p/(cos2~+p2sin2~) [32]
and
- 20~499~
13
h = Im{D}=(l-p2)cos~sin~]/(cos2~+p2sin2~) [33]
Note that for small gain and phase errors, g~p and
h~cos~sin~~(1/2)sin(~/2).
From [32], solving for cos2~ and sin2~ in terms of g
and p:
coS2~,+p2sin2¢~ = p/g
10(l-sin2~)+p2sin2~=l-(l-p2)sin2~=p/g [35]
sin2~=(l-p/g)/(l-p2)=(g-p)/g(l p2) [36]
and
Cos2~+p2(l-cos2~)=p2+ (1-p2)cos2~=p/g [37]
cos2~ = ((p/g)_p2)/(l_p2) = (p_gp2)/g(1_p2) [38]
Substituting equations [36] and [38] into equation [33],
allows us to solve for the gain imbalance, p, in terms of
g and h, which were computed from the FFT outputs.
Hence:
h = (1-p2)cos~sin~/(cos2~+p2sin2~)
= (l-p2) ~ ~ (P-qP2)/g(l-p2) ~ r (g-p)/~(l p2) ] ,,1/2
[ (p_gp2)/g(l_p2) ~+[p2 (g-p)/g(l-p2) ]
30 ~1_p2~ (gp-p2-g2p2+gp3) 1/2 [40]
p_gp2+gp2 _p3
[gp3_ (1+g2) p2+gp] 1/2 [41]
- 14 20S499S
p2h2 = gp3_(1+g2)p2+gp [42]
gp2 _ (1+g2+h2)p+g = 0 [43]
P = (l+g2+h2)-+ r (1+g2+h2)2_4q211/2
2g
Since from equation [32], g ~ p, of the two roots for p,
the one of interest is:
p = (l+g2+h2)- r ( 1+g2+h2)2_4g211/2
2g
Rearranging equation [36] allows us to solve for the
phase imbalance, ~, in terms of the measured value, g,
and the computed value, p, from equation [45]. Hence:
~ = sin~l[(g-p)/g(1-p2)]1/2 [46]
This derivation shows that the phase and gain errors, ~
and p, can be computed from the discriminants g and h
which were computed from the output of the FFT filters 32
containing the target and image signals.
Computation of the Correction Factors m and n:
As shown in equations [14] and [15] above:
m = (l+p)/[(1-p)cos(2~)] [47]
and
n = tan(2~) [48]
20~4~95
_ 15
Thus, the correction factors m and n are computed by the
general purpose computer 34 and fed back to the
compensation circuit 30. The compensation circuit 30
adjusts one of the quadrature signals in accordance with
equation [13]:
Q' = mQ - nI [49]
If there is only a gain error, then
m = (1+p)/[(1-p)cOs(2~)]=(l+p)/(l-p) [49A]
15n = tan(2~) = 0 [49B]
For this special case (or if the effect of phase error is
ignored), compensation can be accomplished using only one
multiply since the coefficient n equals zero.
Calibration:
A target signal is placed in a Doppler filter such
that the target filter and image filter are separated.
The real and imaginary outputs of the FFT are output to
the general purpose computer and equations [32] and [33]
are implemented as follows:
g = Re~Target~xRe~Imagel-Im~Target~xIm~Image~ [50]
Re{Target}2 + Im{Target}
and
h = Re~Tarqet~xIm~Image~+Im~Target~xRe~Imaqe~ [51]
Re{Target}2 + Im{Target}2
20~995
16
The receiver 10 is adapted for use in a radar system
such as that of a missile (not shown).
Thus, the present invention has been described
herein with reference to a particular embodiment for a
particular application. Those having ordinary skill in
the art and access to the present teachings will
recognize additional modifications applications and
embodiments within the scope thereof.
It is therefore intended by the appended claims to
cover any and all such applications, modifications and
embodiments within the scope of the present invention.
Accordingly,
WHAT IS CLAIMED IS:
