下一节:Jacobi-Davidson方法 上一级:二次特征值问题 上一节:QEP的Cayley变换

线性化问题求解的数值方法

如第§9.2.2节所述, 可以使用广义特征值问题 (9.4) 与 (9.5),或者广义特征值 问题 (9.9)与 (9.10) 来求解相应的二次特征值问题 (9.1)。

如果所有矩阵 MCK 都是厄米的,并且 M 是正定的, 正如特殊情况 (9.2)所示, 那么选择就归结为采用本质上非厄米的广义特征值问题 (9.4) 与 (9.5),其中 B 矩阵是厄米正定的, 或者采用广义特征值问题 (9.9)与 (9.10), 其中 AB 矩阵都是厄米的,但它们都不是正定的。

8章讨论的数值方法 可以用于求解这些广义的“线性”特征值问题。例如,在MSC/NASTRAN [274]中, 采用了非厄米形式 (9.4)与 (9.5)。

第§8.6节讨论的对称不定Lanczos方法专门针对 广义对称不定特征值问题 (9.9)。 潜在的问题是在对称不定Lanczos方法中,基向量是相对于不定内积正交的。因此, 这些基向量可能不是线性独立的,算法可能会崩溃并导致数值不稳定。 尽管如此,由于潜在的内存需求和浮点运算节省,这通常是求解原始QEP的一种吸引人的方法。 详见第§8.6节。



下一节:Jacobi-Davidson方法 上一级:二次特征值问题 上一节:QEP的Cayley变换
Susan Blackford 2000-11-20