next up previous


Algorithm of Filipov


If n=1, then the Jordan block coincides with the given matrix and formula (9) is true. Let us suppose that the Jordan form of the matrix A is found by applying the Jordan block construction formula (9) if the order of the matrix A is smaller than n. we will use mathematical induction.

I step. Assuming that A is singular, tex2html_wrap_inline8619 Considering the corresponding tex2html_wrap_inline8621 matrix we find that in this case the construction based on formulas (9) is realizable. Namely, in the space tex2html_wrap_inline8623 there are r independent vectors tex2html_wrap_inline8627 such that the following relations take place
 
equation3364

II step. Let us suppose that tex2html_wrap_inline8629 Every vector of the null space tex2html_wrap_inline8631 is an eigenvector of the matrix A corresponding to the eigenvalue of the matrix A tex2html_wrap_inline8637 Therefore, there must be p chains on the I step which begin with the eigenvectors corresponding to the eigenvalue 0. We are interested in the last vector of each such chain. Since the vectors tex2html_wrap_inline8627 belonging to the subspace tex2html_wrap_inline8645 must also belong to the space tex2html_wrap_inline8071 then they have to be the linear combinations of the column vectors of the matrix A
displaymath8585
with some tex2html_wrap_inline8651. Therefore, the vector tex2html_wrap_inline8651 follows the vector tex2html_wrap_inline8627 in the chain corresponding to the eigenvalue tex2html_wrap_inline8523.

III step. Since tex2html_wrap_inline8659 there must be n-r-p more linearly independent vectors tex2html_wrap_inline8663 of the space tex2html_wrap_inline8631 in the orthogonal complement of the subspace tex2html_wrap_inline8645.

Proposition 5.3.1. The algorithm of Filipov defines r vectors tex2html_wrap_inline8671 p vectors tex2html_wrap_inline8651 and n-r-p vectors tex2html_wrap_inline8679 which determine the Jordan chains. These vectors are linearly independent and suit to be the column-vectors of the matrix X, and J=X-1AX.

Proof. Look Strang (1988, p. 457). tex2html_wrap_inline7853

Example 5.3.1. Let us find the Jordan normal form of the matrix
displaymath8586
using the algorithm of Filipov.

I step. It comes out from the form of the matrix that tex2html_wrap_inline8687 and tex2html_wrap_inline8689 Hence r=1 and there is a vector tex2html_wrap_inline8693 from this subspace tex2html_wrap_inline8623 satisfying the condition (10).

II step. Let us find the basis of the null space tex2html_wrap_inline8631 of the matrix A:
displaymath8587
The vector tex2html_wrap_inline8701 belongs to the subspace tex2html_wrap_inline8645 and tex2html_wrap_inline8705 We solve the system

displaymath8588

III step. We take for the vector tex2html_wrap_inline8707 the vector tex2html_wrap_inline8709 and form the matrix X:
displaymath8589
Now we find the inverse matrix
displaymath8590

displaymath8591
and the Jordan matrix
displaymath8592
The software package ``Maple'' gives for the Jordan decomposition:

displaymath8593

Since the matrix X in the Jordan decomposition of the matrix A is not uniquely defined, then for many problems it is of interest to choose the matrix X so that the conditional number k(X) were the least. Such a problem arose also in example 1.2.9.4.

Problem 5.3.1. Find the Jordan decomposition of the matrix
displaymath8594

Problem 5.3.2.* Find the Jordan decomposition of the matrix
displaymath8595

Problem 5.3.3.* Find the Jordan decomposition of the matrix
displaymath8596

Problem 5.3.4. Let the Jordan decomposition of the matrix tex2html_wrap_inline7067 be A=MJM-1. Show that tex2html_wrap_inline8729


next up previous