從多維矩陣中提取塊

Question

考慮多維矩陣A，其中size(A)具有相同的偶數元素N。應該如何找到矩陣B與size(B)=size(A)/2這樣的：

B(1,1,...,1)=A(1,1,...,1), 
B(1,1,...,2)=A(1,1,...,2), 
... 
B(N/2,N/2,...,N/2)=A(N/2,N/2,...,N/2). 

Answer 1

我一般不喜歡arrayfun（或糊塗的功能），但是如果維數是不是在幾千，那麼這應該是蠻好的：

Nv = size(A)/2; 
S = arrayfun(@(x){1:x},Nv); 
B = A(S{:}); 

也應該使用不同尺寸的尺寸。只需決定如何處理尺寸mod(size(A),2)~=0。