A digital pulse compression filter for use in a radar or sonar transmitting and receiving unit to process sampled and digitized signals. The sampling frequency (fs) is equal to about four times the center frequency (fo) and greater than about twice the bandwidth (Δf) of these signals. The pulse compression filter (5) comprises a time correlation circuit (7) supplied on the one hand with the sampled and digitized signals and on the other hand with signals representing a replica of the radar or sonar transmitter pulse, whereby one of the two types of signals supplied to the correlation circuit (7) is correlated with the orthogonal components of the other type of signals supplied to the correlation circuit (7) to obtain the orthogonal components of the compressed pulse.

Patent
   4673941
Priority
Dec 07 1983
Filed
Nov 19 1984
Issued
Jun 16 1987
Expiry
Nov 19 2004
Assg.orig
Entity
Large
11
8
EXPIRED
1. A digital pulse compression filter for a radar or sonar transmitting and receiving unit arranged for the generation and transmission of frequency-modulated transmitter pulses, the reception of echo signals, and the conversion of these signals into sampled and digitized signals, the sampling frequency (fs) thereof being equal to about four times the center frequency (fo) of the sampled and digitized signals and greater than about twice the bandwidth (Δf) of said signals, characterized in that the pulse compression filter comprises a time correlation circuit supplied with the sampled and digitized signals and with signals representing a replica of the transmitter pulse, one of the two types of signals supplied to the correlation circuit being correlated with the orthogonal components of the other type of signals supplied to the correlation circuit to obtain the orthogonal components of the compressed pulse.
2. A digital pulse compression filter as in claim 1, characterized in that the correlation circuit comprises a non-recursive FIR filter for combining said sampled and digitized signals with two orthogonal pulses which represent the replica of the transmitter pulse, said non-recursive FIR filter including a delay circuit having N elements and having N-1 tappings coupled to corresponding multipliers for multiplying the supplied signals alternately with a value determined by one component of the replica of the transmitter pulse and a value determined by the other component of the replica of the transmitter pulse, the summed output signals of the multipliers alternately representing one of the orthogonal components of the compressed pulse.
3. A digital pulse compression filter as in claim 2, characterized in that the tappings k and N-k, where k=0, 1, 2, . . . 1/2N-1, with increasing values of k are alternately connected to inputs of respective adders and subtractors, the output of each adder and subtractor being connected to the input of a corresponding multiplier, the tapping 1/2N being directly connected to a corresponding multiplier.
4. A digital pulse compression filter as in claim 2, characterized in that the delay circuit comprises first and second simultaneously active parts, each comprising 1/2/N elements, the second part comprising an inverter circuit inserted between the 1/4Nth and (1/4N+1)th element, the tappings k and 1/2N-k of the two parts, where k=0, 1, 2, . . . , 1/4N-1, being connected to inputs of a respective adder, an output of each adder being coupled to an input of a corresponding multiplier, the tapping 1/4N of the first part being directly connected to a corresponding multiplier.

The invention relates to a digital pulse compression filter for a radar or sonar transmitting and receiving unit arranged for the generation and transmission of frequency-modulated transmitter pulses, the reception of echo signals, and the conversion of these signals into sampled and digitised signals.

Various embodiments of digital pulse compression filters are known. U.S. Pat. No. 3,680,105 describes a pulse compression filter utilising a frequency correlation circuit, multiplying the complex conjugate of the discrete Fourier transform of signals representing a replica of the transmitter pulse by the discrete Fourier transform of the sampled and digitised echo signals to obtain the compressed pulse from the obtained product signal through an inverse discrete Fourier transformation. The orthogonal components of the echo signal required for the discrete Fourier transformation are obtainable by quadrature detection in a conventional way. The orthogonal components of the signals representing the replica of the transmitter pulse need be calculated once only. U.S. Pat. No. 4,379,295 describes a digital pulse compression filter, where the orthogonal components of the echo signals obtained by quadrature detection are sampled, digitised and supplied to an FFT unit; the signals supplied over the various frequency output channels are each subjected to a specific delay and combined to form a compressed pulse. These pulse compression filters are however very complicated, while the latter pulse compression filter does not admit of a proper adaptation to the form of the transmitter pulse. The present invention has for its object to obviate these disadvantages to a high degree.

