From octave-sources-request at bevo dot che dot wisc dot edu Thu Dec  5 13:15:28 2002
Subject: Re: Jordan Form
From: A S Hodel <a dot s dot hodel at eng dot auburn dot edu>
To: bebark2 at uky dot edu
Cc: octave-sources at bevo dot che dot wisc dot edu
Date: Thu, 5 Dec 2002 13:15:50 -0600

What you're proposing is a numerically risky procedure; see Golub and 
Van Loan's book Matrix Computations.

Otherwise,  the procedure is (roughly):

[V,Lambda] = eig(A);
Bh = V\B;
Ch =  C*V;
Ah = Lambda;

JordanSys = ss2sys(Ah, Bh, Ch,D);

The dangerous step is V\B.  If A has Jordan blocks, you've got a 
problem.

Take a look at the octave controls toolbox functions is_controllable, 
is_stabilizable, etc. for other control systems subspace issues and 
techniques.

On Thursday, December 5, 2002, at 12:40 PM, Brandon Elam Barker wrote:

> Has anyone written a program to put square matrices into Jordan Form?  
> If not, I may do this as a way to learn Octave.
>
> Thanks,
> Brandon Barker
>
>
>
A. S. Hodel, Assoc. Prof, Dept. Elect & Comp Eng, Auburn University, AL 
36849-5201
(334) 844-1854 200 Broun Hall hodelas at auburn dot edu 
http://ww.eng.auburn.edu/~scotte