Background of the Invention
This invention relates to ground velocity detectors in general and more particularly to a ground velocity detector using the natural radiation from terrain.
It is important in aircraft navigation to determine the velocity of the aircraft with respect to the ground. Prior art devices for determining ground velocity include airborne Doppler navigation radar. However, this system involves the transmission of signals by the aircraft. Such transmissions make a military aircraft vulnerable to detection and destruction by enemy defensive systems. Similarly, the prior art includes velocity sensors which use the pattern of back-scattered power from a radiating antenna to determine ground velocity, as disclosed in U.S. Pat. Nos. 3,838,424 and 3,974,500. Like the Doppler system, this system also is vulnerable to defensive countermeasures because of its use of a transmitting antenna.
Prior art systems also include a passive sensor system which used antennas sited in the front and rear of an aircraft to sense natural terrain radiation. The relative time delay in the sensed signals measured the ratio of ground velocity to altitude. However, another method was needed to determine the aircraft's altitude. The use of active radar to measure altitude would make this system as vulnerable to countermeasures as the Doppler system.
Summary of the Invention
The present invention solves the problems in the prior art by using a passive system to sense ground velocity independent of altitude. Since no transmitting antenna is needed to sense ground velocity or altitude, this system does not betray the aircraft's location and make it vulnerable to attack.
The present invention generates a signal proportional to ground velocity by dividing the output frequency of a radiometric interferometer, which is directly proportional to ground velocity divided by altitude, by the phase difference between two radiometric interferometers, which is inversely proportional to altitude.
Each of the radiometric interferometers analyzes the natural radiation from terrain received by a pair of antennas which are separated by a distance such that their combined radiation pattern would generate a series of interference fringes along the terrain. The two pairs of antennas are spaced apart along the direction in which velocity is to be sensed. Variations in natural radiation from the terrain will cause the sum and difference of the outputs of each pair of antennas to vary periodically. The radiometric interferometer determines the frequency of these variations, which are proportional to the ratio of ground velocity to altitude. The phase difference between the two interferometers is inversely proportional to altitude. By dividing the phase difference by the frequency of the variations, a signal is generated which is proportional to ground velocity and independent of altitude.
By using two triads of antennas and switching between pairs of antennas, signals can be generated which are proportional to along-heading velocity and to drift velocity.
Brief Description of the Drawings
FIG. 1 is a schematic representation of the generation of interference fringes by the radiation patterns of the antennas in two radiometric interferometers.
FIG. 2 is a block diagram of the radiometric interferometer.
FIG. 3 illustrates the variation in power of the sum and difference of the signals received by one pair of antennas.
FIGS. 4a-c illustrate idealized waveforms at several points in the radiometric interferometer.
FIG. 5 illustrates the sensitivity of a pair of antennas to variations in natural radiation from terrain.
FIG. 6 illustrates the phase shift between signals received by the two radiometric interferometers in the passive microwave velocity sensor.
FIG. 7 is a block diagram of a passive microwave velocity sensor using two pairs of antennas to determine ground velocity.
FIG. 7a illustrates the relationship between FIGS. 7a-1 and 7a-2.
FIG. 7a-1 and 7a-2 show idealized waveforms at several points in the passive microwave velocity sensor of FIG. 7.
FIG. 8 shows an antenna configuration for use in a passive microwave velocity sensor that senses velocity in two directions.
FIG. 9 is a block diagram of a velocity sensor that senses velocity in two directions.
Detailed Description of the Invention
FIG. 1 schematically illustrates the generation of interference fringes by a pair of antennas 1 and 3 mounted in the front of an aircraft and a pair of antennas 5 and 7 mounted in the rear of the aircraft. Each pair of antennas is separated by a distance S. The aircraft is flying at height h. If antennas 1 and 3 were radiating at wavelength .lambda., they would generate a fringe pattern with a spacing W between positive maxima (or minima) of
if antennas 1 and 3 were fed in phase a maximum would exist at point A.
Instead of using antennas 1 and 3 as transmitting antennas, the present invention uses antennas 1 and 3 as receiving antennas. If a point P having radiometric passive emissions is located at a maximum of the fringe pattern, the signals received by antennas 1 and 3 will be in phase. On the other hand, if the point P is located at a minimum of the interference pattern, the signals received by antennas 1 and 3 will be 180.degree. out of phase. Thus as the aircraft flies over the terrain at velocity V, points P which have radiometric passive emissions would cause a frequency f to be received by the antennas, where
Equation (1) follows from the fact that the angular directional-sensitivity of antenna pair 1 and 3 is
where k is a constant and .theta. is the angle ASP in FIG. 1. Since .theta. = tan.sup.-1 (x/h), where x is the distance from point A to point P in FIG. 1, this equation can be rewritten as
The function g(x) equals zero when
where n is an odd integer. Thus g(x) is zero where x/h = tan (sin.sup.-1 n.lambda./2S). The separation W between successive nulls (e.g., n=1 and n=3) for small angles (.theta. less than 15 degrees) is nearly constant and is approximately equal to
At small angles .theta. .apprxeq. sin .theta. .apprxeq. tan .theta., and thus
fig. 2 is a block diagram of a radiometric interferometer designed to measure the frequency of passive radiometric emissions from terrain as the aircraft flies over the terrain. Signals received by antennas 1 and 3 are fed into magic Tee 9. The sum output 11 and the difference output 13 of magic Tee 9 are connected to terminals of a radio frequency single pole double throw switch 15. The switch 15 is driven by a 100 Hz. square wave generator 17 so that it alternately samples the sum output 11 and the difference output 13 of the magic Tee 9. The output of the switch 15 is mixed in a mixer 18 with the signal from a 35 GHz. local oscillator 19. The mixer output is fed into a 500 MHz video amplifier 21 and then into a square law detector 23 which tracks the envelope of the video carrier. Phase detector 25, which is driven by a commutating signal from square wave generator 17 removes the DC element from the signal and puts the sampled signal from the difference output 13 in phase with the sum output 11. The output 27 of the radiometric interferometer is a signal with a frequency proportional to the ratio of ground velocity to altitude.
FIG. 3 shows the sum and difference patterns for antennas 1 and 3. This figure indicates that the sum and difference patterns and mirror images of each other with respect to their mean value and thus the sum output 11 and the difference output 13 in FIG. 2 will be 180.degree. out of phase with each other.
FIGS. 4a-c show idealized waveforms at several points in the interferometer. FIG. 4a shows the waveform of the sum output 11 and the difference output 13 as sampled by switch 15. The dots above certain samples indicate that they were taken from the difference output 13.
FIG. 4b shows the output of the video amplifier 21. Finally, FIG. 4c shows the output 27 of the phase detector 25. The detector 25 alternately reverses the polarity of the signal in step with switch 15. Thus FIG. 4c shows the samples of the difference signal inverted in polarity. This subtracts out the DC value of the signal and causes the difference signal to be in phase with the sum signal, thereby cancelling unwanted common signals and noise, while enhancing the desired signal.
The effectiveness of the radiometric interferometer depends on variations in natural radiation from terrain. For example, assume that the interferometer of FIG. 2 is mounted on an aircraft traveling at an altitude of 1000 feet, the separation S between the centers of antennas 1 and 3 is 4.8" and the local oscillator frequency is 35 GHz. From equation (3), the separation between interference fringes is about 70 feet. Thus the interferometer of FIG. 2 will exhibit maximum responses to ground temperature variations with a wavelength of 70 feet. As the aircraft moves over the ground, temperature variations with this wavelength will move into phase with the receiving antennas and 180.degree. out of phase. Ground temperature variations with other wavelengths will not elicit as strong a response from the interferometer. FIG. 5 shows the mean squared response of the interferometer vs ground temperature variation length for a distance between interference fringes of 70 feet.
The radiometric interferometer of FIG. 2 generates a signal with a frequency proportional to the ratio of ground velocity to height. Thus if an aricraft flying at 1000 feet dropped to 500 feet while maintaining the same speed, the frequency of the radiometric interferometer output would double. This occurs because at an altitude of 500 feet the distance between interference maxima is half of the distance between maxima at an altitude of 1000 feet. Thus if a signal can be generated which is inversely proportional to distance between maxima (and hence inversely proportional to height) this signal can be divided into the interferometer output signal to produce a signal which is proportional to ground velocity and independent of height.
A signal proportional to height can be developed by using a second radiometric interferometer connected to a pair of antennas located aft of the antennas connected to the first interferometer. The antennas connected to the second interferometer are spaced apart a distance S which is equal to the distance between the antennas connected to the first interferometer. The second interferometer operates at the same local oscillator and video amplifier frequencies as the first interferometer.
FIG. 1 shows a pair of antennas 5 and 7 spaced apart a distance D from antennas 1 and 3. Since the spacing S between antennas 5 and 7 is the same as the spacing between antennas 1 and 3, the frequency received by antennas 5 and 7 is the same as the frequency received by antennas 1 and 3. However signals received by antennas 5 and 7 will undergo a phase delay of .chi., where
The phase delay is proportional to height but independent of ground velocity, and thus can be used in conjunction with the interferometer output frequency to determine ground velocity.
FIG. 6 illustrates the phase delay between the sum patterns where the antenna pairs are 15 feet apart and the interference stripes are 70 feet apart.
FIG. 7 is a block diagram of a passive microwave velocity sensor which uses two radiometric interferometers 29 and 31. FIGS. 7a-1 and 7-2 illustrate idealized waveforms at several points in the passive microwave velocity sensor of FIG. 7. Antennas 1 and 3 are connected to interferometer 29 and antennas 5 and 7 are connected to interferometer 31. The outputs of the interferometers 29 and 31 are connected to phase-locked loops (PLLs) 33 and 35 which convert the outputs to a square wave of the same frequency. Waveform 101 shows the output of the fore phase-locked loop 33 and waveform 103 shows the output of aft phase-locked loop 35. At the time t.sub.0 the aircraft is at 1000 feet altitude with a speed of 100 knots. At time t.sub.1 the aircraft is at 500 feet altitude with a speed of 100 knots. At time t.sub.1, the frequency of waveforms 101 and 103 has doubled (equation (2)) and the phase delay between 101 and 103 has also doubled (equation (4)). Inverter 37 inverts the output of the aft phase-locked loop 35 to produce waveform 105. The inverted output is combined in an AND gate 39 with the output of the fore phase-locked loop 33 to yield a pulse train 107 with a width proportional to the phase difference between the outputs of the radiometer interferometers 29 and 31. The average voltage of this pulse train is inversely proportional to height (equation (4)). The output pulses of AND gate 39 are inputs to an inverse servo 41. The inverse servo 41 increases the height of the pulses so that their average value equals a preset constant. The output of AND gate 39 is used to control an analog switch 43 which has as an input the output of integrator 45. A circuit 47 subtracts a comparison constant from the analog switch output. The output of circuit 47 is the input of integrator 45, which means that the output of integrator 45 will be driven to a voltage where the average output of circuit 47 is zero. Waveform 109 illustrates the output of analog switch 43 if the output of integrator 45 is 4 volts. Waveform 111 shows the output of circuit 47 if the comparison constant is 1 volt. When the aircraft is at 1000 feet, the average voltage of waveform 111 over a complete cycle (e.g., from t.sub.2 to t.sub.3) will be 0 when the integrator output value 113 is 4 volts. When the aircraft is at 500 feet the average voltage of waveform 111 over a complete cycle (e.g., from t.sub.4 to t.sub.5) will be 0 when the integrator output value 113 is 2 volts.
If the aircraft remained at 1000 feet and doubled its velocity, the frequency of the outputs of phase-locked loops 33 and 35 would double but the phase delay would remain the same. Waveform 101 shows the output of the fore phase-locked loop 33 and waveform 104 shows the output of the aft phase-locked loop 35. Waveform 104 is inverted to waveform 106 and combined with waveform 101 in AND gate 39 to yield waveform 108. When the output of analog switch 43 is waveform 110, then the output waveform 112 of circuit 47 will have an average value of 0. For the output of circuit 47 to have an average value of 0, the output value 115 of integrator 45 must be 4 volts. Thus a change in aircraft velocity does not effect the value of the output of integrator 45. The foregoing analysis shows that the output of integrator 45 is a DC voltage which is directly proportional to height and independent of velocity.
The output of integrator 45 is fed into frequency multiplier 49 which generates a signal with a frequency proportional to ground velocity and independent of height. Frequency multiplier 49 contains a pulse generator 51 which generates a series of fixed width, constant amplitude pulses 116 coincident with the ascending or descending transitions of the output 101 of phase-locked loop 33. Since the output of pulse generator 51 has the same frequency as the output of phase-locked loop 33, the mean value of output of pulse generator 51 over a full cycle is proportional to the ratio of ground velocity to height. For example, when aircraft velocity doubles or aircraft height is halved the output of pulse generator 51 would change from waveform 116 to waveform 117 and the mean value of the output would rise from 1/4 volt to 1/2 volt. Analog switch 53 multiplies the output of pulse generator 51 by the output of integrator 45. Multiplying integrator output value 113 by waveform 116 yields switch output waveform 118. If the aircraft height is halved, the integrator output value 113 becomes 2 volts. Multiplying this value by waveform 117 yields waveform 119. Waveform 119 has the same average value (1 volt) as waveform 118. On the other hand, if aircraft velocity doubles, the integrator output value 115 remains 4 volts. The resulting analog switch output waveform 120 has an average value of 2 volts. Thus the average value of the output of analog switch 53 is directly proportional to ground velocity and independent of height.
The next step is conversion of the output of analog switch 53 to an output waveform with a frequency proportional to ground velocity. Pulse generator 55 generates a pulse train of constant amplitude but variable frequency. Circuit 57 subtracts the output of pulse generator 55 from the output of analog switch 53. The output of circuit 57 is fed into integrator 59. The output of integrator 59 drives voltage controlled oscillator (VCO) 61. The output of VCO 61 sets the frequency of pulse generator 55. If the output of analog switch 53 is waveform 118 and the output of pulse generator 55 is waveform 121, then the input to integrator 59 will be 0, since waveform 118 has the same average value as waveform 121. Thus the output of the pulse generator 55 will remain at the same frequency. If the aircraft drops from 1000 feet to 500 feet, the output of analog switch 53 will be waveform 119. Since waveform 119 has the same average value (1 volt) as waveform 118, the output waveform 122 of the pulse generator will remain at the same frequency. Now assume the aircraft doubles its velocity while remaining at the same altitude. The output of analog switch 53 will be waveform 120, which has an average value of 2 volts. The output of integrator 59 will rise until the output frequency of pulse generator 55 doubles, and the resulting output waveform 123 of pulse generator 55 has an average value of 2 volts. Thus the output frequency of VCO 61 will be directly proportional to ground velocity and independent of height.
The velocity sensor can also generate an output with a polarity that indicates the direction of ground velocity. Such an output would be useful on helicopters which are capable of both forward and backward motion. Examination of the relationship between the leading edges of the pulse train output of AND gate 39 and the ascending transitions of the fore phase-locked loop 33 show that for forward motion the leading edges coincide with the ascending transitions while, for backward motion, the leading edges coincide with descending transitions. By combining the output of pulse generator 51 (which coincides with the ascending transitions) with the output of AND gate 39 in AND gate 63, an output is generated with a polarity that indicates the direction of travel of the helicopter.
The passive microwave velocity sensor shown in FIG. 7 has a bandwidth of 500 MHz. around the 35 GHz. center frequency. Mathmatical analysis shows that this bandwidth does not significantly affect the interference pattern near the center between each pair of antennas. Thus the existence of a finite bandwidth of 500 MHz. does not impair the operation of the radiometric interferometers.
The passive microwave velocity sensor which has just been described senses the along-heading velocity of the aircraft. Drift velocity, i.e., velocity perpendicular to the aircraft's heading, would not be detected by the velocity sensor. By setting up a separate velocity sensor along another axis, (e.g., along the axis of the wing of the aircraft), drift velocity can be sensed. However, the same effect can be achieved more economically by adding a third horn to each interferometer and time sharing the processing electronics.
FIG. 8 shows a configuration for three horn antennas for use in a velocity sensor to determine velocity in two directions. The configuration would be the same for both radiometric interferometers. The three horns may be used as two pairs in three different ways, depending on which horn is common to the two pairs; in this manner any drift angle in a range of 360.degree. may be covered. For the purpose of this description, the horn labelled B in FIG. 8 will be taken to be the common horn and the two pairs will be AB and BC. Assuming that the other interferometer has also been expanded as in FIG. 8, time sharing between pairs would be coordinated so that pairs AB would always be active simultaneously on one half cycle and pairs BC would be active on the other half. At the output of the processing electronics, height would not be affected by the time sharing but velocity, magnitude and direction would apply alternately to the AB and BC components.
Assume that the heading axis of the aircraft bisects angle ABC. Letting
V be the magnitude of the total velocity,
.delta. be the drift angle,
V.sub.1 be the component of velocity directed from A to B,
V.sub.2 be the component of velocity directed from C to B,
V.sub.h be the along heading component of velocity and
V.sub.d be the drift component of velocity, the following equations hold:
thus if velocity along the AB and CB axes is known, along-heading and drift velocities can be computed from equations (5) and (6). If another antenna is used as the common antenna, calculations are done by assuming that the angle bisector is the "heading". Equations (5) and (6) can be generalized to yield V.sub.H ', the component of velocity along the assumed "heading" and V.sub.D ', the drift velocity from the assumed "heading". The actual along-heading and drift velocity are computed from the following equations:
and
where .theta. is 120.degree., 0.degree., or -120.degree., depending on whether A, B, or C is the common horn.
The reason for switching from one common horn to another is that motion along a set of interference fringes must be avoided. For example, if horn B in FIG. 8 were the common horn and the drift angle was 60.degree., then a small, hot source on the ground would not cross any of the interference fringes associated with the BC horn pair. Thus, half of the time (when the BC pair was active) there would be no alternating signal to track and the entire scheme would break down. If the angular deviation of the velocity vector with respect to the angle bisector at the common horn in use were limited to about .+-.30.degree., the foregoing problem would never arise. This can be done by continuously computing the ratio V.sub.D '/V.sub.H '. Whenever .vertline.V.sub.D '/V.sub.H '.vertline. exceeded 5/8, it would signify that the generalized drift angle had exceeded .+-.32.degree. and a new common horn would be selected. The new common horn would be to the right of the previous common horn (e.g. from B to C) if the sign of the ratio were positive and it would move to the left (e.g. from B to A) if the sign were negative. After choosing the common horn, the appropriate value of .theta. would be selected for use in the coordinate rotation equations.
FIG. 9 is a block diagram of a passive microwave velocity sensor designed to implement the foregoing scheme. Many of the components, such as the radiometric interferometers 29 and 31 and the frequency multiplier 49 are the same as in the one-dimensional velocity sensor illustrated in FIG. 7. The input to the radiometric interferometers 29 and 31 is controlled by RF switching matrices 65 and 67. The matrices 65 and 67 designate a common horn and switch between the other two horns. Time sharing control 69 drives matrices 65 and 67 to switch between the other two horns at 25 Hz. The output 71 of frequency multiplier 49 represents the time-shared velocity components V.sub.1 and V.sub.2 along the axis connecting each pair of horns. Circuit 73 uses output 71 and the output of time sharing control 69 to implement equations (5) and (6) with respect to the assumed "heading". This implementation is done with standard analog circuitry, electronic switches, and sample and hold circuits. Sector selector and coordinate rotator 75 then implements equations (7) and (8) to generate outputs representing along-heading velocity V.sub.H and drift velocity V.sub.D. This implementation is also done with standard analog circuitry, electronic switches, and sample and hold circuits. As with the one dimensional velocity sensor, both of these velocity outputs are independent of height. Sector selector and coordinate rotator 75 also checks to see if .vertline.V.sub.D '/V.sub.H '.vertline. > 5/8 and, if so, sends a signal to the switching matrices 65 and 67 to switch common horns. Having chosen a new common horn, the sector selector and coordinate rotator 75 would also choose a new value of .theta. to use in implementing equations (7) and (8).
As previously indicated, the satisfactory operation of the radiometric interferometer depends on the existence of variations in ground temperature with wavelengths comparable to the separation between interference fringes in the combined radiation pattern of the antenna pairs. There are other limitations on the operation of the passive microwave velocity sensors shown in FIGS. 7 and 9. Aircraft flexure which would change the distance between horn pairs could impair the accuracy of the sensors. To accurately determine the center frequency of the radiometric interferometer output at moderate and low signal to noise ratios, the phase locked loop employs a smoothing time constant which is several times as large as the correlation time associated with the input signal spectrum. This introduces a delay in response time which is proportional to the ratio of height to velocity, and hence would be long for high altitude, low speed flight. On the other hand, the response time would be short for the low altitudes and high speed characteristic of "nap of the earth" flights. Furthermore, since navigational accuracy depends on ground velocity over a considerable distance (e.g. 100 nautical miles), the navigational accuracy of the system would be high. A final problem concerns ambiguities in altitude readings. In the velocity sensor the altitude is measured as a function of the phase difference between signals received by the radiometric interferometers. The phase difference function is single-valued only for a shift of less than 180.degree.. At low altitudes, the phase shift could exceed 180.degree., thus introducing an ambiguity into altitude determination. This effect can be ameliorated by mounting additional groups of antennas between the original groups, which would permit a decrease in the distance D between antenna groups, thus decreasing the magnitude of the phase difference. Alternatively, a lower frequency local oscillator (e.g., 18 GHz.) could be employed at low altitudes to increase the separation between interference fringes, which would also decrease the magnitude of the phase difference.
The velocity sensor shown in FIG. 7 and FIG. 9 may be altered by substituting dipole antennas for horn antennas, or changing the frequency of the local oscillator or video amplifier. These and other modifications may be made without departing from the spirit of the invention which is intended to be limited solely by the appended claims.