Sign in

username:

password:



Not a member?

Search Online Books



Search tips

Free Online Books

Sponsor

NEW! TMS320C6474: 3x the performance. 1/3 the cost. Three 1 GHz cores on 1 chip.

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?

  

Multi-Input, Multi-Output (MIMO) Allpass Filters

To generalize lossless filters to the multi-input, multi-output (MIMO) case, we must generalize conjugation to MIMO transfer function matrices:



Theorem: A $ p\times q$ transfer function matrix $ \mathbf{H}(z)$ is lossless if and only if its frequency-response matrix $ \mathbf{H}(e^{j\omega})$ is unitary, i.e.,

$\displaystyle \mathbf{H}^*(e^{j\omega})\mathbf{H}(e^{j\omega}) = \mathbf{I}_q \protect$ (C.2)

for all $ \omega$, where $ \mathbf{I}_q$ denotes the $ q\times q$ identity matrix, and $ \mathbf{H}^\ast(e^{j\omega})$ denotes the Hermitian transpose (complex-conjugate transpose) of $ \mathbf{H}(e^{j\omega})$:

$\displaystyle \mathbf{H}^*(e^{j\omega}) \isdef \overline{\mathbf{H}^T(e^{j\omega})}
$

Let $ \underline{y}_p(n)$ denote the length $ p$ output vector at time $ n$, and let $ {\underline{x}}_q(n)$ denote the input $ q$-vector at time $ n$. Then in the frequency domain we have $ \underline{Y}_p(e^{j\omega})=\mathbf{H}(e^{j\omega})\underline{X}_q(e^{j\omega})$, which implies

$\displaystyle \underline{Y}_p^*\underline{Y}_p = \underline{X}_q^*\underbrace{\...
...H}(e^{j\omega})}_{\mathbf{I}_q}\underline{X}_q
= \underline{X}^*\underline{X},
$

or

$\displaystyle \sum_{i=1}^p\left\vert Y_i(e^{j\omega})\right\vert^2 = \sum_{i=1}^q\left\vert X_i(e^{j\omega})\right\vert^2.
$

Integrating both sides of this equation with respect to $ \omega$ yields that the total energy in equals the total energy out, as required by the definition of losslessness.

We have thus shown that in the MIMO case, losslessness is equivalent to having a unitary frequency-response matrix. A MIMO allpass filter is therefore any filter with a unitary frequency-response matrix.

Note that $ \mathbf{H}^*(e^{j\omega})\mathbf{H}(e^{j\omega})$ is a $ q\times q$ matrix product of a $ q\times p$ times a $ p\times q$ matrix. If $ q>p$, then the rank must be deficient. Therefore, $ p\geq q$. (There must be at least as many outputs as there are inputs, but it's ok to have extra outputs.)



Subsections

Order a Hardcopy of Introduction to Digital Filters

Previous: Paraunitary FiltersC.4
Next: Paraunitary MIMO Filters

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