CN1009519B - Low frequency digital notch filter - Google Patents
Low frequency digital notch filterInfo
- Publication number
- CN1009519B CN1009519B CN 86105367 CN86105367A CN1009519B CN 1009519 B CN1009519 B CN 1009519B CN 86105367 CN86105367 CN 86105367 CN 86105367 A CN86105367 A CN 86105367A CN 1009519 B CN1009519 B CN 1009519B
- Authority
- CN
- China
- Prior art keywords
- filter
- output
- input
- node
- multiplier
- 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.)
- Expired
Links
- 239000002131 composite material Substances 0.000 claims description 4
- 230000001629 suppression Effects 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims 1
- 230000004044 response Effects 0.000 description 4
- 230000003321 amplification Effects 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000000034 method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Images
Landscapes
- Filters That Use Time-Delay Elements (AREA)
- Networks Using Active Elements (AREA)
- Analogue/Digital Conversion (AREA)
- Noise Elimination (AREA)
- Picture Signal Circuits (AREA)
- Oscillators With Electromechanical Resonators (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
低频数字陷波滤波器(图1)包括全通网络滤波器(1)和T型滤波器(3),输入节点(31)与后者反馈连接,并与端输出节点(33)相连接,滤波器(1)的转换函数为A(Z)=(Z-1+K1)/(1+K1Z-1);滤波器(3)的通过量和抽头转换函数为B(Z)=[(K3+K2-K4)Z-1-1)/(1-K3Z-1)和C(Z)=K2Z-1/(1-K3Z-1),这里Z-1是单位延迟算子,K1~K4是乘法器系数,滤波器(3)可用三个乘法器(19,21和23)及一个延迟单元(25)联合实现。The low-frequency digital notch filter (Fig. 1) comprises an all-pass network filter (1) and a T-type filter (3), and the input node (31) is connected with the latter feedback, and is connected with the terminal output node (33), The transfer function of filter (1) is A(Z)=(Z -1 +K 1 )/(1+K 1 Z -1 ); the throughput and tap transfer function of filter (3) is B(Z) =[(K 3 +K 2 −K 4 )Z −1 −1)/(1−K 3 Z −1 ) and C(Z)=K 2 Z −1 /(1−K 3 Z −1 ), Here Z -1 is a unit delay operator, K 1 ~ K 4 are multiplier coefficients, and the filter (3) can be realized jointly by three multipliers (19, 21 and 23) and a delay unit (25).
Description
本发明是关于低频数字陷波滤波器以及包括这种滤波器的复合滤波器的设计的改进。This invention relates to improvements in the design of low frequency digital notch filters and complex filters including such filters.
有许多数字信号处理(DSP)的应用,其中所需的功能元件之一是抑制任何低频信号的滤波器,可能的话还抑制该信号的任何直流成分。这样一个滤波器的典型应用场合是在被处理的数字信号是一个模拟信号的数字表示,且该模拟信号包含不需要的50赫或60赫信号时。这些信号会是从邻近电源设备线路耦合于系统所产生的信号;常常需要除去这些不需的信号及存在的任何直流成分。There are many digital signal processing (DSP) applications where one of the required functional elements is a filter that rejects any low frequency signal and possibly any DC content of that signal. A typical application of such a filter is when the digital signal being processed is a digital representation of an analog signal containing unwanted 50 Hz or 60 Hz signals. These signals will be the result of line coupling into the system from adjacent power equipment; it is often necessary to remove these unwanted signals and any DC components present.
至今已有许多不同的滤波器结构被考虑来实施这样的滤波器,特别是在需要清除的信号的频率大大小于系统采样频率(约1%)的应用场合。在考虑这种结构的适用性时,必须考虑它们的性能的许多方面。譬如为了能保证简易有效的实施,希望滤波器有合理的有规则的结构,并且考虑到价格因素,复杂性要最小。另外一个必要条件是乘法部件要保持最少,最好乘法器系数的字长也为最小,还有导向提高信号字长需要的其他性能考虑。循环滤波器结构倾向于放大任何分层噪声,这样噪声总是跟随在乘法器级后面的舍位处理时产生,通常在信号字长中采用附加位来补偿这种放大,从而使滤波器噪声保持在可接受水平。还必须考虑在滤波器内部结点处的信号大小。一个具有高有效Q-因素的滤波器在其内部结点处会有比输入信号大40dB的信号,在一定的频率上,即使总的增益是1。提高信号字长还可避免限幅。Many different filter structures have been considered to implement such filters, especially in applications where the frequency of the signal to be cleaned is much smaller than the system sampling frequency (about 1%). When considering the suitability of such structures, many aspects of their performance must be considered. For example, in order to ensure simple and effective implementation, it is hoped that the filter has a reasonable and regular structure, and considering the price factor, the complexity should be minimal. Another requirement is that the multiplication components be kept to a minimum, and preferably the word length of the multiplier coefficients is also minimized, along with other performance considerations that lead to the need to increase the signal word length. Loop filter structures tend to amplify any layered noise such that the noise always follows the truncation process following the multiplier stage, usually with additional bits in the signal word length to compensate for this amplification so that the filter noise remains at an acceptable level. The magnitude of the signal at the internal nodes of the filter must also be considered. A filter with a high effective Q-factor will have signals at its internal nodes that are 40dB larger than the input signal at certain frequencies, even though the overall gain is 1. Increasing the signal word length also avoids clipping.
例如,设想这样一种应用情况,在一个采样频率为8赫的系统内50赫和60赫的信号至少需要衰减25dB,并消除所有直流分量,而使200赫的信号衰减小于0.7dB。For example, imagine an application where the 50 Hz and 60 Hz signals need to be attenuated by at least 25 dB in a system with a sampling frequency of 8 Hz, and all DC components should be eliminated, while the 200 Hz signal should be attenuated by less than 0.7 dB.
一个包括双二次方节的滤波器将是很有规则且易于实施的。然而这样一个经选择具有最佳噪声和信号增加长性能的滤波器包括多至14个9位乘法器和6个延迟单元。在内部结点处的噪声增益为0dB情况下,噪声放大约为13dB,如前所述,信号字长需增加附加3位来补偿滤波器。但应注意到存在许多替换结构,例如,用11个乘法器和5个延迟单元实施的滤波器会具有相当可观的坏信号增长特性。A filter including biquadratic sections will be very regular and easy to implement. However such a filter selected for optimum noise and signal adding length performance includes as many as 14 9-bit multipliers and 6 delay elements. With a noise gain of 0dB at the internal nodes, the noise amplification is about 13dB. As mentioned earlier, the signal word length needs to be increased by an additional 3 bits to compensate for the filter. It should be noted, however, that many alternative structures exist, eg a filter implemented with 11 multipliers and 5 delay elements would have considerable bad signal growth characteristics.
本发明的目的是提供一个具有合理通用结构的数字滤波器,但与已有的滤波器相比,结构简单,性能相当或更好。The object of the present invention is to provide a digital filter with a reasonable and general structure, but with a simple structure and equivalent or better performance than the existing filters.
根据本发明提供了一种低频数字陷波滤波器,它包括:A kind of low-frequency digital notch filter is provided according to the present invention, it comprises:
一个输入节点;an input node;
一个连接于上述节点的全通网络滤波器,包括至少一个延迟单元,至少一个具有乘法系数为K1的系数乘法器,该滤波器的转换函数A(Z)由下式给出:An all-pass network filter connected to the above nodes, including at least one delay unit, at least one coefficient multiplier with a multiplication coefficient of K 1 , the transfer function A(Z) of the filter is given by the following formula:
A(Z)=(Z-1+K1)/(1+K1Z-1)式中Z-1是单位延迟算子;A(Z)=(Z -1 +K 1 )/(1+K 1 Z -1 ) where Z -1 is a unit delay operator;
一个T型滤波器,其输入与上述全通网络滤波器的输出连接,其输出与上述输入节点呈负反馈连接,该T型滤波器包括至少一个延迟单元和至少三个系数乘法器,其乘法系数分别为K2、K3和K4。它的通过量转换函数B(Z)由下式给出:A T-shaped filter, whose input is connected to the output of the above-mentioned all-pass network filter, and whose output is connected with the above-mentioned input node in negative feedback, the T-shaped filter includes at least one delay unit and at least three coefficient multipliers, and its multiplication The coefficients are K 2 , K 3 and K 4 , respectively. Its throughput transfer function B(Z) is given by:
B(Z)=〔1+(K2·K4-K3)Z-1〕/(1-K3Z-1),并且它的输入至抽头转换函数C(Z)由下式给出:B(Z)=[1+(K 2 ·K 4 -K 3 )Z -1 ]/(1-K 3 Z -1 ), and its input-to-tap transfer function C(Z) is given by :
C(Z)=K2Z-1/(1-K3Z-1);以及C(Z)=K 2 Z -1 /(1-K 3 Z -1 ); and
一个端输出节点,它与上述T型循环滤波器的输出相连接,并连接到上述输入节点。a terminal output node connected to the output of the above-mentioned T-shaped loop filter and connected to the above-mentioned input node.
上述陷波滤波器根据应用情况可与一个直流抑制滤波器串联在一起。The above notch filter can be connected in series with a DC suppression filter according to the application.
本说明书的附图包括:The accompanying drawings in this manual include:
图1,是根据本发明实施的一个低频数字陷波滤波器的电路图;Fig. 1 is a circuit diagram of a low-frequency digital notch filter implemented according to the present invention;
图2,是与图1所示的陷波滤波器联合使用的直流抑制滤波器的电路图;Fig. 2, is the circuit diagram of the direct current suppressing filter used in combination with the notch filter shown in Fig. 1;
图3,是说明典型的增益频率响应,和使用上述图1、图2所示滤波器的组合所能达到的频率响应的曲线图。Figure 3 is a graph illustrating a typical gain frequency response and the frequency response achievable using the combination of filters shown in Figures 1 and 2 above.
现在参考上述附图,举例描述本发明的实施例。Embodiments of the present invention will now be described by way of example with reference to the above-mentioned drawings.
一个复合带阻滤波器的低频陷波滤波器部件示于图1中,该部件主要包括2个子滤波器,即一个全通网络滤波器1和一个T型滤波器3。A low-frequency notch filter component of a composite band-stop filter is shown in Figure 1, which mainly includes two sub-filters, namely an all-
全通网络滤波器1包括有一个延迟单元5和一个系数乘法器7。该滤波器的转换函数A(Z)由下式给出:The all-
A(Z)=(Z-1+K1)/(1+K1Z-1),其中K1是乘法器7的系数值。滤波器可由标准的循环结构来实现-例如,由RIA.Valerzuela和A.G.Conslantinides合写的题为“用于有效的内插法和取其十分之一的数字信号处理系统”。文章〔见IEE会刊130卷Pt.G No.6(1983.12)第232页〕中所描述的结构,所示的全通网络滤波器1在公共输入延迟单元5和乘法器7之间包括支路节点9和11,还包括连接于延迟单元5和乘法器7的输出处的输出节点13。在延迟单元5的输出和乘法器7的输出之间,以及在乘法器7的输出和延迟单元5之间,分别通过支路节点11和9有交叉连接线15和17。A(Z)=(Z −1 +K 1 )/(1+K 1 Z −1 ), where K 1 is the coefficient value of the multiplier 7 . Filters can be implemented by standard loop structures - eg, "Digital Signal Processing Systems for Efficient Interpolation and Decimation" by RIA. Valerzuela and AG Conslantinides. The structure described in the article [see IEE Transactions Volume 130 Pt.G No.6 (1983.12) page 232], the shown all-
T型滤波器3与全通网络滤波器1的输出节点13相连接,该滤波器包括3个乘法器19、21和23,其乘法系数分别为K2、K3和K4。滤波器3的转换函数B(Z)和C(Z)由下列二式表示:The T-shaped filter 3 is connected to the
B(Z)=〔(K3+K2·K4)Z-1-1〕/(1-K3Z-1)(输入至输出);和C(Z)=K2Z-1/(1-K3Z-1)(输入至抽头)。B(Z)=[(K 3 +K 2 ·K 4 )Z -1 -1]/(1-K 3 Z -1 ) (input to output); and C(Z)=K 2 Z -1 / (1-K 3 Z -1 ) (input to tap).
在如图1所示的该滤波器3的实施例中,乘法器19、支路节点27、延迟单元25和乘法器21为串联,而延迟单元25通过支路节点27被乘法器23分流,乘法器21的输出连接到输出节点29,节点29也连接到T型滤波器的输入端以便与输入信号相减。In the embodiment of this filter 3 as shown in Figure 1, the multiplier 19, the branch node 27, the
加上输入节点31和输出端节点33就完成了低频陷波滤波器部件。反馈信号从T型滤波器的延迟单元25和乘法器21之间的信号通路上一点提取,并与输入节点31处的输入相减。从节点31来的输出信号并联地供给全通网络滤波器33的输入和输出节点33,并在节点33处和T型滤波器3的输出信号相加。Adding an input node 31 and an output node 33 completes the low frequency notch filter block. The feedback signal is taken from a point on the signal path between
图1所示的低频陷波滤波器部件也包括第5个乘法器35,它连接在输出节点33的输出处。这个乘法器的乘法系数K5=1/2。The low frequency notch filter block shown in FIG. 1 also includes a fifth multiplier 35 connected at the output of output node 33 . The multiplication coefficient K 5 =1/2 of this multiplier.
复合滤波器的直流抑制滤波器部件示于图2。该滤波器部件与上述低频陷波滤波器部件串联连接。它主要包括一个全通网络滤波器1和一个输出节点37,该部件的输入信号并联地供给网络滤波器1和输出节点37。在节点37处输入信号与网络滤波器1的输出信号相减,全通网络滤波器1可以与前面参考图1描述的对应滤波器有相同的结构。但这里乘法器7的乘法系数为K6。在图2所示的直流抑制滤波器中,输出节点37与第7个乘法器39相连接。乘法器39的乘法系数K7=1/2。这一直流抑制滤波器的选择,与它的双二次部件对应物相比较具有更好的性能,且在形式上与陷波滤波器十分类似,从而有可能,更有效地实施。The DC rejection filter components of the composite filter are shown in Figure 2. This filter section is connected in series with the above-mentioned low frequency notch filter section. It basically comprises an all-
一个典型的性能说明示于图3。带阻频率为50赫,对200赫以上的频率获得一个水平响应。系数K1~K7的合适的值在下列表1中给出A typical performance illustration is shown in Figure 3. The band stop frequency is 50 Hz, and a horizontal response is obtained for frequencies above 200 Hz. Suitable values for the coefficients K 1 to K 7 are given in Table 1 below
表1Table 1
系数 值coefficient value
K1-(1-3/16)K 1 - (1-3/16)
K21/16
K31-1/128K 3 1-1/128
K41/4
K51/2
K6-(1-1/32)K 6 - (1-1/32)
K71/2
从图3中可以看出,当采用表1中的系数值时就能容易地获得所需的性能。图3说明直流抑制,低频带阻为50赫,在频率为200赫以上有一个接近水平的响应。It can be seen from Figure 3 that the required performance can be easily obtained when using the coefficient values in Table 1. Figure 3 illustrates DC rejection, with a low frequency bandstop of 50 Hz and a nearly flat response at frequencies above 200 Hz.
下列表给出图1和图2所示2部件结合的性能。The table below gives the performance of the combination of the 2 components shown in Figure 1 and Figure 2.
表2Table 2
结构: 合理有规则,用简单的全通网络作为滤波器Structure: Reasonable and regular, using a simple all-pass network as a filter
的基本构成部分。basic components of the .
组成: 7个乘法器,3个延迟单元。Composition: 7 multipliers, 3 delay units.
乘法器: 7位,不过7个乘法器中4个的乘法系数仅Multiplier: 7 bits, but the multiplication coefficient of 4 of the 7 multipliers is only
仅为以2为底数的幂,但能很简单地实施。Only powers of
噪声放大: 约10dB。Noise amplification: about 10dB.
内部节点的最大噪声增益:3dBMaximum noise gain for internal nodes: 3dB
信号字长: 需3位附加位来补偿滤波器Signal word length: 3 additional bits are required to compensate the filter
可以看出上述结构以稍差的规则结构为代价,在较普通的结构上作了许多改进。对于给定的信号字长,它设法获得了与以前研究的双二次部分滤波器有相似的性能,但只用了半数存贮单元和半数乘法器,且在许多情况下乘法系数很简单,例如是1/2等。It can be seen that the above structure provides many improvements over the more general structure at the expense of a less regular structure. For a given signal word length, it manages to obtain performance similar to that of previously studied biquadratic partial filters, but only uses half the number of storage units and half the number of multipliers, and in many cases the multiplication coefficients are simple, For example 1/2 etc.
系数K1~K7可以变化,以满足其他的应用的需要,特别是改变相对于采样频率的陷波滤波器的频率。上述结构的优点之一是许多乘法可以很容易地实现,即使用除以2的级联,甚至在系数改变以符合不同的要求时候也这样。因此,参阅图1和图2,系数K5和K7总是被选择成相等于1/2。如果不是K2和K4都为1/2的幂,以及系数K3与系数K2和K4具有下列关系K3=1-(2·K2K4)的话,通常可以安排K5和K7中的一个值为1/2。Coefficients K 1 -K 7 can be varied to meet other application requirements, in particular changing the frequency of the notch filter relative to the sampling frequency. One of the advantages of the above structure is that many multiplications can be easily implemented, even using cascades of divide-by-2, even when the coefficients are changed to meet different requirements. Therefore, referring to Figures 1 and 2, the coefficients K5 and K7 are always chosen to be equal to 1/2.
Claims (8)
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB8521378 | 1985-08-28 | ||
GB858521378A GB8521378D0 (en) | 1985-08-28 | 1985-08-28 | Digital notch filter |
GB8522643A GB2182513B (en) | 1985-08-28 | 1985-09-12 | Low frequency digital notch filter |
GB8522643 | 1985-09-28 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN86105367A CN86105367A (en) | 1987-02-25 |
CN1009519B true CN1009519B (en) | 1990-09-05 |
Family
ID=26289707
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 86105367 Expired CN1009519B (en) | 1985-08-28 | 1986-08-26 | Low frequency digital notch filter |
Country Status (5)
Country | Link |
---|---|
JP (1) | JPH0693599B2 (en) |
CN (1) | CN1009519B (en) |
BG (1) | BG46609A3 (en) |
CA (1) | CA1288483C (en) |
FI (1) | FI91468C (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015135508A1 (en) | 2014-03-13 | 2015-09-17 | Mediatek Inc. | Feedforward filter using translational filter |
CN111480306B (en) * | 2017-12-13 | 2023-06-30 | 瑞典爱立信有限公司 | Method for estimating propagation delay difference of optical link and apparatus for the same |
-
1986
- 1986-07-21 CA CA000514260A patent/CA1288483C/en not_active Expired - Fee Related
- 1986-08-22 BG BG076252A patent/BG46609A3/en unknown
- 1986-08-26 CN CN 86105367 patent/CN1009519B/en not_active Expired
- 1986-08-26 FI FI863448A patent/FI91468C/en not_active IP Right Cessation
- 1986-08-27 JP JP20118586A patent/JPH0693599B2/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
FI91468C (en) | 1994-06-27 |
JPH0693599B2 (en) | 1994-11-16 |
CN86105367A (en) | 1987-02-25 |
BG46609A3 (en) | 1990-01-15 |
FI863448A (en) | 1987-03-01 |
FI863448A0 (en) | 1986-08-26 |
FI91468B (en) | 1994-03-15 |
JPS6251819A (en) | 1987-03-06 |
CA1288483C (en) | 1991-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bruton | Low-sensitivity digital ladder filters | |
CN1014288B (en) | Adaptive digital filter including non-recursive part and recursive part | |
Ansari et al. | Efficient sampling rate alteration using recursive (IIR) digital filters | |
EP0693235B1 (en) | Decimation filter | |
JP3223188B2 (en) | How to reduce digital filter bank power consumption by reducing the number of products. | |
EP1449298A1 (en) | Time discrete filter comprising upsampling, sampling rate conversion and downsampling stages | |
Laakso et al. | Design and implementation of efficient IIR notch filters with quantization error feedback | |
CN1012408B (en) | Interpolation/Fractional Filter Structure and Its Notch Filter | |
JP2540460B2 (en) | Sampling rate change and filtering circuit | |
CN1009519B (en) | Low frequency digital notch filter | |
Djokic et al. | A new improvement to the Powell and Chau linear phase IIR filters | |
Cabezas et al. | FIR filters using interpolated prefilters and equalizers | |
GB2182513A (en) | Low frequency digital notch filter | |
GB2175167A (en) | Digital filters | |
Krukowski et al. | Almost linear-phase polyphase iir lowpass/highpass filter approach | |
Jovanovic Dolecek et al. | Computationally efficient multiplier-free FIR filter design | |
Bhattacharya et al. | Some observations leading to multiplierless implementation of linear phase FIR filters | |
Petraglia et al. | A switched-capacitor filter having very low sensitivity to capacitance ratio errors and to finite amplifier gains | |
Krukowski et al. | Polyphase filter design with reduced phase non-linearity | |
Bhattacharya et al. | Allpass structures for multiplierless realization of recursive digital filters | |
Kleinmann et al. | Efficient design of low delay IIR QMF banks for speech subband coding | |
Bhattacharaya et al. | Multiplierless realization of recursive digital filters | |
Davis et al. | Design of state-variable SC filters by means of polynomial transformation | |
CN118381486A (en) | A low resource consumption filter based on FPGA and its implementation method | |
CA1269723A (en) | Notch filter for an interpolator/decimator filter structure |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C10 | Entry into substantive examination | ||
C13 | Decision | ||
GR02 | Examined patent application | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C53 | Correction of patent for invention or patent application | ||
CB02 | Change of applicant information |
Address after: Vail County, England, England Applicant after: Plessey Semiconductors Limited Address before: Essex, England, England Applicant before: Plessey Overseas Ltd. |
|
COR | Change of bibliographic data |
Free format text: CORRECT: APPLICANT; FROM: PLESSEY OVERSEAS CO., LTD. TO: PLESSEY SEMICONDUCTORS CO., LTD. |
|
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |