import numpy as np
a = np.matrix([[0, 1, 0, 1, 0], [1, 0, 1, 0 ,0],[0, 1, 0, 0, 1], [1, 0, 0, 0, 1],[0,0, 1, 1,0]])
for i in range(20):
if i == 0:
t = a*a
continue
t = t * a
#t为特征值, Q为特征向量
t, Q = np.linalg.eig(a)
#以t为对角元素建立对角元素
T = np.diag(t)
A2 = Q * (T * T * T) * Q**-1
print(a**3)
print(A2)
Python求矩阵特征值及特征向量和实现对角化矩阵简化矩阵n次方运算
最新推荐文章于 2023-09-02 14:21:35 发布