피크 폭 결정하기
종형 곡선의 합으로 구성된 신호를 생성합니다. 각 곡선의 위치, 높이, 폭을 지정합니다.
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
돌출부의 절반을 기준으로 사용하여 피크의 폭을 측정합니다.
findpeaks(PeakSig,x,'Annotate','extents')
폭을 다시 측정하되, 이번에는 높이의 절반을 기준으로 사용합니다.
findpeaks(PeakSig,x,'Annotate','extents','WidthReference','halfheight') title('Signal Peak Widths')