matlab for循環：最快，最有效的方法來重現大矩陣

Question

我的數據是一個2096x252矩陣的double值。我需要一個for循環或執行以下操作的等效項：
每次矩陣被重現時，第一個數組被刪除，第二個數組成爲第一個。當循環再次運行時，剩餘的矩陣被再現，並且第一個數組被刪除，並且下一個變成第一個，依此類推。 我嘗試過使用repmat，但處理大型矩陣（2096x252）時太慢且乏味。

例輸入：

1 2 3 4 
3 4 5 6 
3 5 7 5 
9 6 3 2 

所需的輸出：

1 2 3 4 
3 4 5 6 
3 5 7 5 
9 6 3 2 
3 4 5 6 
3 5 7 5 
9 6 3 2 
3 5 7 5 
9 6 3 2 
9 6 3 2 

Answer 1

免責聲明：
你可能有大矩陣內存問題，但這不是問題。

現在，對業務：
對於給定矩陣A，與for循環的簡單方法是：

[N, M] = size(A); 
B = zeros(sum(1:N), M); 
offset = 1; 
for i = 1:N 
    B(offset:offset + N - i, :) = A(i:end, :); 
    offset = offset + size(A(i:end, :), 1); 
end 

B是所需的輸出矩陣。
但是，由於for循環，此解決方案預計也會變慢。

編輯：預分配B而不是動態更改大小（此優化應該實現輕微的加速）。

Answer 2

一般來說，使用Matlab預分配大型數組要比增量式構建要快得多。當你事先知道大陣的最終尺寸時，沒有理由不遵循這個一般性建議。

像下面這樣的東西應該做你想做的。假設你有一個數組in(nrows, ncols);那麼

indices = [0 nrows:-1:1]; 
out = zeros(sum(indices),ncols); 
for ix = 1:nrows 
    out(1+sum(indices(1:ix)):sum(indices(1:ix+1)),:) = in(ix:end,:); 
end 

這對您的小測試輸入有效。我希望你能弄清楚發生了什麼。

無論是我不知道的所有可能的方法中最快的方法，但我希望它比增量式構建大矩陣要快得多。