Monopulse processor digital correction circuit

Information

  • Patent Grant
  • 4994810
  • Patent Number
    4,994,810
  • Date Filed
    Monday, March 26, 1990
    34 years ago
  • Date Issued
    Tuesday, February 19, 1991
    33 years ago
Abstract
A method and correction circuit is described incorporating a monopulse receiver, injecting a first signal of known amplitude and arbitrary phase and a second signal of the same amplitude shifted 90.degree. and measuring the output voltages of the channel and using the measurements to generate four coefficients which may be mathematically applied to the ouptut signals of the channel to provide a corrected output. The invention overcomes the problem of compensating for phase and gain drift in the in-phase and quadrature paths of the sum and difference channels of a monopulse receiver.
Description

BACKGROUND OF THE INVENTION
1. Field of the Invention
This invention relates to the calibration of analog circuitry in a monopulse receiver and more particularly to a digital correction circuit for compensating drifts and anomalies in phase and gain between quadrature channels.
2. Description of the Prior Art
A number of airports in the United States have parallel runways for the landing and take off of airplanes. When the weather is inclement so that the runways may not be visually observed by the pilot, one runway is shut down allowing only one runway to be used during periods of inclement weather.
A precision runway monitor may be provided to receive RF signals radiated by airplanes in response to a ground based air traffic control station. A monopulse radar receiver positioned to scan the area of airplanes flying an approach pattern to land on one or the other runway may be tracked from the signals radiated by the respective airplane. A requirement was established that the monopulse radar receiver should have an azimuth accuracy of 60 feet at a range of 60,000 feet which corresponds to 0.057.degree. accuracy. By using the resulting resolution, aircraft approaching the runway may be prevented from colliding into one another and assured of a landing on the correct runway. For example, the runways may be parallel to one another and spaced apart by 3,500 feet. The precision runway monitor may be placed on the ground midway between the two runways.
A typical monopulse receiver would have an antenna which simultaneously generates a sum and difference signal-in response to a received signal emanated by an aircraft. In order to process all the energy that is received indicative of location, the sum signal and difference signal are each processed by in phase and quadrature channels. It is well known that different gain and phase variations in the in-phase and quadrature channels result in an associated error in the sum signal and the difference signal which ultimately results in an error in the azimuth measurements of the incoming signals.
In a publication by F. E. Churchill et al. entitled "The Correction of I and Q Errors In A Coherent Processor", IEEE Transactions on Aerospace and Electronic Systems, Vol. AES-17, No. 1, January, 1981, pp. 131-137, a method is presented for correcting the I and Q channel errors by means of correction coefficients which are derived from measurements of a test signal. The required correction coefficients can be obtained from time samples of a test signal. By taking the discrete Fourier transformation (DFT) of the samples, three-digital filters are formed at dc, the test signal frequency, and the image of the test signal frequency. By using a simple algorithm, errors in gain may be inserted into one channel and errors in phase may be inserted into the other channel.
In U.S. Pat. No. 4,040,055 which issued on Aug. 2, 1977 to T. H. Donahue et al. entitled "Digital Compensator For Transmitter Instability", the in-phase and quadrature coherent phase detector outputs are continuously monitored for phase, amplitude and timing instabilities and arithmetic correction signals are generated and applied to the digitized received signal values to compensate for these instabilities.
In U.S. Pat. No. 4,021,804 which issued on May 3, 1977 to Dounce et al., the use of digital techniques is described for storage and for phase correction processing of the received signals in a coherent-on-receive pulse radar system. The phases of the received video signals are corrected by first converting the video signals to a series of complex digital numbers indicating the amplitude and phase of the received signals with respect to the stable local oscillator. The stored complex digital numbers are then operated on by a digital data processor to provide a new series of digitally corrected complex numbers representing coherent received video signals.
In a publication by A. I. Sinsky and P. C. P. Wang, entitled "Error Analysis of a Quadrature Coherent Detector Processor", IEEE Transactions On Aerospace and Electronic Systems, Vol. AES-10, No. 6, November, 1974, pp. 880-883, the output signal distortion resulting from relative phase and gain errors between two quadrature detector channels of a baseband processing system is computed.
SUMMARY OF THE INVENTION
In accordance with the present invention, a digital correction circuit and method is provided for closed loop digital compensation of phase and gain drift in a channel of a monopulse receiver comprising the steps of injecting a first signal of known amplitude and arbitrary phase into the channel, measuring the uncorrected in-phase and quadrature channel voltages to provide values I.sub.1 and Q.sub.1, inserting a 90.degree. phase shift in the first signal to provide a second signal, injecting the second signal into the channel, measuring the uncorrected in-phase and quadrature channel voltages to provide values I.sub.2 and Q.sub.2, computing four correction coefficients K.sub.1, K.sub.2, K.sub.3 and K.sub.4 from the values I.sub.1, Q.sub.1, I.sub.2 and Q.sub.2, and applying the coefficients to the output signals of the in-phase and quadrature channel voltages to provide third and fourth signals having compensation for phase and gain drift.
It is an object of the invention to increase the angle accuracy of a monopulse receiver by compensating for phase and gain drift of the sum beam in-phase and quadrature channels and the difference beam in phase and quadrature channels.
It is a further object of the invention to provide coefficients which may be multiplied with the output signals to provide values which may be added together to provide a compensated output signal.
It is a further object of the invention to perform computations of correction coefficients periodically which may be then averaged together to form new coefficients, for example, every minute.
It is a further object of the invention to use an analog to digital converter to convert the signals to digital values which may facilitate the speedy computation of the correction coefficients and of the compensated signals of the sum and difference beams.





BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a pictorial view of parallel runways with airplanes in an acceptable approach pattern.
FIGS. 2A and 2B assembled is a block diagram of one embodiment of the invention.
FIG. 3 is one embodiment of a correction circuit.
FIG. 4 is a second embodiment of a correction circuit.
FIG. 5 is a graph of the test signal.
FIG. 6 is a graph of the sum channel with an uncorrected output.
FIG. 7 is a graph of the sum channel with a corrected output.
FIG. 8 is a graph of the difference channel with an uncorrected output.
FIG. 9 is a graph of the difference channel with a corrected output.





DESCRIPTION OF THE PREFERRED EMBODIMENTS
Referring now to the drawings, FIG. 1 shows an airport 10 having parallel runways 12 and 13 positioned on a surface 14 of earth 15. Runways 12 and 13 are spaced apart by a predetermined distance, for example, 1,066.80 meters (3,500 feet) as shown by arrow 16. Airplanes or aircraft 18 and 19 may be on a landing approach along glide paths 20 and 21 to runways 12 and 13 respectively. The path of glide path 21 on earth surface 14 is shown by reference line 22.
An air traffic control signal may be radiated by antenna 29 towards airplanes 18 and 19. Airplanes 18 and 19 may both be carrying air traffic control transponders which upon receiving a predetermined interrogation signal by way of antenna 29 will radiate a return signal which may have encoded therein information to identify the specific aircraft emanating the reply and its altitude. Antenna 29 receives the reply signals from airplanes 18 and 19 and process them in a monopulse receiver to decode the information and determine the angle of their arrival from airplanes 18 and 19 respectively. A signal received from airplane 18 may be resolved within a volume shown by outline 28. The range may be determined by the delayed time from the original interrogation pulse from antenna 29 to the reception of the reply signal and subtracting therefrom a fixed delay that the transponder on the airplane adds prior to transmitting a reply pulse. Antenna 29 may be positioned midway between runways 12 and 13 and towards the end of the runway furthest from the approaching aircraft.
Referring now to FIGS. 2A and 2B, a block-diagram of a monopulse receiver 30 is shown. Reply signals to an air traffic control interrogation from aircraft 18 and 19 is shown by arrows 31 and 32 to antenna 29. Antenna 29 is configured to combine the signals to provide a sum channel over lead 34 and a difference channel over lead 35 to an input of RF amplifiers 38 and 39 respectively. Calibration signal generator 40 functions to provide a first signal of known amplitude and of arbitrary phase over lead 41. Lead 41 is coupled to terminal 42 of switch 43 and to one side of 90 degree phase shifter 44. The other side of phase shifter 44 is coupled over lead 45 to terminal 46 of switch 43. Switch 43 functions as a single pole double throw switch which functions to couple terminal 47 to terminal 42 or to terminal 46. Terminal 47 is coupled over lead 48 to an input of couplers 36 and 37. Calibration signal generator 40, 90.degree. phase shift 44 switch 43 and couplers 36 and 37 function to inject a first test signal of known amplitude and of arbitrary phase into the sum and difference channels 49 and 50 of monopulse receiver 30 during a first time. And at a second time to inject a second signal having a 90.degree. phase shift with respect to the first signal into sum channel 49 and difference channel 50.
The output of RF amplifier 38 is coupled over lead 51 to an input of sum channel detector 52 which may include, for example, two mixers. The output of RF amplifier 39 is coupled over lead 53 to an input of difference channel detector 54 which may include for example two mixers. Local oscillator 55 provides a first frequency signal over lead 56 to an input of sum channel detector 52 and difference channel detector 54. Local oscillator 55 also provides a second signal over lead 57 suitable for mixing. The frequency signal on lead 57 is shifted 90.degree. in phase with respect to the frequency signal on lead 56 which are at the same frequency. Sum channel detector 52 functions to provide an in-phase signal over lead 58 and a quadrature signal over lead 59 to respective inputs of low pass filters 60 and 61. The output of low pass filters 60 and 61 is coupled over leads 62 and 63 to respective inputs of analog to digital converters 64 and 65. Low pass filters 60 and 61 function to select the lower frequency (usually base band video) when mixing the signals on leads 56 and 57 with the frequency on lead 51. Analog to digital converters 64 and 65 function to convert by sampling the analog signals on leads 62 and 63 at a rate of at least twice the frequency of the information components contained on leads 62 and leads 63. The output of analog to digital converter 64 and 65 is coupled over leads 66 and 67 to respective inputs of correction circuit 68. Correction circuit 68 further functions to apply coefficients K.sub.1, K.sub.2, K.sub.3 and K.sub.4 to the output signals of the in-phase and quadrature signals on leads 66 and 67 to provide third and fourth signals over leads 71 and 72 respectively, having compensation for phase and gain drift in sum channel 49 Leads 71 and 72 are coupled to respective inputs of angle measurement 74.
Difference channel 50 has reference numerals with like reference numbers as sum channel 49 except for the numbers being primed. Difference channel 50 and its components function in like manner as the corresponding reference numerals in the sum channel 49. In-phase and quadrature signals on leads 71' and 72' respectively have been compensated for phase and gain drift by correction circuit 68'.
Circuitry 109 for periodic computation of the coefficients K.sub.1 through K.sub.8 for use by correction circuit 68 and 68' is shown in FIG. 2B. Calibration controller 110 provides a signal over lead 111 to calibration signal generator 40 shown in FIG. 2A. Calibration signal generator 40 provides test signals through couplers 36 and 37 to provide the corresponding test voltages on leads 66, 67, 66' and 67' which are coupled to data storage 112 and data storage 114. Calibration controller 110 provides control signals over leads 115 and 116 to data storage 112 and 114 respectively to store the signals during the time respective test signals are inserted by couplers 36 and 37. Calibration controller 110 also provides control signals over lead 117 to processor 118 which functions to compute the coefficients which are provided over leads 119 and 120 to correction circuits 68 and 68'. Samples taken during injection of test signals are coupled from data storage 112 over lead 124 to an input of processor 118 and from data storage 114 over lead 125 to an input of processor 118. Calibration signal generator 40 may provide a signal over lead 128 to solenoid 129 which functions to move switch 43 to its proper position.
Angle measurement 74 functions to determine the angle of the received signal based on the inputs on leads 71 and 72 from sum channel 49 and the inputs of 71' and 72' from difference channel 50. The angle measurement is provided over lead 93 to an input of runway selected for landing 76. An output from angle measurement 74 is passed on to range measurement 77 over lead 78 which may be for example the same signals on leads 71 and 72 from sum channel 49. Range measurements 77 functions to determine the range of the received signals and provides an output over lead 95 to an input of runway selected for landing 94. Lead 78 is coupled to an input of aircraft identification and altitude decoder 97 which functions to decode the signal to determine the identity and altitude of the aircraft. The output of aircraft identification and altitude decoder 97 is coupled over lead 98 to an input of runway selected for landing 94. Runway selected for landing 76 functions to determine from the range measurement and angle measurement the runway selected by the aircraft decoded by aircraft identification and altitude decoder 97. The output of runway selected for landing 94 is coupled over lead 99 to display 101 which functions to display the information. A second output of runway selected for landing 94 is coupled over lead 85 to an input of glide path detector 86 which functions to determine if the aircraft at the respective runway has the appropriate altitude (as determined by aircraft identification and altitude decoder 97), range (as determined by range measurement 7) and azimuth (as determined by angle measurement 4). Glide path detector 86 has an output over lead 87 to an input of display 101 to display appropriate information related to the glide paths of aircrafts 18 and 19 as they approach runways 12 and 13.
FIG. 3 shows one embodiment of correction circuit 68 Lead 66 is coupled to an input of multipliers 69 and 70. Multiplier 69 has a second input over lead 90 having a value representative of coefficient K.sub.1. Multiplier 70 has a second input over lead 91 having a value representative of coefficient K.sub.3. The output of multiplier 69 is coupled over lead 73 to an input of adder 74. The output of multiplier 70 is coupled over lead 75 to an input of adder 76. Lead 67 is coupled to an input of multipliers 79 and 80. A second input of multiplier 79 is coupled over lead 81 which has a value representative of coefficient K.sub.2. A second input of multiplier 80 has a input coupled over lead 82 representative of coefficient K.sub.4. The output of multiplier 79 is coupled over lead 83 to a second input of adder 74. The output of multiplier 80 is coupled over lead 84 to a second input of adder 76. The output of adder 74 is coupled to lead 71 and the output of adder 76 is coupled to lead 72.
FIG. 4 shows one embodiment of correction circuit 68' shown in FIG. 2. In FIG. 4, reference numerals having primes correspond to the function of like reference numerals without primes in FIG. 3. In FIG. 4, coefficient K.sub.5 is coupled to lead 90, lead 91 is coupled to coefficient K.sub.7, lead 81' is coupled to coefficient K.sub.6 and lead 82' is coupled to coefficient K.sub.8.
In operation, the signals on leads 66 and 67 as well as 66' and 67' are normally baseband digitized signals from sum channel 49 and difference channel 50 respectively. Correction circuits 68 and 68' shown in FIGS. 3 and 4 respectively, effectively remove the gain and phase errors from sum channel 49 and difference channel 50. Correction circuit 68 and 68' perform a linear transformation of the uncorrected input signals by use of multipliers 69, 70, 79 and 80 and adders 74 and 76 shown in FIG. 3. The corrected output quadrature components on leads 71 and 72 are related to the uncorrected input quadrature components on leads 66 and 67 by equations 1-4.
V'.sub..SIGMA.I =K.sub.1 V.sub..SIGMA.I +K.sub.2 V.sub..SIGMA.Q (1)
V'.sub..SIGMA.Q =K.sub.3 V.sub..SIGMA.I +K.sub.4 V.sub..SIGMA.Q (2)
V'.sub..DELTA.I =K.sub.5 V.sub..DELTA.I +K.sub.6 V.sub..DELTA.Q (3)
V'.sub..DELTA.Q =K.sub.7 V.sub..DELTA.I +K.sub.8 V.sub..DELTA.Q (4)
The transformation constants K.sub.1 through K.sub.8 used in correction circuits 68 and 68' are chosen so that the magnitudes of the corrected complex channel outputs are equal to G.sub.O V.sub.T and the phase of the sum and difference channel complex outputs are equal to each other. G.sub.O is the precise value of voltage gain referenced from the couplers 36 and 37 for sum channel 49 and difference channel 50 respectively. The value of G.sub.O is selected based on the dynamic range requirements of monopulse receiver 30. Equations 5-8 concisely define the corrected output signal characteristics when the test signal generated by calibration signal generator 40 and coupled by couplers 36 and 37 is V.sub.T expj.phi..
(V'.sub..SIGMA.I).sup.2 +(V'.sub..SIGMA.Q).sup.2 =(G.sub.O V.sub.T).sup.2 ( 5)
tan(.phi.-.phi..sub.o)=V'.sub..SIGMA.Q /V'.sub..SIGMA.I (6)
(V'.sub..DELTA.I).sup.2 +(V'.sub..DELTA.Q).sup.2 =(G.sub.O V.sub.T).sup.2 ( 7)
tan(.phi.-.phi..sub.o)=V'.sub..DELTA.Q /V'.sub..DELTA.I (8)
Note that .phi..sub.0 is the specific value of the test signal phase at couplers 36 and 37 at the time the values of the coefficients K.sub.1 -K.sub.8 are computed. If the K values then were held constant and the test signal phase was allowed to depart from .phi..sub.0, equations 6 and 8 would be satisfied. The values of the K's are such that when the test signal phase is .phi..sub.0, then the corrected outputs on lead 72 and 72' are 0 and lead 71 and 71' are equal to G.sub.O V.sub.T.
The procedure for establishing the correction coefficients K.sub.1 through K.sub.8 may be performed as described below. The input test signal switch 43 is set to terminal 42. The test signal amplitude must be precisely V.sub.T at couplers 36 and 37. The test signal phase is arbitrary and is assumed to be .phi..sub.0 at the time of calibration. The value .phi..sub.0 must remain constant during the period of time it takes to complete the measurements. The uncorrected voltages on leads 66, 67, 66' and 67' are recorded and designated V.sub..SIGMA.I (1), V.sub..SIGMA.Q (1), V.sub..DELTA.I (1) and V.sub..DELTA.Q (1).
The input test signal is then switched by switch 43 to terminal 46 which provides a 90.degree. phase advance path and the uncorrected signals for voltages on leads 66, 67, 66' and 67' are recorded and designated V.sub..SIGMA.I (2), V.sub..SIGMA.Q (2), V.sub..DELTA.I (2), and V.sub..DELTA.Q (2). The coefficients K.sub.1 through K.sub.8 are then computed. To ensure a reliable set of correction coefficients K.sub.1 through K.sub.8, an averaging procedure should be used. A number of sets of uncorrected data must be taken before computing the set of K's to be used for correcting the real world signal inputs. The computation of coefficients K.sub.1 through K.sub.8 is given by equations 9-18.
K.sub.1 =A[V.sub..SIGMA.Q (2)/D.sub..SIGMA. ] (9)
K.sub.2 =-A[V.sub..SIGMA.I (2)/D.sub..SIGMA. ] (10)
K.sub.3 =-A[V.sub..SIGMA.Q (1)/D.sub..SIGMA. ] (11)
K.sub.4 =A[V.sub..SIGMA.I (1)/D.sub..SIGMA. ] (12)
D.sub..SIGMA. =V.sub..SIGMA.I (1)V.sub..SIGMA.Q (2)-V.sub..SIGMA.Q (1)V.sub..SIGMA.I (2) (13)
K.sub.5 =A[V.sub..DELTA.Q (2)/D.sub..DELTA. ] (14)
K.sub.6 =-A[V.sub..DELTA.I (2)/D.sub..DELTA. ] (15)
K.sub.7 =-A[V.sub..DELTA.Q (1)/D.sub..DELTA. ] (16)
K.sub.8 =A[V.sub..DELTA.I (1)/D.sub..DELTA. ] (17)
D.sub..DELTA. =V.sub..DELTA.I (1)V.sub..DELTA.Q (2)-V.sub..DELTA.Q (1)V.sub..DELTA.I (2) (18)
By way of example, gain and phase errors have been selected to illustrate the calibration procedure. In Table I, G.sub..SIGMA. represents the gain common to both the in-phase and quadrature channel while G.sub..SIGMA.I represents the gain variation in the in phase channel and G.sub..SIGMA.Q represents the gain variation or error in the quadrature channel. The above distinction is also applicable to the phase designations in the sum and difference channel and the gains in the difference channel.
TABLE I______________________________________EXAMPLE CHANNEL PHASE & GAIN ERRORS.SIGMA. Channel .DELTA. Channel______________________________________G.sub..SIGMA. = 0.8 G.sub..DELTA. = 1.1G.sub..SIGMA.I = 0.9 G.sub..DELTA.I = 1.2G.sub..SIGMA.Q = 1.2 G.sub..DELTA.Q = 0.8.phi..sub..SIGMA. = 20.degree. .phi..sub..DELTA. = -10.degree..phi..sub..SIGMA.I = 10.degree. .phi..sub..DELTA.I = 20.degree..phi..sub..SIGMA.Q = -25.degree. .phi..sub..DELTA.Q = -15.degree.______________________________________
In Table I, the gain errors are voltage ratios and where the test voltage V.sub.T multiplied by G.sub.0 equals A which in turn equals 1. The signal phase during the calibration procedure will alternate between 40.degree. and 130.degree. at couplers 36 and 37.
Table II lists the analytical expressions for the uncorrected signals at the inputs to correction circuits 68 and 68'.
TABLE II______________________________________UNCORRECTED SIGNAL VOLTAGE EQUATIONS______________________________________ V.sub..SIGMA.I = AG.sub..SIGMA. G.sub..SIGMA.I cos(.phi. + .phi..sub..SIGMA. + .phi..sub..SIGMA.I) V.sub..SIGMA.Q = AG.sub..SIGMA. G.sub..SIGMA.Q sin(.phi. + .phi..sub..SIGMA. + .phi..sub..SIGMA.Q) V.sub..DELTA.I = AG.sub..DELTA. G.sub..DELTA.I cos(.phi. + .phi..sub..DELTA. + .phi..sub..DELTA.I) V.sub..DELTA.Q = AG.sub..DELTA. G.sub..DELTA.Q sin(.phi. + .phi..sub..DELTA. + .phi..sub..DELTA.Q) A = G.sub.o V.sub.T______________________________________
Table III lists the computed values of the uncorrected voltages at the initial phase condition,. .phi..sub.0, of the test signal from calibration signal generator 40 and then at initial phase plus 90.degree. by way of switch 43.
TABLE III______________________________________VALUES OF UNCORRECTED SIGNAL VOLTAGEAT .phi..sub.o and (.phi..sub.o + 90.degree.).phi..sub.o = 40.degree. .phi..sub.o = 130.degree.______________________________________V.sub..SIGMA.I (1) = 0.2463 V.sub..SIGMA.I (2) = -0.6766V.sub..SIGMA.Q (1) = 0.5506 V.sub..SIGMA.Q (2) = 0.7864V.sub..DELTA.I (1) = 0.8485 V.sub..DELTA.I (2) = -1.0112V.sub..DELTA.Q (1) = 0.2278 V.sub..DELTA.Q (2) = 0.8500______________________________________
Table IV lists the values of the eight correction coefficients computed from the values in Table III using equations (9)-(18).
TABLE IV______________________________________COMPUTED CORRECTION COEFFICIENTS______________________________________K.sub.1 = 1.388890 K.sub.5 = 0.893316K.sub.2 = 1.194950 K.sub.6 = 1.062690K.sub.3 = -0.972059 K.sub.7 = -0.239363K.sub.4 = 0.434928 K.sub.8 = 0.891704______________________________________
Table V lists the values of the uncorrected and corrected voltages as a function of the test signal phase (.phi.) when the correction coefficients are fixed at the values in Table IV. When (.phi.) equal 40.degree. the corrected signals are such that V'.sub..SIGMA.Q and V'.sub..DELTA.Q are 0.
TABLE V__________________________________________________________________________UNCORRECTED AND CORRECTED CHANNEL VOLTAGES AS AFUNCTION OF TEST SIGNAL PHASETESTSIGNALPHASE (DEG).phi. V.sub..SIGMA.I V.sub..SIGMA.Q V.sub..DELTA.I V.sub..DELTA.Q V.sub..SIGMA.I ' V.sub..SIGMA.Q ' V.sub..DELTA.I ' V.sub..DELTA.Q '__________________________________________________________________________0 0.6235 -0.0837 1.2999 -0.3719 0.7660 -0.6428 0.7660 -0.642810 0.5516 0.0837 1.2404 -0.2278 0.8660 -0.5000 0.8660 -0.500020 0.4628 0.2485 1.1432 -0.0767 0.9397 -0.3420 0.9397 -0.342030 0.3600 0.4057 1.0112 0.0767 0.9848 -0.1736 0.9848 -0.1736.phi..sub.o .fwdarw. 40 0.2463 0.5506 0.8485 0.2278 1.0000 0.0000 1.0000 0.000050 0.1250 0.6788 0.6600 0.3719 0.9848 0.1736 0.9848 0.173660 0.0000 0.7864 0.4515 0.5047 0.9397 0.3420 0.9397 0.342070 -0.1250 0.8701 0.2292 0.6223 0.8660 0.5000 0.8660 0.500080 -0.2463 0.9273 0.0000 0.7209 0.7660 0.6428 0.7660 0.642890 -0.3600 0.9563 -0.2292 0.7976 0.6428 0.7660 0.6428 0.7660100 -0.4628 0.9563 "0.4515 0.8500 0.5000 0.8660 0.5000 0.8660110 -0.5516 0.9273 -0.6600 0.8767 0.3420 0.9397 0.3420 0.9397120 -0.6235 0.8701 -0.8485 0.8767 0.1736 0.9848 0.1736 0.9848.phi..sub.o +90.degree..fwdarw. 130 -0.6766 0.7864 -1.0112 0.8500 0.0000 1.0000 0.0000 1.0000140 -0.7091 0.6788 -1.1432 0.7976 -0.1736 0.9848 -0.1736 0.9848150 -0.7200 0.5506 -1.2404 0.7209 -0.3420 0.9397 -0.3420 0.9397160 -0.7091 0.4057 -1.2999 0.6223 -0.5000 0.8660 -0.5000 0.8660170 -0.6766 0.2485 -1.3200 0.5047 -0.6428 0.7660 -0.6428 0.7660180 -0.6235 0.0837 -1.2999 0.3719 -0.7660 0.6428 -0.7660 0.6428190 -0.5516 -0.0837 -1.2404 0.2278 -0.8660 0.5000 -0.8660 0.5000200 -0.4628 -0.2485 -1.1432 0.0767 -0.9397 0.3420 -0.9397 0.3420210 -0.3600 -0.4057 -1.0112 -0.0767 -0.9848 0.1737 -0.9848 0.1737220 -0.2463 -0.5506 -0.8485 -0.2278 -1.0000 0.0000 -1.0000 0.0000230 -0.1250 -0.6788 -0.6600 -0.3719 -0.9848 -0.1736 -0.9848 -0.1736240 -0.0000 -0.7864 -0.4515 -0.5047 -0.9397 -0.3420 -0.9397 -0.3420250 0.1250 -0.8701 -0.2292 -0.6223 -0.8660 -0.5000 -0.8660 -0.5000260 0.2463 -0.9273 -0.0000 -0.7209 -0.7660 -0.6428 -0.7660 -0.6428270 0.3600 -0.9563 0.2292 -0.7975 -0.6428 -0.7660 -0.6428 -0.7660280 0.4628 -0.9563 0.4515 -0.8500 -0.5000 -0.8660 -0.5000 -0.8660290 0.5515 -0.9273 0.6600 -0.8767 -0.3420 -0.9397 -0.3420 -0.9397300 0.6235 -0.8701 0.8485 -0.8767 -0.1737 -0.9848 -0.1737 -0.9848310 0.6766 -0.7864 1.0112 -0.8500 -0.0000 -1.0000 -0.0000 -1.0000320 0.7091 -0.6788 1.1432 -0.7976 0.1736 -0.9848 0.1736 -0.9848330 0.7200 -0.5506 1.2404 -0.7209 0.3420 -0.9397 0.3420 -0.9397340 0.7091 -0.4057 1.2999 -0.6223 0.5000 -0.8660 0.5000 -0.8660350 0.6766 -0.2485 1.3200 -0.5048 0.6428 -0.7660 0.6428 -0.7660360 0.6235 -0.0837 1.2999 -0.3719 0.7660 -0.6428 0.7660 -0.6428__________________________________________________________________________
The magnitudes of the corrected signals are all equal and independent of .phi. and the phases of the corrected complex outputs are linearly related to the test signal phase. All corrected signal values obey equations 5-8. To further illustrate the relationships between the uncorrected and corrected output signals, reference is made to FIGS. 5-9, which graphs the complex values given in Table V. As expected, the corrected complex signals shown in FIGS. 7 and 9 are two perfect circles of equal radius and equal phase.
FIG. 5 is a graph of the test signal. In FIG. 5, the ordinate represents the voltage of the test signal quadrature component and the abscissa represents the voltage of the test signal in-phase component. Curve 90 shows the voltage of the test signal with an initial phase of 40.degree. shown by radius 91 and at an initial phase plus 90.degree. shown by radius 92.
FIG. 6 is a graph of the sum channel with an uncorrected output shown by curve 93. In FIG. 6 the ordinate represents V.sub..SIGMA.Q and the abscissa represents V.sub..SIGMA.I. The corresponding voltages with an initial test signal phase equal to 40.degree. is shown by radius 94 and with an additional 90.degree. is shown by radius 95.
FIG. 7 is a graph of the sum channel with a corrected output shown by curve 96. In FIG. 7 the ordinate represents V'.sub..SIGMA.Q and the abscissa represents V'.sub..SIGMA.I. The corrected voltages are now shown for an initial phase of 40.degree. of the test signal by radius 97 and now has a phase of 0.degree.. The Corrected voltages are shown for an initial phase of 130.degree. of the test signal by radius 98 with a phase of 90.degree..
FIG. 8 is a graph of the difference channel with an uncorrected output shown by curve 99. In FIG. 8 the ordinate represents V.sub..DELTA.Q and the abscissa represents V.sub..DELTA.I. The uncorrected output for an initial phase of 40.degree. of the test signal is shown by radius 100. The uncorrected output for an initial phase of 130.degree. is shown by radius 101.
FIG. 9 is a graph of the difference channel with a corrected output shown by curve 102. In FIG. 9 the ordinate represents V'.sub..DELTA.Q and the abscissa represents V'.sub..DELTA.I. The corrected output corresponding to the test signal with initial phase of 40.degree. is shown by radius 103 and the corrected output for an initial phase of 130.degree. is shown by radius 104.
An apparatus and method for closed loop digital compensation of phase and gain drift in the channels of a monopulse receiver has been described comprising the steps of first injecting a signal of known amplitude and arbitrary phase into the quadrature channel, second measuring the uncorrected in-phase and quadrature channel voltages to provide values I.sub.1 and Q.sub.1, third inserting a 90.degree. phase shift in the first signal to provide a second signal, fourth injecting said second signal into the quadrature channel, fifth measuring the uncorrected in-phase and quadrature channel voltages to provide values I.sub.2 and Q.sub.2, sixth computing for correction coefficients K.sub.1, K.sub.2, K.sub.3 and K.sub.4 from the values I.sub.1, Q.sub.1, I.sub.2 and Q.sub.2, and seventh applying the coefficients to the output signals of the in-phase and quadrature channel voltages to provide third and fourth signals having compensation for phase and gain drift.
Claims
  • 1. In a monopulse radar having sum .SIGMA. and difference .DELTA.channels for providing angular measurements of a radar target from the radar,
  • said sum channel including a first quadrature detector,
  • said difference channel including a second quadrature detector,
  • each of said first and second detectors having an in-phase I channel and a quadrature Q channel;
  • a method for digital compensation of phase and gain drift in said first and second detectors comprising the steps of:
  • first, injecting a first signal of known amplitude and arbitrary phase into both said first and second detectors;
  • second, measuring the uncorrected voltage outputs of said in-phase I and quadrature Q channels of said first detector and said second detector to provide values V.sub..SIGMA.I (1) from said in-phase channel of said first detector and V.sub..SIGMA.Q (1) from said quadrature channel of said first detector and V.sub..DELTA.I (1) from said in-phase channel of said second detector and V.sub..DELTA.Q (1) from said quadrature channel of said second detector;
  • third, inserting a 90.degree. phase shift in said first signal to provide a second signal;
  • fourth, injecting said second signal into both said first and said second detectors;
  • fifth, measuring the uncorrected voltage outputs of said in-phase I and quadrature Q channels of said first and second detectors to provide values V.sub..SIGMA.I (2) from said in-phase channel of said first detector and V.sub..SIGMA.Q (2) from quadrature channel of said first detector and V.sub..DELTA.I (2) from said in-phase channel of said second detector and V.sub..DELTA.Q (2) from said quadrature channel of said second detector;
  • sixth, computing directly from said values V.sub..SIGMA.I (1), V.sub..SIGMA.Q (1), V.sub..SIGMA.I (2) and V.sub..SIGMA.Q (2) a first set of four correction coefficients K.sub.1, K.sub.2, K.sub.3 and K.sub.4 ;
  • seventh, computing directly from said values V.sub..DELTA..sub.I (1), V.sub..DELTA.Q (1), V.sub..DELTA.I (2) and V.sub..DELTA.Q (2) a second set of four correction coefficients K.sub.5, K.sub.6, K.sub.7 and K.sub.8 ;
  • eighth, removing said first and second signals from both said first and second quadrature detectors;
  • ninth, applying said first set of coefficients to the output signals of said in-phase and quadrature channels of said first detector to provide third and fourth signals having compensation for phase and gain drift; and
  • tenth, applying said second set of coefficients to the output signals of said in-phase and quadrature channels of said second detector to provide fifth and sixth signals having compensation for phase and gain drift.
  • 2. The method of claim 1 wherein:
  • said sixth step of computing from said values V.sub..SIGMA.I (1), V.sub..SIGMA.Q (1), V.sub.93 I (2) and V.sub..SIGMA.Q (2) a first set of four correction coefficients K.sub.1, K.sub.2, K.sub.3 and K.sub.4 includes:
  • computing said first set of coefficients from the relationships:
  • K.sub.1 =A[V.sub..SIGMA.Q (2)/D.sub..SIGMA. ]
  • K.sub.2 =-A[V.sub..SIGMA.I (2)/D.sub..SIGMA. ]
  • K.sub.3 =-A[V.sub..SIGMA.Q (1)/D.sub..SIGMA. ]
  • K.sub.4 =A[V.sub..SIGMA.I (1)/D.sub..SIGMA. ]
  • where:
  • D.sub..SIGMA. =V.sub..SIGMA.I (1)V.sub..SIGMA.Q (2)-V.sub..SIGMA.Q (1)V.sub..SIGMA.I (2); and
  • A equals the desired corrected signal amplitude.
  • 3. The method of claim 1 wherein:
  • said seventh step of computing from said values V.sub..DELTA.I (1), V.sub..DELTA.Q (1), V.sub..DELTA.I (2) and V.sub..DELTA.Q (2) a second set of four correction coefficients K.sub.5, K.sub.6, K.sub.7 and K.sub.8 includes:
  • computing said second set of coefficients from the relationships:
  • K.sub.5 =A[V.sub..DELTA.Q (2)/D.sub..DELTA. ]
  • K.sub.6 =-A[V.sub..DELTA.I (2)/D.sub..DELTA. ]
  • K.sub.7 =-A[V.sub..DELTA.Q (1)/D.sub..DELTA. ]
  • K.sub.8 =A[V.sub..DELTA.I (1)/D.sub..DELTA. ]
  • where:
  • D.sub..DELTA. =V.sub..DELTA.I (1)V.sub..DELTA.Q (2)-V.sub..DELTA.Q (1)V.sub..DELTA.I (2); and
  • A equals the desired corrected signal amplitude.
  • 4. The method of claim 2 wherein:
  • said seventh step of computing from said values V.sub..DELTA.I (1), V.sub..DELTA.Q (1), V.sub..DELTA.I (2) and V.sub..DELTA.Q (2) a second set of four correction coefficients K.sub.5, K.sub.6,K.sub.7 and K.sub.8 includes:
  • computing said second set of coefficients from the relationships:
  • K.sub.5 =A[V.sub..DELTA.Q (2)/D.sub..DELTA. ]
  • K.sub.6 =-A[V.sub..DELTA.I (2)/D.sub..DELTA. ]
  • K.sub.7 =-A[V.sub..DELTA.Q (1)/D.sub..DELTA. ]
  • K.sub.8 =A[V.sub..DELTA.I (1)/D.sub..DELTA. ]
  • where:
  • D.sub.66 =V.sub..DELTA.I (1)V.sub..DELTA.Q (2)-V.DELTA.Q(1)V.sub..DELTA.I (2); and
  • A equals the desired corrected signal amplitude.
  • 5. The method of claim 4 wherein:
  • said ninth step of applying said first set of coefficients to the output signals of said in-phase and quadrature channels of said first detector includes:
  • applying said first set of coefficients according to the relationships:
  • V'.sub..SIGMA.I =K.sub.1 V.sub..SIGMA.I +K.sub.2 V.sub..SIGMA.Q, and
  • V'.sub..SIGMA.Q =K.sub.4 V.sub..SIGMA.Q +K.sub.3 V.sub..SIGMA.I
  • where:
  • V.sub..SIGMA.I and V.sub..SIGMA.Q are the output signals, respectively, of said in-phase and quadrature channels of said first detector; and
  • V'.sub..SIGMA.I and V'.sub..SIGMA.Q are said third and fourth signals, respectively; and wherein:
  • said tenth step of applying said second set of coefficients to the output signals of said in-phase and quadrature channels of said second detector includes:
  • applying said second set of coefficients according to the relationships:
  • V'.sub..DELTA.I =K.sub.5 V.sub..DELTA.I +K.sub.6 V.sub..DELTA.Q, and
  • V'.sub..DELTA.Q =K.sub.7 V.sub..DELTA.I +K.sub.8 V.sub..DELTA.Q
  • where:
  • V.sub..DELTA.I and V.sub..DELTA.Q and the output signals, respectively, of said in-phase and quadrature channels of said second detector, and
  • V'.sub..DELTA.I and V'.sub..DELTA.Q are said fifth and sixth signals, respectively.
  • 6. In a monopulse radar having sum .SIGMA. and difference .DELTA. channels for providing angular measurements of a radar target from the radar,
  • said sum channel including a first quadrature detector,
  • said difference channel including a second quadrature detector,
  • each of said first and second detectors having an in-phase I channel and a quadrature Q channel;
  • a method for digital compensation of phase and gain drift in said first and second detectors comprising the steps of:
  • first, injecting a first signal of known amplitude and arbitrary phase into both said first and second detectors;
  • second, measuring the uncorrected voltage outputs of said in-phase I and quadrature Q channels of both said first and said second detectors to provide values V.sub..SIGMA.I (1) from said in-phase channel of said first detector and V.sub..SIGMA.Q (1) from said quadrature channel of said first detector and V.sub..DELTA.I (1) from said in-phase channel of said second detector and V.sub..DELTA.Q (1) from said quadrature channel of said second detector;
  • third, inserting a 90.degree. phase shift in said first signal to provide a second signal;
  • fourth, injecting said second signal into both said first and said second detectors;
  • fifth, measuring the uncorrected voltage outputs of said in-phase I and quadrature Q channels of both said first and second detectors to provide values V.sub..SIGMA.I (2) and V.sub..SIGMA.Q (2) from said first detector and V.sub..DELTA.I (2) and V.sub..DELTA.Q (2) from said second detector;
  • sixth, computing directly from said values V.sub..SIGMA.I (1), V.sub..SIGMA.Q (1), V.sub..SIGMA.I (2) and V.sub..SIGMA.Q (2) a first set of four correction coefficients K.sub.1, K.sub.2, K.sub.3 and K.sub.4 ;
  • seventh, computing directly from said values V.sub..DELTA.I (1), V.sub..DELTA.Q (1), V.sub..DELTA.I (2) and V.sub..DELTA.Q (2) a second set of four correction coefficients K.sub.5, K.sub.6, K.sub.7 and K.sub.8 ;
  • eighth, repeating said first, second, third, fourth, fifth, sixth, and seventh steps a plurality of times to provide a plurality of values of said K.sub.1, K.sub.2, K.sub.3 and K.sub.4 of said first set of correction coefficients and K.sub.5, K.sub.6, K.sub.7 and K.sub.8 of said second set of correction coefficients;
  • ninth, computing from said plurality of values of K.sub.1, K.sub.2, K.sub.3 and K.sub.4 a third set of correction coefficients K.sub.1 av., K.sub.2 av., K.sub.3 av. and K.sub.4 av., each of said coefficients of said third set comprising the average value of corresponding coefficients of said plurality of values of said first set;
  • tenth, computing from said plurality of values of K.sub.5, K.sub.6, K.sub.7 and K.sub.8 a fourth set of correction coefficients K.sub.5 av., K.sub.6 av., K.sub.7 av. and K.sub.8 av., each of said coefficients of said fourth set comprising the average value of corresponding coefficients of said plurality of values of said second set;
  • eleventh, removing said first and second signals from both said first and said second detectors;
  • twelfth, applying said third set of correction coefficients to the output signals of said in-phase and quadrature channels of said first detector to provide third and fourth signals having compensation for phase gain and drift; and
  • thirteenth, applying said fourth set of correction coefficients to the output signals of said in-phase and quadrature channels of said second detector to provide fifth and sixth signals having compensation for phase gain and drift.
  • 7. The method as claimed in claim 6 wherein:
  • said sixth step of computing from said values V.sub..SIGMA.I (1), V.sub..SIGMA.Q (1), V.sub..SIGMA.I (2) and V.sub..SIGMA.Q (2) a first set of four correction coefficients K.sub.1, K.sub.2, K.sub.3 and K.sub.4 includes:
  • computing said first set of coefficients from the relationships:
  • K.sub.1 =A[V.sub..SIGMA.Q (2)/D.sub..SIGMA. ]
  • K.sub.2 =-A[V.sub..SIGMA.I (2)/D.sub..SIGMA. ]
  • K.sub.3 =-A[V.sub..SIGMA.Q (1)/D.sub..SIGMA. ]
  • K.sub.4 =A[V.sub..SIGMA.I (1)/D.sub..SIGMA. ]
  • where:
  • D.sub..SIGMA. V.sub..SIGMA.I (1)V.sub..SIGMA.Q (2)-V.sub..SIGMA.Q (1)V.sub..SIGMA.I (2); and
  • A equals the desired corrected signal amplitude;
  • and wherein:
  • said seventh step of computing from said values V.sub..DELTA.I (1), V.sub..DELTA.Q (1), V.sub..DELTA.I (2) and V.sub..DELTA.Q (2) a second set of four correction coefficients K.sub.5, K.sub.6, K.sub.7 and K.sub.8 includes:
  • computing said second set of coefficients from the relationships:
  • K.sub.5 =A[V.sub..DELTA.Q (2)/D.sub..DELTA. ]
  • K.sub.6 =-A[V.sub..DELTA.I (2)/D.sub..DELTA. ]
  • K.sub.7 =-A[V.sub..DELTA.Q (1)/D.sub..DELTA. ]
  • K.sub.8 =A[V.sub..DELTA.I (1)/D.sub..DELTA. ]
  • wherein:
  • D.sub..DELTA. =V.sub..DELTA.I (1)V.sub..DELTA.Q (2)-V.sub..DELTA.Q (1)V.sub..DELTA.I (2); and
  • A equals the desired corrected signal amplitude.
  • 8. The method of claim 7 wherein:
  • said twelfth step of applying said third set of correction signals to the outputs of said in-phase and quadrature channels of said first detector includes:
  • applying said third set of coefficients according to the relationships:
  • V'.sub..SIGMA..sub.I =K.sub.1 av.V.sub..SIGMA.I +K.sub.2 av.V.sub..SIGMA.Q, and
  • V'.sub..SIGMA.Q =K.sub.4 av.V.sub..SIGMA..sub.Q +K.sub.3 av.V.sub..SIGMA.I
  • where:
  • V.sub..SIGMA.I and V.sub..SIGMA.Q are the output signals, respectively, of said in-phase and quadrature channels of said first detector; and
  • V'.sub..SIGMA.I and V'.sub..SIGMA.Q are said third and fourth signals, respectively;
  • and wherein:
  • said thirteenth step of applying said fourth set of correction coefficients to the ptput signals of said in-phase and quadrature channels of said second detector includes:
  • applying said fourth set of coefficients according to the relationships:
  • V'.sub..SIGMA.I =K.sub.5 av.V.sub..SIGMA.I +K.sub.6 av.V.sub..SIGMA.Q, and
  • V'.sub..DELTA.Q =K.sub.7 av.V.sub..DELTA.I +K.sub.8 av.V.sub..DELTA.Q
  • where:
  • V.sub..DELTA.I and V.sub..DELTA.Q are the output signals, respectively, of said in-phase and quadrature channels of said second detector, and
  • V'.sub..DELTA.I and V'.sub..DELTA.Q are said fifth and sixth signals, respectively.
US Referenced Citations (6)
Number Name Date Kind
3950750 Churchill et al. Apr 1976
4003054 Goldstone Jan 1977
4021804 Dounce et al. May 1977
4040055 Donahue et al. Aug 1977
4484194 Arvidsson Nov 1984
4905010 Baechtiger et al. Feb 1990
Non-Patent Literature Citations (2)
Entry
F. E. Churchill et al., "The Correction of I and Q Errors in a Coherent Processor", IEEE, Transactions on Aerospace and Electronic Systems, vol. AES-17, No. 1, Jan. 1981, pp. 131-137.
A. I. Sinksy et al., "Error Analysis of a Quadrature Coherent Detector Processor", IEEE, Transactions on Aerospace and Electronic Systems, vol. AES-10, No. 6, Nov. 1974, pp. 880-883.