accuracy of two methods to solve eignvalues/eignvector problem
显示 更早的评论
I am trying to solve an eignvalue/eignvector problem related to the multi-class LDA. The equation is inv(Sw)*Sb. I have two solutions:
[E,D]=eig(inv(Sw)*Sb)
or
[E,D]=eig(Sb,Sw);
when I try them, I get different eignvectors and eignvalues in both cases. Am I supposed to get the same results and what is more accurate to use the first way or the second one?
采纳的回答
Christine Tobler
2020-4-3
The second way is typically more accurate, because it doesn't incur additional round-off error when conputing inv(Sw) and when multiplying this to Sb.
Can you give an example of the eigenvalues and eigenvectors being different? It's to be expected that the order of the eigenvalues changes between two different ways of calling EIG, and the eigenvectors may be scaled in different ways between the two cases. Any other changes would be due to round-off error.
One way to check the accuracy of the combined eigenvectors and eigenvalues is to compute the residual norm:
For [U, D] = eig(A)
norm(A*U - U*D)
For [U, D] = eig(A, B)
norm(A*U - B*U*D)
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Define Shallow Neural Network Architectures 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
