为什么 np.linalg.eig(L) 不能得到正确的结果?

问题描述

我们都知道拉普拉斯矩阵的最小特征值为零,特征向量为[1 1 ... 1],因为每一行的和为零。 所以对于矩阵:[[ 1. -1. 0.][-1. 2. -1.][ 0. -1. 1.]]

我尝试了 np.linalg.eig(L)np.linalg.eigh(L)scipy.linalg.eig(L),但它们中的每一个都无法得到零和 [1 1 ... 1]

np.linalg.eig(L)中:

(-3.367702055640532e-17,0.9999999999999998,2.999999999999999)
(array([-0.40824829,0.70710678,0.57735027]),array([8.16496581e-01,2.61239546e-16,5.77350269e-01]),array([-0.40824829,-0.70710678,0.57735027]))

np.linalg.eigh(L)中:

(9.996582244115599e-17,1.0,3.0)
(array([-0.57735027,0.40824829]),array([-5.77350269e-01,7.58447699e-17,-8.16496581e-01]),array([-0.57735027,0.40824829]))

scipy.linalg.eig(L)中:

((-3.367702055640532e-17+0j),(0.9999999999999998+0j),(2.999999999999999+0j))
(array([-0.40824829,0.57735027]))

我知道计算机表示错误,但特征向量远非 [1 ... 1]^T

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)