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
of multiplicity
could give
rise to the following
Jordan block:
The ones along the superdiagonal serve to couple the states
corresponding to

and generate polynomial amplitude
envelopes
multiplying the sampled
exponential 
.
G.11Note, however, that a
pole of multiplicity three can also yield two Jordan blocks, such as
or even three Jordan blocks of order 1. The number of Jordan blocks
associated with a single pole

is equal to the number of
linearly
independent eigenvectors of the transition
matrix associated with
eigenvalue

. If all eigenvectors of

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.
Previous:
Repeated PolesNext:
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.