load('TP3_1.mat')
L=length(y);
Fs=1000;  % Frecuencia de muestreo
N=4096*4;   % notar que N > L
Y=fft(y,N);
F=[-Fs/2:Fs/N:Fs/2-Fs/N]';
[picomax,ind]=max(abs(fftshift(Y)));
A=2*picomax/L;  % amplitud de la senoide
Fsen=-F(ind);   % frecuencia de la senoide
figure(1)
plot(t,y)
figure(2)
plot(F,abs(fftshift(Y)))
%============================
% Filtrado en el dominio FFT
%============================
[pico,indi]=max(abs(Y));
filtro=zeros(N,1);
filtro(indi-20:indi+20)=1;
filtro(N-indi-20:N-indi+20)=1;
Yfiltrada=Y.*filtro;
yfiltrada=real(ifft(Yfiltrada,N));
yfiltrada=yfiltrada(1:L);
figure(3)
plot(F,abs(fftshift(Yfiltrada)),F,500*fftshift(filtro))