According to the invention, the sampling frequency (fs) is equal to about four times the centre frequency (fo) of the sampled and digitised signals and greater than about twice the bandwidth (Δf) of the signals, and the pulse compression filter comprises a time correlation circuit supplied on the one hand with the sampled and digitised signals and on the other hand with signals representing a replica of the transmitter pulse, whereby one of the two types of signals supplied to the correlation circuit is correlated with the orthogonal components of the other type of signals supplied to the correlation circuit to obtain the orthogonal components of the compressed pulse. Hence, the supplied single, sampled and digitised echo signals can be correlated with the orthogonal components of the signals representing the replica of the transmitter pulse, and the orthogonal components of the supplied sampled and digitised echo signals to be obtained in the pulse compression filter with the single signals representing the transmitter pulse.

In the latter case, the sampled and digitised echo signals have to provide the orthogonal components thereof. To this effect a Hilbert filter forming part of the pulse compression filter can be used. The Hilbert filter adds an orthogonal component of these signals to the sampled and digitised echo signals. The orthogonality of the thus obtained signals however does not meet the desired accuracy unless this filter, which is designed as a non-recursive FIR Finite Impulse Response filter, comprises a delay circuit of an impractical length. Not until then the ripple value in the passband and that in the stop band of the Hilbert filter, which values are coupled together here, will be sufficiently small, in order that the added component and the ripple-free component already present are orthogonal to some extent. It is better to add the sampled and digitised signals to two filters forming part of the pulse compression filter. These two filters can be obtained as follows: From a non-recursive FIR filter having a low-pass filter characteristic, of which the pulse reponses h(n) are obtained in accordance with the "Remez Exchange Algorithm", see L. R. Rabiner and G. Gold "Theory and Application of Digital Signal Processing" (Prentice Hall, Inc., 1975), § 3.30 and appendix, pp 187-204. Let the pass band of this filter range from (k-1/4)fs to (k+1/4)fs and the stop band from (k+1/4)fs to (k+3/4)fs, where k= . . . , -2, -1, 0, 1, 2, . . . . By multiplying the pulse responses by a factor e1/2jnπ, so that h'(n)=h(n)·e1/2jnπ, the frequency characteristic is shifted in such a way that the pass band ranges from kfs to (k+1/2)fs and the stop band from (k+1/2)fs to (k+1)fs. Of the new pulse responses h'(n) the values for n is even are orthogonal with those for n is odd. The even and odd pulse responses then determine the two filters. The pulse responses of one filter are thus . . . , h'(-2), 0, h'(0), 0, h'(2), 0, h'(4), . . . ; and those of the other filter . . . , h'(-3), 0, h'(-1), 0, h'(1), 0, h'(3), . . . . A more favourable embodiment is obtained by integrating the two filters with an orthogonal pulse response in the correlation circuit. In such a case, the correlation circuit is formed by two non-recursive FIR filters, each of which filters receiving the sampled and digitised signals and correlating these signals with one of the two orthogonal components of the replica of the transmitter pulse. The two filters forming the correlation circuit may be combined to form one single non-recursive FIR filter.

The invention will now be explained with reference to the accompanying figures, of which:

FIG. 1 shows a part of a radar receiving unit having a digital pulse compression filter;

FIG. 2 illustrates a frequency characteristic of the IF-transformed and then sampled FM echo signals;

FIGS. 3 and 4 are two embodiments of the digital pulse compression filter according to the invention;

FIGS. 5A-5F are diagrams useful in explaining the operation of the embodiments shown in FIGS. 3 and 4; and

