Terrestrial and space navigation, rendezvous, and docking efforts require precise navigational data in order to derive accurate knowledge of vehicle relative velocity, range, and altitude. For instance, spacecraft rely on real-time navigational data to control descent and ensure a soft landing at a designated landing site. Radio and/or light wave-based navigation systems may also be used to such ends, with range-to-ground surface measurements performed along multiple lines-of-sight (LOS). The multiple LOS measurements are used by an onboard navigation system to determine vehicle attitude, altitude, vector velocity, and other control parameters. Furthermore, the derived velocity may be used to drive angles of attack and side-slip control actions for real-time guidance and control.
Light Detection and Ranging, better known as lidar, is a particular type of remote sensing technology that functions by illuminating a target-of-interest, typically using continuous wave (CW) or quasi-CW light having an instantaneous frequency that varies linearly with time. A round-trip time-of-flight of each reflected light pulse is determined using a sensor to estimate a range to the target. Lidar systems obtain high-resolution range and relative velocity information from such light beams. In a homodyne receiver configuration common to lidar systems, a portion of the transmitted laser beam serves as a local oscillator for an optical receiver. Mixing of the local oscillator field with a time-delayed received field yields an intermediate frequency (IF) signal that is directly related to the target range, with the IF being fixed for each returned pulse.
A particular embodiment of the lidar system described herein transmits modulated three-section waveforms in different fixed directions to determine relative velocity and range to ground along three laser LOS. The transmitted waveforms are “three-section” in the sense of having separate Doppler, up-chirp/up-ramp, and down-chirp/down-ramp sections. The three LOS measurements are combined to determine three components of a velocity vector, and to accurately measure altitude and attitude of a spacecraft, aircraft, or another host vehicle relative to a ground surface.
A frequency-modulated continuous wave (FMCW) waveform is one in which the laser frequency is modulated linearly with time rather than its amplitude. That is, rather than pulses of light, a frequency-variant continuous wave is detected by the homodyne receiver setup. Linear FMCW waveforms provide the above-noted distance to a target, as well as relative velocity between the sensor and the target. Underlying measurements are obtained in the frequency domain by the homodyne receiver. However, phase information is lost during homodyne detection, and thus conventional approaches to implementing FMCW techniques in a lidar system can produce inaccurate measurements.
An improved frequency-modulated continuous wave (FMCW) Doppler Lidar system is disclosed herein, along with associated methods for using such a system. The Doppler lidar system uses a homodyne receiver to detect a three-section return waveform, i.e., one having a Doppler intermediate frequency (IF) section and separate up-ramp and down-ramp IF sections. Phase information is lost during detection of the three IFs described herein, as noted generally above. Sign ambiguities therefore occur as a result of the lost frequency information. Incorrect sign assignments in turn can lead to inaccurate conversion of the IFs to corresponding engineering units, and thus to navigational errors. However, the nature of the three-section return waveform available at the homodyne receiver provides a way to recover phase information, with the present teachings situationally exploiting qualities of the three-section waveform in a sign estimation and assignment process to improve navigational accuracy as set forth herein.
As will be appreciated by those of ordinary skill in the art, a three-section return waveform within the context of the FMCW lidar includes the above-noted linear up-ramp or “up-chirp” IF section, the linear down-ramp or “down-chirp” IF section, and the constant frequency/Doppler IF section. The present approach uses timing knowledge of the three-section return waveform on each channel of the homodyne receiver to mitigate sign ambiguities, i.e., the particular condition that exists whenever one or more of the three IF sections is not detected due to signal loss, drop out, interference, or other reasons. The method may be encoded as a computer-executable algorithm and used by a control circuit of the Doppler lidar system to estimate and assign a positive or negative sign to each IF section detected by the homodyne receiver, with range and/or velocity determined using the estimated sign. Operation of a host spacecraft, aircraft, or other vehicle may be controlled using the estimated sign and determining correct LOS range and velocity for each beam.
In a particular embodiment, the FMCW Doppler lidar system for use with a host vehicle includes a laser system, a frequency modulator, lenses, a homodyne receiver, and a control circuit/signal processing unit (SPU). The laser system is configured to generate first, second, and third laser beams and first, second, and third local oscillator beams using a seed laser and a beam splitter. The frequency modulator generates a modulated three-section waveform having a Doppler section, an up-ramp section, and a down-ramp section as described above. The lenses transmit the three-section waveform toward a target-of-interest, and receive respective first, second, and third return beams each reflected from the target-of-interest.
The homodyne receiver in this embodiment receives and mixes each of the return beams with a respective one of the local oscillator beams to detect a return waveform having a Doppler IF section, an up-ramp IF section, and a down-ramp IF section. The SPU is coupled to the homodyne receiver and configured to calculate or otherwise determine a respective magnitude and frequency of, and estimate a sign of, each of the up-ramp IF section, the down-ramp IF section, and/or the Doppler IF section, to calculate a range to the target-of-interest and/or velocity of the host vehicle relative to the target-of-interest using the respective magnitude, frequency, and estimated sign, and to thereafter directly or indirectly control an operation of the host vehicle using the range and/or the velocity.
A method is also used for using a FMCW Doppler lidar system aboard a host vehicle. The method includes generating, via a laser system having a seed laser and a beam splitter, first, second, and third laser beams and first, second, and third local oscillator beams. The method also includes generating a modulated three-section waveform via a frequency modulator, the three-section waveform having a Doppler section, an up-ramp section, and a down-ramp section. The method additionally includes transmitting, using first, second, and third lenses, the three-section waveform toward a target-of-interest external to the host vehicle, and receiving, using the lenses, respective first, second, and third return beams reflected from the target-of-interest.
The method in this example embodiment also includes receiving and mixing, via first, second, and third photoreceptors of a homodyne receiver, each of the return beams and a respective one of the local oscillator beams to detect a return waveform having a Doppler IF section, an up-ramp IF section, and a down-ramp IF section. The SPU determines a respective magnitude, frequency, and sign of the up-ramp IF section, the down-ramp IF section, and/or the Doppler IF section, with the sign being estimated as set forth herein, and calculates a range to the target-of-interest and/or velocity of the host vehicle relative to the target-of-interest using the respective magnitude, frequency, and sign. Thereafter, the method includes controlling an operation of the host vehicle using the range and/or the velocity.
These and other features, advantages, and objects of the present disclosure will be further understood and appreciated by those skilled in the art by reference to the following specification, claims, and appended drawings. The present disclosure is susceptible to various modifications and alternative forms, and some representative embodiments have been shown by way of example in the drawings and will be described in detail herein. It should be understood, however, that the novel aspects of this disclosure are not limited to the particular forms illustrated in the appended drawings. Rather, the disclosure is to cover all modifications, equivalents, combinations, subcombinations, permutations, groupings, and alternatives falling within the scope and spirit of the disclosure.
The present disclosure is susceptible to various modifications and alternative forms, and some representative embodiments have been shown by way of example in the drawings and will be described in detail herein. It should be understood, however, that the novel aspects of this disclosure are not limited to the particular forms illustrated in the appended drawings. Rather, the disclosure is to cover all modifications, equivalents, combinations, subcombinations, permutations, groupings, and alternatives falling within the scope and spirit of the disclosure.
For purposes of description herein, the terms “upper,” “lower,” “right,” “left,” “rear,” “front,” “vertical,” “horizontal,” and derivatives thereof shall relate to orientation shown in
The SPU 50 executes the method 100 to estimate the sign of one or more of the IF sections in a case-specific manner, with the case being dependent upon which of the IF sections are received and thus available to the SPU 50 at a given measurement instance. Such sign estimates, along with the respective calculated frequency and magnitude, are then used by the SPU 50 to derive range and/or velocity of the host vehicle 55 relative to a ground surface or target-of-interest, doing so without the use of external state estimators. As will be appreciated, the Doppler lidar system 10 depicted in
The exemplary Doppler lidar system 10 of
As part of such a laser system 11, a frequency modulator 14 is coupled to the laser generator 12. The frequency modulator 14 is configured to receive the laser beam (arrow L) and modulate its frequency to thereby define a modulated waveform/laser beam (arrow LM) having a predetermined three-section “sawtooth” modulated waveform, as described in greater detail below. The frequency modulator 14 may be an electro-optical frequency modulator 14 driven by an electric ramp generator (not shown) using a pre-defined modulation waveform. As will be appreciated, such an electric ramp generator may include a voltage-controlled oscillator controlled by a digital microcontroller to generate pre-defined ramp characteristics, which in turn are provided to the voltage-controlled oscillator to create the pre-defined modulation waveform. An exemplary setup for such an electric ramp generator, with an optical filter used to eliminate all frequency harmonics generated by the modulator 14, is disclosed in U.S. patent application Ser. No. 13/527,638, issued as U.S. Pat. No. 8,897,654 B1, which is hereby incorporated by reference in its entirety.
Within the representative laser system 11 configuration of
The amplifier 18 is configured to amplify power of the first portion of the modulated laser beam (arrow LM) so as to define an amplified laser beam (arrow LAMP). For example, the amplifier 18 may increase power of the laser beam (arrow LM) to a range between 1 W and 10 W (peak power). The amplifier 18 may include, for example, a single-mode fiber amplifier, and may be operated at less than a 100% duty factor, i.e., in a pulsed mode, so as to reduce power consumption. In an embodiment, the amplifier 18 may operate at a duty factor of less than about 10%, and may be enabled for a few milliseconds at a rate of about 30 Hz. Such an embodiment is non-limiting and intended herein to convey exemplary ranges and magnitudes, and therefore other embodiments of the amplifier 18 may be envisioned within the scope of the disclosure.
The Doppler lidar system 10 of
The Doppler lidar system 10 of
An additional beam splitter 19 is disposed between and interconnects the amplifier 18 and the T/R switches 25A, 25B, and 25C. The beam splitter 19 splits the amplified laser beam (arrow LAMP) into equal parts, here three equal parts as indicated by the notation “1:3”, and directs a respective one of the parts L1, L2, and L3 toward each of the T/R switches 25A, 25B, and 25C. The beam splitter 20 noted above receives the second portion of the modulated laser beam (arrow LM) from the primary laser beam splitter 16, and adjusts the polarization of the second portion of the modulated laser beam to define a local oscillator beam (arrow LO) that matches the polarization of each of the return signals L1R, L2R, and L3R.
At the beam splitter 20, the local oscillator beam (arrow LO) likewise splits into three equal parts LO1, LO2, and LO3 that are directed toward a respective one of the photoreceivers 24A, 24B, and 24C. Each of the photoreceivers 24A, 24B, 24C thus receives a respective part LO1, LO2, or LO3 of the local oscillator beam and respective reflected radiation/return signals L1R, L2R, and L3R from a respective T/R switch 25A, 25B, 25C, and then mixes the two together in an optical homodyne detection scheme, as will be appreciated by those of ordinary skill in the art. Thus, the photoreceivers 24A, 24B, 24C collectively form the homodyne receiver 24 within the context of the Doppler lidar system 10. The photoreceivers 24A, 24B, and 24C mix the local oscillator beam with the reflected radiation to produce an intermediate electrical signal with a frequency equal to the difference between frequencies of the reflected signals and the respective parts LO1, LO2, and LO3 of the local oscillator beam (arrow LO).
Within the scope of the disclosure, the SPU 50 noted briefly above is coupled to and in communication with the homodyne receiver 24. The SPU 50 analyzes and stores data received from the individual photoreceivers 24A, 24B, and 24C and uses the data to calculate and/or determine various navigational data as described below with reference to
Referring briefly to
From such data, the SPU 50 derives the LOS velocity and range to the target-of-interest, a velocity vector, altitude, and attitude at a computation node 56, which may be implemented in hardware as another processor, ASIC, chip, or within the FPGA architecture. From such data, the SPU 50 may feed a control node 58 in communication with the host vehicle 55 of
The representative setup of
Referring to
As shown in
In
Thus, a range-to-target is computed as the product of a constant (k1) and half of the frequency difference between the up and down ramp IFs, i.e.,
with k1 being a constant equal to the ratio
IF+ and IF− representing the up-ramp (FUP) and down-ramp (FDN) IFs, respectively, T being the waveform period, c being the speed of light, and B representing the signal bandwidth. Velocity is computed as the product of a second constant (k2) and half of the sum of the two side-band IFs, i.e.,
where
with λ being the wavelength. Thus, one can readily see that errors in range and velocity calculations can result if the sign of any one or more of the IFs is incorrect. The method 100 of
In terms of sign, the frequencies detected by Doppler lidar system 10 have a sign that is either positive or negative depending on the relative direction of motion between the Doppler lidar system 10 and the target/ground, and the relative distance to the target/ground. If used as part of a navigation system aboard a spacecraft during a descent and landing operation, for example, where a positive sign indicates movement of the spacecraft toward the ground and a negative sign indicates movement away from ground, accurate knowledge of the sign is critical. However, homodyne detection capabilities of the Doppler lidar system 10 via the homodyne receiver 24 shown in
In a practical system, there are times when one or more of the three intermediate frequencies (IFs) are not available during a given lidar measurement. Loss of signal can occur when a given frequency crosses from a negative sign to a positive sign, or when the signal-to-noise ratio (SNR) is too low to accurately detect. For such cases, the intermediate frequencies of one or more previous measurements are used to determine the correct signs to each of the available up-ramp IF section (FUP), the down-ramp IF section (FDN), or Doppler IF section (FDOP).
An exemplary method 100 for performing the present teachings is depicted in
In the method 100 of
Beginning with block B102, the method 100 includes defining physical constants or parameters of the lasers used by the Doppler lidar system 10 of
where B is the modulation bandwidth, T is the waveform period, and c is the speed of light. The method 100 proceeds to block B104 when the physical constants have been defined, and thereafter used in interrogating the target-of-interest. Block B102 is a baseline process used in the operation of the Doppler lidar system 10 of
Block B104 entails receiving a return waveform from an irradiated target-of-interest using the homodyne receiver 24 setup depicted in
In the next four blocks, i.e., B106, B108, B110, and B112, the SPU 50 of
At block B108, a state estimation filter is applied to estimate an expected new range and velocity, and to predict the sign of the frequencies for each of the Doppler, Up-ramp, and Down-ramp IF sections noted above. One such filter is the so-called “α, β, γ” tracking filter, as will be understood by those of ordinary skill in the art, or any other suitable filter that, unlike a Kalman filter, does not require access to a state model. Block B108 outputs slopes of Doppler, up-ramp, and down-ramp IFs of previous measurements to be used in determining the sign of the frequency of any one of the Doppler, up-ramp, and down-ramp IFs deemed invalid.
At block B110, a signal threshold is obtained from memory of the SPU 50. Block B112 then determines SNR-based false alarms given the threshold from block B110. The output of block B112 is the set of IF sections, i.e., FDN, FUP, and FDOP of
At block B114, the SPU 50 next evaluates each received IF section from block B112 to determine whether any of the received IF sections is missing. Ideally, all three IF sections FDN, FUP, and FDOP are received and available in a given measurement cycle. However, at times due to factors such as interference or loss of signal, one or more of the IF sections may not be detected by the homodyne receiver, or may be detected at a degraded level that precludes effective use of the received signal. The SPU 50 responds at block B114 by determining which IF sections are detected, available, and suitable for use.
Based on this information, and as set forth in detail below, the SPU 50 selects one of a plurality of possible sign ambiguity cases, which are nominally labeled Case 1 (block B116A described in
CASE 1—NO AMBIGUITY: Referring to
At block B131, the SPU 50 determines if the down-ramp IF section is less than the up-ramp IF section. The SPU 50 assigns a negative sign (“FDOP<0”) to the Doppler IF section in sub-blocks B136 and B139. Otherwise, the SPU 50 assigns a positive sign (“FDOP>0”) to the Doppler IF in sub-blocks B135 and B139.
Still referring to
Sub-block B118 entails computing the range and velocity using the magnitude of the Doppler IF and one of the up-ramp IF or down-ramp IF magnitudes, with the appropriate signs updated in sub-block B139 as noted above.
CASE 2—FDOP and FUP: Referring to
At sub-block B142, the SPU 50 determines if the product of the slope of the Doppler IF section (ΔFDOP) and up-ramp IF section (ΔFUP) of the previous block B108 estimate is positive. If so, block B116B proceeds to sub-block B144. The SPU 50 executes sub-block B146 in the alternative when the product is negative.
At sub-block B143, the SPU 50 determines if the product of the slope of the Doppler IF section (ΔFDOP) and up-ramp IF section (ΔFUP) previously estimated at block B108 is negative. If so, block B116B proceeds to sub-block B146. The SPU 50 executes sub-block B149 in the alternative when the product is positive.
Sub-block B144 includes assigning a negative sign to each of the Doppler IF section and the up-ramp IF section before proceeding to block B118 (
Sub-block B146 respectively entail assigning a positive sign to the Doppler IF section and a negative sign to the up-ramp IF section before proceeding to block B118.
Sub-block B149 entails assigning a positive sign to each of the Doppler IF section and the up-ramp IF section before proceeding to block B118 of
CASE 3—FDOP and FDN: Referring to
At sub-block B153, the SPU 50 determines if the product of the slope of Doppler IF section and the down-ramp IF section is negative, i.e., ΔFDOP*ΔFDN<0. If so, block B116C proceeds to sub-block B158. Sub-block B159 is executed in the alternative when ΔFDOP*ΔFDN>0.
At sub-block B154, the SPU 50 determines if the product of the slope of the Doppler IF section and that of the down-ramp IF section is positive, i.e., ΔFDOP*ΔFDN>0. If so, block B116C proceeds to sub-block B156. Block B116C instead proceeds to sub-block B158 when ΔFDOP*ΔFDN<0.
At sub-block B156, the SPU 50 assigns a negative sign to each of the Doppler IF section and the down-ramp IF section, and then proceeds to sub-block B118 (
At sub-block B158, the SPU 50 assigns a negative sign to the Doppler IF section and a positive sign to the down-ramp IF section, and then proceeds to sub-block B118 (
Sub-block B159 includes assigning, via the SPU 50, a positive sign to each of the Doppler IF section and the down-ramp IF section before proceeding to sub-block B118 (
CASE 4—FDOP only: Referring to
At sub-blocks B163 and B164, the SPU 50 respectively assigns a positive sign (sub-block B163) or a negative sign (sub-block B164) to the Doppler IF section, and then proceeds to block B118A.
CASE 5—FDN and FUP: Referring to
At sub-block B171, the SPU 50 next determines if the product of the slope of the up-ramp IF section (ΔFUP) and the down-ramp IF section (ΔFDN) of the previous estimate (block B108 of
At sub-block B173, the SPU 50 determines if the product of the slope of the up-ramp IF section (ΔFUP) and the down-ramp IF section (ΔFDN) estimate is negative. The SPU 50 proceeds to sub-block B174 when this case is true. The SPU 50 executes sub-block B175 in the alternative when the product is positive.
At sub-blocks B172, B174, and B175, the SPU 50 of
CASE 6—No IF sections, or FDN only, or FUP only: As noted above, Case 6 is an “invalid state” condition. Thus, selection of Case 6 at block B114 (
Referring once again to
At block B120, the SPU 50 obtains slopes of the state estimates from block B118, records the slopes, and then proceeds to block B122. As will be appreciated, a slope is a rate of change of the state estimate, and therefore a predetermined past number of state estimates may be used to determine the slope to provide the desired immunity to noise or errant state estimates.
Block B122 includes setting a data status bit flag in memory indicative of the case number and completion of an iteration of the method 100. The method 100 then proceeds to block B124.
At block B124, the SPU 50 determines if a new measurement by the Doppler lidar system 10 has been received. If so, the method 100 returns to block B104, with the method 100 otherwise being complete (“End”).
The method 100 described above with reference to
As will be appreciated, the Doppler lidar system 10 of
Precision range and velocity data provided by the Doppler lidar system 10 of
Additionally, the disclosed Doppler lidar system 10 can also make a dramatic impact on the ongoing efforts for achieving safe and efficient operation of autonomous ground vehicles. By scanning its laser beam in
While aspects of the present disclosure have been described in detail with reference to the illustrated embodiments, those skilled in the art will recognize that many modifications may be made thereto without departing from the scope of the present disclosure. The present disclosure is not limited to the precise construction and compositions disclosed herein; any and all modifications, changes, and variations apparent from the foregoing descriptions are within the spirit and scope of the disclosure as defined in the appended claims. Moreover, the present concepts expressly include any and all combinations and subcombinations of the preceding elements and features.
This patent application claims the benefit of and priority to U.S. Provisional Patent Application No. 63/030,215, titled “Ambiguity Mitigation for FMCW Sawtooth Waveform,” filed on May 26, 2020, the contents of which are hereby incorporated by reference in their entirety.
The invention described herein was made in the performance of work under a NASA contract and by employees of the United States Government and is subject to the provisions of Public Law 96-517 (35 U.S.C. § 202) and may be manufactured and used by or for the Government for governmental purposes without the payment of any royalties thereon or therefore. In accordance with 35 U.S.C. § 202, the contractor elected not to retain title.
Number | Date | Country | |
---|---|---|---|
63030215 | May 2020 | US |