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?

  

Jordan Canonical Form

The block diagonal system having the eigenvalues along the diagonal and ones in some of the superdiagonal elements (which serve to couple repeated eigenvalues) is called Jordan canonical form. Each block size corresponds to the multiplicity of the repeated pole. As an example, a pole $ p_i$ of multiplicity $ m_i=3$ could give rise to the following $ 3\times 3$ Jordan block:

$\displaystyle D_i = \left[\begin{array}{ccc}
p_i & 1 & 0\\ [2pt]
0 & p_i & 1\\ [2pt]
0 & 0 & p_i
\end{array}\right]
$

The ones along the superdiagonal serve to couple the states corresponding to $ p_i$ and generate polynomial amplitude envelopes multiplying the sampled exponential $ p_i^n$.G.11Note, however, that a pole of multiplicity three can also yield two Jordan blocks, such as

$\displaystyle D_i = \left[\begin{array}{ccc}
p_i & 1 & 0\\ [2pt]
0 & p_i & 0\\ [2pt]
0 & 0 & p_i
\end{array}\right]
$

or even three Jordan blocks of order 1. The number of Jordan blocks associated with a single pole $ p_i$ is equal to the number of linearly independent eigenvectors of the transition matrix associated with eigenvalue $ p_i$. If all eigenvectors of $ A$ are linearly independent, the system can be diagonalized after all, and any repeated roots are ``uncoupled'' and behave like non-repeated roots (no polynomial amplitude envelopes).

Interestingly, neither Matlab nor Octave seem to have a numerical function for computing the Jordan canonical form of a matrix. Matlab will try to do it symbolically when the matrix entries are given as exact rational numbers (ratios of integers) by the jordan function, which requires the Maple symbolic mathematics toolbox. Numerically, it is generally difficult to distinguish between poles that are repeated exactly, and poles that are merely close together. The residuez function sets a numerical threshold below which poles are treated as repeated.


Order a Hardcopy of Introduction to Digital Filters

Previous: Repeated Poles
Next: State-Space Analysis Example: The Digital Waveguide Oscillator

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