第八讲--相似理论

充分条件：1。A是实对称矩阵一定可以相似对角化，可以相似对角化不代表它是实对称矩阵。普通矩阵也有可能出现这种情况，只需要有n个线性无关的特征向量就行了。

                2。A有n个互异的特征值（必然有n个特征向量）好，那么A必然可以相似对角化。但是反过来，A可以相似对角化不代表他有n个互异的特征值，只代表它有n个线性无关的特征向量。（n个线性无关的特征向量对应的特征值是有可能等有可能不等的）

判断A是否可以相似对角化（        A=[1 1 1 1;1 1 -1 -1;1 -1 1 -1;1 -1 -1 1];        ）

1        求A的特征向量和特征值 （        [V,D]=eig(sym(A))        ）

                                                        （   [V,D]=eig(A)             ）

2        A的特征向量的秩（      rank(V)          ）

>> A=[1 1 1 1;1 1 -1 -1;1 -1 1 -1;1 -1 -1 1];
>> [V,D]=eig(sym(A));
>> rank(V)

ans =

     4

>> inv(V)*A*V
 
ans =
 
[-2, 0, 0, 0]
[ 0, 2, 0, 0]
[ 0, 0, 2, 0]
[ 0, 0, 0, 2]

 

>> A=[1 -1 -0.5;2 0.5 -2;-0.5 -1 1];
>> [V,D]=eig(A);
>> V

V =

   -0.4243    0.7071   -0.7071
    0.5657    0.0000    0.0000
   -0.7071    0.7071   -0.7071

>> rank(V)

ans =

     3

>> [V,D]=eig(sym(A));
>> V
 
V =
 
[ 3/5, 1]
[-4/5, 0]
[   1, 1]

poly(A)表示生成A的特征多项式，特征多项式的根是A的特征值。

>> p=poly(A)

p =

    1.0000   -2.5000    1.7500   -0.3750

>> y=poly2sym(p)
 
y =
 
x^3 - (5*x^2)/2 + (7*x)/4 - 3/8
 
>> r=roots(p)

r =

    1.5000
    0.5000
    0.5000

 好用：

A与B相似

 orth()可以用来得到单位正交基（就是线性代数的单位化）