The present invention relates to a system and method for RF Emitter Location with automatic UAS array reconfiguration.
Prior art methods for emitter location are based upon fixed array geometries within a system such as the fixed location of an antenna array mounted on an aircraft fuselage or at fixed terrestrial locations. This provides significant disadvantages. Consequently, there is a need for improving the placement and accuracy of emitter location.
The novel features believed characteristic of the invention are set forth in the appended claims. The invention itself, however, as well as a preferred mode of use, further objectives and advantages thereof, will be best understood by reference to the following detailed description of illustrative embodiments when read in conjunction with the accompanying drawings, wherein:
Several embodiments of Applicant's invention will now be described with reference to the drawings. Unless otherwise noted, like elements will be identified by identical numbers throughout all figures. The invention illustratively disclosed herein suitably may be practiced in the absence of any element which is not specifically disclosed herein.
Prior art methods for array-based emitter location are based upon fixed array geometries within a system such as the fixed location of an antenna array mounted on an aircraft fuselage or at fixed terrestrial locations. For mobile systems, such as a fixed array mounted on an aircraft fuselage, those skilled in the art are aware that there is a degree of freedom with respect to relocating the sensor array during an emitter location mission, but that the relative locations of the array elements are fixed. While emitter location accuracy can be improved through the repositioning of the sensor array, the fixed relative positions of the sensor comprising the system results in limitations. These limitations can degrade the emitter location accuracy and can limit the operable field-of-view (FOV) of the system. For example, an array of sensors in fixed locations on an aircraft fuselage can severely limit the FOV with regard to the elevation angle of the host aircraft and a ground-based emitter while also imposing lesser but significant location accuracies in the azimuthal plane.
This disclosure, in some embodiments, concerns a system comprising one or more unmanned autonomous devices with mobility to provide estimates of the location of a radio frequency (RF) transmitter. Examples include unmanned aerial system (UAS) platforms, unmanned terrestrial-based systems (UTS) that are autonomously mobile, space-based or satellite systems (SBS), combinations of UAS, UTS or SBS platforms operable as a coordinated system, or the augmentation of any of these systems with fixed-location components. While a UAS platform is described below, this is for illustrative purposes only and should not be deemed limiting. Each system component comprises an RF receiver and antenna that is coupled to an onboard computing system, a means for communications among the system components as well as optionally to an external entity, and a local onboard navigation system for each component. Within each system component, the antenna/receiver serves as the RF sensor, the internal navigation unit provides the position of each component, the onboard computing system hosts the emitter location algorithms, and the communications means enables each component to exchange information among other system components and an external entity.
Geolocation of RF emitters external to the said system has applications in areas such as search and rescue, payload delivery, military applications, and others. In these applications, it is useful to provide an accurate estimate of the geolocation of an RF emitter whose actual location is unknown but for which one of more of the system components is capable of detecting the presence of the RF emissions from the external source. For example, a search and rescue mission might require the ability to obtain the geolocation estimate of an emergency transponder beacon.
There are several prior art methods for estimating the geolocation of RF emitters. These include approaches based upon received power density of the emitter, interferometry, TDOA, beamforming, cyclostationary signal processing, and hybrid approaches. Emitter location methods are known to generally provide enhanced emitter location accuracies in proportion to the number of sensors within an array. Emitter location accuracy is further known to depend upon several characteristics of the system including the relative spatial geometry and positioning of the array elements and the RF emitter of interest.
The disclosed system and method solves the problem of reduced emitter location accuracy due to fixed array geometry limitations by replacing each sensor within the fixed array with an independently mobile platform such as UAS or UTS. Furthermore, the system of UAS/UTS components comprises the ability to automatically adapt and reconfigure their relative locations with respect to one another offering the capability to dynamically change the array geometry. Through the capability of dynamic array reconfiguration, subsequent emitter location estimates can occur wherein the array geometry is automatically adapted to enhance subsequent emitter location estimates. The capability to dynamically change the array geometry through repositioning of the array, wherein the array comprises one or more of the said systems, results in increasing the operational FOV and enables accurate emitter location estimates to be acquired with fewer measurement observations.
Additionally, a related method solves the problem of reducing RF emitter location error with a single platform by utilizing signal power density and the precise 3D features of the building or structure the RF emitter is located in or in close proximity to. We refer to this method as Feature Aided Signal Tracking (FAST).
The increased operational FOV provided by this disclosure enables emitter location estimates to be obtained for relative array and emitter geometries that would otherwise be limited due to incompatible overall geometries. In particular three-dimensional geolocation estimates can be obtained since the dynamic array can be oriented in the vertical plane coincident with the emitter and not restricted to a small range of elevation angles. Additionally, limitations due to the antenna pattern peaks and nulls present on each dynamic array element can be overcome by causing each dynamic array element to physically rotate; effectively enabling antenna pattern peaks and nulls to be oriented in any desired direction with respect to the external emitter location.
Since the orientation of the array element with respect to the location of the emitter can be dynamically changed, increased emitter location accuracies can be achieved with fewer samples of the received emissions. Adaptively reconfiguring the array upon subsequent measurements enables an optimal (or estimated optimal) orientation of the array with respect to the position of the emitter. Thus, subsequent measurements increase in accuracy. This is particularly important if the emitter only transmits RF signals in relatively short time periods.
Finally, the use of the dynamic array can enable the geolocation of emitters that may be located in positions for which fixed array methods are infeasible. For example, an emitter that is located beneath a tree canopy or within a large building or subterranean mine can be geolocated with a dynamic array comprised of UAS whereas a fixed array mounted on an aircraft fuselage would be obstructed.
A method may be implemented and present on the onboard computer readable media of the one or more said systems that allows for improving the accuracy of the emitter locations, determining new positions of the one or more said systems, and to classify, predict, compute, or extract secondary information of interest.
Turning to the figures,
The ability to dynamically change the designation of relayed communication links is advantageous and it is generally the case that the overall number of relayed channels should be minimized to balance overall network bandwidth requirements; however, the use of channel relays is important to overcome issues due blockage and multipath or to reduce the power consumption of UAS onboard transmitters. Furthermore, the array also has the ability to dynamically choose or elect a new UAS to serve as the array controller. Thus, the array follower (UAS 111-113, as depicted), can become the new array controller 110, and the previous array controller 110, becomes an array follower. Enabling the designation of the array controller to dynamically change offers advantages that include optimizing the required intra-network SNR, minimizing required intra-network power usage, and minimizing the number of relay UASs.
Within the
Scenario 200 indicates that the emitter location estimates for emitters 130-132 increased in accuracy since the diamond-shaped symbols representing the location estimates are closer to the center of the power density heat maps. Additionally, a new emitter location estimate is computed for emitter 140 and is shown in scenario 200 as the diamond-shaped symbol in the power density heat map of 140.
The computation of the updated array geometry 200 as compared to the original array geometry 100 consisted in part of the array controller 110 receiving information from UAS 111-113. This information includes data regarding the emitter location measurements of each UAS, information concerning the network power usage, and optionally, information regarding potential new emitters in the environment that were not yet geolocated during the scenario 100 measurements. After array controller UAS 110 receives this data from the other UAS, it computes new positions for each array element, including itself, and transmits this information to each UAS in the array. After receiving the new positioning data, the UAS move to their updated positions in scenario 200. Next, the wireless network automatically reconfigures to maximize network overall SNR and to determine if relayed channels should be established. To further optimize the network re-configuration, current array controller 110 determines if another member of the array should serve as the array controller. If it is determined that a new array controller should be chosen, this designation is broadcast to each UAS in the array and the newly chosen array controller assumes the coordinating responsibilities for future dynamic array changes.
UAS mobility control circuitry 303 is coupled to actuators 301 as indicated by line 302. UAS mobility components are well-known to exist in generic UAS platforms and are responsible for enabling the UAS to reposition its physical location.
The navigation system 304 represents any of a plurality of choices. For example, 304 can represent a global positioning system (GPS) receiver coupled to a GPS antenna module 306 via line 305. Block diagram elements 304-306 may likewise comprise a specialized interface to receive navigation information from an external entity. In some embodiments, the navigation system 304 may not require a navigation antenna 306 or coupling line 305; such as the case where the navigation unit 306 is in the form of an initialized inertial navigation system (INS) using onboard gyroscopes. The navigation system comprising blocks 304-306 may represent one or more of the previously mentioned options including the use of customized navigation systems such as those that operate in a GPS-denied environment. In one embodiment the system includes an external source such as one which transmits global coordinate system information to said system. The global coordinate system, such as those provided by a GPS system, can be transmitted to an array controller, array follower, etc. In one embodiment the system converts local coordinates such as described herein to global coordinates. This allows the system to utilize local coordinates when helpful, but convert to global coordinates when necessary. Further, there are situations wherein the GPS provided global coordinates, or global coordinates from some other external system, are not available. The system would then utilize local coordinates. The provided said global coordinates can be for the array, or they can be for another location such as a corner of a building. The array can then move to the corner of the building to correlate its local coordinates with those received by the external system providing global coordinates. The system would then be able to convert global coordinates to local coordinates and vice versa. Thus, in one embodiment, the system needs only one global coordinate which it can match with a local coordinate. Once the system can convert global coordinates to local coordinates and vice versa, the system can convert coordinates for all components of the system including the array controller, array follower, emitters, etc.
Elements 314-316 comprise the UAS infrastructure to enable it to communicate with external entities such as other UAS comprising the dynamic array including the elected array controller UAS. Additional external communications are supported by elements 314-316 and may include communications with a ground station or human operator. The communications protocols supported by elements 314-316 may be standardized wireless protocols including one or more of the WiFi (IEEE 802.11) family, cellular telephone standards including 5G, BlueTooth, or other standards as would be known by one skilled in the art. Furthermore, non-standard communications protocols may be used to support communications elements 314-316.
In one embodiment, each UAS is equipped with an RF intercept receiver 319 and antenna 321 that are coupled via line 320. Although intercept receiver 319 is depicted as a separate component with respect to communications transceiver 314, in some embodiments the communications functions of 316 and the intercept functions of 319 may be shared. In some embodiments, intercept receiver 319 may be in the form of a software-defined radio (SDR). Advantages in using an SDR to fulfill the role of block 319 include flexibility with respect to tuning ranges, demodulation capabilities, received signal bandwidth selection, and others as are known to those skilled in the art. The function of components 319-321 are to detect the presence of RF emissions in the ambient electromagnetic (EM) environment to support the geolocation capability of the dynamic array. Electromagnetic emissions comprise any emissions which contain electromagnetic energy and can include communications, data, etc. Optional features of intercept receiver 319 include the capability to demodulate received EM emissions, to estimate the received emission SNR, and to perform other analyses of the received emissions. These optional features may be controlled through manual, automatic, or combination of manual/automatic means. Automatic detection of signal modulation type and received signal SNR can be accomplished as described in [HL+:20][HS+:21]. Alternatively, the intercept functions of receiver 319 may be in the form of a conventional fixed tunable or non-tunable receiver as would be apparent to those skilled in the art.
Antennas 306, 316 and 321 may be present as separate and distinct antennas, or alternatively one or more of these antennas can be shared with appropriate feed networks (not shown in 300) that comprise filtering, RF distribution, and other components as are known to those skilled in the art. Antennas 306, 316, and 321 may further comprise a single or a plurality of elements. For example, if a beamforming approach is implemented for RF emitter location, intercept antenna 321 will comprise a fixed array of elements onboard the UAS. Likewise, if an amplitude-based geolocation method is implemented, one or more elements comprising the intercept antenna 321 may be supported by a mechanical means to physically rotate 321 at some known RPM (revolutions per minute) angular velocity to enable the intercept antenna pattern to rotate. Likewise, the RF feed network that interfaces the intercept antenna element 321 with the intercept receiver 319 (not shown in
The on-board computing system comprises elements 307-310. In addition to hosting computations to support typical and well-known functionality in support of a UAS, such as interfacing to the mobility control circuitry when the UAS is moving or stationary, additional functionality is present in the system to support the geolocation mission. Any of a variety of emitter location methods are implemented as instructions present in the computer readable memory (CRM) 307. Well-known emitter location methods include amplitude-based, phase interferometric, beamforming, TDOA, beamforming, cyclostationary signal processing, and hybrid approaches. Computer 310 may optionally comprise specialized processing hardware to enhance the emitter location calculations including either commercial off-the-shelf (COTS) components or custom-designed hardware. One example of a COTS hardware component is the inclusion of a graphical processing unit (GPU) that is well-known to increase the performance of certain mathematical computations. An example of a customized component is the inclusion of a high-speed and low-jitter clock signal generation circuit that could be used to time stamp intercepted signal artifacts when a TDOA emitter location method is employed. Those skilled in the art will appreciate and understand the engineering tradeoffs that exist when designing the architecture of the UAS in its supporting role as a dynamic array element for a RF emitter location application. For example, if a cyclostationary emitter location algorithm is implemented, the inclusion of a specialized COTS or custom-designed circuit for computing the covariance matrix and its corresponding eigenvalues of an intercepted signal by an array of elements comprising block 321 would enable enhanced speed in performing emitter location; however, the overall power requirements and weight of the UAS would likewise increase. Alternatively, implementing the covariance matrix formulation and associated eigenvalue computations supporting a cyclostationary emitter location as instructions present in the CRM 307 would enable the UAS to have reduced weight and power requirements while also increasing the delay of processing each observation of an external RF emitter.
CRM 307 is used to store instructions retrieved and executed by computer 310 to enable the UAS to perform a number of tasks. These tasks include support for the generic operation of the UAS such as internal commands to other components and interacting with the mobility control circuitry. In support of the dynamic array geolocation mission, the CRM will contain instructions that implement one or more emitter location methods; either wholly or in part if additional supporting hardware is included as described previously. Additionally, CRM 307 will comprise instructions that enable the UAS to function as an array controller if it is elected to serve in that role. CRM 307 will also contain instructions that implement the functionality to interact with intercept receiver 319 including prioritizing the list of tuning frequencies and other settings such as bandwidth.
Optionally, the CRM may contain instructions to automatically predict received RF signal modulation type and SNR level using methods such as those described in [HL+:20][HS+:21] or others.
Optionally, CRM 307 may contain instructions that enable the determination of specific instances in the intercepted signal wherein time stamps should be applied to generate TOA values. The generation of TOA time stamps is important when the dynamic array for RF emitter location is implementing a method such as the TDOA method. Emitter location accuracy with respect to TDOA methods requires that TOA values be assigned to the same and similar samples of an intercepted RF emission to ensure that calculated TDOA values are as accurate as possible. The spatial separation due to different physical locations of the UAS serving as sensors causes the TOA of an intercepted signal instance to vary due to the finite speed of propagation of EM energy. Therefore, knowledge of the propagation speed, the distances between pairs of UAS, the location of pairs of UAS, and the TOA of intercepted signal samples among pairs of UAS allow the emitter location estimate to be computed. Errors in any of the above values due to noise, clock resolution, or inaccurate TOA time stamps affect the overall emitter location accuracy.
A variety of means can be used to increase the TOA accuracy. One of these means is to use a very high-speed clock generator to reduce error due to resolution limits that also has very low drift and jitter rates to decrease the bias and random error associated with the clock. Clock signal generation and its effects on TDOA-based emitter location accuracy is discussed in more detail in the following paragraphs. Another means is to increase the emitter location estimation accuracy is to account for small differences in UAS location due to hover drift. Yet another means to increase the emitter location estimation accuracy is to ensure that intercepted RF signals have time stamped TOAs applied to the same sample point of the signal as received by different UAS platforms. This latter error source can be minimized in a plurality of ways. One way is to enable each UAS to have knowledge of the signal structure of the received signal. For example, if the intercepted signal is of a form wherein a specific artifact is present, such as the presence of a synchronization or “start” bit in a wireless packet-based data communications signal, each UAS can simply search for the presence of the synchronization bit and provide time stamps to that portion of the signal. This approach requires that the UAS have a priori knowledge, or perhaps automatically recognized knowledge, pertaining to the structure and type of signal that is being received by 319. Another possible means of enhancing the accuracy of time stamping the received signal is to enable a pair of UAS to exchange a sequence of stored samples of a recently received signal. In this manner, each UAS may execute a time-series cross-correlation technique implemented as stored instructions in CRM 307. By maximizing the cross-correlation coefficients, the associated time lag would yield the TDOA directly. Yet another method is to enhance the correlation method or the recognition of signal artifacts with a machine learning (ML) approach such as a neural network or LSTM recurrent network to further enhance the accuracy of choosing stored signal samples representing the same instant in an intercepted signal as received by multiple UAS array elements.
Navigation system 304 may further be augmented with special purpose circuitry as needed to support the RF emitter location mission. As an example, if a TDOA based method is implemented for geolocation, block 304 can supply time stamping information in a periodic or on-demand fashion that is used by computer 310 to determine TOA and TDOA of intercepted RF emission artifacts. Some navigation systems, such as GPS, inherently provide accurate time information. Furthermore, the use of COTS or customized clock signal generators can be incorporated in a variety of different elements within block diagram 300 such as within the computer 310, intercept receiver 321, the communications transceiver 316, or even as a separate standalone element that is interfaced to any of these components.
The dynamic array system that is the subject of
In
When the above three conditions are satisfied, the LOCA emitter location calculation can proceed and is described here with additional details given in [Sch:72]. It is assumed that the ranges from the emitter to the UASs are sufficiently small such that curvature of the earth is negligible, thus all calculations are Euclidean and occur with the assumption that the range values fall along a straight line. Furthermore, it is assumed that the ranges from the UASs to the emitter are long enough that they fall within the far-field of the EM transmissions of the emitter. Far-field propagation is typically considered to be a distance in excess of 20 wavelengths in terms of the emitter's carrier frequency although values other than 20 may be assumed in some cases. Far-field assumptions enable the emissions to be approximated as having planar wavefronts. Furthermore, the effects of in-band signal interference are not considered as it is further assumed that the intercept receivers 319 have sufficiently narrow bandwidths as to reject any transmissions from other emitters that are closely spaced to the emitter of interest in terms of carrier frequency. The wavefront of the transmissions from the emitter of interest is assumed to have a constant propagation speed equal to the speed of light in a vacuum, c. Alternative values for EM phase velocity, c, can be used without loss of generality so long as it is a constant. Distances, or range values r, are easily computed by dividing the EM propagation speed c by a propagation delay value (in units of time).
In keeping with the notation of [Sch:72] UASs 110-112 are mathematically denoted as S1(x1, y1), S2(x2, y2) and S3(x3, y3) respectively and the emitter to be located is denoted as T(x0, y0). Thus, the localized planar coordinates of each station Si are likewise given as (x1, y1), (x2, y2) and (x3, y3) with respect to stations S1, S2 and S3. Likewise, the localized planar coordinates of the emitter to be located is (x0, y0). The inter-station ranges are shown as solid line segments 520, 521 and 522 and are denoted as Δ12, Δ13 and Δ23 where the subscript ij in the variable Δij is the range between UAS Si and UAS Δj in units of distance. Array controller 110 calculates the Δij values as the Euclidean distance among the three (x1, y1), (x2, y2) and (x3, y3) values. The range of each UAS to the arbitrarily chosen origin of the two-dimensional plane where the UASs and the emitter are coplanar is denoted as ai for UAS i. The array controller computes ai2, the square of the distance of UAS i to the origin as ai2=xi2+yi2. From results in analytic geometry, the linear equation for the conic axis in the two-dimensional plane can be computed using the quantities (x1, y1), (x2, y2) and (x3, y3) x1, y1, x2, y2, x3, y3, Δ12, Δ31, Δ23, a12, a22 and a32 as given in Equation (1).
(x1Δ23+x2Δ31+x3Δ12)x+(y1Δ23+y2Δ31+y3Δ12)x=½(Δ12Δ23Δ31+a12Δ23+a22Δ31+a32Δ12) (1)
Since Equation (1) is a line coincident with the conic axis, the emitter location (x0, y0) is a point on the line.
TDOA values. For the sake of notation, we will let TOAi represent the measured TOA for UAS i. Likewise, TDOAij represents the time difference of arrival of an instance of an emitted transmission among UASs i and j.
(x2, y2) and (x3, y3) that correspond to the positions of the three UAS 110-112 the three TDOA values into ranges among the three UASs 110-112 by dividing the propagation speed, typically the speed of light, by each TDOA resulting in the range values it receives TOA measurements from UAS 111 and 112, and computes the three TDOA values as the positive arithmetic difference of pairs of TOA the LOCA method is shown below, from Schmidt [1], where the “Stations” represent known locations of sensors, and the “The Fix” is the location of the target.
While the invention has been particularly shown and described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the invention. Furthermore, it should be noted that the above-mentioned embodiments illustrate rather than limit the invention, and that those skilled in the art will be able to design many alternative embodiments without departing from the scope of the appended claims. In the claims, the word “comprising” does not exclude the presence of elements or steps other than those listed in a claim. The word “a” or “an” preceding an element does not exclude the presence of a plurality of such elements. In the device claim enumerating several means, several of these means can be embodied by one and the same item of hardware. The mere fact that certain measures are recited in mutually different dependent claims does not indicate that a combination of these measures cannot be used to advantage.
The present invention claims priority to U.S. Provisional Application No. 63/369,911 filed Jul. 29, 2022, the entirety of which is hereby incorporated by reference.
Number | Date | Country | |
---|---|---|---|
63369911 | Jul 2022 | US |