FIGS. 6-8 are three detailed embodiments of the digital pulse compression filter according to the invention.

A part of a radar receiving unit is illustrated in FIG. 1. The frequency-modulated, IF-transformed echo signals are obtained from the IF amplifier 1 and the band filter 2 connected thereto. The centre frequency of these signals will then be fIF and their bandwidth Δf. In circuit 3 the IF signals are sampled at a frequency fs and then digitised in circuit 4. The sampling frequency is so chosen that a frequency conversion is performed resulting in sampled and digitised signals having a centre frequency fo =1/4fs. For bandwidth Δf of these signals: Δf≦1/2fs. The signals from the analogue-digital converter 4 are supplied to the digital pulse compression filter 5. Filter 5 performs the time correlation with signal R, representing a replica of the transmitter pulse, and supplies the orthogonal components I and Q of the compressed pulse. FIG. 2 shows a frequency characteristic of the sampled and digitised signals; the fundamental frequencies obtained after the sampling of the intermediate frequency signals lie in the frequency interval [0, 1/2fs ].

FIG. 3 is an embodiment of the pulse compression filter 5, here comprising a Hilbert filter 6 and a time correlator 7. The Hilbert filter adds an orthogonal component to the input signal of filter 5 to correlate the obtained components with the replica R of the transmitter pulse. FIG. 4 shows an embodiment of the pulse compression filter, comprising two non-recursive FIR filters with a mutual orthogonal pulse response 8 and 9, and the time correlator 7. The orthogonal components obtained with the two filters are in turn correlated with signal R. The results, obtained with the use of a Hilbert filter on the one side and with the use of two filters with an orthogonal pulse response on the other side, can be compared with reference to FIGS. 5A-5F.

FIG. 5A shows the ideal amplitude frequency characteristic 10 in the interval [0, fs ] of a digital filter, used to obtain from a signal f(t) the orthogonal component g(t) thereof, i.e. the Hilbert transform of f(t), where g(t)=f(t)*(-1/πt). FIG. 5B illustrates the ideal characteristic for the f(t)+jg(t) signal. The ideal Hilbert filter is of infinite length; the filter used in practice, however, is a non-recursive FIR filter having a limited delay circuit. In such a case, the amplitude-frequency characteristic of the filter is as shown in FIG. 5C, by the numeral 12. In the two intervals [0, 1/2fs ] and [1/2fs, fs ] the ripple is equal to and dependent on the length of the delay circuit applied in the filter. The characteristic for the f(t)+jg(t) signal is then as shown in FIG. 5D by the numeral 13. In the pass and stop bands the ripple is the same. The greater the ripple, the less the orthogonality of the f(t) and g(t) signals with respect to each other. A better orthogonality of the f(t) and g(t) can be obtained by using two filters with orthogonal pulse responses. FIG. 5E illustrates the amplitude-frequency characteristic 14 and 15, of these filters. Since these filters also consist of a non-recursive FIR filter, they also show a ripple depending on the length of the delay circuit; the ripples are however adapted to each other to a high degree, at any rate around the frequency values 1/4fs and 3/4fs and in such a way that the characteristic of the f(t)+jg(t) signal, as indicated in FIG. 5E by the numeral 16, shows a practically rippleless stop band, implying a great mutual orthogonality between the f(t) and g(t) signals.

The embodiments of FIGS. 6-8 all represent a pulse compression filter, where the input signal is correlated with the two orthogonal components IR, QR of the replica R of the transmitter pulse. Like the replica itself, the IR and QR components can be established by one calculation only, and therefore supplied separately to the pulse compression filter, as shown in the embodiment of FIG. 6, as well as entered in memory means forming part of the pulse compression filter, as indicated by the embodiments of FIGS. 7 and 8. Separate filters for obtaining the orthogonal components of signal R are not required in these embodiments; the pulse compression filter here consists only of a time-correlation circuit. Although this circuit may consist of two separate non-recursive FIR filters it is of greater advantage to combine these two filters to form one single non-recursive FIR filter and to accomplish the working of two filters by correct timing. Therefore, in the embodiments depicted in FIGS. 6-8 the pulse compression filter consists of one single non-recursive FIR filter.

