Dear all, how to select the values for k and f, of savitzky-golay filter?

5 次查看（过去 30 天）

显示更早的评论

Ayisha Nayyar 2018-9-15

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/419142-dear-all-how-to-select-the-values-for-k-and-f-of-savitzky-golay-filter

评论： Ayisha Nayyar 2018-9-16

采纳的回答： Image Analyst

Dear all, how to select the values for k and f, of savitzky-golay filter?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

采纳的回答

Image Analyst 2018-9-15

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/419142-dear-all-how-to-select-the-values-for-k-and-f-of-savitzky-golay-filter#answer_336952

You pick them to give the signal you want. There is no one right value for all situations. It's a judgment call. A longer frame width will take more points into consideration when fitting the polynomial and give a smoother curve. A higher polynomial will more accurately "hug" the actual/noisy training data, while a lower polynomial order will give a smoother, less "huggy" fit.

5 个评论
显示 3更早的评论隐藏 3更早的评论

Image Analyst 2018-9-15

在 MATLAB Online 中打开

Try this:

clc;    % Clear the command window.
clearvars;
close all;  % Close all figures (except those of imtool.)
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 20;
% Make a noisy sine wave signal
x = 1 : 3000;
period = 500
y = sin(2*pi*x/period);
noiseAmplitude = 0.8;
y = y + noiseAmplitude * rand(size(y));
subplot(2,1,1);
plot(x, y, 'b-', 'LineWidth', 2);
grid on;
title('Noisy Signal', 'FontSize', fontSize);
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);
% Give a name to the title bar.
set(gcf, 'Name', 'Demo by ImageAnalyst', 'NumberTitle', 'Off')
% Try orders of 1-4, and frame lengths of 3 to about a quarter of your signal length.
frameLengths = 3 : 8 : 201
for polynomialOrder = 2 : 4
  for k = 1 : length(frameLengths)
    % Make sure windowWidth is long enough given the polynomial order we are using.
    windowWidth = max([polynomialOrder + 1, frameLengths(k)]);
    % Make sure windowWidth is odd
    if rem(windowWidth, 2) == 0
      windowWidth = windowWidth + 1;
    end
    caption = sprintf('Smoothed Signal with Polynomial Order = %d, windowWidth = %d', polynomialOrder, windowWidth);
    fprintf('%s\n', caption);
    
    % Now smooth with a Savitzky-Golay sliding polynomial filter
    smoothY = sgolayfilt(y, polynomialOrder, windowWidth);
    subplot(2,1,2);
    plot(x, smoothY, 'b-', 'LineWidth', 2);
    grid on;
    title(caption, 'FontSize', fontSize);
    
    promptMessage = sprintf('Do you want to Continue processing,\nor Quit processing?');
    titleBarCaption = 'Continue?';
    buttonText = questdlg(promptMessage, titleBarCaption, 'Continue', 'Quit', 'Continue');
    if contains(buttonText, 'Quit')
      break;
    end
  end
end

Adapt as needed for your signal.

Ayisha Nayyar 2018-9-16

really thankful....I am trying to implement this discussion on my problem.

请先登录，再进行评论。

更多回答（0 个）

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

类别

Signal Processing Signal Processing Toolbox Signal Generation, Analysis, and Preprocessing Smoothing and Denoising

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

标签

产品

Signal Processing Toolbox

版本

R2015a

Community Treasure Hunt

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

Translated by