[cham,fs] = wavread('chamclip.wav');

%zero pad the sig
cz = zeros(2^17,1);
cz(1:length(cham))=cham;

%take DFT
cz_dft = fft(cz);

freq = 2^13;
%freq = 2^12+2000;

filter = ones(2^17,1);
bassboost = 0*ones(freq,1); % set coeff to zero to kill the low freqs, higher number to boost bass

bb_filt = filter;
%bb_filt = zeros(2^17,1);  % kill the high freqs
bb_filt(1:freq) = bassboost;
bb_filt(2^17-freq+1:2^17) = bassboost;

filt_cham = real(ifft(cz_dft.*bb_filt));