# Patent

## Title
Signal detection

## Abstract

A transmitted signal and a received signal are combined and the combination is expected to determine whether or not a target signal is present. Either the transmitted signal or the received signal is combined with an auxiliary signal containing a range of frequencies corresponding to an anticipated Doppler shift, so that an output of the combined transmitted and received signal will be present only if a target signal exhibiting a Doppler shift within the anticipated range is present. The auxiliary signal preferably comprises finite-duration signal portions of different types so as to provide a substantially uniform frequency response throughout the selected range.

## Background
This is a divisional application of U.S. Patent application No. 10/112,873, filed Apr. 2, 2002 now U.S. Pat. No. 6,753,803


FIELD OF THE INVENTION
This invention relates to a method and apparatus for detecting signals. The invention is especially, but not exclusively, applicable to systems which may be subject to an unknown Doppler frequency shift affecting modulated coherent signals used for communication and ranging purposes, and particularly substantially continuous signals modulated by random or chaotic waveforms utilized by obstacle detection or collision avoidance systems.
BACKGROUND OF THE INVENTION
Such systems are designed to operate in multiuser and often hostile environments, and are intended for a wide range of applications, including automotive applications, industrial robotics and unmanned vehicle navigation.
FIG. 1 is a block diagram of a typical microwave obstacle-detection system. The system comprises a signal generator 1 that produces a substantially continuous waveform x(t) with suitable bandwidth to provide required range resolution. The waveform x(t) may be deterministic (periodic or aperiodic), chaotic or purely random.
The system also has a microwave oscillator 2 that generates a sinusoidal signal with required carrier frequency, a modulator 3 that modulates one or more of the parameters (such as amplitude, phase, or frequency) of the carrier signal with the modulating waveform x(t), a power amplifier (PA) 4 that amplifies the modulated carrier signal to a required level, a microwave transmit antenna (TA) 5 that radiates an electromagnetic wave representing the modulated carrier signal towards an obstacle 6, a microwave receive antenna (RA) 7 that receives an electromagnetic wave reflected back by the obstacle 6, an input amplifier (IA) 8 that amplifies a signal provided by the receive antenna (RA) 7, and a coherent demodulator 9 that processes jointly the reference carrier signal supplied by the oscillator 2 and a signal supplied by the input amplifier (IA) 8 to reconstruct a time-delayed replica y(t) of the modulating waveform x(t).
The modulating waveform x(t) and its time-delayed replica y(t) are then processed jointly during a specified time interval by a suitable processor 10, such as correlator, to produce an estimate of the unknown time delay that is proportional to the distance (range) between the system and the obstacle 6.
FIG. 2 shows an example of a correlation function of a synchronous random binary waveform.
When there occurs a relative movement between a ranging system and an obstacle of interest, an electromagnetic wave reflected back by the obstacle and received by a coherent system will exhibit a Doppler frequency shift. The value ωD0 of this (angular) frequency shift can be determined from:




  
    ω
    D0
  
  =
  
    
      
        2
        ⁢
        
          
        
        ⁢
        
          υ
          0
        
      
      c
    
    ⁢
    
      ω
      0
    
  




where υ0 is the radial speed (i.e., the range rate) of a relative movement between the system and the obstacle and ω0 is the (angular) carrier frequency of a transmitted electromagnetic wave having velocity c.

The signal reflected back by a moving obstacle can be expressed as:

zr(t)=γx(t−τ0)cos [(ω0+ωD0)(t−τ0)+θ]

where γ is the round-trip attenuation, τ0 is the delay corresponding to the range D, ωD0 is the Doppler frequency shift, and θ is an unknown constant phase shift. Strictly speaking, the value of τ0 cannot be constant for a nonzero Doppler frequency ωD0. However, in most practical cases it is assumed that the small changes in range D, of the order of the wavelength of the carrier frequency, cannot be discerned when determining a short-time time-delay estimate. Such assumption justifies decoupling delay and Doppler frequency measurements.

A baseband signal corresponding to the received signal zr(t) has the form

y(t)=x(t−τ0)cos(ωD0t+φ) 

where φ is an unknown phase shift.

The correlator performs the following operation




  
    
      R
      xy
    
    ⁡
    
      (
      τ
      )
    
  
  =
  
    
      1
      
        T
        0
      
    
    ⁢
    
      
        ∫
        0
        
          T
          0
        
      
      ⁢
      
        
          x
          ⁡
          
            (
            
              t
              -
              τ
            
            )
          
        
        ⁢
        
          x
          ⁡
          
            (
            
              t
              -
              
                τ
                0
              
            
            )
          
        
        ⁢
        
          cos
          ⁡
          
            (
            
              
                
                  ω
                  D0
                
                ⁢
                t
              
              +
              φ
            
            )
          
        
        ⁢
        
          
        
        ⁢
        
          d
          t
        
      
    
  




where the integral is evaluated for a plurality of hypothesized time delays τmin<τ<τmax. When the observation interval T0 is much shorter than one period (2π/ωD0) of the Doppler frequency, the value of cos(ωD0tm+φ) is almost constant during T0, so that the correlation integral can be approximated by





  
    
      R
      xy
    
    ⁡
    
      (
      τ
      )
    
  
  ≈
  
    
      1
      
        T
        0
      
    
    ⁢
    
      cos
      ⁡
      
        (
        
          
            
              ω
              DO
            
            ⁢
            
              t
              m
            
          
          +
          φ
        
        )
      
    
    ⁢
    
      
        ∫
        0
        
          T
          0
        
      
      ⁢
      
        
          x
          ⁡
          
            (
            
              t
              -
              τ
            
            )
          
        
        ⁢
        
          x
          ⁡
          
            (
            
              t
              -
              
                τ
                0
              
            
            )
          
        
        ⁢
        
          
        
        ⁢
        
          d
          t
        
      
    
  




where the time instant tm is taken in the middle of the observation interval T0.

