markdown
Diagonalization, Jordan Form, and Linear Systemsmd 976ab59
lecture/math/differential-equations/diagonalization-jordan-form-and-systems.lecture.n.md
Download PDF

Diagonalization, Jordan Form, and Linear Systems

date2026-07-16document_iddoc_1391a6b57f52591851f8c1b3116570c4descriptionJordan 標準形と一般化固有ベクトルの鎖を、定数係数線型系の指数関数×多項式の時間発展へ翻訳する。prerequisites一階連立系と行列指数関数 / ジョルダン標準形の入口type講義content_typelecturestatusactiverelateddata/lecture/math/differential-equations/first-order-systems-and-matrix-exponentials.lecture.n.md / data/lecture/math/linear-algebra/introduction-to-jordan-canonical-form.lecture.n.md / data/lecture/math/differential-equations/linearization-and-eigenvalue-stability.lecture.n.md / data/exercise/math/differential-equations/jordan-form-and-linear-systems.exercise.n.md
mathdifferential-equationslinear-systemslecture

1Purpose of the Lecture

For the constant-coefficient system

x=Ax,x(t0)=x0,

Jordan form translates generalized eigenvector chains into solutions consisting of an exponential multiplied by a finite-degree polynomial. We work over C. A real matrix may be complexified and placed in complex Jordan form; conjugate complex solutions then yield real solutions through their real and imaginary parts. Real Jordan form is outside the scope of this lecture.

data/lecture/math/linear-algebra/introduction-to-jordan-canonical-form.lecture.n.md

2Exponential of a Jordan Block

For a block of size k, write

Jk(λ)=λI+N,N=(0100001000010000).

The matrices λI and N commute. Hence, with τ=t-t0,

eJk(λ)τ=eλτr=0k-1τrr!Nr.

Here Nk=0 and, for k>1, Nk-10. Nilpotency terminates the series and produces polynomial factors through degree k-1. This factorization relies on commutativity and is not valid for arbitrary matrices X,Y.

3Return to the Original Coordinates

If

P-1AP=J,

then

eAτ=PeJτP-1,x(t)=PeJ(t-t0)P-1x0.

Each Jordan block may be exponentiated independently.

4Generalized Eigenvector Chains

Let

(A-λI)v1=0,(A-λI)vj=vj-1(2[PARSE ERROR: Undefined("Command(\"le\")")]j[PARSE ERROR: Undefined("Command(\"le\")")]k)

be a Jordan chain. The solution with initial value vj is

xj(t)=eλτ=0j-1τ!vj-.

In particular,

x1=eλτv1,x2=eλτ(v2+τv1),x3=eλτ(v3+τv2+τ22v1).

Differentiation and the chain relations verify the system directly. The exponential derivative contributes the term with λ, while differentiating the polynomial lowers the chain index by one; (A-λI)vj=vj-1 then gives xj=Axj.

5Example in Non-Jordan Coordinates

For

A=(21-10),v1=(1-1),v2=(10),

we have (A-I)v1=0 and (A-I)v2=v1. Moreover,

det(μI-A)=(μ-1)2,P=(11-10),J=P-1AP=(1101).

Consequently,

eAτ=eτ(1+ττ-τ1-τ).

6Asymptotic Behavior

For τ[PARSE ERROR: Undefined("Command(\"ge\")")]0, a Jordan block satisfies

eJk(λ)τ[PARSE ERROR: Undefined("Command(\"le\")")]C(1+τk-1)eReλτ.

Thus, negative real part implies decay despite polynomial factors, while positive real part produces an exponentially growing direction. When the real part is zero, blocks of size one are bounded, whereas larger blocks admit polynomially unbounded solutions. Both eigenvalue real parts and Jordan block sizes are therefore required at the spectral boundary.

7Numerical Boundary

Jordan structure is sensitive to perturbations and is primarily a theoretical classification. For example,

(λ1ελ)

is one Jordan block when ε=0, whereas its eigenvalues split into λ±ε when ε0. Numerical computation of eAt should not generally assume that an exact Jordan block structure can be identified reliably in floating-point arithmetic.

8Previous Lecture

data/lecture/math/differential-equations/first-order-systems-and-matrix-exponentials.lecture.n.md

9Next Lecture

data/lecture/math/differential-equations/linearization-and-eigenvalue-stability.lecture.n.md

10Exercises

data/exercise/math/differential-equations/jordan-form-and-linear-systems.exercise.n.md
raw .n.md をコピー
loc をコピー (filepath:line ~ line)
copy share link
copy encoded share link
path をコピー
copy share link
copy encoded share link
copy share link
copy encoded share link
タブを全て閉じる