Digital biquad filter

From Wikipedia, the free encyclopedia
Jump to navigation Jump to search

In signal processing, a digital biquad filter is a second-order recursive linear filter, containing two poles and two zeros. "Biquad" is an abbreviation of "biquadratic", which refers to the fact that in the Z domain, its transfer function is the ratio of two quadratic functions:

The coefficients are often normalized such that a0 = 1:

High-order IIR filters can be highly sensitive to quantization of their coefficients, and can easily become unstable. This is much less of a problem with first and second-order filters; therefore, higher-order filters are typically implemented as serially-cascaded biquad sections (and a first-order filter if necessary). The two poles of the biquad filter must be inside the unit circle for it to be stable. In general, this is true for all discrete filters i.e. all poles must be inside the unit circle in the Z-domain for the filter to be stable.

Implementation[edit]

Direct form 1[edit]

The most straightforward implementation is the direct form 1, which has the following difference equation:

or, if normalized:

Here the , and coefficients determine zeros, and , determine the position of the poles.

Flow graph of biquad filter in direct form 1:

Biquad filter DF-I.svg

Direct form 2[edit]

The direct form 1 implementation requires four delay registers. An equivalent circuit is the direct form 2 implementation, which requires only two delay registers:

Biquad filter DF-II.svg

The direct form 2 implementation is called the canonical form, because it uses the minimal amount of delays, adders and multipliers, yielding in the same transfer function as the direct form 1 implementation. The difference equations for direct form 2 are:

where

Transposed direct forms[edit]

Each of the two direct forms may be transposed by reversing the flow graph without altering the transfer function. Branch points are changed to summers and summers are changed to branch points.[1] These provide modified implementations that accomplish the same transfer function which can be mathematically significant in a real-world implementation where precision may be lost in state storage.

The difference equations for Transposed Direct Form 2 are:

where

and

Transposed Direct form 1[edit]

The direct form 1 Digital Biquad Direct Form 1 Untransformed is transposed into Flow diagram of Biquad filter Direct Form 1 Transformed for use in article https://en.wikipedia.org/wiki/Digital_biquad_filter

Transposed Direct form 2[edit]

The direct form 2 Flow diagram of Biquad filter Direct Form 2 for use in article https://en.wikipedia.org/wiki/Digital_biquad_filter is transposed into Flow diagram of Biquad filter Direct Form 2 for use in article https://en.wikipedia.org/wiki/Digital_biquad_filter

Quantizing Noise[edit]

When a sample of n bits is multiplied by a coefficient of m bits, the product has n+m bits. These products are typically accumulated in a DSP register, the addition of five products may need 3 overflow bits; this register is often large enough to hold n+m+3 bits. The z-1 is implemented by storing a value for one sample time; this storage register is usually n bits, the accumulator register is rounded to fit n bits, and this introduced quantizing noise.

In the direct form 1 arrangement, there is a single quantizing/rounding function Flow diagram of Biquad filter Direct Form 1 for use in article https://en.wikipedia.org/wiki/Digital_biquad_filter.

In the direct form 2 arrangement, there is a quantizing/rounding function for an intermediate value. In a cascade, the value may not need rounding between stages, but the final output may need rounding Flow diagram of Biquad filter Direct Form 2 with Quantizing for use in article https://en.wikipedia.org/wiki/Digital_biquad_filter.

Fixed point DSP usually prefers the non transposed forms and has an accumulator with a large number of bits, and is rounded when stored in main memory. Floating point DSP usually prefers the transposed form, each multiplication and potentially each addition are rounded; the additions are higher precision result, when both operands have similar magnitude.

See also[edit]

References[edit]

External links[edit]