Question

I have a question. Suppose I have matrix A =

 1     2     3
 4     5     6
 7     8     9
10    11    12

I need to select n rolling rows from A and transpose elements in new matrix C in rows. The loop that I use is:

n = 3;     %for instance every 3 rows of A

B = []; 

for i = 1:n 

    Btemp = transpose(A(i:i+size(A,1)-n,:)); 

    B = [B;Btemp]; 

end

C=B ;

并产生矩阵C,即:

C =

 1     2     3     4     5     6     7     8     9
 4     5     6     7     8     9    10    11    12

这就是我也想要做的, 但我可以做同样的工作没有循环?

计算一个3280x35大小的A矩阵需要4分钟。

Answer 1

我认为如果你进行初始化,你可以非常快地让它工作。另外一种技巧是首先进行转换,因为MATLAB使用列作为第一索引,而不是行。

tic
A =  reshape(1:3280*35,[3280 35]) ; %# Generate an example A
[nRows, nCols] = size(A);

n = 3; %for instance every 3 rows of A
B = zeros(nRows-n+1,nCols*n);
At = A ;
for i = 1:size(B,1)
    B(i,:) = reshape(At(:,i:i+n-1), [1 nCols*n]); 
end
toc

时间已经过去

Elapsed time is 0.004059 seconds.

Answer 2

我不会在循环圈中使用重塑,而是将A第一行转换为单行(实际上一列也会有效,但无关紧要)

Ar = reshape(A ,1,[]); % the   is important here!

然后从 Ar 中选择元素非常简单 :

[nrows, ncols] = size(A);
new_ncols = ncols*n;
B = zeros(nrows-(n-1),new_ncols);

for ii = 1:nrows-(n-1)
    B(ii,:) = Ar(n*(ii-1)+(1:new_ncols));
end

然而,B 的预先分配使你们有了最大的改进:在http://www.mathworks.nl/help/techdoc/matlab_prog/f8-784135.html

Answer 3

我现在没有马特拉布, 但我认为你可以做到这一点没有这样的循环:

reshape(permute(cat(A(1:end-1,:),A(2:end,:),3),[3,2,1]), [2, size(A,2)*(size(A,1) - 1)]);

这难道不能做你们所欲做的事吗?

A1 = A(1:end-1,:);
A2 = A(2:end,:);
answer = [A1(:) ; A2(:)]

友情链接