%% Script Matlab / Octave
%% Uso de la Fast Fourier Transform
%%
%%                        rhp, me4501 

N=512;                        % Potencia de 2
tmax = 1;                     % en (seg)
t  = linspace(0,tmax,N);  
fs = N/tmax;                  % Frecuencia de Muestreo
f = cos(2*pi*10*t); 
F = fft(f);                   % Fast Fourier Transform
P = abs(fftshift(F)); 

frec = linspace(-fs/2, fs/2, N); % Creamos eje de frecuencias

plot(frec, P);                   % Espectro de f(t)