análisis de fourier -...
TRANSCRIPT
![Page 1: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/1.jpg)
8.- Análisis de
FourierDr. Servando López Aguayo
Agosto-Diciembre 2017
![Page 2: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/2.jpg)
En la clase pasada…
Recordemos:
![Page 3: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/3.jpg)
En esta clase:
Análisis de Fourier
La existencia de dos mundos!
![Page 4: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/4.jpg)
Ah y recordemos hace tiempo… a los super campeones! ☺
![Page 5: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/5.jpg)
Análisis de Fourier
La parte matemática ya la saben… en teoría.
Y la parte física también…. en teoría.
En esta sesión, nos concentraremos en ver
algunos fenómenos que ocurren en la versión
“computacional” del análisis de Fourier:
Series de Fourier
Transformada de Fourier
Transformada rápida de Fourier
![Page 6: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/6.jpg)
Nuestro amigo Fourier Descomposición utilizando sumas de
ondas senoidales.
![Page 7: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/7.jpg)
Recordemos: series de Fourier
¿Cómo obtengo los coeficientes?
¿Qué pasa si la función es par o impar?
![Page 8: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/8.jpg)
Actividad 1 Calcular los coeficientes correspondientes de
la serie de Fourier de la función dada por:
Grafica la aproximación dada por la suma de Fourier usando 5, 10, 50 y 100 Términos.
Comentar y reportar los resultados obtenidos.
![Page 9: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/9.jpg)
La transformada de Fourier
Recordemos el legendario par:
Si pensamos en su versión discreta… ¿qué
podemos concluir en relación
a las series de Fourier?
![Page 10: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/10.jpg)
Extra: percepción visual de
las frecuencias
![Page 11: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/11.jpg)
¿Vemos realmente el mundo
cómo es?
![Page 12: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/12.jpg)
Principio de
incertidumbre
![Page 13: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/13.jpg)
1.- Nuestras simulaciones son hechas en
un espacio discreto.
2.- Necesitamos que nuestras funciones
sean “limitadas en banda”.
Recordar:
![Page 14: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/14.jpg)
Crear el vector x=-1 : .02 : 1;
Crear las funciones:
y1=cos(2*pi*f1*x);
y2=cos(2*pi*f2*x);
A) Graficar y1 & y2 con f1=1 y f2=52.
B) Graficar y1 & y2 con f1 =1 y f2=51
Contestar: ¿Qué se observa? ¿Por qué se da este fenómeno?
Actividad # 2
![Page 15: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/15.jpg)
Errores al realizar el muestreo discreto de
funciones continuas.
En Fourier, se estila que el dominio en
frecuencia angular (kx) sea:
Aliasing
Matlab:
y plot(cos(1:1000),’.’)
![Page 16: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/16.jpg)
Transformada discreta de
Fourier
Consideraciones:
La razón de muestreo:
![Page 17: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/17.jpg)
Transformada discreta de
Fourier
Por lo tanto:
Por lo que se sólo se puede representar
un número FINITO de
frecuencias:
![Page 18: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/18.jpg)
La DFT
Por lo que tenemos como la “DFT”:
![Page 19: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/19.jpg)
La IDFT De manera similar con la transformada
inversa se tiene:
![Page 20: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/20.jpg)
Sin embargo…
Hay varias consecuencias!
El dominio en frecuencias se vuelve
periódico.
Existe un compromiso entre el dominio
temporal y el dominio frecuencial.
Aliasing.
![Page 21: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/21.jpg)
Manera alternativa de la DFT
Podemos reescribir la DFT como:
![Page 22: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/22.jpg)
Vámonos al break!!
Regresamos en 10 minutos!!
![Page 23: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/23.jpg)
Actividad 3 Programar la DFT y la IDFT.
Utilizando su programa, calcular la transformada de Fourier de un pulso Gaussiano que está dado por:
F(t) = exp(-t2)
Comprobar que su transformada de Fourier es otro pulso Gaussiano. Recuperar nuevamente F(t) utilizando la IDFT. Reportar y comentar resultados obtenidos.
![Page 24: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/24.jpg)
Y con ustedes… la fft!
Fast Fourier Transform: es una manera de
calcular la DFT maximizando la velocidad
de cómputo.
Realizado por Cooley y Tuckey en 1965.
![Page 25: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/25.jpg)
Idea de la FFT Utilizar la periodicidad del algoritmo, para
N=8 por ejemplo:
![Page 26: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/26.jpg)
Idea de la FFT Realmente necesitamos calcular tantos
coeficientes?
![Page 27: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/27.jpg)
Idea de la FFT
Por lo que tenemos:
![Page 28: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/28.jpg)
Idea de la FFT
Y reacomodando términos:
![Page 29: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/29.jpg)
Idea de la fft
Y utilizando la “operación mariposa”:
![Page 30: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/30.jpg)
Idea de la FFT
![Page 31: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/31.jpg)
Idea de la FFT Fin: ordenar previamente nuestro vector!
![Page 32: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/32.jpg)
FFT
Eficiencia de la DFT:
Eficiencia de la FFT:
Vale la pena programar dicho algoritmo?
![Page 33: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/33.jpg)
Actividad 4
Comparar el número de operaciones
usando la DFT vs la FFT para N datos,
donde N va de 2 a 4096 datos.
Reportar y comentar los resultados
obtenidos.
![Page 34: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/34.jpg)
Y con Matlab? Tenemos los siguientes comandos: fft, ifft,
pero además fftshift y ifftshift.
Obtener la fft de los datos utilizados en la actividad pasada.
![Page 35: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/35.jpg)
Actividad 4
Calcular la transformada de Fourier de un
pulso Gaussiano utilizando los comandos
de Matlab
Graficar, comentar y reportar los
resultados.
![Page 36: Análisis de Fourier - homepages.mty.itesm.mxhomepages.mty.itesm.mx/slopez/lec_08_fc1_ad_2017.pdf · Análisis de Fourier ... Y con Matlab? Tenemos los ... Actividad 4 Calcular la](https://reader031.vdocumento.com/reader031/viewer/2022012322/5baa350509d3f2c9618b5756/html5/thumbnails/36.jpg)
Y eso es todo!!
Listo mis estimados!!
Nos vemos el siguiente miércoles!!