使用numpy进行类似的矩阵计算

我正在尝试使用随机可逆矩阵P找到与3 X 3矩阵:A相似的矩阵B。 B = P_inv.A.P

import numpy as np
from scipy import linalg as LA
from numpy.linalg import inv

A = np.random.randint(1,10,9).reshape(3,3)
P = np.random.randn(3,3)
P_inv = inv(P)
eig1 = LA.eigvalsh(A)
eig1 = np.sort(eig1)

B1 = P_inv.dot(A)
B = B1.dot(P)
eig2 = LA.eigvalsh(B)
eig2 = np.sort(eig2)

print(np.round(eig1 ,3))
print(np.round(eig2,3))

但是,我注意到eig1和eig2永远不相等。 我缺少什么,还是数字错误?

谢谢

凯达