% Interp_Tempi_DFT
clear
close all
N=20
x=gausswin(N,5);
figure(1)
stem((0:N-1),x,'filled')
title('x: segnale nei tempi')
pause
X=fft(x);
figure(2)
stem((0:N-1),abs(X),"filled")
title('X: DFT del segnale')
pause
Y=zeros(1,N*4);
if round(N/2)==N/2
Y(1:(N/2))=X(1:(N/2));
Y(N*4-(N/2)+1:N*4)=X((N/2)+1:N);
else
Y(1:floor(N/2))=X(1:floor(N/2));
Y(N*4-floor(N/2):N*4)=X(floor(N/2)+1:N); 
end
figure(3)
Y=Y*4;
stem((0:4*N-1),abs(Y),'r','filled')
title("Y:DFT del segnale moltiplicata per 4 e con 3N zero padding")
pause
y=ifft(Y);
figure(4)
stem((0:4*N-1),real(y),"filled",'r')
title("segnale interpolato di un fattore 3")
