11

我有一个矩阵,每列代表一个随时间变化的特征。我需要找到给定窗口大小的这些值的移动平均值。

有没有类似 MATLAB 中的函数

output = tsmovavg(vector, 's', lag, dim)
4

1 回答 1

21

您可以使用过滤器功能。一个例子:

t = (0:.001:1)';                                %#'
vector = sin(2*pi*t) + 0.2*randn(size(t));      %# time series

wndw = 10;                                      %# sliding window size
output1 = filter(ones(wndw,1)/wndw, 1, vector); %# moving average

甚至使用Image Package中的IMFILTERFSPECIAL

output2 = imfilter(vector, fspecial('average', [wndw 1]));

最后一种选择是使用索引(不推荐用于非常大的向量)

%# get indices of each sliding window
idx = bsxfun(@plus, (1:wndw)', 0:length(vector)-wndw);
%'# compute average of each
output3 = mean(vector(idx),1);

请注意padding的区别:output1(wndw:end)对应output3

于 2010-06-24T23:46:07.950 回答