When the correlation integral is calculated repeatedly for successive short processing intervals, each of duration T0, the sequence of observed correlation functions may be represented by the plot of FIG. 3. The rate of change (in time) of the correlation function Rxy(τ) will correspond to the Doppler frequency ωD0. The value of this frequency can be determined by applying some suitable form of spectral analysis.
FIG. 4 is a block diagram of a conventional correlator, comprising variable delay line 11, multiplier 12 and integrator 13, followed by a spectrum analyser 14. When the total number of successive short processing intervals is large enough, the frequency spectrum S(ω), observed at the output of the analyser, will exhibit a pronounced peak at the Doppler frequency ωD0.
FIG. 5 is a block diagram of a multichannel correlator that uses a tapped delay line 15 to cover the entire interval of hypothesized time delays, τmin<τ<τmax, in J steps of unit delay Δ using delay circuits 16. Multipliers 17 process the delayed signals with y(t) and provide the outputs to integrators 18. The required spectral analysis may be performed by a digital processor 19 implementing the Discrete Fourier Transform (DFT). The principle of operation is thus similar to that of FIG. 4.
The systems shown in FIGS. 4 and 5, and also other similar known systems, attempt to approximate the correlation integral with two parameters of interest, delay time τ and Doppler frequency ωD by combining a time sequence of correlation functions determined over a number of relatively short observation intervals T0. Because of this approximation, such systems can only provide suboptimal solutions to the problem of joint estimation of time delay and Doppler frequency.
Another prior art technique involves determining values of the following correlation integral




  
    
      R
      xy
    
    ⁡
    
      (
      
        τ
        ,
        
          
            ω
            D
          
          ;
          φ
        
      
      )
    
  
  =
  
    
      1
      
        T
        0
      
    
    ⁢
    
      
        ∫
        0
        
          T
          0
        
      
      ⁢
      
        
          x
          ⁡
          
            (
            
              t
              -
              τ
            
            )
          
        
        ⁢
        
          cos
          ⁡
          
            (
            
              
                
                  ω
                  D
                
                ⁢
                t
              
              +
              φ
            
            )
          
        
        ⁢
        
          y
          ⁡
          
            (
            t
            )
          
        
        ⁢
        
          
        
        ⁢
        
          d
          t
        
      
    
  




for the entire interval of hypothesized time delays

τmin<τ<τmax,

and also for the entire interval of hypothesized Doppler frequencies

ωDmin<ωD<ωDmax;

the unknown phase φ should also be varied over the interval (0, 2π). In this case, values of the correlation integral, determined for some prescribed range of argument values (τ,ωD), define a two-dimensional correlation function. The specific values of arguments τ and ωD, say τ0 and ωD0, that maximise the two-dimensional correlation function provide estimates of unknown time delay and unknown Doppler frequency.

SUMMARY OF THE INVENTION
FIG. 6 is a block diagram of a suitable system capable of determining values of the correlation integral for some specified time delays τ and Doppler frequencies ωD. The system is formed by a correlator like that of FIG. 4, including variable delay line 11, multiplier 12 and integrator 13. However, the signal x(t) is first applied to a multiplier 20 which multiplies it by a variable Doppler signal cos(ωDt+φ).
The procedure of adjusting the unknown value of phase φ can be avoided by suitably combining the values of the following two integrals:




  
    
      
        
          
            R
            xyc
          
          ⁡
          
            (
            
              τ
              ,
              
                
                  ω
                  D
                
                ;
                φ
              
            
            )
          
        
        =
        
          
            1
            
              T
              0
            
          
          ⁢
          
            
              ∫
              0
              
                T
                0
              
            
            ⁢
            
              
                y
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                x
                ⁡
                
                  (
                  
                    t
                    -
                    τ
                  
                  )
                
              
              ⁢
              
                cos
                ⁡
                
                  (
                  
                    
                      ω
                      D
                    
                    ⁢
                    t
                  
                  )
                
              
              ⁢
              
                
              
              ⁢
              
                d
                t
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  
  
    
      
        
          
            R
            xys
          
          ⁡
          
            (
            
              τ
              ,
              
                
                  ω
                  D
                
                ;
                φ
              
            
            )
          
        
        =
        
          
            1
            
              T
              0
            
          
          ⁢
          
            
              ∫
              0
              
                T
                0
              
            
            ⁢
            
              
                y
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                x
                ⁡
                
                  (
                  
                    t
                    -
                    τ
                  
                  )
                
              
              ⁢
              
                sin
                ⁡
                
                  (
                  
                    
                      ω
                      D
                    
                    ⁢
                    t
                  
                  )
                
              
              ⁢
              
                
              
              ⁢
              
                d
                t
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




evaluated for the entire interval of hypothesized time delays

τmin<τ<τmax 

and the entire interval of hypothesized Doppler frequencies

ωDmin<ωD<ωDmax.

