Вернуться к поиску Вернуться к таблице

 Выделить слова:

 
Патент США N 7852259
Falk December 14, 2010
Clutter filtering

РЕФЕРАТ

Clutter filtering of broadband radar signals is performed by obtaining a measure of the impedance of an echo target. A distance resolution is selected and then a power spectrum is calculated for the necessary transmitting signal. Further a target area profile .rho.(.tau.) is estimated by utilising a correlation between the radiated signal and the received signal in form of a convolution. The target impedance is then calculated, and thereafter, the target can be filtered out using the achieved impedance characterising the echo of the target.

Автор(ы): Falk; Kent Olof (Go{umlaut over ( )}teborg, SE)
Заявитель: Telefonaktiebolaget LM Ericsson (publ) (Stockholm, SE)
Family ID: 42933960
N заявки: 11/487,585
Приоритет: July 17, 2006
Related U.S. Patent Documents
Application NumberFiling DatePatent NumberIssue Date
PCT/SE2004/000077 Jan., 2004
10557466 Nov., 2005 7573419

Current U.S. Class: 342/159 ; 342/118; 342/145; 342/189
Класс МПК: G01S 13/00 (20060101)
Current CPC Class: G01S 7/414 (20130101); H01Q 21/08 (20130101); H01Q 21/205 (20130101); H01Q 21/28 (20130101)
Field of Search: 342/118,145,159,189

Ссылки на источники [Referenced By]
Патентные документы США
3281842 October 1966 Cerer et al.
3346862 October 1967 Raudeep
3441931 April 1969 Shrader
4042928 August 1977 Altes
4099182 July 1978 Ward
4114153 September 1978 Neidell
4150378 April 1979 Barton
4450531 May 1984 Kenyon et al.
4644356 February 1987 Yamano
4719468 January 1988 Jehle et al.
4894660 January 1990 Thomson et al.
5093670 March 1992 Braathen
5140332 August 1992 Martin et al.
5210820 May 1993 Kenyon
5351053 September 1994 Wicks et al.
5446461 August 1995 Frazier
5519400 May 1996 McEwan
5644314 July 1997 Ahmad et al.
5734346 March 1998 Richardson et al.
5877726 March 1999 Kudoh et al.
6014099 January 2000 Bennett et al.
6075479 June 2000 Kudoh
6114987 September 2000 Bjornholt
6115113 September 2000 Flockencier
6131013 October 2000 Bergstrom et al.
6147658 November 2000 Higashi et al.
6397154 May 2002 Jones et al.
6518914 February 2003 Peterson et al.
6539308 March 2003 Monk et al.
7065162 June 2006 Sorrells et al.
7109939 September 2006 Lynch et al.
2002/0085624 July 2002 Lomp
2002/0190915 December 2002 Barnes et al.
2003/0063597 April 2003 Suzuki
2004/0066331 April 2004 Poullin
2006/0087471 April 2006 Hintz
2006/0109173 May 2006 Erikmats et al.
Зарубежные патентные документы
2101717 Oct., 1988 RU
2166769 Oct., 2001 RU

Другие ссылки

