US20030228017A1 - Method and system for waveform independent covert communications - Google Patents

Method and system for waveform independent covert communications Download PDF

Info

Publication number
US20030228017A1
US20030228017A1 US10/400,486 US40048603A US2003228017A1 US 20030228017 A1 US20030228017 A1 US 20030228017A1 US 40048603 A US40048603 A US 40048603A US 2003228017 A1 US2003228017 A1 US 2003228017A1
Authority
US
United States
Prior art keywords
signal
symbols
communication
eigenvalue
noise
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Abandoned
Application number
US10/400,486
Inventor
Edward Beadle
Richard Anderson
John Dishman
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harris Corp
Original Assignee
Beadle Edward Ray
Anderson Richard Hugh
Dishman John Fitzgerald
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beadle Edward Ray, Anderson Richard Hugh, Dishman John Fitzgerald filed Critical Beadle Edward Ray
Priority to US10/400,486 priority Critical patent/US20030228017A1/en
Publication of US20030228017A1 publication Critical patent/US20030228017A1/en
Assigned to HARRIS CORPORATION reassignment HARRIS CORPORATION ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: BEADLE, EDWARD RAY, DISHMAN, JOHN FITZGERALD, ANDERSON, RICHARD HUGH
Abandoned legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/74Multi-channel systems specially adapted for direction-finding, i.e. having a single antenna system capable of giving simultaneous indications of the directions of different signals
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2134Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on separation criteria, e.g. independent component analysis
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04KSECRET COMMUNICATION; JAMMING OF COMMUNICATION
    • H04K1/00Secret communication
    • H04K1/02Secret communication by adding a second signal to make the desired signal unintelligible
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/22Source localisation; Inverse modelling
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7097Interference-related aspects
    • H04B1/7103Interference-related aspects the interference being multiple access interference
    • H04B1/7105Joint detection techniques, e.g. linear detectors
    • H04B1/71052Joint detection techniques, e.g. linear detectors using decorrelation matrix

