Lecture VII
• The Jordan Canonical Form • Examples and Applications
Fall 2020 Prof.Jiang@ECE NYU 295
Review of Canonical Forms
If A is an n n matrix with distinct eigenvalues, then there exists a nonsingular matrix P s.t.
P1 AP diag i .
If A is Hermitian (possibly having repeated eigenvalues),
a unitary matrix U such that U*AU diag .
i
Fall 2020 Prof.Jiang@ECE NYU 296
A Motivating Example
As stated previously, not every matrix can be transformed into a canonical diagonal form. For example,
A1 b, b0 0 1
cannot be transformed into a diagonal matrix. Fall 2020 Prof.Jiang@ECE NYU 297
Jordan Canonical Form
Theorem (Jordan):
Let A be an n n matrix whose different eigenvalues
are , , with multiplicities m ,, m : 1s 1s
i.e., nonsingular P such that
s
i1
m
det I A
Then, A is transformable into a Jordan canonical form.
P1AP blockdiagi J
Fall 2020 Prof.Jiang@ECE NYU 298
i i
Fall 2020
i
Prof.Jiang@ECE NYU 299
Theorem (Jordan), cont’d: P1AP blockdiagi J
where
i 0 0
1 0 i
i 00 1 0
Jordan Block
i 0 0 1
Comments
In some texts, J T is used as Jordan form. Different Jordan blocks, say i , j may be
associated with the same eigenvalues.
s ThetotalnumberofJordanblocks:ss n.
Fall 2020 Prof.Jiang@ECE NYU 300
Fall 2020
Prof.Jiang@ECE NYU 301
Illustration via 3×3 matrices If a 33 matrix A has an eigenvalue 1
of multiplicity three, then it may be reduced into one of the following Jordan forms:
1 0 0 1 0 0 J0 0,J0 0,
1121 00 01
11 1 0 0
J1 0. 31
01 1
Remark 1
The distinct Jordan forms Ji , Jk , i k, are not similar to each other.
Fall 2020 Prof.Jiang@ECE NYU 302
Remark 2
When each Jordan block ( ) in the Jordan ii
form J is one-dimensional (i.e. n 1) and s n, i
the Jordan matrix J becomes diagonal.
Fall 2020 Prof.Jiang@ECE NYU 303
Application to Matrix Analysis of Differential Equations
Given a set of 1st-order differential equations
x(t)Ax(t), x(0)n,
applying the transformation y P1x yields:
y(t)P1APy(t) : Jy(t).
y1
yi(t)yi(t),yimi,y .
Fall 2020 Prof.Jiang@ECE NYU
304
i
ys
Comment
So, with the help of Jordan canonical form, solving differential equations can be reduced down to solving lower-order (disjoint!) differential equations.
(see a forthcoming lecture.)
Fall 2020 Prof.Jiang@ECE NYU 305
Principal Vectors
In order to develop a constructive method for P resulting in Jordan form, let’s introduce the notion of principal vector, or generalized eigenvector,
which is a generalization of eigenvector.
Fall 2020 Prof.Jiang@ECE NYU 306
Principal Vectors
A (possibly zero) vector p is a principal vector of grade g 0 belonging to the eigenvalue i if
iIAg p0,
for which g is the smallest non-negative integer.
Fall 2020 Prof.Jiang@ECE NYU 307
Examples
• The vector p = 0 is the principal vector of grade 0.
• The (nonzero) eigenvectors are the principal vectors of grade 1.
Fall 2020 Prof.Jiang@ECE NYU 308
Motivating Question
In case of transformation to diagonal canonical form, i.e., P1 AP diag i ,
the columns of P are linearly independent eigenvectors.
What about the matrix P in Jordan form?
How to construct P from principal vectors? Fall 2020 Prof.Jiang@ECE NYU 309
gi
i
Linear Spaces
Define the linear space composed of all principal vectors of grade g belonging to i :
Pp|IAp0 g
i.e., the null space of i I Ag . Clearly,
P P P 0i1i2i
Fall 2020 Prof.Jiang@ECE NYU 310
An Interesting Result
Let A be an n n matrix with the distinct
eigenvalues 1, , s , 1 s n, with
multiplicities m , , m . 1s
Then, every vector x n can be written as x p1 p2 ps
where pi is a uniquely defined principal vector associated with i of grade mi .
Fall 2020 Prof.Jiang@ECE NYU 311
Comment 1
A special, but interesting, case is when there are n linearly independent eigenvectors, say,
c1, , cn. In this case, i scalars s.t. xc1cn :p1pn.
1n
Fall 2020 Prof.Jiang@ECE NYU 312
Comment 2
Its proof relies upon the well-known Cayley- Hamilton theorem; see any standard matrix or linear algebra textbook.
Fall 2020 Prof.Jiang@ECE NYU 313
1 0 0 Consider the matrix A 7 1 0.
Example
0 0 2
Compute its eigenvalues and the associated eigenvectors. * Can each column be written as a linear combination
of eigenvectors?
* Show that each column can be written
as a unique representation of principal vectors.
Fall 2020 Prof.Jiang@ECE NYU 314
Fall 2020
Prof.Jiang@ECE NYU 315
Answer
1 m 2, 2 m 1. 1122
The eigenvectors of 1 are of the form col0,1,0, any nonzero scalar.
The eigenvectors of 2 are of the form col 0, 0,1, any nonzero scalar.
P p2|p2col0,0,.
P p1|p1col,,0
21 12
Cayley-Hamilton Theorem Revisited
For any n n matrix A,
(A)An 1An1n1AnIO
where is the characteristic polynomial of A, i.e.,
detIAn n ini. i1
Fall 2020 Prof.Jiang@ECE NYU
316
Fall 2020
Prof.Jiang@ECE NYU
317
() 2 23. A
Example
ConsiderA7 4.Verifythat 8 5
1) The characteristic polynomial A () is:
2)(A)A22A3I00 0.
A
0 0
Another Proof
Define the n n matrix of signed cofactors: C cof I A.
T
Then, usingMcofM (detM)I,
IACT I. In addition,
CT n1C0 Cn1 Cn
for constant matrices Ci ‘ s.
Fall 2020 Prof.Jiang@ECE NYU 318
C0 I
C AC I
101
Proof (cont’d)
By identification of the coefficients of equal powers of gives
Cn1 ACn2 n1I ACn1 nI.
Multiplying the first eq. by An , the second by An1,…,
and then adding them up leads to: O AI.
Fall 2020 Prof.Jiang@ECE NYU 319
Question
How to compute principal vectors for a given matrix?
Fall 2020 Prof.Jiang@ECE NYU 320
A Motivating Example
Consider a 22 Jordan block J 0. 1
Denote P x1 x2 that transforms A into J.
Namely, P1AP J. So, we have AP PJ, or
Ax1 x2x1 x2 0 1
Ax2 x2 , so x2 is an eigenvector;
(AI)x1 x2, so x1 is a principal vector (of grade 2). Fall 2020 Prof.Jiang@ECE NYU 321
Comment
Usually, x , x is called a Jordan Basis for 1 2
this 22 matrix A. In order words,
the JCF transformation matrix P is composed
of a Jordan basis, or a set of linearly independent eigenvectors and principal vectors.
Fall 2020 Prof.Jiang@ECE NYU 322
General Procedure
Step 1: Solve the characteristic equation AIz1 0.
Step 2 : For each independent z1 , solve AIz2 z1
2 wherez2 clearlysolves AI z2 0.
Collect only those z2 which are linearly independent
with the previously found eigenvectors z1.
Fall 2020 Prof.Jiang@ECE NYU 323
General Procedure
Step 3: For each independent z 2 , solve AIz3 z2
3 wherez3 clearlysolves AI z3 0.
Collect only those z3 which are linearly independent with the previously found vectors z1, z2.
Step 4: Continue in this way till the total number of independent eigenvectors and principal vectors equals
to the (algebraic) multiplicity of .
Fall 2020 Prof.Jiang@ECE NYU 324
General Procedure
Step 4 (cont’d): Denote
1 2 mm m1 1
x,x,,xz,z ,,z and
P x1, x2 ,, xm .
eigenvector
Therefore,
P1 AP J (associated with eigenvalue ).
Fall 2020 Prof.Jiang@ECE NYU 325
Comments
• Not any arbitrary choice of linearly independent principal vectors would lead to a correct transformation matrix P.
For example, at Step 2, the linearly independent principal vectors z2 are chosen according to
AIz2 z1 but NOT :
IAz2 z1.
•
Fall 2020 Prof.Jiang@ECE NYU 326
See (the 1960 book of Gantmacher, Vol.1, Chap. VI, Section 8) for another general method of constructing a transformation matrix.
Fall 2020
x4 :AIx3
Prof.Jiang@ECE NYU 327
More on Jordan Basis
Without going into the full details in proving Jordan’s Theorem, let’s illustrate the concept of Jordan basis and its use in the canonical transformation.
Consider a principal vector v of grade g n 4. Define:
x1 : v
x2 :AIx1
x3 :AIx2
Jordan Basis
That is,
eigenvector
Jordan Basis (cont’d)
Then, the 44 matrix A can be transformed into the Jordan canonical form:
0 0 0
100 J 010
Principal vector of grade 4
0 0 1
P1APJ, Px1,x2,x3,x4.
Fall 2020 Prof.Jiang@ECE NYU 328
Comment
If we define P x4,x3,x2,x1 , then A is transformed
into the Jordan canonical form J T , i.e.: 1 0 0
Fall 2020
Prof.Jiang@ECE NYU 329
010 1T
P APJ 0 0 1. 0 0 0
A More Complex Case
If AI has rank n2, i.e. its null space is of dimension 2, then two linearly independent eigenvectorsto AIq0.
Thus, we need n 2 linearly independent principal vectors. In this case, the Jordan basis takes the form
v1,v2,,vk and u1,u2,,ul , kln.
So, A is transformed into the Jordan canonical form P APdiag J1,J2 .
1
Fall 2020 Prof.Jiang@ECE NYU 330
Fall 2020
Prof.Jiang@ECE NYU 331
Exercise 1
Find a transformation matrix P to bring the following matrix
M1 b, b0 0 1
into the Jordan Canonical Form J1 0.
1 1
Fall 2020
Prof.Jiang@ECE NYU 332
Exercise 2
Find a transformation matrix to bring the following matrix into a Jordan form:
1 11 1
3 3 5 4 A 8 4 3 4
15 10 11 11
Solution:
Fall 2020
Prof.Jiang@ECE NYU 333
0 1 0 1 1 5 0 5
P , 0 4 1 5
1110 12
1 1 0 0
0 1 1 0 J .
0 0 1 0 0 0 0 1
I A becomes (after elementary operations on rows and columns: 100 0
0 1 0 0 .
0 0 1 0
3 0 0 0 1
Therefore, the matrix has two elementary divisors:
1 and 13 ,
which give two Jordan blocks, respectively:
1 1 0 J1,J0 1 1.
12 0 0 1
See (the 1960 book of Gantmacher, Vol.1, pp.160-164)
for the details.
Fall 2020 Prof.Jiang@ECE NYU 334
Practicing Problems for Midterm
1. Compute the eigenvalues of the matrix A7 2
41
and transform it to one of the canonical forms.
Fall 2020 Prof.Jiang@ECE NYU 335
Practicing Problems for Midterm
2. Consider the block diagonal matrix A0
Fall 2020
Prof.Jiang@ECE NYU 336
A 1 ,withAnn,nnn. iii12
0A 2
Show that the eigenvalues of A are those of
A and A . 12
Practicing Problems for Midterm
3. Assume A is a nonsingular matrix. If is an eigenvalue of A with eigenvector x,
showthat1 isaneigenvalueofA1.
In addition, give an eigenvector associated
with 1.
Fall 2020 Prof.Jiang@ECE NYU 337
Practicing Problems for Midterm
4. Show that A 0 1 cannot be transformed into 0 0
a diagonal matrix under any similarity transformation.
Fall 2020 Prof.Jiang@ECE NYU 338
Practicing Problems for Midterm
5. For any given 2 2 real orthogonal matrix U , one of the following must hold:
(i) U cos sin for some ; sin cos
(ii)U0 1cos sinforsome.
1 0sin cos
(Only for those who love math proof!)
Fall 2020 Prof.Jiang@ECE NYU 339
Practicing Problems for Midterm
6. Show that JT is similar to J. That is, 01 01
JT J.
10 10
Fall 2020 Prof.Jiang@ECE NYU 340
Practicing Problems for Midterm
7. Assume that A, D are invertible matrices. Show that
A B 1 A1 A1BD1 0D0 D1 .
Fall 2020 Prof.Jiang@ECE NYU 341
Practicing Problems for Midterm
8. Assume that A, D are invertible matrices. Show that
A B1 A1 A1BECA1 A1BE C D ECA1 E
where E is the inverse of the Schur complement 1
ofA: EDCA1B .
Note: A Very Useful Identity.
Fall 2020 Prof.Jiang@ECE NYU 342
Practicing Problems for Midterm
9. Reduce the following matrix into a canonical diagonal form:
Fall 2020
Prof.Jiang@ECE NYU 343
A M 022 0 M
22 where
M0 1 1 0
Practicing Problems for Midterm
10. Reduce the following matrix into a Jordan canonical form:
Fall 2020
Prof.Jiang@ECE NYU 344
3 2 1 A0 3 0
0 0 3
Practicing Problems for Midterm
11. Rank Inequalities (See Horn-Johanson text, page 13)
Sylvester inequality Amk, Bkn, wehave
rankArankBkrankABminrankA, rankB. Frobenius inequality
Amk, Bkp, Cpn, wehave
rankAB rankBC rankB + rankABC
with equality iff there are matrices X and Y such that
B BCX YAB.
Fall 2020 Prof.Jiang@ECE NYU 345
Fall 2020
Prof.Jiang@ECE NYU 346
1. For the matrix
1 0 1
Homework #7
A0 2 0,
0 0 1
identify the spaces P and the principal g
vectors of grade 2.
Fall 2020
Prof.Jiang@ECE NYU 347
Homework #7
2. Express the following vectors as unique representations of principal vectors found in Problem 1:
2 x9
0 , x9.3.
84
0
Fall 2020
Prof.Jiang@ECE NYU 348
Homework #7
3. Can you transform the following matrix into a Jordan form:
A0 , 0?
0 0