Construct Power of Matrix without for loop
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Hi, everyone:
Suppose I have a 2 by 2 matrix A, if I want construct a larger matrix B that is defined as:
B=[A, A^2, A^3, A^4, ... A^N];
is it possible to do it without for loop?
Thanks
0 Kommentare
Akzeptierte Antwort
Azzi Abdelmalek
am 11 Aug. 2013
Bearbeitet: Azzi Abdelmalek
am 11 Aug. 2013
Example
A=magic(2);
N=3;
B=cell2mat(arrayfun(@(x) A^x,1:N,'un',0))
0 Kommentare
Weitere Antworten (2)
Jan
am 12 Aug. 2013
Azzi's suggestion is fine for N=3. If you are talking about larger N, neither the repeated power operator nor arrayfun nor cell2mat are efficient:
N = 10000
A = rand(2, 2);
tic;
B = cell2mat(arrayfun(@(x) A^x,1:N,'un',0));
toc
tic;
B = zeros(2, 2, N);
P = 1;
for k = 1:N
P=P*A;
B(:,:,k) = P;
end
B = reshape(B, 2, N * 2);
toc
Elapsed time is 0.199776 seconds.
Elapsed time is 0.039983 seconds.
So I'd prefer the more efficient FOR loop.
2 Kommentare
Ming
am 11 Aug. 2013
Bearbeitet: Ming
am 11 Aug. 2013
1 Kommentar
Walter Roberson
am 11 Aug. 2013
Note: arrayfun() just hides the "for" loop. If you are willing to use it, then Azzi's example does what you ask.
Siehe auch
Kategorien
Mehr zu Loops and Conditional Statements finden Sie in Help Center und File Exchange
Produkte
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!