How to search for the optimal parameters of a filter

3 次查看（过去 30 天）
Z Liang
Z Liang 2019-7-26
评论： Star Strider 2019-8-2
ep_mean is a 26xN matrix of a single trial, recorded through 26 electrodes, over N time points (roughly 220 points sampled)
I applied the sgolayfilt function to try to smooth the signal. I run into the problem that I don't know what values of the two parameters for sgolayfilt, the order and the frame length, are optimal for the filtering process. I also don't know if this is a technical or a philosophical question.
So I tried a grid sweep of parameter combinations using the following code:
snr_mat = [];
for i1 = 1:99
if mod(i1,2) == 0
c = i1+1;
else
c = i1+2;
end
for i2 = c:2:101
rd = i1;
fl = i2;
filt = sgolayfilt(ep_mean,rd,fl,[],2);
snr_mat(i1,i2) = snr(ep_mean,ep_mean-filt);
end
end
snr_mat(snr_mat>300) = 0; % to remove what I thought were artifacts, but now I think it's my setup/interpretation of the function that's flawed
subplot(4,1,1)
plot(ep_mean(:,ep_window_start:ep_window_end)')
title('unfiltered','fontsize',16)
subplot(4,1,2)
y = sgolayfilt(ep_mean,3,21,[],2);
plot(y(:,ep_window_start:ep_window_end)')
ylim([-10 10])
title('rd = 3; fl = 21','fontsize',16)
subplot(4,1,3)
x = sgolayfilt(ep_mean,11,41,[],2);
plot(x(:,ep_window_start:ep_window_end)')
ylim([-10 10])
title('rd = 11; fl = 41','fontsize',16)
subplot(4,1,4)
[M1,I1] = max(snr_mat);
[M2,I2] = max(M1);
rd = I1(I2);
fl = I2;
z = sgolayfilt(ep_mean,rd,fl,[],2);
plot(z(:,ep_window_start:ep_window_end)')
title(['rd = ' num2str(rd) '; fl = ' num2str(fl)],'fontsize',16)
To produce these figures:
5.PNG
6.PNG
The first figure plots the magnitude of the SNR after using each pair of order (rd) and framelength (fl) parameters. It seems the higher both are, and the closer they are to each other, the better then SNR. But the second figure clarifies what I think to be happening, which is that higher values for rd and fl do not necessarily filter better, they just more closely match the unfiltered signal so the snr function calculates a lower noise and higher ratio. This feels quite stupid to me, what I did I mean.
So my question becomes how to identify the optimal parameters (rd, fl) in the sgolayfilt function to filter a given signal. What metric can I employ to assess the performance of the filtering process? And is a grid sweep the best way to test different parameter combinations?
Thanks for your help!

请先登录，再回答此问题。

采纳的回答

Star Strider
Star Strider 2019-7-26
The Savitzky-Golay filter is very good for some applications and not appropriate for others. I use it as a relatively straightforward way to reduce broad-band noise, since frequency-selective filters are completely ineffective in that situation. (The best way to see if broad-band noise is a problem is to do a Fourier transform of your signal.)
If you simply want to reduce high-frequency noise, then a frequency-selective filter is optimal, and there are several ways to design and implement them in MATLAB.
With respect to sgolayfilt, and the problem you pose, I always use a low-order filter (usually 3) and adjust the frame length until I am happy with the result. I only use it where a frequency-selective filter would not provide the result I want.
  2 个评论
Z Liang
Z Liang 2019-8-2
Thank you, I found a paper here that tries to answer half the question. But for now I'm content with whatever my supervisor is happy with!
Star Strider
Star Strider 2019-8-2
As always, my pleasure!
Thank you for the paper: ‘Optimum window length of Savitzky-Golay filters with arbitrary order’. I saved it to my library directory.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Smoothing and Denoising 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by