function spec = butt_filt(wave,cutoff,n)
f=linspace(0,wave.fs/2,size(wave.spec,1));
lpf = 1 - 1 ./ (1.0 + (f ./ cutoff).^(2*n));
spec=repmat(lpf',1,size(wave.spec,2)).*wave.spec;