Peak finding and measurement: TestPeakfind.m - File Exchange

Code covered by the BSD License

Highlights from
Peak finding and measurement

P=ipeak(DataMatrix,PeakD,AmpT... Version 3.91, FitW now equals number of points, not number of intervals;
AmpTSlider(n,h) Changes AmpThreshold when the AmpThreshold slider is moved.
BGSlider(n,h) Called when the BG slider is clicked.
FitSlider(n,h) Changes FitWidth when the FitWidth slider is moved.
P=findpeaks(x,y,SlopeThreshol... function P=findpeaks(x,y,SlopeThreshold,AmpThreshold,smoothwidth,peakgroup)
P=findpeaks(x,y,SlopeThreshol... function P=findpeaks(x,y,SlopeThreshold,AmpThreshold,smoothwidth,peakgroup,smoothtype)
P=findpeakslidersRedraw(x,y,S... Redraws graph for FindPeakSliders when slider are moved.
PeakSlider(n,h) Changes PeakNumber when the Peak slider is moved.
SlopeTSlider(n,h) Changes SlopeThreshold when the SlopeThreshold slider is moved.
SmoothSlider(n,h) Changes SmoothWidth when the SmoothWidth slider is moved.
V=findvalleys(x,y,SlopeThresh... function P=findvalleys(x,y,SlopeThreshold,AmpThreshold,smoothwidth,peakgroup,smoothtype)
[index,closestval]=val2ind(x,... Returns the index and the value of the element of vector x that is closest to val
d=deriv(a) First derivative of vector using 2-point central difference.
d=secderiv(a) Second derivative of vector using 3-point central difference.
exp(-((x-pos)./(0.6006.*wid))...
lorentzian(x,position,width)
rtslid(fig,f,hh,varargin) RTSLID Slider widget that responds to dragging realtime
smoothwidth. Works well with ...
sy=condense(y,n) Condense y by a factor of n, where n is a non-zero positive integer.
Demo4peaks.m
DemoFindPeak.m
DemoFindPeak.m
DemoFindPeakSliders.m
DemoFindPeakSliders2.m
DemoFindPeakSlidersG.m
FindPeakSliders.m
FindPeakSlidersG.m
RedrawPeak.m
TestPeakfind.m
ipeakdemo.m
ipeakdemo1.m
ipeakdemo2.m
ipeakdemo3.m
ipeakdemo4.m
View all files

from Peak finding and measurement by Tom O'Haver
Function to locate and measure the positive peaks and valleys in noisy data sets.

TestPeakfind.m

x=[1:1000000];
load test y
factor=2;
x=condense(x,factor);
y=condense(y,factor);
starttime=cputime;
P=findpeaks(x,y,0.0000000007,.2,50,50);
ElapsedTime=cputime-starttime;
 ['Found '  num2str(length(P)) ' peaks in a ' num2str(length(x)) '-point signal in ' num2str(ElapsedTime) ' seconds.']
PointsPerSecond=length(y)/ElapsedTime
plot(x(1:2000/factor),y(1:2000/factor))

Highlights from Peak finding and measurement

Highlights from
Peak finding and measurement