一般化固有値問題の部分空間反復法について質問があります。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 ループなしで部分空間法を使用するnnと、( matlab からの) 完全な固有値の最小値と比較すると、最小固有値は同じではありませんeig。と for ループを使用する
nnと、固有値が同じであることがわかります。 アルゴリズムにはこれらの追加の手順は含まれていませんが、これらがないと、matlab の結果は反復法と一致しませんでした。これには何か特別な理由がありますか？
前もって感謝します...