1.
%% TRANSFORMADA DISCRETA DE FOURIER
N = 1024; %definir un valor para N
x = 10 * rand(1,N); %definir una señal con muestras discretas
for k=0:N-1
a = x(1)*exp(-2*pi*j*k*0/N)
for n=1:N-1
a = x(n+1)*exp(-2*pi*j*n*k/N)+a
end
TX(k+1)=a;
end
%Comparando este resultado con el del comando o función "FET".
TX1=fft(x,N);
[TX,"TX1"];
2. %% Transformada Discreta de Fourier de una señal
periódica.
% Figura.1 Señal discreta con Fs = 100Hz
Fs = 100;
n=0:Fs-1;
x=cos(2*pi*n*20/Fs);
stem(n,x)
X_16=fft(x,16);
X_512=fft(x,512);
3. %% Figura.2 Transformada Discreta de Fourier.
Módulo y Fase (N=16).
X_16=fft(x,16);
X_512=fft(x,512);
mX_16 = abs (X_16); % módulo para N=16 puntos
fX_16 = angle (X_16); % fase para N=16 puntos
mX_512 = abs (X_512); % módulo para N=512 puntos
fX_512 = angle (X_512); % fase para N=512 puntos
f_16 = linspace(0,Fs,16);
f_512 = linspace(0,Fs,512);
subplot(2,1,1), stem(f_16 , mX_16)
subplot(2,1,2), stem (f_16 , fX_16)
4. %% Figura.3 Transformada Discreta de Fourier.
Módulo y Fase (N=512).
figure (2)
subplot(2,1,1), stem(f_512, mX_512)
subplot(2,1,2), stem (f_512, fX_512)
5. C%% Figura.4 Transformada Discreta de Fourier. Real e Imaginaria
(N=512).
rX_16 = real (X_16);
iX_16 = imag (X_16);
rX_512 = real (X_512);
iX_512 = imag (X_512);
subplot(2,1,1), stem (f_512, rX_512)
subplot(2,1,2), stem (f_512, iX_512)
6. %% Transformada Discreta de Fourier Inversa IFFT.
% Figura .5 Transformada Discreta de Fourier Inversa
ix = ifft (X_512, 512);
rix = real (ix);
stem(0:Fs-1, rix(1:Fs))
7. %% GRÁFICA DE LA FUNCIÓN FFT
% Figura .6 Módulo de la TDF de una señal periódica
Fs = 120;
n = 0:Fs-1;
x = cos(2*pi*n*10/Fs) + 2*cos(2*pi*n*30/Fs) + 4*cos (2*pi*n*43/Fs);
tX = fft(x,1024);
tX = fftshift( tX );
mtX = abs(tX);
f = linspace( -Fs/2,Fs/2 ,1024);
plot( f, mtX), grid
8. %% Graficando la FFT de una Señal Coseno sumado a un ruido
% Figura 7 Módulo de la TDF de una señal periódica con ruido
Fs=100;
n = 0:Fs-1;
x = 2*cos(2*pi*n*10/Fs);
r = randn(1,Fs) / 1.5;
xr = x + r;
figure(1) , plot(xr)
tXR = fft (xr, 1024);
mtXR = abs (tXR);
f = linspace(0,Fs, 1024);
plot(f,mtXR)
9. %% Graficando la FFT con N = 1024 para una Señal Cuadrada de F = 10
Hz.
Fs=200;
n=0:Fs-1;
x=square(2*pi*n*10/Fs);
figure (1);
stem(x)
10. %% Figura 8 Módulo de la TDF de una señal periódica (cuadrada)
tX= fft ( x , 1024 );
mtX= abs (tX);
mtX = fftshift( mtX );
f = linspace( -Fs/2 ,Fs/2 ,1024);
plot(f,mtX)