European communication mailed Mar. 9, 2010 in corresponding European Application No. 04 704 783.2--2220. cited by other .
E.K. Walton, "Future Concepts for Ground Penetrating Noise Radar", in: Ann Franchois (Ed.): PIERS Workshop on Advances in Radar Methods--A Collection of Slides, Jul. 1998, Joint Research Centre of the European Commission, Ispra (Italy), XP007912036. cited by other .
B.M. Horton: "Noise-Modulated Distance Measuring Systems", Proceedings of the I.R.E., vol. 46, No. 5, May 1, 1959, pp. 821-828, XP011156578, ISSN: 0096-8390. cited by other .
Albrecht Moschwitzer: "Formeln der Electrotechnik und Elektronik", 1989, VEB Verlag Technik, Berlin, ISBN: 3341004807, pp. 242-243, XP007912037. cited by other .
International Search Report for PCT/SE2004/000077 dated Aug. 26, 2004. cited by other .
U.S. Patent Application to Falk, U.S. Appl. No. 10/557,466, filed Nov. 18, 2005. cited by other .
Translation of Decision on Grant mailed Feb. 29, 2008 in corresponding Russian Application No. 2006119635/09(021336). cited by other .
International Search Report for PCT/SE03/00833 dated Dec. 8, 2003. cited by other .
International Search Report for PCT/SE2003/001616 dated Feb. 18, 2004. cited by other .
Strifors et al., Scattering of short EM-pulses by simple and complex targets in the combined time-frequency domain using impulse radar, Record of the 1993 IEEE National Radar Conference, Apr. 20-22, 1993, pp. 70-75. cited by other .
Anderson et al., Ultra-wideband beamforming in sparse arrays, Microwaves, Antennas and Propagation, IEEE Proceedings H, vol. 138, No. 4, Aug. 1991, pp. 342-346. cited by other .
U.S. Appl. No. 10/557,768, filed Nov. 21, 2005, Inventor: Kent Olof Faulk. cited by other .
Office Action mailed Jun. 21, 2007 in co-pending U.S. Appl. No. 10/557,768. cited by other .
Office Action mailed Dec. 7, 2007 in co-pending U.S. Appl. No. 10/557,768. cited by other .
Office Action mailed Jul. 2, 2008 in co-pending U.S. Appl. No. 10/557,768. cited by other .
Office Action mailed Feb. 26, 2008 in co-pending U.S. Appl. No. 10/557,466. cited by other .
Office Action mailed Jul. 17, 2008 in co-pending U.S. Appl. No. 10/577,466. cited by other.

Primary Examiner: Tarcza; Thomas H
Assistant Examiner: Brainard; Timothy A
Attorney, Agent or Firm: Nixon & Vanderhye P.C.
Parent Case Text


RELATED APPLICATIONS

This application is a continuation of commonly-assigned PCT/SE2004/000077, filed on Jan. 23, 2004 and is a continuation in part of commonly-assigned U.S. patent application Ser. No. 10/557,466, filed on Nov. 18, 2005 now U.S. Pat. No. 7,573,419. The entire contents of these applications are incorporated herein by reference.
ФОРМУЛА ИЗОБРЕТЕНИЯ


The invention claimed is:

1. A method for creating a tool for clutter filtering of broadband radar signals by obtaining a measure of an impedance of a target, wherein the method comprises the steps of: selecting a distance resolution (p(.tau.)); calculating a power spectrum (P(.omega.)) for a necessary transmit signal; estimating a target area profile (.rho.(.tau.)) by using a correlation (g(.tau.)) between a radiated signal and received signal; calculating a target impedance Z(.tau.) as a function of time using a relation .function..tau..times..rho..function..tau..rho..function..tau. ##EQU00013## where Z.sub.0 is the impedance of free space; and filtering out the target using the target impedance.

2. The method according to claim 1, further comprising: calculating the correlation (g(.tau.)) between the radiated signal and received signal as a convolution of the target area profile (.rho.(.tau.)) and the distance resolution (p(.tau.)).

3. The method according to claim 1, further comprising: recalculating the time (.tau.) into distance (R) using the fact that R=c.sub.0.tau.1/2, where c.sub.0 is the velocity of light.

4. The method according to claim 1, further comprising: visualizing an angle of the target impedance being the sum of partial reflections within a resolution cell using a complex presentation plane, wherein a dominating metallic contribution of the detected target echo results in an angle close to -.pi..

5. A radar system comprising: means for selecting a distance resolution (p(.tau.)); means for calculating a power spectrum (P(.omega.)) for a necessary transmit signal; means for estimating a target area profile (.rho.(.tau.)) by using a correlation (g(.tau.)) between a radiated signal and received signal; means for calculating a target impedance Z(.tau.) as a function of time using a relation .function..tau..times..rho..function..tau..rho..function..tau. ##EQU00014## where Z.sub.0 is the impedance of free space; and a filter for filtering out the target using the target impedance.

