function cleannoise()
 sound([0]);
 one = wavread('one.wav');
 figure;
 subplot(2,2,1); plot(one);
 axis([0 3000 -0.5 0.5]);

[x,y,but]=ginput(1);
while (but==1)
   sound(one);
   [x,y,but]=ginput(1);
end
 
 
 
 ff = fft(one);
 subplot(2,2,3); plot(fftshift(abs(ff))); 
 
 for t=300:800
    rr=rand;
  ff(1501+t)=60*rr;
  ff(1501-t)=60*rr;
 end    
 
 withnoise = real(ifft(ff));
 subplot(2,2,2);
 plot(withnoise);
 axis([0 3000 -0.5 0.5]);
 
 
[x,y,but]=ginput(1);
while (but==1)
   sound(withnoise); 
   [x,y,but]=ginput(1);
end


 noiseff = fft(withnoise);
 subplot(2,2,4); plot(fftshift(abs(noiseff)));

[x,y,but]=ginput(1);
 
 subplot(2,2,4); hold on; line([800 800],[-20 40]);
		          line([2200 2200],[-20 40]);

  newff = bandpassfilter(noiseff,1,650);
 
 subplot(2,2,4); hold on; plot(fftshift(abs(newff)),'g');

[x,y,but]=ginput(1);

 cleanednoise = real(ifft(newff));

 subplot(2,2,2); hold on; plot(cleanednoise,'g');
 
[x,y,but]=ginput(1);
while (but==1)
   sound(cleanednoise);
   [x,y,but]=ginput(1);
end

  