1、给出一个3*3的可逆矩阵：import numpy as npa=[[1,2,3],[0,2,1],[3,1,3]]
2、计算a的特征向量：c=np.linalg.eig(a)d=c[1]
3、假设d的逆矩阵是d'，那么d'.a.d就是一个对角矩阵：e=np.dot(np.dot(np.linalg.inv(d),a),d)
4、运行结果看起来不是对角矩阵，这是为何？原来，这是精确度的问题。如果矩阵中的每个数字都保留两位有效数字，对角矩阵就显出了原形。e=np.around(e, decimals=2, out=None)
5、矩阵a=[[1,2,3,3],[0,2,1,2],[3,1,3,1],[3,5,6,9]]的对角化是：[缪梨痤刻[12.56687 0. -0. -0. ][ 0. -0.82664 0. 0. ][-0. -0. 2.36388 -0. ][-0. 0. -0. 0.89589]]
