Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Infinite impulse response

From Wikipedia, the free encyclopedia
Property of many linear time-invariant (LTI) systems
icon
This articleneeds additional citations forverification. Please helpimprove this article byadding citations to reliable sources. Unsourced material may be challenged and removed.
Find sources: "Infinite impulse response" – news ·newspapers ·books ·scholar ·JSTOR
(April 2015) (Learn how and when to remove this message)

Infinite impulse response (IIR) is a fundamental property applying to manylinear time-invariant systems that are distinguished by having animpulse responseh(t){\displaystyle h(t)} that does not become exactly zero past a certain point but continues indefinitely.[1] This is in contrast to afinite impulse response (FIR) system, in which the impulse responsedoes become exactly zero at timest>T{\displaystyle t>T} for some finiteT{\displaystyle T}, thus being of finite duration. Common examples of linear time-invariant systems are mostelectronic anddigital filters. Systems with this property are known asIIR systems orIIR filters.

In practice, the impulse response, even of IIR systems, usually approaches zero and can be neglected past a certain point. However the physical systems which give rise to IIR or FIR responses are dissimilar, and therein lies the importance of the distinction. For instance, analog electronic filters composed of resistors, capacitors, and/or inductors (and perhaps linear amplifiers) are generally IIR filters. On the other hand,discrete-time filters (usually digital filters) based on a tapped delay lineemploying no feedback are necessarily FIR filters. The capacitors (or inductors) in the analog filter have a "memory" and their internal state never completely relaxes following an impulse (assuming the classical model of capacitors and inductors where quantum effects are ignored). But in the latter case, after an impulse has reached the end of the tapped delay line, the system has no further memory of that impulse and has returned to its initial state; its impulse response beyond that point is exactly zero.

Implementation and design

[edit]

Although almost allanalog electronic filters are IIR, digital filters may be either IIR or FIR. The presence of feedback in the topology of a discrete-time filter (such as the block diagram shown below) generally creates an IIR response. Thez domaintransfer function of an IIR filter contains a non-trivial denominator, describing those feedback terms. The transfer function of an FIR filter, on the other hand, has only a numerator as expressed in the general form derived below. All of theai{\displaystyle a_{i}} coefficients withi>0{\displaystyle i>0} (feedback terms) are zero and the filter has no finitepoles.

The transfer functions pertaining to IIR analog electronic filters have been extensively studied and optimized for their amplitude and phase characteristics. These continuous-time filter functions are described in theLaplace domain. Desired solutions can be transferred to the case of discrete-time filters whose transfer functions are expressed in the z domain, through the use of certain mathematical techniques such as thebilinear transform,impulse invariance, orpole–zero matching method. Thus digital IIR filters can be based on well-known solutions for analog filters such as theChebyshev filter,Butterworth filter, andelliptic filter, inheriting the characteristics of those solutions.

Transfer function derivation

[edit]

Digital filters are often described and implemented in terms of thedifference equation that defines how the output signal is related to the input signal:

y[n]=b0x[n]+b1x[n1]++bPx[nP]+a1y[n1]+a2y[n2]++aQy[nQ]{\displaystyle {\begin{aligned}y[n]{}=&b_{0}x[n]+b_{1}x[n-1]+\cdots +b_{P}x[n-P]\\&{}+a_{1}y[n-1]+a_{2}y[n-2]+\cdots +a_{Q}y[n-Q]\end{aligned}}}

where:

A more condensed form of the difference equation is:

 y[n]=i=0Pbix[ni]+i=1Qaiy[ni]{\displaystyle \ y[n]=\sum _{i=0}^{P}b_{i}x[n-i]+\sum _{i=1}^{Q}a_{i}y[n-i]}

To find thetransfer function of the filter, we first take theZ-transform of each side of the above equation to obtain:

 Y(z)=X(z)i=0Pbizi+Y(z)i=1Qaizi{\displaystyle \ Y(z)=X(z)\sum _{i=0}^{P}b_{i}z^{-i}+Y(z)\sum _{i=1}^{Q}a_{i}z^{-i}}

