clear all
s1 = wavread('for_notes.wav');
[s2, fs] = wavread('taaa.wav');

f2 = fft(s2(1:2^15));
a2 = abs(f2);

f3 = fft(a2);
f3(50:end-50) = 0;
r3 = real(ifft(f3));

f1 = fft(s1(1:2^15))';
a1 = abs(f1);

f4 = fft(a1);
f4(50:end-50) = 0;
r4 = real(ifft(f4));

r = real(ifft(f1 .* r3./r4));
soundsc(r, fs)