6. The radar system according to claim 5, further comprising: means for calculating the correlation (g(.tau.)) between the radiated signal and received signal as a convolution of the target area profile (.rho.(.tau.)) and the distance resolution (p(.tau.)).

7. The radar system according to claim 5, further comprising: means for recalculating the time (.tau.) into distance (R) using the fact that R=c.sub.0.tau.1/2, where c.sub.0 is the velocity of light.

8. The radar system according to claim 5, further comprising: means for visualizing an angle of the target impedance being the sum of partial reflections within a resolution cell using a complex presentation plane, wherein a dominating metallic contribution of the detected target echo results in an angle close to -.pi..

9. A radar system comprising: processing circuitry configured to perform the following: select a distance resolution (p(.tau.)); calculate a power spectrum (P(.omega.)) for a necessary transmit signal; estimate a target area profile (.rho.(.tau.)) by using a correlation (g(.tau.)) between a radiated signal and received signal; calculate a target impedance Z(.tau.) as a function of time using a relation .function..tau..times..rho..function..tau..rho..function..tau. ##EQU00015## where Z.sub.0 is the impedance of free space, and a filter for filtering out the target using the target impedance.

10. The radar system according to claim 9, wherein the processing circuitry is further configured to calculate the correlation (g(.tau.)) between the radiated signal and received signal as a convolution of the target area profile (.rho.(.tau.)) and the distance resolution (p(.tau.)).

11. The radar system according to claim 9, wherein the processing circuitry is further configured to recalculate the time (.tau.) into distance (R) using the fact that R=c.sub.0.tau.1/2, where c.sub.0 is the velocity of light.

12. The radar system according to claim 9, wherein the processing circuitry is further configured to visualize an angle of the target impedance being the sum of partial reflections within a resolution cell using a complex presentation plane, wherein a dominating metallic contribution of the detected target echo results in an angle close to -.pi..
ОПИСАНИЕ ИЗОБРЕТЕНИЯ


TECHNICAL FIELD

The technical field relates to a method for differentiating echoes from valid radar targets of interest from targets of no interest as well as disturbances considered as clutter.

BACKGROUND

A desire to be able to distinguish targets from clutter has always has been an objective with radar sensors. Traditionally this has been based on the fact that the velocity of a moving target will be different from background disturbances such as present clutter. Clutter echoes may, for instance, result from small gradients in a surrounding medium like air or a surrounding water surface. Therefore different Doppler filters have been developed to solve this task.

Thus, all present solutions generally seem to be based on the Doppler signal information extracted from received echo signals. However typical solutions of today utilise coded long pulses, for instance using a binary phase code or spread spectrum signals operating at large receiver bandwidths. Some of these signals resemble signal noise or clutter signals or a jamming signal source.

Such broadband signal radar signals generally do not exhibit a pronounced carrier, which means that a regular velocity filtering can not take place using ordinary Doppler methods. However, with new broadband radar systems there will be further possibilities to measure other echo parameters and to construct filters which, for instance, may directly filter out reflections from objects having metallic surfaces and separate those for example from those targets with generally non-metallic surfaces.

SUMMARY

Clutter filtering of broadband radar signals is performed by obtaining a measure of the impedance of an echo target. A distance resolution is then selected, and a power spectrum is calculated for the necessary transmitting signal. Further, a target area profile .rho.(.tau.) is estimated by utilising a correlation between the radiated signal and the received signal in form of a convolution. The target impedance thereafter is calculated as a function of time using the relation

.function..tau..times..rho..function..tau..rho..function..tau. ##EQU00001## where Z.sub.0 is the impedance of free space. Thereafter, by means of the achieved impedance characterising the echo of the target, the target can be filtered out.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 illustrates an example radar system.

