python – 如何计算最小的非零特征值

发布时间：2020-12-20 13:10:55 所属栏目：Python 来源：网络整理

导读：给定正半定矩阵M,我想找到其最小的非零特征值.在 python中,这段代码看起来很诱人 import numpy as np(w,v) = np.linalg.eigh(M)minw = np.amin(w)if (np.isclose(minw,0) and minw 0): print M,minw 这是一个小输入矩阵示例. [ 6 2 -4 -2][ 2 6 0 -6][-4 0 6

给定正半定矩阵M,我想找到其最小的非零特征值.在 python中,这段代码看起来很诱人

import numpy as np
(w,v) = np.linalg.eigh(M)
minw = np.amin(w)
if (np.isclose(minw,0) and minw > 0):
    print M,minw

这是一个小输入矩阵示例.

[ 6  2 -4 -2]
[ 2  6  0 -6]
[-4  0  6  0]
[-2 -6  0  6]

不幸的是,如果你试试这个,你会得到8.90238403828e-16.我一般不知道如何判断非常小的数字是否为零.

如何找到矩阵的最小非零特征值(并确保它确实非零)？

解决方法

这可能是函数中的浮动问题.解决方案将取决于问题的严重性.您可以尝试循环结果并查找最小值,例如,函数( numpy)：

np.isclose(a,b)

,给出容差内两个值之间的相似性.这不是一个干净的解决方案,但它通常被认为是一个安全的比较.至于“eigh”代码本身或许有些收敛问题,但我实在说不出来.

（编辑：李大同）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!