我有一个关于广义特征值问题的子空间迭代方法的问题。我正在使用 MATLAB 使用子空间迭代方法求解一些最低特征值。求解简化系统后，我们是否对质量矩阵的特征向量进行归一化？
例如：
[eve,eva] = eig(full(kred),full(mred));
x = xbar*eve % update eigen vector
在x更新之前，我们应该执行
nn = eve'*mred*eve
for i = 1:min(2*m,m+8)
evec(:,i) = evec(:,i)/sqrt(i,i);
end
然后更新x？当我使用没有和 for 循环的子空间方法时，与完整特征值中的最低值（来自 matlab ）nn相比，最低特征值是不一样的。eig当我使用
nnand 循环时，我看到特征值是相同的。 该算法不包括这些额外的步骤，但没有这些，matlab 的结果与我的迭代方法不匹配。这背后有什么特别的原因吗？
提前致谢...