我正在解决具有逆复数移位的广义特征值问题：

$$(M-\sigma J)^{-1}J \boldsymbol{x} = \boldsymbol{x} \nu \enspace .$$

我的矩阵是通过 Getfem++ 的有限元离散化获得的。雅可比取决于一些参数。我在计算矩阵和的同一 C++ 代码中使用 Matlab 和直接使用 Arpack 计算了特征值。当我使用 Matlab 计算特征值时，我使用 Getfem++ 加载先前以市场矩阵格式保存和 ，并使用函数. 现在，在某些情况下，Arpack 特征值和 Matlab 特征值之间的差异非常大。特别是，当我的参数值使得我的特征值几乎为零时。使用 Arpack，实部大约为$J$$M$$J$$M$$J$eigs$10^{-10}$与 Matlab 的数量级为。此外，在这些情况下，Matlab 给了我一些关于矩阵病态的警告。你有没有遇到过这样的问题？$10^{-4}$

PS：我也尝试更改 eigs 的默认选项，但我总是得到相同的结果！