Sign in

username:

password:



Not a member?

Search Online Books



Search tips

Free Online Books

Chapters

Chapter Contents:

Search Introduction to Digital Filters

  

Book Index | Global Index


Would you like to be notified by email when Julius Orion Smith III publishes a new entry into his blog?

  

FIR Part of a PFE

When $ M\geq N$ in Eq.$ \,$(6.7), we may perform a step of long division of $ B(z)/A(z)$ to produce an FIR part in parallel with a strictly proper IIR part:

$\displaystyle H(z) \isdefs \frac{B(z)}{A(z)} \eqsp F(z) + \sum_{i=1}^{N} \frac{r_i}{1-p_iz^{-1}} \protect$ (7.10)

where

\begin{eqnarray*}
B(z) &=& b_0 + b_1 z^{-1}+ b_2z^{-2}+ \cdots + b_M z^{-M}\\
A...
...=& f_0 + f_1z^{-1}+ f_2z^{-2}+ \cdots + f_K z^{-K}, \quad K=M-N.
\end{eqnarray*}

When $ M<N$, we define $ F(z)=0$. This type of decomposition is computed by the residuez function (a matlab function for computing a complete partial fraction expansion, as illustrated in §6.8.8 below).

An alternate FIR part is obtained by performing long division on the reversed polynomial coefficients to obtain

$\displaystyle H(z) \eqsp F(z) + z^{-(K+1)}\sum_{i=1}^{N} \frac{r_i}{1-p_iz^{-1}}, \protect$ (7.11)

where $ K=M-N\geq 0$ is again the order of the FIR part. This type of decomposition is computed (as part of the PFE) by residued, described in §J.6 and illustrated numerically in §6.8.8 below.

We may compare these two PFE alternatives as follows: Let $ A_N$ denote $ A(z)$, $ F_K\isdeftext F(z)$, and $ B_M\isdeftext B(z)$. (I.e., we use a subscript to indicate polynomial order, and `$ (z)$' is omitted for notational simplicity.) Then for $ K=M-N\geq 0$ we have two cases:

\begin{eqnarray*}
(1) && H(z) \eqsp F_K + \frac{B^\prime_{N-1}}{A_N} \eqsp \frac...
..._N} \eqsp \frac{F_K A_N + z^{-(K+1)}B^{\prime\prime}_{N-1}}{A_N}
\end{eqnarray*}

In the first form, the $ B^\prime_{N-1}$ coefficients are ``left justified'' in the reconstructed numerator, while in the second form they are ``right justified''. The second form is generally more efficient for modeling purposes, since the numerator of the IIR part ( $ B^{\prime\prime}_{N-1}(z)$) can be used to match additional terms in the impulse response after the FIR part $ F_K(z)$ has ``died out''.

In summary, an arbitrary digital filter transfer function $ H(z)$ with $ N$ distinct poles can always be expressed as a parallel combination of complex one-pole filters, together with a parallel FIR part when $ M\geq N$. When there is an FIR part, the strictly proper IIR part may be delayed such that its impulse response begins where that of the FIR part leaves off.

In artificial reverberation applications, the FIR part may correspond to the early reflections, while the IIR part provides the late reverb, which is typically dense, smooth, and exponentially decaying [85]. The predelay (``pre-delay'') control in some commercial reverberators is the amount of pure delay at the beginning of the reverberator's impulse response. Thus, neglecting the early reflections, the order of the FIR part can be viewed as the amount of predelay for the IIR part.



Subsections

Order a Hardcopy of Introduction to Digital Filters

Previous: Inverting the Z Transform
Next: Example: The General Biquad PFE

written by Julius Orion Smith III
Julius Smith's background is in electrical engineering (BS Rice 1975, PhD Stanford 1983). He is presently Professor of Music and Associate Professor (by courtesy) of Electrical Engineering at Stanford's Center for Computer Research in Music and Acoustics (CCRMA), teaching courses and pursuing research related to signal processing applied to music and audio systems. See http://ccrma.stanford.edu/~jos/ for details.


Comments


No comments yet for this page


Add a Comment
You need to login before you can post a comment (best way to prevent spam). ( Not a member? )