计算实矩阵的特征值和特征向量,然后将复共轭特征值对变换为实数块。
求实矩阵的特征值和特征向量。
X = [1 1 1 1 1
0 4 5 1 1
0 -5 4 1 1
0 0 2 3 1
0 0 -3 -2 1];
[V,D] = eig(X)
V = 5×5 complex
1.0000 + 0.0000i -0.0179 - 0.1351i -0.0179 + 0.1351i 0.1593 - 0.4031i 0.1593 + 0.4031i
0.0000 + 0.0000i 0.0130 - 0.6214i 0.0130 + 0.6214i 0.0704 - 0.0267i 0.0704 + 0.0267i
0.0000 + 0.0000i 0.6363 + 0.0000i 0.6363 + 0.0000i -0.1261 + 0.1032i -0.1261 - 0.1032i
0.0000 + 0.0000i 0.1045 - 0.2087i 0.1045 + 0.2087i -0.2279 - 0.4161i -0.2279 + 0.4161i
0.0000 + 0.0000i -0.1156 + 0.3497i -0.1156 - 0.3497i 0.7449 + 0.0000i 0.7449 + 0.0000i
D = 5×5 complex
1.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 3.8801 + 5.1046i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 3.8801 - 5.1046i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 2.1199 + 0.7018i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 2.1199 - 0.7018i
D 包含一个实数特征值和两对复共轭特征值。
变换 V 和 D,使 Dnew 为实数块对角形,并且 Vnew 满足 X*Vnew = Vnew*Dnew
[Vnew,Dnew] = cdf2rdf(V,D)
Vnew = 5×5
1.0000 -0.0179 -0.1351 0.1593 -0.4031
0 0.0130 -0.6214 0.0704 -0.0267
0 0.6363 0 -0.1261 0.1032
0 0.1045 -0.2087 -0.2279 -0.4161
0 -0.1156 0.3497 0.7449 0
Dnew = 5×5
1.0000 0 0 0 0
0 3.8801 5.1046 0 0
0 -5.1046 3.8801 0 0
0 0 0 2.1199 0.7018
0 0 0 -0.7018 2.1199
Dnew 仍有实数特征值,但复共轭特征值被替换为 2×2 的块。