FIG. 2 illustrates an example step of procedures for implementing clutter filtering.

DETAILED DESCRIPTION one skilled in the art that other embodiments may be practiced apart from the specific details disclosed below. In other instances, detailed descriptions of well-known methods, devices, techniques, etc. are omitted so as not to obscure the description with unnecessary detail. Individual functions are shown in blocks in the figures. Those skilled in the art will appreciate that the functions of those blocks may be implemented using individual hardware circuits, using software programs and data in conjunction with a suitably programmed microprocessor or general purpose computer, using applications specific integrated circuitry (ASIC), and/or using one or more digital signal processors (DSPs).

FIG. 1 illustrates an example radar system for detecting a target. The radar system includes conventional transmission, reception, and processing circuitry. As indicated above, new broadband radar systems open up possibilities for measuring also other parameters, for example an indication of impedance of the targets and then constructing filters which filter out, for instance, reflections from metallic objects (Z=0).

To simplify calculations, suppose a bandwidth limited broadband noise signal having an average value equal to zero in accordance with: s(t)=s.sub.1(t)cos(2.pi.f.sub.0t)-s.sub.Q(t)sin(2.pi.f.sub.0t) (1) where s.sub.1(t) and s.sub.Q(t) are two Gaussian noise processes and f.sub.0 is the centre frequency. Gaussian noise has a large amplitude variation, which means that an amplitude-limited transmitter will be ineffectively used. Phase modulated noise of the type s(t)=cos [.omega..sub.0t+.theta.(t)] having some suitable statistical distribution may therefore in practice be more preferable. Equation (1) may then alternatively be written: s(t)=A(t)cos [2.pi.f.sub.0t+.phi.(t)] (2)

The envelope function A(t) is Rayleigh distributed and given by: A(t)= {square root over (s.sub.1.sup.2(t)+s.sub.Q.sup.2(t))}{square root over (s.sub.1.sup.2(t)+s.sub.Q.sup.2(t))} (3)

The phase function .phi.(t) is equally distributed over the interval [-.pi., .pi.], and given by:

.PHI..function..function..function..function. ##EQU00002##

We also may represent the signal in a complex analytical form according to:

.function..function.e.pi..function.e.pi. ##EQU00003##

Suppose that the radar station is stationary and transmits an ergodic wave form s(t) with an average value of zero. Ergodic theory can best be described as the statistical and qualitative behaviour of measurable group and semi-group action on a measurable space possessing a non-negative measure. (Ergodic theory has its origins in the work of Boltzman in statistical mechanics problems where time- and space-distribution averages are equal.) This adoption simplifies the calculations, but implies still no limitation. A moving radar station having the velocity .nu. implies that time is scaled according to t(.tau.)=t.sub.0+2.nu..tau./c.sub.0. The received signal from the reflection profile .rho.(.tau.) for the distance 0 to R.sub.max can be written according to:

.function..intg..times..rho..function..zeta..function..zeta..times.d.zeta- ..function. ##EQU00004## where c.sub.0 is the velocity of light and n(t) represents additive noise and possible disturbances. The reflection profile .rho.(.tau.) lies within the unity circle in the complex plane and represents the sum of all subset reflections within the resolution cell. For a resolution cell for instance constituting a mix of air and metal the reflection profile .rho.(.tau.) will be a sum of contributions having a reflection coefficient of 0 and -1, respectively. Thus, the resulting reflection coefficient will be positioned along the line between 0 and -1.

According to the theory of ergodic processes, r(t) also becomes stationary and ergodic. Therefore time correlation may be utilised to approximate the cross correlation between a time delayed copy of the conjugate of the transmit signal s*(t-.zeta.) and the received signal r(t) according to:

