1

我对 matlab 很陌生，我正在尝试找到一种不使用 for 循环来完成以下任务的方法：

我有一个如下所示的数据集：

data = [1 5; 1 3; 1 8; 2 1; 2 2; 2 5; 3 3; 3 8; 3 4]

第一列是一个组（将来是一个月年的组合）

现在我想计算具有给定窗口大小的第二列的总和，但前提是组索引相同 - 如果不是，则应计算该组中的最大总和。

使用 window size=2 我想创建以下结果：

summed_data = [1 8; 1 11; 1 8; 2 3; 2 7; 2 5; 3 11; 3 12; 3 4]

窗口大小为 3 时，结果将如下所示：

summed_data = [1 16; 1 11; 1 8; 2 8; 2 7; 2 5; 3 15; 3 12; 3 4]

等等。

我曾考虑通过创建足够的子索引来使用 accumarray - 但我对窗口大小有疑问，并且总和是重叠的。

有没有人知道如何在不使用循环的情况下实现它？

在此先感谢和最好的问候斯蒂芬

4

1 回答 1

4

这似乎有效：

ws = 2; k = [ones(ws,1);zeros(mod(ws,2),1)];
C = accumarray(data(:,1),data(:,2),[],@(v){conv(v,k,'same')})

您似乎锚定到窗口中的当前像素并向前看。那是对的吗？

不确定这是否涵盖所有极端情况，但它可能会引导您朝着正确的方向前进。

测试：ws = 2

ws = 2; k = [ones(ws,1);zeros(mod(ws,2),1)];
C = accumarray(data(:,1),data(:,2),[],@(v){conv(v,k,'same')});
summed_data = [data(:,1) vertcat(C{:})]
summed_data =

     1     8
     1    11
     1     8
     2     3
     2     7
     2     5
     3    11
     3    12
     3     4

测试：ws = 3

summed_data = [data(:,1) vertcat(C{:})]
summed_data =

     1    16
     1    11
     1     8
     2     8
     2     7
     2     5
     3    15
     3    12
     3     4
于 2014-10-31T17:24:52.993 回答