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

Ready1998

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

SixPy

>>> import numpy as np
>>> A=np.mat([[2,-2,0],
              [-2,1,-2],
              [0,-2,0]])
>>> A
matrix([[ 2, -2,  0],
        [-2,  1, -2],
        [ 0, -2,  0]])
>>> np.all(A==A.T)
True

>>> D,V=np.linalg.eig(A) 
>>> D #D 对角线上的元素为A的特征值
array([ 4.,  1., -2.])
>>> V # V为相对应的特征向量所构成的矩阵。
matrix([[-0.66666667, -0.66666667,  0.33333333],
        [ 0.66666667, -0.33333333,  0.66666667],
        [-0.33333333,  0.66666667,  0.66666667]])

>>> np.diag(D) # 单位化
array([[ 4.,  0.,  0.],
       [ 0.,  1.,  0.],
       [ 0.,  0., -2.]])