SMOOTHING FUNCTION FOR FOURIER AMPLITUDE SPECTRUM

SMOOTHING FUNCTION FOR FOURIER AMPLITUDE SPECTRUM

smoothSpectra offers various different window smoothing options for Fourier amplitude spectrum (FAS) including boxcar, triangle, Parzen, Hann, Hanning, Hamming, Gaussian. Default window function is Konno-Ohmachi (see Konno and Ohmachi (1998), page 234), which symmetric in log space. The smoothing is performed by convolving the window function with the FAS. Three examples are provided in demo.m file.

USAGE:

[smoothFAS] = smoothSpectra(Y,varargin)

STATIC INPUT:

Y = vector of Fourier amplitudes (1xn or nx1)

VALID PROP_NAME / PROP_VAL PAIRS:
-----------------------------------------
'w' --> (1x1)-[numeric]-[default:40]
'b' --> (1x1)-[numeric]-[default: 20]
'method' --> [text]-[default: konno-ohmachi]
'debug' --> [text]-[default: False]

NOTES:

w = width of window function (e.g., 100)

b = bandwidth coefficient of konno-ohmachi window (e.g., 20)

method = window function (e.g., boxcar, gaussian, hamming, hann,
hanning, konno-ohmachi, parzen, triang)

debug = 'True' for print debug messages

OUTPUT:

smoothFAS = smoothed vector of Fourier amplitudes (1xn)

EXAMPLES:

see demo.m file

REQUIREMENTS:

smoothSpectra function does not require any MatLAB toolbox

ACKNOWLEDGEMENT:

In preparing this function, I benefitted from konv.m function written by Ali Jadoon.

If you find this code useful for your application, please don't forget to rate it.
For questions / suggestions / comments and bug reports: kalkan76@gmail.com