After rearranging:

 Y(z)[1i=1Qaizi]=X(z)i=0Pbizi{\displaystyle \ Y(z)\left[1-\sum _{i=1}^{Q}a_{i}z^{-i}\right]=X(z)\sum _{i=0}^{P}b_{i}z^{-i}}

We then define the transfer function to be:

H(z)=Y(z)X(z)=i=0Pbizi1i=1Qaizi{\displaystyle H(z)={\frac {Y(z)}{X(z)}}={\frac {\sum _{i=0}^{P}b_{i}z^{-i}}{1-\sum _{i=1}^{Q}a_{i}z^{-i}}}}
Simple IIR filter block diagram
An example of a block diagram of an IIR filter. Thez1{\displaystyle z^{-1}} block is a unit delay.

Stability

[edit]

The transfer function allows one to judge whether or not a system isbounded-input, bounded-output (BIBO) stable. To be specific, the BIBO stability criterion requires that theROC of the system includes the unit circle. For example, for acausal system to be stable, allpoles of the transfer function have to have an absolute value smaller than one. In other words, all poles must be located within a unit circle in thez{\displaystyle z}-plane.

The poles are defined as the values ofz{\displaystyle z} which make the denominator ofH(z){\displaystyle H(z)} equal to 0:

 0=j=0Qajzj{\displaystyle \ 0=\sum _{j=0}^{Q}a_{j}z^{-j}}

Clearly, ifaj0{\displaystyle a_{j}\neq 0} then the poles are not located at the origin of thez{\displaystyle z}-plane. This is in contrast to theFIR filter where all poles are located at the origin, and is therefore always stable.

IIR filters are sometimes preferred over FIR filters because an IIR filter can achieve a much sharper transition regionroll-off than an FIR filter of the same order.

Example

[edit]

Let the transfer functionH(z){\displaystyle H(z)} of adiscrete-time filter be given by:

H(z)=B(z)A(z)=11az1{\displaystyle H(z)={\frac {B(z)}{A(z)}}={\frac {1}{1-az^{-1}}}}

governed by the parametera{\displaystyle a}, a real number with0<|a|<1{\displaystyle 0<|a|<1}.H(z){\displaystyle H(z)} is stable and causal with a pole ata{\displaystyle a}. The time-domainimpulse response can be shown to be given by:

h(n)=anu(n){\displaystyle h(n)=a^{n}u(n)}

whereu(n){\displaystyle u(n)} is theunit step function. It can be seen thath(n){\displaystyle h(n)} is non-zero for alln0{\displaystyle n\geq 0}, thus an impulse response which continues infinitely.

IIR filter example

Advantages and disadvantages

[edit]

The main advantage digital IIR filters have over FIR filters is their efficiency in implementation, in order to meet a specification in terms of passband, stopband, ripple, and/or roll-off. Such a set of specifications can be accomplished with a lower order (Q in the above formulae) IIR filter than would be required for an FIR filter meeting the same requirements. If implemented in a signal processor, this implies a correspondingly fewer number of calculations per time step; the computational savings is often of a rather large factor.

On the other hand, FIR filters can be easier to design, for instance, to match a particular frequency response requirement. This is particularly true when the requirement is not one of the usual cases (high-pass, low-pass, notch, etc.) which have been studied and optimized for analog filters. Also FIR filters can be easily made to belinear phase (constantgroup delay vs frequency)—a property that is not easily met using IIR filters and then only as an approximation (for instance with theBessel filter). Another issue regarding digital IIR filters is the potential forlimit cycle behavior when idle, due to the feedback system in conjunction with quantization.

Design Methods

[edit]

Impulse Invariance

[edit]

Impulse invariance is a technique for designing discrete-time infinite-impulse-response (IIR) filters from continuous-time filters in which the impulse response of the continuous-time system is sampled to produce the impulse response of the discrete-time system. Impulse invariance is one of the commonly used methods to meet the two basic requirements of the mapping from the s-plane to the z-plane. This is obtained by solving the T(z) that has the same output value at the same sampling time as the analog filter, and it is only applicable when the inputs are in a pulse.
Note that all inputs of the digital filter generated by this method are approximate values, except for pulse inputs that are very accurate. This is the simplest IIR filter design method. It is the most accurate at low frequencies, so it is usually used in low-pass filters.