In many industrial and automotive applications, systems capable of detecting moving obstacles and employing substantially continuous modulated microwave (or other coherent radiation), developed in accordance with the prior art, will be too complicated and also too expensive.
Furthermore, obstacle-detection systems designed for multiuser environments will preferably exploit some suitable form of random or chaotic modulation which may preclude the use of conventional Doppler signal processors constructed in accordance with the prior art.
Aspects of the present invention are set out in the accompanying claims.
In accordance with a further aspect of the invention, representations of a transmitted signal and a received signal are processed together in such a manner as to determine whether or not a target signal is present (this target signal possibly representing an object reflecting the transmitted signal). An auxiliary representation, for example formed of one or more auxiliary signals, is introduced into the process. The auxiliary representation includes frequencies distributed throughout a range which corresponds to anticipated frequency modifications (for example Doppler shifts) in the target signal, as a result of which the system has response characteristics such that a target signal exhibiting frequency modifications within this range will give rise to a significant system output, whereas frequency modifications outside this range will not.
Such an arrangement avoids the need for determining the value of the frequency modification (e.g. Doppler shift) but enables the rapid detection of any target signal incorporating a frequency modification of interest, using a simple structure. The above-mentioned arrangement is thus particularly useful in applications where it is sufficient to discriminate between stationary obstacles and obstacles moving with specified velocities, and for applications in which obstacles are placed into a number of classes according to obstacle velocity, for example zero/low, moderate/average, high/very high.
Preferred aspects of the present invention are directed to the nature of the auxiliary representation and give rise to a number of advantages as explained below.
Preferably, the auxiliary representation is in the form of finite-duration signal portions comprising multiple frequencies within the range of interest, and preferably distributed throughout this range. The duration of the signal portions is selected so as to maintain orthogonality of the components of the respective frequencies. By using finite-duration portions and selecting the shapes of these portions appropriately, it is possible to extend the system response characteristics such that a significant output is produced for not only target signal frequencies equal to those of the signal components used to construct the auxiliary representation, but also adjacent frequencies.
Preferably, each signal portion is created using a window function which results in a significant frequency response for all frequencies within the range of interest. However, it is difficult to arrange for a signal portion to produce a flat response throughout the frequency range of interest.
For this reason, in accordance with a preferred embodiment of the invention, the auxiliary representation comprises different finite-duration signal portions, each comprising components having a respective set of frequencies, and arranged such that they give rise to complementary frequency response characteristics. Accordingly, by combining the results achieved using the respective different auxiliary signal portions, it is possible to get an overall frequency response characteristic which is significantly flatter than if only one type of signal portion were used.
Preferably, each signal portion comprises components of frequencies interleaved with those of the other signal portion or portions.
An indicated by FIG. 3, at any given instant the correlation function may have a low level depending upon the phase of the Doppler shift. In order to ensure a significant output irrespective of the unknown phase, the auxiliary representation preferably comprises a first representation and a second representation, the components of respective frequencies in the first representation having a quadrature phase relationship with the corresponding frequency components in the second representation. By combining the responses achieved as a result of theses two auxiliary signals, the output is substantially unaffected by the phase of the Doppler shift.
Accordingly, in the preferred embodiment, there are at least four types of auxiliary signal portions: signal portions which differ in that they comprise frequencies interleaved with frequencies of other portions, and signal portions which differ in that their frequency components are in quadrature relationship with those of other portions.
The preferred embodiment uses four types of signal portions, but it would be possible to produce a system which uses more types, while maintaining the advantages indicated above. For example, the frequency components may be split into three or more sub-sets.
It would be desirable to synthesise digitally the auxiliary representation. In order to facilitate this, the auxiliary representation preferably exhibits small values of the peak factor, i.e., the peak value divided by its root-mean square value. Minimising this enables the waveform to be represented digitally with a minimum number of bits for a specified reconstruction error.
A small peak factor can be achieved by appropriate selection of the initial phase values of the signal components for each signal portion.
In the preferred embodiment described below, the transmitted and received representations are signals which are combined by correlation. However, other arrangements are possible. For example, the received signal may be modified (e.g. multiplied) by the auxiliary representation and then delivered to a filter matched to the characteristics of the transmitted signal. This use of matched filters, which utilise a representation of the transmitted signal rather than the signal itself, is, per se, known in the art.
In another alternative embodiment, the correlator is replaced by a suitably modified version of the time-delay discriminator described in WO-A-00/39643, the contents of which are incorporated herein.
Although the present invention is being described primarily in the context of detection of obstacles which produce a Doppler shift in a reflected signal, the invention is applicable also in other areas. For example, the invention is applicable to communication systems in which the local oscillators in the transmitter and receiver may not be accurately matched; a frequency offset between the oscillators would give rise to a frequency shift similar to a Doppler shift. The invention thus permits use of equipment which has wider tolerances, and is thus less expensive, than prior art arrangements.
The invention relates both to a method of signal processing and to apparatus operating according to such a method.
Arrangements embodying the invention will now be described by way of example with reference to the accompanying drawings.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a block diagram of a typical microwave obstacle-detection system.
FIG. 2 is a correlation function of a synchronous random binary waveform.
FIG. 3 shows the effects of different values of phase of the Doppler shift in a target signal on the correlation function of a synchronous random binary waveform.
FIG. 4 is a block diagram of a conventional correlator followed by a spectrum analyser.
FIG. 5 A block diagram of a multichannel correlator followed by a digital processor implementing the Discrete Fourier Transform.
FIG. 6 is a block diagram of a system for determining values of the correlation integral for various time delays and various Doppler frequencies.
FIG. 7 is a representation of the time-frequency plane divided into range-Doppler resolution cells.
FIG. 8 is a representation of the time-frequency plane divided into range resolution cells for a specified interval of Doppler frequencies.
FIG. 9 is a block diagram of a two-channel Doppler-tolerant signal processor in accordance with the present invention.
FIG. 10 is a block diagram of a multichannel correlator constructed in accordance with the prior art.
FIG. 11 is an example of an auxiliary representation in the form of signal portions for use in the processor of FIG. 9.
FIG. 12 schematically illustrates an example of operations performed by a channel combiner of the processor of FIG. 9.
FIG. 13 shows an alternative configuration of a two-channel Doppler-tolerant signal processor in accordance with the present invention.
FIG. 14 is a block diagram of a four-channel Doppler-tolerant signal processor in accordance with the present invention.
FIG. 15 is a block diagram of a quadrature Doppler generator for use in a Doppler-tolerant signal processor in accordance with the present invention.
FIGS. 16(a) and (b) and FIGS. 17(a) and (b) show examples of synthesized signal portions used to construct the auxiliary representation.
FIG. 18 shows the partial frequency response of a two-channel Doppler-tolerant signal processor employing the signal portions of FIGS. 16(a) and (b).
FIG. 19 shows the complementary frequency response of a two-channel Doppler-tolerant signal processor employing the signal portions of FIGS. 17(a) and (b).
FIG. 20 shows the combined frequency response of a two-channel Doppler-tolerant signal processor employing all four signal portions.
FIG. 21 shows ripples of the combined frequency response of FIG. 20.
FIG. 22 shows the combined frequency response of a two-channel Doppler-tolerant signal processor employing a squared cosine window.
FIG. 23 shows ripples of the combined frequency response of FIG. 22.
FIG. 24 shows a modified Kaiser window (solid line) and a squared cosine window (broken line).
FIG. 25 shows two sets, A and B, of auxiliary signals employed by two two-channel Doppler-tolerant signal processors operating in a tandem configuration.



