UNIVERSIDAD AUTONOMA DE CHIHUAHUA

FACULTAD DE INGENIERIA INGENIERIA CIVIL

ECUACIONES DIFERENCIALES PARCIALES

"SERIES DE FOURIER, TRANSFORMADA DE FOURIER Y APLICACIONES"

ALUMNOS: MARIA MAGDALENA ISLAS PAREDES JOSE HOMAR SANZABAS RIVERA LUIS ANDRES OSUNA TIRADO JOSE DE JESUS GARCIA LOPEZ 282126 256832 259543 256986

FECHA: 28 DE NOVIEMBRE DE 2013

OBJETIVO
Demostrar la importancia de la Transformada de Fourier y aplicaciones en la ingeniera y otras ciencias.

INTRODUCCION

La idea bsica de las series de Fourier es que toda funcin peridica de perodo T puede ser expresada como una suma trigonomtrica de senos y cosenos del mismo perodo T. El problema aparece naturalmente en astronoma, de hecho Neugebauer (1952) descubri que los Babilonios utilizaron una forma primitiva de las series de Fourier en la prediccin de ciertos eventos celestiales. La historia moderna de las series de Fourier comenz con DAlembert (1747) y su tratado de las oscilaciones de las cuerdas del violn. El desplazamiento u = u(t, x) de una cuerda de violn, como una funcin del tiempo t y de la posicin x, es solucin de la ecuacin diferencial

sujeto a las condiciones iniciales u(t,0) = u (t,1) = 0 para t 0, ( ) para 0 < x < 1. La solucin de este problema es la

superposicin de dos ondas viajando en direcciones opuestas a la velocidad 1, como lo expresa la formula de DAlembert: ( ) ( ) ( )

en la cual f es una funcin impar de perodo 2 que se anula en los puntos x = 0,1,2,... Euler en 1748 propuso que tal solucin poda ser expresada en una serie de la forma ( ) Y como consecuencia ( ) ,

( )

Las mismas ideas fueron luego expuestas por D. Bernoulli (1753) y LaGrange (1759). La formula ( ) para calcular los ( ) coecientes apareci por primera vez en un artculo escrito por Euler en 1777.

