d*****y
发帖数: 140 | 1 在用eigs(A,10,'la')计算A的最大的10的特征值的时候，A是对称矩阵，结果中会有许
多0，而且0明显不是A的特征向量，这是咋回事？ | t***s
发帖数: 4666 | 2 show us your matrix.
my guess is that your matrix has many negative eigenvalues and is almost
sigular.
【在 d*****y 的大作中提到】
: 在用eigs(A,10,'la')计算A的最大的10的特征值的时候，A是对称矩阵，结果中会有许
: 多0，而且0明显不是A的特征向量，这是咋回事？
| d*****y
发帖数: 140 | 3 数据比较大，大部分是negative eigenvalues, 用eig算出来最接近于0的eigenvalues是
-0.0139， 0.0128， 0.0305 and 0.0443 算是almost sigular么？
【在 t***s 的大作中提到】
: show us your matrix.
: my guess is that your matrix has many negative eigenvalues and is almost
: sigular.
| d*****y
发帖数: 140 | 4 我觉得可能不是negative和sigular的问题，
我的矩阵A最负的eigenvalue有个理论的lower bound (假设是b),我用eigs分解 2b*I+A
, (I is identity matrix) 这样2b*I+A non-sigular,所有eigenvalues都是正的.
eigs分解出来的前10个最大的eigenvalues还是有很多0。
【在 t***s 的大作中提到】
: show us your matrix.
: my guess is that your matrix has many negative eigenvalues and is almost
: sigular.
| t***s
发帖数: 4666 | 5 don't know what's wrong. can you provide a minimum example that
exhibits your problem?
the following example works as expected for me.
a = randn(100);
b = a*a.';
c = -300*eye(100)+b;
eigs(c,10,'la')
+A
【在 d*****y 的大作中提到】
: 我觉得可能不是negative和sigular的问题，
: 我的矩阵A最负的eigenvalue有个理论的lower bound (假设是b),我用eigs分解 2b*I+A
: , (I is identity matrix) 这样2b*I+A non-sigular,所有eigenvalues都是正的.
: eigs分解出来的前10个最大的eigenvalues还是有很多0。
|
|
