真的不应该有任何区别这两个电话。你能把一些输入矩阵这发生在哪里?
eigs使用smallestabs vs标量
8视图(30天)
显示旧的评论
你好,
我注意到,在某些情况下使用
eigs
命令来解决一个广义特征值问题,获得的最小非零特征值及其对应的特征函数在使用“smallestabs”是复杂的。但是同样的问题(),当目标的最小非零特征值使用真正的标量,得到的特征值和特征向量是真实的。之间有不一致的原因eigs (a, B, k, smallestabs)和eigs (a, B, k,标量)针对同一eignevalue吗?
谢谢
接受的答案
克里斯汀Tobler
2022年2月18日
嗨,杰克,
我最初被误解,你得到不同的结果,当传入的smallestabs与传入数字标量0,本来eigs bug。现在我的问题是了解不同特征值计算有“smallestabs”(相当于0)之间通过与1.248(接近目标特征值)。
我不能重现看到1.248 + / - 0.0003我自己,但这可能是因为一些计算机有关舍入。σ的选择影响舍入误差,因为EIGS的第一步是计算的LU分解-σ* B,然后使用这个矩阵和线性系统解决的算法来确定特征值。
正因为如此,一个小小的改变会影响舍入错误,这可以相对较大的转变时接近一个特征值(0,因为一个特征值计算约e-12)。所以近奇异矩阵像,任何选择σ可能有所改善,甚至2为例,因为这意味着矩阵a -σ* B使用的不是接近奇异了。
没有一个真正返回真实和复杂数字系统的原因。特殊情况的对称输入定义的和对称的积极B,特征值总是真实的。但是对于一个真正的非对称像这里,特征值是真实的或复杂的共轭双。因为有两个特征值接近1.2483,有可能他们会分裂成一对共轭复数,这仍然是一个短的距离两个真正的价值。
