PRACTICA No.
3
“DESARROLLO DE LA SERIE DE FURIER"
I. OBJETIVO
Haciendo uso de OCTAVE u otro software, desarrollar la transformada de funciones no
periódicas y la transformada Rápida de Fourier FFT de señales muestreadas y mostrar las
gráficas correspondientes en el dominio del tiempo y la frecuencia.
II. EQUIPOS Y MATERIALES
Matlab
Octave
PC Pentium II – actual
III. INTRODUCCION TEORICA
Matlab
Matlab es una herramienta interactiva basada en matrices para cálculos científicos y de
ingeniería (de hecho, el termino matlab procede de matrix laboratory). Desde el punto de vista
del control, matlab se puede considerar un entorno matemático de simulación que puede
utilizarse para modelar y analizar sistemas. Permitir a el estudio de sistemas continuos,
discretos, lineales y no lineales, mediante descripción interna y externa, en el dominio temporal
y frecuencia.
Octave
Octave o GNU Octave es un programa y lenguaje de programación para realizar cálculos numéricos.
Como su nombre indica, Octave es parte del proyecto GNU. Es considerado el
equivalente libre de MATLAB. Entre varias características que comparten, se puede destacar que
ambos ofrecen un intérprete, permitiendo ejecutar órdenes en modo interactivo. Nótese que Octave
no es un sistema de álgebra computacional, como lo es Maxima, sino que está orientado al análisis
numérico.
IV. PROCEDIMIENTO
1. Desarrolle la transformada de fourier usando octve u otro software cuya expresión
N=128;
t=linspace(0,3,N);
f=2*exp(-20*t);
figure(1)
plot(t,f)
xlabel('Time,seg'),ylabel('f(t)'),grid
axis([0 0.3 0 2]);
Ts=t(2)-t(1);
Ws=2*pi/Ts;
F=fft(f);
Fp=F(1:N/2+1)*Ts;
W=Ws*(0:N/2)/N;
figure(2)
plot(W,abs(Fp),'+')
xlabel('Frequency,Rad/s'),ylabel('|F(W)|')
Simulación
2. Desarrolle la gráfica de tranformada de fourier desarrollada:
N=128;
t=linspace(0,3,N);
Ts=t(2)-t(1);
Ws=2*pi/Ts;
W=Ws*(0:N/2)/N;
Fa=2./(20+j*W);
figure(3)
plot(W,abs(Fa))
xlabel('Frequency,Rad/s'),ylabel('|F(W)|')
Simulación:
3. Desarrolle la transformada rápida de fourier de una señal muestreada
𝟐𝝅
𝑿(𝒌) = ∑ 𝒙(𝒏)𝒆− 𝟒 𝒏𝒌𝒋 , Donde k = 0,1,2,3,.....
Cuyo desarrollo esta dada por el siguiente programa:
m=[0,1,2,3,4,5];
Xn=[1,2,3,4,5,6];
Xk=fft(Xn);
Xmag=abs(Xk);
Xphase=angle(Xk);
figure(1)
plot(m,Xmag) ,axis([0 5 0 23]);
figure(2)
stem(m,Xmag)
figure(3)
stem(m,Xphase)
Simulación:
4. Para la suma de dos señales senoidales contaminado con ruido desarrolle la gráfica en el
dominio del tiempo y su respectiva transformada de fourier.
t=0:0.001:0.6;
x=sin(2*pi*50*t)+sin(2*pi*120*t);
y=x+2*randn(size(t));
figure(4)
plot(y(1:50))
Y=fft(y,512);
Pyy=Y.*conj(Y)/512;
f=1000*(0:255)/512;
figure(5)
plot(f,Pyy(1:256))
Simulación:
5. Desarrolle la transformada de fourier de la suma de tres señales senoidales:
Fs=100;
t=(1:100)/Fs;
s1=5*sin(2*pi*t*5); s2=10*sin(2*pi*t*15); s3=7*sin(2*pi*t*30);
s=s1+s2+s3;
figure(1)
plot(t,s);
S=fft(s,512);
w=(0:255)/256*(Fs/2);
figure(2)
plot(w,abs([S(1:256)]));
Simulación:
6. Desarrolle la gráfica de la función de muestreo Sa(x):
fplot('5*sin(x)./x',[-30 30 -.2 6])
title('Fplot of f(x)=5.sin(x)/x')
xlabel ('x')
ylabel('f(x)')
Simulación:
V. CUESTIONARIO
1. Desarrolle la transformada rápida de fourier de la función Sa(t).
Función Sa(t)
Transformada:
𝒆𝒋𝒘𝒕 +𝒆−𝒋𝒘𝒕
2. Si 𝒇(𝒕) = 𝟐
= 𝒄𝒐𝒔(𝒘𝒕). Determine su transformada rápida de fourier.
Función: cos(wt)
Transformada:
3. Dado 𝒇(𝒕) = 𝟐𝒔𝒆𝒏(𝒘𝒕). Desarrolle su transformada rápida de Fourier.
Función: 2sin(wt):
Transformada:
4. Desarrolle la transformada de Fourier de la señal muestreada m=[0, 1, 2, 3] y Xm=[2,
3,4, 5].
Unión de los puntos muestreados:
Grafica de la magnitud
Grafica de la fase
VI. CONCLUSIONES
La transformada rápida de Fourier nos sirve para graficar más rápido y computacionalmente sea
mucho más fácil. También sirve para el tratamiento digital de señales y filtrado digital en general a
la resolución de ecuaciones en derivadas parciales o los algoritmos de multiplicación rápida de
grandes enteros.