Coefficient Matrix for Savitzky-Golay Data Filtering

Starting from a polynomial coefficient vector such as: [-3/35 12/35 17/35 12/35 -3/35] I would like to create a coefficient matrix of the form:

by iteration of this vector in order to perfrom the Savitzky-Golay signal smoothing. How is this possible without using a for ... loop?

[Jan]

Why do you want to avoid a loop? Are you aware that this must be done by a loop, even if it is hidden inside a Matlab command? Why do you want to create a potentially huge matrix with such a low density of information? Note that you have 4 independent values only in your example. Do you want at least a sparse array?

The explicite creation of the matrix is such a waste of time, that a FOR loop will not matter at all. If you want to compute the matrix product also, you loose even more time, because the matrix consists almost only of zeros. Use filter instead.

Answer 1

[Image Analyst]

Why do you want this? Who cares what the coefficients are? Just use sgolayfilt() and be done with it.

[Saeid]

Thanks for the answer! I didn't have the Matlab module for this fitting so I wanted to program from scratch. I still haven't used the program you have sent but that could solve the problem and save me buying the Matlab module!

[Image Analyst]

You can do the Savitzky-Golay filter yourself. It's essentially just a polynomial fit to a window that slides across the signal. So just have a for loop to have a window move across and use polyfit() to get the polynomial, then use polyval() to get the value at the center of the window.