密集广义厄米特不定特征值问题

计算科学 本征系统 拉帕克 特征值
2021-11-26 21:03:58

Lapack 包含一个驱动例程,用于解决以下形式的密集广义 Hermitian 正定特征值问题Ax=λBx, 在哪里AB都是厄米特式的,并且B是肯定的。我想知道是否有一种方法(代码或修改zhegv)时B是不确定的(具体来说,B几乎可以肯定是非奇异的)。我知道驱动程序首先执行正定的 Cholesky 分解B, 并将其转换为标准特征值问题,但这仍然有效B是无限期的吗?我宁愿使用专门的方法而不是完全通用的求解器,因为我希望保证特征向量的特征值和正交关系的真实性。

1个回答

什么时候B是不确定的,特征值可能很复杂,利用对称性几乎没有优势。相反,通常使用 QZ 算法。

编辑:如果您先验地知道所有特征值都是真实的,那可能是因为您可以先验地建立一些线性组合C=sA+tB是肯定的。在这种情况下,您应该考虑您的问题A或者B取而代之C. 它具有相同的特征向量,并且特征值通过莫比乌斯变换相关。