回 帖 发 新 帖 刷新版面

主题:跪求求解稀疏矩阵本征值和本征矢量的代码

我发现imsl库里面的子程序evcrg
CALL EVCRG (N, A, LDA, EVAL, EVEC, LDEVEC)  我运算的结果,其中有几个本征矢量和math.pdf上印出来的结果不一样,怎么回事?

回复列表 (共1个回复)

沙发

首先应该清楚求解是精确还是近似的、全部解还是子空间的。如果是近似的则有可能精度问题,高阶越高精度越底。如果是子空间的有丢失一个或几个特征值的可能性,尤其是后面的特征值。
其次,特征向量不唯一,特征向量乘以任意一个数均是正确解。通常人们为了得到“唯一”的解会做“正交归一”处理,但处理方法不唯一,须采用同一种方法处理才可得到唯一特征向量
此外,特征值有重根时特征向量不唯一(即使正交归一处理也不唯一),对应于相同特征值的所有特征向量的线性组合都是解

我来回复

您尚未登录,请登录后再回复。点此登录或注册