每n行矩阵求和

Question

有什么方法可以总结矩阵中每组三行的列值？
我可以手动总结三行。

例如

% matrix is the one I wanna store the new data.
% data is the original dataset.
matrix(1,1:end) = sum(data(1:3, 1:end))
matrix(2,1:end) = sum(data(4:6, 1:end))
...

但如果数据集很大，这将行不通。
有没有什么方法可以在没有循环的情况下自动执行此操作？

Answer 1

Prashant 之前回答得很好，但我有一个简单的修改：

fl = filterLength;
A = yourVector (where mod(A,fl)==0)
sum(reshape(A,fl,[]),1).'/fl;

即使在 fl==1（原始值）时，也有使行运行的“,1”。 我在这样的 for 循环中运行它时发现了这一点：

... read A ...
% Plot data
hold on;

averageFactors = [1 3 10 30 100 300 1000];
colors = hsv(length(averageFactors));
clear legendTxt;


for i=1:length(averageFactors)
% ------ FILTERING ----------
clear Atrunc;
clear ttrunc;
clear B;
fl = averageFactors(i); % filter length
Atrunc = A(1:L-mod(L,fl),:);
ttrunc = t(1:L-mod(L,fl),:);

B = sum(reshape(Atrunc,fl,[]),1).'/fl;
tB = sum(reshape(ttrunc,fl,[]),1).'/fl;
length(B)
plot(tB,B,'color',colors(i,:) )
%kbhit ()
endfor

Answer 2

还有其他四种方法：

1. 必填for-loop：

   % for-loop over each three rows
   matrix = zeros(size(data,1)/3, size(data,2));
   counter = 1;
   for i=1:3:size(data,1)
       matrix(counter,:) = sum(data(i:i+3-1,:));
       counter = counter + 1;
   end
   

2. 使用mat2cell 进行平铺：

   % divide each three rows into a cell
   matrix = mat2cell(data, ones(1,size(data,1)/3)*3);
   
   % compute the sum of rows in each cell
   matrix = cell2mat(cellfun(@sum, matrix, 'UniformOutput',false));
   

3. 使用第三维度（基于this）：

   % put each three row into a separate 3rd dimension slice
   matrix = permute(reshape(data', [], 3, size(data,1)/3), [2 1 3]);
   
   % sum rows, and put back together
   matrix = permute(sum(matrix), [3 2 1]);
   

4. 使用accumarray：

   % build array of group indices [1,1,1,2,2,2,3,3,3,...]
   idx = floor(((1:size(data,1))' - 1)/3) + 1;
   
   % use it to accumulate rows (appliead to each column separately)
   matrix = cell2mat(arrayfun(@(i)accumarray(idx,data(:,i)), 1:size(data,2), ...
       'UniformOutput',false));
   

当然，到目前为止的所有解决方案都假设行数可以被3 整除。

Answer 3

这个单行reshapes 使特定单元格所需的所有值都在一个列中，然后sum，然后reshapes 回到预期的形状。

reshape(sum(reshape(data, 3, [])), [], size(data, 2))

如果您想将不同数量的行加在一起，可以更改裸露的3。您需要确保每组中的行数平均分配。

Answer 4

将矩阵分成三部分并将它们相加：

matrix = data(1:3:end, :) + data(2:3:end, :) + data(3:3:end, :);

如果size(data,1) 不是三的倍数，这将产生错误，因为这三个部分的大小不同。如果适合您的数据，您可以通过截断 data 或在末尾附加一些零来解决此问题。

您还可以使用 reshape 和 3D 数组来做一些花哨的事情。但我更喜欢上面的（除非你需要用变量替换3...）