MATLAB学习笔记：方阵的相似对角化

>> A=[11 -6 4 -10 -4;-3 5 -2 4 1;-8 12 -3 12 4;1 6 -2 3 -1;8 -18 8 -14 -1];
>> [V D]=eig(A)

V =

-0.3244   -0.4983   -0.7759   -0.2343   -0.5752
0.1622    0.1878    0.0887   -0.0186    0.3890
0.6489    0.5633    0.2660    0.5059    0.3724
0.1622    0.0650   -0.5098    0.2903   -0.5919
-0.6489   -0.6284    0.2438   -0.7775   -0.1695

D =

3.0000         0         0         0         0
0    5.0000         0         0         0
0         0    5.0000         0         0
0         0         0    1.0000         0
0         0         0         0    1.0000

>> r=rank(V)

r =

5

A有5个线性无关的特征向量，所以矩阵A可以相似对角化，取P=V，则

P^(-1)*A*P=D

>> B=[-2 1 -2;-5 3 -3;1 0 2];
>> [V d]=eig(B)

V =

0.5774 - 0.0000i   0.5774 + 0.0000i  -0.5774
0.5773 - 0.0000i   0.5773 + 0.0000i  -0.5774
-0.5774            -0.5774             0.5773

d =

1.0000 + 0.0000i        0                  0
0             1.0000 - 0.0000i        0
0                  0             1.0000

>> r=rank(V)

r =

3

当A为实对称矩阵时，使用schur和eig可以求得正交的相似变换矩阵Q，使得Q^(-1)*A*Q为对角矩阵。

schur的调用格式为：

[Q,D]=schur(A) D是由A的特征值构成的对角矩阵，Q是正交矩阵，满足Q^(-1)*A*Q=D。

>> A=[0 0 0.32;0.2 0 0;0 0.72 0.95];
>> [P1 D1]=eig(A)

P1 =

0.6570             0.6570             0.3052
-0.0659 - 0.6075i  -0.0659 + 0.6075i   0.0613
-0.0476 + 0.4390i  -0.0476 - 0.4390i   0.9503

D1 =

-0.0232 + 0.2138i        0                  0
0            -0.0232 - 0.2138i        0
0                  0             0.9964

线性无关，可以为一组基。



求m趋于无穷时候，xm的极限：

>> syms c1 c2 c3 m;
>> xm=c1*D1(1,1)^m*P1(:,1)+c2*D1(2,2)^m*P1(:,2)+c3*D1(3,3)^m*P1(:,3);
>> limit(xm,m,inf)

ans =

0
0
0

因此，m趋于无穷时候，xm的极限为0.