Programación en Matlab: Algoritmo FFT N=2^L basado en extracción de tiempo y extracción de frecuencia.

Calculado en base a la fórmula de transformada discreta de Fourier de secuencia de longitud finita: f1=50;%10Hz

f2=100;?%100Hz

% frecuencia de muestreo

Fs=1000;?%100Hz

% puntos de muestreo N

L=10;

N=2^L;

p>

% de secuencia de pulsos de muestreo

n?=?0:N-1;

t?=?n./Fs;

% ?f2?Número de muestras en un ciclo

M?=?floor(Fs/f2);

% señal muestreada

x ?=?cos( 2*pi*f1.*t)+sin(2*pi*f2.*t);

% secuencia de muestreo

subplot(311);

stem(t(1:2*M),x(1:2*M));

hold?off;

% transformada de Fourier

%Calcular DFT basado en la fórmula de transformada discreta de Fourier de secuencias de longitud finita

n?=?0:N-1;

k?=?0: N-1;

F?=?x?*?exp(-j*2*pi/N).^(n'*k);

subtrama(312) ;

plot(n,abs(F));

subplot(313);

plot(k,angle(F));