[x,fs,nbits] = wavread('cntry.wav');    %use any audio file you like
%sound(x,fs)                             %listen to original audio file

h = zeros(1,4000);
h(1) = 1;
h(1000) = 1/4;
h(2000) = 1/16;
h(3000) = 1/4;
h(4000) = 1/16;

y = conv(x,h);
y = y/max(y);                             %normalize filter output
sound(y,Fs);
