function seefelder(inFile,outFile) x = wavread(inFile); a = x(1:11025); lp = fix(length(x)/330750)-1 lastSample = 330750; if(lp>0) for i=1:lp i b = x((lastSample+1):(lastSample+11025)); xcor = xcorr(a,b); [mx ix] = max(xcor); delta = 11025 - ix; sig = x(((lastSample+1)-delta):((lastSample)-delta+330750)); lastSample = ((lastSample)-delta+330750) t = [1:length(sig)]; ct = x(t); aux(t) = sig(t) - ct(t); diff(i,t) = aux(t); lVec(i)=length(aux); clear b; clear xcor; clear sig; clear aux; end rest = rem(length(x),330750) if(rest>11025) b = x((lastSample+1):(lastSample+11025)); xcor = xcorr(a,b); [mx ix] = max(xcor); delta = 11025 - ix; sig = x(((lastSample+1)-delta):length(x)); t = [1:(length(x)-((lastSample+1)-delta))]; ct = x(t); aux(t) = sig(t) - ct(t); diff(lp+1,t) = aux(t); lVec(lp+1)=length(aux); lastSample = 0; for i = 1:lp+1 lDiff = lVec(i); totalSig(lastSample+1:lastSample+lDiff) = diff(i,1:lDiff); lastSample = lastSample + lDiff; end else lastSample = 0; for i = 1:lp lDiff = lVec(i); totalSig(lastSample+1:lastSample+lDiff) = diff(i,1:lDiff); lastSample = lastSample + lDiff; end end else rest = rem(length(x)/330750) if(rest>11025) b = x(330751:341775); xcor = xcorr(a,b); [mx ix] = max(xcor); delta = 11025 - ix; sig = x((330751-delta):length(x)); t = [1:(length(x)-(330751-delta))]; ct = x(t); totalSig(t) = sig(t) - ct(t); end end wavwrite(totalSig,11025,16,outFile); clear all; %pack;