В настоящее время я ищу, чтобы найти FWHM сигнала. Пик сигнала составляет около 1,0, но самое низкое значение составляет всего около 0,6. Итак, на самом деле у меня нет половины максимального значения.Попытка найти полную ширину с половиной максимума для шумового сигнала
Как я мог продолжить анализ кривой аналогичным образом?
Вот изображение кривой:
Предполагая, что минимальные и максимальный сигнал y составляет 0,6 и 1 соответственно вы можете найти полувысот:
idx1 и idx2 каждого возвращение индексы из двух точек до и после нужной точки в половине. Мы можем использовать эти точки для интерполяции значения x в половине y.
%height at half
h=(0.6+1)/2;
idx1=find(y>h,1) +[-1 0];
idx2=find(y>h,1,'last') +[0 1];
x1 = interp1(y(idx1),x(idx1),h);
x2 = interp1(y(idx2),x(idx2),h);
w = x2 - x1;
Я предполагаю, что вы должны были бы компенсировать сигнал, центрированный в точке 0 до расчета FWHM. Однако кажется, что FWHM - это значение, которое обычно вычисляется для gaussians, и ваши данные очень не гауссовы –
, если я вас правильно понимаю, вы принимаете смещение по оси x? кривая уже центрирована на x = 0, это видно на изображении моего сигнала, который я прикреплл к своему сообщению. или вы имели в виду что-то еще? –
Я имею в виду ось y, ось x - это то, что вы ищете. Кроме того, он выглядит так, как обычно, чтобы получить FWHM по негауссовским данным, извинения. –