DETAILED DESCRIPTION
FIG. 7 shows a representation of the time-frequency plane divided into range-Doppler resolution cells and a hypothetical obstacle occupying a single cell. Prior art arrangements are designed to establish the coordinates τ0,ωD0 of the cell.
When obstacles are to be classified into only two alternative classes: obstacles moving with specified velocities and all other obstacles, the time-frequency plane will be divided into a number of resolution cells equal to the number of required range cells. An obstacle detected in any of those resolution cells will have the range specified by that cell, but its velocity may be one of the plurality of admissible velocities, corresponding to an entire interval of Doppler frequencies ωD, say, ωDmin<ωD<ωDmax.
FIG. 8 is a representation of the time-frequency plane divided into a number of range resolution cells for a specified interval of Doppler frequencies. An obstacle located at the range cell To may have a radial velocity within an interval of velocities corresponding to the entire Doppler interval ωDmin<ωD<ωDmax.
In order to detect obstacles moving with any of the specified velocities, the two correlation integrals, given by (1) and (2), can be replaced by




  
    
      
        
          
            R
            xyc
          
          ⁡
          
            (
            
              τ
              ,
              
                
                  ω
                  D
                
                ;
                φ
              
            
            )
          
        
        =
        
          
            1
            
              T
              0
            
          
          ⁢
          
            
              ∫
              0
              
                T
                0
              
            
            ⁢
            
              
                y
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                x
                ⁡
                
                  (
                  
                    t
                    -
                    τ
                  
                  )
                
              
              ⁢
              
                
                  d
                  C
                
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                
              
              ⁢
              
                d
                t
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  
  
    
      
        
          
            R
            xys
          
          ⁡
          
            (
            
              τ
              ,
              
                
                  ω
                  D
                
                ;
                φ
              
            
            )
          
        
        =
        
          
            1
            
              T
              0
            
          
          ⁢
          
            
              ∫
              0
              
                T
                0
              
            
            ⁢
            
              
                y
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                x
                ⁡
                
                  (
                  
                    t
                    -
                    τ
                  
                  )
                
              
              ⁢
              
                
                  d
                  S
                
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                
              
              ⁢
              
                d
                t
              
            
          
        
      
    
    
      
        (
        4
        )
      
    
  




where dC(t) and dS(t) are especially designed auxiliary signals comprising a plurality of suitable frequency components.
The auxiliary signals dC(t) and dS(t) should be designed in such a way as:

    
    
        to contain all, mutually orthogonal within the interval T0, frequency components covering the entire interval of Doppler frequencies of interest;
        to contain two quadrature components at each frequency of interest, thus providing independence of unknown phase shift;
        to exhibit small values of the peak factor (i.e., the peak value of a waveform divided by its root-mean-square value) required to represent an analogue waveform digitally with a minimum number of bits for a specified value of a reconstruction error;
        to provide a signal processor with a flat frequency response over the entire interval of frequencies of interest, and also a required attenuation at all other Doppler frequencies.
    
    


