Hankel transform
Efficient implementations of the Hankel transform and the inverse Hankel transform, respectively.
Author: Marcel Leutenegger

Hi Xuesong. The "Undefined function" is because Matlab will not add folders titled "private" to the search path. Perhaps the name of that folder should be changed.

22 Jun 2013

Peak Fitter
Command-line peak fitter for time-series signals. Version 5.5, July, 2014
Author: Tom O'Haver

A couple of errors in the post below:

sqrt(2 log(2)) = 1.177

sqrt(log(2)) = 0.832

22 Jun 2013

Peak Fitter
Command-line peak fitter for time-series signals. Version 5.5, July, 2014
Author: Tom O'Haver

Hi Tom-

This is a great tool.

I have a question about the constant '0.600516' that appears in the function 'gaussian'. In the comment for that function, you write "half-width=wid". If you define your gaussian as:

exp(-x^2/(2*sigma)^2)

then the half width = sqrt(2log(2))*sigma = 0.776*sigma. If you define your gaussian as:

exp(-x^2/sigma^2)

then the half width = sqrt(log(2))*sigma = 0.549*sigma.

It appears that you chose the latter definition but with "wid=full width at half max". That way wid = 2*sqrt(log(2))*sigma or alternatively sigma = 0.600516*wid.