systems and methods are presented for passive location of transmitters in which two or more receivers time stamp received signals from target transmitters and the time stamped data for each target signal of interest is isolated to identify a peak power time of arrival for the signal at each transmitter from which differential scan observation values are derived, and for each signal of interest a line of position curve is computed based on the differential scan observation value and corresponding receiver locations, and for each signal of interest an estimated target transmitter location is determined based on an intersection of two corresponding line of position curves.
|
10. A system for locating a radar transmitter comprising:
a plurality of receivers operative to receive and time stamping signals from one or more target transmitters; and
a central processing facility coupled with the receivers and comprising a processor and a memory storing differential scan observation values derived from peak power time of arrival data from a pair of the receivers for each signal of interest, the processor being operative for each signal of interest to compute a line of position curve for each differential scan observation value based on the differential scan observation value and corresponding receiver locations, wherein the line of position curve comprises a function of a rotational speed of the radar transmitter, and for each signal of interest to determine an estimated target transmitter location based on an intersection of two corresponding line of position curves.
1. A method for passive location of a radar transmitter, the method comprising:
at each of a plurality of receivers, receiving and time stamping signals from one or more target transmitters;
for each receiver, isolating amplitude and time data for each signal of interest;
for each signal of interest for each receiver, identifying a peak power time of arrival;
for each signal of interest, merging peak power time of arrival data from multiple receivers:
for each signal of interest, computing a differential scan observation value for each pair of receivers with peak power time of arrival data for the signal of interest;
for each signal of interest, computing a line of position curve for each differential scan observation value based on the differential scan observation value and corresponding receiver locations, wherein the line of position curve comprises a function of a rotational speed of the radar transmitter; and
for each signal of interest, determining an estimated target transmitter location based on an intersection of two corresponding line of position curves.
2. The method of
where tb−ta is the differential scan observation value computed as the difference between peak power time of arrival values for receivers a and b, ω is the transmitter scanning frequency in radians per second ra is a multi-dimensional vector location of receiver a, rb is a multi-dimensional vector location of receiver b, and r is the line of position curve.
3. The method of
4. The method of
5. The method of
6. The method of
7. The method of
8. The method of
9. The method of
11. The system of
12. The system of
13. The system of
where tb−ta is the differential scan observation value computed as the difference between peak power time of arrival values for receivers a and b,ω is the transmitter scanning frequency in radians per second, ra is a multi-dimensional vector location of receiver a, rb is a multi-dimensional vector location of receiver b, and r is the line of position curve.
14. The system of
where tb−ta is the differential scan observation value computed as the difference between peak power time of arrival values for receivers a and b, ω is the transmitter scanning frequency in radians per second, ra is a multi-dimensional vector location of receiver a, rb is a multi-dimensional vector location of receiver b, and r is the line of position curve.
15. The system of
where tb−ta is the differential scan observation value computed as the difference between peak power time of arrival values for receivers a and b, ω is the transmitter scanning frequency in radians per second, ra is a multi-dimensional vector location of receiver a, rb is a multi-dimensional vector location of receiver b, and r is the line of position curve.
16. The system of
17. The system of
18. The system of
19. The system of
20. The system of
|
The present disclosure relates generally to locating the sources of radar transmissions and more particularly to methods and apparatus for passively locating radar emissions from transmitters.
A variety of applications exist for precisely locating the source or radio waves such as radar transmitters. Conventional location techniques employ signal arrival time-difference and phase-difference at spatially separated low gain or so-called main beam receivers and require simultaneous illumination of the receivers by the target radar signal. However, the receivers must have sufficiently small angular separation relative to the radar to achieve simultaneous illumination, which is not always practical. In addition, the geometric dilution of precision increases as the separation of the receivers decreases. Thus, there remains a need for improved techniques and apparatus for detecting the source of radar and other radio wave transmission.
Various details of the present disclosure are hereinafter summarized to facilitate a basic understanding, where this summary is not an extensive overview of the disclosure, and is intended neither to identify certain elements of the disclosure, nor to delineate the scope thereof. Rather, the primary purpose of this summary is to present some concepts of the disclosure in a simplified form prior to the more detailed description that is presented hereinafter. Systems and methods are presented hereinafter for passive location of radar transmitters which do not require simultaneous illumination of the receivers as was the case in conventional techniques. The system includes fixed or moving receivers that receive and time stamp signals from one or more target transmitters and a central processing facility with a processor and a memory storing differential scan observation values derived from peak power time of arrival data from a pair of the receivers for each signal of interest. For each signal of interest, the processor computes a line of position curve based on the differential scan observation value and corresponding receiver locations, and determines an estimated target transmitter location based on an intersection of two corresponding line of position curves for each signal of interest.
The receivers and/or the central processing facility in certain embodiments identify peak power time of arrival values for each signal of interest and each receiver, and merge peak power time of arrival data from multiple receivers for each signal of interest, where the central processing facility computes a differential scan observation value for each pair of receivers with peak power time of arrival data for the signal of interest. The central processing facility in one embodiment solves an equation for each differential scan observation value to compute the line of position curves, and may linearize the equation about an a priori location estimate and iteratively update the location estimate so as to minimize a root mean square (RMS) of the difference between actual and modeled differential scan observations. In certain embodiments, moreover, the central processing facility curve-fits received and time stamped signals using a second-order polynomial to identify the peak power arrival time at a maximum value of the fitted curve. In certain embodiments, three or more stationary receivers are used, and the location estimate is determined from the intersection of two or more line of position curves computed from two or more differential scan observation values. In other implementations, two or more receivers are used, with at least one receiver moving, where two differential scan observation values and corresponding line of position curves are computed to determine the estimated transmitter location.
Passive radar transmitter location methods are provided, which include receiving and time stamping signals from one or more target transmitters at each of a plurality of receivers, isolating amplitude and time data for each signal of interest and identifying a peak power time of arrival for each signal of interest for each receiver, as well as merging peak power time of arrival data from multiple receivers for each signal of interest. The methods further include computing a differential scan observation value for each pair of receivers with peak power time of arrival data and computing a line of position curve for each differential scan observation value based on the differential scan observation value and corresponding receiver locations for each signal of interest. In addition, the methods include determining an estimated target transmitter location based on an intersection of two corresponding line of position curves for each signal of interest. The line of position computation in certain embodiments includes solving an equation for each differential scan observation value and linearizing the equation about an a priori location estimate and iteratively updating the location estimate so as to minimize a root mean square of the difference between actual and modeled differential scan observations. The identification of the peak power time of arrival in some implementations includes curve-fitting the received and time stamped signals using a second-order polynomial and identifying a time value corresponding to a maximum value of the fitted curve.
The following description and drawings set forth certain illustrative implementations of the disclosure in detail, which are indicative of several exemplary ways in which the various principles of the disclosure may be carried out. The illustrated examples, however, are not exhaustive of the many possible embodiments of the disclosure. Other objects, advantages and novel features of the disclosure will be set forth in the following detailed description of the disclosure when considered in conjunction with the drawings, in which:
Referring to the figures, several embodiments or implementations are hereinafter described in conjunction with the drawings, where like reference numerals are used to refer to like elements throughout, and where the various features are not necessarily drawn to scale. Referring initially to
As best shown in
In the illustrated system 2, the receivers 6 receive signals emitted by the target transmitters, such as radar or other radio waves, and quantify the received signal strength and the time of signal receipt, such as by time stamping. Referring also to
The received signal power levels and arrival time (time stamp) information may be forwarded to the central processing facility 40 and/or the receivers 6 may perform one or more signal processing/analysis tasks prior to forwarding data to the central facility 40. In particular, the central processing facility 40 or the receivers 6 are configured to identify a peak power time of arrival value for each signal of interest and each receiver. As further shown in graphs 211, 214, and 217 in
Referring to
In operation, the central processing facility 40 merges the peak power time of arrival data from multiple receivers 6 and uses these compute a differential scan observation value 45 for each pair of receivers 6 with peak power time of arrival data for each signal of interest. In the illustrated embodiments, the differential scan observation values 45 are computed as the difference between corresponding peak power time of arrival values from pairs of receivers 6, for example, tb−ta computed as the difference between peak power time of arrival values for two receivers a and b.
In one embodiment, moreover, the central processing facility 40 solves the following equation (1) for each differential scan observation value 45 to compute a line of position curve 46:
where tb−ta is the differential scan observation value 45, ω is the transmitter scanning frequency in radians per second, ra is a multi-dimensional vector location of receiver a, rb is a multi-dimensional vector location of receiver b, and r is the line of position curve 46.
To improve the estimate of the target location, moreover, the central processing facility 40 in certain embodiments linearizes the equation (1) about an a priori location estimate and iteratively updates the estimated location so as to minimize a root mean square (RMS) of the difference between actual and modeled differential scan observations 45. Least squares data fitting (weighted least squares differential correction) may thus be used to identify a model estimated location for which the sum of squared residuals is minimal, where the residuals are computed as the difference between an observed value and the value given by the model.
The method 100 begins at 102 with receipt and time stamping of signals received at the plurality of receivers 6 from one or more target transmitters 10, 14. For each receiver 6, amplitude (e.g., power) and time data are isolated at 104 for each signal of interest. For each signal of interest and for each receiver 6, a peak power time of arrival is identified at 106, such as by curve-fitting the received and time stamped signals using a second-order polynomial and identifying a time value corresponding to a maximum value of the fitted curve (e.g., as illustrated and described above in connection with
The above examples are merely illustrative of several possible embodiments of various aspects of the present disclosure, where equivalent alterations and/or modifications will occur to others skilled in the art upon reading and understanding this specification and the annexed drawings. In particular regard to the various functions performed by the above described components (assemblies, devices, systems, circuits, and the like), the terms (including a reference to a “means”) used to describe such components are intended to correspond, unless otherwise indicated, to any component, such as hardware, software, or combinations thereof, which performs the specified function of the described component (i.e., that is functionally equivalent), even though not structurally equivalent to the disclosed structure which performs the function in the illustrated implementations of the disclosure. In addition, although a particular feature of the disclosure may have been illustrated and/or described with respect to only one of several implementations, such feature may be combined with one or more other features of the other implementations as may be desired and advantageous for any given or particular application. Also, to the extent that the terms “including”, “includes”, “having”, “has”, “with”, or variants thereof are used in the detailed description and/or in the claims, such terms are intended to be inclusive in a manner similar to the term “comprising”.
D'Agostino, Anthony, Middour, Jay W, Bynum, Kelia, Huffine, Christopher, Chrisman, Charles, Ellis, C Lane, Nichols, Randolph L.
Patent | Priority | Assignee | Title |
8447252, | Jan 21 2009 | AVAGO TECHNOLOGIES INTERNATIONAL SALES PTE LIMITED | Adaptive channel scanning for detection and classification of RF signals |
8849213, | Jan 21 2009 | AVAGO TECHNOLOGIES INTERNATIONAL SALES PTE LIMITED | Integrated circuit for signal analysis |
Patent | Priority | Assignee | Title |
5179573, | Feb 13 1992 | Bae Systems Information and Electronic Systems Integration INC | Amplitude measurement of received pseudonoise sequence using digital correlation |
5890068, | Oct 03 1996 | SAMSUNG ELECTRONICS CO , LTD | Wireless location system |
7307915, | Mar 09 2006 | Gula Consulting Limited Liability Company | Seismic modem |
20070247368, | |||
20090303130, |
Executed on | Assignor | Assignee | Conveyance | Frame | Reel | Doc |
Apr 30 2009 | MIDDOUR, JAY W | USA AS REPRESENTED BY THE SECRETARY OF THE NAVY, THE | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 022943 | /0858 | |
Apr 30 2009 | ELLIS, C LANE | USA AS REPRESENTED BY THE SECRETARY OF THE NAVY, THE | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 022943 | /0858 | |
Apr 30 2009 | CHRISMAN, CHARLES | USA AS REPRESENTED BY THE SECRETARY OF THE NAVY, THE | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 022943 | /0858 | |
May 01 2009 | BYNUM, KELLA | USA AS REPRESENTED BY THE SECRETARY OF THE NAVY, THE | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 022943 | /0858 | |
May 01 2009 | HUFFINE, CHRISTOPHER | USA AS REPRESENTED BY THE SECRETARY OF THE NAVY, THE | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 022943 | /0858 | |
May 04 2009 | AGOSTINO, ANTHONY D | USA AS REPRESENTED BY THE SECRETARY OF THE NAVY, THE | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 022943 | /0858 | |
May 04 2009 | NICHOLS, RANDOLPH L | USA AS REPRESENTED BY THE SECRETARY OF THE NAVY, THE | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 022943 | /0858 | |
Jun 09 2009 | The United States of America as represented by the Secretary of the Navy | (assignment on the face of the patent) | / |
Date | Maintenance Fee Events |
Sep 22 2014 | M1551: Payment of Maintenance Fee, 4th Year, Large Entity. |
Sep 18 2018 | M1552: Payment of Maintenance Fee, 8th Year, Large Entity. |
Jan 16 2023 | REM: Maintenance Fee Reminder Mailed. |
Jul 03 2023 | EXP: Patent Expired for Failure to Pay Maintenance Fees. |
Date | Maintenance Schedule |
May 31 2014 | 4 years fee payment window open |
Dec 01 2014 | 6 months grace period start (w surcharge) |
May 31 2015 | patent expiry (for year 4) |
May 31 2017 | 2 years to revive unintentionally abandoned end. (for year 4) |
May 31 2018 | 8 years fee payment window open |
Dec 01 2018 | 6 months grace period start (w surcharge) |
May 31 2019 | patent expiry (for year 8) |
May 31 2021 | 2 years to revive unintentionally abandoned end. (for year 8) |
May 31 2022 | 12 years fee payment window open |
Dec 01 2022 | 6 months grace period start (w surcharge) |
May 31 2023 | patent expiry (for year 12) |
May 31 2025 | 2 years to revive unintentionally abandoned end. (for year 12) |