According to a preferred embodiment of the present invention, the two auxiliary signals, dC(t) and dS(t), are obtained by suitably combining four finite-duration chirps, dCE(t), dSE(t), dCO(t) and dSO(t), defined by:




  
    
      
        
          
            d
            CE
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            
              w
              E
            
            ⁡
            
              (
              t
              )
            
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              K
            
            ⁢
            
              
            
            ⁢
            
              cos
              ⁡
              
                [
                
                  
                    2
                    ⁢
                    
                      π
                      ⁡
                      
                        (
                        
                          
                            2
                            ⁢
                            k
                          
                          +
                          u
                        
                        )
                      
                    
                    ⁢
                    
                      f
                      p
                    
                    ⁢
                    t
                  
                  +
                  
                    α
                    k
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  
  
    
      
        
          
            d
            SE
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            -
            
              
                w
                E
              
              ⁡
              
                (
                t
                )
              
            
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              K
            
            ⁢
            
              
            
            ⁢
            
              sin
              ⁡
              
                [
                
                  
                    2
                    ⁢
                    
                      π
                      ⁡
                      
                        (
                        
                          
                            2
                            ⁢
                            k
                          
                          +
                          u
                        
                        )
                      
                    
                    ⁢
                    
                      f
                      p
                    
                    ⁢
                    t
                  
                  +
                  
                    α
                    k
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  
  
    
      
        
          
            d
            CO
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            
              w
              O
            
            ⁡
            
              (
              t
              )
            
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              K
            
            ⁢
            
              
            
            ⁢
            
              cos
              ⁡
              
                [
                
                  
                    2
                    ⁢
                    
                      π
                      ⁡
                      
                        (
                        
                          
                            2
                            ⁢
                            k
                          
                          +
                          1
                          +
                          u
                        
                        )
                      
                    
                    ⁢
                    
                      f
                      p
                    
                    ⁢
                    t
                  
                  +
                  
                    β
                    k
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  
  
    
      
        
          
            d
            SO
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            -
            
              
                w
                O
              
              ⁡
              
                (
                t
                )
              
            
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              K
            
            ⁢
            
              
            
            ⁢
            
              sin
              ⁡
              
                [
                
                  
                    2
                    ⁢
                    
                      π
                      ⁡
                      
                        (
                        
                          
                            2
                            ⁢
                            k
                          
                          +
                          1
                          +
                          u
                        
                        )
                      
                    
                    ⁢
                    
                      f
                      p
                    
                    ⁢
                    t
                  
                  +
                  
                    β
                    k
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




where wE(t) and w0(t) are suitable window (or tapering) functions, K is the number of even or odd frequency components, u is the integer frequency shift, fp is the fundamental frequency, and αk and βk are phase angles so selected as to obtain the minimum, or at least low, value of the peak factor.

The chirp dSE(t) is a quadrature copy of the corresponding chirp dCE(t); similarly, the chirp dSO(t) is a quadrature copy of the corresponding chirp dCO(t). All harmonics of the fundamental frequency fp are shifted (by u) and split into two sets, each comprising either even or odd harmonics. Those sets are then employed to synthesize two orthogonal groups of chirps: {dCE(t), dSE(t)}, and {dCO(t), dSO(t)}.
Window functions wE(t) and w0(t) belong to a large class of known and widely used window functions, such as Blackman windows, Kaiser-Bessel windows or Dolph-Chebyshev windows. The duration Tw of a selected window is so chosen as to exploit the orthogonality of the frequency components, and the window type and shape are mainly determined by a required frequency resolution at each Doppler frequency of interest.
In order to simplify the optimisation procedure and also to obtain symmetrical chirps, the values of phases αk and βk can be restricted to only two values, 0 and π. In such a case, the four chirps, dCE(t), dSE(t), dCO(t) and dSO(t) will be synthesized as follows:




  
    
      
        
          
            d
            CE
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            
              w
              E
            
            ⁡
            
              (
              t
              )
            
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              K
            
            ⁢
            
              
            
            ⁢
            
              
                a
                k
              
              ⁢
              
                
              
              ⁢
              
                cos
                ⁡
                
                  [
                  
                    2
                    ⁢
                    
                      π
                      ⁡
                      
                        (
                        
                          
                            2
                            ⁢
                            k
                          
                          +
                          u
                        
                        )
                      
                    
                    ⁢
                    
                      f
                      p
                    
                    ⁢
                    t
                  
                  ]
                
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  
  
    
      
        
          
            d
            SE
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            -
            
              
                w
                E
              
              ⁡
              
                (
                t
                )
              
            
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              K
            
            ⁢
            
              
            
            ⁢
            
              
                a
                k
              
              ⁢
              
                
              
              ⁢
              
                sin
                ⁡
                
                  [
                  
                    2
                    ⁢
                    
                      π
                      ⁡
                      
                        (
                        
                          
                            2
                            ⁢
                            k
                          
                          +
                          u
                        
                        )
                      
                    
                    ⁢
                    
                      f
                      p
                    
                    ⁢
                    t
                  
                  ]
                
              
            
          
        
      
    
    
      
        (
        10
        )
      
    
  
  
    
      
        
          
            d
            CO
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            
              w
              O
            
            ⁡
            
              (
              t
              )
            
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              K
            
            ⁢
            
              
            
            ⁢
            
              
                b
                k
              
              ⁢
              
                
              
              ⁢
              
                cos
                ⁡
                
                  [
                  
                    2
                    ⁢
                    
                      π
                      ⁡
                      
                        (
                        
                          
                            2
                            ⁢
                            k
                          
                          +
                          1
                          +
                          u
                        
                        )
                      
                    
                    ⁢
                    
                      f
                      p
                    
                    ⁢
                    t
                  
                  ]
                
              
            
          
        
      
    
    
      
        (
        11
        )
      
    
  
  
    
      
        
          
            d
            SO
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            -
            
              
                w
                O
              
              ⁡
              
                (
                t
                )
              
            
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              K
            
            ⁢
            
              
            
            ⁢
            
              
                b
                k
              
              ⁢
              
                
              
              ⁢
              
                sin
                ⁡
                
                  [
                  
                    2
                    ⁢
                    
                      π
                      ⁡
                      
                        (
                        
                          
                            2
                            ⁢
                            k
                          
                          +
                          1
                          +
                          u
                        
                        )
                      
                    
                    ⁢
                    
                      f
                      p
                    
                    ⁢
                    t
                  
                  ]
                
              
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




where the coefficients ak and bk may assume only two values +1 or −1, and are so selected as to minimize the value of the peak factor.
A signal processor determining and combining the values of the two correlation integrals (3) and (4) will be referred to as a “Doppler-tolerant” signal processor. The resulting wideband Doppler compensation will make all obstacles of interest detectable (‘visible’) at the output of a Doppler-tolerant signal processor, while all other obstacles will be strongly attenuated.
FIG. 9 is a block diagram of a two-channel Doppler-tolerant signal processor, incorporating a quadrature Doppler generator QDG, according to the present invention. The processor utilizes two identical signal processing blocks, CRI and CRQ, either of which may be a conventional correlator or another appropriate signal processor. In the following, it will be assumed that each of the blocks CRI and CRQ is a conventional multichannel correlator.
FIG. 10 is a block diagram of a multichannel correlator, constructed in accordance with the prior art (see FIG. 5), which is an example of a suitable structure for each of blocks CRI and CRQ. The correlator provides J discrete values of the correlation function determined at the multiples of unit delay Δ.
In the two-channel Doppler-tolerant signal processor of FIG. 9, the function of Doppler compensation is implemented jointly by a quadrature Doppler generator QDG, two analogue multipliers, AMI and AMQ, and a channel combiner IQC.
A complete set of operations is performed by the two-channel signal processor during a basic processing time interval 2Tw split into two equal subintervals Tw, referred to as an even cycle and an odd cycle.
A baseband signal y(t), supplied by the coherent demodulator 9 in FIG. 1 (shown here as comprising a mixer AM for mixing the input signal with a local oscillator signal, followed by a low pass filter LPF), is applied simultaneously to inputs I1 and Q1 of the multipliers AMI and AMQ, respectively. The other two inputs, I2 and Q2, of the multipliers are driven by auxiliary signals produced by the quadrature Doppler generator QDG.
During an even cycle, the generator QDG produces simultaneously two chirps: dCE(t) at output I2, and dSE(t) at output Q2. During an odd cycle, another two chirps, dCO(t) and dSO(t), are produced simultaneously at outputs I2 and Q2, respectively. When the processing interval 2Tw is replicated, the alternate cycles contain chirps dCE(t) and dCO(t) which will thus form a periodic auxiliary signal dC(t). Similarly, a periodic auxiliary signal dS(t) will be obtained by periodic extension of combined chirps dSE(t) and dSO(t). Each of these two resulting quadrature signals, dC(t) and dS(t), contains the same frequency components covering the entire hypothesized Doppler frequency range of interest.
The generator QDG also produces a pulse train EOC to indicate the time instants of cycle transitions, but not necessarily a particular even or odd cycle.
FIG. 11 shows an example of two auxiliary quadrature signals, dC(t) and dS(t), and a pulse train EOC produced by the quadrature Doppler generator QDG.
Each correlator, CRI or CRQ, processes jointly a reference signal x(t) and a signal supplied by the corresponding multiplier, AMI or AMQ, during even and odd cycles. These cycles, each of duration Tw, are initiated by consecutive pulses of the pulse train EOC. Results presented at the output of the correlator CRI, after completing an even cycle and an odd cycle, are denoted by CE and CO, respectively. Similarly, SE and SO denote results obtained at the output of the correlator CRQ after completing an even cycle and an odd cycle, respectively
During an even cycle, the channel combiner IQC determines a representation of

RE=√{square root over (CE2+SE2)}

and during an odd cycle, IQC determines a representation of

RO=√{square root over (CO2+SO2)}

The above representations can be determined with the use of a conventional ‘Pythagoras processor’ or another suitable processor constructed according to the prior art.
In order to simplify the structure of the channel combiner IQC, the values of RE and RO can also be determined from some suitable approximations. For example,

RE≈max{|CE|,|SE|}+0.5 min{|CE|,|SE|}

RO≈max{|CO|,|SO|}+0.5 min{|CO|,|SO|}

Other approximations can also be constructed as known in the art.

Representations of RE and RO, determined for a prescribed number L of the processing intervals 2Tw, are suitably combined together, e.g., added or averaged, and then the resulting combined representation REO is provided at the output of the combiner IQC. For example, for a single observation interval, L=1 and REO=RE+RO. All operations performed by the channel combiner can be implemented either by a special-purpose digital processor or by a standard programmable digital signal processor.
For obstacle detection purposes, the combined representation REO is to be compared to a predetermined decision threshold. Obstacle detection will be declared in those delay (range) cells, where the threshold has been exceeded by the corresponding value of the combined representation REO.
FIG. 12 is a schematic illustration of an example of operations performed by the channel combiner IQC and the resulting detection decisions, which employ a decision threshold DT.
FIG. 13 shows an alternative configuration of a two-channel Doppler-tolerant signal processor with Doppler compensation. In this case, it is the reference signal x(t), and not a demodulated signal y(t), that is multiplied by the chirps: dCE(t), dSE(t), dCO(t) and dSO(t) produced by the quadrature Doppler generator QDG. Except for this modification, all operations and functions performed by the processor remain identical to those performed by the configuration shown in FIG. 9.
FIG. 14 is a block diagram of yet another configuration of a Doppler-tolerant signal processor utilizing Doppler compensation. The system comprises four channels: a pair of in-phase channels and a pair of quadrature channels, each with their own correlator CRIE, CRQE, CRIO or CRQO. In this case, the chirps: dCE(t), dSE(t), dCO(t) and dSO(t) are supplied simultaneously during a single cycle Tw. Consequently, the quadrature Doppler generator has now four signal outputs, in addition to an output for a processing interval signal PI indicative of the interval over which correlation should take place. Using additional separate channels for even-cycle chirps and odd-cycle chirps reduces by half the time required for determining the combined representation REO by a suitably modified channel combiner IQC.
It is also possible to devise another configuration of the quad-channel Doppler-tolerant signal processor processor shown in FIG. 14, by exploiting the modification presented in FIG. 13.
In some applications, it may be convenient to employ a single-channel Doppler tolerant signal processor, in which the four chirps, dCE(t), dSE(t), dCO(t) and dSO(t) will be produced and used sequentially.
All the above configurations of a Doppler-tolerant signal processor, and also other similar structures, will utilize an appropriate form of a quadrature Doppler generator QDG to provide suitable auxiliary signals for Doppler frequency compensation in accordance with the present invention.
FIG. 15 is a block diagram of a suitable quadrature Doppler generator QDG. The generator comprises an N-bit binary counter BCT, driven by a clock generator CLK, two memories, CVM and SVM, two buffer registers, CBR and SBR, two digital-to-analogue converters, CDA and SDA, and a trigger circuit TRG.
The operations performed by the quadrature Doppler generator can be summarised as follows. The most significant bit (MSB) of the output of the free-running counter BCT identifies even cycles, when MSB=0, and odd cycles, when MSB=1. The MSB and remaining (N−1) bits are used jointly to define addresses of 2N cells of each memory, CVM and SVM. The memory CVM contains digital amplitude samples of a chirp dCE(t) at addresses with MSB=0, whereas samples of a chirp dCO(t) are stored at addresses with MSB=1. Similarly, the memory SVM contains digital samples of a chirp dSE(t) at addresses with MSB=0, whereas samples of a chirp dSO(t) are stored at addresses with MSB=1. Consequently, each of the four chirps is represented by 2N−1 binary words, each comprising M bits, where the number M is chosen to attain a specified accuracy of the resulting digital representation of chirp samples.
Successive digital words appearing at the outputs of memories CVM and SVM are transferred to the corresponding buffer registers, CBR and SBR, at the time instants determined by a suitably delayed pulse train available at output 2 of the clock generator CLK. Digital words held in the buffer registers CBR and SBR are then converted respectively into analogue samples by two corresponding digital-to-analogue converters, CDA and SDA. As a result, the converter CDA produces repeatedly: 2N−1 analogue samples of chirp dCE(t) followed by 2N−1 analogue samples of chirp dCO(t). Similarly, the converter SDA produces repeatedly: 2N−1 analogue samples of chirp dSE(t) followed by 2N−1 analogue samples of chirp dSO(t).
A trigger circuit, TRG, generates a short pulse at the time instants coinciding with each transition of a binary waveform representing values of the most significant bit (MSB) of successive states of the counter BCT. A sequence of such short pulses, denoted by EOC, is used to indicate the beginning of each (even or odd) cycle. The pulse train EOC is utilized in the Doppler-tolerant signal processor, shown in FIG. 9 and FIG. 13, to define integration intervals Tw of the correlators, CRI and CRQ, and also to initiate various operations performed by the channel combiner IQC.
For a quad-channel Doppler-tolerant signal processor, the quadrature Doppler generator will be modified as follows. The number of bits of the binary counter BCT will be reduced to (N−1), and successive counter states will determine addresses of four memories containing digital samples of corresponding four chirps, dCE(t), dSE(t), dCO(t) and dSO(t). The output of each of the four memories will then drive a separate digital-to-analogue converter via a suitable buffer register. As a result, analogue representations of the four chirps will be available simultaneously at four outputs of such modified quadrature Doppler generator.
The modifications of the quadrature Doppler generator required for a single-channel operation will be obvious to those skilled in the art.
In order to illustrate some basic aspects of the present invention, and also to demonstrate its potential benefits, consider the following example.
Let the number of frequency components K=31 and the integer frequency shift u=4; hence, the minimum frequency and the maximum frequency are equal to 6fp and 67fp, respectively. Assume that each of the four chirps is to be represented by 2048 discrete samples and that wE(t)=w0(t)=w(t); assume also that ak=bk for k=1, 2, . . . , 31.
It has been found that the peak factor achieves its minimum value, when

    a1=+1 a2=+1 a3=+1 a4=+1 a5=+1 a6=+1 a7=+1
    a8=−1 a9=−1 a10=−1 a11=+1 a12=+1 a13=+1
    a14=−1 a15=−1 a16=−1 a17=+1 a18=+1 a19=−1
    a20=+1 a21=+1 a22=−1 a23=−1 a24=+1 a25=−1
    a26=−1 a27=+1 a28=−1 a29=+1 a30=−1 a31=+1.


Application of the above coefficients has resulted in the reduction of the peak factor value from 31/√(31·10.5)≈7.9 (when all components are simply added with unit coefficients) to only 1.9. It is desirable to ensure the peak factor value is less than 2.
To achieve a flat frequency response over the selected frequency range, and provide a required attenuation at all other frequencies, the discrete window function w(j) is a modified Kaiser window of the form




  
    w
    ⁡
    
      (
      j
      )
    
  
  =
  
    
      
        
          I
          0
        
        ⁡
        
          (
          
            ξ
            ⁢
            
              
                1
                -
                
                  m
                  2
                
              
            
          
          )
        
      
      -
      1
    
    
      
        
          I
          0
        
        ⁡
        
          (
          ξ
          )
        
      
      -
      1
    
  




where m=(j−1024)/1024 for j=0, 1, . . . , 2047, I0(·) is the modified Bessel function, and the shape parameter ξ is equal to 6.24.

Samples of the four discrete chirps, dCE(j), dSE(j), dCO(j) and dSO(j), are determined from



  
    
      
        d
        CE
      
      ⁡
      
        (
        j
        )
      
    
    =
    
      2
      ⁢
      
        w
        ⁡
        
          (
          j
          )
        
      
      ⁢
      
        
          ∑
          
            k
            =
            1
          
          31
        
        ⁢
        
          
        
        ⁢
        
          
            a
            k
          
          ⁢
          
            
          
          ⁢
          
            cos
            ⁡
            
              [
              
                2
                ⁢
                
                  π
                  ⁡
                  
                    (
                    
                      
                        2
                        ⁢
                        k
                      
                      +
                      4
                    
                    )
                  
                
                ⁢
                
                  j
                  /
                  1024
                
              
              ]
            
          
        
      
    
  




  
    
      
        d
        SE
      
      ⁡
      
        (
        j
        )
      
    
    =
    
      
        -
        2
      
      ⁢
      
        w
        ⁡
        
          (
          j
          )
        
      
      ⁢
      
        
          ∑
          
            k
            =
            1
          
          31
        
        ⁢
        
          
        
        ⁢
        
          
            a
            k
          
          ⁢
          
            
          
          ⁢
          
            sin
            ⁡
            
              [
              
                2
                ⁢
                
                  π
                  ⁡
                  
                    (
                    
                      
                        2
                        ⁢
                        k
                      
                      +
                      4
                    
                    )
                  
                
                ⁢
                
                  j
                  /
                  1024
                
              
              ]
            
          
        
      
    
  




  
    
      
        d
        CO
      
      ⁡
      
        (
        j
        )
      
    
    =
    
      2
      ⁢
      
        w
        ⁡
        
          (
          j
          )
        
      
      ⁢
      
        
          ∑
          
            k
            =
            1
          
          31
        
        ⁢
        
          
        
        ⁢
        
          
            a
            k
          
          ⁢
          
            
          
          ⁢
          
            cos
            ⁡
            
              [
              
                2
                ⁢
                
                  π
                  ⁡
                  
                    (
                    
                      
                        2
                        ⁢
                        k
                      
                      +
                      5
                    
                    )
                  
                
                ⁢
                
                  j
                  /
                  1024
                
              
              ]
            
          
        
      
    
  




  
    
      
        d
        SO
      
      ⁡
      
        (
        j
        )
      
    
    =
    
      
        -
        2
      
      ⁢
      
        w
        ⁡
        
          (
          j
          )
        
      
      ⁢
      
        
          ∑
          
            k
            =
            1
          
          31
        
        ⁢
        
          
        
        ⁢
        
          
            a
            k
          
          ⁢
          
            
          
          ⁢
          
            sin
            ⁡
            
              [
              
                2
                ⁢
                
                  π
                  ⁡
                  
                    (
                    
                      
                        2
                        ⁢
                        k
                      
                      +
                      5
                    
                    )
                  
                
                ⁢
                
                  j
                  /
                  1024
                
              
              ]
            
          
        
      
    
  



for j=0, 1, . . . , 2047.

The chirps dCE(j), dSE(j), dCO(j) and dSO(j), synthesized above, can be produced by a quadrature Doppler generator QDG, with parameters N=12 and M=12.
The synthesized chirps dCE(j), dSE(j), dCO(j) and dSO(j) are shown in FIG. 16 and FIG. 17. As seen, the chirps exhibit various forms of symmetry which can be exploited to reduce the memory of the quadrature Doppler generator QDG, at the expense of more complicated timing and control functions.
The chirps dCE(j) and dSE(j) are used for correlation during an even cycle to produce the results CE and SE, respectively. Those results are then combined by the channel combiner IQC as follows

RE=√{square root over (CE2+SE2)}

The resulting partial frequency response of a two-channel Doppler-tolerant signal processor is shown in FIG. 18.
Similarly, the chirps dCO(j) and dSO(j) are used for correlation during an odd cycle to produce the results CO and SO, respectively. Those results are combined by the channel combiner IQC as follows

RO=√{square root over (CO2+SO2)}

The resulting complementary frequency response of a two-channel Doppler-tolerant signal processor is shown in FIG. 19.
The two frequency response characteristics are superimposed (added) to produce the combined frequency response, shown in FIG. 20. As seen, the relative attenuation at out-of-band frequencies can exceed 50 dB.
The plot of the combined frequency response in FIG. 21 reveals that the ripples of the response do not exceed ±0.25 dB at all frequencies of interest.
FIG. 22 shows the combined frequency response when the modified Kaiser window is replaced a simple cosine window of the form




  
    w
    ⁡
    
      (
      j
      )
    
  
  =
  
    
      cos
      2
    
    ⁡
    
      (
      
        
          
            
              j
              -
              1024
            
            ⁢
            
              
            
          
          1024
        
        ⁢
        
          π
          2
        
      
      )
    
  




for j=0, 1, 2047.

The cosine window provides greater attenuation of stationary obstacles (when Doppler frequency is zero); however, the frequency response is not as flat as that corresponding to the modified Kaiser window. The shape and location of ripples are shown in FIG. 23.
In all the above configurations of a Doppler-tolerant signal processor, the received signal y(t) is multiplied by chirps with nonuniform envelopes determined by a selected window. Because of this loss of signal energy, the detection performance of the system can be significantly reduced. The examination of various windows utilized in a Doppler-tolerant signal processor shows that their shape is close to that of a squared cosine function. For example, FIG. 24 shows that a modified Kaiser window and a squared cosine window almost overlap. This observation and the trigonometric identity sin2φ+cos2φ=1 suggest that the loss due to the use of windows can be almost eliminated when a received signal is processed simultaneously by two Doppler-tolerant signal processors operating in tandem configuration. While either processor employs the same set of chirps, all chirps of a duplicated set are shifted in time by Tw/2, where Tw is the window duration.
FIG. 25 shows two sets of chirps employed by a tandem arrangement comprising two two-channel Doppler-tolerant signal processors, each arranged in accordance with FIG. 9.
A more efficient system will utilize a tandem configuration comprising two quad-channel Doppler-tolerant signal processors, each developed in accordance with FIG. 14.
This technique, whereby auxiliary signals are used to investigate a primary signal, and wherein each auxiliary signal comprises finite-duration portions overlapping portions of the other auxiliary signal, is considered independently inventive. It can be applied to other areas of signal processing beyond the specific areas disclosed herein, including but not limited to arrangements in which the auxiliary signals contain components of different frequencies. It can also be extended to cover the use of more than two trains of auxiliary signals.
The correlators of the embodiments described above could be replaced by other arrangements such as a suitably modified time-delay discriminator as described in WO-A-00/39643.

## Claims

1. A signal processing method comprising: modifying a primary signal using, at least, first and second auxiliary signals; and
investigating the modified primary signal, wherein each auxiliary signal comprises successive finite-duration shaped portions having relatively low-amplitude leading and trailing parts, the portions being interleaved with signal portions of the other auxiliary signal such that the leading parts of the portions of each signal overlap the trailing parts of the portions of the other signal.

2. The method as claimed in claim 1, wherein each shaped portion has a shape substantially similar to that of a squared cosine.

3. The method as claimed in claim 2, wherein each shaped portion is produced using a modified Kaiser window function.

4. A signal processing method comprising;
modifying a primary signal using, at least, first and second auxiliary signals; and
investigating the modified primary signal, wherein each auxiliary signal comprises successive finite-duration shaped portions having relatively low-amplitude leading and trailing parts, the portions being interleaved with, and overlapping, signal portions of the other auxiliary signal, and
wherein each finite-duration shaped portion contains multiple different predetermined frequencies.

5. The method as claimed in claim 4, wherein each auxiliary signal comprises finite-duration shaped portions containing a first set of frequencies, and intervening finite-duration shaped portions containing a second set of frequencies, the frequencies of the first set being interleaved with the frequencies of the second set.

6. The method as claimed in claim 4, wherein the primary signal is additionally modified by third and fourth auxiliary signals, the third auxiliary signal having finite-duration shaped portions which are produced simultaneously with the finite-duration shaped portions of the first auxiliary signal, and the fourth auxiliary signal containing finite-duration shaped portions which are produced simultaneously with the finite-duration shaped portions of the second auxiliary signal, and wherein the simultaneous finite-duration shaped portions of the first and third auxiliary signals contain frequency components in quadrature relationship with each other, and wherein the simultaneous finite duration shaped portions of the second and fourth auxiliary signals contain frequency components in quadrature relationship with each other.

7. The method as claimed in claim 5, wherein the primary signal is additionally modified by third and fourth auxiliary signals, the third auxiliary signal having finite-duration shaped portions which are produced simultaneously with the finite-duration shaped portions of the first auxiliary signal, and the fourth auxiliary signal containing finite-duration shaped portions which are produced simultaneously with the finite-duration shaped portions of the second auxiliary signal, and wherein the simultaneous finite-duration shaped portions of the first and third auxiliary signals contain frequency components in quadrature relationship with each other, and wherein the simultaneous finite duration shaped portions of the second and fourth auxiliary signals contain frequency components in quadrature relationship with each other.