find peak

findpeaksDzone.m

function [resultPeak,resultIndex,strmzeroWoPeak] = findpeaksDzone(InDat)

Ndat = length(InDat);
[sizeC,sizeR] = size(InDat);

if Ndat~=(sizeC*sizeR)
	errMess = sprintf('stream not single vector [%d,%d] ;_;',[sizeC,sizeR]);
	error(errMess)
end

Dat = reshape(InDat,Ndat,1);

DiffFoward = [0; Dat(2:(end)) - Dat(1:(end-1))];
%DiffBackward = [flipud( InDat((end-1):-1:1) - InDat((end):-1:2) );0];
DiffBackward = [ Dat(1:(end-1)) - Dat(2:(end));0];

peaksFlag = (DiffFoward>0).*(DiffBackward>0);

[resultIndex]=find(peaksFlag);
resultPeak=Dat(resultIndex);

strmZeroOrPeak =reshape(peaksFlag.*Dat,sizeC,sizeR);