The present invention relates generally to radar systems, and more particularly to frequency modulated continuous wave (FMCW) radar systems.
Radar is a system that uses electromagnetic waves to identity the range, direction, and/or velocity of both moving objects and stationary targets. For example, radar is often used to detect weather conditions, ships, aircraft, motor vehicles, geological formations, as well as many other applications.
To facilitate this functionality, radar systems include a radar transmitter that transmits electromagnetic waves, such as radio waves, which are scattered or reflected by a target. A radar receiver, which is typically in approximately the same location as the transmitter, then receives the scattered wave and analyzes it to determine the range, direction, and/or velocity of the target. Because radio frequency signals are relatively easy to amplify, even though the scattered wave is usually very weak, the receiver can amplify the received scattered wave so it can be suitably processed. Therefore, radar systems can detect objects at ranges where other electromagnetic waves, such as sound or visible light, would be too weak to detect.
As will be appreciated from the above discussion and embodiments described herein, there is an on-going need for improvements in radar systems.
The following presents a simplified summary in order to provide a basic understanding of one or more aspects of the invention. This summary is not an extensive overview of the invention, and is neither intended to identify key or critical elements of the invention, nor to delineate the scope thereof. Rather, the primary purpose of the summary is to present some concepts of the invention in a simplified form as a prelude to the more detailed description that is presented later.
One embodiment of the present invention relates to a system and method for detecting a range and a velocity of a target. In this method, an electromagnetic wave is transmitted over a frequency range, where a period of the wave comprises a number of consecutive ramps. A first ramp in the period is transmitted over a first portion of the frequency range, and a second ramp in the period is transmitted over a second portion of the frequency range that differs from the first portion. The second ramp is offset by a frequency shift relative to the first ramp. A scattered wave is received from the target and processed to determine the range and the velocity of the target.
The following description and annexed drawings set forth in detail certain illustrative aspects and implementations of the invention. These are indicative of only a few of the various ways in which the principles of the invention may be employed.
One or more implementations of the present invention will now be described with reference to the attached drawings, wherein like reference numerals are used to refer to like elements throughout.
Referring to
In one embodiment, radar system 100, and the transmitter 102 and the receiver 104, comprises a continuous wave radar system. Continuous wave radar works by detecting a phase shift between the transmitted wave 108 and the scattered wave 112, and then analyzing the phase shift to determine information about a target 114. To facilitate this detection capability, the transmitter 102 changes the frequency of the transmitted wave 108 as a function of time. After being transmitted, there is a delay before the transmitted wave 108 reflects off the target 114 and the scattered wave 112 is received by the receiver 104. For FMCW radar systems a frequency ramp is generated and transmitted. The scattered signal 112 is down-converted in a mixer with the help of the transmitted signal. The received signal is a delayed frequency ramp having a phase shift that is a function of the time of flight. Multiplication of the transmitted and received signal, which is performed by a mixer, provides a baseband signal having a frequency that is directly proportional to the distance of the target from the radar system 100. During this delay, the transmitter 102 shifts to transmit at some other frequency. Generally, the amount of shift between the transmitted wave 108 and the scattered wave 112 is greater over longer times, so a greater phase shift may tend to correspond to longer distances, for example. In this way, the range to the target 114 can be determined in one embodiment. Unambiguously determining both the velocity and range of a given target, however, can be difficult.
As shown, the electromagnetic wave 300 spans a total frequency range, Ft, and has a period, Tt, which in some embodiments corresponds to a period of a baseband frequency. Within the illustrated period, Tt, the electromagnetic wave 300 includes a series of N consecutive ramps that each has a duration of Ts, where the ramps are transmitted in ascending or descending order with a frequency shift, Δf, therebetween. Thus, the N consecutive ramps cover the frequency range, Ft, but each consecutive ramp covers only a fraction of the frequency range Ft. In the illustrated embodiment, the ramps are triangular ramps. However, in other embodiments the ramps can also be substantially triangular and/or have other geometries.
More specifically, the wave 300 includes a first consecutive ramp (Ramp 1) transmitted over a first portion FP1 of the frequency range FT. The first consecutive ramp (Ramp 1) includes a first ramp segment 302 that is piece-wise continuous with a second ramp segment 304. A second consecutive ramp (Ramp 2) is transmitted over a second portion FP2 of the frequency range Ft and includes a third ramp segment 306 that is piece-wise continuous with both the second ramp segment 304 and a fourth ramp segment 308. Often, the first and third ramp segments 302, 306 have one slope (e.g., a positive slope), and the second and fourth ramp segments 304, 308 have another slope (e.g., a negative slope). Again, the first and second consecutive ramps (Ramp 1, Ramp 2) are offset by the frequency shift, Δf, relative to one another. Additional consecutive ramps are also formed within the period Tt in a similar manner as shown.
Along each ramp there are, in various embodiments, L sampling points 310. In essence, the radar system takes measurements (samples) of both the transmitted wave and received scattered wave at these times, and then measures a phase difference between the transmitted and received signals. This phase difference can then be used to unambiguously determine the range and velocity of a target. As shown, consecutive sampling times are separated by a sampling interval, TA. There is a frequency step, fstep, between frequencies associated with consecutive sampling times.
Referring now to
To transmit the desired electromagnetic wave 300, the memory 402 stores data related to the ramps of the electromagnetic wave 300. More specifically, an increment memory 420 stores the slope associated with each ramp, a cycle memory 422 stores the duration for which each slope is to be transmitted, and a delay memory 424 determines when to start the control state machine 404 and analog-to-digital (A/D) converter 416.
The control state machine 404 periodically reads from the memory 402 by incrementing a memory address to facilitate the desired ramp shaping. For example, to create the first ramp, the control state machine 404 reads the increment memory 420 to find the slope of the first ramp (e.g., “2”), reads the cycle memory 422 to ascertain the duration of the first ramp (e.g., “10 clock cycles”), and reads the delay memory 424 for the delay (e.g., “100 clock cycles”). Thus, after receiving the “start” signal, the control state machine 404 waits (delay) 100 clock cycles, and then provides a ramp slope of 2 for the next 10 clock cycles to the accumulator 406.
In one embodiment the accumulator 406 adds the value of the present slope to a present count value stored therein. In this manner, the accumulator 406 outputs a series of un-interpolated frequencies to the interpolation control 408. Over time, the un-interpolated frequencies are used to form the ramps as shows in
The interpolation control circuit 408 may work in conjunction with the calibration memory tuning law circuit 410 to ensure that the ramp segments that are actually transmitted are linear. To this end, before the ramp segment generation starts (e.g., during the delay cycles), the calibration memory tuning law circuitry 410 characterizes the behavior of the VCO 508 over the frequency at which the electromagnetic waves will be transmitted. If the tuning law of the VCO is nonlinear, the calibration memory tuning law circuitry 410 stores data indicative of this non-linearity. During transmission, the interpolation control circuit 408 receives the present frequency to be transmitted from the accumulator 406, and calculates the necessary D/A-converter samples for a linear ramp. based on the data provided by the calibration memory tuning law circuit 410.
In one embodiment the linear interpolation circuit 412 includes subtractors 426, 428; mixers 430, 432; and an adder 434 that provide a final interpolated frequency value to the D/A converter 414. The D/A converter 414 converts this final interpolated frequency value into analog form and passes it to the VCO 508.
On the receiver-side, the circuit 400 receives the baseband signal (which has been down-converted), and converts it to a digital signal at the A/D converter 416. The A/D controller 418 works in conjunction with the A/D converter 416 to provide the sampled data that is representative of the scattered wave.
Referring now to
Now that some example embodiments of systems that can achieve advantageous radar techniques have been discussed, reference is made to
In
At 604, a scattered wave is received, wherein the scattered wave is reflected from a target.
At 606, a phase shift is measured between the ramp being presently transmitted and the received scattered wave presently received.
At 608, the phase shift and frequency of the waves are analyzed to determine the target's unambiguous range and velocity. Therefore, both the range and velocity of the target can be measured.
In
In 704, the scattered wave is received after reflecting off a target, and is down-converted to a down-converted signal.
In 706, sampled data is obtained by sampling the down-converted signal.
In 708, an inverse two-dimensional inverse Fast Fourier-transform (2D-IFFT) is performed on the sampled data.
In 710, a constant false alarm rate (CFAR) algorithm is applied to the result of the 2D-IFFT to unambiguously determine both the range and velocity of the target.
For purposes of understanding some implementations of the invention, what follows below is a mathematical explanation of a measurement principle in accordance with a ramp sequence such as set for forth in
The transmitter frequency for the first ramp in
The phase φT(t) of the transmitted signal cos(φT(t)) becomes after integration:
where fc is the frequency of a carrier wave, and fstep and TA are as shown and described in
τ is the delay between the transmitted and received signal of one target. The last term in the equation above can be neglected because τ/T<<1. For calculation of the delay τ=2(R+υ·t)/c, a target at distance R with constant velocity υ is assumed. This leads to:
The last term is called Range-Doppler-Coupling and can be neglected again:
As described in
With τ=2(R+υ·t)/c the phase difference of the equation above can be written as:
The phase difference is evaluated at the discrete time t=n·Ts+l·TA. Sorting the terms according to n and l leads finally to:
If an I/Q-Demodulator is present in the system, the received signal is:
sIF=cos(Δφ(n·Ts+l·TA))−i sin(Δφ(n·Ts+l·TA))=e−iΔφ(n·T
In order to detect targets, an inverse two-dimensional Fast Fourier transform (2D-IFFT) can be applied:
N is the number of ramps within a given period, and L is the number of sampling points on one ramp. NZ and LZ are the number of points in the 2D-IFFT after zero-padding. The window function w(l,n) is used to reduce the side lobes after the IFFT. w(l,n) is set to one to simplify the following equations without restricting the generality. The equation above can be re-ordered to that two-dimensional spectrum can calculated from one-dimensional IFFTs.
A peak of a target with distance R and velocity υ occurs therefore at the following position:
The position of one peak, that has been detected at k and p, can be used to calculate the distance R and velocity v by solving the equations above:
The range resolution along the k-axis can be derived from equation 12:
The velocity resolution is given by:
The range resolution along the p-axis is related to the bandwidth of a short ramp:
The velocity resolution is determined by the duration of the short ramp:
In some embodiments, FMCW-radar system comprise a single mixer instead of an I/Q-Demodulator in order to reduce current consumption and space for a second mixer and the 90°-phaseshift module. As a consequence, half of the unambiguous range is lost. The presented principle also applies to a single mixer embodiment. In equation 9 the Q-signal is lost and cos(Δφ) remains in case of a single mixer, which leads to:
SIF=cos(Δφ(n·Ts+l·TA) (20)
The 2D-IFFT becomes:
The term cos(Δφ) can be written as ei·Δφ−i·ei·Δφ:
A peak of one target occurs therefore not only at k and p, but also at −k and −p
This means that a rotation of 180° around the center maps the 2D-IFFT-spectrum onto itself again. The unambiguous range is therefore reduced by a factor of 2 compared to the embodiment with an I/Q-Demodulator:
This is also true for the unambiguous velocity range:
The peak position p is nearly independent of the velocity and is positive, because the duration Ts is rather short in a practical system. The slope of the short ramp allows application of the principle without the I/Q-Demodulator. The parameter Δf is chosen according to equation 26 so that the whole distance range of interest is covered by Lz/2. If the peak position occurs in pε[0, Lz/2], the true position lies at point P=(k,|p|). Otherwise the true position is P=(−k,|p|) for pε[−Lz/2+1,0] or pε[Lz/2+1,Lz]. The decision rule is finally:
p≧0P=(k,|p|) (28)
p<0P=(−k+r·Nz|−p|)rε{0,1} (29)
For the case p<0 the velocity v can be calculated according to equation 27 for the two possibilities rε{0,1}. Thus, a velocity for each of the two possibilities will be calculated. The velocity that fulfills the unambiguous range from equation 27 is the true velocity for the target and the other velocity can be discarded. The case r=0 occurs if k becomes negative due to a negative velocity. The other case r=1 covers situations in which k is greater than Nz/2 due to a positive velocity. Other cases of aliasing, where the velocity is higher than the unambiguous range, are often not detected.
A simulation example is given in
R=1.5 m, υ=50 km/h,
k would become negative, so the peak appears at k=21, p=249
Case: p<(0, r=0P=(−20, 256−248)
R=5 m, υ=20 km/h,
Peak Position: k=88, p−27
Case: p≧0P=(88,27)
R=20 m, υ=−220 km/h
Peak Position: k=355, p=147
Case: p<0,r=1P=(−355+Nz=669,109)
If the number of sampling points on a ramp segment is low, the peak along the p-axis becomes very broad and the peak position p is hard to estimate. This influences the accuracy of the calculated distance and velocity of one target. In order to overcome this problem, two 2D-IFFTs can be calculated: the first one without the last point of every ramp segment, the second one without the first point of every ramp segment.
The peak of one target occurs at the same position k, p of the two 2D-IFFTs. But there is a phase difference between the two peaks according to the frequency shift fstep and time shift TA.
Together with equation 23 the phase difference ΔφR can be used to calculate the distance and velocity more precisely:
The phase difference ΔφR can also be used instead of equation 28 to determine the sign of k:
ΔφR≧0P=(k,|p|) (35)
ΔφR<0P=(−k+r·Nz,|−p|)rε{0,1} (36)
DC-Offset: Unwanted signal due to imperfect system
Moving corner reflector at R=1.82 m, υ=12 km/h
Door frame at R=4.3 m, υ=0.48 km/h≈0
The distance and velocity were calculated with equations 33 and 34. A surface plot of the scene is depicted in
Some methods and corresponding features of the present disclosure can be performed by hardware modules, software routines, or a combination of hardware and software. To the extent that software is employed, for example by a baseband processor or other processor associated with the power amplifier, the software may be provided via a “computer readable medium”, which includes any medium that participates in providing instructions to the processor. Such a computer readable medium may take numerous forms, including but not limited to, non-volatile media, volatile media, and transmission media. Non-volatile media includes, for example, optical disks (such as CDs, DVDs, etc.) or magnetic disks (such as floppy disks, tapes, etc.). Volatile media includes dynamic memory, such as ferroelectric memory, SRAM, or DRAM. Transmission media includes coaxial cables, copper wire, fiber optics, etc. that could deliver the instructions over a network or between communication devices. Transmission media can also include electromagnetic waves, such as a voltage wave, light wave, or radio wave.
Although the invention has been illustrated and described with respect to one or more implementations, alterations and/or modifications may be made to the illustrated examples without departing from the spirit and scope of the appended claims. For example, although several examples above are discussed with regards to triangular or substantially triangular ramps, other non-triangular ramps are also contemplated as constituting legal equivalents. In particular regard to the various functions performed by the above described components or structures (assemblies, devices, circuits, systems, etc.), the terms (including a reference to a “means”) used to describe such components are intended to correspond, unless otherwise indicated, to any component or structure which performs the specified function of the described component (e.g., that is functionally equivalent), even though not structurally equivalent to the disclosed structure which performs the function in the herein illustrated exemplary implementations of the invention. In addition, while a particular feature of the invention may have been disclosed with respect to only one of several implementations, such feature may be combined with one or more other features of the other implementations as may be desired and advantageous for any given or particular application. Furthermore, to the extent that the terms “including”, “includes”, “having”, “has”, “with”, or variants thereof are used in either the detailed description and the claims, such terms are intended to be inclusive in a manner similar to the term “comprising”.
Number | Name | Date | Kind |
---|---|---|---|
3750171 | Faris | Jul 1973 | A |
4388622 | Fletcher, Jr. | Jun 1983 | A |
4503433 | Tomasi | Mar 1985 | A |
5111206 | Madni et al. | May 1992 | A |
5349358 | Canal | Sep 1994 | A |
5442359 | Rubin | Aug 1995 | A |
5483242 | Grein et al. | Jan 1996 | A |
5619208 | Tamatsu et al. | Apr 1997 | A |
5963163 | Kemkemian et al. | Oct 1999 | A |
6011507 | Curran et al. | Jan 2000 | A |
6104336 | Curran et al. | Aug 2000 | A |
6392588 | Levanon | May 2002 | B1 |
6864832 | Mende et al. | Mar 2005 | B2 |
7737885 | Winkler | Jun 2010 | B2 |
7821443 | Winkler et al. | Oct 2010 | B2 |
7932855 | Uebo | Apr 2011 | B2 |
20050046607 | Volman | Mar 2005 | A1 |
20050156780 | Bonthron et al. | Jul 2005 | A1 |
20050225481 | Bonthron | Oct 2005 | A1 |
20100085233 | Wintermantel et al. | Apr 2010 | A1 |
20100289692 | Winkler | Nov 2010 | A1 |
20110037641 | Steffens | Feb 2011 | A1 |
20110050484 | Nakanishi | Mar 2011 | A1 |
Number | Date | Country |
---|---|---|
19602437 | Jul 1996 | DE |
100 50 278 | Apr 2002 | DE |
Number | Date | Country | |
---|---|---|---|
20100289692 A1 | Nov 2010 | US |