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>
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 p>
%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));