Removal of Directwave High Frequency Signal for Ionospheric Sounder Return Processing
20200018841 ยท 2020-01-16
Inventors
- James H. Nixon (Helotes, TX)
- Brock A. Martin (Helotes, TX)
- John J. Signorotti (Fair Oaks Ranch, TX)
- Robert G. Fenske, Jr. (San Antonio, TX)
Cpc classification
G01S19/07
PHYSICS
G01S13/0218
PHYSICS
Y02A90/10
GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
International classification
Abstract
A method of operating a vertical incidence sounder. The sounder transmits a random modulated sounder signal, and receives a return signal, comprising a directwave component and a reflected component. A first cross-ambiguity process between the transmitted signal and the return signal is used to time and frequency align the return signal with the transmitted signal. The direct wave signal is removed, thereby producing a processed return signal. A second cross-ambiguity process between the transmitted signal and the processed return signal is performed to determine a time offset, which is used to estimate an ionospheric height.
Claims
1. A method of operating a vertical incidence sounder, comprising: generating a random modulated sounder waveform; transmitting a random modulated signal based on the sounder waveform; receiving a return signal, comprising a directwave component and a reflected component; performing a first cross-ambiguity process between the transmitted signal and the return signal; using a first peak value from the results of the first cross-ambiguity process to time and frequency align the return signal with the transmitted signal; removing the direct wave signal, thereby producing a processed return signal; performing a second cross-ambiguity process between the transmitted signal and the processed return signal; and using a second peak value from the result of the cross-ambiguity process to determine a time offset; and using the time offset to determine an ionospheric height.
2. The method of claim 1, further comprising the steps of using the second peak value to determine a frequency offset, and using the frequency offset to determine a Doppler shift.
3. The method of claim 1, wherein the step of removing the direct wave signal is performed with a nulling process.
4. The method of claim 1, wherein the step of using a second peak value to determine a time offset is performed by determining peaks above a predetermined time delay threshold.
5. The method of claim 4, wherein the time delays correspond to 100 to 700 kilometer heights.
6. A vertical incidence sounder, comprising: two dipole antennas; a transceiver for transmitting and receiving a transmit signal and a return signal; a direct wave removal process for performing the following tasks: generating a random modulated sounder waveform and delivering a random modulated signal based on the sounder waveform to the transceiver; receiving the return signal, comprising a directwave component and a reflected component; performing a first cross-ambiguity process between the transmitted signal and the return signal; using a first peak value from the results of the first cross-ambiguity process to time and frequency align the return signal with the transmitted signal; removing the direct wave signal, thereby producing a processed return signal; performing a second cross-ambiguity process between the transmitted signal and the processed return signal; and using a second peak value from the result of the cross-ambiguity process to determine a time offset; and using the time offset to determine an ionospheric height.
7. The sounder of claim 6, wherein the direct wave removal process further uses the second peak value to determine a frequency offset, and using the frequency offset to determine a Doppler shift.
8. The sounder of claim 6, wherein removing the direct wave signal is performed with a nulling process.
9. The sounder of claim 6, wherein using a second peak value to determine a time offset is performed by determining peaks above a predetermined time delay threshold.
10. The sounder of claim 9, wherein the time delays correspond to 100 to 700 kilometer heights.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0008] A more complete understanding of the present embodiments and advantages thereof may be acquired by referring to the following description taken in conjunction with the accompanying drawings, in which like reference numbers indicate like features, and wherein:
[0009]
[0010]
[0011]
[0012]
[0013]
DETAILED DESCRIPTION OF THE INVENTION
[0014] The following description is directed to a vertical incidence ionospheric sounder that has processing for removing the unwanted direct wave component from the return signal. As explained below, the time of transmission of the direct wave is intrinsically measured. A novel signal nulling technique removes the direct wave component without the need for T/R switching or gating. This method simplifies the sounder hardware, and eliminates the processing gain penalty of conventional sounder designs.
[0015] The sounder may be operated at very low power levels (on the order of 10 milliwatts) as compared to the 10's or 100's of watts used by conventional sounders. It can be implemented using standard, off-the-shelf, software defined radio (SDR) HF transceiver hardware. Software-defined radio (SDR) is a radio communication system where components that have been traditionally implemented in hardware (e.g. mixers, filters, amplifiers, modulators/demodulators, detectors, etc.) are instead implemented by means of software on a personal computer or embedded processing system.
[0016]
[0017]
[0018]
[0019] Transceiver 207 sends and receives sounding signals, as described below. Processing system 209 performs signal processing tasks, including removal of the directwave from the return signal, as described below.
[0020] Unlike conventional vertical incidence sounders, processing system 209 that does not use T/R switching. As stated in the Background, vertical incidence sounders conventionally avoid overloading their receiver with the unwanted direct signal by using T/R switching.
[0021] Processing system 209 is assumed to have appropriate processing hardware and software for performing the tasks described herein. Processing system 209 may be part of a more comprehensive sounder processing system for performing additional tasks related to ionospheric sounding.
[0022]
[0023] Step 201 is generating a random modulated digital I/Q sounder waveform. Random modulation is a known modulation technique, which creates a new signal from two other signals by the process of quadrature amplitude modulation. In I/Q modulation, I is the in-phase component of the waveform and Q represents the quadrature component.
[0024] In Step 201, the generated waveform is a fixed frequency transmission that is transmitted for a short time. Then, another frequency is selected and transmitted, and so on. This signal is used herein as a reference signal during processing, and is thus referred to below as a reference signal.
[0025] Step 202 is transmitting the random modulated signal generated in Step 201.
[0026] Step 203 is receiving the return signal. Prior to processing, the return signal comprises both the directwave component and the ionospheric reflected signal component.
[0027] Step 204 is digitizing both the transmitted (reference) signal and return signal, and performing a cross-ambiguity function between them. The digitized data from these signals, respectively, may be referred to as the reference data and the collected data for processing purposes.
[0028] Step 205 is determining the peak value of the cross-ambiguity function. This peak value is used to time and frequency align the return waveform received in Step 203 with the reference waveform generated in Step 201. For frequency alignment, the collected data is mixed by the frequency of the peak to remove the frequency offset relative to the reference data. The frequency-adjusted data is aligned in time. The result is the best estimate of the time offset of the directwave.
[0029] Step 206 is applying a nulling process to remove the directwave signal from the return waveform as aligned in Step 205. The directwave signal is canceled from the collected data across time and frequency dimensions.
[0030] Step 207 is performing a cross-ambiguity process between the waveform resulting from Step 206 (the return signal after nulling the directwave component) and the reference waveform. To remove unwanted correlation artifacts, the negative cross-ambiguity plane is subtracted from the positive plane.
[0031] Step 208 is finding peaks in data resulting from Step 207 that are above desired thresholds and are within expected time delays. The expected time delays are those that correspond to 100-700 kilometer heights. These results correspond to ionospheric return heights and Doppler shift.
[0032]
[0033]
[0034] Thus, the vector A is a time-varying vector of measured and reference waveforms, and when all time samples are collected, the result is a matrix A.
[0035] The matrix R is a correlation (covariance) matrix between the measured waveform and reference waveforms.
[0036] The vector a is a constraint vector, which forces nulling of all reference waveform energy matching the u.sub.mn references.
[0037] The vector W is a normalized weight vector, which is used to optimally combine the measured and reference waveforms to remove energy matching the reference waveforms.
[0038] The values {dot over (v)}=W.sup.HA represent the subsequent signal with all localized reference waveform energy removed. These values contain the weak returned reference signal, as it does not match the filtered versions.