GB2030410A - A Digital Time Domain Transversal Filter - Google Patents
A Digital Time Domain Transversal Filter Download PDFInfo
- Publication number
- GB2030410A GB2030410A GB7928969A GB7928969A GB2030410A GB 2030410 A GB2030410 A GB 2030410A GB 7928969 A GB7928969 A GB 7928969A GB 7928969 A GB7928969 A GB 7928969A GB 2030410 A GB2030410 A GB 2030410A
- Authority
- GB
- United Kingdom
- Prior art keywords
- samples
- signal
- output
- input signal
- multiplier circuits
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000005070 sampling Methods 0.000 claims abstract description 8
- 230000001131 transforming effect Effects 0.000 claims abstract description 4
- 238000001914 filtration Methods 0.000 abstract description 2
- 238000001228 spectrum Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 238000000034 method Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/06—Non-recursive filters
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0202—Two or more dimensional filters; Filters for complex signals
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0211—Frequency selective networks using specific transformation algorithms, e.g. WALSH functions, Fermat transforms, Mersenne transforms, polynomial transforms, Hilbert transforms
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Complex Calculations (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
The filter comprises two samplers (5 and 6), the first of which receives an input signal x(t) and the second of which receives the Hilbert transform signal @(t) of the input signal x(t), the samplers operating synchronously at the rate 1/T and respectively providing samples x(m) of the input signal x(t) and @(m) of the signal @(t), together with a calculator unit for transforming the sample pairs (x(m), @(m)) into output samples y(m) of an output signal, according to the relationship: <IMAGE> where m and k are integers and b(2k) and c(2k) are real weighting coefficients. It is convenient to use Hilbert transform when filtering single sideband signals, and the present filter makes use of both the input signal and its Hilbert transform, but without doubling the number of samples on which calculations are performed during any one sampling period. <IMAGE>
Description
SPECIFICATION
A Digital Time Domain Transversal Filter
The present invention relates to a digital time domain transversal filter for filtering electrical signals.
A digital time domain transversal filter is, in generai, constituted by a coupler receiving an input signal x(t), and providing signal samples x(m) at a rate 1/T, and by a calculator unit for transforming the samples x(m) into samples y(m) according to the relationship:
where k is an integer and a(k) are real weighting coefficients.
The non-recursive part of the calculator unit is often considered to be in the form of a circuit comprising: a delay line with intermediate taps which is connected to the output of the sampler and in which the samples x(m) progress from one tap to the next at the sampling rate 1/T; multiplier circuits which are connected to the input, to the intermediate taps and to the output of the delay line and which weight the samples available at the said input, intermediate taps and outputs with the weighting coefficients; and a summing circuit adding the weighted samples of the input signal available at the outputs of the multiplier circuits and providing samples of the output signal.
The present invention provides a digital time domain transversal filter which makes use both an input signal and of its Hilbert transform.
The present invention provides a digital filter comprising two samplers, the first of which receives an input signal x(t) and the second of which receives the Hilbert transform signal x(t) of the input signal x(t), the samplers operating synchronously at the rate 1/T and respectively providing samples x(m) of the input signal x(t) and x(m) of the signal x(t), together with a calculator unit for transforming the sample pairs (x(m), x(m)) into output samples y(m) of an output signal, according to the relationship:
where m and k are integers and b(2k) and c(2k) are real weighting coefficients.
The calculator unit may be considered to be a circuit comprising: a first delay line with intermediate taps which is connected in series with the sampler receiving the signal x(t) and in which the samples x(m) progress from one tap to the next at a rate 1/2T (i.e. half the sampling rate) with a sample stored between each tap, a second delay line identical to the first connected in series with the sampler receiving the signal Q(t), multiplier circuits connected to the inputs, the intermediate taps and the outputs of the delay lines to weight the samples available at the said inputs, intermediate taps and output using the weighting coefficients; and a summing circuit which is connected to the outputs of the multiplier circuits to sum the weighted samples and provides the output samples.
An embodiment of the invention is described by way of example with reference to the accompanying drawing, in which: Figure 1 is a circuit diagram of a non-recursive digital time domain transversal filter according to the prior art; and
Figure 2 is a non-recursive digital time domain transversal filter in accordance with the invention.
The non-recursive digital time domain transversal filter of the prior art shown in Figure 1 comprises, as is conventional, a sampler 1 which is connected at its input and which delivers samples x(m) of an input signal x(t) at regular intervals T, a delay line 2 with intermediate taps connected to the output of the sampler 1 and in which the successive samples (x(m), x(m-1),... x(m-2N)) progress from one tap to the next at the rate 11T, multiplier circuits 3 which are connected to the input, to the intermediate taps and to the output of the delay line 2 and which deliver at their outputs samples of the input signal which have been weighted with weighting coefficients (a(O), a(l),... a(2N)), and a summing circuit 4 which delivers the samples y(m) of the output signal.The samples y(m) obey the relationship:
where k is an integer and a(k) are real weighting coefficients.
When the filter is considered to be of infinite length, this relationship may be written in the form:
As is well known, this relationship shows that a time domain transversal filter of infinite length provides a correlation between the samples of x(m) of the function of x(t) and sample a(m) of the impulse response of the function. The filter is defined by the coefficients a(m) which may be used to define a function a(t) of limited spectrum [-1/2T, 1/2T]:
which is used below and which is called herein the analog impulse response of the filter. Its samples a(kT) define the filter coefficients.
Figure 2 is the circuit diagram of a non-recursive digital time domain transversal filter in accordance with the invention. This filter comprises:
two samplers 5 and 6, the first of which 5 receives an input signal x(t) and the second of which 6 receives the Hilbert transform Q(t) of the input signal x(t), the samplers operate in synchronism at the rate 1/T and the first sampler 5 provides samples x(m) and the second sampler 6 provides samples
Q(m);
a first delay line 7 with intermediate taps connected to the output of the first sampler 5 and in which the samples x(m) progress from one tap to the next at a rate 1/2T (i.e. half the sampling rate) with a sample x(m) being stored between each pair of adjacent taps;;
a second delay line 8 identical to the first and connected to the output of the sampler 6 and in which the samples Q(m) progress from one tap to the next at a rate 1/2T (i.e. half the sampling rate) with a sample Q(m) being stored between each pair of adjacent taps;;
first multiplier circuits 9 connected respectively to the input, intermediate taps and the output of the first delay line 7, and second multiplier circuits 10 are similarly connected to the second delay line 8, the first multipliers 9 delivering samples of the input signal x(t) weighted by respective weighting coefficients (b(O), b(2), b(4),... b(2N)), the second multiplier circuits 10 delivering samples of the
Hilbert transform signal Q(t) of the input signal x(t) weighted by respective weighting coefficients (c(O), c(2), c(4), . . . c(2N)); and
a summing circuit 11 connected to the output of the multiplier circuits to sum the output signals of the multiplier circuits 9 and 10 and deliver samples y(m) of the output signal.
The samples y(m) of the output signal are defined on the basis of the samples x(m) of the input signal x(t) and Q(m) of the signal Q(t), i.e. the Hilbert transform of the signal x(t), by relationship:
where k is an integer and b(2k), c(2k) are real weighting coefficients.
The delay lines 7 and 8 may be constituted by delay lines identical to the delay lines 2 of the filter shown in Figure 1 with only every other tap being used.
The coefficients of a filter in accordance with the invention may, in a manner analogous to that of the filter shown in Figure 1, be determined from the analog impulse response as defined above. To show this it is sufficient to observe for both filters, i.e. that shown in Figure 1 and that shown in Figure 2, and supposing then to be of infinite lengths will give the same output signal provided certain relationships hold between their weighting coefficients.
The relationship defining the output samples of the filter shown in Figure 1 , supposed to be of infinite length, is:
and may be written with the odd and even indices separated:
According to the sampling theorem, samples x(m) defining a real signal x' (t) of limited spectrum [--1/2T, 1/2T] unambiguously obey the relationship:
Likewise, samples Q(m) defining a real signal x'(t) of limited spectrum [--1/2T, 1/2To unambiguously, representative of the Hilbert transform of the signal x'(t), obey the relationship:
Now it can be shown that the analytic signal X'(t) corresponding to the real frequency limited signal x'(t) where: : X' (t)=x' (t) +job' (t) whose spectrum is limited to the range [0, 1/2T] is completely defined by the pairs of samples (x' (2KT), Q'(2KT)) which are none other than (x(2m), Q(2m)) and the samples x(m) and Q(m) are linked by the relationships:
To demonstrate these properties reference can be made to an article by J. Oswald entitiled "les signaux à spectre limited et leurs transformations" (Spectrum limited signals and their transformations) which appeared in the journal "Cables et transmission" year 4 July 1 950.
The relationships 2 make it possible to express the samples x(m-2p-1) as a function of the samples Q(m) as follows:
Relationship 1 then becomes:
By changing the index p in the first part of relationship 3 for an index u defined by the relationship: u=p-i the following relationship is obtained:
Considering, as before, the coefficients a(m) as samples at regular intervals T of the limited spectrum function a(t), it is possible to define coefficients â(m) as being the samples at regular intervals
T of the function â(t) (the Hilbert transform of the function a(t)) and by applying the relationships (2), the following relationship is obtained::
The relationship (4) then becomes
The expression for the output signal of the filter of Figure 2, supposing it to be of infinite length, can then be found by putting:
These relationships show that the coefficients b(2p) may be considered to be samples a(2pT) of the analog impulse response a(t) while the coefficients c(2pT) may be considered as negative versions of the samples â(2pT) of the Hilbert transform of the function a(t). They also give a simple method for determining the weighting coefficients of the filter shown in Figure 2 from those of the filter shown in
Figure 1.
The filter which has just been described with reference to Figure 2 is better adapted than the prior art filter shown in Figure 1 for cases where common information is to be drawn from a signal x(t) and its Hilbert transform, e.g. in single sideband amplitude modulation transmission systems.
Without going beyond the scope of the invention, it is possible to modify various arrangements or to replace various means by other equivalent means. In particular, the transversal filter in accordance with the invention may be used in conjunction with a conventional recursive portion.
Claims (3)
1. A digital time domain transversal filter comprising two samplers, the first of which receives an input signal x(t) and the second of which receives the Hilbert transform signal Q(t) of the input signal x(t), the samplers operating synchronously at the rate 1/T and respectively providing samples x(m) of the input signal x(t) and Q(m) of the signal Q(t), together with a calculator unit for transforming the sample pairs (x(m), x(m)) into output samples y(m) of an output signal, according to the relationship:
where m and k are integers and b(2k) and c(2k) are real weighting coefficients.
2. A filter according to claim 1, wherein the calculator unit comprises:
a first delay line with intermediate taps which is connected to the output of the first sampler and in which the samples x(m) progress from one tap to the next at a rate 1/2T (i.e. half the sampling rate) with a sample x(m) being stored in between each pair of adjacent taps;
a second delay line identical to the first, connected to the output of the second sampler and in which the samples Q(m) progress from one tap to the next at the rate 1/2T (i.e. half the sampling rate) with a sample Q(m) being stored in between each pair of adjacent taps::
first multiplier circuits connected respectively to the input, the intermediate taps and the output of the first delay line, and second multiplier circuits connected respectively to the input, the intermediate taps and the output of the second delay line, the first multiplier circuits delivering samples of the input signal x(t) weighted by respective weighting coefficients (b(O), b(2), b(4),... b(2N)) and the second multiplier circuits delivering samples of the Hilbert transform signal Q(t) of the input signal x(t) weighted by respective weighting coefficients (c(O), c(2), c(4),... c(2N)); and
a summing circuit connected to the outputs of the multiplier circuits and delivery samples y(m) of the output signal.
3. A digital time domain transversal filter substantially as herein described with reference to
Figure 2 of the accompanying drawings.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR7825839A FR2435858A1 (en) | 1978-09-08 | 1978-09-08 | TIME DOMAIN TRANSVERSE DIGITAL FILTER |
Publications (1)
Publication Number | Publication Date |
---|---|
GB2030410A true GB2030410A (en) | 1980-04-02 |
Family
ID=9212445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
GB7928969A Withdrawn GB2030410A (en) | 1978-09-08 | 1979-08-20 | A Digital Time Domain Transversal Filter |
Country Status (7)
Country | Link |
---|---|
BE (1) | BE878383A (en) |
DE (1) | DE2933003A1 (en) |
FR (1) | FR2435858A1 (en) |
GB (1) | GB2030410A (en) |
IT (1) | IT1119359B (en) |
LU (1) | LU81648A1 (en) |
NL (1) | NL7906619A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0290336A1 (en) * | 1987-05-06 | 1988-11-09 | Societe De Prospection Electrique Schlumberger | Seismic signal velocity filtering |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3144456A1 (en) * | 1981-11-09 | 1983-05-19 | Siemens AG, 1000 Berlin und 8000 München | Transversal filter for converting digital signals |
-
1978
- 1978-09-08 FR FR7825839A patent/FR2435858A1/en active Granted
-
1979
- 1979-08-14 DE DE19792933003 patent/DE2933003A1/en not_active Withdrawn
- 1979-08-20 GB GB7928969A patent/GB2030410A/en not_active Withdrawn
- 1979-08-23 BE BE1/9501A patent/BE878383A/en unknown
- 1979-09-04 NL NL7906619A patent/NL7906619A/en not_active Application Discontinuation
- 1979-09-05 LU LU81648A patent/LU81648A1/en unknown
- 1979-09-07 IT IT68785/79A patent/IT1119359B/en active
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0290336A1 (en) * | 1987-05-06 | 1988-11-09 | Societe De Prospection Electrique Schlumberger | Seismic signal velocity filtering |
FR2614995A1 (en) * | 1987-05-06 | 1988-11-10 | Schlumberger Prospection | METHOD FOR SPEED FILING SEISMIC SIGNALS AND INSTALLATION FOR ITS IMPLEMENTATION |
US4852068A (en) * | 1987-05-06 | 1989-07-25 | Schlumberger Technology Corporation | Seismic signal velocity filtering |
Also Published As
Publication number | Publication date |
---|---|
BE878383A (en) | 1980-02-25 |
FR2435858A1 (en) | 1980-04-04 |
IT1119359B (en) | 1986-03-10 |
NL7906619A (en) | 1980-03-11 |
FR2435858B1 (en) | 1981-01-30 |
LU81648A1 (en) | 1980-04-21 |
IT7968785A0 (en) | 1979-09-07 |
DE2933003A1 (en) | 1980-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Allen | Short term spectral analysis, synthesis, and modification by discrete Fourier transform | |
Rabiner et al. | The chirp z‐transform algorithm and its application | |
Diniz et al. | Digital signal processing: system analysis and design | |
Claasen et al. | On stationary linear time-varying systems | |
Gallagher et al. | A theoretical analysis of the properties of median filters | |
Erseghe et al. | Unified fractional Fourier transform and sampling theorem | |
Chen | One-dimensional digital signal processing | |
Kiya et al. | A development of symmetric extension method for subband image coding | |
Gazi | Understanding digital signal processing | |
DE69012164T2 (en) | Digital fir filter for high speed communication systems. | |
US3883726A (en) | Redundant fast fourier transform data handling computer | |
EP0180989B1 (en) | Rate conversion digital filter | |
US5265042A (en) | Device and method for a nonlinear comb filter | |
GB2030410A (en) | A Digital Time Domain Transversal Filter | |
EP0123278B1 (en) | Digital wave filter with bridged digital wave filters | |
EP0791242B1 (en) | Improved digital filter | |
Huang et al. | Frequency-domain considerations of LSV digital filters | |
Rabiner | Digital techniques for changing the sampling rate of a signal | |
Angus | Modern Sampling: A Tutorial | |
Butzer et al. | Herbert Raabe's work in multiplex signal transmission and his development of sampling methods | |
Ismail | New z-domain continued fraction expansions | |
EP0258574A2 (en) | Filter device | |
Petraglia et al. | On the use of filter banks in adaptive filtering | |
Matiu-Iovan et al. | A cubic B-spline interpolation algorithm that uses the first derivative values of the input function in the knots | |
Kittel | New general approach to commensurate TEM transmission line networks using state space techniques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WAP | Application withdrawn, taken to be withdrawn or refused ** after publication under section 16(1) |