# Gaussian fit to xy data and extracting FWHM

286 visualizaciones (últimos 30 días)
Matlab User el 8 de En. de 2020
Respondida: Fynn Reinbacher el 5 de Nov. de 2020
Hello,
I am hoping someone could help. I am going round in circles although i'm sure it's a straightforward answer.
I have an xy data set (see attached for example). I want to fit a gaussian to this, from which the FWHM of the gaussian and the goodness of fit extracted (those that do not have good gaussian fits will be then rejected on the basis of r squared values.
I tried using cftool, and from what i have read, I can use FWHM= 2*sqrt(2*log(2))*c1 to get the FWHM. But the value I get (~30), doesn't match up to what i expect from looking at the curve at the half-maximum (~20). Thankyou
##### 0 comentariosMostrar -1 comentarios más antiguosOcultar -1 comentarios más antiguos

Iniciar sesión para comentar.

Daniel M el 8 de En. de 2020
This is because of the slightly different way cftool has defined the gaussian equation for the fit, and it ends up multipling the c1 coefficient by a factor of sqrt(2) from the true value of the standard deviation.
The equation for FWHM is
FWHM = 2*sqrt(2*log(2))*sigma
%%% sigma, NOT c1!
For fitting, MATLAB uses
f(x) = a1*exp(-((x-b1)/c1)^2)
And wherever you're reading the FWHM equation defined the gaussian as
f(x) = (1/sigma/sqrt(2*pi))*exp(-(x-mu)^2/(2*sigma^2))
By equating the two, you'll find that
sigma = c1/sqrt(2)
Therefore the FWHM equation becomes
FWHM = 2*sqrt(2*log(2)) * (c1/sqrt(2))
FWHM = 2*sqrt(log(2)) * c1
And in fact, when you use this new equation, you get FWHM = 21.78, which is accurate.
##### 1 comentarioMostrar NingunoOcultar Ninguno
Matlab User el 8 de En. de 2020
Thank-you so much!

Iniciar sesión para comentar.

### Más respuestas (1)

Fynn Reinbacher el 5 de Nov. de 2020
##### 0 comentariosMostrar -1 comentarios más antiguosOcultar -1 comentarios más antiguos

Iniciar sesión para comentar.

### Categorías

Más información sobre Point Cloud Processing en Help Center y File Exchange.

### Community Treasure Hunt

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

Start Hunting!