Hola amigos, quisiera que me ayuden con lo siguiente:

Necesito graficar en matlab la transformada de fourier de la funcion impulso(delta de dirac) y de la función escalón unitario:

¿alguien sabe el código para poder hacerlo?

He intentado esto pero no me sale:

Fs=6400;
T=1/Fs;
N=Tfinal/T;
n=1:N;

s=ones(1,100);
t=1:length(s);

figure;
plot(t,s);
grid;

transformada=fft(s);
espectro_magnitud=abs(transformada);
espectro_fase=(180/pi)*angle(round(transformada));

resf=Fs/N;
f=-Fs/2:resf:Fs/2-resf;

figure;

subplot(2,1,1);
stem(f,fftshift(espectro_magnitud));grid;

subplot(2,1,2);
stem(f,fftshift(espectro_fase));
grid;