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?

  

Example State Space Filter Transfer Function

In this example, we consider a second-order filter ($ N = 2$) with two inputs ($ p=2$) and two outputs ($ q=2$):

\begin{eqnarray*}
A &=& g\left[\begin{array}{rr} c & -s \\ [2pt] s & c \end{arra...
... \left[\begin{array}{cc} 0 & 0 \\ [2pt] 0 & 0 \end{array}\right]
\end{eqnarray*}

so that

\begin{eqnarray*}
\left[\begin{array}{c} x_1(n+1) \\ [2pt] x_2(n+1) \end{array}\...
...left[\begin{array}{c} x_1(n) \\ [2pt] x_2(n) \end{array}\right].
\end{eqnarray*}

From Eq.$ \,$(G.5), the transfer function of this MIMO digital filter is then

\begin{eqnarray*}
H(z) &=& C(zI-A)^{-1}B = (zI-A)^{-1} = \left[\begin{array}{cc}...
...z^{-2}}{\displaystyle 1-2gcz^{-1}+g^2z^{-2}} \end{array}\right].
\end{eqnarray*}

Note that when $ g=1$, the state transition matrix $ A$ is simply a 2D rotation matrix, rotating through the angle $ \theta$ for which $ c=\cos(\theta)$ and $ s=\sin(\theta)$. For $ g<1$, we have a type of normalized second-order resonator [51], and $ g$ controls the ``damping'' of the resonator, while $ \theta =
2\pi f_r/f_s$ controls the resonance frequency $ f_r$. The resonator is ``normalized'' in the sense that the filter's state has a constant $ L2$ norm (``preserves energy'') when $ g=1$ and the input is zero:

$\displaystyle \left\Vert\,{\underline{x}}(n+1)\,\right\Vert \isdef \sqrt{x_1^2(...
...x}}(n)\,\right\Vert \equiv \left\Vert\,{\underline{x}}(n)\,\right\Vert \protect$ (G.6)

since a rotation does not change the $ L2$ norm, as can be readily checked.

In this two-input, two-output digital filter, the input $ u_1(n)$ drives state $ x_1(n)$ while input $ u_2(n)$ drives state $ x_2(n)$. Similarly, output $ y_1(n)$ is $ x_1(n)$, while $ y_2(n)$ is $ x_2(n)$. The two-by-two transfer-function matrix $ H(z)$ contains entries for each combination of input and output. Note that all component transfer functions have the same poles. This is a general property of physical linear systems driven and observed at arbitrary points: the resonant modes (poles) are always the same, but the zeros vary as the input or output location are changed. If a pole is not visible using a particular input/output pair, we say that the pole has been ``canceled'' by a zero associated with that input/output pair. In control-theory terms, the pole is ``uncontrollable'' from that input, or ``unobservable'' from that output, or both.


Order a Hardcopy of Introduction to Digital Filters

Previous: Transfer Function of a State Space Filter
Next: Transposition of a State Space Filter

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? )