A lot of eigenvalues

2 次查看(过去 30 天)
tomas
tomas 2011-8-23
Hello everybody!
I have many 3x3 matrixes (e.g. 3 000 000 matrixes) and I need to compute their eigenvalues. Nowadays, I'm doing it in loop having 3 000 000 counts, but this way is very time consuming. I'd like to know if there is any possibility how to replace the loop by some matrix operation. For example, I can make a 3x3x3000000 matrix from which I would get the eigenvalues, somehow?
Can anybody help with this?
Thank you in advance,
Tom
  1 个评论
Oleg Komarov
Oleg Komarov 2011-8-23
c(1:3e6) = {rand(3)};
tic
egv = cellfun(@eig,c,'un',0);
toc
23 seconds to me it seems fair.

请先登录,再进行评论。

采纳的回答

Doug Hull
Doug Hull 2011-8-23
This from Oleg seems like an answer, not a comment. I prefer to see questions in the 'answered' state.
c(1:3e6) = {rand(3)};
tic
egv = cellfun(@eig,c,'un',0);
toc
23 seconds to me it seems fair.
  3 个评论
tomas
tomas 2011-8-24
Olegs answer works. Although I don't understand last two input parameters of cellfun function - 'un',0. Can you tell me what they mean, please?
Thank you guys.
Daniel Shub
Daniel Shub 2011-8-24
They are shorthand for 'uniform' and false

请先登录,再进行评论。

更多回答(1 个)

Daniel Shub
Daniel Shub 2011-8-24
While the loop is ~20% slower on my machine then using cellfun, the answer is in a double array and not a cell array. There may be hidden costs of getting the information from the cell array into a usable form.
N = 3e6;
x = rand(3, 3, N);
tic
egv = zeros(3, N);
for ii = 1:N
egv(:, ii) = eig(x(:, :, ii));
end
toc

类别

Help CenterFile Exchange 中查找有关 Loops and Conditional Statements 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by