A = np.mat(np.zeros((26, 6)))
for i in range(6):
A[:,i] = (sqrt(eig_value[i]) * eigvector[:,i])
运行之后出现错误,查了一下发现是等式右边的数组形式为(26,)
因此需要用reshape()重新定义数组形式
即:
A = np.mat(np.zeros((26, 6)))
for i in range(6):
A[:,i] = (sqrt(eig_value[i]) * eigvector[:,i]).reshape(26,1)```