Full Width at Half Max Peaks

148 views (last 30 days)
Naif Alsalem
Naif Alsalem on 5 Nov 2020
Commented: Image Analyst on 10 Mar 2022
Dear all,
I have 3 peaks close to eachother and I would like to find the FWHM and draw some lines like the graph produced by the Image Analyst here: https://au.mathworks.com/matlabcentral/answers/510144-peak-width-at-half-height
I have loaded my data and normalized already then I utilized this following code (with very slight modification) produced by the Image Analyat here:
% Find the half height - midway between min and max y values.
halfHeight = (min(normy) + max(normy)) / 2;
hold on;
yline(halfHeight, 'Color', 'g', 'LineWidth', 2);
% Find left edge
index1 = find(normy >= halfHeight, 1, 'first');
x1 = x(index1)
line([x1, x1], [0, normy(index1)], 'Color', 'r', 'LineWidth', 2);
text(x1+1, 5, sprintf('x1 = %.2f', x1), 'FontSize', 1, 'Color', 'r');
% Find right edge
index2 = find(normy >= halfHeight, 1, 'last');
x2 = x(index2)
line([x2, x2], [0, normy(index2)], 'Color', 'r', 'LineWidth', 2);
text(x2+1, 5, sprintf('x2 = %.2f', x2), 'FontSize', 1, 'Color', 'r');
% Compute the full width, half max.
fwhm = x2 - x1
text( -0.02688, halfHeight+2, sprintf('width = %.2f', fwhm), 'FontSize', 1, 'Color', 'r', 'HorizontalAlignment', 'center');
caption = sprintf('Full Width, Half Max = %.2f', x2 - x1);
title(caption, 'FontSize', 1);
The "normy" is the y-axis values normalized. The value -0.02688 in the "text" command is the c1 value which is used here in the fwhm equation:
f(x)=a1*exp(-((x-b1)/c1)^2). The produced graph shows the right edge is not accurate and it rather picked up the centre of the peak itself.
Is there anything I am missing here?
Thank you
  2 Comments
Naif Alsalem
Naif Alsalem on 5 Nov 2020
Sindar,
Thank you very much for your answer. It really looks like that it is the last point the code is picking. I was just hoping it picks the half point at that side (right side) just like in the left side.

Sign in to comment.

Answers (1)

Star Strider
Star Strider on 5 Nov 2020
You have not posted the release you are using, however the Signal Processing Toolbox findpeaks function will return the FWHM for every peak it identifies (that you want it to identify).
  2 Comments
Image Analyst
Image Analyst on 10 Mar 2022
@Francesco Digeronimo if it doesn't mention it int he documentation you'd have to get it manually knowing the left and right location and the signal.

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!