Zen Pharaohs of Ultimate Love mathematics, music, science, wood....what else could you ask?    Next: Bidiagonal Representation of TIB Up: Triangular Input Balanced Realizations Previous: Meixner Functions and TIB

## Gauss-Meixner Summation

The impulse response of a TIB realization with a single eigenvalue is given by where and, as usual, c is the vector of coefficients. In order to compute c we exploit Gaussian quadrature (in this case summation) properties of the orthogonal polynomials. The normalized three-term recurrence for the Meixner polynomials results in the Jacobi matrix and the associated Gauss-Meixner summation formula where p is any polynomial of degree . The are the eigenvalues of J and the are proportional to the squares of the first components of the corresponding orthonormal eigenvectors. We fix the by choosing the proportionality constant so that for consistency, (i.e. consider p=1).

We put this in terms of a summation formula for the Meixner functions: This allows us to compute Meixner coefficients, that is for of the form by the formula Now we can compute recursively without explicitly computing as follows. then using the recursion for : and where we put . A check that this has worked is that for k=1,...,d.

Gauss-Meixner summation works very well when we can bound the order d of the model in advance, and it is one attractive way of computing the TIB filter realization of the maximally flat impulse response.

Fri Jun 27 03:10:38 EDT 1997
 Copyright © 1998,1999 Zen Pharaohs of Ultimate Love. All rights reserved.