For Laplace transform or z-transform, the output after the transformation is just the input multiplied by the corresponding transformation function, T(s) or T(z). Y(s) and Y(z) are the converted output of input X(s) and input X(z), respectively.

Y(s)=T(s)X(s){\displaystyle Y(s)=T(s)X(s)}
Y(z)=T(z)X(z){\displaystyle Y(z)=T(z)X(z)}

When applying the Laplace transform or z-transform on the unit impulse, the result is 1. Hence, the output results after the conversion are

Y(s)=T(s){\displaystyle Y(s)=T(s)}
Y(z)=T(z){\displaystyle Y(z)=T(z)}

Now the output of the analog filter is just the inverse Laplace transform in the time domain.

y(t)=L1[Y(s)]=L1[T(s)]{\displaystyle y(t)=L^{-1}[Y(s)]=L^{-1}[T(s)]}

If we use nT instead of t, we can get the output y(nT) derived from the pulse at the sampling time. It can also be expressed as y(n)

y(n)=y(nT)=y(t)|t=sT{\displaystyle y(n)=y(nT)=y(t)|_{t=sT}}

This discrete time signal can be applied z-transform to get T(z)

T(z)=Y(z)=Z[y(n)]{\displaystyle T(z)=Y(z)=Z[y(n)]}
T(z)=Z[y(n)]=Z[y(nT)]{\displaystyle T(z)=Z[y(n)]=Z[y(nT)]}
T(z)=Z{L1[T(s)]t=nT}{\displaystyle T(z)=Z\left\{L^{-1}[T(s)]_{t=nT}\right\}}

The last equation mathematically describes that a digital IIR filter is to perform z-transform on the analog signal that has been sampled and converted to T(s) by Laplace, which is usually simplified to

T(z)=Z[T(s)]T{\displaystyle T(z)=Z[T(s)]*T}

Pay attention to the fact that there is a multiplier T appearing in the formula. This is because even if the Laplace transform and z-transform for the unit pulse are 1, the pulse itself is not necessarily the same. For analog signals, the pulse has an infinite value but the area is 1 at t=0, but it is 1 at the discrete-time pulse t=0, so the existence of a multiplier T is required.

Step Invariance

[edit]

Step invariance is a better design method than impulse invariant. The digital filter has several segments of input with different constants when sampling, which is composed of discrete steps. The step invariant IIR filter is less accurate than the same input step signal to the ADC. However, it is a better approximation for any input than the impulse invariant.
Step invariant solves the problem of the same sample values when T(z) and T(s) are both step inputs. The input to the digital filter is u(n), and the input to the analog filter is u(t). Apply z-transform and Laplace transform on these two inputs to obtain the converted output signal.
Perform z-transform on step inputZ[u(n)]=zz1{\displaystyle Z[u(n)]={\dfrac {z}{z-1}}}
Converted output after z-transformY(z)=T(z)U(z)=T(z)zz1{\displaystyle Y(z)=T(z)U(z)=T(z){\dfrac {z}{z-1}}}
Perform Laplace transform on step inputL[u(t)]=1s{\displaystyle L[u(t)]={\dfrac {1}{s}}}
Converted output after Laplace transformY(s)=T(s)U(s)=T(s)s{\displaystyle Y(s)=T(s)U(s)={\dfrac {T(s)}{s}}}
The output of the analog filter is y(t), which is the inverse Laplace transform of Y(s). If sampled every T seconds, it is y(n), which is the inverse conversion of Y(z).These signals are used to solve for the digital filter and the analog filter and have the same output at the sampling time.
The following equation points out the solution of T(z), which is the approximate formula for the analog filter.

