b=11025;
clear a
points = 30;
focus = 880;

for i=1:points
    v = fastconv(w1(1:b),4096*2,focus-points/2+i, 14.7); %subplot(1,1,1);plot(v)
    [av1,fv1] = desa(v);%plot(av1);r=find(av1==max(av1));
    above = time_above( av1 );
    maxA  = max( av1 );
    favg  = sum( fv1( above(1):above(2) ) );
    favg  = favg/ ( above(2) - above(1) );
    favg  = favg*44100/2/pi;
    a(:,i) = av1;
    f(:,i) = favg;
end

plot( sum( a ) )
find(sum(a)==max(sum(a)))