步骤 1:求特征值与约旦标准型 J
计算特征方程:
∣λI−A∣=λ0−2−1λ50−1λ−4=λ3−4λ2+5λ−2=(λ−1)2(λ−2)=0
解得特征值为:λ1=λ2=1(二重根),λ3=2(单根)。
由于重根 λ=1 的几何重数为 3−rank(I−A)=3−2=1<2,系统无法对角化,对应一个二阶约旦块。
约旦标准型及其指数矩阵为:
J=100110002⟹eJt=et00tetet000e2t
步骤 2:构造变换矩阵 T 与 T−1
设 T=[v1v2v3],满足 AT=TJ:
- Av1=1⋅v1⟹(I−A)v1=0⟹v1=111(主特征向量)
- Av2=v1+1⋅v2⟹(A−I)v2=v1⟹v2=−101(广义特征向量)
- Av3=2v3⟹(2I−A)v3=0⟹v3=124(特征值 2 对应的特征向量)
求得相应的变换矩阵及其逆矩阵:
T=111−101124,T−1=−2−2153−2−2−11
步骤 3:矩阵连乘求 eAt
先计算 TeJt:
TeJt=111−101124et00tetet000e2t=etetet(t−1)ettet(t+1)ete2t2e2t4e2t
再乘以 T−1:
eAt=Φ(t)=TeJtT−1=etetet(t−1)ettet(t+1)ete2t2e2t4e2t−2−2153−2−2−11=−2tet+e2t2(e2t−tet−et)−2tet−4et+4e2t3tet+2et−2e2t3tet+5et−4e2t3tet+8et−8e2t−tet−et+e2t−tet−2et+2e2t−tet−3et+4e2t
Discussion
Comments
Thoughts, corrections, and follow-up notes are welcome here.