US 5241313 A Abstract To allow measurement of the angle-of-arrival (AOA) of the radar pulses from an uncooperative ground-based emitter, the method exploits the time doppler shift resulting from the velocity of a high performance aircraft and the high maneuver ability available to such aircraft to form initial angle calculations. These initial angle calculations along with inherent radar stability are then used for subsequent AOA measurements. The aircraft is flown at a constant velocity V
_{B} along successive legs with an angle φ_{K} between successive legs, and the velocity V_{B} and the angles φ_{K} are found using a navigation system such as GPS or inertial navigation. The time difference T_{n} '-T_{1} ' is measured between the arrival of the first and last of n pulses for each sample. A general approximation equation of T_{n} '-T_{1} ' is normalized such that for any K samples the following equation applies,e wherein θ
_{K} is the angle from the line of flight for a leg K to the emitter, where e_{K} =(T'_{n} -T'_{1})/(n-1), s=V_{B} /c(n-1), c is the velocity of light, and R=e_{K} /(1-s cos θ_{2}). The value of the angle θ_{K} is calculated using the following equation. ##EQU1## where p=-e_{1} +e_{2} cosφr=(e
_{2} ^{2} -2e_{1} e_{2} cosφ+e_{2} ^{2})^{1/2}.Claims(1) 1. A method of measuring the angle of arrival (AOA) of radar pulses from a ground-based emitter via doppler shift, using apparatus including a pulse collection radar receiver, navigation means and computing means deployed on a high performance aircraft;
wherein the aircraft is flown at a constant velocity V _{B} along successive legs with an angle φ_{K} between successive legs, wherein the method includes determining the velocity V_{B} and the angles φ_{K} using the navigation means, using the radar receiver for receiving samples of radar pulses from said emitter, with each sample comprising n pulses, n being constant for all samples, and measuring the time difference T_{n} '-T_{1} ' between the arrival of the first and last of the n pulses for each sample;assuming that a general approximation equation of T _{n} '-T_{1} ' is normalized such that for any K samples the following equation applies,e wherein θ _{K} is the angle from the line of flight for a leg K to the emitter, and e_{K} and s are defined as ##EQU31## where c is the velocity of light; wherein the method of measuring angles of arrival θ_{K} comprises the steps:a) while the aircraft is flown in a straight line along a first leg K1 and the aircraft velocity V _{B} is determined using the navigation means, measuring a train of pulses 1 to n inclusive to determine the time interval T'_{n} -T'_{1} for that sample, calculating the value of e_{1} as ##EQU32## b) after the aircraft then executes a turn until it is flying in a straight path along leg K2, and the angle φ_{1} between legs K1 and K2 is determined using the navigation means, taking another sample by measuring a train of pulses 1 to n inclusive to determine a new time interval T'_{n} -T'_{1}, calculating the value of e_{2} as ##EQU33## c) calculating a value of θ_{2} using the following equation with the value of φ_{1} for φ. ##EQU34## where p=-e_{1} +e_{2} cosφr=(e _{2} ^{2} -2e_{1} e_{2} cos φ+e_{2} ^{2})^{1/2},d) calculating an initial value for R using the following equation, ##EQU35## e) after the aircraft then executes a turn until it is flying in a straight path along a leg K3, and the angle φ _{2} between legs K2 and K3 is determined using the navigation means, taking another sample by measuring a train of pulses 1 to n inclusive to determine a new time interval T'_{n} -T'_{1}, calculating the value of e_{3} as ##EQU36## f) calculating a value of θ_{3} using the following equation, with the value of φ_{2} for φ: ##EQU37## where p=-e_{2} +e_{3} cos φr=(e _{3} ^{2} -2e_{2} e_{3} cos φ+e_{3} ^{2})^{1/2},g) calculating another value for R using the equation ##EQU38## h) continuing the iteration repeating the above steps (e), (f) and (g) along successive legs until the value of R is refined to the point that it remains substantially constant from leg to leg, and calculating a value of the angle of arrival θ _{K} on any leg.Description The invention described herein may be manufactured and used by or for the Government of the United States for all governmental purposes without the payment of any royalty. The present invention relates generally to a method of measuring the angle of arrival (AOA) of radar pulses from a ground-based emitter via doppler shift. There are many methods employed to measure an uncooperative radar's AOA ranging from widely spaced receivers (i.e., multiple aircraft) to single aircraft approaches. Since this invention is for a single aircraft application the following discussion will be limited to only single aircfaft approaches. One method employs multiple antennas and measures angle via power difference among antennas. This method requires multiple antennas and multiple receivers and accuracies available are usually limited to plus or minus 15 degrees. Another method employs closely spaced multiple antenna elements where phase differences between an incoming pulse are measured and angles estimated. This method will allow accurate angle measurement (in the order of plus or minus 1 degree) however, it also requires a considerable amount of hardware for implementation. The following United States patents are of interest. U.S. Pat. No. 4,825,213--Smrek U.S. Pat. No. 4,746,924--Lightfoot U.S. Pat. No. 4,549,184--Boles et al U.S. Pat. No. 4,172,255--Barrick et al The above patents relate to apparatus and methods for detecting and tracking targets. In particular, the Smerk patent describes a radar system and technique for use in detecting and tracking targets. Return radar signals are doppler processed, phase shifted and compared in a manner which preserves the angle of arrival of a moving target irrespective of the boresight direction. Stationary target return signals are constructively combined so as to augment the target signal gain independent of the antenna boresight direction. The Lightfoot patent is directed to apparatus and methods for locating a target relative to a receiver aircraft utilizing electromagnetic emissions. The receiver aircraft includes a pair of ring tip synchronized scanning interferometer antennas for receiving non-reflected emissions from the target to determine the range from the receiver aircraft to the target. Other variables are calculated for locating and displaying the position of the target on a display in the receiver aircraft. To generate bearing information, the signals are received by a receiver aircraft lens antenna and split into a number of separate beams which are scanned at separate parts of the lens antenna in a predetermined manner. The bearing of the target is determined by comparing the amplitudes of the reflected signals received at lens ports. The Boles et al patent relates to a method and apparatus for moving target ordnance control. The azimuth angle of a moving target having a radial component of velocity relative to the electrical boresight of an interferometer antenna is determined by first measuring the radial component of velocity relative of the moving target. The observed net doppler frequency of the moving target is then shifted by a fixed quantity so that the moving target resides in a doppler cell whose frequency corresponds to that of clutter in the immediate vicinity of the moving target. By plotting observable data, the azimuth angle relative to the electrical boresight of the antenna can be obtained. The azimuth angle is then used for relative range and azimuth angle weapon guidance. The Barrick et al patent describes a system for radar remote sensing of near surface ocean currents in a coastal regions. The system comprises a pair of radar units to scatter signals from the shore off to the ocean waves. Underlying surface currents impart a slight change in velocity to the ocean waves, which is detected by the radar units. Each radar unit can determine the angular direction of arrival of the radar echo signals by comparing the phase of the signals received at three short antennas on the shore. Signals scattered from the same point on the ocean by each of the two geographically separated radar units are used to construct a complete current vector for that point. Although the prior art patents relate to tracking targets, they do not describe a method of exploiting the time doppler shift resulting from the relative velocity of the tracking system and the target to obtain a mathematical relationship to which an angle approximation technique is applied to obtain the angle of arrival of the radar pulses. An objective of the invention is to provide a method to be deployed on a high performance aircraft to allow measurement of the angle-of-arrival (AOA) of the radar pulses from an uncooperative ground-based emitter. The method according to the invention exploits the time doppler shift resulting from the velocity of the high performance aircraft and the high maneuver ability available to such aircraft to form initial angle calculations. These initial angle calculations along with inherent radar stability are then used for subsequent AOA measurements. ADVANTAGES AND NEW FEATURES: The potential exists for getting a relatively accurate angle measurement on a transmitter by exploiting the time doppler shift. This can be accomplished with less hardware than other approaches. FIG. 1 is a diagram showing the relationship of a stationary emitter and a receiver; FIG. 2 is a diagram showing the relationship of a stationary emitter, and a receiver in motion directly toward the emitter; FIG. 3 is a diagram showing the relationship of a stationary emitter, and a receiver in motion directly toward the emitter, with a sequence of pulses; FIG. 4 is a diagram showing the relationship of a stationary emitter, and a receiver in motion; FIG. 5 is a diagram showing the relationship of a stationary emitter, and a receiver in motion, with a sequence of pulses; FIG. 6 is a diagram showing an angle estimation geometry; FIG. 6a is a diagram showing an angle estimation geometry over a plurality of legs; and FIG. 7 is a bar graph showing count versus angle error. In order to understand this invention some familiarity with time doppler shift must be established. In the following analysis, the emitter is a ground based radar and the receiver is an airborne collector utilizing an omnidirectional antenna. TIME DOPPLER SHIFT: First consider the situation shown in FIG. 1, where an emitter is at point A, a receiver at point B, the distance between A and B is d In this case the time difference of arrival at the receiver is the same as the time difference of transmission at the transmitter. FIG. 2. Now consider relative motion between the emitter and receiver as shown in FIG. 2. Specifically consider the case where A is stationary and B is travelling toward A with constant velocity V In general, any pulse transmitted at time T FIG. 3. However, A and B are always moving toward each other and d d Using the relationship shown in equation (4), the time of arrival difference, T Under the previous assumptions and for V FIG. 4. Now consider the situation where the transmitter is stationary and the receiver is moving with a constant velocity, V In general, any pulse transmitted at time T when the transmitter and receiver are d FIG. 5. Now again consider a sequence of pulses where the first pulse is transmitted at time T Now assuming the approximations shown in equations (10a), (10b) and (10c), the time difference of arrival, T So there is a time shift occurring at the receiver dependent upon the receiver velocity and the relative angle between the transmitter and receiver. To get some idea of this time shift, let the emitter be transmitting a stable (constant rate) pulse train where the time of transmission difference between any two pulses is one msec. After 1001 pulses, the time difference between pulse number 1001 and the first pulse will be one second, and if there was no relative velocity between the transmitter and receiver the time difference of arrival at the receiver would also be one second. Now assume B is moving directly toward A with a constant (over the measurement time) velocity of 332 meters/sec (mach 1). The results of this calculation are shown in equation (12). ##EQU12## or there has been a 1.107 μsec shift in the received pulses relative to no movement or relative to the receiver flying normal to the emitter (i.e., θ Time Doppler Exploitation: As can be seen from the previous equations, the measured time shift is a function of the velocity of the aircraft and the angle of aircraft relative to the emitter. If the actual time T Angle Approximation Technique (FIG. 6): Assume the general approximation equation of T
e In equation (13), e One normalization technique for a stable emitter would be to require the receiver to collect the same number of pulses on each k. In this case equation (16) would result. ##EQU15## Equation (16) may be considered as a definition of the quantities e and R. Recall that T is the time of transmittal of pulses at the emitter at point A, and T' is the time that the pulses are received at the aircraft at point B. The quantity s is defined as ##EQU16## where V
e The above equation is in the form shown in equation (18) (Reference, "Table of Integrals and Other Mathematical Data", Herbert Dwight, page 79, equation 401.2). ##EQU17## Then using equation (18), equation (17) and solving for e Equation (19) can readily be solved for θ While the above equation was generated by assuming both e Iteration over several legs (FIG. 6a): The diagram of FIG. 6a shows the extension of the situation of FIG. 6 over a plurality of legs. The aircraft may fly along a first leg K1, make a turn and fly along a second leg K2, make another turn and fly along a third leg K3, and so on. The angle φ may have a positive or negative value, depending on the direction of a turn. The calculations are based on the assumption that the lines from points along the flight path to the emitter at point A are all substantially parallel. For example, if the aircraft were to fly along each leg for one second at mach 1, each leg would be 332 meters long, while the distance to the emitter may be of the order of 100 to 200 kilometers. For an initial calculation of a value for R, the angles θ In conclusion, one method of measuring angles of arrival assumes that the general approximation equation of T
e The method of measuring angles of arrival θ The aircraft then executes a turn until it is flying in a straight path along leg K2. The angle φ The value of θ r=(e An initial value for R can now be calculated using equation (20). ##EQU24## The aircraft then executes another turn until it is flying in a straight path along leg K3. The angle φ The value of θ r=(e Another value for R can now be calculated using equation (20). ##EQU27## The iteration then continues repeating the above steps along successive legs K4, etc. until the value of R is refined to the point that it remains substantially constant from leg to leg. The value of the angle of arrival θ Error Affects (FIG. 6): Up to this point the affects of transmission and measurement errors have been ignored. The following discussion addresses two types of potential errors. 1) Transmission jitter and 2) measurement quantization. The basic time of transmission error (commonly called jitter) can primarily be traced back to the transmitter not emitting a pulse exactly when directed to do so. Jitter is a function of the internal mechanics and design of the transmitter and varies for the various classes of transmitters and even among transmitter of the same class. Jitter is usually considered as ranging randomly between 0 and some maximum value. Jitter affects can easily be incorporated into the previous equations by changing the time of transmissions from T Equating equation (21), solving for T For a small j the above equation can be approximated as shown in equation (23). ##EQU29## Receiver time quantization error is a direct result of the measurement scheme of the receiver which in essence overlaps a continuous running clock of some period p with the detected incoming pulse. Typically a time measurement is recorded when an edge of a clock becomes active and a pulse has been detected. The time of arrival recorded then becomes the actual time of arrival t Some observations can now be made concerning the capability of our angle measurement scheme. First of all, the error affects of both jitter and quantization will become less as a longer time interval is measured or it may be possible to mitigate the affects by taking several short snap shots and averaging the result. Secondly, we typically have no control over the transmission jitter but on new systems random jitter is minimized. Thirdly, we do have control over our own measurement accuracy and reasonable state-of-the-art in this measurement can easily be less than 3 nsec which also minimizes this effect. Experimental Results: To test this invention a simulation was generated in which a receiver is flown in an environment consisting of one transmitter located approximately 140 km away initially. The receiver has some initial flight path which is randomly selected. Initially a collection of radar pulses is made and the receiver then jinked (defined as a change in angle in a short period of time) and another collection made. Jinks are randomly selected but not less than 10 degrees or larger than 45 degrees. During the jink time, data is not collected. The aircraft is flying at a constant velocity of mach 1 and the turn rate used is 5 degrees/sec. Each collection is made over a 5 sec period. Where data is collected for about 250 msec each sec and an average (over the 5 sec), collection value is used in each leg. A total of 5 jinks are collected and the angle and R value calculated on each jink. The R values for all jinks are averaged and this is used throughout the rest of the simulation to measure angle. The results of several runs are shown in FIG. 7. These are error counts on the y axis, and angle error value on the x axis. In one case, the transmitter has zero jitter and the receiver clock is 0.1 nanoseconds. As can be seen a total of 40 angle measurements were made (2 simulation runs) and the angle error ranged from 0.3 to 0.5 degrees. In the other, the emitter had 30 nanoseconds of jitter and the receiver clock period was 10 nanoseconds. In this FIG. 80 angle measurements were made (4 simulated runs) and the error was primarily centered around -1°±1°. It is understood that certain modifications to the invention as described may be made, as might occur to one with skill in the field of the invention, within the scope of the appended claims. Therefore, all embodiments contemplated hereunder which achieve the objects of the present invention have not been shown in complete detail. Other embodiments may be developed without departing from the scope of the appended claims. Patent Citations
Referenced by
Classifications
Legal Events
Rotate |