每日值到多年月均值 Matlab

Question

我观察到了需要与生成的每月数据进行比较的每日数据，因此我需要获得三十年期间每个月的平均值。

我观察到的数据集目前是 365x31 的，每一天都是行（没有闰年！），额外的一列是月份数 (1-12)。

我遇到的问题是我似乎只能得到一个脚本来获得所有年份的平均值。 IE。我无法弄清楚如何让脚本分别为每一列执行此操作。数据示例如下：

1    12    14
1    -15   10
2    13    3
2    2     37
...all the way to 12 for 365 rows

SO：回顾一下，我需要得到 [12; -15； 13; 2] 然后 [14; 10个； 3； 37]等等。

我一直在尝试使用 unique() 函数进行循环，该函数可用于使行数平均但方法不正确。现在我需要它每个月（28-31 行）和列单独执行。结果应该是一个 12x30 矩阵。我觉得我错过了一些简单的东西。代码：

u = unique(m); %get unique values of m (months) ie) 1-12

for i=1:length(u)
   month(i) = mean(obatm(u(i), (2:31)); % the average for each month of each year
end

欣赏任何想法！谢谢！

Answer 1

您可以简单地过滤每个月的行，然后应用 mean，如下所示：

month = zeros(12, size(obatm, 2));
for k = 1:12
    month(k, :) = mean(obatm(obatm(:, 1) == k, :));
end

编辑：
如果你想要一些花哨的东西，你也可以这样做：

cols = size(obatm, 2) - 1;
subs = bsxfun(@plus, obatm(:, 1), (0:12:12 * (cols - 1)));
vals = obatm(:, 2:end);
month = reshape(accumarray(subs(:), vals(:), [12 * cols, 1], @mean), 12, cols)

看，妈，没有循环！