The pulse compression filter of FIG. 6 comprises a delay circuit consisting of N elements 17. Each of the N+1 tappings of this delay circuit is connected to a corresponding multiplying circuit 18. Each of the multiplying circuits 18 is supplied with a certain digital number through the switch 19. These numbers represent either the value of the IR component, or the value of the QR component. Switch 19 is operated at a frequency twice as large as the shifting frequency of the signals through the delay circuit; i.e the pulse repetition frequency of the radar or sonar transmitting and receiving unit. This is to achieve that all signals of the tappings are multiplied both by a series of numbers, representing the value of the IR component, and by a series of numbers representing the value of the QR component. In this way, two series of product values are produced in a time corresponding with the delay time of one of the elements 17, i.e. the shift period. From the two series of product values the adder 20 produces two sum values. Each of the sum values appearing after a half and after a full shift period are alternately fed to the registers 21 and 22. The values loaded in these registers represent the two orthogonal components I and Q of the compressed pulse.

The frequency-modulated signals to be sampled, digitised, and fed to the pulse compression circuit can be expressed by the relationship:

A(t)=A0 sin[2π∫(fo -1/2Δf+Δft/T)dt],

where a linear relation between frequency and time is assumed and where, exclusive of the values fo and Δf already defined, A0 is a constant and T the pulse duration. Since the pulse is sampled at a frequency fs =4fo and the total number of samples per pulse equals N+1, while β=(Δf/2fo), t=(g/4fo, and T=(N/4fo), the amplitude of the signal sample A(k) can be expressed by:

A(k)=A0 sin 2π{1/4k-1/4βk·(N-k/N}

and, hence, that of signal sample A(N=k) by:

A(N-k)=A0 sin 2π{1/4N-1/4k-1/4βk·(N-k/N).

From these two relationships it follows that, provided N is an even number,

A(N-K)=(-1)1/2N+k A(k).

If now N is chosen to be divisible by 4, then:

A(N-k)=(-1)k A(k).

A similar relationship is obtained when the frequency modulation is not linear in time; in such a case, the deviation from the linearity should however be skew-symmetric with respect to the frequency fo. The latter relationship offers the possibility to simplify the embodiment of the pulse compression filter in FIG. 6 to the embodiment shown in FIG. 7. The embodiment of FIG. 7 comprises a delay circuit consisting of N elements. The tappings k and N-k, where k=0, 1, 2, . . . , 1/2N-1, with increasing values of k, are alternately connected to the inputs of an adder 23 and subtractor 24, respectively. The outputs of adder 23 and subtractor 24 are connected to multipliers 25. The tapping 1/2N is directly coupled to a corresponding multiplier. The numbers forming the IR and QR signals to multiply the output signals of the adder and the subtractor are here already incorporated permanently in the multipliers, or the multipliers contain memory means for storing these numbers. The number of multipliers is then practically halved with respect to the number in the embodiment of FIG. 6. The operation of the embodiment of FIG. 7 is otherwise identical to that in FIG. 6.

A further simplification is provided in the embodiment of FIG. 8. Here the delay circuit is built up of two simultaneously active parts, comprising delay elements 26 and 27. The number of delay elements in each of the parts is 1/2N. With circuit 28 it is achieved that the added series of signal samples is divided into two parallel series of signal samples to be shifted, viz. a series A0, A2, A4, . . . and a series A1, A3, A5, . . . , where A(N-k)=(-1) A(k). This division into two separate series of signal samples enables to supply the signals from the tappings of the two delay circuit parts to adders 29 only, as distinct from the embodiment of FIG. 7. It is true that this requires the insertion of an inverter circuit 30 between the 1/4Nth and (1/4N+1)th delay elements. Another consequence of the division of the signal sample series into two separate series is that the frequency at which the I and Q components are obtained can be halved. The operation of this embodiment is otherwise the same as that of the embodiments of FIGS. 6 and 7.

VAN DER Mark, Jacobus

Patent Priority Assignee Title
10677901, Apr 08 2016 Samsung Medison Co., Ltd.; SOGANG UNIVERSITY RESEARCH FOUNDATION Ultrasound diagnosis apparatus and method of operating the ultrasound diagnosis apparatus
4768156, May 06 1986 TULSA DENTAL PRODUCTS INC Imaging system
5583826, May 29 1992 CGGVERITAS SERVICES SA Method of processing signals for geophysical prospecting, the method making use of an improved wave field extrapolation operator
7443337, Sep 25 2006 AAI Corporation Synthesizing arbitrary waveforms using convolution processors
9280315, Oct 27 2011 Intel Corporation Vector processor having instruction set with vector convolution function for fir filtering
9362977, Oct 27 2011 Intel Corporation Incremental preamble detection
9363068, Apr 17 2013 Intel Corporation Vector processor having instruction set with sliding window non-linear convolutional function
9507017, Mar 15 2013 SRC, INC Simultaneous multi-frequency signal processing method
9923595, Apr 17 2013 Intel Corporation Digital predistortion for dual-band power amplifiers
9935761, Apr 17 2013 Intel Corporation Modeling of a target volterra series using an orthogonal parallel wiener decomposition
9960900, Apr 17 2013 Intel Corporation Modeling of a physical system using two-dimensional look-up table and linear interpolation
Patent Priority Assignee Title
3156914,
3217324,
4028700, Nov 16 1972 Raytheon Company Pulse compression radar and method of operating such a radar
4156876, Dec 16 1976 International Standard Electric Corporation Autocorrelation sidelobe suppression device for a continuous periodic phase-coded signal
4416016, Jun 11 1981 Bell Telephone Laboratories, Incorporated Differential phase shift keyed receiver
4556010, Mar 20 1981 Method of controlling noxious insects
4591857, Jul 11 1983 The United States of America as represented by the Secretary of the Air Programmable LFM signal processor
GB1046318,
//
Executed onAssignorAssigneeConveyanceFrameReelDoc
Nov 08 1984VAN DER MARK, JACOBUSHOLLANDSE SIGNAALAPPARATEN B V ASSIGNMENT OF ASSIGNORS INTEREST 0043470064 pdf
Nov 19 1984Hollandse Signaalapparaten B.V.(assignment on the face of the patent)
Date Maintenance Fee Events
Dec 03 1990M173: Payment of Maintenance Fee, 4th Year, PL 97-247.
Jan 09 1991ASPN: Payor Number Assigned.
Sep 23 1994M184: Payment of Maintenance Fee, 8th Year, Large Entity.
Oct 05 1994ASPN: Payor Number Assigned.
Oct 05 1994RMPN: Payer Number De-assigned.
Jan 05 1999REM: Maintenance Fee Reminder Mailed.
Jun 13 1999EXP: Patent Expired for Failure to Pay Maintenance Fees.


Date Maintenance Schedule
Jun 16 19904 years fee payment window open
Dec 16 19906 months grace period start (w surcharge)
Jun 16 1991patent expiry (for year 4)
Jun 16 19932 years to revive unintentionally abandoned end. (for year 4)
Jun 16 19948 years fee payment window open
Dec 16 19946 months grace period start (w surcharge)
Jun 16 1995patent expiry (for year 8)
Jun 16 19972 years to revive unintentionally abandoned end. (for year 8)
Jun 16 199812 years fee payment window open
Dec 16 19986 months grace period start (w surcharge)
Jun 16 1999patent expiry (for year 12)
Jun 16 20012 years to revive unintentionally abandoned end. (for year 12)