我将此函数(https://stackoverflow.com/a/31207520/8948867)应用于3D数组(5x5x5),以获得该体积的3x3粗麻布矩阵。所得的粗麻布矩阵的形状为3x3x5x5x5。
我的目的是计算2维hessian矩阵的3个特征值,但是由于我的hessian矩阵具有2个以上的维度(3x3x5x5x5),因此计算出的特征值的形状为3x3x5x5。
我不知道问题出在粗麻布矩阵中还是在估计特征值的函数中,或者我可能错过了两者之间的某些步骤。
这是我使用的代码(hessian是上面链接中显示的函数):
hss = hessian (arr)
from numpy import linalg as LA
e_val, e_vec = LA.eigh (hss)