Definitions

  • a message is desirably sent covertly from a mobile or fixed transmitter to a remote receiver or, more importantly, a friendly or cooperative receiver with out the presence of a signal being detected by an unintended receiver.
  • Covertly transmitting a signal may be necessary to reflect the fact that the transmission is conducted in such a way as to avoid intentional or inadvertent detection by systems monitoring the electromagnetic (EM) environment local to the covert transmitter.
  • EM electromagnetic
  • any transmission, covert or not may be desired to be received only by intended receivers.
  • Intentional detection of the signal or message can be accomplished in military systems that use specially designed electronic support measures (ESM) receivers. These ESM receivers are often found in signal intelligence (SIGINT) applications.
  • SIGINT signal intelligence
  • devices employed by service providers i.e. spectral monitors, error rate testers
  • IGFINT signal intelligence
  • devices employed by service providers i.e. spectral monitors, error rate testers
  • Inadvertent detection can also occur, such as when a user or service provider notices degradation in link performance (e.g., video quality, audio quality, or increased bit error rate).
  • covert also implies the additional goals of evading interception and exploitation by unintended receivers.
  • Interception is the measurement of waveform features or parameters useful for classifying/identifying a transmitter and/or the waveform type and/or deriving information useful for denying (i.e. jamming) the communication.
  • Exploitation is processing a signal by an unintended receiver in the attempt to locate the transmitter and/or recover the message content.
  • LPD low probability of detection
  • LPI low probability of intercept
  • LPE low probability of exploitation
  • an LPD communication system shown in FIG. 1 is designed by minimizing the radiated energy in all directions.
  • the goal is to have the emitter 102 radiate zero energy in all directions, save for the direction of the intended receiver 103 .
  • This approach limits the geographic locations in which the signal would likely be detected.
  • the area 110 receives little radiation and thus the unintended receiver 105 would have greater difficulty detecting the signal than the unintended receiver 106 in area 111 .
  • LPD systems use only enough energy in the direction of the intended receiver to be reliably detected and the message recovered by the intended receiver.
  • the radiation in non-preferred directions is primarily controlled by the aperture (i.e. antenna system) and additional performance improvement for LPD can be gained by employing a directional antenna.
  • the aperture i.e. antenna system
  • An unintended receiver such as the receiver 103 in FIG. 1 with the goal of detecting a covert communication must reliably differentiate between the binary noise-only and signal-plus-noise hypothesis.
  • the signal detection process is typically based on an energy threshold.
  • the power P avg is the received covert signal power S plus internal receiver noise power N.
  • E tot (S+N)T xmit .
  • the bandwidth of the message signal is artificially broadened to reduce the energy loading on the band, as is done in direct sequence spread spectrum (DSSS) communication, the intended receiver has a processing gain over the unintended receivers that use the same co-channel receiver bandwidth to capture the signal.
  • DSSS direct sequence spread spectrum
  • the intended receiver 103 must have knowledge of the spreading code and be able to acquire it and remain synchronized to the covert transmitter code, in addition to the usual communication receiver functions of carrier and timing recovery. This requirement, for some applications, is not a problem but for other applications, covert message transmission functionality is greatly diminished.
  • Minimizing transmit power has two direct system benefits. First, the total signal power used will be a small fraction of the total noise power present in the same band. Thus, if the message is limited in time duration, the total energy measured by an unintended receiver 106 , which may be an ESM receiver, is indistinguishable from a noise-only environment. Since ESM receivers are often of energy threshold type, there is an obvious trade-off of average power for time duration in order for a signal to remain undetectable. Second, the low transmit power scenario enables usage by transmitters with very limited power supplies (i.e. batteries).
  • Embodiments of the present inventive system and method address the above needs while requiring only an extremely low power signal.
  • FIG. 1 is a general representation of a prior art approach to LPD.
  • FIG. 2 depiction of covert scenarios.
  • FIG. 3 is a representation of an embodiment of a waveform independent covert communication system.
  • FIG. 4 is a depiction of a binary symbol message 101011 according to an embodiment of the invention.
  • FIG. 5 a is a representation of a spectrogram of a transmission of the binary symbol message depicted in FIG. 4 and revealing the symbol boundaries of the message according to an embodiment of the invention.
  • FIG. 5 b is a representation of a spectrogram of a transmission of the binary symbol message depicted in FIG. 4 as might been seem by a typical ESM/SIGINT spectral monitor.
  • FIG. 6 a is a generalized schematic representation of embodiments of the invention.
  • FIG. 6 b is a schematic representation of an embodiment of the invention with dual signal sources.
  • FIG. 6 c is a schematic representation of an embodiment of the invention with dual temporal filters.
  • FIG. 6 d is a schematic representation of an embodiment of the invention with dual non-Gaussian noise generators.
  • FIG. 6 e is a schematic representation of an embodiment of the invention with dual non-Gaussian noise generator and temporal filter branches.
  • FIG. 6 f is a schematic representation of an embodiment of the invention with dual temporal filters and dual signal sources.
  • FIG. 7 is a schematic representation of a Laplacian noise generator with multiple signal sources according to an embodiment of the invention.
  • FIG. 8 is a representation of eigenvalue tracks for a frequency shifter laplacian noise waveform for message “101011” according to an embodiment of the invention.
  • FIG. 9 is a representation of block-to-block eigenvalue correlations.
  • FIG. 10 is a schematic representation of a message recovery system with spatial information according to an embodiment of the invention.
  • FIG. 11 is a flow diagram for covert communication via message recovery with spatial information according to an embodiment of the invention.
  • FIG. 12 is a schematic representation of a message recovery system without spatial information according to an embodiment of the invention.
  • FIG. 13 is a flow diagram for covert communication via message recovery without spatial information according to an embodiment of the invention.
  • FIG. 14 is a flow diagram for covert communication via encoded eigenvalues according to an embodiment of the invention.
  • FIG. 15 is a representation of a binary message sequence 101011 encoded using different carrier waveforms.
  • a useful feature of embodiments described herein is the use of time duration to convey information, while the individual waveforms used to convey the message (in the time duration of the fourth-order characteristic) are in a sense superfluous or independent of the message to be conveyed.
  • This is a significant advantage for ubiquitous application, allowing for parasitic use of present communication infrastructure and devices.
  • covert communication methods and apparatus because of the independence of the information transfer on the “carrier waveform”.
  • embodiments of the present invention by contrast can function equally well for any waveform, and the receiver does not require any a priori knowledge of the “carrier waveform”.
  • embodiments of the covert transmitter can be waveform agile without informing the intended receiver.
  • the embodiments herein are predicated on selecting and transmitting carrier waveforms with unique higher order spatial statistics.
  • Such higher order statistics include 2 nd order spatial moments and 4 th order spatial cumulants.
  • the number of type of encoding schemes mapping an M-ary alphabet to time durations is limitless, but in point applications there will be preferred methods since they ease hardware and software realization of the system.
  • Specific embodiments are identified herein. One method is by design, such as using frequency shifted Laplacian noise.
  • Another method is to use the inherent properties of conventional waveforms (e.g., BPSK, QPSK, GMSK, QAM, DBPSK MFSK, FSK, and DQPSK). Yet another is to alter channel filtering in the transmitter. Common to all these approaches is that a high order statistic, specifically, but not limited to, the fourth-order characteristic, of the transmitted signal is recoverable with a spatial fourth order cumulant matrix pencil (SFOCMP). Furthermore, the message symbols map to time durations of these higher order statistics and not to specific values of the SFOCMP. Hence, from the receiver's point of view, the communication is independent of waveform selection by the transmitter. The primary restriction is that the receiver and transmitter must use the same “codebook” of time durations and alphabet. However, the “codebook” need not be known a priori by the receiver in all instances. As shown herein, the receiver will be able to discern the different time durations and be able to extract some useful information therefrom.
  • BPSK, QPSK e.g.,
  • the waveform is not conventionally demodulated. Rather, a straightforward block or batch estimation algorithm estimates the generalized eigenvalues of the SFOCMP for each signal in the receiver field-of-view (FOV) over time using only the array output.
  • FOV receiver field-of-view
  • the sizing for the block processing e.g., the block of contiguous array observations, sometimes known as “snapshots” is dependent on several factors. Chiefly we must ensure that each block has enough sample support so that the eigenvalue estimates from the GEVD of the SFOCMP in each block over a symbol duration nominally match. This means that estimation error is negligible.
  • changes in the eigenstructure can be reliably detected, and this change indicates a symbol boundary.
  • the degree to which a nominal match is required within a block depends on the complexity of the signal environment (e.g., extraneous co-channel signals), the communication errors (e.g., partially received messages) tolerable in a given application, and the receiver processing resources to recover the message in a timely manner.
  • the source 300 produces a digital data stream of a covert message to be transmitted by the covert transmitter 302 .
  • the digital data can either be inherently digital (e.g., Morse coded messages) or analog (e.g., voice) and subsequently digitized.
  • the digital data can also include framing or formatting of the message. Typical of the framing would be start/stop and data fields. Though other fields can be used as needs dictate. This framing structure can ensure that the receiver can reliably find the beginning of a message for synchronization. Source coding or compression could also be applied to the incoming data stream to reduce the required bandwidth.
  • the user may also encrypt the data to protect it by an optional encryption device of suitable complexity.
  • This data is then output to the forward error correction (FEC) module, which currently is envisioned as applying block coding.
  • FEC forward error correction
  • the coding is useful to aid the receiver in resolving message ambiguities say caused by fades or other unresolved time coincident measurements, which in this system could be processed as “erasures” up to the correction capability of the code.
  • the error correction code would find application when the symbol stream is improperly decoded due to statistical variations in the time duration measurement process.
  • the potential for a message protocol using automated repetition of the message might be advantageous as error patterns in each transmission will likely be different. This type of protocol would require the framing as stated above, and possibly the inclusion of other fields such as the number of times the message has been or will be repeated.
  • Message formatting, compression, encryption and/or FEC is performed by the processor 304 .
  • the processing in processor 304 is not necessary for practicing the present invention but is included to indicate the ability to add additional robustness and reliability to the transmission of a message.
  • the binary symbol stream is from source 300 is then mapped into suitably selected M-ary symbols (or time durations) in the symbol mapping device 306 , via look-up tables or other codebook like mapping mechanisms.
  • the alphabet size is limited only by system implementation and design considerations for particular applications.
  • the encoded data is then applied to the “carrier waveform” by carrier waveform generator 308 which is a signal generator, implementations of which are further described below.
  • carrier waveform generator 308 which is a signal generator, implementations of which are further described below.
  • the “carrier waveform” has its fourth-order characteristic modified according to the control by the M-ary alphabet. There is no limit to the strategies potentially adopted by the covert transmitter for this operation, so long as the characteristic is measurable reliably by the receiver and it conforms to the time duration and alphabet size assumed by the receiver.
  • the baseband data from the carrier waveform generator 308 is fed to a suitable RF subsystem 310 and antenna 312 for transmission.
  • the minimum duration and duration increment must be such that synchronizing the data block boundaries used in the receiver to that of the symbol timing in the transmitter 302 is not relied upon. It is desirable that the covert transmitter 302 use a fundamental signaling period of several (e.g., 10) “receiver block” durations for the minimum signal, and may have a signaling duration increment of the same size to define the alphabet. However other choices are applicable depending on the particular implementation and application. The goal is to make the time duration alphabet as disparate as possible while meeting performance objectives (e.g., data throughput). Sample data rate computations can be determined as shown below.
  • Achievable data rates are in principle limited by operating conditions (received SNR, tolerable BER, cumulant estimation errors, etc.).
  • the incremental time duration for the alphabet should be at least ES(1/R), where E is the number of blocks desired by the designer to provide a balance between adequate safety margin in the time duration decision process and required throughput rate.
  • E can be as low as unity which would enhance the achievable signaling rate for a fixed alphabet size.
  • the receiver 303 uses an N-element (or port) receive array 327 and an RF processor 305 to obtain the transmitted signal.
  • the array data is first sampled and digitized at some rate suitable for the application. Each array output is digitized simultaneously producing a vector observation in the vector digitizer and buffer 307 .
  • the array output data is buffered and subdivided into non-overlapping blocks in 307 . Block-wise across signal samples (i.e.
  • the operation of the BSS requires the selection of a triplicate of time lags provided by the time lags selection device 311 .
  • the superscript b is used as a block counter in the receiver.
  • the two fourth-order spatial cumulant matrices required to form the SFOCMP are formed using pre-selected delay triplets.
  • the delays can be either pre-selected, or subjected to online modification using a programmed search routine (if necessary). This search routine might be necessary when certain conditions, such as repeated eigenvalues for different signals are encountered. Processing of the received covert message this difficulty is likely of no consequence.
  • provisions are made for signals whose eigenstructure match at the delays selected to be repressed at different delays to provide improved discrimination if desirable.
  • the GEVD is computed. From the GEVD, the eigenvalues and eigenvectors are used to determine the signal environment over time block b. Subsequently, the eigenvectors are used to determine the signal steering vectors and then the eigenstructure is correlated block wise in the Blockwise Eigenvalue Correlator 313 to determine any changes in the signal environment. A change, such as symbol boundary, in the communication signal will alter its contribution to the signal environment eigenstructure, measured by the SFOCMP, in a detectable manner. This means a “significant” movement in the complex plane of eigenvalues. As signal changes are detected, those signals are cued for storage in the signal history database 317 .
  • the spatial location of the signal is determined (i.e., either angle-of-arrival (AOA) or geolocation, whichever applies given the specific application) by the AOA/Geolocation processor 315 .
  • the steering vector can be recorded, which is useful when refined spatial information is unavailable and the relative motion of the transmitter and receiver is negligible.
  • the eigenvalues no longer correlating with the present signal structure are also written to the database.
  • the temporal support (i.e. durations) of the eigenvalues no longer correlating with the current signal structure is measured and stored. All this data is formed and recorded in the signal history database 317 along with other ancillary data that may be useful for signal post-processing applications such as data mining or covert message recovery.
  • An important function of a tracker is the track initiation and deletion logic.
  • An embodiment of the tracks uses a fixed distance and a fixed number of consecutive “good associations” for initiation and a single “no association” for a track deletion.
  • a “good association” is any measurement that is “close enough” to track.
  • a “no association” condition occurs when all the measurements are “Too far” from a particular track. the distance indicative of an good association can be set empirically or experimentally.
  • Track initiation and track deletion strategies can also be used to adapt to various situation.
  • a Kalman-like approach to association gates can be adapted as the number of observation for a track are accumulated. Such an approach also has the advantage of replacing fixed averaging of the measurements.
  • the covert communication receive processor 319 To recover the covert message, in the covert communication receive processor 319 , first the time durations are time-gated and those passing the time-gate are then mapped to the M-ary symbols of the “codebook” by the Symbol Demapper 321 .
  • the Message Sequence Correlator 323 is an association engine that “stitches” all the signal eigenvalues durations and symbol values corresponding to a “consistent” spatial location into a time-ordered stream or signal track.
  • the data is supplied to the optional data processor 325 , where the data can then be FEC decoded, decrypted (if necessary), decompressed (if necessary) and then the message can be checked for synchronization, as well as, the data provided to the sink 301 with the appropriate format.
  • the design also allows for multiple access for communications.
  • multiple remote covert emitters are sending data. It is unlikely that they would have exactly the same fourth-order cumulant representation, even if they are using the same base waveform. This is because any deviation from nominal waveform implementation (e.g., frequency change, waveform change, matrix pencil eigenvalue change, phase noise, I/Q imbalance, timing jitter, phase jitter, symbol rate change, pulse shape change, a fourth-order statistic change, relative rotational alignment of a signal constellation change, power amplifier rise/fall time change, and Doppler shift change) causes the fourth-order statistics of these signals to differ. Further, the multiple access signals are assumed distinguishable by spatial location.
  • FIG. 4 is a representation of a binary message sequence 101011 where the symbols are discriminated by duration.
  • the binary sequence is mapped to binary time durations. It shall be readily appreciated by those of skill in the art that the process may be extended to M-ary time durations.
  • the carrier frequency for the embodiment shown is offset at the time boundaries, thus shifting the fourth-order cumulant (i.e., the eigenvalue) for that symbol. Ultimately the time duration of the eigenvalue is recovered to reconstruct the message.
  • the frequency shift for the repeated “1” symbol at the end of is needed (in this case) to differentiate a repeated “1” from a single “0” and visa versa in the case of a double “0” (not shown).
  • the temporal correlation of the SFOCMP eigenstructure must be broken at a message symbol boundary. This break of temporal correlation in FIG. 5 a is accomplished by altering the fourth-order cumulant with a frequency shift.
  • FIG. 5 a is a spectrogram of the emitter message shown in FIG. 4 when using the frequency shifts of the preferred embodiment discussed herein. The symbol boundaries defined by frequency offsets at the transmitter are clearly evident in the spectrogram.
  • FIG. 5 b is a representation of a typical Electronic Support Measures (ESM) receiver/Signal Intelligence Spectral Monitor output viewing the received signal of FIG. 5 a.
  • ESM Electronic Support Measures
  • FIG. 5 b is a representation of a typical Electronic Support Measures (ESM) receiver/Signal Intelligence Spectral Monitor output viewing the received signal of FIG. 5 a.
  • ESM Electronic Support Measures
  • FIG. 5 b the time varying spectral pattern of the message shown in 5 a is buried in interference and noise (SNIR ⁇ 10 db) making the covert message very difficult to detect, and thus even more difficult to intercept and exploit.
  • SNIR stands for Signal to Noise and Interference Ratio.
  • a mathematical element of the invention is the use of spatial high order statistics to separate signal sources, such as a blind source separation algorithm that utilizes a normalized spatial fourth-order cumulant matrix pencil and its generalized eigenvalue decomposition (GEVD).
  • GSVD generalized eigenvalue decomposition
  • the equations presented herein use the following subscripting convention. Quantities relating to the array observations available to the system are denoted with a boldface subscript x. However, the subscript should not be confused with the representation of the vector observation from the array output, also denoted as a boldface x. From the context the meanings shall be clear to those of skill in the art. Further, quantities relating to the propagating signals impinging on a receive array are denoted with a boldface subscript r.
  • C x 4 represent the spatial fourth-order autocumulant matrices.
  • the arguments of the terms indicate the triplet of time delays used to form the matrices.
  • the matrix is N ⁇ N
  • the subscript rc indicates the element in the r th row and the c th column.
  • the quantity V shown in equation 2 is a N ⁇ M s matrix composed of the steering vectors for each signal impinging on the array, where N is the number of array ports available to the user and M s , M s ⁇ N, is the number of signals.
  • the well-known array propagation vector is a steering vector (i.e., the time delay is represented as phase).
  • the matrix V is of full rank. This guarantees an equivalence between the eigen structure of the pencils P r and P x .
  • [0059] contributes one finite eigenvalue, and it is expressed as the inverse normalized fourth-order autocumulant for that signal as expressed by equation 3.
  • ⁇ m c r m 4 ⁇ ( 0,0,0 )
  • c r m 4 ⁇ ( ⁇ 1 , ⁇ 2 , ⁇ 3 ) ⁇ ⁇ for ⁇ ⁇ m 1 , 2 , ... ⁇ , M ( 3 )
  • [0061] represent the individual fourth-order cumulant terms for each signal. These terms are actually the diagonal terms of the pencil P r as shown in equation (4).
  • P r ⁇ ( ⁇ , ⁇ ) [ c r 1 4 ⁇ ( 0,0,0 ) - ⁇ ⁇ ⁇ c r 1 4 ⁇ ( ⁇ 1 , ⁇ 2 , ⁇ 3 ) 0 ⁇ ⁇ ⁇ ⁇ 0 ⁇ ⁇ ⁇ c j 4 ⁇ ( 0,0,0 ) - ⁇ ⁇ ⁇ c j 4 ⁇ ( ⁇ 1 , ⁇ 2 , ⁇ 3 ) ⁇ ⁇ ⁇ ⁇ ⁇ 0 ⁇ ⁇ c r M 4 ⁇ ( 0,0,0 ) - ⁇ ⁇ ⁇ c r M 4 ⁇ ( 0,0,0 ) - ⁇ ⁇ ⁇ c r M 4 ⁇ ( ⁇ 1 , ⁇ 2 , ⁇ 3 ) ] (
  • eigenvalues are available to an analysis system, and in theory are independent of system Gaussian noise level given sufficient length data records.
  • the eigenvalues are implicit characteristics of the signals carrying the emitter's covert message in each symbol duration.
  • the receiver will typically form blocks or batches of received data for the purpose of correlating the eigenstructure over time to determine patterns of persistent values (FIG. 3) augmented by the availability of spatial data. It is important to recall that only the time duration of the emitter's statistical characteristic as measured by the SFOCMP is relevant, and not the exact values. Hence, the emitter is completely free to choose the “carrier waveforms” at will (FIGS. 6 a - f ).
  • the steering vectors can be estimated from the cumulant data for each signal in the FOV of the receiver.
  • the Blind Source Separation processor 309 forms and applies a separation Matrix and enumerates the number of sources. As described above, from an array output the spatial 4 th order cumulant matrices are estimated and the estimates are used to determine the eigen analysis for the first-order matrix pencil. Signal detection and enumeration providing the number of source is performed and the separation matrix from the pencil eigenvectors is accomplish. The separation matrix is also used to determine spatial variables, such as steering vectors, AoA or geolocation in block 315 of FIG. 3.
  • covert transmitter can control the desired characteristic of its emitted waveform, some of which leave lower-order statistics unchanged. For example altering the channel filter (i.e. Nyquist pulse shaping) between maximum phase and minimum phase realizations is undetectable in the second-order domain (i.e. power spectra), but evident in the fourth-order domain as measured by our SFOCMP. Also one could conceive on signaling with kurtosis, a fourth-order statistic, applied in the transmitter. Or, one could simply shift between classic waveforms, for example, BPSK, GMSK, QPSK, QAM, or potentially even just variants (i.e. constellation rotations, different pulse shape filters) of a fixed modulation type. There also exists the possibility that the “carrier waveforms” might be chosen as chaotic to appear more noise-like or designed using numerical techniques and generated using direct synthesis in a transmitter.
  • the information transfer from the transmitter to the receiver is contained in the duration of the change in the eigenvalues of the SFOCMP, not the particular eigenvalues. Since our technique is independent of the particular eigenvalue, it is independent of the waveforms used by the emitter. Which allows, in principle, any transmitter to make use of the receiver having the capability to exploit fourth-order cumulants. The degree to which a specific emitter wishes to “hide” from say conventional ESM receivers holds the implications on the implementation details of the “carrier waveforms”.
  • FIG. 6 a is a generalized schematic diagram of an embodiment of a noise signal generator for a waveform independent covert communication system.
  • a noise generator 602 generates temporal dependent non-Gaussian noise.
  • the output of the noise generator 602 is combined by combiner 606 with a carrier signal or waveform source 604 .
  • the combiner can be an adder, a mixer, a multiplier, a non-linear device or other type of combiner that facilitate a change in a higher order signal statistic. This list of combiners in not exhaustive and should not be construed to limit the scope of the invention.
  • the combined signal forms a type of baseband signal that is processed (eg. amplified, upconverted, etc.) and transmitted by transmitter 608 .
  • FIG. 6 a form a basic framework and facilitates several different embodiment for noise signal generation and transmission of a covert message.
  • FIG. 6 b is a specific implementation of FIG. 6 a with dual signal sources.
  • the noise generator 602 is implemented with a Non-Gaussian Noise Generator (NGNG) 610 and a temporal filter 612 that creates an output with temporal dependence.
  • the signal source 604 is implemented with two unique signal sources 614 a and 614 b.
  • Unique signal sources 614 a and 614 b are selected to create unique spatial high order statistics when combined by the combiner 606 with the output of the noise generator 602 .
  • the unique signal sources 614 a and 614 b are connected to the combiner by a switch 618 .
  • the switch is driven by the symbol bit stream of the covert message.
  • the switch 618 connects alternate unique signal sources at the conclusion of each successive symbol duration.
  • FIG. 6 c is another specific implementation of FIG. 6 a having dual temporal filters.
  • the noise generator 602 is implement with a non-Gaussian noise generator 610 connected to two unique temporal filters 612 a and 612 b.
  • the unique temporal filters 612 a and 612 b are selected to create unique spatial high order statistics of the transmitted combined signal (e.g. maximum and minimum phase filters).
  • the unique temporal filters 612 a and 612 b are connected by switch 618 to the combiner 606 and combined with the output of the signal source 604 . Again the switch is driven by the symbol bit stream of the covert message in the same manner described above for the dual signal source implementation.
  • the signal or waveform transmitted from transmitter 608 of FIG. 6 c provides alternating unique spatial high order statistics for each successive message symbol independent of the symbol transmitter.
  • FIG. 6 d is a specific implementation of the noise generator 602 including two unique non Gaussian noise generators 610 a and 610 b connected by switch 618 to the temporal filter 612 . Again the output of the noise generator combined with the signal source 604 in combiner 606 and transmitted by emitter or transmitter 608 .
  • the unique NGNGs 610 a and 610 b are likewise selected for their effect on spatial high order statistic of the transmitter signal.
  • FIG. 6 e is an implementation similar to the noise generator 602 of FIG. 6 d , where a non-Gaussian noise generator 610 a and temporal filter 612 a form a unique branch and NGNG 610 b and temporal filter 612 b form another unique branch.
  • the branches are unique in the sense of their respective spatial-high order statistics.
  • the branches are connected to the combiner 606 by switch 618 .
  • FIG. 6F is but one of the many possible combinations.
  • the implementation of the noise generator 602 of FIG. 6 f includes the noise generator implementation of FIG. 6 d containing two unique NGNGs 610 a and 610 b coupled with the signal source 604 of FIG. 6 b containing two unique signal sources 614 a and 614 b. While the switches 618 a and 618 b are constrained to switching at only at the duration boundaries, their operation advantageously would operate independently. The result of this specific implementation allows for switching between four unique signals.
  • FIG. 7 is an embodiment using a Laplacian generator 710 and an infinite impulse response low pass filter 712 (e.g. an autoregressive single-pole filter) as the noise generator 602 .
  • the Laplacian generator 710 is chosen for leptokurtosis property. We have chosen to form the necessary sequence of random Laplacian variates by selecting a pair of random variants, say x and y, uniformly distributed between 0 and 1 and forming log(x/y) for each sample required. For complex baseband symbols a total of four uniformly distributed variates must be sampled providing a stream of complex samples where the real and imaginary parts of each are independently Laplacian distributed.
  • the output of the noise generator 602 is combined by combiner 606 with the signal source 604 .
  • the signal source 604 is implemented with multiple unique signal sources 714 1 , 714 2 , . . . 714 1 that are used to switch spatial high order statistics (i.e. matrix pencil eignevalues for the present embodiment) for each symbol of the covert messages.
  • the unique signal sources 714 shown are offset in frequency to provide unique matrix pencil eigenvalues. In this particular embodiment, L must be greater than or equal to the alphabet size M.
  • the basic waveform need not be altered from QPSK to say BPSK, but uniqueness can be achieve by controlling the parameters of, for example, timing jitter, “carrier waveform” symbol rate, phase noise, pulse shaping and the like to provide the necessary separation in the fourth-order domain or other high order domains.
  • the “carrier waveforms” can contain completely worthless symbols, hence anyone trying to conventionally exploit the “carrier waveform” information would be wasting their time, alternatively the “carrier waveforms” may carry information with a secondary message.
  • a secondary message contained within the carrier waveform would not necessarily be LPE, LPI or LPD.
  • FIG. 8 An example of a potential message recovery embodiment is shown in FIG. 8.
  • the receiver requires three parameters relative to the time durations of the message symbols. The receiver must know the minimum symbol duration 801 , maximum symbol duration 803 , and some other information regarding the other M-2 symbols in the alphabet. In the embodiment shown in FIG. 8 only a binary symbol alphabet is used. This information can be provided, say by policy where the symbol duration increments are integer-related to the minimum symbol duration 801 . In this case the receiver need only know the maximum duration, the minimum duration and alphabet sizes. Otherwise the receiver can just know a-priori the pre-selected durations which can then be arbitrarily selected.
  • the symbol durations should be selected in an application for easing the decision process in the receiver to map the measured durations to one of the alphabet elements. In this case we used a symbol twice as long as another, but such widely disparate time durations may impact throughput.
  • a minimum symbol duration threshold 807 and a maximum symbol duration threshold 811 can be used as a time gate to remove spurious or interfering tracks.
  • Spatially correlated signal tracks 823 , 825 , 827 and 829 are examples of eigenvalue durations which exceed the threshold in the time gate and thus are ignored or discarded.
  • Spatially correlated eigenvalue time durations 805 9 , 805 16 , 805 36 , 805 53 , 805 68 and 805 80 are sequenced in signal track 821 .
  • the eigenvalue time duration of signal track 821 are compared to a decision threshold 809 to map and recover the encoded message “101011” of signal track 821 .
  • M ⁇ 1 decision thresholds are required.
  • FIG. 9 An illustration of a portion of the block-to-block eigenvalue correlation result is shown in FIG. 9.
  • the large complex plane diagram in the top portion of the FIG. 9 shows the complex eigenvalue locations of the SFOCMP (GEVD) results and the predicted locations of the blockwise eigenvalue correlator.
  • the legend identifies the four levels of eigenvalue correlation confidence, (“New,” “Tentative,” “Candidate,” and “Confirmed”).
  • the five consistent signal eigenvalues of the four co channel interferers and of the covert emitter are indicated by the smaller rectangle.
  • the covert emitter eigenvalue lies below the imaginary axis.
  • Blocks 30 and 31 correspond to two different symbols.
  • the lower portion of FIG. 9 illustrates the blockwise changes in eigenvalue locations over blocks 30 to 34 where the message symbol “boundary” occurs between blocks 30 and 31 .
  • the covert emitter eigenvalue lies in a consistent location above the imaginary axis of the complex plane.
  • the stronger co channel interference source eigenvalues, within the oval for block 34 lie in consistent locations over all five blocks.
  • the consistency of the spatial locations of the interference emitters and covert emitter allow message recovery embodiment to “stitch” together eigenvalues that fall within a “time-gate” and come from the same spatial location.
  • FIG. 10 An embodiment of the receiver 303 is shown in FIG. 10.
  • the receiver 303 augments the message recovery process by using spatial information regarding the locations of the transmitters in the field of view.
  • signal detections are performed in processor 1001 which is composed of a BSS processor 309 and (eigenvalue) tracking processor 313 previously described.
  • the results of the blockwise processing and eigenvalue correlation previously described are stored in the Signal History Database 317 .
  • the data envisioned for the database are signal characteristics measured which include but are not limited to eigenvalue identification number, time initialized, time deleted, number of observations, an estimate of the eigenvalue, and importantly an estimate of the spatial location associated with the eigenvalue.
  • the detections surviving the time gate test are stored in the signal history database 317 along with other associated eigenvalue characteristics, including but not limited to Signal track identification, Symbol estimate, Time initialized, time deleted, number of observations and importantly spatial location.
  • the detections are assembled into strings of symbols (tracks) using the criteria that a message must associate with a “consistent” spatial location.
  • the preferred spatial information is angle-of-arrival (AoA) or geolocation for each signal detected in each block which are determined in the AoA/Geo computation processor 315 .
  • spatial location refers to AoA if only calibration is known, or geolocation if both receiver position and array calibration are known. It is preferable that the spatial relationship between the transmitter 302 and receive array 327 be fixed over the message length, however, a slowly varying spatial property can be accommodated by introducing a method to correlate a sequence of spatial variables across contiguous blocks of data in the receiver 303 specifically in the Message Sequence Correlator 323 of the Covert Communication Processor 319 .
  • Spatial correlation can be broaden to include simply steering vectors. This is useful when the array and transmitter have a stable geometry. Relative motion between the transmitter and sensing array causes the steering vectors to have a detrimental time dependency. Again, if the spatial variable for correlating the message data is “slowly” varying then small incremental changes can be tolerated.
  • the covert messages is indented to be recovered using a “consistency” of the spatial domain information of the computed eigenstructure from the SFOCMP for the signals of interest. But to account for the possibility of “fixed” location emitters and other emitter who are not of interest to the communication process, a time-gate decision process as noted earlier is advantageously applied. This way the receiver need only attempt to decode “message strings” that emanate from “consistent” spatial locations with the appropriate time character.
  • the steering vector will be estimated using the data available from the GEVD of the SFOCMP.
  • the estimation of the other spatial variables, AoA and/or geolocation can be determine by methods well-known in the art.
  • the characteristic that allows this computation is that the eigenvalues and eigenvectors of the GEVD of the SFOCMP have a 1:1 correspondence as in standard eigenanalysis. So, when the eigenvalues are used to measure the high-order statistical properties of the received signals, an index relating directly to where that signal with that characteristic emanated from is available. Again, the spatial dimension for signal association can be exploited to remove any ambiguity of the temporal decomposition, since we assume that no two emitters are identically located.
  • Signal tracks 1021 , 1023 and 1025 containing associating matrix pencil eigenvalues from consistent sources as evident from common spatial variables are used to recover the message by correlating track information against spatial variables and time gates. This recovery process in performed by the covert communication receive processor 319 .
  • the recovered messages are composed of sequentially ordered matrix pencil eigenvalues with a duration within the time gate originating from the same location as determined by the steering vector estimate, AoA or geolocation.
  • the “signal sequence 1” 1021 includes eigenvalues 1, 4, 7, 10, and 12 with durations mapping message “01011”, these eigenvalues and their associated signals all originated from a “consistent” spatial location, namely, Geo 1 , AoA 1 or had the same steering vector SV 1 .
  • “signal sequence 2” 1023 includes eigenvalues and durations with common steering vectors, AoA, or geolocation and the durations pass the time gate filter thus recovering a “000100” message.
  • “Signal sequence 3” 1025 however had two eigenvalues 3 and 9 with high spatial correlation, however the signal is not decoded since the durations exceed the time gate allowable. In short, the receive processor recognized that the time durations in sequence 1025 do not correspond to any symbols of the expected alphabet.
  • the spatial variables aid in rejecting extraneous pulsed emitters based on their spatial locations being anti-correlated over time to the persistent spatial locations of the covert emitter(s).
  • spatial variables allow the basic signaling approach to support multiple access of covert emitters without undue burden in the receiver for properly assembling the pulsed message sequences. This is because the additional covert emitters will very likely emanate from resolvable spatial locations, and the receiver can use the consistency of the spatial locations over time to associate the proper message sequence.
  • the message sequence is represented by the time durations of the eigenvalues of the appropriately selected matrix pencil, where we have preferred the SFOCMP approach.
  • any emitter is independent of the exact value of its corresponding eigenvalues available from the GEVD of the SFOCMP.
  • the spatial variables provide additional “distance” in the recovery process, since it is now multi-dimensional. For example, two signals may have very similar eigenvalues. But, if their spatial locations are resolvable by the receiver, and fairly constant, then the eigenvalues corresponding to those spatial locations can be easily assigned. Then the message can be recovered using the time duration of each eigenvalue in the sequence assigned a given spatial location using the same technique as previously described when only a single covert signal was in view.
  • FIG. 11 is a flow chart for covert communication with a transmitter 1102 and a receiver 1103 using message recovery with spatial information as shown in FIG. 10.
  • the source 1100 provides a data stream to the covert transmitter 1102 as a stream of M-ary alphabet symbols.
  • each M-ary alphabet symbol is assigned a unique duration.
  • Plural waveforms with determinable high order statistic that are constant and non zero are generated in block 1106 .
  • One of the plural waveforms is selected in block 1108 and is transmitted for the symbol assigned duration in block 1110 .
  • the receiver 1103 receives the transmitted waveform along with environmental and random noise in a multi element array as shown in block 1105 .
  • matrix pencil eigenvalues MPE
  • MPE matrix pencil eigenvalues
  • a spatial variable or location is determined from the MPE in block 1111 , again the spatial variable can be selected from a Steering vector, AoA or geolocation.
  • the MPE detection durations are time gate filtered in block 1113 to reject detections outside of the duration minimum and maximum thresholds.
  • the passed MPE and durations are sequentially sorted by common or consistent spatial variable in to signal tracks in block 1115 .
  • the signal tracks are maps are mapped (i.e. their durations are correlated to symbols and their respective assigned durations)to recover the message in block 1117 .
  • the message may be processed otherwise transformed to provide the message at the sink 1101 .
  • FIG. 12 presents a receiver 303 for processing covert messages without spatial location or spatial variable information.
  • Processor 1201 contains BSS processing 309 and eigenvalue correlation 313 operating in the same manner as previously described, however without transferring or developing spatial variables. Therefore the eigenvalue characteristic and time stamps for determining eigenvalue duration are stored in the signal history database 317 , without a spatial variable.
  • the tracks are not sorted by spatial location by rather eigenvalues.
  • Signal tracks contain only detections of the same eigenvalue with a detected duration that passes the time gate.
  • Signal 1221 includes only one eigenvalue with 4 detection durations that pass the time gate. Signal 1221 thus is mapped to recover message “0010”.
  • the durations of the “OFF” period of pulse is irrelevant as long as its duration is greater that 1 block to enable the previous detection to be terminated the tracking correlator 313 .
  • Signal tracks 1223 and 1225 are both ignored or discarded for failing the time gate.
  • FIG. 13 is a flow chart for covert communication with a transmitter 1302 and a receiver 1303 using message recovery without spatial information as shown in FIG. 12. The process is similar to that shown in FIG. 11.
  • the source 1300 provides a data stream to the covert transmitter 1302 as a M-ary alphabet symbols.
  • each M-ary alphabet symbol is assigned a unique duration.
  • a signal waveform with a determinable high order statistic that is constant and nonzero is generated in block 1306 .
  • the waveform is transmitted for the symbol assigned duration in block 1310 .
  • the receiver 1303 receives the transmitted waveform along with environmental and random noise in a multi element array as shown in block 1305 .
  • matrix pencil eigenvalues are obtained in block 1307 and the MPE are tracked to determine their respective durations in block 1309 .
  • the MPE detection durations are time gate filtered in block 1313 to reject detections outside of the duration minimum and maximum thresholds.
  • the passed MPE and durations are sequentially sorted by specific eigenvalue into signal tracks in block 1321 .
  • the eigenvalue durations forming the “signal sequences” are mapped (i.e., their durations are correlated to symbols and their respective assigned durations) to recover the message in block 1317 and delivered to the sink 1301 .
  • FIG. 14 is a flow diagram for covert communication using encoded eigenvalues and spatial information.
  • the source 1400 provides a data stream to the covert transmitter 1402 as a M-ary alphabet symbols.
  • the M-ary symbols of the messages are encoded as MPE with unique durations.
  • the eigenvalue is transmitted for the symbol assigned duration as a respective source waveform in block 1410 .
  • the source waveform is the low order signal whose high order statistic creates the respective eigenvalue.
  • the receiver 1403 receives the transmitted waveform along with environmental and random noise in a multi element array as shown in block 1405 .
  • the matrix pencil eigenvalues encoding the symbol are recovered and tracked in block 1407 to determine their respective durations and spatial information in block 1409 .
  • the MPEs are decoded based on MPED and Spatial information as previously described in block 1417 to recover the message and provide the message symbols to the sink in block 1401 .
  • the noise is white Gaussian noise, color noise or interferer signals.
  • FIG. 15 is a representation of a binary message sequence 101011 transmitted via data carrying waveforms according to an embodiment.
  • the symbol are transmitted by alternating waveforms of BPSK, QPSK and GMSK.
  • the carrier waveforms be modulated with random data at a rate greatly exceeding the covert message symbol rate.
  • Other waveforms including to DQPSK, DBPSK, FSK, QAM, DPCM are equally applicable and are also envisioned, however embodiments of the invention should not be construed to be limited to the particular waveforms listed.
  • the waveform duration discriminates the message symbol.
  • the symbol 1 has a duration of 20 blocks and the symbol “0” has a duration of 10 blocks.
  • the message is transmitted as a GMSK waveform for a duration of 20 blocks indicating a symbol “1”.
  • the next symbol is transmitted as a different waveform QPSK for a duration of 10 block indicating a symbol “0”.
  • the next symbol is transmitted by the GMSK waveform for a duration of 20 blocks again indicating a “1” symbol and a QPSK waveform is used to transmit the “0” symbol.
  • the symbol “1” is then transmitted as a BPSK waveform for a 20 block duration and the next symbol “1” is transmitted by the QPSK waveform for a duration of 20 block.
  • the symbol “1” is then transmitted as a BPSK waveform for a 20 block duration and the next symbol “1” is transmitted by the QPSK waveform for a duration of 20 block.
  • the QPSK waveform for a duration of 20 block.
  • identical waveforms can communicate different symbols while maintaining waveform content and message independence.

Abstract

A transmitted signal exploits unique higher order statistics of temporally dependent waveforms to encode message symbols with unique durations enabling low power transmission of the covert message. The message is independent of the characteristics or encoded data of the transmitted waveform. The method uses spatial fourth order cumulants or spatial second order moments in a Blind Source Separation and generalized eigenvalue decomposition to determine unique matrix pencil eigenvalues. Sequential detection in successive blocks determine the duration of the eigenvalue. The durations of the detected eigenvalues are sorted sequentially and can be sorted spatially by Steering vectors, AoA or geolocation into signal tracks that are mapped to recover the transmitted cover message. Generation of the transmitted signals with unique high order statistics may be accomplished using combination of noise generators, temporal filters, signal sources, combiners and switches. The receiver includes a multi-element array and does not need a priori knowledge of the transmitted signal source to recover the message. The methods and apparatus for covert communication does not require typical demodulation. The covert communication system is multiple access.

Description

    RELATED APPLICATIONS
  • The present application is related to and co-pending with commonly-assigned U.S. patent application Ser. No. 10/360,631 entitled “Blind Source Separation Utilizing A Spatial Fourth Order Cumulant Matrix Pencil”, filed on Feb. 10, 2003, the disclosure of which is hereby incorporated herein by reference. [0001]
  • The present application is related to and co-pending with U.S. Provisional Patent Application Serial No. 60/374,149 filed Apr. 22, 2002 entitled “Blind Source Separation Using A Spatial Fourth Order Cumulant Matrix Pencil”, the entirety of which is hereby incorporated herein by reference. [0002]
  • The present application is related to and filed concurrently with U.S. Provisional Patent Application Serial No. ______ entitled “Cooperative SIGINT for Covert Communication and Location Provisional”, the entirety of which is hereby incorporated herein by reference.[0003]
  • BACKGROUND
  • In the advent of globalization, information is a fundamental and valuable commodity. Information and intelligence regarding national defense and security comes at an even a higher premium. In may instances this information and intelligence can only be obtained covertly so as to not reveal the source. Information regarding location of a source such as for surveillance or combat search and rescue can be degraded in value if detected by unfriendly entities, such as enemy forces in the case of a downed pilot or a marked terrorist under surveillance. Additionally, command and control among conventional and/or special operation forces can betray missions by alerting hostile governments or organization of their presence. Embedded voice and data bugs can also expose the existence of these forces if their transmissions are detected resulting in the desired information moving out of reach of the bugging device. These are but a view of the many possible scenarios where a message is desirably sent covertly from a mobile or fixed transmitter to a remote receiver or, more importantly, a friendly or cooperative receiver with out the presence of a signal being detected by an unintended receiver. Covertly transmitting a signal may be necessary to reflect the fact that the transmission is conducted in such a way as to avoid intentional or inadvertent detection by systems monitoring the electromagnetic (EM) environment local to the covert transmitter. However, any transmission, covert or not, may be desired to be received only by intended receivers. [0004]
  • Intentional detection of the signal or message can be accomplished in military systems that use specially designed electronic support measures (ESM) receivers. These ESM receivers are often found in signal intelligence (SIGINT) applications. In commercial applications, devices employed by service providers (i.e. spectral monitors, error rate testers) can be used to detect intrusion on their spectral allocation. Inadvertent detection can also occur, such as when a user or service provider notices degradation in link performance (e.g., video quality, audio quality, or increased bit error rate). [0005]
  • The term covert also implies the additional goals of evading interception and exploitation by unintended receivers. Interception is the measurement of waveform features or parameters useful for classifying/identifying a transmitter and/or the waveform type and/or deriving information useful for denying (i.e. jamming) the communication. Exploitation is processing a signal by an unintended receiver in the attempt to locate the transmitter and/or recover the message content. In the broad literature on covert communications these characteristics as applied to transmitted information signals are referred to as low probability of detection (LPD), low probability of intercept (LPI), and/or low probability of exploitation (LPE) by an unintended receiver. [0006]
  • Typically an LPD communication system shown in FIG. 1 is designed by minimizing the radiated energy in all directions. The goal is to have the [0007] emitter 102 radiate zero energy in all directions, save for the direction of the intended receiver 103. This approach limits the geographic locations in which the signal would likely be detected. The area 110 receives little radiation and thus the unintended receiver 105 would have greater difficulty detecting the signal than the unintended receiver 106 in area 111. Typically, LPD systems use only enough energy in the direction of the intended receiver to be reliably detected and the message recovered by the intended receiver. The radiation in non-preferred directions is primarily controlled by the aperture (i.e. antenna system) and additional performance improvement for LPD can be gained by employing a directional antenna. However, there remains a need to minimize the radiated energy necessary for the intended receiver to detect and recover the message and thus reduce the ability of an unintended receiver in area 111 to detect the signal.
  • Given the desirability to transmit messages covertly, it is helpful to understand considerations that enhance or degrade LPD, LPI and LPE. An unintended receiver such as the receiver [0008] 103 in FIG. 1 with the goal of detecting a covert communication must reliably differentiate between the binary noise-only and signal-plus-noise hypothesis. As is known to those of skill in the art, for an unintended receiver the signal detection process is typically based on an energy threshold. The energy the receiver measures is given by Etot=PavgTxmit. Where under general conditions the power Pavg is the received covert signal power S plus internal receiver noise power N. Hence, Etot=(S+N)Txmit. If the signal power used to communicate is only a small fraction of the receiver noise, S<<N, it is extremely difficult for the unintended receiver to reliably detect the presence of the covert signal because the total energy detected will only be marginally greater than the noise-only (S=0) case.
  • In prior art LPD systems, the bandwidth of the message signal is artificially broadened to reduce the energy loading on the band, as is done in direct sequence spread spectrum (DSSS) communication, the intended receiver has a processing gain over the unintended receivers that use the same co-channel receiver bandwidth to capture the signal. However, for prior art DSSS, the intended receiver [0009] 103 must have knowledge of the spreading code and be able to acquire it and remain synchronized to the covert transmitter code, in addition to the usual communication receiver functions of carrier and timing recovery. This requirement, for some applications, is not a problem but for other applications, covert message transmission functionality is greatly diminished.
  • Minimizing transmit power has two direct system benefits. First, the total signal power used will be a small fraction of the total noise power present in the same band. Thus, if the message is limited in time duration, the total energy measured by an unintended receiver [0010] 106, which may be an ESM receiver, is indistinguishable from a noise-only environment. Since ESM receivers are often of energy threshold type, there is an obvious trade-off of average power for time duration in order for a signal to remain undetectable. Second, the low transmit power scenario enables usage by transmitters with very limited power supplies (i.e. batteries).
  • Therefore, as naturally arise in military environments such as depicted in FIG. 2, there is a need for a low power message system and method, covert or otherwise, such as covert communications for Intel or Special Forces, “stealth” IFF for low observable ground vehicles, and combat search and rescue (CSAR). There is also such a need in a number of civilian or public safety applications as well, such as asset tracking/location or “lost child” detection/location and surveillance. In particular in these latter-described applications it may be particularly desirable to receive both a message and location the source of the message. [0011]
  • Embodiments of the present inventive system and method address the above needs while requiring only an extremely low power signal.[0012]
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • FIG. 1 is a general representation of a prior art approach to LPD. [0013]
  • FIG. 2 depiction of covert scenarios. [0014]
  • FIG. 3 is a representation of an embodiment of a waveform independent covert communication system. [0015]
  • FIG. 4 is a depiction of a binary symbol message 101011 according to an embodiment of the invention. [0016]
  • FIG. 5[0017] a is a representation of a spectrogram of a transmission of the binary symbol message depicted in FIG. 4 and revealing the symbol boundaries of the message according to an embodiment of the invention.
  • FIG. 5[0018] b is a representation of a spectrogram of a transmission of the binary symbol message depicted in FIG. 4 as might been seem by a typical ESM/SIGINT spectral monitor.
  • FIG. 6[0019] a is a generalized schematic representation of embodiments of the invention.
  • FIG. 6[0020] b is a schematic representation of an embodiment of the invention with dual signal sources.
  • FIG. 6[0021] c is a schematic representation of an embodiment of the invention with dual temporal filters.
  • FIG. 6[0022] d is a schematic representation of an embodiment of the invention with dual non-Gaussian noise generators.
  • FIG. 6[0023] e is a schematic representation of an embodiment of the invention with dual non-Gaussian noise generator and temporal filter branches.
  • FIG. 6[0024] f is a schematic representation of an embodiment of the invention with dual temporal filters and dual signal sources.
  • FIG. 7 is a schematic representation of a Laplacian noise generator with multiple signal sources according to an embodiment of the invention. [0025]
  • FIG. 8 is a representation of eigenvalue tracks for a frequency shifter laplacian noise waveform for message “101011” according to an embodiment of the invention. [0026]
  • FIG. 9 is a representation of block-to-block eigenvalue correlations. [0027]
  • FIG. 10 is a schematic representation of a message recovery system with spatial information according to an embodiment of the invention. [0028]
  • FIG. 11 is a flow diagram for covert communication via message recovery with spatial information according to an embodiment of the invention. [0029]
  • FIG. 12 is a schematic representation of a message recovery system without spatial information according to an embodiment of the invention. [0030]
  • FIG. 13 is a flow diagram for covert communication via message recovery without spatial information according to an embodiment of the invention. [0031]
  • FIG. 14 is a flow diagram for covert communication via encoded eigenvalues according to an embodiment of the invention. [0032]
  • FIG. 15 is a representation of a binary message sequence 101011 encoded using different carrier waveforms.[0033]
  • DETAILED DESCRIPTION
  • A useful feature of embodiments described herein is the use of time duration to convey information, while the individual waveforms used to convey the message (in the time duration of the fourth-order characteristic) are in a sense superfluous or independent of the message to be conveyed. This is a significant advantage for ubiquitous application, allowing for parasitic use of present communication infrastructure and devices. Thus there are few restrictions on the pairing between potential covert transmitters and the intended receiver using the disclosed covert communication methods and apparatus because of the independence of the information transfer on the “carrier waveform”. This is unlike prior art systems where the receivers designed or instantiated for a certain signal type cannot accurately recover the message if the receiver is presented with another signal type. However, embodiments of the present invention by contrast can function equally well for any waveform, and the receiver does not require any a priori knowledge of the “carrier waveform”. In fact, embodiments of the covert transmitter can be waveform agile without informing the intended receiver. [0034]
  • The embodiments herein are predicated on selecting and transmitting carrier waveforms with unique higher order spatial statistics. Such higher order statistics include 2[0035] nd order spatial moments and 4th order spatial cumulants. There are a number of ways that the fourth-order statistics of high order statistic of a signal may be altered. The number of type of encoding schemes mapping an M-ary alphabet to time durations is limitless, but in point applications there will be preferred methods since they ease hardware and software realization of the system. Specific embodiments are identified herein. One method is by design, such as using frequency shifted Laplacian noise. Another method is to use the inherent properties of conventional waveforms (e.g., BPSK, QPSK, GMSK, QAM, DBPSK MFSK, FSK, and DQPSK). Yet another is to alter channel filtering in the transmitter. Common to all these approaches is that a high order statistic, specifically, but not limited to, the fourth-order characteristic, of the transmitted signal is recoverable with a spatial fourth order cumulant matrix pencil (SFOCMP). Furthermore, the message symbols map to time durations of these higher order statistics and not to specific values of the SFOCMP. Hence, from the receiver's point of view, the communication is independent of waveform selection by the transmitter. The primary restriction is that the receiver and transmitter must use the same “codebook” of time durations and alphabet. However, the “codebook” need not be known a priori by the receiver in all instances. As shown herein, the receiver will be able to discern the different time durations and be able to extract some useful information therefrom.
  • To recover the message information, the waveform is not conventionally demodulated. Rather, a straightforward block or batch estimation algorithm estimates the generalized eigenvalues of the SFOCMP for each signal in the receiver field-of-view (FOV) over time using only the array output. For this discussion we assume the data has been digitized appropriately. The sizing for the block processing (e.g., the block of contiguous array observations, sometimes known as “snapshots”) is dependent on several factors. Chiefly we must ensure that each block has enough sample support so that the eigenvalue estimates from the GEVD of the SFOCMP in each block over a symbol duration nominally match. This means that estimation error is negligible. Accordingly, changes in the eigenstructure can be reliably detected, and this change indicates a symbol boundary. The degree to which a nominal match is required within a block depends on the complexity of the signal environment (e.g., extraneous co-channel signals), the communication errors (e.g., partially received messages) tolerable in a given application, and the receiver processing resources to recover the message in a timely manner. [0036]
  • Further, in practical situations, as power sources become impaired (e.g., batteries running low on power), transmitted waveforms become increasingly distorted. This situation limits the effectiveness of matched filtering as used in prior art systems, since the concept of matched filter relies on knowledge of the transmitted waveform in the receiver. Embodiments of the inventive technique are impervious to such distortion since it is the duration and not the actual value of the eigenvalues of the SFOCMP that matter. So as long as the eigenstructure characteristic of the distorted signals is nominally constant during a message symbol, the inventive system and method is robust as to degraded transmitter performance. Therefore, the present inventive system and method will operate successfully under conditions that would normally be detrimental to conventional systems. The use of lower order matrix pencils are also contemplated by the present inventive system and method. [0037]
  • In FIG. 3, the [0038] source 300 produces a digital data stream of a covert message to be transmitted by the covert transmitter 302. The digital data can either be inherently digital (e.g., Morse coded messages) or analog (e.g., voice) and subsequently digitized. The digital data can also include framing or formatting of the message. Typical of the framing would be start/stop and data fields. Though other fields can be used as needs dictate. This framing structure can ensure that the receiver can reliably find the beginning of a message for synchronization. Source coding or compression could also be applied to the incoming data stream to reduce the required bandwidth. The user may also encrypt the data to protect it by an optional encryption device of suitable complexity. This data is then output to the forward error correction (FEC) module, which currently is envisioned as applying block coding. The coding is useful to aid the receiver in resolving message ambiguities say caused by fades or other unresolved time coincident measurements, which in this system could be processed as “erasures” up to the correction capability of the code. Additionally, the error correction code would find application when the symbol stream is improperly decoded due to statistical variations in the time duration measurement process. Thus, the potential for a message protocol using automated repetition of the message might be advantageous as error patterns in each transmission will likely be different. This type of protocol would require the framing as stated above, and possibly the inclusion of other fields such as the number of times the message has been or will be repeated. Message formatting, compression, encryption and/or FEC is performed by the processor 304. The processing in processor 304 is not necessary for practicing the present invention but is included to indicate the ability to add additional robustness and reliability to the transmission of a message.
  • The binary symbol stream is from [0039] source 300 is then mapped into suitably selected M-ary symbols (or time durations) in the symbol mapping device 306, via look-up tables or other codebook like mapping mechanisms. The alphabet size is limited only by system implementation and design considerations for particular applications. The encoded data is then applied to the “carrier waveform” by carrier waveform generator 308 which is a signal generator, implementations of which are further described below. The “carrier waveform” has its fourth-order characteristic modified according to the control by the M-ary alphabet. There is no limit to the strategies potentially adopted by the covert transmitter for this operation, so long as the characteristic is measurable reliably by the receiver and it conforms to the time duration and alphabet size assumed by the receiver. To this end, there will be some minimum duration and maximum duration for a symbol, and a preferred duration increment for each symbol in the alphabet. The exact choices depend on application, however, making the durations too disparate can negatively affect data rate limiting this technique to lower data rate applications. It is desirable to provide durations that are easy to resolve into the M-ary symbols. The baseband data from the carrier waveform generator 308 is fed to a suitable RF subsystem 310 and antenna 312 for transmission.
  • The minimum duration and duration increment must be such that synchronizing the data block boundaries used in the receiver to that of the symbol timing in the [0040] transmitter 302 is not relied upon. It is desirable that the covert transmitter 302 use a fundamental signaling period of several (e.g., 10) “receiver block” durations for the minimum signal, and may have a signaling duration increment of the same size to define the alphabet. However other choices are applicable depending on the particular implementation and application. The goal is to make the time duration alphabet as disparate as possible while meeting performance objectives (e.g., data throughput). Sample data rate computations can be determined as shown below.
  • Defining S as the array snapshots/block, “b” blocks for the minimum length symbol, “B” blocks for the maximum length symbol, and “R” as ADC (analog to digital converter) conversion rate in the receiver, the minimum and maximum symbol durations for a binary alphabet are:[0041]
  • T sym(min)=S*b*1/R=5,000*10*10−9=50×10−6 sec
  • T sym(max)=S*B*1/R=5,000*20*10−9=100×10−6 sec
  • The values b=10 and B=20 along with S=5,000 and R=1 Gsamples/sec are subject to implementation choice, and used here for illustration only. Assuming that a system would have an equal number of binary symbols of each type, the average (over the long-term) data rate is nominally 13 kbps. If M-ary signaling is implemented with the same maximum and minimum symbol durations, the data rate can be improved by factor of log[0042] 2(M), but at potentially increased channel BER. Achievable data rates are in principle limited by operating conditions (received SNR, tolerable BER, cumulant estimation errors, etc.).In addition to reliably detecting a change in the SFOCMP eigenvalues using a basic correlation technique, a minimum b consecutive blocks are required (currently b=4) for each of S vector samples from the receive array, thereby making the theoretical minimum symbol duration equal to bS(1/R). Similarly, the incremental time duration for the alphabet should be at least ES(1/R), where E is the number of blocks desired by the designer to provide a balance between adequate safety margin in the time duration decision process and required throughput rate. In theory, E can be as low as unity which would enhance the achievable signaling rate for a fixed alphabet size. However, this is likely not a practical choice since numerous errors can occur due to the receive block processing not being time aligned to transmitter symbol boundaries. Conversely, selecting a large value for E for a fixed alphabet minimizes the effect of block misalignment, but limits the available data throughput. In the present binary example E=10 thereby yielding symbol durations of 10 blocks and 20 blocks for the binary alphabet. For M-ary alphabets, one might construct an alphabet with uniformly spaced symbol durations yielding durations of (b+qE)S(1/R) where q=0, 1, 2, 3, . . . , M−1. An alphabet of equally-spaced symbol durations might be selected for convenience of implementation, but non-uniform spacing of the symbol durations is not precluded by the teachings of the present invention. It shall be understood by those of skill in the art that the values given above are exemplary only and are not to be construed as limiting the invention in any way.
  • The [0043] receiver 303 uses an N-element (or port) receive array 327 and an RF processor 305 to obtain the transmitted signal. In order to capture the temporal character (i.e. the time duration modulation of the SFOCMP eigenvalues) of the covert signal, the array data is first sampled and digitized at some rate suitable for the application. Each array output is digitized simultaneously producing a vector observation in the vector digitizer and buffer 307. The array output data is buffered and subdivided into non-overlapping blocks in 307. Block-wise across signal samples (i.e. the vector observations) are then collected from an array at the intended receiver aperture and the cumulants are block estimated, the matrix pencil is formed, and the generalized eigenvalue decomposition (GEVD) is performed by the Blind Source Separation processor 309. The operation of the BSS requires the selection of a triplicate of time lags provided by the time lags selection device 311. The GEVD provides a set of N eigenvalues λk (b) and N eigenvectors Vk (b), where k=1, 2, 3, . . . , N (i.e. assuming an N-port array is used) for each block of data. The superscript b is used as a block counter in the receiver. We assume there are Ms generalized eigenvalues representing the SFOCMP properties for each of the Ms signals in the field of view (FOV) of the receive array 327, M≦N. The remaining N−Ms eigenvalues are of the indeterminate (i.e. 0/0 type). Thus when using a sequence of block estimates for the SFOCMP eigenvalues of the Ms consistent signals will be apparent. Further processing is required and performed in the communication receive processor 319 to assemble valid messages. A major part of this processing is to use spatial information available from the GEVD processing. The spatial information regarding the transmitter location and the message content are linked in a 1:1 fashion by the generalized eigenvalues produced by the processing in 309. Using the “side information” of the available spatial variable greatly eases message recovery since we assume that a transmitter spatial location will be “slowly” varying (i.e., changing at a rate much less than the symbol rate of the message), hence a message can be reconstructed in part by looking for message symbols represented by eigenvalues and their durations, associated with a “consistent” location. The designer must ensure that the symbol duration alphabet has sufficient minimum support and increments such that the practical time duration recovery issues where ambiguous results can be obtained do not adversely affect the system performance.
  • As may be apparent to those of skill in the art, there may be some advantage to overlapping blocks of the data. However, the following discussion deals with non-overlapping blocks. On each block, the two fourth-order spatial cumulant matrices required to form the SFOCMP are formed using pre-selected delay triplets. The delays can be either pre-selected, or subjected to online modification using a programmed search routine (if necessary). This search routine might be necessary when certain conditions, such as repeated eigenvalues for different signals are encountered. Processing of the received covert message this difficulty is likely of no consequence. However, provisions are made for signals whose eigenstructure match at the delays selected to be repressed at different delays to provide improved discrimination if desirable. After the matrix pencil is formed, the GEVD is computed. From the GEVD, the eigenvalues and eigenvectors are used to determine the signal environment over time block b. Subsequently, the eigenvectors are used to determine the signal steering vectors and then the eigenstructure is correlated block wise in the [0044] Blockwise Eigenvalue Correlator 313 to determine any changes in the signal environment. A change, such as symbol boundary, in the communication signal will alter its contribution to the signal environment eigenstructure, measured by the SFOCMP, in a detectable manner. This means a “significant” movement in the complex plane of eigenvalues. As signal changes are detected, those signals are cued for storage in the signal history database 317. As part of the storage procedure, the spatial location of the signal is determined (i.e., either angle-of-arrival (AOA) or geolocation, whichever applies given the specific application) by the AOA/Geolocation processor 315. Additionally, the steering vector can be recorded, which is useful when refined spatial information is unavailable and the relative motion of the transmitter and receiver is negligible. The eigenvalues no longer correlating with the present signal structure are also written to the database. The temporal support (i.e. durations) of the eigenvalues no longer correlating with the current signal structure is measured and stored. All this data is formed and recorded in the signal history database 317 along with other ancillary data that may be useful for signal post-processing applications such as data mining or covert message recovery.
  • An important function of a tracker is the track initiation and deletion logic. An embodiment of the tracks uses a fixed distance and a fixed number of consecutive “good associations” for initiation and a single “no association” for a track deletion. A “good association” is any measurement that is “close enough” to track. A “no association” condition occurs when all the measurements are “Too far” from a particular track. the distance indicative of an good association can be set empirically or experimentally. Track initiation and track deletion strategies can also be used to adapt to various situation. A Kalman-like approach to association gates can be adapted as the number of observation for a track are accumulated. Such an approach also has the advantage of replacing fixed averaging of the measurements. [0045]
  • To recover the covert message, in the covert communication receive [0046] processor 319, first the time durations are time-gated and those passing the time-gate are then mapped to the M-ary symbols of the “codebook” by the Symbol Demapper 321. The Message Sequence Correlator 323 is an association engine that “stitches” all the signal eigenvalues durations and symbol values corresponding to a “consistent” spatial location into a time-ordered stream or signal track. After the data is assembled into streams, the data is supplied to the optional data processor 325, where the data can then be FEC decoded, decrypted (if necessary), decompressed (if necessary) and then the message can be checked for synchronization, as well as, the data provided to the sink 301 with the appropriate format.
  • The design also allows for multiple access for communications. Consider the case where multiple remote covert emitters are sending data. It is unlikely that they would have exactly the same fourth-order cumulant representation, even if they are using the same base waveform. This is because any deviation from nominal waveform implementation (e.g., frequency change, waveform change, matrix pencil eigenvalue change, phase noise, I/Q imbalance, timing jitter, phase jitter, symbol rate change, pulse shape change, a fourth-order statistic change, relative rotational alignment of a signal constellation change, power amplifier rise/fall time change, and Doppler shift change) causes the fourth-order statistics of these signals to differ. Further, the multiple access signals are assumed distinguishable by spatial location. Of course this requires enough data to be collected to resolve the location, and the array must also provide such resolving power. But, if automated location is not possible at the receiver, due to, for instance, no calibration, the covert transmitters may still have multiple access if the multiple access signals can be assured uniqueness amongst themselves and the environment of sufficient degree in the SFOCMP eigenvalues. The receiver need not know the exact eigenvalues that will be used, but in this mode it is incumbent on the individual transmitters to use one and only one eigenvalue and not switch waveforms. In principle, correlation algorithms to properly sort this data are readily imaginable, though the details depend specifically on the signal designs. And of course if location of the covert emitter is still is necessary, the remote covert emitter could encode his position, say from his GPS, into the data stream. This might make the encryption more important, and certainly elevates the need for LPD/LPI capability. [0047]
  • Returning to the concept of the SFOCMP, it is not correct to think of the possible eigenvalues for a given transmission to constitute a signaling constellation. While signaling using the eigenvalues could be approached this way, it is not a requirement. This would restrict the choices and implementation and performance for the system. The full flexibility of the proposed approach is reached when we allow a non-constant mapping of the baseband symbols (i.e. time durations) to the potential plurality “carrier waveform” options available in the transmitter. The number of “carrier waveforms” can exceed the M-ary alphabet. For example, some alternative “carrier waveforms” might be encountered as distortions of the canonical set such as when power systems or components degrade over time or for some reason the system designer desires a larger signal set be available. [0048]
  • A specific example of waveform design to be feature-less is shown in FIG. 4 and FIGS. 5[0049] a and 5 b. FIG. 4 is a representation of a binary message sequence 101011 where the symbols are discriminated by duration. In the transmitter the binary sequence is mapped to binary time durations. It shall be readily appreciated by those of skill in the art that the process may be extended to M-ary time durations. In response to a new symbol, the carrier frequency for the embodiment shown is offset at the time boundaries, thus shifting the fourth-order cumulant (i.e., the eigenvalue) for that symbol. Ultimately the time duration of the eigenvalue is recovered to reconstruct the message. The frequency shift for the repeated “1” symbol at the end of is needed (in this case) to differentiate a repeated “1” from a single “0” and visa versa in the case of a double “0” (not shown). In particular the temporal correlation of the SFOCMP eigenstructure must be broken at a message symbol boundary. This break of temporal correlation in FIG. 5a is accomplished by altering the fourth-order cumulant with a frequency shift.
  • FIG. 5[0050] a is a spectrogram of the emitter message shown in FIG. 4 when using the frequency shifts of the preferred embodiment discussed herein. The symbol boundaries defined by frequency offsets at the transmitter are clearly evident in the spectrogram. FIG. 5b is a representation of a typical Electronic Support Measures (ESM) receiver/Signal Intelligence Spectral Monitor output viewing the received signal of FIG. 5a. With reference to the spectrogram shown in FIG. 5b, the time varying spectral pattern of the message shown in 5 a is buried in interference and noise (SNIR<−10 db) making the covert message very difficult to detect, and thus even more difficult to intercept and exploit. As is know to those of skill in the art the acronym SNIR stands for Signal to Noise and Interference Ratio.
  • Given an environment with several interferers and the already negative received SNR an unintended receiver (even using a front-end filter) will likely not reliably detect the presence of the covert signal. But even if a machine detects the presence of the signal energy, it would likely not be acted upon since it would fail all modulation recognition tests and show no exploitable temporal structure. The signal represented in FIG. 4[0051] b is frequency shifted Laplacian (double-exponential) Noise. Viewed by a casual observer the signal would mimic additional thermal noise, hence even if detected by an energy detector of suitable design, the detection would likely be discarded.
  • A mathematical element of the invention is the use of spatial high order statistics to separate signal sources, such as a blind source separation algorithm that utilizes a normalized spatial fourth-order cumulant matrix pencil and its generalized eigenvalue decomposition (GEVD). The equations presented herein use the following subscripting convention. Quantities relating to the array observations available to the system are denoted with a boldface subscript x. However, the subscript should not be confused with the representation of the vector observation from the array output, also denoted as a boldface x. From the context the meanings shall be clear to those of skill in the art. Further, quantities relating to the propagating signals impinging on a receive array are denoted with a boldface subscript r. Following this convention, the matrix pencil of the array output data is given as is given as [0052] equation 1. An assumption is made that the received signals r comprising the vector observation of the array output x are independent. Therefore the spatial fourth-order cumulant matrix pencil (SFOCMP) of the array output Px can be written as: P x ( λ , τ ) = C x 4 ( 0 , 0 , 0 ) - λ C x 4 ( τ 1 , τ 2 , τ 3 ) ( 1 )
    Figure US20030228017A1-20031211-M00001
  • where the arguments of the pencil P[0053] x represent a generalized eigenvalue, λ, and a triplet of time delays, τ. The theoretical set of finite generalized eigenvalues turns out to be the inverse of the normalized fourth-order autocumulants of the M signals, { r i ( t ) } i = 1 M
    Figure US20030228017A1-20031211-M00002
  • in the field of view (FOV) during the observation interval. The terms C[0054] x 4 represent the spatial fourth-order autocumulant matrices. The arguments of the terms indicate the triplet of time delays used to form the matrices. The explicit computation is given as [ C x 4 ( τ 1 , τ 2 , τ 3 ) ] rc = i = 1 N Cum [ x i * ( t - τ 1 ) , x i ( t - τ 2 ) , x r ( t ) , x c * ( t - τ 3 ) ]
    Figure US20030228017A1-20031211-M00003
  • where the matrix is N×N, and the subscript rc indicates the element in the r[0055] th row and the cth column. The subscript on the function x in the argument on the right-hand side, indicates which array port, i,r,c=1,2 . . . , N, is being used.
  • Because of the unique definition of the pencil of the array output data P[0056] x is related to the pencil of the impinging signals Pr as given in equation 2 P x ( λ , τ ) = C x 4 ( 0,0,0 ) - λ C x 4 ( τ 1 , τ 2 , τ 3 ) = V [ C r 4 ( 0,0,0 ) - λ C r 4 ( τ 1 , τ 2 , τ 3 ) ] V H = VP r ( λ , τ ) V H ( 2 )
    Figure US20030228017A1-20031211-M00004
  • The quantity V shown in [0057] equation 2 is a N×Ms matrix composed of the steering vectors for each signal impinging on the array, where N is the number of array ports available to the user and Ms, Ms≦N, is the number of signals. In a very simplistic and idealized case the well-known array propagation vector is a steering vector (i.e., the time delay is represented as phase). However, in general if the array is well-designed (i.e., no grating lobes) and the signals are emitted from non-identical locations, then the matrix V is of full rank. This guarantees an equivalence between the eigen structure of the pencils Pr and Px.
  • Since P[0058] r is a pencil solely of the received signals, and the signals are assumed independent, then by virtue of the properties of cumulants, the pencil Pr is diagonal. This property does not hold true for the pencil formed with the array output data x. However, because of “equivalence” finite eigenvalues of Px are the finite eigenvalues of Pr, access to an exploitable high-order statistical property the eigenstructure of the SFOCMP is available. As introduced here these eigenvalues represent the fourth-order characteristics of each received signal. Specifically, each signal in { r i ( t ) } i = 1 M
    Figure US20030228017A1-20031211-M00005
  • contributes one finite eigenvalue, and it is expressed as the inverse normalized fourth-order autocumulant for that signal as expressed by [0059] equation 3. λ m = c r m 4 ( 0,0,0 ) c r m 4 ( τ 1 , τ 2 , τ 3 ) for m = 1 , 2 , , M ( 3 )
    Figure US20030228017A1-20031211-M00006
  • where the terms [0060] c r m 4
    Figure US20030228017A1-20031211-M00007
  • represent the individual fourth-order cumulant terms for each signal. These terms are actually the diagonal terms of the pencil P[0061] r as shown in equation (4). P r ( λ , τ ) = [ c r 1 4 ( 0,0,0 ) - λ c r 1 4 ( τ 1 , τ 2 , τ 3 ) 0 0 0 c r j 4 ( 0,0,0 ) - λ c r j 4 ( τ 1 , τ 2 , τ 3 ) 0 c r M 4 ( 0,0,0 ) - λ c r M 4 ( τ 1 , τ 2 , τ 3 ) ] ( 4 )
    Figure US20030228017A1-20031211-M00008
  • Thus the GEVD of the two pencils P[0062] x and Pr have the same set of finite solutions for the eigenvalues. The eigenvalues are the terms where the rank of the pencil is reduced. It should be readily apparent that values given by equation (3) are the eigenvalues of the pencil equation (1).
  • These eigenvalues are available to an analysis system, and in theory are independent of system Gaussian noise level given sufficient length data records. The eigenvalues are implicit characteristics of the signals carrying the emitter's covert message in each symbol duration. To exploit this property, as mentioned before, the receiver will typically form blocks or batches of received data for the purpose of correlating the eigenstructure over time to determine patterns of persistent values (FIG. 3) augmented by the availability of spatial data. It is important to recall that only the time duration of the emitter's statistical characteristic as measured by the SFOCMP is relevant, and not the exact values. Hence, the emitter is completely free to choose the “carrier waveforms” at will (FIGS. 6[0063] a-f).
  • The steering vectors can be estimated from the cumulant data for each signal in the FOV of the receiver. A cumulant matrix formed by the receive data, say C[0064] x 4(0,0,0) and for each eigenvector available from the pencil Px forms, C x 4 ( 0,0,0 ) e x ( j ) = [ i = 1 M c r i ( 0,0,0 ) v i v i H ] e x ( j ) = β v i
    Figure US20030228017A1-20031211-M00009
  • The last equality follows directly from the fact that each eigenvector of the SFOCMP P[0065] x is orthogonal to each signals steering vector, vi Hex (j)=0 when i≠j. This fact is generated by the unique construction of the SFOCMP and the definitions of the cumulants.
  • In FIG. 3, the Blind [0066] Source Separation processor 309 forms and applies a separation Matrix and enumerates the number of sources. As described above, from an array output the spatial 4th order cumulant matrices are estimated and the estimates are used to determine the eigen analysis for the first-order matrix pencil. Signal detection and enumeration providing the number of source is performed and the separation matrix from the pencil eigenvectors is accomplish. The separation matrix is also used to determine spatial variables, such as steering vectors, AoA or geolocation in block 315 of FIG. 3.
  • There are numerous ways the covert transmitter can control the desired characteristic of its emitted waveform, some of which leave lower-order statistics unchanged. For example altering the channel filter (i.e. Nyquist pulse shaping) between maximum phase and minimum phase realizations is undetectable in the second-order domain (i.e. power spectra), but evident in the fourth-order domain as measured by our SFOCMP. Also one could conceive on signaling with kurtosis, a fourth-order statistic, applied in the transmitter. Or, one could simply shift between classic waveforms, for example, BPSK, GMSK, QPSK, QAM, or potentially even just variants (i.e. constellation rotations, different pulse shape filters) of a fixed modulation type. There also exists the possibility that the “carrier waveforms” might be chosen as chaotic to appear more noise-like or designed using numerical techniques and generated using direct synthesis in a transmitter. [0067]
  • In any scheme adopted, the information transfer from the transmitter to the receiver is contained in the duration of the change in the eigenvalues of the SFOCMP, not the particular eigenvalues. Since our technique is independent of the particular eigenvalue, it is independent of the waveforms used by the emitter. Which allows, in principle, any transmitter to make use of the receiver having the capability to exploit fourth-order cumulants. The degree to which a specific emitter wishes to “hide” from say conventional ESM receivers holds the implications on the implementation details of the “carrier waveforms”. [0068]
  • FIG. 6[0069] a is a generalized schematic diagram of an embodiment of a noise signal generator for a waveform independent covert communication system. A noise generator 602 generates temporal dependent non-Gaussian noise. The output of the noise generator 602 is combined by combiner 606 with a carrier signal or waveform source 604. The combiner can be an adder, a mixer, a multiplier, a non-linear device or other type of combiner that facilitate a change in a higher order signal statistic. This list of combiners in not exhaustive and should not be construed to limit the scope of the invention. The combined signal forms a type of baseband signal that is processed (eg. amplified, upconverted, etc.) and transmitted by transmitter 608. The generalized elements in FIG. 6a form a basic framework and facilitates several different embodiment for noise signal generation and transmission of a covert message. The generalized diagram in 6 a as well as embodiment shown in FIGS. 6a-6 e presented to illustrate variations of noise signal generators and not intended to limit the scope of the disclosure.
  • FIG. 6[0070] b is a specific implementation of FIG. 6a with dual signal sources. The noise generator 602 is implemented with a Non-Gaussian Noise Generator (NGNG) 610 and a temporal filter 612 that creates an output with temporal dependence. The signal source 604 is implemented with two unique signal sources 614 a and 614 b. Unique signal sources 614 a and 614 b are selected to create unique spatial high order statistics when combined by the combiner 606 with the output of the noise generator 602. The unique signal sources 614 a and 614 b are connected to the combiner by a switch 618. The switch is driven by the symbol bit stream of the covert message. The switch 618 connects alternate unique signal sources at the conclusion of each successive symbol duration. It is important to recall that only the time duration of the emitter's statistical characteristic as measured by the SFOCMP is relevant, and not the exact values. Therefore the emitter is completely free to choose the carrier waveforms. While two unique signal sources are shown in FIG. 6b, multiple unique signal source can likewise be applied with the switch alternating between the signal sources.
  • FIG. 6[0071] c is another specific implementation of FIG. 6a having dual temporal filters. The noise generator 602 is implement with a non-Gaussian noise generator 610 connected to two unique temporal filters 612 a and 612 b. The unique temporal filters 612 a and 612 b are selected to create unique spatial high order statistics of the transmitted combined signal (e.g. maximum and minimum phase filters). The unique temporal filters 612 a and 612 b are connected by switch 618 to the combiner 606 and combined with the output of the signal source 604. Again the switch is driven by the symbol bit stream of the covert message in the same manner described above for the dual signal source implementation. The signal or waveform transmitted from transmitter 608 of FIG. 6c provides alternating unique spatial high order statistics for each successive message symbol independent of the symbol transmitter.
  • FIG. 6[0072] d is a specific implementation of the noise generator 602 including two unique non Gaussian noise generators 610 a and 610 b connected by switch 618 to the temporal filter 612. Again the output of the noise generator combined with the signal source 604 in combiner 606 and transmitted by emitter or transmitter 608. The unique NGNGs 610 a and 610 b are likewise selected for their effect on spatial high order statistic of the transmitter signal.
  • FIG. 6[0073] e is an implementation similar to the noise generator 602 of FIG. 6d, where a non-Gaussian noise generator 610 a and temporal filter 612 a form a unique branch and NGNG 610 b and temporal filter 612 b form another unique branch. The branches are unique in the sense of their respective spatial-high order statistics. The branches are connected to the combiner 606 by switch 618.
  • Combinations of the specific implementation described and others that should be readily apparent from an understanding of this disclosure are likewise envisioned. FIG. 6F is but one of the many possible combinations. The implementation of the [0074] noise generator 602 of FIG. 6f includes the noise generator implementation of FIG. 6d containing two unique NGNGs 610 a and 610 b coupled with the signal source 604 of FIG. 6b containing two unique signal sources 614 a and 614 b. While the switches 618 a and 618 b are constrained to switching at only at the duration boundaries, their operation advantageously would operate independently. The result of this specific implementation allows for switching between four unique signals.
  • FIG. 7 is an embodiment using a [0075] Laplacian generator 710 and an infinite impulse response low pass filter 712 (e.g. an autoregressive single-pole filter) as the noise generator 602. The Laplacian generator 710 is chosen for leptokurtosis property. We have chosen to form the necessary sequence of random Laplacian variates by selecting a pair of random variants, say x and y, uniformly distributed between 0 and 1 and forming log(x/y) for each sample required. For complex baseband symbols a total of four uniformly distributed variates must be sampled providing a stream of complex samples where the real and imaginary parts of each are independently Laplacian distributed. The output of the noise generator 602 is combined by combiner 606 with the signal source 604. The signal source 604 is implemented with multiple unique signal sources 714 1, 714 2, . . . 714 1 that are used to switch spatial high order statistics (i.e. matrix pencil eignevalues for the present embodiment) for each symbol of the covert messages. The unique signal sources 714 shown are offset in frequency to provide unique matrix pencil eigenvalues. In this particular embodiment, L must be greater than or equal to the alphabet size M.
  • As shall be understood by those of skill in the art, the specific example discussed above may be extended to use random mappings of frequency offsets over time. Also, we could alter the channel filter. There is no requirement that the filter be IIR as shown in the Figure. A number of alternative implementations could be chosen depending on the application. The key feature of the filter is to introduce a temporal dependence of the input noise waveform. Further one could also consider altering the input noise generator. However, a consideration is to select a source with suitable fourth-order properties. Any or all of these parameters can be modified to control the fourth-order properties for the transmitted waveform so long as the “codebook” constraints (time duration and alphabet size) are maintained. A natural alternative to frequency shifting, would be to pulse the carrier on/off. However, this approach reduces the number of signal samples available for geolocation given a fixed observation interval as discussed hereafter. [0076]
  • If one wished to use standard waveforms as the “carrier waveforms” this mode of operation is also possible with this invention. The transmitter shown in FIG. 7 could be modified to look like alternatives shown in FIGS. 6[0077] b-6 f. It is also important to notice that any deviations from a nominal waveform type, such as a QPSK waveform without phase noise or timing or I/Q imbalance, will cause a detectable shift in a signal's fourth-order statistic. In principle, the basic waveform need not be altered from QPSK to say BPSK, but uniqueness can be achieve by controlling the parameters of, for example, timing jitter, “carrier waveform” symbol rate, phase noise, pulse shaping and the like to provide the necessary separation in the fourth-order domain or other high order domains. The “carrier waveforms” can contain completely worthless symbols, hence anyone trying to conventionally exploit the “carrier waveform” information would be wasting their time, alternatively the “carrier waveforms” may carry information with a secondary message. However, a secondary message contained within the carrier waveform would not necessarily be LPE, LPI or LPD.
  • An example of a potential message recovery embodiment is shown in FIG. 8. To recover a message the receiver requires three parameters relative to the time durations of the message symbols. The receiver must know the [0078] minimum symbol duration 801, maximum symbol duration 803, and some other information regarding the other M-2 symbols in the alphabet. In the embodiment shown in FIG. 8 only a binary symbol alphabet is used. This information can be provided, say by policy where the symbol duration increments are integer-related to the minimum symbol duration 801. In this case the receiver need only know the maximum duration, the minimum duration and alphabet sizes. Otherwise the receiver can just know a-priori the pre-selected durations which can then be arbitrarily selected. The symbol durations should be selected in an application for easing the decision process in the receiver to map the measured durations to one of the alphabet elements. In this case we used a symbol twice as long as another, but such widely disparate time durations may impact throughput. A minimum symbol duration threshold 807 and a maximum symbol duration threshold 811 can be used as a time gate to remove spurious or interfering tracks. Spatially correlated signal tracks 823, 825, 827 and 829 are examples of eigenvalue durations which exceed the threshold in the time gate and thus are ignored or discarded.
  • Spatially correlated eigenvalue time durations [0079] 805 9, 805 16, 805 36, 805 53, 805 68 and 805 80 are sequenced in signal track 821. The eigenvalue time duration of signal track 821 are compared to a decision threshold 809 to map and recover the encoded message “101011” of signal track 821. Where a M-ary alphabet is used M−1 decision thresholds are required.
  • An illustration of a portion of the block-to-block eigenvalue correlation result is shown in FIG. 9. For [0080] block 30 of the covert emitter example in FIG. 4, the large complex plane diagram in the top portion of the FIG. 9 shows the complex eigenvalue locations of the SFOCMP (GEVD) results and the predicted locations of the blockwise eigenvalue correlator. The legend identifies the four levels of eigenvalue correlation confidence, (“New,” “Tentative,” “Candidate,” and “Confirmed”). The five consistent signal eigenvalues of the four co channel interferers and of the covert emitter are indicated by the smaller rectangle. The inconsistent non-signal eigenvalue outside this box tends to move about the complex plane in an erratic/unstable fashion from one block to the next. In this block 30 of this example, the covert emitter eigenvalue lies below the imaginary axis. Blocks 30 and 31 correspond to two different symbols. The lower portion of FIG. 9 illustrates the blockwise changes in eigenvalue locations over blocks 30 to 34 where the message symbol “boundary” occurs between blocks 30 and 31. For blocks 31 to 34, the covert emitter eigenvalue lies in a consistent location above the imaginary axis of the complex plane. The stronger co channel interference source eigenvalues, within the oval for block 34, lie in consistent locations over all five blocks. In a case involving pulsed interference emitters, the consistency of the spatial locations of the interference emitters and covert emitter allow message recovery embodiment to “stitch” together eigenvalues that fall within a “time-gate” and come from the same spatial location.
  • As mentioned above, using a simple time-gating operation in the receiver, it is possible to determine which eigenvalues are potentially information carrying. By correlating the GEVD over many blocks of data the persistence of the eigenvalues can be measured. The persistence of eigenvalues of the SFOCMP over time from the covert transmitter provides the signaling mechanism. However, there may be a number of extraneous pulsed signals in the FOV time coincident with the desired communication signal. This makes message recovery complicated, though with proper message construction and error recovery/correction (i.e. FEC), the system is robust to several types of errors such as “erasures” when ambiguous results may be obtained in the decoding and symbol recovery errors. These results can be encountered due to signal fades (i.e. erasures) or symbol recovery errors in the receiver due to statistical fluctuations in time duration measurements exceeding a tolerable threshold. We can correct improper decisions regarding the detection of a symbol in the message in the receiver using typical error control coding. [0081]
  • An embodiment of the [0082] receiver 303 is shown in FIG. 10. The receiver 303 augments the message recovery process by using spatial information regarding the locations of the transmitters in the field of view. In FIG. 10 signal detections are performed in processor 1001 which is composed of a BSS processor 309 and (eigenvalue) tracking processor 313 previously described. The results of the blockwise processing and eigenvalue correlation previously described are stored in the Signal History Database 317. The data envisioned for the database are signal characteristics measured which include but are not limited to eigenvalue identification number, time initialized, time deleted, number of observations, an estimate of the eigenvalue, and importantly an estimate of the spatial location associated with the eigenvalue. To recover the message in the communication receive processor 319, it is advantageous to initially time gate (filter) all the detections (i.e. durations of the eigenvalues) to limit the scope of the message recovery process. The detections surviving the time gate test are stored in the signal history database 317 along with other associated eigenvalue characteristics, including but not limited to Signal track identification, Symbol estimate, Time initialized, time deleted, number of observations and importantly spatial location. The detections are assembled into strings of symbols (tracks) using the criteria that a message must associate with a “consistent” spatial location. The preferred spatial information is angle-of-arrival (AoA) or geolocation for each signal detected in each block which are determined in the AoA/Geo computation processor 315. This requires that the receiver platform have knowledge of the N-port sensor array calibration and the platform position during receipt of the covert message, such knowledge would normally be available. For simplicity of presentation, from hereon the term “spatial location” or spatial variable refers to AoA if only calibration is known, or geolocation if both receiver position and array calibration are known. It is preferable that the spatial relationship between the transmitter 302 and receive array 327 be fixed over the message length, however, a slowly varying spatial property can be accommodated by introducing a method to correlate a sequence of spatial variables across contiguous blocks of data in the receiver 303 specifically in the Message Sequence Correlator 323 of the Covert Communication Processor 319.
  • Spatial correlation can be broaden to include simply steering vectors. This is useful when the array and transmitter have a stable geometry. Relative motion between the transmitter and sensing array causes the steering vectors to have a detrimental time dependency. Again, if the spatial variable for correlating the message data is “slowly” varying then small incremental changes can be tolerated. The covert messages is indented to be recovered using a “consistency” of the spatial domain information of the computed eigenstructure from the SFOCMP for the signals of interest. But to account for the possibility of “fixed” location emitters and other emitter who are not of interest to the communication process, a time-gate decision process as noted earlier is advantageously applied. This way the receiver need only attempt to decode “message strings” that emanate from “consistent” spatial locations with the appropriate time character. [0083]
  • Although access to the spatial variables using only the receive array output data has been previously described. It is useful to note a blind source separation algorithm based on a fourth-order cumulant matrix pencil produces eigenvectors that are orthogonal all but one signal's steering vector. Thus using the eigenvectors it is possible to estimate each corresponding signal's steering vector. Two methods are possible. The first method is to use the blockwise estimates directly available from the BSS process as described in relation to FIG. 3 such as averaging across the blocks. A second method is to use the time stamps available from the [0084] blockwise correlation process 313 of particular eigenvalues as indexes into the digitizer buffer memory. Subsequently, the raw data so indexed is used to directly make an estimate of the steering vector. In both cases the steering vector will be estimated using the data available from the GEVD of the SFOCMP. Once steering vector estimates are available, the estimation of the other spatial variables, AoA and/or geolocation, can be determine by methods well-known in the art. The characteristic that allows this computation is that the eigenvalues and eigenvectors of the GEVD of the SFOCMP have a 1:1 correspondence as in standard eigenanalysis. So, when the eigenvalues are used to measure the high-order statistical properties of the received signals, an index relating directly to where that signal with that characteristic emanated from is available. Again, the spatial dimension for signal association can be exploited to remove any ambiguity of the temporal decomposition, since we assume that no two emitters are identically located.
  • Signal tracks [0085] 1021, 1023 and 1025 containing associating matrix pencil eigenvalues from consistent sources as evident from common spatial variables are used to recover the message by correlating track information against spatial variables and time gates. This recovery process in performed by the covert communication receive processor 319.
  • The recovered messages are composed of sequentially ordered matrix pencil eigenvalues with a duration within the time gate originating from the same location as determined by the steering vector estimate, AoA or geolocation. The “[0086] signal sequence 1” 1021 includes eigenvalues 1, 4, 7, 10, and 12 with durations mapping message “01011”, these eigenvalues and their associated signals all originated from a “consistent” spatial location, namely, Geo1, AoA1 or had the same steering vector SV1. Similarly, “signal sequence 2” 1023 includes eigenvalues and durations with common steering vectors, AoA, or geolocation and the durations pass the time gate filter thus recovering a “000100” message. “Signal sequence 3” 1025 however had two eigenvalues 3 and 9 with high spatial correlation, however the signal is not decoded since the durations exceed the time gate allowable. In short, the receive processor recognized that the time durations in sequence 1025 do not correspond to any symbols of the expected alphabet.
  • The advantages of incorporating spatial variables into the message recovery process warrant explanation. First, the spatial variables aid in rejecting extraneous pulsed emitters based on their spatial locations being anti-correlated over time to the persistent spatial locations of the covert emitter(s). By the same token, spatial variables allow the basic signaling approach to support multiple access of covert emitters without undue burden in the receiver for properly assembling the pulsed message sequences. This is because the additional covert emitters will very likely emanate from resolvable spatial locations, and the receiver can use the consistency of the spatial locations over time to associate the proper message sequence. For each transmitted signal, the message sequence is represented by the time durations of the eigenvalues of the appropriately selected matrix pencil, where we have preferred the SFOCMP approach. [0087]
  • The spatial location of any emitter is independent of the exact value of its corresponding eigenvalues available from the GEVD of the SFOCMP. Lastly, the spatial variables provide additional “distance” in the recovery process, since it is now multi-dimensional. For example, two signals may have very similar eigenvalues. But, if their spatial locations are resolvable by the receiver, and fairly constant, then the eigenvalues corresponding to those spatial locations can be easily assigned. Then the message can be recovered using the time duration of each eigenvalue in the sequence assigned a given spatial location using the same technique as previously described when only a single covert signal was in view. [0088]
  • The ability to resolve spatial location has system implications that are interrelated. Some top-level practical design issues that must be reconciled, are desired proximity of transmitters, expected noise environment, block processing sample support for estimating spatial locations and the eigenvalues available from the SFOCMP, digitizer sample rates, signal bandwidths and center frequencies, aperture design (i.e. element type, size, number of ports, operating frequency), and the like. This is of course in addition to having an appropriate level of calibration and positional knowledge of the receive platform. Many of these considerations are direct carry-overs from standard array-based signal processing systems. [0089]
  • FIG. 11 is a flow chart for covert communication with a [0090] transmitter 1102 and a receiver 1103 using message recovery with spatial information as shown in FIG. 10. The source 1100 provides a data stream to the covert transmitter 1102 as a stream of M-ary alphabet symbols. In block 1104 each M-ary alphabet symbol is assigned a unique duration. Plural waveforms with determinable high order statistic that are constant and non zero are generated in block 1106. One of the plural waveforms is selected in block 1108 and is transmitted for the symbol assigned duration in block 1110. In block 1112 it is determined if the transmitted symbol is the last in the message, if the last symbol has be transmitted the process is returned in block 1116 to standby or other quasi active state pending initiation of a new message. If other symbols remain to be transmitted, a different waveform from that previously transmitted is selected in block 1114 and transmitted for the respective symbols assigned duration in block 1110.
  • The [0091] receiver 1103 receives the transmitted waveform along with environmental and random noise in a multi element array as shown in block 1105. Using BSS and GEVD, matrix pencil eigenvalues (MPE) are obtained in block 1107 and the MPE are blockwise correlated to determine their respective durations in block 1109. A spatial variable or location is determined from the MPE in block 1111, again the spatial variable can be selected from a Steering vector, AoA or geolocation. The MPE detection durations are time gate filtered in block 1113 to reject detections outside of the duration minimum and maximum thresholds. The passed MPE and durations are sequentially sorted by common or consistent spatial variable in to signal tracks in block 1115. The signal tracks are maps are mapped (i.e. their durations are correlated to symbols and their respective assigned durations)to recover the message in block 1117. The message may be processed otherwise transformed to provide the message at the sink 1101.
  • If spatial data is unavailable, say because calibration of the sensing array has been degraded, the communication process can still operate. However, the freedom of waveform selection by the transmitter is reduced. In this case the transmitter must select a specific waveform type and use it exclusively (in a pulsed fashion) over the entire message (FIG. 4). Unfortunately with this implementation option the achievable data rate is reduced because of the need to introduce “deadtime” to define symbol boundaries. Thus it is possible to use a single “carrier waveform” that can be pulsed “ON” for each symbol followed by a period of “OFF” time. In this way the covert transmitter need only use a single waveform and need not modify it's fourth-order cumulant signature. This could be an advantage in systems where additional spatial correlation variables preferred to aid unambiguous assignment of the received eigenvalues are unavailable, since we identify one and only one eigenvalue. But as mentioned this is disadvantageous in a multi-emitter environment. For such an environment we prefer a waveform agile emitter where the “carrier waveform” sequencer logic can be designed to select the “carrier waveform” for a specific duration, controlled by the message symbol, in either a fixed map or in some other manner. The mapping choice would be up to the transmitter designer and need not be known to the receiver. The receiver processing is shown in FIG. 12. [0092]
  • FIG. 12 presents a [0093] receiver 303 for processing covert messages without spatial location or spatial variable information. Processor 1201 contains BSS processing 309 and eigenvalue correlation 313 operating in the same manner as previously described, however without transferring or developing spatial variables. Therefore the eigenvalue characteristic and time stamps for determining eigenvalue duration are stored in the signal history database 317, without a spatial variable. For the embodiment of FIG. 12, the tracks are not sorted by spatial location by rather eigenvalues. Signal tracks contain only detections of the same eigenvalue with a detected duration that passes the time gate. Signal 1221 includes only one eigenvalue with 4 detection durations that pass the time gate. Signal 1221 thus is mapped to recover message “0010”. The durations of the “OFF” period of pulse is irrelevant as long as its duration is greater that 1 block to enable the previous detection to be terminated the tracking correlator 313. Signal tracks 1223 and 1225 are both ignored or discarded for failing the time gate.
  • FIG. 13 is a flow chart for covert communication with a [0094] transmitter 1302 and a receiver 1303 using message recovery without spatial information as shown in FIG. 12. The process is similar to that shown in FIG. 11. The source 1300 provides a data stream to the covert transmitter 1302 as a M-ary alphabet symbols. In block 1304 each M-ary alphabet symbol is assigned a unique duration. However, now a signal waveform with a determinable high order statistic that is constant and nonzero is generated in block 1306. The waveform is transmitted for the symbol assigned duration in block 1310. In block 1312 it is determined if the transmitted symbol is the last in the message, if the last symbol has been transmitted the process is returned in block 1316 to standby or other quasi active state pending initiation of a new message. If other symbols remain to be transmitted, the transmitter is paused, or another unique waveform is transmitter for an small duration t=k in block 1320. The waveform generated in block 1306 is again transmitted for the respective symbol assigned duration in block 1310.
  • The [0095] receiver 1303 receives the transmitted waveform along with environmental and random noise in a multi element array as shown in block 1305. Using BSS and GEVD, matrix pencil eigenvalues are obtained in block 1307 and the MPE are tracked to determine their respective durations in block 1309. The MPE detection durations are time gate filtered in block 1313 to reject detections outside of the duration minimum and maximum thresholds. The passed MPE and durations are sequentially sorted by specific eigenvalue into signal tracks in block 1321. The eigenvalue durations forming the “signal sequences” are mapped (i.e., their durations are correlated to symbols and their respective assigned durations) to recover the message in block 1317 and delivered to the sink 1301.
  • FIG. 14 is a flow diagram for covert communication using encoded eigenvalues and spatial information. The [0096] source 1400 provides a data stream to the covert transmitter 1402 as a M-ary alphabet symbols. In block 1404 the M-ary symbols of the messages are encoded as MPE with unique durations. The eigenvalue is transmitted for the symbol assigned duration as a respective source waveform in block 1410. The source waveform is the low order signal whose high order statistic creates the respective eigenvalue. In block 1412 it is determined if the transmitted symbol is the last in the message, if the last symbol has be transmitted the process is returned in block 1416 to standby or other quasi active state pending initiation of a new message. If other symbols remain to be transmitted, a different MPE is selected in block 1414 and transmitted with a respective symbol assigned duration as a respective source waveform in block 1410.
  • The receiver [0097] 1403 receives the transmitted waveform along with environmental and random noise in a multi element array as shown in block 1405. Using BSS and GEVD, the matrix pencil eigenvalues encoding the symbol are recovered and tracked in block 1407 to determine their respective durations and spatial information in block 1409. The MPEs are decoded based on MPED and Spatial information as previously described in block 1417 to recover the message and provide the message symbols to the sink in block 1401. Generally the noise is white Gaussian noise, color noise or interferer signals.
  • FIG. 15 is a representation of a binary message sequence 101011 transmitted via data carrying waveforms according to an embodiment. In the embodiment shown the symbol are transmitted by alternating waveforms of BPSK, QPSK and GMSK. As the message is independent of the encoded data, for LPI it is preferred that the carrier waveforms be modulated with random data at a rate greatly exceeding the covert message symbol rate. The random data may be modulated for a BPSK waveform at 5 sample/per random symbol with a frequency offset of 0.1 f[0098] s, for QPSK wave at 10 samples/random symbol with a frequency offset of −0.05 fs and for a GMSK waveform (GSM:hmod=½, Btb=⅓) 2 samples/random symbol with a frequency offset of 0.0005 fs. Other waveforms including to DQPSK, DBPSK, FSK, QAM, DPCM are equally applicable and are also envisioned, however embodiments of the invention should not be construed to be limited to the particular waveforms listed.
  • As discussed previously the waveform duration discriminates the message symbol. As shown in FIG. 15 the [0099] symbol 1 has a duration of 20 blocks and the symbol “0” has a duration of 10 blocks. The message is transmitted as a GMSK waveform for a duration of 20 blocks indicating a symbol “1”. The next symbol is transmitted as a different waveform QPSK for a duration of 10 block indicating a symbol “0”. The next symbol is transmitted by the GMSK waveform for a duration of 20 blocks again indicating a “1” symbol and a QPSK waveform is used to transmit the “0” symbol. As illustrative of the independent of the waveform and message content the symbol “1” is then transmitted as a BPSK waveform for a 20 block duration and the next symbol “1” is transmitted by the QPSK waveform for a duration of 20 block. Evident in FIG. 15 is that identical waveforms can communicate different symbols while maintaining waveform content and message independence.
  • While preferred embodiments of the present inventive system and method have been described, it is to be understood that the embodiments described are illustrative only and that the scope of the embodiments of the present inventive system and method is to be defined solely by the appended claims when accorded a full range of equivalence, many variations and modifications naturally occurring to those of skill in the art from a perusal hereof. [0100]

Claims (150)

We claim:
1. A communication device for providing a communication signal from a data signal comprised of a sequence of bits, comprising:
a noise generator for generating a first signal comprised of a bit stream;
a filter for receiving said first signal and generating a temporal dependence between immediately adjacent bits in said first signal to thereby provide a third signal;
a plurality of signal sources each providing a source signal with a characteristic unique from the others;
a switch for selectively operatively connecting as a function of said data signal one of said plurality of signal sources in sequence to a combining means to thereby provide a second signal; and
said combining means for combining said third signal with said second signal to thereby provide the communication signal from said data signal.
2. The communication device of claim 1 wherein said noise generator is a Laplacian noise generator.
3. The communication device of claim 1 wherein said noise generator generates a signal with a distribution function for which the kurtosis is greater than the kurtosis of a Gaussian distribution function.
4. The communication device of claim 1 wherein said noise generator generates a signal with a distribution function for which the kurtosis is different than the kurtosis of a Gaussian distribution function.
5. The communication device of claim 1 wherein said filter is a spectral limiting filter.
6. The communication device of claim 1 wherein said filter is selected from the group consisting of an autoregressive filter, an infinite impulse response filter combined with a low pass filter, and a finite impulse response filter with a moving average.
7. The communication device of claim 1 wherein said characteristic is from the group consisting of frequency, waveform, matrix pencil eigenvalue, phase noise, I/Q imbalance, timing jitter, phase jitter, symbol rate, pulse shape, a fourth-order statistic, relative rotational alignment of a signal constellation, power amplifier rise/fall time, and Doppler shift.
8. The communication device of claim 1 wherein said characteristic is a waveform from the group consisting of BPSK, QPSK, GMSK, QAM, DBPSK MFSK, FSK, and DQPSK waveforms.
9. The communication device of claim 1 wherein said characteristic is a waveform with a matrix pencil eigenvalue that is a function of the spatial second order moment of said waveform.
10. The communication device of claim 1 wherein said characteristic is a waveform with a matrix pencil eigenvalue that is a function of the spatial fourth order cumulant of said waveform.
11. The communication device of claim 1 wherein at least one of said plurality of signal sources is from the group consisting of an oscillator, a random signal generator, and a Laplacian noise generator.
12. The communication device of claim 1 wherein said combining means is a logic circuit.
13. The communication device of claim 1 wherein said combining means is from the group consisting of a mixer, a multiplier, a divider, and an adder.
14. The communication device of claim 1 wherein said communication signal has a received SNR less than 0 dB.
15. The communication device of claim 1 wherein said communication signal has a SNIR less than 0 dB at the intended receiver.
16. The communication device of claim 1 further comprising circuitry for forming symbols of an M-ary alphabet by assigning sequentially to each of said symbols a predetermined number of bits of said data signal.
17. The communication device of claim 16 wherein said predetermined number of bits is from the group consisting of positive integers.
18. The communication device of claim 16 wherein said communication signal is comprised of symbols of said M-ary alphabet.
19. The communication device of claim 18 wherein each of said symbols in said M-ary alphabet has a time duration unique among said alphabet.
20. The communication device of claim 19 wherein for each of said symbols said switch operatively connects one of said plural signal sources other than the signal source operatively connected for the immediately preceding symbol.
21. The communication device of claim 20 wherein said switch operatively connects one of said plural signal sources for the duration of the respective symbol.
22. The communication device of claim 1 wherein said first signal includes an I signal and a Q signal.
23. The communication device of claim 1 further comprising radio frequency processing circuitry for processing said communication signal for transmission.
24. The communication device of claim 23 further comprising a radio frequency transmitter and antenna for transmitting said communication signal.
25. A communication device for providing a communication signal from a sequence of bits comprising:
first means for providing a first signal;
second means for providing a second signal comprising, in sequence, as a function of said sequence of bits, one of a plurality of third signals each with a matrix pencil eigenvalue unique among said third signals; and
third means for combining said first and second signals to thereby provide the communication signal.
26. The communication device of claim 25 further comprising circuitry for forming symbols of an M-ary alphabet by assigning sequentially to each of said symbols a predetermined number of bits in said sequence of bits.
27. The communication device of claim 26 wherein said predetermined number of bits is from the group consisting of positive integers.
28. The communication device of claim 26 wherein said communication signal is comprised of symbols of said M-ary alphabet.
29. The communication device of claim 28 wherein each of said symbols in said M-ary alphabet has a time duration unique among said alphabet.
30. The communication device of claim 29 wherein said first means comprises a noise generator for providing a noise signal of which said first signal is a function thereof.
31. The communication device of claim 30 wherein said first means further comprises a filter for receiving said noise signal and producing therefrom a signal with temporal dependence between adjacent bits of said noise signal to thereby provide said first signal.
32. The communication device of claim 31 wherein said second means comprises a plurality of signal sources each providing one of said plurality of third signals.
33. The communication device of claim 32 wherein said second means further comprises a first switch for operatively connecting for each of said symbols one of said plural signal sources other than the signal source operatively connected for the immediately preceding symbol to thereby provide said second signal.
34. The communication device of claim 33 wherein said noise generator is a plurality of noise generators the outputs of which are operatively connected to a second switch which operates synchronously with said first switch to thereby provide said noise signal.
35. The communication device of claim 33 wherein said noise generator is a plurality of noise generators the outputs of which are operatively connected to a second switch which operates asynchronously with said first switch to thereby provide said noise signal.
36. The communication device of claim 33 wherein said filter is a plurality of filters the inputs of which are operatively connected to a second switch which receives said noise signal and which operates synchronously with said first switch to thereby provide said first signal.
37. The communication device of claim 33 wherein said filter is a plurality of filters the inputs of which are operatively connected to a second switch which receives said noise signal and which operates asynchronously with said first switch to thereby provide said first signal.
38. The communication device of claim 33 wherein said noise generator is a plurality of noise generators and said filter is a plurality of filters whereby each noise generator is operatively connected to one of said plurality of filters, wherein the outputs of said filters are operatively connected to a second switch which operates synchronously with said first switch to thereby provide said first signal.
39. The communication device of claim 33 wherein said noise generator is a plurality of noise generators and said filter is a plurality of filters whereby each noise generator is operatively connected to one of said plurality of filters, wherein the outputs of said filters are operatively connected to a second switch which operates asynchronously with said first switch to thereby provide said first signal.
40. The communication device of claim 29 wherein said first means comprises a signal source for providing said first signal.
41. The communication device of claim 40 wherein said second means comprises a plurality of noise generators each providing one of said plurality of third signals.
42. The communication device of claim 41 wherein said second means further comprises a first switch for operatively connecting for each of said symbols one of said plural noise generators other than the noise generator operatively connected for the immediately preceding symbol to thereby provide said second signal.
43. The communication device of claim 41 wherein each of said noise generators includes a filter for providing temporal dependence between immediately adjacent bits in the associated said third signal.
44. The communication device of claim 29 further comprising a filter for receiving said communication signal and producing therefrom a second communication signal with temporal dependence between adjacent symbols of said second communication signal.
45. The communication device of claim 29 wherein said first signal comprises a noise signal with a Laplacian distribution with temporal dependence between adjacent bits of said noise signal.
46. The communication device of claim 29 wherein said first signal comprises a noise signal with temporal dependence between adjacent bits of said noise signal and with a distribution function for which the kurtosis is greater than the kurtosis of a Gaussian distribution function.
47. The communication device of claim 29 wherein said first signal comprises a noise signal with temporal dependence between adjacent bits of said noise signal and with a distribution function for which the kurtosis is different than the kurtosis of a Gaussian distribution function.
48. In a receiver including an antenna with a plurality of antenna elements for receiving a communication signal comprised of a plurality of symbols and a digitizer for providing a bit stream from the received symbols, the improvement comprising a filter for generating a temporal dependence between immediately adjacent bits of said bit stream.
49. The receiver of claim 48 wherein the improvement further comprises means for determining a matrix pencil eigenvalue for at least one of said symbols.
50. The receiver of claim 49 wherein the improvement further comprises means for determining the generalized eigenvalue decomposition of said matrix pencil eigenvalue.
51. The receiver of claim 50 wherein the improvement further comprises means for determining the duration of said matrix pencil eigenvalue.
52. The receiver of claim 51 wherein the improvement further comprises means for filtering said symbols as a function of the determined duration of the associated matrix pencil eigenvalue.
53. The receiver of claim 52 wherein the improvement further comprises means for spatially correlating the filtered matrix pencil eigenvalues.
54. The receiver of claim 53 wherein the improvement further comprises means for sequentially mapping the correlated matrix pencil eigenvalues.
55. A communication system comprising a transmitter and a receiver geographically separated from the transmitter for communicating a communication signal formed from a sequence of bits wherein said communication signal is comprised of a plurality of symbols, comprising:
a transmitter comprising:
first means for providing a first signal;
second means for providing a second signal comprising in sequence as a function of said sequence of bits one of a plurality of third signals each with a matrix pencil eigenvalue unique among said third signals;
third means for combining said first and second signals to thereby provide said communication signal; and
transmitting means for transmitting said communication signal; and
a receiver comprising:
receiving means for receiving and digitizing the symbols of said communication signal to thereby produce a first received signal;
means for determining a matrix pencil eigenvalue for at least one of said first received signal symbols;
means for determining the generalized eigenvalue decomposition of said matrix pencil eigenvalue;
means for determining the duration of said matrix pencil eigenvalue; [ditto]
means for filtering said symbols as a function of the determined duration of the associated matrix pencil eigenvalue; [ditto]
means for spatially correlating the filtered matrix pencil eigenvalues; and
means for sequentially mapping the correlated matrix pencil eigenvalues to thereby determine at the receiver the sequence of bits in the communication signal.
56. The communication system of claim 55 further comprising at the receiver means for determining a spatial variable of said matrix pencil eigenvalue.
57. The communication system of claim 56 wherein said spatial variable is a steering vector.
58. The communication system of claim 57 further comprising at the receiver means for determining the angle of arrival of the communication signal as a function of said steering vector.
59. The communication system of claim 58 further comprising at the receiver means for determining the geolocation of said transmitter as a function of said steering vector.
60. A communication system comprising a transmitter and a receiver geographically separated from the transmitter for communicating a communication signal formed from a sequence of bits wherein said communication signal is comprised of a plurality of symbols, comprising:
a transmitter comprising:
first means for providing a first signal;
second means for providing a second signal comprising in sequence as a function of said sequence of bits one of a plurality of third signals each with a matrix pencil eigenvalue unique among said third signals;
third means for combining said first and second signals to thereby provide said communication signal; and
transmitting means for transmitting said communication signal; and
a receiver comprising:
receiving means for receiving and digitizing the symbols of said communication signal to thereby produce a first received signal;
means for determining a matrix pencil eigenvalue for at least one of said first received signal symbols;
means for determining the generalized eigenvalue decomposition of said matrix pencil eigenvalue;
means for determining the duration of said matrix pencil eigenvalue;
means for filtering said symbols as a function of the determined duration of the associated matrix pencil eigenvalue;
means for sorting the filtered matrix pencil eigenvalues into at least one group of like eigenvalues; and
means for sequentially mapping the correlated matrix pencil eigenvalues to thereby determine at the receiver the sequence of bits in the communication signal.
61. A method for providing a communication signal from a data signal comprised of a sequence of bits, comprising the steps of:
providing a first signal comprised of a bit stream;
receiving said first signal and generating a temporal dependence between immediately adjacent bits in said first signal to thereby provide a third signal;
providing a plurality of signal sources each providing a source signal with a characteristic unique from the others;
providing a combining means;
selectively operatively connecting as a function of said data signal one of said plurality of signal sources in sequence to said combining means to thereby provide a second signal; and
combining at said combining means said third signal with said second signal to thereby provide the communication signal from said data signal.
62. The method of claim 61 wherein said first signal has a Laplacian distribution function.
63. The method of claim 61 wherein said first signal has a distribution function for which the kurtosis is greater than the kurtosis of a Gaussian distribution function.
64. The method of claim 61 wherein said first signal has a distribution function for which the kurtosis is different than the kurtosis of a Gaussian distribution function.
65. The method of claim 61 wherein said characteristic is from the group consisting of frequency, waveform, matrix pencil eigenvalue, phase noise, I/Q imbalance, timing jitter, phase jitter, symbol rate, pulse shape, a fourth-order statistic, relative rotational alignment of a signal constellation, power amplifier rise/fall time, and Doppler shift.
66. The method of claim 61 wherein said characteristic is a waveform from the group consisting of BPSK, QPSK, GMSK, QAM, DBPSK MFSK, FSK, and DQPSK waveforms.
67. The method of claim 61 wherein said characteristic is a waveform with a matrix pencil eigenvalue that is a function of the spatial second order moment of said waveform.
68. The method of claim 61 wherein said characteristic is a waveform with a matrix pencil eigenvalue that is a function of the spatial fourth order cumulant of said waveform.
69. The method of claim 61 wherein said communication signal has a received SNR less than 0 dB.
70. The method of claim 61 wherein said communication signal has a SNIR less than 0 dB at the intended receiver.
71. The method of claim 61 further comprising the step of forming symbols of an M-ary alphabet by assigning sequentially to each of said symbols a predetermined number of bits of said data signal.
72. The method of claim 71 wherein said predetermined number of bits is from the group consisting of positive integers.
73. The method of claim 71 wherein said communication signal is comprised of symbols of said M-ary alphabet.
74. The method of claim 73 wherein each of said symbols in said M-ary alphabet has a time duration unique among said alphabet.
75. The method of claim 74 wherein for each of said symbols the step of selectively operatively connecting is accomplished so that one of said plural signal sources other than the signal source for the immediately preceding symbol is operatively connected.
76. The method of claim 75 wherein the step of selectively operatively connecting is accomplished such that one of said plural signal sources is operatively connected for the duration of the respective symbol.
77. The method of claim 71 wherein said first signal includes an in-phase signal and a quadrature signal.
78. The method of claim 71 further comprising the step of processing said communication signal for radio frequency transmission.
79. The method of claim 78 further comprising the step of transmitting said communication signal.
80. A method for providing a communication signal from a sequence of bits comprising:
providing a first signal;
providing a second signal comprising, in sequence, as a function of said sequence of bits, one of a plurality of third signals each with a matrix pencil eigenvalue unique among said third signals; and
combining said first and second signals to thereby provide the communication signal.
81. The method of claim 80 further comprising the step of forming symbols of an M-ary alphabet by assigning sequentially to each of said symbols a predetermined number of bits in said sequence of bits.
82. The method of claim 81 wherein said predetermined number of bits is from the group consisting of positive integers.
83. The method of claim 81 wherein said communication signal is comprised of symbols of said M-ary alphabet.
84. The method of claim 83 wherein each of said symbols in said M-ary alphabet has a time duration unique among said alphabet.
85. The method of claim 84 wherein said first signal is comprised of a noise signal of which said first signal is a function thereof.
86. The method of claim 85 wherein said noise signal is produced with temporal dependence between adjacent bits to thereby provide said first signal.
87. The method of claim 86 wherein said second signal is comprised of a sequential combination of a plurality of third signals provided by a corresponding plurality of signal sources.
88. The method of claim 87 wherein the step of providing said second signal is accomplished by operatively connecting said third signals to a first switch whereby for each of said symbols said switch operatively connects one of said third signals other than the third signal that was operatively connected for the preceding symbol.
89. The method of claim 88 wherein said noise signal is a plurality of noise signals which are operatively connected to a second switch which operates synchronously with said first switch to thereby provide said noise signal.
90. The method of claim 88 wherein said noise signal is a plurality of noise signals which are operatively connected to a second switch which operates asynchronously with said first switch to thereby provide said noise signal.
91. The method of claim 88 wherein said temporal dependence between adjacent bits of said noise signal is provided by a plurality of filters arranged in parallel the inputs of which are operatively connected to a second switch which receives said noise signal and which operates synchronously with said first switch to sequentially operatively connect one of said plurality of filters to thereby provide said first signal.
92. The method of claim 88 wherein said temporal dependence between adjacent bits of said noise signal is provided by a plurality of filters arranged in parallel the inputs of which are operatively connected to a second switch which receives said noise signal and which operates asynchronously with said first switch to sequentially operatively connect one of said plurality of filters to thereby provide said first signal.
93. The method of claim 88 wherein said noise signal is a plurality of noise signals and said temporal dependence between adjacent bits of said noise signal is provided by a corresponding plurality of filters arranged in parallel whereby each said noise signal is operatively connected to one of said plurality of filters wherein the outputs of said filters are operatively connected to a second switch which operates synchronously with said first switch to thereby provide said first signal.
94. The method of claim 88 wherein said noise signal is a plurality of noise signals and said temporal dependence between adjacent bits of said noise signal is provided by a corresponding plurality of filters arranged in parallel whereby each said noise signal is operatively connected to one of said plurality of filters wherein the outputs of said filters are operatively connected to a second switch which operates asynchronously with said first switch to thereby provide said first signal.
95. The method of claim 84 wherein said first signal is provided by a signal source.
96. The method of claim 95 wherein said second signal is provided by a first switch for operatively connecting for each of said symbols one of said third signals other than the noise signal operatively connected for the immediately preceding symbol.
97. The method of claim 96 wherein for each of said third signals the unique matrix pencil eigenvalue is a function of the temporal dependence between immediately adjacent bits provided by a filter.
98. The method of claim 84 further comprising the step of providing temporal dependence between adjacent bits of said communication signal and producing therefrom a second communication signal.
99. The method of claim 98 wherein the temporal dependence is provided by a filter.
100. The method of claim 84 wherein said first signal comprises a noise signal with a Laplacian distribution with temporal dependence between adjacent bits of said noise signal.
101. The method of claim 84 wherein said first signal comprises a noise signal with temporal dependence between adjacent bits of said noise signal and with a distribution function for which the kurtosis is greater than the kurtosis of a Gaussian distribution function.
102. The method of claim 84 wherein said first signal comprises a noise signal with temporal dependence between adjacent bits of said noise signal and with a distribution function for which the kurtosis is different than the kurtosis of a Gaussian distribution function.
103. In a method of receiving a communication signal comprised of a plurality of symbols in a receiver including an antenna with a plurality of antenna elements and a digitizer for providing a bit stream from the received symbols, the improvement comprising the step of providing a temporal dependence between immediately adjacent bits of said bit stream.
104. The method of claim 103 wherein the improvement further comprises the step of determining a matrix pencil eigenvalue for at least one of said symbols.
105. The method of claim 104 wherein the improvement further comprises the step of determining the generalized eigenvalue decomposition of said matrix pencil eigenvalue.
106. The method of claim 105 wherein the improvement further comprises the step of determining the duration of said matrix pencil eigenvalue.
107. The method of claim 106 wherein the improvement further comprises the step of filtering said symbols as a function of the determined duration of the associated matrix pencil eigenvalue.
108. The method of claim 107 wherein the improvement further comprises the step of spatially correlating the filtered matrix pencil eigenvalues.
109. The method of claim 108 wherein the improvement further comprises the step of sequentially mapping the correlated matrix pencil eigenvalues.
110. A method of communicating between a transmitter and a receiver including a multi-element array antenna which is geographically spaced apart from said transmitter a communication signal formed from a sequence of bits wherein said communication signal is comprised of a plurality of symbols, comprising the steps of:
in the transmitter:
providing a first signal;
providing a second signal comprising in sequence as a function of said sequence of bits one of a plurality of third signals each with a matrix pencil eigenvalue unique among said third signals;
combining said first and second signals to thereby provide said communication signal; and
transmitting said communication signal; and
in the receiver:
receiving and digitizing the symbols of said communication signal to thereby produce a first received signal;
determining a matrix pencil eigenvalue for at least one of said first received signal symbols;
determining the generalized eigenvalue decomposition of said matrix pencil eigenvalue;
determining the duration of said matrix pencil eigenvalue;
filtering said symbols as a function of the determined duration of the associated matrix pencil eigenvalue;
spatially correlating the filtered matrix pencil eigenvalues; and
sequentially mapping the correlated matrix pencil eigenvalues to thereby determine at the receiver the sequence of bits in the communication signal.
111. The method of claim 110 wherein the step of determining the matrix pencil eigenvalue for a received symbol includes the step of determining a spatial second-order moment.
112. The method of claim 110 wherein the step of determining the matrix pencil eigenvalue for a received symbol includes the step of determining a spatial fourth-order cumulant.
113. The method of claim 110 further comprising the step of determining a spatial variable of said matrix pencil eigenvalue.
114. The method of claim 113 wherein said spatial variable is a steering vector of the communication signal.
115. The method of claim 114 further comprising at the receiver the step of determining the angle of arrival of the communication signal as a function of said steering vector.
116. The method of claim 114 further comprising at the receiver the step of determining the geolocation of said transmitter as a function of said steering vector.
117. The method of claim 110 wherein said first signal is a noise signal comprised of a bit stream with a Laplacian distribution with temporal dependence between immediately adjacent bits.
118. The method of claim 110 wherein each of said third signals has a unique frequency among said third signals.
119. The method of claim 110 wherein each of said third signals has a unique waveform among said third signals.
120. The method of claim 119 wherein said waveforms are from the group consisting of BPSK, QPSK, GMSK, QAM, DBPSK MFSK, FSK, and DQPSK waveforms.
121. The method of claim 119 further comprising the step of modulating the waveform of at least one of said third signals with random data at a rate substantially exceeding the rate of said sequence of bits.
122. The method of claim 119 further comprising the step of modulating the waveform of at least one of said third signals with data from a third communication signal.
123. The method of claim 110 wherein the step of receiving and digitizing includes the step of organizing the digitized bits into blocks of data comprised of sequential snapshots wherein said snapshots are each comprised of a bit from each element of said multi-element array antenna
124. The method of claim 123 wherein the matrix pencil eigenvalue is determined from said blocks of data.
125. The method of claim 124 wherein said block of data includes approximately 5000 snapshots.
126. The method of claim 125 wherein said duration of matrix pencil eigenvalues is the number of sequential blocks of data with approximately the same matrix pencil eigenvalue.
127. The method of claim 110 wherein the communication signal has a received SNR less than 0 dB.
128. The method of claim 110 wherein the communication signal has a SNIR less than 0 dB at the receiver.
129. The method of claim 110 wherein the communication signal has a SNIR less than −6 dB at the receiver.
130. The method of claim 110 wherein said sequence of third signals is determined at random.
131. The method of claim 110 wherein said sequence of third signals is predetermined.
132. A method of communicating between a transmitter and a receiver including a multi-element array antenna which is geographically spaced apart from said transmitter a communication signal formed from a sequence of bits wherein said communication signal is comprised of a plurality of symbols, comprising the steps of:
in the transmitter:
providing a first signal;
providing a second signal comprising in sequence as a function of said sequence of bits one of a plurality of third signals each with a matrix pencil eigenvalue unique among said third signals;
combining said first and second signals to thereby provide said communication signal; and
transmitting said communication signal; and
in the receiver:
receiving and digitizing the symbols of said communication signal to thereby produce a first received signal;
determining a matrix pencil eigenvalue for at least one of said first received signal symbols;
determining the generalized eigenvalue decomposition of said matrix pencil eigenvalue;
determining the duration of said matrix pencil eigenvalue;
filtering said symbols as a function of the determined duration of the associated matrix pencil eigenvalue;
sorting the filtered matrix pencil eigenvalues into at least one group of like eigenvalues; and
sequentially mapping the correlated matrix pencil eigenvalues to thereby determine at the receiver the sequence of bits in the communication signal.
133. A method of communicating information comprised of a sequence of bits from a transmitter to a receiver comprising the steps of:
modulating a carrier wave with a data signal to thereby produce a first signal wherein the data signal does not represent the information to be communicated;
transmitting the first signal from the transmitter;
receiving the first signal at the receiver;
determining at the receiver the communicated information independently of the data signal.
134. The method of claim 133 wherein the step of determining the communicated information is performed without demodulating the first signal.
135. The method of claim 133 further comprising the step of providing a second signal with temporal dependence between immediately adjacent bits and combining said second signal with said modulated carrier wave to thereby produce said first signal.
136. The method of claim 135 further comprising the steps of:
sequentially grouping a predetermined number of bits in said sequence of bits into symbols to thereby provide symbols of an M-ary alphabet;
providing a plurality of third signals each with a matrix pencil eigenvalue unique among said third signals;
providing a switch whereby for each of said symbols said switch operatively connects one of said third signals other than the third signal that was operatively connected for the preceding symbol to thereby form said second signal.
137. The method of claim 136 wherein each of said symbols in said M-ary alphabet has a time duration unique among said alphabet.
138. The method of claim 137 wherein said information is communicated to the receiver as a function of the durations of said symbols in said M-ary alphabet.
139. In a method of communicating information from a transmitter to a receiver using a communication signal comprised of a carrier wave modulated by a data signal, the improvement comprising determining the information at the receiver independent of the data signal.
140. In a method of communicating information from a transmitter to a receiver using a communication signal containing a plurality of symbols comprised of a carrier wave modulated by a data signal, the improvement comprising determining the information at the receiver as a function of a higher-order statistic of the received symbols.
141. The method of claim 140 wherein said symbols comprise an M-ary alphabet and said higher-order statistic for each symbol in the M-ary alphabet is unique among said alphabet.
142. The method of claim 140 wherein said information includes a spatial variable.
143. In a method of communicating a first information stream from a transmitter to a receiver using a communication signal containing a plurality of symbols comprised of a carrier wave modulated by a data signal, the improvement comprising communicating a second information stream from the transmitter to the receiver by selectively altering a characteristic of said communication signal wherein said second information stream is different than said first information stream.
144. The method of claim 143 wherein said characteristic is from the group consisting of frequency, waveform, matrix pencil eigenvalue, phase noise, I/Q imbalance, timing jitter, phase jitter, symbol rate, pulse shape, a fourth-order statistic, relative rotational alignment of a signal constellation, power amplifier rise/fall time, and Doppler shift.
145. In a method of communication using plural waveforms where the duration of the transmission of each waveform represents a symbol, the improvement wherein the duration of transmission is determined by the state of a high-order statistic of a characteristic of the waveforms.
146. The method of claim 145 wherein said high-order statistic is a spatial second-order moment.
147. The method of claim 145 wherein said high-order statistic is a spatial fourth-order cumulant.
148. A method for communication of a message comprising a plurality of sequenced M-ary alphabet symbols, said method comprising the steps of:
assigning each M-ary alphabet symbol a unique symbol duration;
generating a waveform for which both the spatial 2nd order moment or the spatial 4th order cumulant are constant and not zero; and
for each symbol, transmitting the waveform for the assigned symbol duration immediately followed by a pause in waveform transmission.
149. A method for communication of a message comprising a plurality of sequenced M-ary alphabet symbols, said method comprising the steps of:
assigning each M-ary alphabet symbol a unique symbol duration;
generating plural waveforms for which the spatial 2nd order moment and the spatial 4th order cumulant are constant, nonzero, and unique; and
for each symbol, transmitting one of the plural waveforms for the assigned symbol duration so that immediately adjacent symbols are transmitted by different waveforms.
150. A method for communicating a message comprising a plurality of sequenced M-ary alphabet symbols represented by M-ary durations, said method comprising the steps of:
receiving a signal in a multiple element array;
determining a high order statistic of a signal characteristic of each symbol, and determining the duration of each symbol by the state of the determined high order statistic.
US10/400,486 2002-04-22 2003-03-28 Method and system for waveform independent covert communications Abandoned US20030228017A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US10/400,486 US20030228017A1 (en) 2002-04-22 2003-03-28 Method and system for waveform independent covert communications

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US37414902P 2002-04-22 2002-04-22
US10/400,486 US20030228017A1 (en) 2002-04-22 2003-03-28 Method and system for waveform independent covert communications

Publications (1)

Publication Number Publication Date
US20030228017A1 true US20030228017A1 (en) 2003-12-11

Family

ID=29251149

Family Applications (2)

Application Number Title Priority Date Filing Date
US10/360,631 Expired - Fee Related US6711528B2 (en) 2002-04-22 2003-02-10 Blind source separation utilizing a spatial fourth order cumulant matrix pencil
US10/400,486 Abandoned US20030228017A1 (en) 2002-04-22 2003-03-28 Method and system for waveform independent covert communications

Family Applications Before (1)

Application Number Title Priority Date Filing Date
US10/360,631 Expired - Fee Related US6711528B2 (en) 2002-04-22 2003-02-10 Blind source separation utilizing a spatial fourth order cumulant matrix pencil

Country Status (8)

Country Link
US (2) US6711528B2 (en)
EP (1) EP1500007B1 (en)
JP (1) JP4271041B2 (en)
CN (1) CN100386764C (en)
AT (1) ATE487187T1 (en)
AU (1) AU2003223605A1 (en)
DE (1) DE60334790D1 (en)
WO (1) WO2003090127A1 (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050118963A1 (en) * 2003-12-02 2005-06-02 Mao-Ching Chiu Method and apparatus for I/Q imbalance calibration of a transmitter system
US20080001821A1 (en) * 2004-09-14 2008-01-03 Akira Tanaka Signal Arrival Direction Deducing Device, Signal Arrival Direction Deducing Method, and Signal Direction Deducing Program
US20090003418A1 (en) * 2005-06-22 2009-01-01 Karabinis Peter D Systems, Methods, Devices And/Or Computer Program Products For Providing Communications Deviod Of Cyclostationary Features
US20090268783A1 (en) * 2005-06-22 2009-10-29 Eices Research, Inc. Wireless communications systems and/or methods providing low interference, high privacy and/or cognitive flexibility
US20130096699A1 (en) * 2010-06-21 2013-04-18 Optimized Systems And Solutions Limited Asset health monitoring
US8537916B2 (en) 2010-03-29 2013-09-17 Eices Research, Inc. Increased capacity communications for OFDM-based wireless communications systems/methods/devices
US8576940B2 (en) 2005-06-22 2013-11-05 Eices Research, Inc. Systems/methods of adaptively varying a bandwidth and/or frequency content of communications
WO2014066279A1 (en) * 2012-10-22 2014-05-01 Gassi Donatello Apelusion Information security based on eigendecomposition
US9374746B1 (en) 2008-07-07 2016-06-21 Odyssey Wireless, Inc. Systems/methods of spatial multiplexing
US9806790B2 (en) 2010-03-29 2017-10-31 Odyssey Wireless, Inc. Systems/methods of spectrally efficient communications
US20180007045A1 (en) * 2016-06-30 2018-01-04 Mehdi Arashmid Akhavain Mohammadi Secure coding and modulation for optical transport
CN108702219A (en) * 2016-03-08 2018-10-23 骁阳网络有限公司 Reduce the adaptive constellation that phase distortion influences
CN109547371A (en) * 2019-01-17 2019-03-29 西安电子科技大学 Signal processing method based on quasi- high power spectrum
USRE47633E1 (en) 2005-06-22 2019-10-01 Odyssey Wireless Inc. Systems/methods of conducting a financial transaction using a smartphone
CN114095045A (en) * 2022-01-24 2022-02-25 成都航天通信设备有限责任公司 Communication method, device and system based on constellation pseudo-shape modulation

Families Citing this family (101)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE521024C2 (en) * 1999-03-08 2003-09-23 Ericsson Telefon Ab L M Method and apparatus for separating a mixture of source signals
US7952511B1 (en) 1999-04-07 2011-05-31 Geer James L Method and apparatus for the detection of objects using electromagnetic wave attenuation patterns
US6711528B2 (en) * 2002-04-22 2004-03-23 Harris Corporation Blind source separation utilizing a spatial fourth order cumulant matrix pencil
JP4309110B2 (en) * 2002-09-27 2009-08-05 パナソニック株式会社 Adaptive antenna wireless communication device
US7187326B2 (en) * 2003-03-28 2007-03-06 Harris Corporation System and method for cumulant-based geolocation of cooperative and non-cooperative RF transmitters
US6931362B2 (en) * 2003-03-28 2005-08-16 Harris Corporation System and method for hybrid minimum mean squared error matrix-pencil separation weights for blind source separation
FR2853427B1 (en) * 2003-04-01 2005-06-17 Thales Sa METHOD OF BLINDLY IDENTIFYING MIXTURES OF SOURCES WITH HIGHER ORDERS
FR2853480B1 (en) * 2003-04-01 2005-06-17 METHOD AND DEVICE FOR AUTOMATICALLY IDENTIFYING A SUBDEFINED MIXTURE OF SOURCES IN THE FOURTH ORDER
US7672834B2 (en) * 2003-07-23 2010-03-02 Mitsubishi Electric Research Laboratories, Inc. Method and system for detecting and temporally relating components in non-stationary signals
JP4000095B2 (en) * 2003-07-30 2007-10-31 株式会社東芝 Speech recognition method, apparatus and program
WO2005040739A2 (en) * 2003-10-22 2005-05-06 Softmax, Inc. System and method for spectral analysis
FR2862173B1 (en) * 2003-11-07 2006-01-06 Thales Sa METHOD FOR BLINDLY DEMODULATING HIGHER ORDERS OF A LINEAR WAVEFORM TRANSMITTER
GB0326539D0 (en) 2003-11-14 2003-12-17 Qinetiq Ltd Dynamic blind signal separation
US20100265139A1 (en) * 2003-11-18 2010-10-21 Harris Corporation System and method for cumulant-based geolocation of cooperative and non-cooperative RF transmitters
US8204149B2 (en) 2003-12-17 2012-06-19 Qualcomm Incorporated Spatial spreading in a multi-antenna communication system
US7336746B2 (en) * 2004-12-09 2008-02-26 Qualcomm Incorporated Data transmission with spatial spreading in a MIMO communication system
US8169889B2 (en) 2004-02-18 2012-05-01 Qualcomm Incorporated Transmit diversity and spatial spreading for an OFDM-based multi-antenna communication system
US20050180312A1 (en) * 2004-02-18 2005-08-18 Walton J. R. Transmit diversity and spatial spreading for an OFDM-based multi-antenna communication system
KR100600313B1 (en) 2004-02-26 2006-07-14 남승현 Method and apparatus for frequency domain blind separation of multipath multichannel mixed signal
US20050238111A1 (en) * 2004-04-09 2005-10-27 Wallace Mark S Spatial processing with steering matrices for pseudo-random transmit steering in a multi-antenna communication system
US8923785B2 (en) * 2004-05-07 2014-12-30 Qualcomm Incorporated Continuous beamforming for a MIMO-OFDM system
US8285226B2 (en) * 2004-05-07 2012-10-09 Qualcomm Incorporated Steering diversity for an OFDM-based multi-antenna communication system
US7356075B2 (en) * 2004-05-24 2008-04-08 General Dynamics C4 Systems, Inc. Method and apparatus for signal separation
US7110463B2 (en) 2004-06-30 2006-09-19 Qualcomm, Incorporated Efficient computation of spatial filter matrices for steering transmit diversity in a MIMO communication system
US7978649B2 (en) 2004-07-15 2011-07-12 Qualcomm, Incorporated Unified MIMO transmission and reception
US7978778B2 (en) 2004-09-03 2011-07-12 Qualcomm, Incorporated Receiver structures for spatial spreading with space-time or space-frequency transmit diversity
US7860182B2 (en) * 2004-09-23 2010-12-28 Interdigital Technology Corporation Receiver hardware reduction for spatially independent signals and associated methods
US7414579B2 (en) * 2004-09-23 2008-08-19 Interdigital Technology Corporation Blind signal separation using correlated antenna elements
US7190308B2 (en) * 2004-09-23 2007-03-13 Interdigital Technology Corporation Blind signal separation using signal path selection
FR2887379B1 (en) * 2005-06-17 2007-08-31 Thales Sa METHOD FOR BLINDLY DEMODULATING UPPER ORDERS OF SEVERAL LINEAR WAVEFORM TRANSMITTERS
US7330801B2 (en) * 2005-07-29 2008-02-12 Interdigital Technology Corporation Signal separation using rank deficient matrices
JP2007034184A (en) * 2005-07-29 2007-02-08 Kobe Steel Ltd Device, program, and method for sound source separation
US7599714B2 (en) * 2005-09-30 2009-10-06 Alcatel-Lucent Usa Inc. Increasing the range of access point cells for a given throughput in a downlink of a wireless local area network
US7769118B2 (en) * 2006-02-10 2010-08-03 Interdigital Technolgy Corporation Method and apparatus for equalizing received signals based on independent component analysis
JP2009529699A (en) * 2006-03-01 2009-08-20 ソフトマックス,インコーポレイテッド System and method for generating separated signals
US8543070B2 (en) 2006-04-24 2013-09-24 Qualcomm Incorporated Reduced complexity beam-steered MIMO OFDM system
US8290089B2 (en) * 2006-05-22 2012-10-16 Qualcomm Incorporated Derivation and feedback of transmit steering matrix
JP2007324754A (en) * 2006-05-30 2007-12-13 Ntt Docomo Inc Signal receiving section detector
US7603245B2 (en) * 2006-06-20 2009-10-13 Southwest Research Institute Blind estimation of bandwidth and duration parameters of an incoming signal
WO2008002205A1 (en) * 2006-06-27 2008-01-03 Telefonaktiebolaget Lm Ericsson (Publ) A radio frequency emitter detection and location method and system
CN100433046C (en) * 2006-09-28 2008-11-12 上海大学 Image blind separation based on sparse change
JP5034469B2 (en) * 2006-12-08 2012-09-26 ソニー株式会社 Information processing apparatus, information processing method, and program
US7555064B2 (en) * 2007-08-27 2009-06-30 Harris Corporation System and method for estimating noise power level in a multi-signal communications channel
US7499515B1 (en) 2007-08-27 2009-03-03 Harris Corporation System and method for automated link quality measurement for adaptive modulation systems using noise level estimates
US7447276B1 (en) 2007-08-29 2008-11-04 Harris Corporation Receiver and method for blind adaptive thresholding using noise estimator
US20090060008A1 (en) * 2007-08-29 2009-03-05 Harris Corporation System and method for blind source separation of signals using noise estimator
US7884602B2 (en) * 2007-09-18 2011-02-08 Baker Hughes Incorporated Nuclear magnetic resonance evaluation using independent component analysis (ICA)-based blind source separation
US8195118B2 (en) 2008-07-15 2012-06-05 Linear Signal, Inc. Apparatus, system, and method for integrated phase shifting and amplitude control of phased array signals
JP5305841B2 (en) * 2008-11-04 2013-10-02 三菱電機株式会社 Signal processing device
CN102461035B (en) * 2009-05-28 2015-03-11 骁阳网络有限公司 Method and arrangement for blind demultiplexing a polarisation diversity multiplex signal
EP2476008B1 (en) 2009-09-10 2015-04-29 Rudjer Boskovic Institute Underdetermined blind extraction of components from mixtures in 1d and 2d nmr spectroscopy and mass spectrometry by means of combined sparse component analysis and detection of single component points
CA2693012C (en) * 2009-09-18 2015-06-23 Her Majesty The Queen In Right Of Canada As Represented By The Ministeof National Defence Signal detection in fading environments
US8340588B2 (en) * 2009-09-18 2012-12-25 Her Majesty The Queen In Right Of Canada As Represented By The Minister Of National Defence Handheld synthetic antenna array
US8666734B2 (en) * 2009-09-23 2014-03-04 University Of Maryland, College Park Systems and methods for multiple pitch tracking using a multidimensional function and strength values
US8872719B2 (en) 2009-11-09 2014-10-28 Linear Signal, Inc. Apparatus, system, and method for integrated modular phased array tile configuration
CN101729157B (en) * 2009-12-11 2016-02-17 南京航空航天大学 Method for separating vibration signal blind sources under a kind of strong noise environment
US8548067B2 (en) * 2010-01-29 2013-10-01 Goran Ivkovic Single sensor radio scene analysis for packet based radio signals using 2nd and 4th order statistics
CN101807301B (en) * 2010-03-17 2012-11-07 北京航空航天大学 High spectral image target detection method based on high order statistic
US8775341B1 (en) 2010-10-26 2014-07-08 Michael Lamport Commons Intelligent control with hierarchical stacked neural networks
US9015093B1 (en) 2010-10-26 2015-04-21 Michael Lamport Commons Intelligent control with hierarchical stacked neural networks
JP5675285B2 (en) * 2010-11-10 2015-02-25 富士通テン株式会社 Radar equipment
US9015007B2 (en) * 2010-12-13 2015-04-21 Southwest Research Institute Sensor array processor with multichannel reconstruction from random array sampling
US8878726B2 (en) 2011-03-16 2014-11-04 Exelis Inc. System and method for three-dimensional geolocation of emitters based on energy measurements
US8878725B2 (en) * 2011-05-19 2014-11-04 Exelis Inc. System and method for geolocation of multiple unknown radio frequency signal sources
CN102288285B (en) * 2011-05-24 2012-11-28 南京航空航天大学 Blind source separation method for single-channel vibration signals
US8615190B2 (en) 2011-05-31 2013-12-24 Exelis Inc. System and method for allocating jamming energy based on three-dimensional geolocation of emitters
US8723730B2 (en) 2011-07-27 2014-05-13 Exelis Inc. System and method for direction finding and geolocation of emitters based on line-of-bearing intersections
AU2014201436A1 (en) * 2013-03-22 2014-10-09 Cgg Services Sa System and method for interpolating seismic data
US10042037B2 (en) * 2014-02-20 2018-08-07 Nestwave Sas System and method for estimating time of arrival (TOA)
CN103973411B (en) * 2014-04-10 2017-12-22 沈阳理工大学 A kind of differential space-time coding signal blind checking method
CN104267265B (en) * 2014-09-30 2017-02-08 中国科学院新疆天文台 Evaluating system and method based on radio astronomical instrument device electromagnetic radiation
CN104796360B (en) * 2015-04-09 2018-04-17 西安交通大学 A kind of channel estimation methods based on ICA filtering methods
US20180284312A1 (en) * 2015-10-19 2018-10-04 Halliburton Energy Services, Inc. T1 distribution-based logging systems and methods using blind source separation independent component analysis
CN105609112A (en) * 2016-01-15 2016-05-25 苏州宾果智能科技有限公司 Sound source positioning method and apparatus and time delay estimation method and apparatus
GB2550581A (en) * 2016-05-23 2017-11-29 Vodafone Ip Licensing Ltd Dynamic cyclic prefix configuration
CN106226729B (en) * 2016-07-15 2018-08-31 西安电子科技大学 Relatively prime array direction of arrival angle method of estimation based on fourth-order cumulant
US10324168B2 (en) 2016-09-12 2019-06-18 The Boeing Company Systems and methods for spatial filtering using data with widely different error magnitudes
US10324167B2 (en) 2016-09-12 2019-06-18 The Boeing Company Systems and methods for adding functional grid elements to stochastic sparse tree grids for spatial filtering
US10429491B2 (en) * 2016-09-12 2019-10-01 The Boeing Company Systems and methods for pulse descriptor word generation using blind source separation
US10332530B2 (en) * 2017-01-27 2019-06-25 Google Llc Coding of a soundfield representation
JP6741864B2 (en) * 2017-04-13 2020-08-19 日本電信電話株式会社 Signal separation device and signal separation method
EP3596484B1 (en) * 2017-05-10 2022-06-01 Huawei Technologies Co., Ltd. Method and device for estimating a time of arrival of a radio signal
CN108366034B (en) * 2018-02-11 2020-02-14 哈尔滨工程大学 High-order single-channel digital-analog integrated communication method and system
CN108872928A (en) * 2018-03-19 2018-11-23 长江大学 The evaluation method of planar array angle based on four linear decompositions
CN108594194A (en) * 2018-03-19 2018-09-28 长江大学 Bistatic MIMO radar angle evaluation method based on four linear decompositions
CN109188473B (en) * 2018-07-26 2020-12-15 广东工业大学 Beidou satellite weak signal high-precision rapid capturing method based on blind separation technology
CN109192220B (en) * 2018-08-02 2022-11-11 福建师范大学福清分校 Mixed signal self-adaptive rapid separation method based on proportional-differential control
KR102192564B1 (en) * 2019-01-23 2020-12-17 국방과학연구소 Apparatus and method for controlling performance of receiver for sub-device in mimo cognitive radio systems
CN110236589B (en) * 2019-06-03 2022-04-29 苏州美糯爱医疗科技有限公司 Real-time heart-lung sound automatic separation method of electronic stethoscope
US11108457B2 (en) * 2019-12-05 2021-08-31 Bae Systems Information And Electronic Systems Integration Inc. Spatial energy rank detector and high-speed alarm
CN111027645B (en) * 2019-12-27 2023-05-26 北京工业大学 Microbial fermentation process fault monitoring method based on high-order statistic construction in sub-stage
CN111025424B (en) * 2020-01-06 2021-09-07 中国石油化工股份有限公司 Pre-stack time migration response operator quantification judgment method
CN111239718B (en) * 2020-01-17 2022-11-15 电子科技大学 Multi-base-station target positioning method based on single-satellite radiation source
CN112307927A (en) * 2020-10-26 2021-02-02 重庆邮电大学 BP network-based identification research for MPSK signals in non-cooperative communication
CN112881975B (en) * 2021-01-08 2023-09-08 电子科技大学 Single pulse and difference beam angle measurement method based on subarray feature matrix joint diagonalization
CN112885366A (en) * 2021-01-18 2021-06-01 宇龙计算机通信科技(深圳)有限公司 Active noise reduction method and device, storage medium and terminal
CN113376604A (en) * 2021-05-17 2021-09-10 西安电子科技大学 Airborne three-dimensional heterogeneous array clutter and interference suppression method based on blind source separation
CN113345456B (en) * 2021-05-31 2023-06-06 北京小米移动软件有限公司 Echo separation method, device and storage medium
CN114039824B (en) * 2021-11-09 2022-08-09 北京邮电大学 Data transmission method and system based on bit operation and layered enumeration spherical forming
CN114371145A (en) * 2022-03-21 2022-04-19 武汉工程大学 Detection method and device for milk oil mixed pigment, electronic equipment and storage medium
CN116866124A (en) * 2023-07-13 2023-10-10 中国人民解放军战略支援部队航天工程大学 Blind separation method based on baseband signal time structure

Citations (44)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3031527A (en) * 1959-09-08 1962-04-24 Manson Lab Inc Ultra stable frequency shift keying system
US3223925A (en) * 1962-01-29 1965-12-14 Ibm Digital data modulation device
US4166980A (en) * 1977-08-25 1979-09-04 Sanders Associates, Inc. Method and apparatus for signal recognition
US4368439A (en) * 1978-08-24 1983-01-11 Fujitsu Limited Frequency shift keying system
US4506225A (en) * 1981-12-28 1985-03-19 Barringer Research Limited Method for remote measurement of anomalous complex variations of a predetermined electrical parameter in a target zone
US4607305A (en) * 1983-02-08 1986-08-19 Ampex Corporation Monolithic multichannel multistack magnetic transducer assembly and method of manufacturing thereof
US4825396A (en) * 1986-02-14 1989-04-25 Siemens Aktiengesellschaft Digital circuit for sampling rate variation and signal filtering and method for constructing the circuit
US4965732A (en) * 1985-11-06 1990-10-23 The Board Of Trustees Of The Leland Stanford Junior University Methods and arrangements for signal reception and parameter estimation
US5283813A (en) * 1991-02-24 1994-02-01 Ramat University Authority For Applied Research & Industrial Development Ltd. Methods and apparatus particularly useful for blind deconvolution
US5381450A (en) * 1993-08-20 1995-01-10 Hitachi America, Ltd. Technique for automatically detecting the constellation size of a quadrature amplitude modulated (QAM) signal
US5459668A (en) * 1993-06-04 1995-10-17 University Of Southern California Method and apparatus for signal analysis employing a virtual cross-correlation computer
US5493612A (en) * 1962-03-27 1996-02-20 The United States Of America As Represented By The Secretary Of The Navy Secure communication keying system
US5539832A (en) * 1992-04-10 1996-07-23 Ramot University Authority For Applied Research & Industrial Development Ltd. Multi-channel signal separation using cross-polyspectra
US5640419A (en) * 1994-07-19 1997-06-17 Grumman Aerospace Corporation Covert communication system
US5651030A (en) * 1991-05-06 1997-07-22 Motorola, Inc. Receiver with signal classifier
US5706402A (en) * 1994-11-29 1998-01-06 The Salk Institute For Biological Studies Blind signal processing system employing information maximization to recover unknown signals through unsupervised minimization of output redundancy
US5848160A (en) * 1996-02-20 1998-12-08 Raytheon Company Digital synthesized wideband noise-like waveform
US5883922A (en) * 1993-05-12 1999-03-16 Olivetti Telemedia S.P.A. FSK communication device for digital information
US5893031A (en) * 1996-06-27 1999-04-06 Cellular Technical Services Company, Inc. System and method for collection of transmission characteristics
US5909646A (en) * 1995-02-22 1999-06-01 U.S. Philips Corporation System for estimating signals received in the form of mixed signals
US5959966A (en) * 1997-06-02 1999-09-28 Motorola, Inc. Methods and apparatus for blind separation of radio signals
US5996406A (en) * 1997-05-20 1999-12-07 Motorola, Inc. Advanced signal process for a material storage measuring device
US6021334A (en) * 1996-11-07 2000-02-01 France Telecom Method for transmission by a base station equipped with a multi-element antenna to a mobile
US6177906B1 (en) * 1999-04-01 2001-01-23 Arraycomm, Inc. Multimode iterative adaptive smart antenna processing method and apparatus
US6185309B1 (en) * 1997-07-11 2001-02-06 The Regents Of The University Of California Method and apparatus for blind separation of mixed and convolved sources
US6208295B1 (en) * 1995-06-02 2001-03-27 Trw Inc. Method for processing radio signals that are subject to unwanted change during propagation
US6252962B1 (en) * 1994-04-22 2001-06-26 Raytheon Company Featureless covert communication system
US6310704B1 (en) * 1995-06-02 2001-10-30 Trw Inc. Communication apparatus for transmitting and receiving signals over a fiber-optic waveguide using different frequency bands of light
US6426977B1 (en) * 1999-06-04 2002-07-30 Atlantic Aerospace Electronics Corporation System and method for applying and removing Gaussian covering functions
US6449560B1 (en) * 2000-04-19 2002-09-10 Schlumberger Technology Corporation Sonic well logging with multiwave processing utilizing a reduced propagator matrix
US20020153891A1 (en) * 1999-07-06 2002-10-24 Smith John Alec Sydney Methods of and apparatus for analysing a signal
US6522691B1 (en) * 1999-06-28 2003-02-18 Nortel Networks Limited Method of and apparatus for determining the capacity of a telecommunications system
US6535666B1 (en) * 1995-06-02 2003-03-18 Trw Inc. Method and apparatus for separating signals transmitted over a waveguide
US20030204378A1 (en) * 2001-06-26 2003-10-30 Cai Khiem V. System and method for forming a beam and creating nulls with an adaptive array antenna using orthogonal eigen-weighting
US6658234B1 (en) * 1995-06-02 2003-12-02 Northrop Grumman Corporation Method for extending the effective dynamic range of a radio receiver system
US6681054B1 (en) * 2000-05-01 2004-01-20 Eastman Kodak Company Noise reduction method utilizing probabilistic weighting, apparatus, and program for digital image processing
US20040012525A1 (en) * 2000-05-12 2004-01-22 Matsushita Electric Industrial Co., Ltd Direction of arrival estimator and direction of arrival estimation method
US6711528B2 (en) * 2002-04-22 2004-03-23 Harris Corporation Blind source separation utilizing a spatial fourth order cumulant matrix pencil
US20040259565A1 (en) * 2003-06-20 2004-12-23 Nortel Networks Limited Method and device for selecting parameters for a cellular radio communication network based on occurrence frequencies
US6845164B2 (en) * 1999-03-08 2005-01-18 Telefonaktiebolaget Lm Ericsson (Publ) Method and device for separating a mixture of source signals
US6941123B2 (en) * 2001-12-08 2005-09-06 Sas Technologies Co., Ltd. Signal processing method and apparatus for obtaining diversity gain in wireless communication system equipping antenna array
US20050212703A1 (en) * 2003-03-27 2005-09-29 Fiore Paul D Direction-finding for multiple cochannel sources
US7068977B1 (en) * 2002-10-11 2006-06-27 Navini Networks, Inc. Method and system for interference assessment and reduction in a wireless communication system
US7068715B2 (en) * 2001-10-10 2006-06-27 Genral Electric Company Ultra-wideband communications system and method using a delay hopped, continuous noise transmitted reference

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6167417A (en) * 1998-04-08 2000-12-26 Sarnoff Corporation Convolutive blind source separation using a multiple decorrelation method
EP1017253B1 (en) * 1998-12-30 2012-10-31 Siemens Corporation Blind source separation for hearing aids

Patent Citations (44)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3031527A (en) * 1959-09-08 1962-04-24 Manson Lab Inc Ultra stable frequency shift keying system
US3223925A (en) * 1962-01-29 1965-12-14 Ibm Digital data modulation device
US5493612A (en) * 1962-03-27 1996-02-20 The United States Of America As Represented By The Secretary Of The Navy Secure communication keying system
US4166980A (en) * 1977-08-25 1979-09-04 Sanders Associates, Inc. Method and apparatus for signal recognition
US4368439A (en) * 1978-08-24 1983-01-11 Fujitsu Limited Frequency shift keying system
US4506225A (en) * 1981-12-28 1985-03-19 Barringer Research Limited Method for remote measurement of anomalous complex variations of a predetermined electrical parameter in a target zone
US4607305A (en) * 1983-02-08 1986-08-19 Ampex Corporation Monolithic multichannel multistack magnetic transducer assembly and method of manufacturing thereof
US4965732A (en) * 1985-11-06 1990-10-23 The Board Of Trustees Of The Leland Stanford Junior University Methods and arrangements for signal reception and parameter estimation
US4825396A (en) * 1986-02-14 1989-04-25 Siemens Aktiengesellschaft Digital circuit for sampling rate variation and signal filtering and method for constructing the circuit
US5283813A (en) * 1991-02-24 1994-02-01 Ramat University Authority For Applied Research & Industrial Development Ltd. Methods and apparatus particularly useful for blind deconvolution
US5651030A (en) * 1991-05-06 1997-07-22 Motorola, Inc. Receiver with signal classifier
US5539832A (en) * 1992-04-10 1996-07-23 Ramot University Authority For Applied Research & Industrial Development Ltd. Multi-channel signal separation using cross-polyspectra
US5883922A (en) * 1993-05-12 1999-03-16 Olivetti Telemedia S.P.A. FSK communication device for digital information
US5459668A (en) * 1993-06-04 1995-10-17 University Of Southern California Method and apparatus for signal analysis employing a virtual cross-correlation computer
US5381450A (en) * 1993-08-20 1995-01-10 Hitachi America, Ltd. Technique for automatically detecting the constellation size of a quadrature amplitude modulated (QAM) signal
US6252962B1 (en) * 1994-04-22 2001-06-26 Raytheon Company Featureless covert communication system
US5640419A (en) * 1994-07-19 1997-06-17 Grumman Aerospace Corporation Covert communication system
US5706402A (en) * 1994-11-29 1998-01-06 The Salk Institute For Biological Studies Blind signal processing system employing information maximization to recover unknown signals through unsupervised minimization of output redundancy
US5909646A (en) * 1995-02-22 1999-06-01 U.S. Philips Corporation System for estimating signals received in the form of mixed signals
US6535666B1 (en) * 1995-06-02 2003-03-18 Trw Inc. Method and apparatus for separating signals transmitted over a waveguide
US6658234B1 (en) * 1995-06-02 2003-12-02 Northrop Grumman Corporation Method for extending the effective dynamic range of a radio receiver system
US6208295B1 (en) * 1995-06-02 2001-03-27 Trw Inc. Method for processing radio signals that are subject to unwanted change during propagation
US6310704B1 (en) * 1995-06-02 2001-10-30 Trw Inc. Communication apparatus for transmitting and receiving signals over a fiber-optic waveguide using different frequency bands of light
US5848160A (en) * 1996-02-20 1998-12-08 Raytheon Company Digital synthesized wideband noise-like waveform
US5893031A (en) * 1996-06-27 1999-04-06 Cellular Technical Services Company, Inc. System and method for collection of transmission characteristics
US6021334A (en) * 1996-11-07 2000-02-01 France Telecom Method for transmission by a base station equipped with a multi-element antenna to a mobile
US5996406A (en) * 1997-05-20 1999-12-07 Motorola, Inc. Advanced signal process for a material storage measuring device
US5959966A (en) * 1997-06-02 1999-09-28 Motorola, Inc. Methods and apparatus for blind separation of radio signals
US6185309B1 (en) * 1997-07-11 2001-02-06 The Regents Of The University Of California Method and apparatus for blind separation of mixed and convolved sources
US6845164B2 (en) * 1999-03-08 2005-01-18 Telefonaktiebolaget Lm Ericsson (Publ) Method and device for separating a mixture of source signals
US6177906B1 (en) * 1999-04-01 2001-01-23 Arraycomm, Inc. Multimode iterative adaptive smart antenna processing method and apparatus
US6426977B1 (en) * 1999-06-04 2002-07-30 Atlantic Aerospace Electronics Corporation System and method for applying and removing Gaussian covering functions
US6522691B1 (en) * 1999-06-28 2003-02-18 Nortel Networks Limited Method of and apparatus for determining the capacity of a telecommunications system
US20020153891A1 (en) * 1999-07-06 2002-10-24 Smith John Alec Sydney Methods of and apparatus for analysing a signal
US6449560B1 (en) * 2000-04-19 2002-09-10 Schlumberger Technology Corporation Sonic well logging with multiwave processing utilizing a reduced propagator matrix
US6681054B1 (en) * 2000-05-01 2004-01-20 Eastman Kodak Company Noise reduction method utilizing probabilistic weighting, apparatus, and program for digital image processing
US20040012525A1 (en) * 2000-05-12 2004-01-22 Matsushita Electric Industrial Co., Ltd Direction of arrival estimator and direction of arrival estimation method
US20030204378A1 (en) * 2001-06-26 2003-10-30 Cai Khiem V. System and method for forming a beam and creating nulls with an adaptive array antenna using orthogonal eigen-weighting
US7068715B2 (en) * 2001-10-10 2006-06-27 Genral Electric Company Ultra-wideband communications system and method using a delay hopped, continuous noise transmitted reference
US6941123B2 (en) * 2001-12-08 2005-09-06 Sas Technologies Co., Ltd. Signal processing method and apparatus for obtaining diversity gain in wireless communication system equipping antenna array
US6711528B2 (en) * 2002-04-22 2004-03-23 Harris Corporation Blind source separation utilizing a spatial fourth order cumulant matrix pencil
US7068977B1 (en) * 2002-10-11 2006-06-27 Navini Networks, Inc. Method and system for interference assessment and reduction in a wireless communication system
US20050212703A1 (en) * 2003-03-27 2005-09-29 Fiore Paul D Direction-finding for multiple cochannel sources
US20040259565A1 (en) * 2003-06-20 2004-12-23 Nortel Networks Limited Method and device for selecting parameters for a cellular radio communication network based on occurrence frequencies

Cited By (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050118963A1 (en) * 2003-12-02 2005-06-02 Mao-Ching Chiu Method and apparatus for I/Q imbalance calibration of a transmitter system
US7139536B2 (en) * 2003-12-02 2006-11-21 Mediatek Inc. Method and apparatus for I/Q imbalance calibration of a transmitter system
US20080001821A1 (en) * 2004-09-14 2008-01-03 Akira Tanaka Signal Arrival Direction Deducing Device, Signal Arrival Direction Deducing Method, and Signal Direction Deducing Program
US7436358B2 (en) * 2004-09-14 2008-10-14 National University Corporation Hokkaido University Signal arrival direction deducing device, signal arrival direction deducing method, and signal direction deducing program
US8811502B2 (en) 2005-06-22 2014-08-19 Eices Research, Inc. Systems and/or methods of wireless communications
US9641202B2 (en) 2005-06-22 2017-05-02 Odyssey Wireless, Inc. Systems/methods of carrier aggregation
US7876845B2 (en) 2005-06-22 2011-01-25 Eices Research, Inc. Wireless communications systems and/or methods providing low interference, high privacy and/or cognitive flexibility
US7881393B2 (en) 2005-06-22 2011-02-01 Eices Research, Inc. Waveforms comprising a plurality of elements and transmission thereof
US20110206086A1 (en) * 2005-06-22 2011-08-25 Eices Research, Inc. Systems/methods of sequential modulation of a single carrier frequency by a plurality of elements of a waveform
US8050337B2 (en) * 2005-06-22 2011-11-01 Eices Research, Inc. Systems, methods, devices, and/or computer program products for providing communications devoid of cyclostationary features
US8199837B2 (en) 2005-06-22 2012-06-12 Eices Research, Inc. Systems/methods of sequential modulation of a single carrier frequency by a plurality of elements of a waveform
USRE47633E1 (en) 2005-06-22 2019-10-01 Odyssey Wireless Inc. Systems/methods of conducting a financial transaction using a smartphone
US8537910B2 (en) 2005-06-22 2013-09-17 Eices Research, Inc. Private, covert and/or cognitive communications systems and/or methods based upon pseudo-randomly generated communications alphabets
US9705535B2 (en) 2005-06-22 2017-07-11 Odyssey Wireless, Inc. Systems/methods of carrier aggregation
US8576940B2 (en) 2005-06-22 2013-11-05 Eices Research, Inc. Systems/methods of adaptively varying a bandwidth and/or frequency content of communications
US8660169B1 (en) 2005-06-22 2014-02-25 Eices Research, Inc. Systems/methods of adaptively varying a bandwidth and/or frequency content of communications
US8670493B2 (en) 2005-06-22 2014-03-11 Eices Research, Inc. Systems and/or methods of increased privacy wireless communications
US20090268783A1 (en) * 2005-06-22 2009-10-29 Eices Research, Inc. Wireless communications systems and/or methods providing low interference, high privacy and/or cognitive flexibility
US20090003418A1 (en) * 2005-06-22 2009-01-01 Karabinis Peter D Systems, Methods, Devices And/Or Computer Program Products For Providing Communications Deviod Of Cyclostationary Features
US8855230B1 (en) 2005-06-22 2014-10-07 Eices Research, Inc. Systems/methods of transmitting information via baseband waveforms comprising frequency content agility and an orthogonality therebetween
US8879606B2 (en) 2005-06-22 2014-11-04 Eices Research, Inc. Systems/methods of transmitting information via baseband waveforms comprising agility in frequency content and an orthogonality therebetween
US8891645B2 (en) 2005-06-22 2014-11-18 Eices Research, Inc. Systems/methods of carrier aggregation providing increased capacity communications
US20140348211A1 (en) * 2005-06-22 2014-11-27 Eices Research, Inc. Systems/methods of adaptively varying a spectral content of communications
US9392451B2 (en) 2005-06-22 2016-07-12 Odyssey Wireless, Inc. Systems/methods of conducting a financial transaction using a smartphone
US9124381B2 (en) 2005-06-22 2015-09-01 Odyssey Wireless, Inc. Systems/methods of carrier aggregation
US9185553B2 (en) 2005-06-22 2015-11-10 Odyssey Wireless, Inc. Systems/methods of preferential communications
US9332429B2 (en) * 2005-06-22 2016-05-03 Odyssey Wireless, Inc. Systems/methods of adaptively varying a spectral content of communications
US9374746B1 (en) 2008-07-07 2016-06-21 Odyssey Wireless, Inc. Systems/methods of spatial multiplexing
US8537916B2 (en) 2010-03-29 2013-09-17 Eices Research, Inc. Increased capacity communications for OFDM-based wireless communications systems/methods/devices
US9806790B2 (en) 2010-03-29 2017-10-31 Odyssey Wireless, Inc. Systems/methods of spectrally efficient communications
US20130096699A1 (en) * 2010-06-21 2013-04-18 Optimized Systems And Solutions Limited Asset health monitoring
US9054870B2 (en) 2012-10-22 2015-06-09 Donatello Apelusion Gassi Information security based on eigendecomposition
WO2014066279A1 (en) * 2012-10-22 2014-05-01 Gassi Donatello Apelusion Information security based on eigendecomposition
CN108702219A (en) * 2016-03-08 2018-10-23 骁阳网络有限公司 Reduce the adaptive constellation that phase distortion influences
US20180007045A1 (en) * 2016-06-30 2018-01-04 Mehdi Arashmid Akhavain Mohammadi Secure coding and modulation for optical transport
CN109547371A (en) * 2019-01-17 2019-03-29 西安电子科技大学 Signal processing method based on quasi- high power spectrum
CN114095045A (en) * 2022-01-24 2022-02-25 成都航天通信设备有限责任公司 Communication method, device and system based on constellation pseudo-shape modulation

Also Published As

Publication number Publication date
CN1656485A (en) 2005-08-17
ATE487187T1 (en) 2010-11-15
US6711528B2 (en) 2004-03-23
EP1500007A4 (en) 2009-09-02
DE60334790D1 (en) 2010-12-16
JP4271041B2 (en) 2009-06-03
US20030204380A1 (en) 2003-10-30
EP1500007A1 (en) 2005-01-26
CN100386764C (en) 2008-05-07
JP2005523475A (en) 2005-08-04
EP1500007B1 (en) 2010-11-03
WO2003090127A1 (en) 2003-10-30
AU2003223605A1 (en) 2003-11-03

Similar Documents

Publication Publication Date Title
US20030228017A1 (en) Method and system for waveform independent covert communications
US7187326B2 (en) System and method for cumulant-based geolocation of cooperative and non-cooperative RF transmitters
US6993440B2 (en) System and method for waveform classification and characterization using multidimensional higher-order statistics
US8699971B1 (en) System for specific emitter identification
US7606524B1 (en) Integrated monitoring and communications receiver architecture
US8755469B1 (en) Method of spectrum mapping and exploitation using distributed sensors
US5422913A (en) High frequency multichannel diversity differential phase shift (DPSK) communications system
US5265121A (en) Spread spectrum coherent processor
US9057772B1 (en) Referenced-based radio frequency feature estimation
Zeng et al. Robustness of the cyclostationary detection to cyclic frequency mismatch
KR20110004443A (en) Passive radar for presence and motion detection
JPH07105777B2 (en) Method and apparatus for improving anti-jam performance with conversion gain
US8917802B1 (en) Modulation scheme determination through use of multiple sensors
US20100265139A1 (en) System and method for cumulant-based geolocation of cooperative and non-cooperative RF transmitters
US6993460B2 (en) Method and system for tracking eigenvalues of matrix pencils for signal enumeration
US6263031B1 (en) Method and apparatus for signal burst classification
US6973142B2 (en) Timing synchronization for M-DPSK channels
CN109085613B (en) Satellite deception jamming identification method and device based on constellation trajectory
US8488578B1 (en) Identifying a CDMA scrambling code
Pearce et al. Signal discrimination and exploitation of ads-b transmission
Parmaksız et al. A review of recent developments on secure authentication using RF fingerprints techniques
Guo et al. A Mode 5 signal detection method based on phase and amplitude correlation
Venturino et al. Radar-enabled ambient backscatter communications
Xiao et al. An investigation of non-data-aided SNR estimation techniques for analog modulation signals
US10560254B1 (en) Picosecond LPD time alignment

Legal Events

Date Code Title Description
AS Assignment

Owner name: HARRIS CORPORATION, FLORIDA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:BEADLE, EDWARD RAY;DISHMAN, JOHN FITZGERALD;ANDERSON, RICHARD HUGH;REEL/FRAME:016764/0976;SIGNING DATES FROM 20051107 TO 20051110

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO PAY ISSUE FEE