Matlab 中向量的高斯滤波器

Question

我有一个 n 维向量（1xn 数据集，它不是图像数据），我想对其应用高斯滤波器。我有图像处理工具包和其他一些工具包（询问您是否需要列表）。

据推测，我可以将 fspecial 函数的 hsize 参数设置为类似于 [1 n] 的形式。 作为下一步，我仍然可以使用 imfilter 将其应用到我的向量中，还是应该使用其他东西？

我已经看过很多关于如何在 Matlab 中将高斯滤波器应用于二维图像数据的示例，但我对 Matlab 作为平台还比较陌生，因此一个示例非常好。

注意：我目前无法尝试一下，看看会发生什么（目前不在安装了 Matlab 的机器上），否则我会先尝试一下，然后只询问我是否在使用 fspecial< 时遇到问题/code> 和 imfilter。

Answer 1

为什么不自己创建高斯滤波器？您可以查看 fspecial 中的公式（或高斯的任何其他定义）：

sigma = 5;
sz = 30;    % length of gaussFilter vector
x = linspace(-sz / 2, sz / 2, sz);
gaussFilter = exp(-x .^ 2 / (2 * sigma ^ 2));
gaussFilter = gaussFilter / sum (gaussFilter); % normalize

为了应用它，您可以使用过滤器：

y = rand(500,1);
yfilt = filter (gaussFilter,1, y);

并且不要忘记过滤器有延迟，这意味着滤波后的信号与输入信号相比发生了偏移。由于此过滤器是对称的，因此您可以使用 conv 而不是 filter 获得非移位输出，并使用 same 选项：

yfilt = conv (y, gaussFilter, 'same');