The present invention relates generally to the field of systems and methods for determining an angle-of-arrival for a detected signal. More particularly, the present invention relates to a system and method for resolving phase ambiguity of a long baseline phase interferometer using weighted unambiguous frequency differences in the detected signal.
Phase interferometry is the technique of using patterns of interference created by the superposition of two or more waves to diagnose the properties of the aforementioned waves. A phase interferometer may be used to provide precise angle-of-arrival information for radio frequency (RF) signals of interest. Phase interferometry has historical application in SIGINT and COMINT.
SIGINT is intelligence-gathering by interception of signals, whether between people (i.e., COMINT or communications intelligence) or between machines (i.e., ELINT or electronic intelligence), or mixtures of the two. Communications intelligence (COMINT) is a sub-category of signals intelligence that engages in dealing with messages or voice information derived from the interception of foreign communications. Electronic signals intelligence (ELINT) is intelligence gathering by use of electronic sensors.
The data that is available from SIGINT may reveal information about an opponent's defense network, especially the electronic parts such as radars, surface-to-air missile systems, aircraft, etc. SIGINT can be used to detect ships and aircraft by their radar and other electromagnetic radiation. SIGINT can be collected from ground stations near the opponent's territory, ships off their coast, aircraft near or in their airspace, or by satellite.
In addition to the content of signals, SIGINT provides information based on both type of electronic transmission and its originating location. Triangulation and more sophisticated radiolocation techniques, such as time of arrival methods, require multiple receiving points at different locations. These receivers send location-relevant information to a central point, or perhaps to a distributed system in which all participate, such that the information can be correlated and a location computed.
What is needed is a system and method for resolving the phase ambiguity created when using two or more antenna that are widely spaced relative to the wavelength of the signals being measured to provide an accurate direction finding measurement.
It would be desirable to provide a system and/or method that provides one or more of these or other advantageous features. Other features and advantages will be made apparent from the present specification. The teachings disclosed extend to those embodiments which fall within the scope of the appended claims, regardless of whether they accomplish one or more of the aforementioned needs.
One embodiment of the invention relates to a method of determining an angle-of-arrival for a detected signal. The method includes receiving the signal at a first antenna, receiving the signal at a second antenna, and generating an angle-of-arrival for the detected signal by resolving a phase ambiguity using frequency differences between the signal from the first antenna and the signal from the second antenna.
Another embodiment of the invention relates to a system for determining an angle-of-arrival for a detected signal. The system includes a first radio frequency receiver configured to receive the detected signal at a first antenna, a second radio frequency receiver configured to receive the detected signal at a second antenna, and an interferometer configured to generate an angle-of-arrival value for the detected signal by resolving a phase ambiguity using frequency differences between the signal from the first receiver and the signal from the second receiver.
Yet another embodiment of the invention relates to a method for determining an angle-of-arrival for a detected signal. The method includes the steps of receiving the signal at a first antenna, receiving the signal at a second antenna, generating a plurality of spectral lines based on a complex multiplication of a complex conjugation of the signal received at the first antenna and the signal received at the second antenna, generating a delta phase difference for pairs of the spectral lines using frequency differences, and generating an angle-of-arrival for the detected signal based on an average of the delta phase difference.
Alternative examples and other exemplary embodiments relate to other features and combinations of features as may be generally recited in the claims.
The invention will become more fully understood from the following detailed description, taken in conjunction with the accompanying drawings, wherein like reference numerals refer to like elements, in which:
Before describing in detail the particular improved system and method, it should be observed that the invention includes, but is not limited to, a novel structural combination of conventional data/signal processing components and communications hardware and software, and not in particular detailed configurations thereof. Accordingly, the structure, methods, functions, control, and arrangement of conventional components and circuits have, for the most part, been illustrated in the drawings by readily understandable block representations and schematic diagrams, in order not to obscure the disclosure with structural details which will be readily apparent to those skilled in the art, having the benefit of the description herein. Further, the invention is not limited to the particular embodiments depicted in the exemplary diagrams, but should be construed in accordance with the language in the claims.
Referring to
Since only two antenna and RF front-ends are needed, interferometer 110 uses fewer antennas and less hardware than traditional interferometers, reducing size, weight, power consumption, and cost. Further, interferometer 110 may be implemented in platforms that do not currently have direction finding capability. For example, many aircraft in use today only have top and bottom antenna.
Although shown and described above with reference to a single aircraft, interferometer 110 could also be implemented by multiple aircraft with single antennas forming a networked interferometer between them. Further, although shown according to a specific embodiment, interferometer 110 may alternatively be configured to include more, fewer, and/or different components configured to implement functions described herein.
Referring now to
Referring now to
In a step 302, signal processing begins by using a Fast Fourier Transform on the incoming signal received and both of antennas 120 and 130 and receivers 125 and 135 to generate the discrete Fourier transform representing the frequency domain of the signal received separately on each antenna. Although described and shown herein as being processed using Fast Fourier transforms, the incoming signal may alternatively be processed used a number of other methods including, but not limited to, delay discriminators, an acousto-optic modulator or Bragg cell, etc.
A first channel may be designated as a reference channel, shown as antenna 130 and receiver 135 in
Although antenna 120 and receiver 125 are shown and described herein as the difference channel and antenna 130 and receiver 135 as the reference channel, it should be understood that these designation may be reversed dependent on which receiver path contains the complex conjugation. The receiver path containing the complex conjugation is referred to as the reference channel.
The resultant discrete Fourier transform on the reference channel may be transformed by complex conjugation in a step 304 such that the signal from the reference channel is given by the equation:
SR=ARej(ωt) (2)
In a step 306, the resultant equation from both channels may be complex multiplied together to generate a magnitude, represented by the equation:
and an angle giving the phase, represented by the equation:
Where in the term ND in equations 3 and 4 represents the number of times the phase signal has passed through 2π when passing from the reference channel to the difference channel. This term, as used herein, represents the phase ambiguity.
In a step 308, following the complex multiplication of the reference and difference channels, interferometer 110 may be configured to threshold spectral lines. If Fast Fourier transforms are used to compute the frequency components on each channel, this threshold is applied to each of the frequency bins to eliminate noise and interference from the calculation and to eliminate spectral artifacts induced by the receivers' front end (as well as windowing and filtering in the digital domain).
Referring now to
Referring again to
Equation (5) is generated by expanding equation (4) for two frequencies at wavelengths of λ0 and λ1. The difference of the resultant two equations provides equation (5), referred to hereinafter as the delta phase differencing equation.
In the second instance, the signal wavefront hits the difference antenna 120 first. Referring now to
When the signal wavefront hits the reference antenna 130 first and the difference antenna 120 last, the difference equation (3), increases in the positive direction in phase as the difference (1/λ1−1/λ0) increases. When the signal wavefront hits the difference antenna 120 first the difference equation (3) increases in the negative direction in phase as the difference (1/λ1−1/λ0) increases. Determination of the signal arrival orientation may be done by computing the slope of the phase difference computed at two different frequencies.
Referring again to
Referring again to
Alternatively, by substituting λ with c/f, where c is the speed of light, the equation may be rewritten to provide the maximum unambiguous frequency difference as a function of the interferometer baseline, d, as:
According to another embodiment, the actual electromagnetic propagation velocity may be used in place of the speed of light.
If the spectral lines are not within the limit of a frequency constraint for ambiguity resolution in step 314, step 312 may be repeated to examine the next two spectral lines for first and second frequency. In the spectral lines are within the limit of a frequency constraint for ambiguity resolution, interferometer 110 may compute the delta phase difference using equation (5) in a step 316.
In a step 318, a determination is made whether N1−N2 is equal to 1. If so, in a step 320, Frequency 1, f1, is adjusted by 2π. If not, or following step 320, the average delta phase difference may be weighted by magnitude in a step 322. Once the spectral lines for a signal have been thresholded and selected, they can be processed by weighting the phase difference solution, equation (5) by the corresponding amplitude or signal-to-noise ratio. The leading to be accomplished using a sliding window constrained to f1−f2≦c/d (equation (9)) where the weighted phase difference is computed by the equation:
Following step 320, a determination is made in a step 322 whether all frequencies have been processed. If not, processing on the next frequency is performed beginning with step 312.
If all of the frequencies have been processed, an angle-of-arrival is computed in a step 324. Following the complex multiplication of the complex conjugated reference channel with the unconjugated difference channel, a real, Rs and an imaginary Is part is provided. Depending on the signal arrival orientation, as previously described, the phase may be computed as either tan−1(Is/Rs) or tan−1(Rs/Is). Once the phase difference is computed for the correct signal arrival orientation, the angle-of-arrival, α, may be computed using the equation:
The angle-of-arrival, α, applies to the angle the signal wavefront makes with the interferometer baseline as well as the angle-of-arrival relative to an interferometer boresight.
Referring to
While the detailed drawings, specific examples and particular formulations given describe preferred and exemplary embodiments, they serve the purpose of illustration only. The inventions disclosed are not limited to the specific forms shown. For example, the methods may be performed in any of a variety of sequence of steps. The hardware and software configurations shown and described may differ depending on the chosen performance characteristics and physical characteristics of the computing devices. For example, the type and number of antenna us, the method in which frequency difference is used to resolved the phase ambiguity, the position of the antenna, etc. may differ. The systems and methods depicted and described are not limited to the precise details and conditions disclosed. Furthermore, other substitutions, modifications, changes, and omissions may be made in the design, operating conditions, and arrangement of the exemplary embodiments without departing from the scope of the invention as expressed in the appended claims.
Number | Name | Date | Kind |
---|---|---|---|
5477230 | Tsui | Dec 1995 | A |
5708443 | Rose | Jan 1998 | A |
6700536 | Wiegand | Mar 2004 | B1 |
6791493 | Rose | Sep 2004 | B1 |