Optimizing calculation of eigenvectors and eigenvalues
显示 更早的评论
I have a quadratic matrix A with a size of about 2000x2000. I want to calculate its eigenvectors and eigenvalues. The eigenvalues must be in a vector (not a diagonal matrix as usual). I use the following code:
eigenvals = abs(real(eig(A)));
[eigenvecs, ~] = eig(A);
eigenvecs = real(eigenvecs);
The problem is that it takes a lot of time (about 15 seconds) and I need to repeat this process a lot of times in a loop. So I tried to optimize it and changed it to this:
[eigenvecs, eigenvals]=eig(A);
eigenvecs=real(eigenvecs);
eigenvals=abs(real(diag(eigenvals)));
This runs about 5 seconds faster. The calculated eigenvectors are the same for both code snippets, however the eigenvalues are not and differ a bit. How could I further optimize the code in a way that it delivers equivalent results to the first code snippet?
Thank you very much in advance!
采纳的回答
David Goodmanson
2021-1-16
编辑：David Goodmanson
2021-1-16
Hi Guiseppi,
there is nothing guaranteed about the order of the eigenvalues that are produced by eig. And if you do the first method, which is basically
lambda = eig(A)
% and
[v, ~] = eig(A)
the order of eigenvalues may not be the same in each case, meaning that the order of eigenvalues and the order of rows of the eigenvector matrix may not match up. In fact for some examples I tried, they do not match up, which is similar to your experience. This means the first method does not work.
If you would like a vector rather than a diagonal matrix for lambda, there is the option
[v lambda] = eig(A,'vector')
which only calls eig once and so is faster than the first method, in addition to being correct.
1 个评论
which only calls eig once and so is faster than the first method
Also, with
[v, ~] = eig(A)
a 2000x2000 matrix memory allocation for the second output argument is still executed by eig, even though it is promptly discarded in the caller workspace because of '~'. Using the form eig(A,'vector') avoids that extra memory allocation step.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Linear Least Squares 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
