在 Scipy 中，我有一个 7x7 矩阵，其中有一些条目大到 10^22。
[[7.87720699e+21 5.81984000e+12 4.06283195e+18 1.36914426e+15
6.84296262e+21 1.18920842e+20 2.30532710e+21]
[5.81984000e+12 6.32869618e+03 3.12159108e+09 8.49240505e+05
5.05790011e+12 8.25020510e+10 1.71012096e+12]
[4.06283195e+18 3.12159108e+09 2.10517370e+15 6.90360153e+11
3.52958137e+18 6.08623401e+16 1.18959479e+18]
[1.36914426e+15 8.49240505e+05 6.90360153e+11 2.68568446e+08
1.18907574e+15 2.15022766e+13 3.99729760e+14]
[6.84296262e+21 5.05790011e+12 3.52958137e+18 1.18907574e+15
5.94451366e+21 1.03298010e+20 2.00265814e+21]
[1.18920842e+20 8.25020510e+10 6.08623401e+16 2.15022766e+13
1.03298010e+20 1.81926760e+18 3.47747090e+19]
[2.30532710e+21 1.71012096e+12 1.18959479e+18 3.99729760e+14
2.00265814e+21 3.47747090e+19 6.74706499e+20]]
当我打电话给spicy.linalg.eigvalsh时，我得到
scipy.linalg.eigvalsh(square_matrix)
array([-1.65239967e+05, 3.55247340e+04, 2.64944833e+06, 2.26542682e+09,n 2.01743752e+14, 5.56910661e+16, 1.44981926e+22])
这不好，因为矩阵是对称的正半定的(它的形式是AA^T(，所以它应该具有所有非负特征值。当我只问它最小的特征值时，它得到的答案与上次完全不同。
scipy.linalg.eigvalsh(square_matrix, eigvals = (0,0))
array([-464577.85826165])
这是怎么回事？只是数字太大以至于在计算中造成了巨大的错误吗？
感谢您的帮助，我对这样的数值实验有些陌生。
编辑：这是另一个矩阵，以列表形式发布：
[[178484429459288.62, 1.262534539362581e+16, 5756113437609551.0, 1.5274842899247696e+17, 6343247145.960138, 3119899812227.452, 2.6847705451184886e+17], [1.262534539362581e+16, 9.311154944318847e+17, 4.25182007604871e+17, 1.0968899948471126e+19, 418362473621.39294, 236566750583994.28, 1.929426966567391e+19], [5756113437609551.0, 4.25182007604871e+17, 1.9416578714719302e+17, 5.00372260140452e+18, 190262613997.92053, 108142802848362.1, 8.801791711005543e+18], [1.5274842899247696e+17, 1.0968899948471126e+19, 5.00372260140452e+18, 1.3144272477543575e+20, 5288261203200.205, 2737268307603631.5, 2.3109381789197917e+20], [6343247145.960138, 418362473621.39294, 190262613997.92053, 5288261203200.205, 258730.2206841113, 99079200.39025305, 9282686524313.389], [3119899812227.452, 236566750583994.28, 108142802848362.1, 2737268307603631.5, 99079200.39025305, 61292039604.753365, 4817296126969346.0], [2.6847705451184886e+17, 1.929426966567391e+19, 8.801791711005543e+18, 2.3109381789197917e+20, 9282686524313.389, 4817296126969346.0, 4.0629951073761985e+20]]
scipy.linalg.eigvalsh(square_matrix, eigvals = (0,0))
> array([-9362.07065027])
scipy.linalg.eigvalsh(square_matrix)
> array([3.04005337e+01, 1.06920150e+04, 9.94140064e+06, 3.42415026e+09,n 2.27074109e+13, 1.99821497e+16, 5.38847690e+20])
我已经确认我在 Scipy 1.1.0 中。这些是在脚本中间调试期间获得的，如果这有什么不同的话。使用 VS Code 作为 IDE。
矩阵的条目上升到1e22但其最小的特征值可能要小得多，可能是1000左右。当使用双精度数字(每一步1e-16相对误差(执行复杂计算时，尤其是对近似特征值的迭代计算，得到大小1e-12的相对误差也就不足为奇了。从绝对值来看，它变得1e10，这完全使最小的特征值相形见绌，并且很容易使其为负。
搜索 k 个最小或最大的特征值是一项不同的任务：SciPy 不仅计算所有特征值并截断输出。在源代码中，根据参数的存在，可以看到对基于 LAPACK 的代码的不同调用eigvals我相信该参数在 LAPACK 内部也会有所不同。所有这些都意味着您得到的随机错误可能会有所不同，具体取决于参数eigvals.
我正在尝试估计条件数
我建议为此使用np.linalg.cond。
一般来说，处理 A 的奇异值比处理 A 的特征值乘以其转置(条件数几乎平方
(更好。