View License

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video

Highlights from

5.0 | 4 ratings Rate this file 10 Downloads (last 30 days) File Size: 20 KB File ID: #5661 Version: 1.3
image thumbnail



Jan Simon (view profile)


09 Aug 2004 (Updated )

Fast Savitzky Golay filter as multi-threaded C-Mex

| Watch this File

File Information

fSGolayFilt is a fast polynomial smoothing filter for uniformly spaced signals. This implementation offers less flexibility than other versions found in the MATLAB Central (eg. from Jianwen Luo, YangQuan Chen or Vassili Pastushenko), but it is optimized for speed by using a C-Mex function for the calculations.

Y = fSGolayFilt(X, K, F, [W], [Dim])
  X: Noisy data, real DOUBLE or SINGLE vector or array.
  K: Polynomial order.
  F: Length of the moving window.
  W: Weighting factors, optional.
  Dim: Dimension to operate on, optional.
  Y: Smoothed data, same size and type as X.

Compared to SGOLAYFILT of the Signal Processing Toolbox of Matlab 2009a, this implementation is 70% to 85% faster for short [100 x 1] and large [1E7 x 1] vectors (see screenshot taken on a single threaded machine). Multi-threading can increase the speed remarkably, but Matlab's filter method benefits from multiple cores also.

The C-mex have to be compiled as usual: mex -O fSGolayCore.c
Unfortunately multi-threading is implemented for Windows only. For OS-X and Linux a single-thread version is included. For futher help see fSGolayCore.c, e.g.: how to set the maximum number of threads at compile time and limit the number of used threads at run time.
After compiling run the unit-test TestfSGolayFilt to test validity and speed on your machine. The test does not need the Signal-Processing-Toolbox anymore.

Mar-2009: This new version is compatible with 32 and 64 bit addressing. Now the initial and final phases are computed in the Mex also to improve the performance by about 10%.
Jun-2010: Optimized for MSVC 2008 (30% faster). Multi-threading increases the speed in addition. Input [Dim]. Computes DOUBLE and SINGLE arrays.

I'd appreciate suggestions for improvements and bug reports sent through email - thanks.


Savitzky Golay Smoothing Filter, Savitzky Golay Smoothing And Differentiation Filter, 2 D Savitzky Golay Smoothing Filter, and Savfilt inspired this file.

This file inspired Z Transform Of 1 D & 2 D Savitzky Golay Smoothing And Differentiation Filter.

MATLAB release MATLAB 7.8 (R2009a)
Other requirements Tested: Matlab 6.5, 7.7, 7.8, Compiler: LCC2.4/3.8, BCC5.5, Open Watcom 1.8, MSVC 2008 (preferred) Pre-compiled Mex:
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (7)
18 Jun 2016 Jan Simon

Jan Simon (view profile)

@Taro Ichimura: As mentioned in the description, the C-mex have to be compiled at first: mex -O fSGolayCore.c

Comment only
16 Jun 2016 Taro Ichimura

look nice but when using
Y = fSGolayFilt(rawdata, 3, 25) on a 1x800 X matrix, I have this error:

Undefined function or variable 'fSGolayCore'.

Error in fSGolayFilt (line 201)
Y = fSGolayCore(X, F, B);

Error in testsmoothing (line 7)
Y = fSGolayFilt(rawdata, 3, 25)

Comment only
25 Jul 2012 AM

AM (view profile)

Works great for smoothing, but have a question. How would I get the actually coefficients out to be able to calculate derivatives?

Comment only
25 Jul 2012 AM

AM (view profile)

06 May 2011 Rhys

Rhys (view profile)

Fantastic! Works ~80x times faster than the script I was using previously. (Matlab 2010a Windows XP 32bit Intel C2D)

01 Dec 2005 Olivier Berné

Fast, Simple, and accurate ! (Works fine on matlab 7.1)

14 Aug 2004 Jianwen Luo

Excellent! The speedup is really dramatic.

09 Jun 2009 1.1

Inserted a notes for compatibility with Matlab 2009a.

05 Mar 2010 1.2

32 and 64 bit addressing, improved performance

09 Jun 2010 1.3

SINGLE support, multi-threading, faster execution, Dim can be specified, unit-test does not need SPT anymore.

Contact us