The present invention relates to spectrum analysis and management for radio frequency signals, and more particularly for automatically detecting devices operating in white space.
Generally, it is known in the prior art to provide wireless communications spectrum management for detecting devices for managing the space. Spectrum management includes the process of regulating the use of radio frequencies to promote efficient use and gain net social benefit. A problem faced in effective spectrum management is the various numbers of devices emanating wireless signal propagations at different frequencies and across different technological standards. Coupled with the different regulations relating to spectrum usage around the globe effective spectrum management becomes difficult to obtain and at best can only be reached over a long period of time.
Another problem facing effective spectrum management is the growing need from spectrum despite the finite amount of spectrum available. Wireless technologies have exponentially grown in recent years. Consequently, available spectrum has become a valuable resource that must be efficiently utilized. Therefore, systems and methods are needed to effectively manage and optimize the available spectrum that is being used.
Most spectrum management devices may be categorized into two primary types. The first type is a spectral analyzer where a device is specifically fitted to run a ‘scanner’ type receiver that is tailored to provide spectral information for a narrow window of frequencies related to a specific and limited type of communications standard, such as cellular communication standard. Problems arise with these narrowly tailored devices as cellular standards change and/or spectrum use changes impact the spectrum space of these technologies. Changes to the software and hardware for these narrowly tailored devices become too complicated, thus necessitating the need to purchase a totally different and new device. Unfortunately, this type of device is only for a specific use and cannot be used to alleviate the entire needs of the spectrum management community.
The second type of spectral management device employs a methodology that requires bulky, extremely difficult to use processes, and expensive equipment. In order to attain a broad spectrum management view and complete all the necessary tasks, the device ends up becoming a conglomerate of software and hardware devices that is both hard to use and difficult to maneuver from one location to another.
While there may be several additional problems associated with current spectrum management devices, at least four major problems exist overall: 1) most devices are built to inherently only handle specific spectrum technologies such as 900 MHZ cellular spectrum while not being able to mitigate other technologies that may be interfering or competing with that spectrum, 2) the other spectrum management devices consist of large spectrum analyzers, database systems, and spectrum management software that is expensive, too bulky, and too difficult to manage for a user's basic needs, 3) other spectrum management devices in the prior art require external connectivity to remote databases to perform analysis and provide results or reports with analytics to aid in management of spectrum and/or devices, and 4) other devices of the prior art do not function to provide real-time or near real-time data and analysis to allow for efficient management of the space and/or devices and signals therein.
Examples of relevant prior art documents include the following:
However, none of the prior art references provide solutions to the limitations and longstanding unmet needs existing in this area for detecting signal emitting devices. Thus, there remains a need for automated device sensing in white space for wireless communications.
The present invention addresses the longstanding, unmet needs existing in the prior art and commercial sectors to provide solutions to the at least four major problems existing before the present invention. The present invention relates to systems, methods, and devices of the various embodiments enable spectrum management by identifying, classifying, and cataloging signals of interest based on radio frequency measurements. In an embodiment, signals and the parameters of the signals may be identified and indications of available frequencies may be presented to a user. In another embodiment, the protocols of signals may also be identified. In a further embodiment, the modulation of signals, data types carried by the signals, and estimated signal origins may be identified.
It is an object of this invention is to provide an apparatus for identifying signal emitting devices including: a housing, at least one processor and memory, and sensors constructed and configured for sensing and measuring wireless communications signals from signal emitting devices in a spectrum associated with wireless communications: and wherein the apparatus is operable to automatically analyze the measured data to identify at least one signal emitting device in near real time from attempted detection and identification of the at least one signal emitting device.
The present invention further provides systems for identifying signal emitting devices including at least one apparatus, wherein the at least one apparatus is operable for network-based communication with at least one server computer including a database, and/or with at least one other apparatus, but does not require a connection to the at least one server computer to be operable for identifying signal emitting devices: wherein each of the apparatus is operable for identifying signal emitting devices including: a housing, at least one processor and memory, and sensors constructed and configured for sensing and measuring wireless communications signals from signal emitting devices in a spectrum associated with wireless communications: and wherein the apparatus is operable to automatically analyze the measured data to identify at least one signal emitting device in near real time from attempted detection and identification of the at least one signal emitting device.
The present invention is further directed to a method for identifying signal emitting devices including the steps of: providing a device for measuring characteristics of signals from signal emitting devices in a spectrum associated with wireless communications, with measured data characteristics including frequency, power, bandwidth, duration, modulation, and combinations thereof: the device including a housing, at least one processor and memory, and sensors constructed and configured for sensing and measuring wireless communications signals within the spectrum: and further including the following steps performed within the device housing: assessing whether the measured data includes analog and/or digital signal(s): determining a best fit based on frequency, if the measured power spectrum is designated in an historical or a reference database(s) for frequency ranges; automatically determining a category for either analog or digital signals, based on power and sideband combined with frequency allocation: determining a TDM/FDM/CDM signal, based on duration and bandwidth: and identifying at least one signal emitting device from the composite results of the foregoing steps.
These and other aspects of the present invention will become apparent to those skilled in the art after a reading of the following description of the preferred embodiment when considered with the drawings, as they support the claimed invention.
The accompanying drawings, which are incorporated herein and constitute part of this specification, illustrate exemplary embodiments of the invention, and together with the general description given above and the detailed description given below, serve to explain the features of the invention.
Referring now to the drawings in general, the illustrations are for the purpose of describing at least one preferred embodiment and/or examples of the invention and are not intended to limit the invention thereto. Various embodiments are described in detail with reference to the accompanying drawings. Wherever possible, the same reference numbers are used throughout the drawings to refer to the same or like parts. References made to particular examples and implementations are for illustrative purposes, and are not intended to limit the scope of the invention or the claims.
The word “exemplary.” is used herein to mean “serving as an example, instance, or illustration.” Any implementation described herein as “exemplary” is not necessarily to be construed as preferred or advantageous over other implementations.
The present invention provides systems, methods, and devices for spectrum analysis and management by identifying, classifying, and cataloging at least one or a multiplicity of signals of interest based on radio frequency measurements and location and other measurements, and using near real-time parallel processing of signals and their corresponding parameters and characteristics in the context of historical and static data for a given spectrum.
The systems, methods and apparatus according to the present invention preferably have the ability to detect in near real time, and more preferably to detect, sense, measure, and/or analyze in near real time, and more preferably to perform any near real time operations within about 1 second or less. Advantageously, the present invention and its real time functionality described herein uniquely provide and enable the apparatus units to compare to historical data, to update data and/or information, and/or to provide more data and/or information on the open space, on the device that may be occupying the open space, and combinations, in the near real time compared with the historically scanned (15 min to 30 days) data, or historical database information.
The systems, methods, and devices of the various embodiments enable spectrum management by identifying, classifying, and cataloging signals of interest based on radio frequency measurements. In an embodiment, signals and the parameters of the signals may be identified and indications of available frequencies may be presented to a user. In another embodiment, the protocols of signals may also be identified. In a further embodiment, the modulation of signals, data types carried by the signals, and estimated signal origins may be identified.
Embodiments are directed to a spectrum management device that may be configurable to obtain spectrum data over a wide range of wireless communication protocols. Embodiments may also provide for the ability to acquire data from and sending data to database depositories that may be used by a plurality of spectrum management customers.
In one embodiment, a spectrum management device may include a signal spectrum analyzer that may be coupled with a database system and spectrum management interface. The device may be portable or may be a stationary installation and may be updated with data to allow the device to manage different spectrum information based on frequency, bandwidth, signal power, time, and location of signal propagation, as well as modulation type and format and to provide signal identification, classification, and geo-location. A processor may enable the device to process spectrum power density data as received and to process raw I/Q complex data that may be used for further signal processing, signal identification, and data extraction.
In an embodiment, a spectrum management device may comprise a low noise amplifier that receives a radio frequency (RF) energy from an antenna. The antenna may be any antenna structure that is capable of receiving RF energy in a spectrum of interest. The low noise amplifier may filter and amplify the RF energy. The RF energy may be provided to an RF translator. The RF translator may perform a fast Fourier transform (FFT) and either a square magnitude or a fast convolution spectral periodogram function to convert the RF measurements into a spectral representation. In an embodiment, the RF translator may also store a timestamp to facilitate calculation of a time of arrival and an angle of arrival. The In-Phase and Quadrature (I/Q) data may be provided to a spectral analysis receiver or it may be provided to a sample data store where it may be stored without being processed by a spectral analysis receiver. The input RF energy may also be directly digital down-converted and sampled by an analog to digital converter (ADC) to generate complex I/Q data. The complex I/Q data may be equalized to remove multipath, fading, white noise and interference from other signaling systems by fast parallel adaptive filter processes. This data may then be used to calculate modulation type and baud rate. Complex sampled I/Q data may also be used to measure the signal angle of arrival and time of arrival. Such information as angle of arrival and time of arrival may be used to compute more complex and precise direction finding. In addition, they may be used to apply geo-location techniques. Data may be collected from known signals or unknown signals and time spaced in order to provide expedient information. I/Q sampled data may contain raw signal data that may be used to demodulate and translate signals by streaming them to a signal analyzer or to a real-time demodulator software defined radio that may have the newly identified signal parameters for the signal of interest. The inherent nature of the input RF allows for any type of signal to be analyzed and demodulated based on the reconfiguration of the software defined radio interfaces.
A spectral analysis receiver may be configured to read raw In-Phase (I) and Quadrature (Q) data and either translate directly to spectral data or down convert to an intermediate frequency (IF) up to half the Nyquist sampling rate to analyze the incoming bandwidth of a signal. The translated spectral data may include measured values of signal energy, frequency, and time. The measured values provide attributes of the signal under review that may confirm the detection of a particular signal of interest within a spectrum of interest. In an embodiment, a spectral analysis receiver may have a referenced spectrum input of 0 Hz to 12.4 GHz with capability of fiber optic input for spectrum input up to 60 GHZ.
In an embodiment, the spectral analysis receiver may be configured to sample the input RF data by fast analog down-conversion of the RF signal. The down-converted signal may then be digitally converted and processed by fast convolution filters to obtain a power spectrum. This process may also provide spectrum measurements including the signal power, the bandwidth, the center frequency of the signal as well as a Time of Arrival (TOA) measurement. The TOA measurement may be used to create a timestamp of the detected signal and/or to generate a time difference of arrival iterative process for direction finding and fast triangulation of signals. In an embodiment, the sample data may be provided to a spectrum analysis module. In an embodiment, the spectrum analysis module may evaluate the sample data to obtain the spectral components of the signal.
In an embodiment, the spectral components of the signal may be obtained by the spectrum analysis module from the raw I/Q data as provided by an RF translator. The I/Q data analysis performed by the spectrum analysis module may operate to extract more detailed information about the signal, including by way of example, modulation type (e.g., FM, AM, QPSK, 16QAM, etc.) and/or protocol (e.g., GSM, CDMA, OFDM, LTE, etc.). In an embodiment, the spectrum analysis module may be configured by a user to obtain specific information about a signal of interest. In an alternate embodiment, the spectral components of the signal may be obtained from power spectral component data produced by the spectral analysis receiver.
In an embodiment, the spectrum analysis module may provide the spectral components of the signal to a data extraction module. The data extraction module may provide the classification and categorization of signals detected in the RF spectrum. The data extraction module may also acquire additional information regarding the signal from the spectral components of the signal. For example, the data extraction module may provide modulation type, bandwidth, and possible system in use information. In another embodiment, the data extraction module may select and organize the extracted spectral components in a format selected by a user.
The information from the data extraction module may be provided to a spectrum management module. The spectrum management module may generate a query to a static database to classify a signal based on its components. For example, the information stored in static database may be used to determine the spectral density, center frequency, bandwidth, baud rate, modulation type, protocol (e.g., GSM, CDMA, OFDM, LTE, etc.), system or carrier using licensed spectrum, location of the signal source, and a timestamp of the signal of interest. These data points may be provided to a data store for export. In an embodiment and as more fully described below, the data store may be configured to access mapping software to provide the user with information on the location of the transmission source of the signal of interest. In an embodiment, the static database includes frequency information gathered from various sources including, but not limited to, the Federal Communication Commission, the International Telecommunication Union, and data from users. As an example, the static database may be an SQL database. The data store may be updated, downloaded or merged with other devices or with its main relational database. Software API applications may be included to allow database merging with third-party spectrum databases that may only be accessed securely.
In the various embodiments, the spectrum management device may be configured in different ways. In an embodiment, the front end of system may comprise various hardware receivers that may provide In-Phase and Quadrature complex data. The front end receiver may include API set commands via which the system software may be configured to interface (i.e., communicate) with a third party receiver. In an embodiment, the front end receiver may perform the spectral computations using FFT (Fast Fourier Transform) and other DSP (Digital Signal Processing) to generate a fast convolution periodogram that may be re-sampled and averaged to quickly compute the spectral density of the RF environment.
In an embodiment, cyclic processes may be used to average and correlate signal information by extracting the changes inside the signal to better identify the signal of interest that is present in the RF space. A combination of amplitude and frequency changes may be measured and averaged over the bandwidth time to compute the modulation type and other internal changes, such as changes in frequency offsets, orthogonal frequency division modulation, changes in time (e.g., Time Division Multiplexing), and/or changes in I/Q phase rotation used to compute the baud rate and the modulation type. In an embodiment, the spectrum management device may have the ability to compute several processes in parallel by use of a multi-core processor and along with several embedded field programmable gate arrays (FPGA). Such multi-core processing may allow the system to quickly analyze several signal parameters in the RF environment at one time in order to reduce the amount of time it takes to process the signals. The amount of signals computed at once may be determined by their bandwidth requirements. Thus, the capability of the system may be based on a maximum frequency Fs/2. The number of signals to be processed may be allocated based on their respective bandwidths. In another embodiment, the signal spectrum may be measured to determine its power density, center frequency, bandwidth and location from which the signal is emanating and a best match may be determined based on the signal parameters based on information criteria of the frequency.
In another embodiment, a GPS and direction finding location (DF) system may be incorporated into the spectrum management device and/or available to the spectrum management device. Adding GPS and DF ability may enable the user to provide a location vector using the National Marine Electronics Association's (NMEA) standard form. In an embodiment, location functionality is incorporated into a specific type of GPS unit, such as a U.S. government issued receiver. The information may be derived from the location presented by the database internal to the device, a database imported into the device, or by the user inputting geo-location parameters of longitude and latitude which may be derived as degrees, minutes and seconds, decimal minutes, or decimal form and translated to the necessary format with the default being ‘decimal’ form. This functionality may be incorporated into a GPS unit. The signal information and the signal classification may then be used to locate the signaling device as well as to provide a direction finding capability.
A type of triangulation using three units as a group antenna configuration performs direction finding by using multilateration. Commonly used in civil and military surveillance applications, multilateration is able to accurately locate an aircraft, vehicle, or stationary emitter by measuring the “Time Difference of Arrival” (TDOA) of a signal from the emitter at three or more receiver sites. If a pulse is emitted from a platform, it will arrive at slightly different times at two spatially separated receiver sites, the TDOA being due to the different distances of each receiver from the platform. This location information may then be supplied to a mapping process that utilizes a database of mapping images that are extracted from the database based on the latitude and longitude provided by the geo-location or direction finding device. The mapping images may be scanned in to show the points of interest where a signal is either expected to be emanating from based on the database information or from an average taken from the database information and the geo-location calculation performed prior to the mapping software being called. The user can control the map to maximize or minimize the mapping screen to get a better view which is more fit to provide information of the signal transmissions. In an embodiment, the mapping process does not rely on outside mapping software. The mapping capability has the ability to generate the map image and to populate a mapping database that may include information from third party maps to meet specific user requirements.
In an embodiment, triangulation and multilateration may utilize a Bayesian type filter that may predict possible movement and future location and operation of devices based on input collected from the TDOA and geolocation processes and the variables from the static database pertaining to the specified signal of interest. The Bayesian filter takes the input changes in time difference and its inverse function (i.e., frequency difference) and takes an average change in signal variation to detect and predict the movement of the signals. The signal changes are measured within 1 ns time difference and the filter may also adapt its gradient error calculation to remove unwanted signals that may cause errors due to signal multipath, inter-symbol interference, and other signal noise.
In an embodiment the changes within a 1 ns time difference for each sample for each unique signal may be recorded. The spectrum management device may then perform the inverse and compute and record the frequency difference and phase difference between each sample for each unique signal. The spectrum management device may take the same signal and calculates an error based on other input signals coming in within the 1 ns time and may average and filter out the computed error to equalize the signal. The spectrum management device may determine the time difference and frequency difference of arrival for that signal and compute the odds of where the signal is emanating from based on the frequency band parameters presented from the spectral analysis and processor computations, and determines the best position from which the signal is transmitted (i.e., origin of the signal).
The signal processor 214 may include a signal detection module 216, a comparison module 222, a timing module 224, and a location module 225. Additionally, the signal processor 214 may include an optional memory module 226 which may include one or more optional buffers 228 for storing data generated by the other modules of the signal processor 214.
In an embodiment, the signal detection module 216 may operate to identify signals based on the RF energy measurements received from the RF receiver 210. The signal detection module 216 may include a Fast Fourier Transform (FFT) module 217 which may convert the received RF energy measurements into spectral representation data. The signal detection module 216 may include an analysis module 221 which may analyze the spectral representation data to identify one or more signals above a power threshold. A power module 220 of the signal detection module 216 may control the power threshold at which signals may be identified. In an embodiment, the power threshold may be a default power setting or may be a user selectable power setting. A noise module 219 of the signal detection module 216 may control a signal threshold, such as a noise threshold, at or above which signals may be identified. The signal detection module 216 may include a parameter module 218 which may determine one or more signal parameters for any identified signals, such as center frequency, bandwidth, power, number of detected signals, frequency peak, peak power, average power, signal duration, etc. In an embodiment, the signal processor 214 may include a timing module 224 which may record time information and provide the time information to the signal detection module 216. Additionally, the signal processor 214 may include a location module 225 which may receive location inputs from the location receiver 212 and determine a location of the spectrum management device 202. The location of the spectrum management device 202 may be provided to the signal detection module 216.
In an embodiment, the signal processor 214 may be connected to one or more memory 230. The memory 230 may include multiple databases, such as a history or historical database 232 and characteristics listing 236, and one or more buffers 240 storing data generated by signal processor 214. While illustrated as connected to the signal processor 214 the memory 230 may also be on chip memory residing on the signal processor 214 itself. In an embodiment, the history or historical database 232 may include measured signal data 234 for signals that have been previously identified by the spectrum management device 202. The measured signal data 234 may include the raw RF energy measurements, time stamps, location information, one or more signal parameters for any identified signals, such as center frequency, bandwidth, power, number of detected signals, frequency peak, peak power, average power, signal duration, etc., and identifying information determined from the characteristics listing 236, In an embodiment, the history or historical database 232 may be updated as signals are identified by the spectrum management device 202, In an embodiment, the characteristic listing 236 may be a database of static signal data 238. The static signal data 238 may include data gathered from various sources including by way of example and not by way of limitation the Federal Communication Commission, the International Telecommunication Union, telecom providers, manufacture data, and data from spectrum management device users. Static signal data 238 may include known signal parameters of transmitting devices, such as center frequency, bandwidth, power, number of detected signals, frequency peak, peak power, average power, signal duration, geographic information for transmitting devices, and any other data that may be useful in identifying a signal. In an embodiment, the static signal data 238 and the characteristic listing 236 may correlate signal parameters and signal identifications. As an example, the static signal data 238 and characteristic listing 236 may list the parameters of the local fire and emergency communication channel correlated with a signal identification indicating that signal is the local fire and emergency communication channel.
In an embodiment, the signal processor 214 may include a comparison module 222 which may match data generated by the signal detection module 216 with data in the history or historical database 232 and/or characteristic listing 236. In an embodiment the comparison module 222 may receive signal parameters from the signal detection module 216, such as center frequency, bandwidth, power, number of detected signals, frequency peak, peak power, average power, signal duration, and/or receive parameter from the timing module 224 and/or location module 225. The parameter match module 223 may retrieve data from the history or historical database 232 and/or the characteristic listing 236 and compare the retrieved data to any received parameters to identify matches. Based on the matches the comparison module may identify the signal. In an embodiment, the signal processor 214 may be optionally connected to a display 242, an input device 244, and/or network transceiver 246. The display 242 may be controlled by the signal processor 214 to output spectral representations of received signals, signal characteristic information, and/or indications of signal identifications on the display 242. In an embodiment, the input device 244 may be any input device, such as a keyboard and/or knob, mouse, virtual keyboard or even voice recognition, enabling the user of the spectrum management device 202 to input information for use by the signal processor 214. In an embodiment, the network transceiver 246 may enable the spectrum management device 202 to exchange data with wired and/or wireless networks, such as to update the characteristic listing 236 and/or upload information from the history or historical database 232.
In block 312 the processor 214 may identify one or more signal above a threshold. In an embodiment, the processor 214 may analyze the spectral representation data to identify a signal above a power threshold. A power threshold may be an amplitude measure selected to distinguish RF energies associated with actual signals from noise. In an embodiment, the power threshold may be a default value. In another embodiment, the power threshold may be a user selectable value. In block 314 the processor 214 may determine signal parameters of any identified signal or signals of interest. As examples, the processor 214 may determine signal parameters such as center frequency, bandwidth, power, number of detected signals, frequency peak, peak power, average power, signal duration for the identified signals. In block 316 the processor 214 may store the signal parameters of each identified signal, a location indication, and time indication for each identified signal in a history database 232. In an embodiment, a history database 232 may be a database resident in a memory 230 of the spectrum management device 202 which may include data associated with signals actually identified by the spectrum management device.
In block 318 the processor 214 may compare the signal parameters of each identified signal to signal parameters in a signal characteristic listing. In an embodiment, the signal characteristic listing may be a static database 238 stored in the memory 230 of the spectrum management device 202 which may correlate signal parameters and signal identifications. In determination block 320 the processor 214 may determine whether the signal parameters of the identified signal or signals match signal parameters in the characteristic listing 236. In an embodiment, a match may be determined based on the signal parameters being within a specified tolerance of one another. As an example, a center frequency match may be determined when the center frequencies are within plus or minus 1 kHz of each other. In this manner, differences between real world measured conditions of an identified signal and ideal conditions listed in a characteristics listing may be accounted for in identifying matches. If the signal parameters do not match (i.e., determination block 320)=“No”), in block 326 the processor 214 may display an indication that the signal is unidentified on a display 242 of the spectrum management device 202. In this manner, the user of the spectrum management device may be notified that a signal is detected, but has not been positively identified. If the signal parameters do match (i.e., determination block 320)=“Yes”), in block 324 the processor 214 may display an indication of the signal identification on the display 242. In an embodiment, the signal identification displayed may be the signal identification correlated to the signal parameter in the signal characteristic listing which matched the signal parameter for the identified signal. Upon displaying the indications in blocks 324 or 326 the processor 214 may return to block 302 and cyclically measure and identify further signals of interest.
If the average of the sample block is equal to or greater than the noise floor average (i.e., determination block 510=“Yes”), the sample block may potentially include a signal of interest and in block 512 the processor 214 may reset a measurement counter (C) to 1. The measurement counter value indicating which sample within a sample block is under analysis. In determination block 516 the processor 214 may determine whether the RF measurement of the next frequency sample (C) is greater than the signal power threshold. In this manner, the value of the measurement counter (C) may be used to control which sample RF measurement in the sample block is compared to the signal power threshold. As an example, when the counter (C) equals 1, the first RF measurement may be checked against the signal power threshold and when the counter (C) equals 2 the second RF measurement in the sample block may be checked, etc. If the C RF measurement is less than or equal to the signal power threshold (i.e., determination block 516=“No”), in determination block 517 the processor 214 may determine whether the cross block flag is set. If the cross block flag is not set (i.e., determination block 517=“No”), in determination block 522 the processor 214 may determine whether the end of the sample block is reached. If the end of the sample block is reached (i.e., determination block 522=“Yes”), in block 506 the processor 214 may load the next available sample block and proceed in blocks 508, 510, 514, and 512 as discussed above. If the end of the sample block is not reached (i.e., determination block 522=“No”), in block 524 the processor 214 may increment the measurement counter (C) so that the next sample in the sample block is analyzed.
If the C RF measurement is greater than the signal power threshold (i.e., determination block 516=“Yes”), in block 518 the processor 214 may check the status of the cross block flag to determine whether the cross block flag is set. If the cross block flag is not set (i.e., determination block 518=“No”), in block 520 the processor 214 may set a sample start. As an example, the processor 214 may set a sample start by indicating a potential signal of interest may be discovered in a memory by assigning a memory location for RF measurements associated with the sample start. Referring to
In determination block 530 the processor 214 may determine whether the C RF measurement (e.g., the next RF measurement because the value of the RF measurement counter was incremented) is greater than the signal power threshold. If the C RF measurement is greater than the signal power threshold (i.e., determination block 530)=“Yes”), in determination block 532 the processor 214 may determine whether the end of the sample block is reached. If the end of the sample block is not reached (i.e., determination block 532=“No”), there may be further RF measurements available in the sample block and in block 526 the processor 214 may store the C RF measurement in the memory location for the sample. In block 528 the processor may increment the measurement counter (C) and in determination block 530 determine whether the C RF measurement is above the signal power threshold and in block 532 determine whether the end of the sample block is reached. In this manner, successive sample RF measurements may be checked against the signal power threshold and stored until the end of the sample block is reached and/or until a sample RF measurement falls below the signal power threshold. If the end of the sample block is reached (i.e., determination block 532=“Yes”), in block 534 the processor 214 may set the cross block flag. In an embodiment, the cross block flag may be a flag in a memory available to the processor 214 indicating the signal potential spans across two or more sample blocks. In a further embodiment, prior to setting the cross block flag in block 534, the slope of a line drawn between the last two RF measurement samples may be used to determine whether the next sample block likely contains further potential signal samples. A negative slope may indicate that the signal of interest is fading and may indicate the last sample was the final sample of the signal of interest. In another embodiment, the slope may not be computed and the next sample block may be analyzed regardless of the slope.
If the end of the sample block is reached (i.e., determination block 532=“Yes”) and in block 534 the cross block flag is set, referring to
If the end of the sample block is reached (i.e., determination block 532=“Yes”) and in block 534 the cross block flag is set, referring to
If the mean is greater than the signal threshold (i.e., determination block 544=“Yes”), in block 546 the processor 214 may identify the sample as a signal of interest. In an embodiment, the processor 214 may identify the sample as a signal of interest by assigning a signal identifier to the signal, such as a signal number or sample number. In block 548 the processor 214 may determine and store signal parameters for the signal. As an example, the processor 214 may determine and store a frequency peak of the identified signal, a peak power of the identified signal, an average power of the identified signal, a signal bandwidth of the identified signal, and/or a signal duration of the identified signal. In block 552 the processor 214 may clear the cross block flag (or verify that the cross block flag is unset). In block 556 the processor 214 may determine whether the end of the sample block is reached. If the end of the sample block is not reached (i.e., determination block 556=“No”) in block 558 the processor 214 may increment the measurement counter (C), and referring to
In block 614 the processor 214 may compare the signal parameters of the identified signal to signal parameters in a signal characteristic listing 236. In an embodiment, the characteristic listing 236 may be a static database separate from the history database 232, and the characteristic listing 236 may correlate signal parameters with signal identifications. In determination block 616 the processor 214 may determine whether the signal parameters of the identified signal match any signal parameters in the signal characteristic listing 236. In an embodiment, the match in determination 616 may be a match based on a tolerance between the signal parameters of the identified signal and the parameters in the characteristic listing 236. If there is a match (i.e., determination block 616=“Yes”), in block 618 the processor 214 may indicate a match in the history database 232 and in block 622 may display an indication of the signal identification on a display 242. As an example, the indication of the signal identification may be a display of the radio call sign of an identified FM radio station signal. If there is not a match (i.e., determination block 616=“No”), in block 620 the processor 214 may display an indication that the signal is an unidentified signal. In this manner, the user may be notified a signal is present in the environment, but that the signal does not match to a signal in the characteristic listing.
If there are matches (i.e., determination block 706=“Yes”), in optional block 708 the processor 214 may display a plot of one or more of the signals matching the current location. As an example, the processor 214 may compute the average frequency over frequency intervals across a given spectrum and may display a plot of the average frequency over each interval. In block 712 the processor 214 may determine one or more open frequencies at the current location. As an example, the processor 214 may determine one or more open frequencies by determining frequency ranges in which no signals fall or at which the average is below a threshold. In block 714 the processor 214 may display an indication of one or more open frequency on a display 242 of the spectrum management device 202.
The protocol module 806 may identify the communication protocol (e.g., LTE, CDMA, etc.) associated with a signal of interest. In an embodiment, the protocol module 806 may use data retrieved from the characteristic listing, such as protocol data 804 to help identify the communication protocol. The symbol detector module 816 may determine symbol timing information, such as a symbol rate for a signal of interest. The protocol module 806 and/or symbol module 816 may provide data to the comparison module 222. The comparison module 222 may include a protocol match module 814 which may attempt to match protocol information for a signal of interest to protocol data 804 in the characteristic listing to identify a signal of interest. Additionally, the protocol module 806 and/or symbol module 816 may store data in the memory module 226 and/or history database 232. In an embodiment, the protocol module 806 and/or symbol module 816 may use protocol data 804 and/or other data from the characteristic listing 236 to help identify protocols and/or symbol information in signals of interest.
The optimization module 818 may gather information from the characteristic listing, such as noise figure parameters, antenna hardware parameters, and environmental parameters correlated with an identified signal of interest to calculate a degradation value for the identified signal of interest. The optimization module 818 may further control the display 242 to output degradation data enabling a user of the spectrum management device 802 to optimize a signal of interest.
In block 1110 the processor 214 may compare the signal parameters and protocol data of the identified signal to signal parameters and protocol data in the signal characteristic listing 236. In determination block 1112 the processor 214 may determine whether the signal parameters and protocol data of the identified signal match any signal parameters and protocol data in the signal characteristic listing 236. If there is a match (i.e., determination block 1112=“Yes”), in block 1114 the processor 214 may indicate a match in the history database and in block 1118 may display an indication of the signal identification and protocol on a display. If there is not a match (i.e., determination block 1112=“No”), in block 1116 the processor 214 may display an indication that the signal is an unidentified signal. In this manner, the user may be notified a signal is present in the environment, but that the signal does not match to a signal in the characteristic listing.
In an embodiment, based on signal detection, a time tracking module, such as a TDOA/FDOA module 1204, may track the frequency repetition interval at which the signal is changing. The frequency repetition interval may also be tracked for a burst signal. In an embodiment, the spectrum management device may measure the signal environment and set anchors based on information stored in the historic or static database about known transmitter sources and locations. In an embodiment, the phase information about a signal be extracted using a spectral decomposition correlation equation to measure the angle of arrival (“AOA”) of the signal. In an embodiment, the processor of the spectrum management device may determine the received power as the Received Signal Strength (“RSS”) and based on the AOA and RSS may measure the frequency difference of arrival. In an embodiment, the frequency shift of the received signal may be measured and aggregated over time. In an embodiment, after an initial sample of a signal, known transmitted signals may be measured and compared to the RSS to determine frequency shift error. In an embodiment, the processor of the spectrum management device may compute a cross ambiguity function of aggregated changes in arrival time and frequency of arrival. In an additional embodiment, the processor of the spectrum management device may retrieve FFT data for a measured signal and aggregate the data to determine changes in time of arrival and frequency of arrival. In an embodiment, the signal components of change in frequency of arrival may be averaged through a Kalman filter with a weighted tap filter from 2 to 256 weights to remove measurement error such as noise, multipath interference, etc. In an embodiment, frequency difference of arrival techniques may be applied when either the emitter of the signal or the spectrum management device are moving or when then emitter of the signal and the spectrum management device are both stationary. When the emitter of the signal and the spectrum management device are both stationary the determination of the position of the emitter may be made when at least four known other known signal emitters positions are known and signal characteristics may be available. In an embodiment, a user may provide the four other known emitters and/or may use already in place known emitters, and may use the frequency, bandwidth, power, and distance values of the known emitters and their respective signals. In an embodiment, where the emitter of the signal or spectrum management device may be moving, frequency deference of arrival techniques may be performed using two known emitters.
The processor 214 of spectrum management devices 202, 802 and 1202 may be any programmable microprocessor, microcomputer or multiple processor chip or chips that can be configured by software instructions (applications) to perform a variety of functions, including the functions of the various embodiments described above. In some devices, multiple processors may be provided, such as one processor dedicated to wireless communication functions and one processor dedicated to running other applications. Typically, software applications may be stored in the internal memory 226 or 230 before they are accessed and loaded into the processor 214. The processor 214 may include internal memory sufficient to store the application software instructions. In many devices the internal memory may be a volatile or nonvolatile memory, such as flash memory, or a mixture of both. For the purposes of this description, a general reference to memory refers to memory accessible by the processor 214 including internal memory or removable memory plugged into the device and memory within the processor 214 itself.
The present invention provides for systems, methods, and apparatus solutions for device sensing in white space, which improves upon the prior art by identifying sources of signal emission by automatically detecting signals and creating unique signal profiles. Device sensing has an important function and applications in military and other intelligence sectors, where identifying the emitter device is crucial for monitoring and surveillance, including specific emitter identification (SEI).
At least two key functions are provided by the present invention: signal isolation and device sensing. Signal Isolation according to the present invention is a process whereby a signal is detected, isolated through filtering and amplification, amongst other methods, and key characteristics extracted. Device Sensing according to the present invention is a process whereby the detected signals are matched to a device through comparison to device signal profiles and may include applying a confidence level and/or rating to the signal-profile matching. Further, device sensing covers technologies that permit storage of profile comparisons such that future matching can be done with increased efficiency and/or accuracy. The present invention systems, methods, and apparatus are constructed and configured functionally to identify any signal emitting device, including by way of example and not limitation, a radio, a cell phone, etc.
Regarding signal isolation, the following functions are included in the present invention: amplifying, filtering, detecting signals through energy detection, waveform-based, spectral correlation-based, radio identification-based, or matched filter method, identifying interference, identifying environmental baseline(s), and/or identify signal characteristics.
Regarding device sensing, the following functions are included in the present invention: using signal profiling and/or comparison with known database(s) and previously recorded profile(s), identifying the expected device or emitter, stating the level of confidence for the identification, and/or storing profiling and sensing information for improved algorithms and matching. In preferred embodiments of the present invention, the identification of the at least one signal emitting device is accurate to a predetermined degree of confidence between about 80 and about 95 percent, and more preferably between about 80 and about 100 percent. The confidence level or degree of confidence is based upon the amount of matching measured data compared with historical data and/or reference data for predetermined frequency and other characteristics.
The present invention provides for wireless signal-emitting device sensing in the white space based upon a measured signal, and considers the basis of license(s) provided in at least one reference database, preferably the federal communication commission (FCC) and/or other defined database including license listings. The methods include the steps of providing a device for measuring characteristics of signals from signal emitting devices in a spectrum associated with wireless communications, the characteristics of the measured data from the signal emitting devices including frequency, power, bandwidth, duration, modulation, and combinations thereof: making an assessment or categorization on analog and/or digital signal(s): determining the best fit based on frequency if the measured power spectrum is designated in historical and/or reference data, including but not limited to the FCC or other database(s) for select frequency ranges: determining analog or digital, based on power and sideband combined with frequency allocation: determining a TDM/FDM/CDM signal, based on duration and bandwidth: determining best modulation fit for the desired signal, if the bandwidth and duration match the signal database(s): adding modulation identification to the database: listing possible modulations with best percentage fit, based on the power, bandwidth, frequency, duration, database allocation, and combinations thereof: and identifying at least one signal emitting device from the composite results of the foregoing steps.
According to methods of the present invention, the following steps are performed automatically by the apparatus unit(s): based on the measured signal(s), input the basis of the license provided in the FCC and/or user defined database: measure the frequency, power, bandwidth, and/or duration of the measured signal(s); determine the best method of modulation identification: perform an assessment on analog or digital signal: if power spectrum is stored in designated FCC, historical, and/or user database frequency ranges, determine best fit based on frequency: based on power and sideband combined with frequency allocation, determine if analog or digital signal(s): based on duration and bandwidth determine a TDM/FDM/CDM signal: if bandwidth and duration match signal database then determine best modulation fit for the desired signal: and add modulation identification to the database and list possible modulations with best percentage fit based on the power, bandwidth, frequency, and/or duration and database allocation.
In embodiments of the present invention, an apparatus is provided for automatically identifying devices in a spectrum, the apparatus including a housing, at least one processor and memory, and sensors constructed and configured for sensing and measuring wireless communications signals from signal emitting devices in a spectrum associated with wireless communications: and wherein the apparatus is operable to automatically analyze the measured data to identify at least one signal emitting device in near real time from attempted detection and identification of the at least one signal emitting device. The characteristics of signals and measured data from the signal emitting devices include frequency, power, bandwidth, duration, modulation, and combinations thereof.
The present invention systems including at least one apparatus, wherein the at least one apparatus is operable for network-based communication with at least one server computer including a database, and/or with at least one other apparatus, but does not require a connection to the at least one server computer to be operable for identifying signal emitting devices: wherein each of the apparatus is operable for identifying signal emitting devices including: a housing, at least one processor and memory, and sensors constructed and configured for sensing and measuring wireless communications signals from signal emitting devices in a spectrum associated with wireless communications: and wherein the apparatus is operable to automatically analyze the measured data to identify at least one signal emitting device in near real time from attempted detection and identification of the at least one signal emitting device.
The present invention provides for systems, methods, and apparatus solutions for automatically identifying open space, including open space in the white space of a wireless communication spectrum. Importantly, the present invention identifies the open space as the space that is unused and/or seldomly used (and identifies the owner of the licenses for the seldomly used space, if applicable), including unlicensed spectrum, white space, guard bands, and combinations thereof. Method steps of the present invention include: automatically obtaining a listing or report of all frequencies in the frequency range: plotting a line and/or graph chart showing power and bandwidth activity: setting frequencies based on a frequency step and/or resolution so that only user-defined frequencies are plotted: generating a .csv or .pdf file showing the average and/or aggregated values of power, bandwidth and frequency for each derived frequency step: and showing an activity report over time, over day vs. night, over frequency bands if more than one, in white space if requested, in ISM space if requested; and if frequency space seldomly used in area list frequencies and license holders. Additional steps include: scanning the frequency span, wherein a default scan includes a frequency span between about 54 MHZ and about 804 MHZ: an ISM scan between about 900 MHz and about 2.5 GHZ: an ISM scan between about 5 GHZ and about 5.8 GHZ: and/or a frequency range based upon inputs provided by a user. Also, method steps include scanning for an allotted amount of time between a minimum of about 15 minutes up to about 30 days: preferably scanning for allotted times selected from the following: a minimum of about 15 minutes; about 30 minutes; about 1 hour increments; about 5 hour increments; about 10 hour increments; about 24 hours; about 1 day: and about up to 30 days: and combinations thereof. In preferred embodiments, if the apparatus is configured for automatically scanning for more than about 15 minutes, then the apparatus is preferably set for updating results, including updating graphs and/or reports for an approximately equal amount of time (e.g., every 15 minutes). The systems, methods, and apparatus also provide for automatically calculating a percent activity on each frequency band.
Various reports for describing and illustrating with visualization the data and analysis of the device, system and method results from spectrum management activities include at least reports on power usage. RF survey, and variance.
The systems, methods, and devices of the various embodiments enable spectrum management by identifying, classifying, and cataloging signals of interest based on radio frequency measurements. In an embodiment, signals and the parameters of the signals may be identified and indications of available frequencies may be presented to a user. In another embodiment, the protocols of signals may also be identified. In a further embodiment, the modulation of signals, devices or device types emitting signals, data types carried by the signals, and estimated signal origins may be identified.
Referring again to the drawings,
By way of example, and not limitation, the computing devices 3820, 3830, 3840 are intended to represent various forms of digital devices and mobile devices, such as a server, blade server, mainframe, mobile phone, a personal digital assistant (PDA), a smart phone, a desktop computer, a netbook computer, a tablet computer, a workstation, a laptop, and other similar computing devices. The components shown here, their connections and relationships, and their functions, are meant to be exemplary only, and are not meant to limit implementations of the invention described and/or claimed in this document.
In one embodiment, the computing device 3820 includes components such as a processor 3860, a system memory 3862 having a random access memory (RAM) 3864 and a read-only memory (ROM) 3866, and a system bus 3868 that couples the memory 3862 to the processor 3860. In another embodiment, the computing device 3830 may additionally include components such as a storage device 3890 for storing the operating system 3892 and one or more application programs 3894, a network interface unit 3896, and/or an input/output controller 3898. Each of the components may be coupled to each other through at least one bus 3868. The input/output controller 3898 may receive and process input from, or provide output to, a number of other devices 3899, including, but not limited to, alphanumeric input devices, mice, electronic styluses, display units, touch screens, signal generation devices (e.g., speakers) or printers.
By way of example, and not limitation, the processor 3860 may be a general-purpose microprocessor (e.g., a central processing unit (CPU)), a graphics processing unit (GPU), a microcontroller, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA), a Programmable Logic Device (PLD), a controller, a state machine, gated or transistor logic, discrete hardware components, or any other suitable entity or combinations thereof that can perform calculations, process instructions for execution, and/or other manipulations of information.
In another implementation, shown in
Also, multiple computing devices may be connected, with each device providing portions of the necessary operations (e.g., a server bank, a group of blade servers, or a multi-processor system). Alternatively, some steps or methods may be performed by circuitry that is specific to a given function.
According to various embodiments, the computer system 3800 may operate in a networked environment using logical connections to local and/or remote computing devices 3820, 3830, 3840 through a network 3810. A computing device 3820 may connect to a network 3810 through a network interface unit 3896 connected to the bus 3868. Computing devices may communicate communication media through wired networks, direct-wired connections or wirelessly such as acoustic, RF or infrared through a wireless communication antenna 3897 in communication with the network's wireless communication antenna 3812 and the network interface unit 3896, which may include digital signal processing circuitry when necessary. The network interface unit 3896 may provide for communications under various modes or protocols.
In one or more exemplary aspects, the instructions may be implemented in hardware, software, firmware, or any combinations thereof. A computer readable medium may provide volatile or non-volatile storage for one or more sets of instructions, such as operating systems, data structures, program modules, applications or other data embodying any one or more of the methodologies or functions described herein. The computer readable medium may include the memory 3862, the processor 3860, and/or the storage device 3890 and may be a single medium or multiple media (e.g., a centralized or distributed computer system) that store the one or more sets of instructions 3900. Non-transitory computer readable media includes all computer readable media, with the sole exception being a transitory, propagating signal per se. The instructions 3900 may further be transmitted or received over the network 3810 via the network interface unit 3896 as communication media, which may include a modulated data signal such as a carrier wave or other transport mechanism and includes any delivery media. The term “modulated data signal” means a signal that has one or more of its characteristics changed or set in a manner as to encode information in the signal.
Storage devices 3890 and memory 3862 include, but are not limited to, volatile and non-volatile media such as cache, RAM, ROM, EPROM, EEPROM, FLASH memory, or other solid state memory technology: discs (e.g., digital versatile discs (DVD), HD-DVD, BLU-RAY, compact disc (CD), or CD-ROM) or other optical storage: magnetic cassettes, magnetic tape, magnetic disk storage, floppy disks, or other magnetic storage devices; or any other medium that can be used to store the computer readable instructions and which can be accessed by the computer system 3800.
It is also contemplated that the computer system 3800 may not include all of the components shown in
The present invention further provides for aggregating data from at least two apparatus units by at least one server computer and storing the aggregated data in a database and/or in at least one database in a cloud-based computing environment or virtualized computing environment, as illustrated in
In other embodiments of the present invention, which include the base invention described hereinabove, and further including the functions of machine “learning”, modulation detection, automatic signal detection, FFT replay, and combinations thereof.
Automatic modulation detection and machine “learning” includes automatic signal variance determination by at least one of the following methods: date and time from location set, and remote access to the apparatus unit to determine variance from different locations and times, in addition to the descriptions of automatic signal detection and threshold determination and setting. Environments vary, especially where there are many signals, noise, interference, variance, etc., so tracking signals automatically is difficult, and a longstanding, unmet need in the prior art. The present invention provides for automatic signal detection using a sample of measured and sensed data associated with signals over time using the at least one apparatus unit of the present invention to provide an automatically adjustable and adaptable system. For each spectrum scan, the data is automatically subdivided into “windows”, which are sections or groups of data within a frequency space. Real-time processing of the measured and sensed data on the apparatus unit(s) or devices combined with the windowing effect provides for automatic comparison of signal versus noise within the window to provide for noise approximation, wherein both signals and noise are measured and sensed, recorded, analyzed compared with historical data to identify and output signals in a high noise environment. It is adaptive and iterative to include focused windows and changes in the window or frequency ranges grouped. The resulting values for all data are squared in the analysis, which results in signals identified easily by the apparatus unit as having significantly larger power values compared with noise: additional analytics provide for selection of the highest power value signals and review of the original data corresponding thereto. Thus, the at least one apparatus automatically determines and identifies signals compared to noise in the RF spectrum.
The apparatus unit or device of the present invention further includes a temporal anomaly detector (or “learning channel”). The first screen shot illustrated in
In a similar capacity, the user can load a spreadsheet that they have constructed on their own to describe the channels that they expect to see in a given environment, as illustrated in
Automatic signal detection of the present invention eliminates the need for a manual setting of a power threshold line or bar, as with the prior art. The present invention does not require a manual setting of power threshold bar or flat line to identify signals instead of noise, instead it uses information on the hardware parameters of the apparatus unit or device, environment parameters, and terrain data to derive the threshold bar or flatline, which are stored in the static database of the apparatus unit or device. Thus, the apparatus unit or device may be activated and left unattended to collect data continuously without the need for manual interaction with the device directly. Furthermore, the present invention allows remote viewing of live data in real time on a display of a computer or communications device in network-based connection but remotely positioned from the apparatus unit or device, and/or remote access to device settings, controls, data, and combinations thereof. The network-based communication may be selected from mobile, satellite, Ethernet, and functional equivalents or improvements with security including firewalls, encryption of data, and combinations thereof.
Regarding FFT replay, the present invention apparatus units are operable to replay data and to review and/or replay data saved based upon an unknown event, such as for example and not limitation, reported alarms and/or unique events, wherein the FFT replay is operable to replay stored sensed and measured data to the section of data nearest the reported alarm and/or unique event. By contrast, prior art provides for recording signals on RF spectrum measurement devices, which transmit or send the raw data to an external computer for analysis, so then it is impossible to replay or review specific sections of data, as they are not searchable, tagged, or otherwise sectioned into subgroups of data or stored on the device.
The previous approach to ASD was to subtract a calibration vector from each FFT sample set (de-bias), then square each resulting value and look for concentrations of energy that would differentiate a signal from random baseline noise. The advantages of this approach are that, by the use of the calibration vector (which was created using the receiver itself with no antenna), we are able to closely track variations in the baseline noise that are due to the characteristics of the receiver, front end filtering, attenuation and A/D converter hardware. On most modern equipment, the designers take steps to keep the overall response flat, but there are those that do not.
The drawbacks to this approach are: 1) It requires the use of several “tuning” variables which often require the user to adjust and fiddle with in order to achieve good signal recognition. A fully automatic signal detection system should be able to choose values for these parameters without the intervention of an operator. 2) It does not take into account variations in the baseline noise floor that are introduced by RF energy in a live environment. Since these variations were not present during calibration, they are not part of the calibration vector and cannot be “canceled out” during the de-bias phase. Instead they remain during the square and detect phase, often being mistakenly classified as signal. An example of this is
In order to solve these two problems, and provide a fully automatic signal detection system, a new approach has been taken to prepare the calibration vector. The existing square and detect algorithm works well if the data are de-biased properly with a cleverly chosen calibration vector, it's just that the way we were creating the calibration vector was not sufficient.
So, rather than attempt to make the calibration vector fit the data, the new approach examines the data itself in an attempt to use parts of it as the correction vector. This is illustrated by the light purple and baby blue lines in the
The new Gradient Detection algorithm is applied to the smoothed data to detect locations where the slope of the line changes quickly. In places where the slope changes quickly in a positive direction, the algorithm marks the start of a signal. On the other side of the signal the gradient again changes quickly to become more horizontal. At that point the algorithm determines it is the end of a signal. A second smoothing pass is performed on the smoothed data, but this time, those values that fall between the proposed start and end of signal are left out of the average. The result is line 4 (baby blue) in
One of the other user-tunable parameters in the existing ASD system was called “Sensitivity”. This was a parameter that essentially set a threshold of energy, above which each FFT bin in a block of bins averaged together must exceed in order for that block of bins to be considered a signal. In this way, rather than a single horizontal line to divide signal from noise, each signal can be evaluated individually, based on its average power. The effect of setting this value too low was that tiny fluctuations of energy that are actually noise would sometimes appear to be signals. Setting the value too high would result in the algorithm missing a signal. In order to automatically choose a value for this parameter, the new system uses a “Quality of Service” feedback from the Event Compositor, a module that processes the real-time events from the ASD system and writes signal observations into a database. When the sensitivity value is too low, the random bits of energy that ASD mistakenly sees as signal are very transient. This is due to the random nature of noise. The Event Compositor has a parameter called a “Pre-Recognition Delay” that sets the minimum number of consecutive scans that it must see a signal in order for it to be considered a candidate for a signal observation database entry (in order to catch large fast signals, an exception is made for large transients that are either high in peak power, or in bandwidth). Since the random fluctuations seldom persist for more than 1 or 2 sweeps, the Event Compositor ignores them, essentially filtering them out. If there are a large number of these transients, the Event Compositor provides feedback to the ASD module to inform it that its sensitivity is too low. Likewise, if there are no transients at all, the feedback indicates the sensitivity is too high. Eventually, the system arrives at an optimal setting for the sensitivity parameter.
The result is a fully automated signal detection system that requires no user intervention or adjustment. The black brackets at the top of
Because the system relies heavily upon averaging, a new algorithm was created that performs an N sample average in fixed time: i.e. regardless of the width of the average, N, each bin requires 1 addition, 1 subtraction, and 1 division. A simpler algorithm would require N additions and 1 division per bin of data. A snippet of the code is probably the best description:
Automatic Signal Detection (ASD) with Temporal Feature Extraction (TFE)
The system in the present invention uses statistical learning techniques to observe and learn an RF environment over time and identify temporal features of the RF environment (e.g., signals) during a learning period.
A knowledge map is formed based on learning data from a learning period. Real-time signal events are detected by an ASD system and scrubbed against the knowledge map to determine if the real-time signal events are typical and expected for the environment, or if there is any event not typical nor expected.
The knowledge map consists of an array of normal distributions, where each distribution column is for each frequency bin of the FFT result set provided by a software defined radio (SDR). Each vertical column corresponds to a bell-shaped curve for that frequency. Each pixel represents a count of how many times that frequency was seen, detected, or observed at that power level.
A learning routine takes power levels of each frequency bin, uses the power levels as an index into each distribution column corresponding to each frequency bin, and increments the counter in a location corresponding to a power level.
The TFE function monitors its operation and produces a “settled percent.” The settled percent is the percentage of the values of the incoming FFT result set that the system has seen, detected, or observed before. In this way, the system can know if it is ready to interpret the statistical data that it has obtained. Once it reaches a point where most of the FFT values have been seen, detected, or observed before (99.95% or better), it can then perform an interpretation operation.
Gradient detection is then applied to the profile to identify areas of transition. An algorithm continues to accumulate a gradient value as long as the “step” from the previous cell to this cell is always non-zero and the same direction. When it arrives at a zero or different direction step, it evaluates the accumulated difference to see if it is significant, and if so, considers it a gradient. A transition is identified by a continuous change (from left to right) that exceeds the average range between the high and low boundaries of power levels shown as lines b (white) in
To a prior art receiver, the narrow band signal hidden within a wideband signal is not distinguishable or detectable. The systems and methods and devices of the present invention are operable to scan a wideband with high resolution or high definition to identify channel divisions within a wideband, and identify narrowband signals hidden within the wideband signal, which are not a part of the wideband signal itself, i.e., the narrow band signals are not part of the bundled channels within the wideband signal.
Also, at the red cursor in
The ASD system has the ability to distinguish between large eruptions of energy that increase the baseline noise and the narrow band signals that could normally be swamped by the additional energy because it generates its baseline from the spectrum itself and looks for relative gradients rather than absolute power levels. This baseline is then subtracted from the original spectrum data, revealing the signals, as displayed by the brackets at the top of the screen. Note that the narrow-band signals are still being detected (tiny brackets at the top that look more like dots) even though there is a hump of noise super-imposed on them.
TFE is a learning process that augments the ASD feature in the present invention. The ASD system enhanced with TFE function in the present invention can automatically tune parameters based on a segmented basis, the sensitivity within an area is changeable. The TFE process accumulates small differences over time and signals become more and more apparent. In one embodiment, the TFE takes 40 samples per second over a 5-minute interval. The ASD system in the present invention is capable of distinguishing signals based on gradients from a complex and moving noise floor without a fixed threshold bar when collecting data from an environment.
The ASD system with TFE function in the present invention is unmanned and water resistant. It runs automatically 24/7, even submerged in water.
The TFE is also capable of detecting interferences and intrusions. In the normal environment, the TFE settles, interprets and identifies signals. Because it has a statistical knowledge of the RF landscape, it can tell the difference between a low power, wide band signal that it normally sees and a new higher power narrow band signal that may be an intruder. This is because it “scrubs” each of the FFT bins of each event that the ASD system detects against its knowledge base. When it detects that a particular group of bins in a signal from ASD falls outside the statistical range that those frequencies normally are observed, the system can raise an anomaly report. The TFE is capable of learning new knowledge, which is never seen, detected, or observed before, from the signals identified by a normal detector. In one embodiment, a narrow band signal (e.g., a pit crew to car wireless signal) impinges on an LTE wideband signal, the narrow band signal may be right beside the wideband signal, or drift in and out of the wideband signal. On display, it just looks like an LTE wideband signal. For example, a narrow band signal with a bandwidth of 12 kHz or 25-30 kHz in a wideband signal with a bandwidth of 5 MHz over a 6 GHz spectrum just looks like a spike buried in the middle. But, because signals are characterized in real time against learned knowledge, the proposed ASD system with TFE function is able to pick out narrow band intruder immediately.
The present invention is able to detect a narrow band signal with a bandwidth from 1-2 kHz to 60 KHz inside a wideband signal (e.g., with a bandwidth of 5 MHz) across a 6 GHz spectrum. In
Statistical learning techniques are used for extracting temporal feature, creating a statistical knowledge map of what each frequency is and determining variations and thresholds and etc. The ASD system with TFE function in the present invention is capable of identifying, demodulating and decoding signals, both wideband and narrow band with high energy.
If a narrow band signal is close to the end of wideband LTE signal, the wideband LTE signal is distorted at the edge. If multiple narrowband signals are within a wideband signal, the top edge of the wideband signal is ragged as the narrow band signal is hidden within the wide band signal. If one narrow band signal is in the middle of a wideband signal, the narrow band signal is usually interpreted as a cell within the wideband signal. However, the ASD system with TFE function in the present invention learns power levels in a spectrum section over time, and is able to recognize the narrow band signal immediately.
The present invention is operable to log the result, display on a channel screen, notify operator and send alarms, etc. The present invention auto records spectrum, but does not record all the time. When a problem is identified, relevant information is auto recorded in high definition.
The ASD system with TFE in the present invention is used for spectrum management. The system in the present invention is set up in a normal environment and starts learning and stores at least one learning map in it. The learning function of the ASD system in the present invention can be enabled and disabled. When the ASD system is exposed to a stable environment and has learned what is normal in the environment, it will stop its learning process. The environment is periodically reevaluated. The learning map is updated at a predetermined timeframe. After a problem is detected, the learning map will also be updated.
The ASD system in the present invention can be deployed in stadiums, ports, airports, or on borders. In one embodiment, the ASD system learns and stores the knowledge in that environment. In another embodiment, the ASD system downloads prior knowledge and immediately displays it. In another embodiment, an ASD device can learn from other ASD devices globally.
In operation, the ASD system then collects real time data and compares to the learning map stored for signal identification. Signals identified by the ASD system with TFE function may be determined to be an error by an operator. In that situation, an operator can manually edit or erase the error, essentially “coaching” the learning system.
The systems and devices in the present invention create a channel plan based on user input, or external databases, and look for signals that are not there. Temporal Feature Extraction not only can define a channel plan based on what it learns from the environment, but it also “scrubs” each spectrum pass against the knowledge it has learned. This allows it to not only identify signals that violate a prescribed channel plan, but it can also discern the difference between a current signal, and the signal that it has previously seen, detected, or observed in that frequency location. If there is a narrow band interference signal where there typically is a wide band signal, the system will identify it as an anomaly because it does not match the pattern of what is usually in that space.
The device in the present invention is designed to be autonomous. It learns from the environment, and, without operator intervention, can detect anomalous signals that either were not there before, or have changed in power or bandwidth. Once detected, the device can send alerts by text or email and begin high resolution spectrum capture, or IQ capture of the signal of interest.
The identification and classification algorithms that the system uses to identify Temporal Features are optimized to be used in real time. Notice that, even though only fragments of the low level wide band signal are detected on each sweep, the system still matches them with the signal that it had identified during its learning phase.
Also as the system is running, it is scrubbing each spectral sweep against its knowledge map. When it finds coherent bundles of energy that are either in places that are usually quiet, or have higher power or bandwidth than it has seen, detected, or observed before, it can automatically send up a red flag. Since the system is doing this in Real Time, it has critical relevance to those in harm's way—the first responder, or the war fighter who absolutely must have clear channels of communication or instant situational awareness of imminent threats. It's one thing to geolocate a signal that the user has identified. It's an entirely different dimension when the system can identify the signal on its own before the user even realizes it's there. Because the device in the present invention can pick out these signals with a sensitivity that is far superior to a simple threshold system, the threat does not have to present an obvious presence to be detected and alerted.
Devices in prior art merely make it easy for a person to analyze spectral data, both in real time and historically, locally or remotely. But the device in the present invention operates as an extension of the person, performing the learning and analysis on its own, and even finding things that a human typically may miss.
The device in the present invention can easily capture signal identifications, match them to databases, store and upload historical data. Moreover, the device has intelligence and the ability to be more than a simple data storage and retrieval device. The device is a watchful eye in an RF environment, and a partner to an operator who is trying to manage, analyze, understand and operate in the RF environment.
The prior art is dependent upon a synchronized receiver for power, phase, frequency, angle, and time of arrival, and an accurate clock for timing, and significantly, requires three devices to be used, wherein all are synchronized and include directional antennae to identify a signal with the highest power. Advantageously, the present invention does not require synchronization of receivers in a multiplicity of devices to provide geolocation of at least one apparatus unit or device, thereby reducing cost and improving functionality of each of the at least one apparatus in the systems described hereinabove for the present invention. Also, the present invention provides for larger frequency range analysis, and provides database(s) for capturing events, patterns, times, power, phase, frequency, angle, and combinations for the at least one signal of interest in the RF spectrum. The present invention provides for better measurements and data of signal(s) with respect to time, frequency with respect to time, power with respect to time, and combinations thereof. In preferred embodiments of the at least one apparatus unit of the present invention, geolocation is provided automatically by the apparatus unit using at least one anchor point embedded within the system, by power measurements and transmission that provide for “known” environments of data. The known environments of data include measurements from the at least one anchorpoint that characterize the RF receiver of the apparatus unit or device. The known environments of data include a database including information from the FCC database and/or user-defined database, wherein the information from the FCC database includes at least maximum power based upon frequency, protocol, device type, and combinations thereof. With the geolocation function of the present invention, there is no requirement to synchronize receivers as with the prior art: the at least one anchorpoint and location of an apparatus unit provide the required information to automatically adjust to a first anchorpoint or to a second anchorpoint in the case of at least two anchorpoints, if the second anchorpoint is easier to adopt. The known environment data provide for expected spectrum and signal behavior as the reference point for the geolocation. Each apparatus unit or device includes at least one receiver for receiving RF spectrum and location information as described hereinabove. In the case of one receiver, it is operable with and switchable between antennae for receiving RF spectrum data and location data: in the case of two receivers, preferably each of the two receivers are housed within the apparatus unit or device. A frequency lock loop is used to determine if a signal is moving, by determining if there is a Doppler change for signals detected.
Location determination for geolocation is provided by determining a point (x, y) or Lat Lon from the at least three anchor locations (x1, y1): (x2, y2): (x3, y3) and signal measurements at either of the node or anchors. Signal measurements provide a system of non-linear equations that must be solved for (x, y) mathematically: and the measurements provide a set of geometric shapes which intersect at the node location for providing determination of the node.
For trilateration methods for providing observations to distances the following methods are used:
Also, in another embodiment of the present invention, a geolocation application software operable on a computer device or on a mobile communications device, such as by way of example and not limitation, a smartphone, is provided. Method steps are illustrated in the flow diagram shown in
The equations referenced in
Note that check if RF propagation distances form circles where one or more circles are Fully Enclosed if it is based upon Mod Type and Power Measured, then Set Distance 1 of enclosed circle to Distance 2 minus the distance between the two points. Also, next, check to see if some of the RF Propagation Distances Form Circles, if they do not intersect, then if so based on Mod type and Max RF power Set Distance to each circle to Distance of Circle+(Distance between circle points−Sum of the Distances)/2 is used. Note that deriving z component to convert back to known GPS lat lon coordinate is provided by: z=sqrt(Dist2−x2−y2).
Accounting for unknowns using Differential Received Signal Strength (DRSS) is provided by the following equation when reference or transmit power is unknown:
And where signal strength measurements in dBm are provided by the following:
For geolocation systems and methods of the present invention, preferably two or more devices or units are used to provide nodes. More preferably, three devices or units are used together or “joined” to achieve the geolocation results. Also preferably, at least three devices or units are provided. Software is provided and operable to enable a network-based method for transferring data between or among the at least two device or units, or more preferably at least three nodes, a database is provided having a database structure to receive input from the nodes (transferred data), and at least one processor coupled with memory to act on the database for performing calculations, transforming measured data and storing the measured data and statistical data associated with it: the database structure is further designed, constructed and configured to derive the geolocation of nodes from saved data and/or from real-time data that is measured by the units: also, the database and application of systems and methods of the present invention provide for geolocation of more than one node at a time. Additionally, software is operable to generate a visual representation of the geolocation of the nodes as a point on a map location.
Errors in measurements due to imperfect knowledge of the transmit power or antenna gain, measurement error due to signal fading (multipath), interference, thermal noise, no line of sight (NLOS) propagation error (shadowing effect), and/or unknown propagation model, are overcome using differential RSS measurements, which eliminate the need for transmit power knowledge, and can incorporate TDOA and FDOA techniques to help improve measurements. The systems and methods of the present invention are further operable to use statistical approximations to remove error causes from noise, timing and power measurements, multipath, and NLOS measurements. By way of example, the following methods are used for geolocation statistical approximations and variances; maximum likelihood (nearest neighbor or Kalman filter): least squares approximation; Bayesian filter if prior knowledge data is included: and the like. Also, TDOA and FDOA equations are derived to help solve inconsistencies in distance calculations. Several methods or combinations of these methods may be used with the present invention, since geolocation will be performed in different environments, including but not limited to indoor environments, outdoor environments, hybrid (stadium) environments, inner city environments, etc.
Typically, prior art arrays are more localized and deployed in a symmetrical fashion to reduce the complexity of mathematics and the equipment. The problem with localized fixed arrays are twofold: they require a large footprint for assembly and operation to gain accuracy in directional measurements. Conversely, smaller footprint arrays of geometric antenna systems can lose significant accuracy of the directional measurements. To avoid these limitations, a large variable array is used with fixed or mobile sites to allow greater accuracy.
In one embodiment of the present invention, geolocation using angle of arrival is provided by a fixed position antenna system constructed and configured with a four-pole array in a close proximity to each other. The antenna system is a unique combination of a half (½) Adcock antenna array positioned at each unit. The antenna system is fixed and is operable to be deployed with a switching device to a low-cost full Adcock system. The use of a phase difference on the dual receiver input allows the local unit to determine a hemisphere of influence in a full Adcock configuration or a group of the deployed units as a full space diversity Adcock antenna system. This embodiment advantageously functions to eliminate directions in the vector-based math calculation, thereby eliminating a large group of false positives.
The antenna system used with the geolocation systems and methods of the present invention includes three or more deployed units where none of the units is a full-time master nor slave. Each unit can be set to scan independently for target profiles. Once acquisition is obtained from one unit, the information is automatically disseminated to the other units within the cluster, i.e., the information is communicated wirelessly through a network. Preferably, the unit array is deployed in an asymmetrical configuration.
The antenna system in the present invention utilizes Normalized Earth Centered Earth Fixed vectors. Two additional vector attributes of the monitoring station are selected from the following: pitch, yaw, velocity, altitude (positive and negative) and acceleration.
Once a target acquisition from a single unit is acquired, a formatted message is broadcast to the deployed monitoring array stations. The formatted message includes but is not limited to the following: center frequency, bandwidth, modulation schema, average power and phase lock loop time adjustment from the local antenna system.
The monitoring units include a GPS receiver to aid in high resolution clocks for timing of signal processing and exact location of the monitoring unit. This is key to determine an exact location of the monitoring units, either fixed or mobile, to simulate mathematically the variable large scale antenna array. The phased-locked inputs determine the orientation of the incoming target signal into hemispheres of influence.
For this example,
The next step in the process is to determine for each target measurement the delays of arrival at each location. This will further reveal the direction of travel to the target or additionally if the target is within the large-scale variable array's own footprint.
Once the unit processing the data has received information from the other units in the array, processing of the information begins. First, the unit automatically sorts the array time of arrival at each location of the at least three units to construct mathematically a synthesis of the array. This is crucial to the efficiency and accuracy of the very large scale array, since no single monitoring unit is the point of reference. The point of reference is established by mathematical precedence involving time of arrival and the physical location of each monitoring unit at that point in time.
An aperture is synthesized between any two points on the array using the difference in the arrival time. Establishing a midpoint between two monitoring units establishes a locus for the bearing measurement along the synthesized aperture.
The aperture is given in radians by the following equation, where A is the wavelength in meters, and Distance is the arc length in meters.
Distance is calculated by the following equations, where R is the radius of the earth in kilometers, and Lat and Lon refer to the points on installation for latitude and longitude in radians.
The radial distance directly related to the angle of arrival across the aperture is given by the equation representing the radial time between monitor unit 1 and monitor unit 2 divided by Aperture Length:
Using fundamental logic, two possible angles of arrival between the units defining the synthetic aperture for a bearing from the midpoint as illustrated in
The use of a second component to establish a synthetic aperture yields another bearing as illustrated in
The foregoing method descriptions and the process flow diagrams are provided merely as illustrative examples and are not intended to require or imply that the steps of the various embodiments must be performed in the order presented. As will be appreciated by one of skill in the art the order of steps in the foregoing embodiments may be performed in any order. Words such as “thereafter,” “then,” “next,” etc. are not intended to limit the order of the steps: these words are simply used to guide the reader through the description of the methods. Further, any reference to claim elements in the singular, for example, using the articles “a,” “an” or “the” is not to be construed as limiting the element to the singular.
The various illustrative logical blocks, modules, circuits, and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both. To clearly illustrate this interchangeability of hardware and software, various illustrative components, blocks, modules, circuits, and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the overall system. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The hardware used to implement the various illustrative logics, logical blocks, modules, and circuits described in connection with the aspects disclosed herein may be implemented or performed with a general purpose processor, a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field programmable gate array (FPGA) or other programmable logic device, discrete gate or transistor logic, discrete hardware components, or any combination thereof designed to perform the functions described herein. A general-purpose processor may be a microprocessor, but, in the alternative, the processor may be any conventional processor, controller, microcontroller, or state machine. A processor may also be implemented as a combination of computing devices, e.g., a combination of a DSP and a microprocessor, a plurality of microprocessors, one or more microprocessors in conjunction with a DSP core, or any other such configuration. Alternatively, some steps or methods may be performed by circuitry that is specific to a given function.
In one or more exemplary aspects, the functions described may be implemented in hardware, software, firmware, or any combination thereof. If implemented in software, the functions may be stored as one or more instructions or code on a non-transitory computer-readable medium or non-transitory processor-readable medium. The steps of a method or algorithm disclosed herein may be embodied in a processor-executable software module which may reside on a non-transitory computer-readable or processor-readable storage medium. Non-transitory computer-readable or processor-readable storage media may be any storage media that may be accessed by a computer or a processor. By way of example but not limitation, such non-transitory computer-readable or processor-readable media may include RAM, ROM, EEPROM, FLASH memory, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other medium that may be used to store desired program code in the form of instructions or data structures and that may be accessed by a computer. Disk and disc, as used herein, includes compact disc (CD), laser disc, optical disc, digital versatile disc (DVD), floppy disk, and BLU-RAY disc where disks usually reproduce data magnetically, while discs reproduce data optically with lasers. Combinations of the above are also included within the scope of non-transitory computer-readable and processor-readable media. Additionally, the operations of a method or algorithm may reside as one or any combination or set of codes and/or instructions on a non-transitory processor-readable medium and/or computer-readable medium, which may be incorporated into a computer program product.
Certain modifications and improvements will occur to those skilled in the art upon a reading of the foregoing description. The above-mentioned examples are provided to serve the purpose of clarifying the aspects of the invention and it will be apparent to one skilled in the art that they do not serve to limit the scope of the invention. All modifications and improvements have been deleted herein for the sake of conciseness and readability but are properly within the scope of the present invention.
This application is related to and claims priority from the following U.S. patents and patent applications: this application is a continuation of U.S. application Ser. No. 18/374,395, filed Sep. 28, 2023, which is a continuation of U.S. application Ser. No. 18/201,295, filed May 24, 2023, which is a continuation of U.S. application Ser. No. 17/490,736, filed Sep. 30, 2021, which is a continuation of U.S. application Ser. No. 16/795,875, filed Feb. 20, 2020, which is a continuation of U.S. application Ser. No. 16/295,691, filed Mar. 7, 2019, which is a continuation of U.S. application Ser. No. 15/596,756, filed May 16, 2017, which is a continuation-in-part of U.S. application Ser. No. 15/412,982, filed Jan. 23, 2017, a continuation-in-part of U.S. application Ser. No. 14/940,299, filed Nov. 13, 2015, and a continuation-in-part of U.S. application Ser. No. 15/478,916, filed Apr. 4, 2017, U.S. application Ser. No. 14/940,299, filed Nov. 13, 2015, is a continuation of U.S. application Ser. No. 14/504,784, filed Oct. 2, 2014, which is a continuation of U.S. application Ser. No. 14/329,820, filed Jul. 11, 2014, which is a continuation of U.S. application Ser. No. 14/086,861, filed Nov. 21, 2013, which is a continuation-in-part of U.S. application Ser. No. 14/082,873, filed Nov. 18, 2013, which is a continuation of U.S. application Ser. No. 13/912,683, filed Jun. 7, 2013, which claims the benefit of U.S. application Ser. No. 61/789,758, filed Mar. 15, 2013, U.S. application Ser. No. 14/086,861, filed Nov. 21, 2013, is also a continuation-in-part of U.S. application Ser. No. 14/082,916, filed Nov. 18, 2013, which is a continuation of U.S. application Ser. No. 13/912,893, filed Jun. 7, 2013, which claims the benefit of U.S. application Ser. No. 61/789,758, filed Mar. 15, 2013, U.S. application Ser. No. 14/086,861, filed Nov. 21, 2013, is also a continuation-in-part of U.S. application Ser. No. 14/082,930, filed Nov. 18, 2013, which is a continuation of U.S. application Ser. No. 13/913,013, filed Jun. 7, 2013, which claims the benefit of U.S. application Ser. No. 61/789,758, filed Mar. 15, 2013, U.S. application Ser. No. 15/478,916, filed Apr. 4, 2017, is a continuation-in-part of U.S. application Ser. No. 14/934,808 filed Nov. 6, 2015, which is a continuation of U.S. application Ser. No. 14/504,836 filed Oct. 2, 2014, which is a continuation of U.S. application Ser. No. 14/331,706 filed Jul. 15, 2014, which is a continuation-in-part of U.S. application Ser. No. 14/086,875 filed Nov. 21, 2013, which is a continuation-in-part of U.S. application Ser. No. 14/082,873 filed Nov. 18, 2013, which is a continuation of U.S. application Ser. No. 13/912,683 filed Jun. 7, 2013, which claims the benefit of U.S. application Ser. No. 61/789,758 filed Mar. 15, 2013. U.S. application Ser. No. 14/086,875 is also a continuation-in-part of U.S. application Ser. No. 14/082,916 filed Nov. 18, 2013, which is a continuation of U.S. application Ser. No. 13/912,893 filed Jun. 7, 2013, which claims the benefit of U.S. application Ser. No. 61/789,758 filed Mar. 15, 2013. U.S. application Ser. No. 14/086,875 is also a continuation-in-part of U.S. application Ser. No. 14/082,930 filed Nov. 18, 2013, which is a continuation of U.S. application Ser. No. 13/913,013 filed Jun. 7, 2013, which claims the benefit of U.S. application Ser. No. 61/789,758 filed Mar. 15, 2013. Each of the U.S. Applications mentioned above is incorporated herein by reference in its entirety.
Number | Date | Country | |
---|---|---|---|
61789758 | Mar 2013 | US | |
61789758 | Mar 2013 | US | |
61789758 | Mar 2013 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 18374395 | Sep 2023 | US |
Child | 18432739 | US | |
Parent | 18201295 | May 2023 | US |
Child | 18374395 | US | |
Parent | 17490736 | Sep 2021 | US |
Child | 18201295 | US | |
Parent | 16795875 | Feb 2020 | US |
Child | 17490736 | US | |
Parent | 16295691 | Mar 2019 | US |
Child | 16795875 | US | |
Parent | 15596756 | May 2017 | US |
Child | 16295691 | US | |
Parent | 14504836 | Oct 2014 | US |
Child | 14934808 | US | |
Parent | 14331706 | Jul 2014 | US |
Child | 14504836 | US | |
Parent | 13912683 | Jun 2013 | US |
Child | 14082873 | US | |
Parent | 13912893 | Jun 2013 | US |
Child | 14082916 | US | |
Parent | 13913013 | Jun 2013 | US |
Child | 14082930 | US | |
Parent | 14504784 | Oct 2014 | US |
Child | 14940299 | US | |
Parent | 14329820 | Jul 2014 | US |
Child | 14504784 | US | |
Parent | 14086861 | Nov 2013 | US |
Child | 14329820 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 15478916 | Apr 2017 | US |
Child | 15596756 | US | |
Parent | 14934808 | Nov 2015 | US |
Child | 15478916 | US | |
Parent | 14086875 | Nov 2013 | US |
Child | 14331706 | US | |
Parent | 14082873 | Nov 2013 | US |
Child | 14086875 | US | |
Parent | 14082916 | Nov 2013 | US |
Child | 14086875 | US | |
Parent | 14082930 | Nov 2013 | US |
Child | 14086875 | US | |
Parent | 15412982 | Jan 2017 | US |
Child | 15596756 | US | |
Parent | 14940299 | Nov 2015 | US |
Child | 15412982 | US | |
Parent | 14082873 | Nov 2013 | US |
Child | 14086861 | US | |
Parent | 14082916 | Nov 2013 | US |
Child | 14082873 | US | |
Parent | 14082930 | Nov 2013 | US |
Child | 14082916 | US |