Transformada de fourier

%Aplicar transformada de fourier discreta
%Veo dos señales
fftsignal1 = fft(M(;,1))
fftsignal2 = fft(M(;,2))

%Pintar
close all
clf
figure

subplot(2,2,1)
plot(abs(fftsignal1))
axis ([0,1,0,800])
axis ("auto x")
title ("Diagrama de magnitud señal 1");
xlabel("Frecuencia (Hz)")
ylabel("Amplitud")

subplot(2,2,2)
plot(angle(fftsignal1))
title ("Diagrama de fase señal 1");
axis ("auto x")
xlabel("Frecuencia (Hz)")
ylabel("Fase (radianes)")

subplot(2,2,3)
plot(abs(fftsignal2))
axis ([0,1,0,800])
title ("Diagrama de magnitud señal 2");
axis ("auto x")
xlabel("Frecuencia (Hz)")
ylabel("Amplitud")

subplot(2,2,4)
plot(angle(fftsignal2))
title ("Diagrama de fase señal 2");
axis ("auto x")
xlabel("Frecuencia (Hz)")
ylabel("Fase (radianes)")