特征值NaN和inf

时间:2014-02-06 08:09:26

标签: python numpy linear-algebra

假设我有一个系统AX = nBX,其中A和B是已知的martrices,X是系数矩阵。

我使用切比雪夫多项式来解决这个问题。

BC是u(-1)= 0 = u(1)

我正在为矩阵A和B的第一行和最后一行强加bc。

e=solve(A,B)

e[1]=0

e[-1]=0

x=solve(A,e)

这有什么问题?

1 个答案:

答案 0 :(得分:2)

你似乎在问的问题是:为什么广义特征值问题有特征值inf和nan?

你的广义特征值问题是奇异的,并且具有特征值lambda = alpha / beta,使得(alpha = 0,beta = 0)和(alpha!= 0,beta = 0)。由于eigvals报告特征值,因此它们相应地为0/0 = nan或x / 0 = inf。

http://www.netlib.org/lapack/lug/node35.html

如果你的问题不应该有这样的特征值,那么矩阵的构造可能会出错。