那位大神知道，如何使用python将一个实对称矩阵对角化呀？

Ready1998

那位大神知道，如何使用python将一个实对称矩阵对角化呀？求出的对角矩阵就是该实对称矩阵的特征值组成的矩阵

SixPy

1. >>> import numpy as np
   
2. >>> A=np.mat([[2,-2,0],
   
3.               [-2,1,-2],
   
4.               [0,-2,0]])
   
5. >>> A
   
6. matrix([[ 2, -2,  0],
   
7.         [-2,  1, -2],
   
8.         [ 0, -2,  0]])
   
9. >>> np.all(A==A.T)
   
10. True
    
11. 
    
12. >>> D,V=np.linalg.eig(A)
    
13. >>> D #D 对角线上的元素为A的特征值
    
14. array([ 4.,  1., -2.])
    
15. >>> V # V为相对应的特征向量所构成的矩阵。
    
16. matrix([[-0.66666667, -0.66666667,  0.33333333],
    
17.         [ 0.66666667, -0.33333333,  0.66666667],
    
18.         [-0.33333333,  0.66666667,  0.66666667]])
    
19. 
    
20. >>> np.diag(D) # 单位化
    
21. array([[ 4.,  0.,  0.],
    
22.        [ 0.,  1.,  0.],
    
23.        [ 0.,  0., -2.]])
    
24. 
    

复制代码