How to do multiple gaussian profiles summation, peak detection and fwhm analysis?

Question

Tusharkumar Sorathiya 2018-9-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/419328-how-to-do-multiple-gaussian-profiles-summation-peak-detection-and-fwhm-analysis

评论： Akira Agata 2018-10-11

Gaussian summation and analysis.txt

I have 10 gaussians profiles f1...f10(see the text file)and "f" is sum of all. QUE_1)I wanted to find fwhm and plot fwhm lines for each 10 curves. QUE_2) I have used findpeaks function to detect the peaks of all curves,but problem i am facing now is ,it can detect the peak but peak does not move when i change the mean values of one or more curves. QUE_3) How to calculate the fwhm and plot it, when I see two or more curves as a final output (f) QUE_4) how I can change the x and y axis values to some wavelength(x) and reflectivity(y) values. I tried to change X= -1550:0.1:-1545(to see x axis as a nano meter values)but it changes the curves,is there any way to relate x axis as wavelength values without affecting bell curves? QUE_5) How to calculate the peak shift in "f"(in x and y direction) when mean values are changed?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Akira Agata 2018-9-19

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/419328-how-to-do-multiple-gaussian-profiles-summation-peak-detection-and-fwhm-analysis#answer_337391

在 MATLAB Online 中打开

Regarding QUE_1 and QUE_3, please try the following:

findpeaks(f1,'WidthReference','halfheight','Annotate','extents')

Without setting output variables, findpeaks function generates graph with peaks and some additional information if you set 'Annotate' option to 'extents'.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Tusharkumar Sorathiya 2018-10-9

Hello Akira, again thanks a lot. Any idea how i can do the same(which you answered) in OCTAVE?

Akira Agata 2018-10-11

Hi Tusharkumar-san,

I regret to say that I'm not so familiar with OCTAVE. Looking at the following help page, it seems that OCTAVE's findpeaks function does not have exactly the same option used in my above post. But, using "parabol", "height" and "baseline" values stored in the 3rd output argument, maybe you can calculate approximate FWHM value for each pulse.

https://octave.sourceforge.io/signal/function/findpeaks.html

请先登录，再进行评论。