matlab - 如何实现匹配过滤器

Question

我有一个模板。我通过对模板的傅里叶变换的共轭进行傅里叶逆变换来计算匹配滤波器的脉冲响应。我想使用 Matlab 中的“过滤器”命令对我可用的 EEG channel 之一执行匹配过滤操作。使用滤波器命令，系数“b”是我的脉冲响应？此外，我想实现 Matlab 代码来对匹配滤波器的输出设置阈值以检测峰值。我该如何实现？

Answer 1

这是你的开始，

% A template is given
temp = randn(100,1);

% Create a matched filter based on the template
b = flipud(temp(:));

% For testing the matched filter, create a random signal which
% contains a match for the template at some time index
x = [randn(200,1); temp(:); randn(300,1)];
n = 1:length(x);

% Process the signal with the matched filter
y = filter(b,1,x);

% Set a detection threshold (exmaple used is 90% of template)
thresh = 0.9

% Compute normalizing factor
u = temp.'*temp;

% Find matches
matches = n(y>thresh*u);

% Plot the results
plot(n,y,'b', n(matches), y(matches), 'ro');

% Print the results to the console
display(matches);

正如安德烈亚斯在他的回答中提到的，不需要傅立叶变换。如果你有一个时域模板，那么它的匹配过滤器就是它自身的时间反转版本(我用 flipud 实现)。随着您的进行，您会发现许多细微差别需要解决。这段代码效果很好，因为我从头到尾都在掌控之中，但是一旦你开始处理真实数据，事情就会变得复杂得多。例如，选择合适的阈值需要对您将要使用的数据有一些了解。

事实上，峰值检测可能是一项非常重要的任务，具体取决于信号的性质等。在我的例子中，峰值检测很容易，因为我的信号与模板完全不相关，除了在中间，我也确切地知道我期望看到的振幅。所有这些假设都是对我用来演示概念的问题的过度简化。