我有一个矩阵,每列代表一段时间内的特征。我需要找到给定窗口大小的这些值的移动平均值。
MATLAB中是否有类似于one的函数?
output = tsmovavg(vector, 's', lag, dim)
最佳答案
您可以使用FILTER函数。一个例子:
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中的IMFILTER和FSPECIAL
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);
请注意填充的区别:
output1(wndw:end)
对应于output3