假设我有一个系统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)
这有什么问题?
答案 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
如果你的问题不应该有这样的特征值,那么矩阵的构造可能会出错。