Создайте сигнал, который состоит из суммы кривых нормального распределения. Задайте местоположение, высоту и ширину каждой кривой.
x = linspace(0,1,1000); Pos = [1 2 3 5 7 8]/10; Hgt = [4 4 2 2 2 3]; Wdt = [3 8 4 3 4 6]/100; for n = 1:length(Pos) Gauss(n,:) = Hgt(n)*exp(-((x - Pos(n))/Wdt(n)).^2); end PeakSig = sum(Gauss);
Постройте отдельные кривые и их сумму.
plot(x,Gauss,'--',x,PeakSig)
grid
Измерьте ширины peaks с помощью половины выдающегося положения в качестве ссылки.
findpeaks(PeakSig,x,'Annotate','extents')
Измерьте ширины снова, на этот раз с помощью половины высоты в качестве ссылки.
findpeaks(PeakSig,x,'Annotate','extents','WidthReference','halfheight') title('Signal Peak Widths')