This disclosure generally pertains to a method of determining the location of an emitter using a single receiver unit.
The following references are to various papers describing localization methods and related methodologies. Throughout this disclosure, these references are cited using an abbreviated form the citation including square brackets and the bracketed one-digit or two-digit number preceding the full citation in the list below. E.g., in this disclosure ‘[4]’ is a citation to E. Xu, Z. Ding, and S. Dasgupta, “Reduced Complexity Semidefinite Relaxation Algorithms for Source Localization Based on Time Difference of Arrival,” IEEE Transactions on Mobile Computing, vol. 10, no. 9, pp. 1276-1282, 2010.
Emitter localization has been an area of interest in wireless communications because it is important for civil and military purposes and for detection of foreign bodies that could be jammers or unknown radio frequency (RF) signal sources. Emitter localization has been studied over the years, and various algorithms and techniques have been proposed to that effect.
For example, time of arrival (TOA), time difference of arrival (TDOA), angle of arrival (AOA), and received signal strength (RSS) methods have been studied in the literature for emitter localization. TOA-based algorithms use the propagation time of the signal from an emitter to a receiver or sensor for localization estimation. The signal is processed at different receivers, and localization is performed using the TOA measurements from these receivers. Typically, multiple receivers (RXs) are required for the TOA-based techniques. The work in [1], [2] shows algorithms based on TOA. TDOA is the difference between TOA measurements from the same emitter and the same transmission time [3]—[9].
Localization of the emitters can also be performed using angle of arrival (AOA) estimation, where several antenna arrays are used at several receivers [10], [11]. There have also been algorithms that use combinations of these methods. The work in [12] and [13] uses joint TDOA and AOA for emitter localization, and the work in [14] presents an algorithm based on RSS.
A Doppler frequency at an RX is a function of a relative velocity vector and a relative position vector between an RX and an electromagnetic interference (EMI) source transmitter (TX) or emitter. Therefore, Doppler frequency has been measured in practice. Furthermore, multiple RXs have been employed because an EMI TX position is a function of multiple variables, e.g., EMI TX position vector (x0, y0, z0) and velocity (vx0, vy0, vz0).
Most existing algorithms employ multiple RXs that measure and exchange the necessary information extracted from the received EMI source signal. A few algorithms use just one RX. References [5], [7]—[9] use more than one receiver, while [1] and [14] use just one receiver.
Investigators in the literature have assumed that RX positions of significant distances and RX velocities of significantly different directions can create more effective Doppler frequencies. For example, airplanes are flying at a minimum 9,260 meters (5 nm) apart, drones at 4.827 km (3 statute miles) apart, and satellites at 700 km altitude and 122 km apart.
The other differential Doppler (DD) method and its variants have been presented in the literature for many years. For example, in 1984, Torrieri analyzed a DD method that employs multiple physically separated RXs that intercept an EMI emitter signal at different positions and measure their respective Doppler frequencies [29]. Then, each RX sends the measured Doppler frequency information to a signal processing center RX. At each interception interval along a trajectory, the difference in frequency shifts among the multiple RXs is computed in the first step. Then, the position of the emitter is estimated based on the results of the previous step. These are the two steps involved in the DD method.
In 2008, Amar and Weiss improved upon this DD method with the direct position determination (DPD) approach [30]. This approach uses a single step without Doppler frequency measurement, whereas the DD method involves two steps: Step 1—measurement of Doppler frequency at each RX; and Step 2—measurement of the Doppler frequency differences among the RXs. In their approach, a fading coefficient, e.g., Rayleigh fading, was also considered in the RX signal at every interception interval, and an additive white Gaussian noise (AWGN) was added at each sampling time. Then, the maximum likelihood (ML) function was used as the objective cost function for a given fading coefficient. The EMI source position vector that maximizes the ML cost function among possible grid points was selected in a single step from all available intercepted data. In other words, no other steps were involved. It was shown that the DPD approach has less error than the conventional DD method at lower signal-to-noise-ratio (SNR) values. Both DD and DPD assumed at least two physically separated RXs; therefore, the RXs were positioned many kilometers apart, e.g., 10 km in [30].
Also, the DD method in [29] measures the Doppler frequencies at the multiple dislocated RXs at the same interception interval. Hence, calculating the difference between the two measured frequencies at two different RXs cancels out the unstable TX frequency component. Therefore, the Doppler frequency measured using the DD method in [29] will be free of the EMI TX frequency instability, whereas, the DPD in [30] does not measure the Doppler frequency but rather takes samples of the Doppler-shifted received signal for localization of the EMI TX.
In a recent study, an RF geolocation problem was also investigated using a single RX [31]. Doppler frequency and Doppler frequency rate were measured at different sampling times, and constrained unscented Kalman filtering (cUKF) was applied by converting almost noiseless nonlinear equations into system state vector-variable equations [31]. For example, a high SNR was assumed, e.g., 20 dB to 40 dB at the single RX, or zero AWGN was considered in the received baseband signal (equation (10) of [31]). In addition, no fading signal or EMI TX oscillator frequency drift (i.e., instability) was assumed when the Doppler frequency and Doppler rate were measured. The last assumption in [31] may prohibit the cUKF application in practice. This is because most EMI clocks have frequency instability, and the EMI TX frequency estimation has nonzero errors. Since a single RX is employed in [31], the Doppler frequencies should be measured at two contiguous sample time points and then subtracted to estimate the Doppler frequency rate. However, since most EMI TX oscillators are unstable, two different random frequencies will be added to the carrier frequency at the two contiguous interception points and then added to the true Doppler-shifted EMI carrier frequency. Hence, the Doppler frequency rate measurement will be in error when two Doppler frequencies are contaminated by the EMI TX clock's instability. Hence, the method in [31], which requires both Doppler frequency and Doppler frequency rate measurements, will fail under the TX frequency instability conditions.
Furthermore, the cUKF method in [31] requires an Earth surface constraint for convergence and accuracy. In other words, the cUKF method requires that the location of the EMI TX is restricted to the Earth's surface, whereas the proposed method in this current paper does not require such an Earth surface constraint. To reduce search time when the EMI TX location is in three-dimensional (3D) space, the arrival angle of EMI TX signal is exploited for the proposed method. In other words, the search range is restricted to a certain area instead of the entire 3D space.
In summary, the conventional DD method, e.g., in [29], requires multiple RXs and a two-step signal processing to find the location of an EMI emitter source and shows worse performance than the existing DPD method [30]. The DPD method requires only one-step processing and shows better performance than DD [30]. However, DPD requires multiple RXs separated by significant distances. The development, launching, and maintenance of multiple RXs is expensive. The cUKF in [31] considers only the known-signal case.
In one aspect, this disclosure provides an emitter localization method that uses only a single RX instead of multiple RXs to localize an EMI TX. Therefore, the proposed method can save costs significantly by a factor of the number of RXs used in traditional methods [5], [7]-[9], [29], [30].
In another aspect, this disclosure provides a random memoryless search and random memory search, which can reduce the searching time in [30] significantly, e.g., 27 times, and enhances its resolution.
In another aspect, this disclosure provides for the creation of multiple Doppler frequency effects via frequency-mixer (FMx) antennas at a single RX, even when the EMI TX uses a single carrier frequency.
In another aspect, this disclosure provides for the creation of a chirp signal at the RX to enhance the localization resolution by using multiple FMx antennas at the RX.
In another aspect, this disclosure provides an emitter localization method that applies for both narrowband and wideband EMI signals. This is because the maximum propagation delay spread between multiple RX antennas in a single RX unit can be Δd/c=3.33×10−9 seconds for Δd=1 meter separation. Hence, the bandwidth of the intercepted signal, which is the inverse of the multipath delay spread, is equal to c/Δd=3×108=0.3 GHz, where c is the speed of light. In contrast, the DPD method in [30] is applicable for only narrow-band emitter signals because of the larger separation distance, e.g., c/(Δd=10 km)=30 kHz.
In another aspect, this disclosure provides an EMI TX localization method that considers fading signal and TX frequency instability, and does not assume that the RX receives a high-power line-of-sight (LOS) EMI signal with 20-40 dB SNR, contra [31].
Other aspects will be in part apparent and in part pointed out hereinafter.
Corresponding parts are given corresponding reference characters throughout the drawings.
This disclosure aims to challenge the assumption that RXs physically separated by a larger distance can create better Doppler frequency effects and hence higher resolution in localization estimation. This paper claims that a single RX unit with multiple physically separated RX antennas can create different Doppler frequency effects and can locate the EMI emitter effectively, even if the EMI TX uses a single carrier frequency. To achieve that effect, this paper proposes using frequency-mixing (FMx) antennas or frequency-mixing intelligent reflecting surface (FMx-IRS) antennas as they appear in [15]. These antennas can create multiple carrier frequency signals at a single receiver, even if a single carrier frequency signal is transmitted by the EMI TX. The IRS has been investigated intensively and recently as a future wireless communication technology to improve channel capacity efficiency as well as secrecy channel capacity performance in wireless communications. For example, IRS studies are featured in [15]—[28].
This disclosure assumes a stationary emitter with zero velocity in the x, y, and z directions. Only the RX is mobile. The proposed method assumes no Global Positioning System (GPS) signal. Furthermore, the proposed method assumes a Rayleigh fading channel and includes the frequency instability of the EMI TX oscillator with the following uniform random variable distribution: uniform[−100, 100][Hz]. Therefore, the proposed method is considered under more practical and severe RF environments than the cUKF in [31].
The proposed method assumes that a modulated signal of only a single carrier frequency is transmitted by the EMI TX. Then, (L−1) subcarrier frequency signals are created in addition to the received main carrier by the (L−1) FMx [15] antennas at a single RX unit (e.g., L=2) and are processed with the signal received directly from the TX to the RX. L is the total number of antennas at the RX unit, and all antennas are in proximity. The separation distances between the main RX antenna and FMx antennas are nonzero, (e.g., Δd is about 1 m or less than 10 m). Therefore, the RX antenna and (L−1) FMx antennas can be located within a single RX unit but not co-located. A massive number of antennas can be co-located for millimeter wavelength carrier frequencies in the fifth generation (5G) and future wireless communication systems. The velocity of the main RX antenna and the (L−1) FMx antennas are all equal. The objective is to localize an EMI source. The received complex equivalent lowpass signal at the lth FMx or main antenna for the kth interception interval at time t is written [30] as Equation (1) below, where Tis the observation time interval, blk, is an unknown complex scalar representing the path attenuation (e.g., Rayleigh fading) at the kth interception interval observed by the lth FMx antenna, sk(t), which may be known or unknown depending on the application (e.g., known for search and rescue (SAR) and unknown for an intentional EMI emitter), is the observed signal envelope during the kth interception interval, flk is the down converted frequency, and wlk(t) is a complex AWGN∈CN(0, σn{circumflex over ( )}2), l=1, . . . , L, k=1, . . . , K.
r
lk(t)=blksk(t)ej2π
Let vlk, p0, and plk denote, respectively, the relative velocity vector between the stationary EMI TX and the mobile RX, the stationary EMI TX position vector, and the mobile lth FMx or the mobile main RX antenna position vector within a single RX at interception interval k.
The key idea of the proposed system is that the arrival frequencies at the RX are multiple frequencies, e.g., fc and fc+fr from the direct link main antenna and the indirect link via the FMx antenna, respectively, for L=2 in
The output frequency observed by the lth receiver antenna before the down conversion during the kth interception interval is given by Equation (2):
f
lk=[fcl+νk][1+μlk(p0)] (2)
where:
f
cl
=f
c+(l−1)fr, l=1, . . . , L (3)
fc is the known transmitted signal nominal frequency, which is also the input frequency at the lth antenna, fcl is the output frequency at the lth antenna, fr=fcΔf is the FMx frequency separation, νk is the unknown transmitted frequency shift due to source instability during the kth interception interval, and μlk (p0) is the Doppler frequency normalized by the carrier frequency fa and written as Equation (4):
The down-converted frequency flk−fcl is approximated according to Equation (5) because νkμlk(p0)<<fcl. The emitter frequency instability νk is modeled by a uniform random variable as Equation (6) and collects N samples per interception interval T (seconds). Hence, the sampling interval is represented by Equation (7).
lk≃νk+fclμlk(p0) (5)
νk=Unif[−100,100]Hz, (6)
T
s
=T/(N−1),0≤t≤T. (7)
For RX signal processing, the same direct position determination approach in [30] can be used by changing fc in [30] to fcl in Equation (3) throughout. The analog RX signal in Equation (1) is sampled as rlk[n]=rlk(nTs), n={0, 1 . . . , N−1} and represented by Equation (8), where Alk is the N×N Doppler shift diagonal matrix written as Equation (9), Ck is an N×N emitter frequency instability diagonal matrix written as Equation (10), sk is the N×1 TX signal vector, and wlk˜CN (0, σn2IN) is the N×1 complex Gaussian noise vector. In Equation (9), the μlk(p0) will be replaced by μlk(p) for the grid search or random search, where p is a position vector and candidate for EMI TX localization, and p0 is the true location.
r
lk
=b
lk
A
lk
C
k
s
k
+w
lk (8)
A
lk
diag{1,ej2πf
C
k
diag{1,ej2πν
The proposed algorithm maximizes the conditional likelihood function in (4) of [29] or the correlation cost function between the received vector including noise and the received vector excluding noise, given the observation vector rlk. Equivalently, the proposed algorithm minimizes the exponent in the likelihood function, called the minimum distance cost function, between the received signal vector including noise and the received vector excluding noise, as done in (5) of [29], [30]:
A Minimum Distance Cost function is defined by Equation (11).
A Maximum Correlation Cost function is defined by Equation (12).
Superscript H denotes the Hermitian. Equation (12) can be simplified further. The proposed algorithm will employ two simplified cost functions below from [30] for the unknown signal case and known signal case.
The objective function in Equation (12) is rewritten for the unknown signal case, and the best estimate of the EMI emitter location is found, respectively, as shown in Equations (13) and (14), where p is a possible grid point in a given search area, Vk and Qk are defined by Equations (15) and (16), and λmax(Qk) is the maximum eigenvalue of Qk, which is ≥0 since Qk is positive semi-definite.
For the known signal case, the objective function in (12) is rewritten, and the best estimate of the EMI emitter location is found, respectively, as Equations (17) or (18) and Equation (19), where ck is defined according to Equation (20) and Equation (21), zk is defined according to Equation (22), Ck is defined according to Equation (23), and Sk is defined according to Equation (24).
Sk and ck are transmitted together by the EMI. So, when Skck in (14) is treated as known, the EMI TX unstable frequency νk and the transmitted sequence vector sk (e.g., quadrature phase shift keying (QPSK)) are known as well.
However, the known-signal case means that only the modulated symbol sequence vector sk in the received signal vector rk is known because ck is unintentionally transmitted by the EMI. The frequency instability νk in the EMI TX signal is unknown at the RX in general, and νk should be unknown even for the known signal case. Therefore, Lks (p) in Equation (18) should be maximized with respect to the possible grid point vector p as well as the frequency instability parameter νk as in Equations (25)-(27), wherein Bk is defined by Equation (28), Gk is defined by Equation (29), and αm is the sum of elements on the m-th diagonal of matrix Gk. Hence, the cost function can be rewritten as Equation (30), which finds the νk, for each k, that maximizes the expression for βm in Equation (31). Therefore, the estimated EMI position for the known signal case, excluding the frequency instability νk, is given by Equation (32).
The same notations as the ones in [30] for Equations (1) to (33) are used in this paper because they are directly applicable for the 2D system model, 2D grid search, and even for the proposed method. The next section of this disclosure describes the proposed method for a 3D grid search. Since this paper considers localization of a stationary EMI emitter, it focuses on the unknown signal case rather than the known case.
For the existing 2D grid search in [30], the z-axis altitude is set to zero, as shown in
Since the altitude, 200 km, is much higher than one side of the search area, 10×10 (km×km), the cost functions will not be sensitive to the x-axis but only to the y-axis when the RX moves along a parallel line to the y-axis, and vice versa.
If a meter-level RMSE is desirable, then the grid size should be small, e.g., in meters because the grid search time increases exponentially. For example, if the grid size in the x-axis and y-axis, respectively, are Δdx=Δdy=0.06 km, then the number of grid points in the x-axis and y-axis would be Nx=Ny=10 km/Δdx=167, and the total number of grid points would be NxNy=167{circumflex over ( )}2=27,889. At each grid point, the cost function in Equation (13) is calculated. The correlation is computed between the actual intercepted signal rlk and a hypothetical received signal (AlkCksk) in Equation (12), which becomes Equation (13), assuming that the EMI TX is located at a hypothetical grid point, and the grid point p corresponding to the maximum correlation cost function in Equation (14) is selected. It takes 85.54 seconds for ten trials using a powerful workstation, e.g., Precision 5820 Tower. However, if Δdx=Δdy=0.001 km, then it will take more than a month for the simulation. One trial means K=10 interceptions in
To address these limitations, this disclosure contemplates two computationally efficient methods without sacrificing the RMSE performance: (a) a random memoryless search and (b) a random memory search method to reduce the grid search computation time significantly.
The QPSK symbol transmission rate is assumed to be 10 ksps. The RX does not need to have the transmitted signal information for the unknown signal case, and it takes N=100 samples of the received signal every interception interval for K=10 interception intervals and L=2 RX antennas. In the existing 2D method, the grid point vector that maximizes the cost function is selected out of Ngrid=NxNy grid points per trial. Then, the average RMSE is computed over Nexp=100 trials as Equation (33).
For a random memoryless search, a hypothetical position vector p is randomly generated in the specified search zone, e.g., 10×10 (km×km) in
For a random memory search, the final hypothetical position vector and its corresponding cost function from the previous trial are used for the current trial (unlike a random memoryless search where a hypothetical position vector is randomly generated for every trial). Once the initial starting position vector is established, the same procedure as that of the random memoryless search is conducted for the rest of the trial. Hence, there is a memory between adjacent trials.
The complexity in a grid search will be proportional to the total number of grid points, which is the square of the number of grid points per axis in a 2D search. The complexity in the proposed random memoryless or memory search is proportional to the number of randomly generated points, which can be significantly smaller than the total number of grid points. Hence, the random search can be much faster than the grid search at the same resolution. In addition, the resolution can be improved by the random search because the step size in a grid search determines the resolution; however, there is no step size in the random search. A search point in an axis can be any real number represented by a high number of bits.
An example Random Memory Search Algorithm can be conducted in the following steps:
STEP 0: Define all parameters, e.g., total number of trials=itotal=10.
STEP 1: Trial i=1.
STEP 2: For each SNR value and each trial value, set the cost function Lus(p)init to zero. For each SNR value and each trial value, generate the signal sk(t), noise wk(t), unknown complex scalar path attenuation blk, and channel phase ϕlk, using the given corresponding probability density functions.
STEP 3: Generate rlk using Equation (8).
STEP 4: Define trial initial position (px,init, py,init) as the center position, e.g., (px,init=5, py,init5).
STEP 5: Define a sufficient number of iterations (e.g., iiteration=10,000) for the ith trial search run.
STEP 6: Generate a new position (px,new, py,new) with Gaussian distributions: px,new˜N(px,init=5,√{square root over (px,init=5)}) and py,new˜N(py,init=5,√{square root over (py,init=5)}).
STEP 7: Compute μ1(p0) and μ2(p0) using Equation (4).
STEP 8: Compute A1k and A2k using Equation (9).
STEP 9: Compute V1k and V2k using Equation (15).
STEP 10: Compute Qk using Equation (16), and find the max eigenvalue of the Qk matrix.
STEP 11: Compute a new cost function Lus(p)new.
STEP 12: If Lus(p)new>Lus(p)init, then update (px,init, py,init) as (px,new, py,new). Also update the cost function with the latest value. Else, do not update position vector and cost function value. Go to STEP 6 until the iteration is complete for each trial.
STEP 13: Find the position estimate {circumflex over ( )}p(i) that gives the final value of Lus(p)new at the ith trial.
STEP 14: Increase i to i=i+1. Go to STEP 2 with Lus(p)init=Lus(p)new and (px,init, py,init)=(px,new, py,new) if i≤itotal. Else stop.
STEP 15: Compute the distance between the last trial position estimate {circumflex over ( )}p(itotal) and the true position p0.
Table 1 lists an example of computational complexity and running time using the Precision 5820 Tower workstation and MATLAB for the existing grid search, random memoryless search, and random memory search, when the scenarios in
Table 2 lists examples of the RMSE values for the existing grid search, random memoryless search, and random memory search when the same scenarios in
Observe from Table 2 that the random memory search, the random memoryless search, and the grid search with grid size 0.06 km show a similar RMSE resolution. The quality of the random memory search and the random memoryless search over the grid search is not degraded, although the search time can be significantly reduced, as indicated in Table 1.
To improve the performance further, the inventors contemplate using chirp RX processing with a finite number of FMx antennas. This approach requires a time-dependent RX frequency at each sample time. For the case without an RX chirp, L=2 number of RX frequencies are created using L=2 number of antennas (i.e., main RX antenna and one FMx antenna). For the chirp case, multiple sample time-dependent RX frequencies are created, e.g., a different RX frequency at each sample time, by switching from an FMx antenna to another FMx antenna in a set of (N−1) FMx antennas. If N samples are taken per interception interval, then it is desirable to employ N−1 number of FMx antennas and one main antenna at the RX side for each group l=1, . . . , L. Here, l represents the RX group index for the chirp case, whereas l represents the RX frequency index for the no-chirp case. Each group consists of n=1, N number of different sample time-dependent frequencies. Here, at the first sample time n=1, the RX frequency is generated via the main RX antenna, and at sample time n, 2≤n≤N, the RX frequency is generated via the nth FMx antenna, Hence, L(N−1) number of FMx antennas are necessary. If a millimeter wave is used by the EMI TX, then the wavelength is a few millimeters, and hence, (L=2)(N−1)=254 FMx antennas can be placed at one RX unit in an area less than 10×10 (cm×cm), assuming a wavelength separation between adjacent antenna elements.
The nth element of a diagonal Doppler shift matrix in Equation (9) is rewritten for the chip case as Equation (34), where l=1, L, n=1, N, and μlk(p0) denotes the normalized Doppler frequency fraction in Equation (4). Hence, different Doppler frequency effects can be created at each sampling time and can be applicable for the proposed method. Note again that the proposed method does not measure Doppler frequency as was done in [31] but rather uses only the received samples.
Accordingly, it can be seen that the present disclosure provides a system for determining a location of an emitter. A system according to the present disclosure generally comprises a single mobile frame configured for movement relative to the emitter (which can emit a known emitter signal (e.g., SAR application) or an unknown emitter signal (e.g., military reconnaissance application). For example, the frame can be an air frame of a SAR aircraft or reconnaissance aircraft, the body of a satellite, a ground vehicle chassis, or a water vessel, etc.
As shown in
As explained above, this system enables localization of the emitter with only one mobile frame because the main receiver receives two Doppler frequencies at offset locations. The first Doppler frequency is derivable from the direct emitter signal, and the second Doppler frequency is derivable from the frequency-mixed signal. Accordingly, the localization methods of the prior art that rely on two Doppler frequencies to establish emitter location can be adapted for use with the above described single-frame localization system.
In an exemplary embodiment, the system further comprises a processor connected to the main receiver and configured to determine the first Doppler frequency and the second Doppler frequency from the direct emitter signal and the frequency-mixed signal, respectively. The processor can be adapted use any suitable multi-Doppler frequency localization method to determine the location of the emitter in a two-dimensional grid or a three-dimensional space. In an exemplary embodiment, the processor is configured to determine the location of the emitter based on the first Doppler frequency and the second Doppler frequency using a direct position determination approach as described above and adapted from [30]. In certain embodiments, the processor is configured to use one of a random memoryless search algorithm or a random memory search algorithm as described above to reduce the grid search computation time for a defined search area. In one or more embodiments, the processor is configured to execute chirp RX processing techniques as described above to further improve localization performance.
The processor mentioned above may reside in a computer including a variety of computer hardware, including memory for storing processor-executable instructions. In general, the processor is configured to execute the programs stored in memory for carrying out the above functions. Components of the programs may reside at various times in different storage components of a computing device, and are executed by the processor of the device.
This disclosure is not limited to any particular computing environment. Examples of computing systems, environments, and/or configurations that may be suitable for use with aspects of the invention include, but are not limited to, personal computers, server computers, hand-held or laptop devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, mobile telephones, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.
Embodiments of the aspects of the present disclosure may be described in the general context of data and/or processor-executable instructions, such as program modules, stored one or more tangible, non-transitory storage media and executed by one or more processors or other devices. Generally, program modules include, but are not limited to, routines, programs, objects, components, and data structures that perform particular tasks or implement particular abstract data types. Aspects of the present disclosure may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote storage media including memory storage devices.
In operation, processors, computers and/or servers may execute the processor-executable instructions (e.g., software, firmware, and/or hardware) such as those illustrated herein to implement aspects of the invention.
Embodiments may be implemented with processor-executable instructions. The processor-executable instructions may be organized into one or more processor-executable components or modules on a tangible processor readable storage medium. Also, embodiments may be implemented with any number and organization of such components or modules. For example, aspects of the present disclosure are not limited to the specific processor-executable instructions or the specific components or modules illustrated in the figures and described herein. Other embodiments may include different processor-executable instructions or components having more or less functionality than illustrated and described herein.
The order of execution or performance of the operations in accordance with aspects of the present disclosure illustrated and described herein is not essential, unless otherwise specified. That is, the operations may be performed in any order, unless otherwise specified, and embodiments may include additional or fewer operations than those disclosed herein. For example, it is contemplated that executing or performing a particular operation before, contemporaneously with, or after another operation is within the scope of the invention.
When introducing elements of the present disclosure or the preferred embodiment(s) thereof, the articles “a”, “an”, “the” and “said” are intended to mean that there are one or more of the elements. The terms “comprising”, “including” and “having” are intended to be inclusive and mean that there may be additional elements other than the listed elements.
In view of the above, it will be seen that the several objects of the disclosure are achieved and other advantageous results attained.
As various changes could be made in the above products and methods without departing from the scope of the disclosure, it is intended that all matter contained in the above description shall be interpreted as illustrative and not in a limiting sense.
This application claims priority to U.S. Provisional Patent Application No. 63/234,439, filed Aug. 18, 2021, and to 63/371,645, filed Aug. 17, 2022 each of which is hereby incorporated by reference in its entirety.
This invention was made with government support under Grant Nos. FA9453-21-2-0011, FA9453-19-2-0003, and FA9453-21-0565 awarded by the Air Force Research Laboratory. The government of the United States has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
63371645 | Aug 2022 | US | |
63234439 | Aug 2021 | US |