Deconvolution of fluorescence spectra

I want to make a deconvolution of a fluorescence spectra into its Gaussian components. I am not sure how many components there are, and I want to do it in Matlab. I tried "fitgmdist" but I did not manage to do it so i tried to do it manually generating three gaussian distribution and I adjusted the mean and standard deviation to get the image image below. It is just that I don't know how to adjust the parameter to get a correct estimaton, and I want to implement a algorithm that could do it automatically.

Answer 1

[Jyothis Gireesh]

I am assuming that you want to create an algorithm which can automatically determine the optimal number of Gaussian components needed to fit the given data.

Here are a few suggestions which may help in resolving your issue

• To determine the number of Gaussian components needed to fit the data optimally you may make use of the Akaike Information Criteria (AIC) statistic. AIC gives the optimized loglikelihood values. You may try fitting the spectrum with different models by varying the number of Gaussian components in the “fitgmdist()” function. Select the model with the lowest value of AIC.
• In addition to the above, you may also try regularization by setting the “RegularizationValue” parameter in the “fitgmdist()” function.

You may also make use of the following documentation for further clarification.

https://www.mathworks.com/help/stats/fitgmdist.html