More non-zero generalized eigenvalues than in theory

Question

Fernando García-García 2014-10-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/160587-more-non-zero-generalized-eigenvalues-than-in-theory

评论： Fernando García-García 2014-10-29

采纳的回答： Fernando García-García

在 MATLAB Online 中打开

Dear all,

I'm trying to compute generalized eigenvalues. First I've tried:

[V,D]=eig(A,B);

Other attempts

    [V,D]=eig(A,B,'chol');
    [V,D]=eig(A,B,'qz');
    [V,D]=eig(A/B);

Well, my problem is that in all cases I obtain more non-zero eigenvalues than in theory. A,B are both of 42x42 size and rank(A)=2, rank(B)=41. (By the way, in case it's useful: A and B are between- and within-scatter matrices from a k=3 multiclass Linear Discriminant Analysis, i.e. three classes/groups).

Theoretically, since rank(A)=2, I would get two non-zero eigenvalues... but results are: [3.1055, 0.9127, 0.7718, etc.] All other are indeed <1e-12. I've also tried computing the scatter matrices of my data having been first transformed to zero-mean unit-stdev, but yet I get three non-negligible eigenvalues!

I don't know where the problem is. I've verified my code (mainly to check if my computation of matrices is correct) against a given example of size 6x6, and I do get two non-zero eigenvalues as expected.

Thank you very much in advance! Best regards,

Fernando

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Fernando García-García 2014-10-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/160587-more-non-zero-generalized-eigenvalues-than-in-theory#answer_157040

在 MATLAB Online 中打开

I think I've found the solution for this issue by using the intermediate output of the QZ algorithm.

[AA,BB,Q,Z,V,W] = qz(A,B,'real'); % in my case I think 'real' is safe, still needs a check

Right generalized eigenvectors are rows of W transposed. Eigenvales are:

D(i,i) = AA(i,i)/BB(i,i);

I guess that the code in eig() doesn't check if BB(i,i) has an aprox. 0 value, and numerical precision due to a 0/0 made an extra non-zero eigenvalue appear in my example data.

Let's see if I'm wrong or if anyone has a different explanation. Hope this helps!

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Fernando García-García 2014-10-29

在 MATLAB Online 中打开

Or most problably, eig() does check if approx. 0... but in my problem (with a considerable matrix dimension but low rank), I suppose I needed a greater tolerance than the default. So I implemented by hand the

D(i,i) = AA(i,i)/BB(i,i);

checking if B(i,i)>TOL. I took TOL as big as 1e-6, because values in the order of 1e-7 appear in B's diagonal for my case.

请先登录，再进行评论。