T(z)=z1zY(z){\displaystyle T(z)={\dfrac {z-1}{z}}Y(z)}
T(z)=z1zZ[y(n)]{\displaystyle T(z)={\dfrac {z-1}{z}}Z[y(n)]}
T(z)=z1zZ[Y(s)]{\displaystyle T(z)={\dfrac {z-1}{z}}Z[Y(s)]}
T(z)=z1zZ[T(s)s]{\displaystyle T(z)={\dfrac {z-1}{z}}Z[{\dfrac {T(s)}{s}}]}

Bilinear Transform

[edit]

The bilinear transform is a special case of a conformal mapping, often used to convert a transfer functionHa(s){\displaystyle H_{a}(s)} of a linear, time-invariant (LTI) filter in the continuous-time domain (often called an analog filter) to a transfer functionHd(z){\displaystyle H_{d}(z)} of a linear, shift-invariant filter in the discrete-time domain.The bilinear transform is a first-order approximation of the natural logarithm function that is an exact mapping of thez-plane to thes-plane. When the Laplace transform is performed on a discrete-time signal (with each element of the discrete-time sequence attached to a correspondingly delayed unit impulse), the result is precisely the Z transform of the discrete-time sequence with the substitution of

z=esT=esT/2esT/21+sT/21sT/2{\displaystyle {\begin{aligned}z&=e^{sT}\\&={\frac {e^{sT/2}}{e^{-sT/2}}}\\&\approx {\frac {1+sT/2}{1-sT/2}}\end{aligned}}}

whereT{\displaystyle T} is the numerical integration step size of the trapezoidal rule used in the bilinear transform derivation; or, in other words, the sampling period. The above bilinear approximation can be solved fors{\displaystyle s} or a similar approximation fors=(1/T)ln(z){\displaystyle s=(1/T)\ln(z)} can be performed.

The inverse of this mapping (and its first-order bilinear approximation) is

s=1Tln(z)=2T[z1z+1+13(z1z+1)3+15(z1z+1)5+17(z1z+1)7+]2Tz1z+1=2T1z11+z1{\displaystyle {\begin{aligned}s&={\frac {1}{T}}\ln(z)\\&={\frac {2}{T}}\left[{\frac {z-1}{z+1}}+{\frac {1}{3}}\left({\frac {z-1}{z+1}}\right)^{3}+{\frac {1}{5}}\left({\frac {z-1}{z+1}}\right)^{5}+{\frac {1}{7}}\left({\frac {z-1}{z+1}}\right)^{7}+\cdots \right]\\&\approx {\frac {2}{T}}{\frac {z-1}{z+1}}\\&={\frac {2}{T}}{\frac {1-z^{-1}}{1+z^{-1}}}\end{aligned}}}

This relationship is used in the Laplace transfer function of any analog filter or the digital infinite impulse response (IIR) filter T(z) of the analog filter.
The bilinear transform essentially uses this first order approximation and substitutes into the continuous-time transfer function,Ha(s){\displaystyle H_{a}(s)}

s2Tz1z+1.{\displaystyle s\leftarrow {\frac {2}{T}}{\frac {z-1}{z+1}}.}

That is

Hd(z)=Ha(s)|s=2Tz1z+1=Ha(2Tz1z+1). {\displaystyle H_{d}(z)=H_{a}(s){\bigg |}_{s={\frac {2}{T}}{\frac {z-1}{z+1}}}=H_{a}\left({\frac {2}{T}}{\frac {z-1}{z+1}}\right).\ }

which is used to calculate the IIR digital filter, starting from the Laplace transfer function of the analog filter.

See also

[edit]

References

[edit]
  1. ^Oppenheim, Alan; Schafer, Ronald (2009-08-18).Discrete-Time Signal Processing. Prentice Hall Signal Processing (3rd ed.). Pearson.ISBN 978-0131988422.

External links

[edit]
Retrieved from "https://en.wikipedia.org/w/index.php?title=Infinite_impulse_response&oldid=1338682813"
Categories:
Hidden categories:

[8]ページ先頭

©2009-2026 Movatter.jp