load handel
Fmi=220*2^(7/12);
n=[0:1:length(y)-1]';
ruido_mi=sin(2*pi*Fmi*n/Fs);
yruidosa=y+0.3*ruido_mi;
N=2^17;
Y=fft(y,N);
Yruidosa=fft(yruidosa,N);
F=[-Fs/2:Fs/N:Fs/2-Fs/N]';
Fresos=220*2^(6/12);
Ffa=220*2^(8/12);
kresos=floor(Fresos*N/Fs);
kfa=ceil(Ffa*N/Fs);
filtro=ones(N,1);
filtro(kresos:kfa)=0;
filtro(N-kfa:N-kresos)=0;
Yfiltrada=filtro.*Yruidosa;
yfiltrada=real(ifft(Yfiltrada,N));
k=[0:1:N-1]';
figure(1)
plot(k,filtro)
figure(2)
plot(F,abs(fftshift(Y)))
figure(3)
plot(F,abs(fftshift(Yruidosa)))
figure(4)
plot(F,abs(fftshift(Yfiltrada)))
figure(5)
plot(F,abs(fftshift(Yfiltrada)),F,500*fftshift(filtro))
soundsc(y,Fs)
soundsc(yruidosa,Fs)
soundsc(yfiltrada,Fs)