.function..tau..intg..times..intg..times..rho..function..zeta..function..- zeta..function..tau..times.d.zeta..times.d.intg..times..function..function- ..tau..times.d ##EQU00005## where T represents the correlation of the integration interval. For a large T the second term in Eq. (7) will go towards zero since n(t) is not correlated to the transmit signal s(t). This gives:

.function..tau..fwdarw..infin..times..function..intg..times..rho..functio- n..zeta..function..tau..zeta..times.d.zeta..rho..function..tau..function..- tau. ##EQU00006## {circle around (x)} symbolises a convolution. Besides it is valid that:

.function..tau..fwdarw..infin..times..intg..times..function..function..ta- u..times.d ##EQU00007## where p(.tau.) is the auto correlation function of the radiated noise signal s(t). It comes from Eq. (8) that the received signal is the convolution between the reflection profile and the auto correlation function of the radiated signal. Thus the distance resolution is given by the auto correlation function p(.tau.) of the transmit signal, which depends on the bandwidth of the transmit signal and the form of its power spectrum P(.omega.). There is a simple relation between the auto correlation function p(.tau.) of the signal and its power spectrum P(.omega.), since these two expressions form a pair of Fourier transforms:

.function..omega..intg..infin..infin..times..function..tau.e.omega..tau..- times.d.tau..function..tau..pi..intg..infin..infin..times..function..omega- .e.omega..tau..times.d.omega. ##EQU00008##

Consequently it is possible to select a distance resolution p(.tau.) and calculate the power spectrum P(.omega.) for the necessary transmit signal. A convolution in the time domain then corresponds to a multiplication in the frequency domain. Therefore the target area profile .rho.(.tau.) can be estimated according to:

.rho..function..tau..times..function..function..tau..function..function..- tau..times..function..function..tau..function..omega. ##EQU00009## and the impedance as a function of time (distance) can be calculated according to:

.function..tau..times..rho..function..tau..rho..function..tau. ##EQU00010## where Z.sub.0=(.mu..sub.0/.epsilon..sub.0).sup.1/2 is the impedance of free space (120.pi. or .apprxeq.377 ohms). The time .tau. can be recalculated into distance R according to:

.tau. ##EQU00011##

As already mentioned, metallic objects are characterised by .rho.=-1 and Z=0. Thus, a target having another typical reflection coefficient and impedance can be filtered out in the same way as in the example with metallic objects. However it should be noted that the reflection coefficient of each individual resolution cell represents the sum of all partial reflections within the cell. In practice a threshold value representing a certain distance from the origin of coordinates need to be contribution will result in an angle close to -.pi..

Generally it is valid that the impedance for an object is given by its conductivity .sigma., dielectric constant .epsilon. and the permeability .mu. according to:

.function..omega..sigma..mu..omega..mu..mu..sigma..omega. ##EQU00012##

The detection of metallic objects is extremely simple since the frequency dependency of Eq. (14) vanish as Z(.omega., .sigma., .epsilon..sub.r, .mu..sub.r)=0.

FIG. 2 illustrates in flowchart form example procedures for clutter filtering broadband radar signals. The distance resolution and power spectrum are selected for transmitting a signal. The target area profile is then estimated. A target impedance is calculated, and the target is filtered out using that target impedance.

Filtering the target impedance using this approach permits echoes from a metallic target to be extracted even in a situation with a stationary target against a stationary background clutter. A standard Doppler detection would not be able to distinguish the target from the background clutter in this situation.

For further separation of an achieved target impedance from another detected impedance, other conventional filtering methods may be used, which are well known to a person skilled in the art.

Although various embodiments have been shown and described in detail, the claims are not limited to any particular embodiment or example. None of the above description should be read as implying that any particular element, step, range, or function is essential such that it must be included in the claims scope. The scope of patented subject matter is defined only by the to claims. The extent of legal protection is defined by the words recited in the allowed claims and their equivalents. No claim is intended to invoke paragraph 6 of 35 USC .sctn.112 unless the words "means for" are used.

* * * * *