SERIES DE FOURIER Se llama serie de Fourier de una funcion f(x) en el intervalo [,] a: f(x) = a0 2 + X n=1 (an cosnx + bn sennx) () A los coecientes a0,a1,,an,b0,b1,,bn se les llama coecientes de Fourier de f(x) en [,]. Debido a que Z senmxsennxdx = ( 0 si n 6= m 6= 0 si n = m Z cosnxdx = 0 Z sennxdx = 0 1 2 Tema 7. Series de Fourier Z

cosmxcosnxdx = ( 0 si n 6= m 6= 0 si n = m Z senmxcosnxdx = 0 e integrando termino a termino en la igualdad () obtenemos: Z f(x)cosnxdx = an Z cos2 xdx = an an = 1 Z f(x)cosnxdx Z f(x)dx = a0 2 2 a0 = 1 Z f(x)dx

Z f(x)sennxdx = bn Z sen2 xdx = bn bn = 1 Z f(x)sennxdx Las anteriores propiedades de las funciones sennx,cosmx se pueden resumir en que el sistema {1,senx,sen2x,,cosx,cos2x} es un sistema ortogonal de funciones respecto del producto escalar (f(x),g(x)) = Z f(x)g(x)dx y la serie de Fourier no es mas que la expresion de un vector f(x) como combinacion lineal de los vectores de la anterior base ortogonal. Denicion 7.2 Se llama serie de Fourier de una funcion f(x) en el intervalo [L,L] a: f(x) a0 2 + X n=1 an cos n L x + bn sen n L x donde a0 = 1L ZL L f(x)dx an = 1L ZL L f(x)cos n L

xdx bn = 1L ZL L f(x)sen n L xdx Este hecho se basa en que el sistema de vectores 1,sen x L ,sen 2x L ,,cos x L ,cos 2x L , es un sistema ortogonal de funciones respecto del producto escalar (f(x),g(x)) = ZL L f(x)g(x)dx Analogamente se puede denir la serie de Fourier de una funcion f(x) denida en un intervalo [a,b] haciendo una traslacion del punto medio a + b 2 al origen. 7.1. Series de Fourier 3 Tomo L = b a+b2 =

ba2 L=a a+b2 = ab2 Denicion 7.3 Se llama serie de Fourier de una funcion f(x) en el intervalo [a,b] a f(x) = a0 2 + X n=1 an cos n ba 2 x + bn sen n ba 2 x ! donde a0 = 2ba Zb a f(x)dx an = 2ba

Zb a f(x)cos 2n b a xdx bn = 2ba Zb a f(x)sen 2n b a xdx Las series anteriores tambien se podran haber escrito de la forma: f(x) C0 + X n=1 Cn cos(n0t n) donde Cn = q a2 n + b2 n, cosn = an q a2 n + b2 n senn = bn q

a2 n + b2 n n = arctang bn an siendo 0 = 1, L , 2 b a segun hayamos utilizado una de las tres formulas anteriores.

TRASFORMADA DE FOURIER

La transformada de Fourier es una operacin que se realiza sobre funciones. Es decir, vamos a coger dos variables, la una dependiente de la otra -una funcin- y la vamos a convertir en otra variable que depende de una nueva -otra funcin-. A la primera funcin la llamamos transformada de Fourier de es: , o si preferimos, . Por definicin, la

Qu nos dice esta expresin? Es una especie de suma infinita -una integral- de la funcin de partida por el nmero elevado a la unidad imaginaria que multiplica a , a la nueva variable que decamos, y a

. Ya tenemos otra funcin, que depende de . Si aplicamos una conocida equivalencia que debemos a Euler , se puede transformar lo anterior en:

Lo que hemos hecho es mirar de otra manera la exponencial elevada a un complejo. Ahora tenemos trigonometra, y es ms fcil entender lo que se le hace a la funcin. Lo que Fourier demostr en realidad es que cualquier funcin que cumpla una serie de condiciones razonables (las condiciones de Dirichlet) es equivalente a su transformada. Esta condiciones razonables se cumplen para la mayora de las funciones con las que vamos a topar, por lo que se puede decir que la TF se puede aplicar casi siempre en el anlisis de series temporales en Astronoma. Y qu tiene de bueno que una funcin sea equivalente a su TF? Lo veremos con unos ejemplos. Supongamos que tenemos la funcin , es decir, un polinomio sencillo, definido entre 0 y 1 por simplicidad. Si aplicamos la TF y dibujamos

el resultado en el mismo plano en que dibujamos 1.

, nos queda el grfico

Figura 1.

Figura 2.

Hemos transformado una funcin en una suma de senos y cosenos que encaja bastante bien con ella. Lo interesante viene cuando aplicamos la transformada a funciones ms complejas. Para una funcin como en la grfica 2, que ya no es tan trivial, tenemos que se sigue realizando un buen ajuste, transformando una expresin polinmica en una simple suma de senos y cosenos.

Y qu pasa si no conocemos la expresin analtica de la funcin de partida?. Llegamos al quid de la cuestin. No siempre vamos a tener polinomios, o a conocer como vara en funcin de . En las series temporales lo que tenemos es una tabla de valores, un nmero, a veces muy grande, de pares de nmeros que se pueden reflejar en una grfica. Pero no sabemos qu sucede si tomamos un valor para la que no est en la tabla. Y precisamente, eso es lo que nos interesa: saber qu pasar entre los valores que tenemos -interpolar- o, ms interesante, ms all de ellos extrapolar-.

APLICACIONES DE LA TRANSFORMADA DE FOURIER En relacin a la Astronoma, hay muchos campos en los que es til aplicarla.

Figura 3. Curva de luz de una estrella variable peridica.

Si hemos realizado 200 observaciones de una estrella variable irregular, queremos saber qu brillo va a tener el mircoles de la semana que viene. Si hemos calculado el nmero de Wolf, querramos saber cul va a ser el mes o el ao que viene. Si hemos accedido a la web de la NASA y nos hemos trado una tabla de temperaturas en Marte, querramos ver qu pauta sigue su variacin para establecer un modelo. Hay que apresurarse a decir que la TF no resuelve del todo estos problemas: es un medio de aproximarse a ellos, mucho ms eficaz que otros. Por ejemplo, en el caso de una variable como la del grfico 3 podramos intentar realizar un ajuste polinmico. El error que se comete en el mejor de los casos, con un polinomio grandecito de grado 6, es muy superior al que se comete calculando la serie de Fourier. Si lo intentamos con U Geminorum, con una grfica como la de la figura 4, no vamos a tener tanta suerte para encontrar un buen ajuste polinmico (ntese que las dos grficas son diferentes). Necesitamos algo operativo y ms o menos automtico.

Figura 4. Curva de luz de U Gem.

Temperatura de la tierra Un problema sencillo pero muy interesante es el de calcular la temperatura de la tierra a una profundidad x a partir de la temperatura de la supercie. Describamos la temperatura de la supercie terrestre como una funcion f periodica en el tiempo t y de perodo 1 (un ano). La temperatura u(t,x) en el tiempo t 0 y profundidad x 0 es tambien periodica en t y es n atural asumir que |u|kfk. Bajo estas circunstancias u(t,x) puede ser expandida mediante una serie de Fourier para cada 0 x < jo como sigue: u(t,x) =X nZ cn(x)e2int, con coecientes de Fourier cn(x) =Z 1 0 u(t,x)e2int dt. Sabemos que la funcion u satisface la ecuacion diferencial parcial u t = 1 22u x2(Ecuacion del calor).

Por lo tanto, c00 n =Z 1 0 2u x2e2int dt = 2Z 1 0 u te2int dt = 4incn. En otras palabras, los coecientes cn satisfacen la ecuacion c00 n = [(2|n|)1/2(1i)]2cn, tomando el signo positivo o negativo de acuerdo a si n > 0 o n < 0. Por otra parte, sabemos que cn(0) = R1 0 f(t)e2int dt = f(n). Resolviendo la ecuacion, obtenemos que: cn(x) = f(n)exp[(2|n|)1/2(1i)x], 54 Genaro Gonzalez y por lo tanto resulta nalmente u(t,x) =X nZ f(n)exp[(2|n|)1/2x]exp[2int(2|n|)1/2ix]. Supongamos por ejemplo que la temperatura de la supercie viene dada por una funcion sinusoidal simple f(t) = sin(2t) (lo cual signica que la temperatura anual media f(0) =R1 0 f es cero). En este caso, la funcion u vendra dada por: u(t,x) = exp( 2x)sin(2t2x). Esta formula nos dice que la temperatura a la profundidad x =p 2 queda afectada por el factor e y esta completamente fuera de fase con respecto a las estaciones como lo indica la siguiente gura. invierno verano u(t,p/2) 1/2 1/2 t

CONCLUSIONES

La importancia de la transformada de Fourier radica en que permite representar funciones complicadas de forma que presenten propiedades muy tiles, que a menudo facilitan el tratamiento de la funcin original. La visualizacin simultanea de una funcin y su transformada de Fourier es, en muchas ocasiones, la clave del xito para solucionar problemas.

REFERENCIAS http://fundacion-eticotaku.org/2010/11/05/senales-convergentes-latransformada-de-fourier/ http://www.astro-digital.com/8/fourier.html www.emis.de/journals/DM/v5/art6.pdf