This disclosure relates to quadrature signal processing of local oscillator and Doppler frequency-shifted signals in a lidar or other coherent optical systems.
A primary obstacle of fiber lidar is assumed to be the birefringent depolarization of the local oscillator (LO) signal from the transmitted carrier after splitting from the lidar output path. The effect can destroy the heterodyne efficiency at the detector and hence lidar operation unless polarization preserving fiber is utilized in the system past the split point in homodyne systems. This effect is assumed worse in heterodyne systems utilizing different LO and transmitter sources. The only form of the optical fiber lidar “immune” from this effect utilizes a local oscillator signal taken from the Fresnel reflection at the end of the transmit fiber immediately preceeding the output telescope. However, this latter mode of operation is not required as conventionally assumed. Laboratory tests have shown that phonon modulation of the birefringence in the local oscillator path gives rise to AM modulation of the detected signals within the dynamic range required of the lidar to perform its basic task. This provides a statistically detectable signal.
Furthermore, in conventional lidar systems, a frequency offset between a local oscillator signal and a transmitted beam has been traditionally required. This has traditionally been achieved in homodyne operation via a frequency shifting device such as an expensive acousto-optic (A/O) modulator, or in heterodyne operation by maintaining a fixed offset between the frequencies of the two coherent sources. It is desirable to perform such heterodyning or homodyning without the use of such acousto-optic modulators.
The disclosed invention can be used in free-space lidar systems, fiber lidar systems, and other systems based upon coherent mixing to eliminate the costly A/O cell used for offset homodyne operation or the difficult to stabilize offset heterodyne source. These elements are replaced with inexpensive detectors and couplers with savings of several thousands of dollars. The use of the disclosed invention allows the effective use of non-polarized or polarization preserving fibers, depending on the coherent system design requirements. The disclosed invention can be utilized effectively in the presence of birefringent de-polarization.
Signal to noise ratio for the disclosed technique is within 3 dB of that engendered by the use of the typical A/O cell, but alignment and temperature sensitivities are considerably reduced. Further, the bandwidth requirements necessary in the processing electronics are cut in half relative to the A/O modulator or offset heterodyne systems. Lastly, the electronic support components required for the other system forms are eliminated with considerable savings in volume and electronic power. The use of multiple coherent wavelengths can be achieved with this disclosed invention
The disclosed technique enables considerably more compact systems to be fabricated and cost effectively extends the applicability of the typical fiber lidar into a wider range of applications that require fall signed Doppler spectrum (vector velocity). Typical applications that will see substantial benefit include vibration sensing, turbulence sensing and velocity lidars (e.g. police radar applications, relative motion sensing applications, optical air data systems, etc.) of any type (e.g. linear velocity, tangential velocity, spin sensing, etc.)
Applications for coherent Doppler lidars include velocity sensing applications (platforms and objects), volumetric/fluidic flow sensing, vibration monitoring, range to target and other related standoff sensing applications. The lidar detects the Doppler frequency shift imposed on coherent light scattered from a moving target by mixing the scattered, frequency shifted light with a reference beam of light (local oscillator) which is not shifted in frequency on the detector. A difference frequency results from this mixing process which is proportional to the velocity of the scattering medium. It is the Doppler frequency shift imposed on the light scattered from the target that provides the mechanism used for velocity detection. The reference beam can be either derived from the transmit beam (homodyne operation) or derived from another stable coherent source (heterodyne operation). By measuring the Doppler shift from three (or more) spatially separated lidar beams a complete vector velocity can be computed along with statistical velocity information.
In general, fiber lidar systems utilize the same optical functions to perform the lidar mission, except the optical elements are created by guided-wave optics (i.e. optical fiber devices). The laser source is generally a combination of a suitable solid state, DFB laser diode and one or more cascaded optical fiber amplifiers of the appropriate wavelength, although fiber or free-space lasers could be used as the source elements. For the most part, the amplifier of choice is the erbium-doped fiber amplifier (EDFA) operating at a wavelength of 1.54 mu.m. In one embodiment of an offset homodyne fiber lidar 100 shown in
If the optical fiber quadrature processing assembly 200 shown in
(sin(.+−..omega..sub.d)t−j cos(.+−..omega..sub.d)t)
to develop a signed velocity spectrum. While the equations below are used in RF spectrum analysis and are standard in communications textbooks for illustrating Fourier transform theory, heretofore it has not been connected to optical lidar signal processing using the phase characteristics of the coupled-mode equation.
Signals in a single mode, directional optical fiber coupler (fused, integrated optics, etc.) have a −90.degree. phase shift in a transferred evanescent wave arm relative to the “straight through” fiber path due to the requirements of the wave equations for coupled waveguide solutions. This fact can be used as to develop in-quadrature signals for the spectrum analysis process that resolves the Doppler frequency and directional ambiguity in a Doppler based LIDAR (fiber or free-space based) used for velocity measurements. A shift in frequency is imposed on the transmitted light beam of a LIDAR (lidar) by the velocity of any object from which the light is reflected (i.e. the Doppler effect). However, a velocity magnitude toward or away from the lidar beam will generate the same differential frequency in the standard heterodyne process. This “directional ambiguity” must be resolved from the sign change in the axial vector velocity (i.e. change of velocity direction along a given axis) by use of the absolute frequency of the optical wave, by use of an offset frequency or via phase information relative to the carrier. The absolute carrier frequency is too high to work with in the electronic domain and the use of an offset frequency via an expensive acousto-optic cell (or other frequency shifting device), though conventionally used, is not to be preferred. The disclosed technique therefore develops the required information from the phase domain of signals.
The Doppler frequency shift in a lidar is related to the velocity according to the equation: 1d=−4 V s (rad/sec) or (1a) f d=−2V s(Hz) (1b)
where V is the target velocity in meters per second and .lambda..sub.s is the laser source wavelength in the medium.
The network or array of signal couplers 200 illustrated in
E.sub.1=−E.sub.s cos [(.omega..sub.c.+−..omega..sub.d)t]+E.sub.lo sin(.omega..sub.lo)t (2)
where E.sub.s and E.sub.lo are the vector magnitudes of the signal and local oscillator field strengths respectively, .omega..sub.c is the radian frequency of the transmitted optical carrier beam and .omega..sub.d is the radian frequency of the Doppler shift imposed on the light by moving target. The sign of omega..sub.d is dependent on the direction of the velocity vector and is positive if the target is moving toward the beam (or lidar) and negative if it is moving away from the beam (or lidar). In general, omega..sub.d is a spectrum of frequencies with a bandwidth determined by the target velocity, surface figure, etc. Similarly, omega..sub.c and omega..sub.lo likewise have a finite bandwidth that is dependent on the laser source(s) being used in the lidar. For the purposes of the current development, omega..sub.d, .omega..sub.c, and omega..sub.lo may be assumed to be radian frequencies of zero bandwidth. The total signal content after processing is then simply the sum of the power spectral densities of each signal's bandwidth after mixing in the optical detectors. Likewise, at the second optical detector, 216:
E.sub.2=−E.sub.s cos [(.omega..sub.c.+−.omega..sub.d)t]+E.sub.lo cos(.omega..sub.lo)t (3)
The detected signal currents are proportional to the power in the field and therefore, proportional to the square of the total field vector on each detector 214, 216. This fact is what causes the frequencies on the detectors to mix or “heterodyne.” It is assumed that the polarizations of E.sub.s and E.sub.lo have been adjusted to achieve linear addition of the field vectors (essentially a heterodyne efficiency of unity). This is usually achieved by the use of polarization preserving waveguide structures, but birefringent structures associated with normal optical fiber guides will work well under most conditions where some compromise in signal to noise ratio may be offset with temporal averaging of the results. Returning to the signal current, under the given assumptions the intensity of the signals detected is, for example at the first detector 214:
I.sub.s.varies..vertline.E.sub.1.vertline..sup.2 (4)
Therefore, working with detector 214, the in-phase signal is:
I.sub.P.varies..vertline.E.sub.s.vertline..sup.2 cos .sup.2[(.omega..sub.c.+−..omega..sub.d)t]+.vertline.E.sub.lo.vertline..su− p.2 sin(.omega..sub.lo)t)−2.vertline.E.sub.s.parallel.E.sub.lo.vertline.c− os[(.omega..sub.c.+−..omega..sub.d)t] sin(.omega..sub.lo)t (5)
The first two terms in proportionality (5) comprise the DC current term in the equation, which are removed by filters in the processing system 328 (
2l p=−2 E s E lo [12 sin(cd+lo)t−12 sin(cd−lo)t] (6)
In equation (6).omega..sub.d is very small in comparison to .omega..sub.c or .omega..sub.lo and the average radian frequencies of these two terms are essentially equal as they are derived by splitting a single laser source (homodyne operation), i.e. .omega..sub.lo=.omega..sub−.c. If these two terms are derived from separate sources (heterodyne operation), the theory of the calculations will not change, however the measured Doppler frequency will deviate from the assumed condition by an offset equal to the frequency difference between the carrier and local oscillator laser.frequencies (.omega..sub.d=.omega..sub.d,true+.omega.su− b.offset). This issue can be ignored in the current calculations as the offset can be later added to the result. Therefore, provided sufficient coherence length is available in the laser source(s) such that .omega..sub.lo (t)=.omega..sub.c (t), the sum frequencies are absorbed by the detector material as loss, leaving I.sub.P=+E.sub.sE.sub.lo sin(.+−..omega..sub.d)t (7)
Similarly, the signal current in detector 216, the in-quadrature signal, may be calculated as: 3IQ−2E s E lo cos [(cd)t] cos(lo)t=−E s E lo cos (d)t(8)
It can be seen from equations (7) and (8) that the two Doppler, photo signal currents are separated by 90 degrees in phase and are therefore in-quadrature. To process the Doppler velocity then the signals are summed and the complex Fourier Transform is taken as follows: 4 F( )=E s E lo−.infin.+.infin. [sin(d)t−j cos (d)t]exp{−jt}t (9)
Using Euler's identity: exp{.+−.jX}=sin(x).+−.j cos(x), then: 5 F ( )=lim a/2->.infin. E s E lo−a 2+a/2exp{−j(d)t}exp{−jt}=lim a/2->.infin.E s E lo−a/2+a/2exp{−j(d)t}t=lim a/2->infin. E s E lo−j(d)t exp{−j(d)t}−a/2a/2-lim a/2->infin. Es E lo-j (d) a/2exp{−j(d)t}—a 2 a/2 (10)=lim a/2->.infin.E s E lo−j(d)a/2 j2 sin [(d)]a/2(11)
Mathematically, equation (11) then describes a frequency magnitude spectrum that is a zero bandwidth delta function with magnitude proportional to the product of E.sub.sE.sub.lo and a power spectral density proportional to .vertline.E.sub.sE.sub.lo.vertline..sup.2 at a radian frequency of .omega.=+.omega..sub.d or−.omega..sub.d according to the vector direction of the target moving toward or away from the lidar respectively. The final equation is then:
F(.omega.)−2.pi.E.sub.sE.sub.lo.delta.(.omega..+−..omega..sub.d) (12)
As was previously noted, if a finite bandwidth is associated with the laser source, local oscillator and/or target motion, the delta function of equation (12) is repeated over a power spectral density function whose width is equal to the sum of source bandwidth, local oscillator bandwidth and any additional bandwidth resulting from the target modulation effects. The center frequency of the distribution however, is still .omega..sub.d and its sign is either positive or negative in accordance with the direction of the Doppler shift. Thus analysis of the Fourier spectrum computed from the quadrature signals and equation (11) will yield both the magnitude spectrum of the Doppler signals (which may be further processed for velocity magnitude according to the equations 1a or 1b) and the sign of the velocity vector (inherent in the positive or negative sign of the frequency in the Fourier plane).
Referring to
Referring to
Referring to
In
A fifth signal coupler 206 is receptive of the unphase-shifted Doppler frequency-shifted signal 224 and the unphase-shifted local oscillator signal 226. The unphase-shifted Doppler frequency-shifted signal 224 and the unphase-shifted local oscillator signal 226 are mixed in the fifth signal coupler 206 and the unphase-shifted Doppler frequency-shifted signal 224 is phase-shifted by −90 degrees. The fifth signal coupler 206 provides as output an in-phase signal 212 which includes the unphase-shifted local oscillator signal 226 and the −90 degree phase-shifted Doppler frequency-shifted signal.
The in-phase signal 212 and the in-quadrature signal 218 are provided as input to optical detectors 214, 216 which provide as output electrical signals 220, 222 indicative of the intensities, I.sub.P and I.sub.Q, of the in-phase and in-quadrature signals 220, 222. Fourier transforming the complex sum of the in-phase and in-quadrature signals yields a frequency spectrum centered around zero with the sign of the power spectral density components representing the sign of the vector velocity in the lidar beam 140 axis. The processing bandwidth is effectively one half of that which is required using a conventional A/O cell.
This method and apparatus can be achieved in the electronic domain under conditions in which tracking the Doppler frequency through zero velocity (zero frequency) is not necessary, i.e. a velocity scenario in which the Doppler frequency is unipolar and sufficiently displaced from zero at all times. However, the dynamic range and simplicity of the optical system disclosed is superior under all conditions and is therefore to be preferred under most circumstances supported by the photonics of the lidar itself. It should also be noted that this technique can be implemented in free-space optics with optical analogs (beam splitters and waveplates) to the fused waveguide couplers originally intended and to a limited degree in multi-mode optical waveguides. In this regard, it has not been obvious to the user community that the phase shift of the waveguide coupler may be used in manner disclosed.
In the case of polarized fiber systems, coupler 208 in
Also shown in
Thus, based upon the foregoing description, a quadrature processed lidar system is disclosed with application general coherent optical systems. While preferred embodiments have been shown, and described, various modifications and substitutions may be made thereto without departing from the spirit and scope of the invention. Accordingly, it is to be understood that the present invention has been described by way of illustration only, and such illustrations and embodiments as have been disclosed herein are not to be construed as limiting the claims.
Number | Date | Country | |
---|---|---|---|
Parent | 10969964 | Oct 2004 | US |
Child | 11882280 | Jul 2007 | US |
Parent | 10323756 | Dec 2002 | US |
Child | 10969964 | Oct 2004 | US |