3

我没有这个包nlfilter，我也没有完全按照这个例子。

我有一个非常简单的函数fun，我想将它应用于数组的移动窗口。数组是Nx1，我想看看长度k间隔，比如说。所以对于N=10和我会k=3得到fun = @(x) min(x);

A = [13 14 2 14 10 3 5 9 15 8];

filter(A,k,fun) = [2 2 2 3 3 3 5 8];

在这里，我只想查看索引 1,2,3 然后是 2,3,4 然后 ... 然后是 8,9,10，因此最终序列的长度为 7。我可以使用 for 循环轻松完成此操作，但我不知道如何为 Matlab 向量化它。请帮忙。谢谢。

4

2 回答 2

6

这是一种非常简单快捷的方法：

>> min([A(1:(end-2)); A(2:(end-1)); A(3:end)], [], 1)

ans =

     2     2     2     3     3     3     5     8

编辑：既然你想要一个完整的功能......

function running_min = running_min(x, k)

xrep = repmat(x, 1, k);
xrep = reshape([xrep zeros(1, k)], length(x)+1, k);
running_min = min(xrep, [], 2)';
running_min = running_min(1:end-k);
于 2011-11-02T01:28:16.697 回答
4

您提到的帖子提供了构建滑动窗口的通用解决方案（您可以控制：重叠与不同，滑动步骤，重叠量，窗口大小）

在您的情况下，它要简单得多，并且可以使用HANKEL函数轻松执行：

x = [13 14 2 14 10 3 5 9 15 8];
idx = hankel(1:3, 3:length(x))
min( x(idx) )

如果您想构建可重用的解决方案：

function y = myFilter(x,k,fcn)
    idx = hankel(1:k, k:length(x));
    y = cellfun(fcn, num2cell(x(idx),1));
end

我们用作：

x = [13 14 2 14 10 3 5 9 15 8];
y = myFilter(x, 3, @(x)min(x))

注意我正在使用 CELLFUN 以防fcn无法以矢量化方式跨维度操作......

于 2011-11-09T00:30:44.623 回答