Revision fec1958a3f9502fd8e3054cb0b95f55fc77ba527 authored by Max Turner on 12 April 2016, 00:46:24 UTC, committed by Max Turner on 12 April 2016, 00:46:24 UTC
1 parent ce75ca2
getRebounds.m
function r = getRebounds(peaks_ind,trace,searchInterval)
%get rebound as fraction of peak amplitude
%trace = abs(trace);
peaks = trace(peaks_ind);
r = zeros(size(peaks));
for i=1:length(peaks)
endPoint = min(peaks_ind(i)+searchInterval,length(trace));
nextMin = getPeaks(trace(peaks_ind(i):endPoint),-1);
if isempty(nextMin), nextMin = peaks(i);
else nextMin = nextMin(1); end
nextMax = getPeaks(trace(peaks_ind(i):endPoint),1);
if isempty(nextMax), nextMax = 0;
else nextMax = nextMax(1); end
if nextMin<peaks(i) %not the real spike min
r(i) = 0;
else
r(i) = nextMax;
end
end
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...