UNIVERSIDAD NACIONAL DE CAJAMARCA
Norte de la universidad peruana
FACULTAD DE INGENIERÍA
ESCUELA ACADÉMICO PROFESIONAL DE INGENIERÍA CIVIL
Filial Jaén
Tema: SERIES DE FOURIER
Curso: ANÁLISIS MATEMÁTICO III
Profesor: ELADIO SÁNCHEZ CULQUI
Jaén – Perú
INDICE
1. DEFINICIÓN ............................................................................................................ 3
2. APLICACIONES ...................................................................................................... 3
2.1. El problema isoperimétrico ................................................................................ 3
2.2. Temperatura de la tierra ..................................................................................... 4
2.3. Solución de ecuaciones diferenciales ................................................................ 5
2.4. Flujo del calor .................................................................................................... 5
2.5. Ecuación de ondas ............................................................................................. 6
2.6. Identidad de Jacobi ............................................................................................ 6
3. FUNCIONES PARES E IMPARES ......................................................................... 7
4. SIMETRIA DE MEDIA ONDA ............................................................................... 9
5. FUNCIONES PERIÓDICAS .................................................................................. 10
6. ANÁLISIS ARMÓNICO ........................................................................................ 13
6.1. EJERCICIO APLICATIVO ............................................................................... 13
SERIES DE FOURIER
1. DEFINICIÓN
Las series de Fourier son una familia de funciones ortogonales y periódicas.
Sea 𝑓(𝑥) una función de una variable real. Supongamos que dicha función es integrable
en un determinado intervalo de longitud T. Se define la serie de Fourier de 𝑓(𝑥) como:
𝑎0 ∞
𝑓(𝑥 ) = + ∑ (𝑎𝑛 cos(𝑛𝜔𝑥 ) + 𝑏𝑛 sin(𝑛𝜔𝑥 ))
2 𝑛=1
Donde:
2𝜋
𝜔= frecuencia fundamental.
𝑇
𝑎0 , 𝑎𝑛 , 𝑏𝑛 son los coeficientes de Fourier. Debemos tener en cuenta que
𝑎𝑛 𝑦 𝑏𝑛 hacen referencia a infinitos términos.
2. APLICACIONES
Las series de Fourier tiene una variedad de aplicaciones en las diferentes ramas de la
ingeniería, la matemática y la fisica-matematica, además de ser una herramienta
sumamente útil en la teoría matemática abstracta.
Áreas de aplicación incluyen análisis vibratorio, acústica, óptica, procesamiento de
imágenes y señales. En ingeniería, para el caso de los sistemas de telecomunicaciones y
a través del uso de los componentes espectrales de frecuencia de una señal dada, se puede
optimizar el diseño de un sistema para la señal portadora del mismo.
2.1. El problema isoperimétrico
Teorema 1. Si C es una curva cerrada simple 𝐶 1 y de longitud 1, entonces el área A
1
encerrada por C satisface la desigualdad 𝐴 ≤ 4𝜋. La se satisface si y sólo si C es una
circunferencia. En consecuencia, entre todas las curvas cerradas simples de longitud 1
la que encierra mayor área es la circunferencia.
2.2. Temperatura de la tierra
Un problema sencillo pero muy importante es calcular la temperatura de la tierra a una
profundidad x a partir de la temperatura de la superficie.
Describamos la temperatura de la superficie terrestre como una función 𝑓 periódica en el
tiempo t y de periodo 1 (un año). La temperatura 𝑢(𝑡, 𝑥) en el tiempo 𝑡 ≥ 0 y profundidad
𝑥 ≥ 0 es también periódica en t y es natural asumir que |𝑢| ≤ ‖𝑓‖∞ . Bajo estas
circunstancias 𝑢(𝑡, 𝑥) puede ser expandida mediante una serie de Fourier para cada
0 ≤ 𝑥 < ∞ fijo como sigue:
𝑢(𝑡, 𝑥) = ∑ 𝐶𝑛 (𝑥)𝑒 2𝜋𝑖𝑛𝑡
𝑛∈ℤ
Con coeficientes de Fourier
1
𝐶𝑛 (𝑥) = ∫ 𝑢(𝑡, 𝑥)𝑒−2𝜋𝑖𝑛𝑡 𝑑𝑡
0
Sabemos que la función u satisface la ecuación diferencial parcial
𝜕𝑢 1 𝜕2 𝑢
= ( 2 ) Ecuación del calor
𝜕𝑡 2 𝜕𝑥
Por lo tanto,
1 2 1 𝜕𝑢
𝜕 𝑢 −2𝜋𝑖𝑛𝑡
𝐶𝑛′′ = ∫ ( 2) 𝑒 𝑑𝑡 = 2 ∫ ( ) 𝑒−2𝜋𝑖𝑛𝑡 𝑑𝑡 = 4𝜋𝑖𝑛𝐶𝑛
0 𝜕𝑥 0 𝜕𝑡
En otras palabras, los coeficientes 𝐶𝑛 satisfacen la ecuación
2
𝐶𝑛′′ = [(2𝜋|𝑛|)1⁄2 (1 ± 𝑖)] 𝐶𝑛
Tomando el signo positivo o negativo de acuerdo a si 𝑛 > 0 ó 𝑛 < 0. Por otra parte,
1
sabemos, sabemos que 𝐶𝑛 (0) = ∫0 𝑓(𝑡) 𝑒−2𝜋𝑖𝑛𝑡 𝑑𝑡, se obtiene:
𝐶𝑛 (𝑥) = 𝑓̂(𝑛)𝑒𝑥𝑝[−(2𝜋|𝑛|)1⁄2 (1 ± 𝑖)𝑥]
Y por lo tanto resulta
𝑢(𝑡, 𝑥) = ∑ 𝑓̂(𝑛)𝑒𝑥𝑝[−(2𝜋|𝑛|)1⁄2 𝑥]𝑒𝑥𝑝[2𝜋𝑖𝑛𝑡 ∓ (2𝜋|𝑛|)1⁄2 𝑖𝑥]
𝑛∈ℤ
Supongamos que la temperatura de la superficie viene dad función sinusoidal simple
𝑓(𝑡) = sin(2𝜋𝑡)
𝑢(𝑡, 𝑥) = 𝑒𝑥𝑝(−√2𝜋𝑥) sin(2𝜋𝑡 − √2𝜋𝑥)
𝜋
Esta fórmula nos dice que la temperatura a la profundidad 𝑥 = √ 2 queda afectada por el
factor 𝑒 −𝜋 y está completamente fuera de fase con respecto a las estaciones.
2.3. Solución de ecuaciones diferenciales
Tal vez una de las propiedades más importantes de las integrales de Fourier es que
transforma operadores diferenciales con coeficientes constantes en multiplicación por
polinomios de acuerdo con la fórmula 𝑓̂ ′ = 2𝜋𝑖𝛾𝑓̂.
Veamos el siguiente ejemplo como resolver la ecuación diferencial
𝑢′′ − 𝑢 = −𝑓
En la cual f es una función conocida y debemos encontrar u. Aplicando el operador ^ en
ambos lados de la ecuación, obtenemos que (4𝜋 2 𝛾 2 + 1)𝑢̂ = 𝑓̂, equivalente a :
𝑢̂ = (4𝜋 2 𝛾 2 + 1)−1 𝑓̂
1
Pero (4𝜋 2 𝛾 2 + 1)−1 es la transformada de Fourier de la función 𝑒 −|𝑥| y por
2
consiguiente, tenemos que:
1
𝑢 = ∫ 𝑒 −|𝑥−𝑦| 𝑓(𝑦)𝑑𝑦
2
2.4. Flujo del calor
El problema del flujo del calor se describe mediante la ecuacion:
𝜕𝑢 1 𝜕 2 𝑢
= ( ), 𝑡 > 0, 𝑥 ∈ ℝ
𝜕𝑡 2 𝜕𝑥 2
Con condicion de borde lim 𝑢 = 𝑓. La solucion de este problema es similar al anterior,
𝑡=0
primero aplicamos la transformada de Fourier en ambos lados de la acuacion:
̂ ⁄𝜕𝑡 = −2(𝜋 2 𝛾 2 𝑢̂)
𝜕𝑢
Luego calculamos 𝑢̂ :
𝑢̂ = 𝑓̂𝑒𝑥𝑝(−2𝜋 2 𝛾 2 𝑡)
Finalmente invertimos y obtenemos
∞
𝑒𝑥𝑝[−(𝑥 − 𝑦)2 ⁄2𝑡]
𝑢(𝑡, 𝑥) = [𝑒𝑥𝑝](−2𝜋 2 𝛾 2 𝑡)𝑓̂ = 𝑝𝑡 𝑜 𝑓 = ∫ 𝑓(𝑦)𝑑𝑦
−∞ (2𝜋𝑡)1⁄2
(𝑥)2 ⁄2𝑡
Donde 𝑝𝑡 (𝑥) = es el llamado Kernel de Gauss.
(2𝜋𝑡)1⁄2
2.5. Ecuación de ondas
La ecuación de ondas viene dada por:
𝜕2𝑢 𝜕2𝑢
= , 𝑡 > 0, 𝑥∈ℝ
𝜕𝑡 2 𝜕𝑥 2
Con condicion de borde lim 𝑢 = 𝑓. Primero aplicaremos la transformada de Fourier:
𝑡=0
𝜕 2 𝑢̂
2
= −4𝜋2 𝛾2 𝑢̂
𝜕𝑡
Luego calculamos 𝑢̂ :
sin 2𝜋𝛾𝑡
𝑢̂(𝑡, 𝛾) = cos 2𝜋𝛾𝑡𝑓̂(𝛾) + 𝑔̂(𝛾)
2𝜋𝛾
1 1 𝑡
= [𝑒 2𝜋𝑖𝛾𝑡 + 𝑒 −2𝜋𝑖𝛾𝑡 ]𝑓̂(𝛾) + ∫ 𝑒 2𝜋𝑖𝛾𝑡 𝑑𝑦 𝑔̂(𝛾)
2 2 −𝑡
Luego invertimos para finalmente obtener:
1 1 𝑥+𝑡
𝑢(𝑡, 𝑥) = [𝑓(𝑥 + 𝑡) + 𝑓(𝑥 − 𝑡)] + ∫ 𝑔(𝑦)𝑑𝑦
2 2 𝑥−𝑡
Esta es la llamada fórmula de D’ Alembert.
2.6. Identidad de Jacobi
La función θ es una función trascendente que aparece en diferentes áreas de la matemática
como teoría de números, funciones elípticas, ecuaciones de la física matemática y
mecánica estadística. Se define la siguiente suma.
∞
𝜃(𝑡) = ∑ exp(−𝜋𝑛2 𝑡), para 𝑡 > 0
𝑛=−∞
La Identidad de Jacobi establece que:
𝜃(𝑡) = 𝑡 −1⁄2 𝜃(1⁄𝑡), para 𝑡 > 0
Demostración: Consideremos el Kernel de Gauss
1 2 ⁄2𝑡
𝑝𝑡 (𝑥) = 𝑒 −𝑥
√2𝜋𝑡
Entonces:
2 𝛾2 𝑡
𝑝̂𝑡 (𝛾) = 𝑒 −2𝜋
Y aplicado la formula de Poisson, obtenemos que:
1 2 ⁄2𝑡 2 𝑛2 𝑡
∑ 𝑒 −𝑛 = ∑ 𝑒 −2𝜋
𝑛∈ℤ
√2𝜋𝑡 𝑛𝜖ℤ
Luego, haciendo el cambio 𝑡 → 𝑡⁄2𝜋, obtenemos que:
1 2 ⁄𝑡 2𝑡
∑ 𝑒 −𝜋𝑛 = ∑ 𝑒 −𝜋𝑛
√𝑡 𝑛𝜖ℤ 𝑛∈ℤ
Es decir:
1 1
𝜃(𝑡) = 𝜃( )
√𝑡 𝑡
3. FUNCIONES PARES E IMPARES
Diremos que la función f : R R , es par si satisface la condición f ( t ) f ( t ), t R
y es impar si satisface la condición f ( t ) f ( t ), t R .
OBSERVACIÓN: Observamos que una función par, es simetría respecto del eje vertical
en el origen, mientras que una función impar es antisimétrica respecto del eje vertical en
el origen.
3.1. PROPIEDADES DE LAS FUNCIONES PARES E IMPARES.
1. Si f, g : R R son funciones pares entonces f.g es función par.
Demostración
Sea h ( t ) f ( t ).g( t ) , por demostrar que h ( t ) es función par
h ( t ) f ( t ).g( t ) f ( t ).g( t ) h ( t ) Por lo tanto
Como h ( t ) h ( t ) entonces f ( t ).g( t ) es función par
2. Si f, g : R R son funciones impares entonces f.g es una función par.
Demostración
f , g Son funciones impares f ( t ) f ( t ) y g( t ) g( t )
Sea h ( t ) f ( t ).g( t ) , por demostrar que h ( t ) es función par
h ( t ) f ( t ).g( t ) (f ( t ).( g( t )) f ( t ).g( t ) h ( t )
Como h ( t ) h ( t ) entonces f ( t ).g( t ) es función par.
3. Si f, g : R R es una función par y g : R R es una función impar entonces f.g
es una función impar.
Demostración
Como f es una función par, entonces f ( t ) f ( t ) y como g es una función impar,
entonces g( t ) g( t ) .
Sea h ( t ) f ( t ).g( t ) por demostrar que h ( t ) es función impar.
h ( t ) f ( t ).g( t ) f ( t ).( g( t )) f ( t ).g( t ) h ( t )
Como h ( t ) h ( t ) , entonces f ( t ).g( t ) es función impar.
4. Toda función f ( t ) , de puede expresar como la suma de dos funciones
componentes, de las cuales la una es par y la otra es impar.
Demostración
Se conoce que cualquier función f ( t ) se puede expresar así:
1 1 1 1 1 1
f ( t ) f ( t ) f ( t ) f ( t ) f ( t ) [f ( t ) f ( t )] [f ( t ) f ( t )] (1)
2 2 2 2 2 2
f e ( t ) 2 f ( t ) f ( t )
1
Sean
f ( t ) 1 [f ( t ) f ( t )]
0 2
1 1
Donde f e ( t ) [f ( t ) f ( t )] [f ( t ) f ( t )] f e ( t )
2 2
Por lo tanto f e ( t ) es función par
1 1
f 0 ( t ) [f ( t ) f ( t )] [f ( t ) f ( t )] f 0 ( t )
2 2
Por lo tanto f 0 ( t ) es función impar, entonces f (t ) f e (t ) f 0 (t ) , donde f e (t ) es
la componente par y f 0 ( t ) es la componente impar de la función f ( t ) .
4. SIMETRIA DE MEDIA ONDA
Si la función f ( t ) es periódica con periodo T, entonces se dice que la función f ( t ) tiene
T
simetría de media onda si satisface la condición: f ( t ) f ( t )
2
OBSERVACIÓN:
Mostraremos en la figura una forma de onda con simetría de media onda se debe observar
que la porción negativa de la onda es el reflejo de la porción positiva desplazando
horizontalmente medio periodo.
NOTA. - Simetría de media onda no es par ni impar.
TEOREMA.- Si una función negativa f ( t ) tiene simetría de media onda. Demostrar que:
1
f ( t ) f ( t T )
2
Demostración
1
Como f (t) tiene simetría de media onda, entonces se tiene: f ( t ) f ( t T) (1)
2
Como (t) es periódica con periodo T, entonces:
1 1 T
f ( t T) f ( t T T) f ( t ) (2)
2 2 2
1
Por lo tanto de (1) y (2) se tiene: f ( t ) f ( t T)
2
5. FUNCIONES PERIÓDICAS
Una función 𝑓(𝑡) tiene un período T o es periódica con período T si para todo t, 𝑓(𝑡 +
𝑇) = 𝑓(𝑡), donde T es una constante positiva. El valor más pequeño de 𝑇 > 0 se llama
el período principal o período fundamental o simplemente el período de 𝑓(𝑡).
Observación:
1) Si f y g son periódicas con periodo T, entonces h = af + bg, a y b constantes,es
también periódica con periodo T.
2) Si f y g son periódicas con periodo T, entonces f. g es también periódica con
perido T.
TEOREMA:
Si f t es una función periódica con periodo T, entonces se cumple:
T T
a
2 2
I) f t dt f t dt
T T
a
2 2
T t t
II) f t dt f t dt
T 0
Caso como en la función 𝑠𝑒𝑛𝑥 y cosx tiene períodos 2π, 4π, 6π ya que:
𝑠𝑒𝑛(𝑥 + 2𝜋), 𝑠𝑒𝑛(𝑥 + 4𝜋), 𝑠𝑒𝑛(𝑥 + 6𝜋), … Todas son iguales a 𝑠𝑒𝑛𝑥.
𝑐𝑜𝑠(𝑥 + 2𝜋), 𝑐𝑜𝑠(𝑥 + 4𝜋), 𝑐𝑜𝑠(𝑥 + 6𝜋), … Todas son iguales a 𝑐𝑜𝑠𝑥.
6. SERIE DE FOURIER PARA UNA FUNCIÓN DE PERIODO 2L
En muchas ocasiones es deseable adaptar la forma de una serie de Fourier a funciones
periódicas de periodo P 2L 0 en el intervalo [−L.L[. Esto se consigue gracias a un
cambio de variable.
Sea f (t ) una función periódica de periodo 2L. Para desarrollar en serie de Fourier en
[−L L[ hacemos un cambio de variable, poniendo.
t L
x
Entonces
L
f (t ) f ( x)
Si definimos
L
g ( x) f ( x)
claramente la función g es una función periódica de x de periodo 2 ya que
L L L
g ( x 2 ) f ( ( x 2 )) f ( x 2 L) f ( x) g ( x)
De esta forma si el desarrollo en serie de Fourier de la función g(x) es
L a0
g ( x) f ( x) (an cos(nx) bn sin(nx))
2 n 1
con coeficientes de Fourier
x t
entonces, como L sustituyendo
w0
Generalmente se escribe L y por lo tanto
7. ANÁLISIS ARMÓNICO
El punto de partida del análisis armónico fue el desarrollo en serie de Fourier. Afirma que
cualquier función 1-periódica razonable 𝑓: ℝ → ℂ se puede realizar como
∞ 1
𝑓(𝑥) = ∑ 𝑎𝑛 𝑒(𝑛𝑥) 𝑐𝑜𝑛 𝑎𝑛 = ∫ 𝑓(𝑥)𝑒(−𝑛𝑥)𝑑𝑥
𝑛=−∞ 0
donde se ha empleado la notación sintética 𝑒(𝑥) = 𝑒 2𝜋𝑖𝑥 para englobar los senos y
cosenos que representan los tonos puros. A los 𝑎𝑛 se les llama coeficientes de Fourier.
Las funciones {𝑒(𝑛𝑥)}∞
𝑛=−∞ son los armónicos en esta descomposición.
Por la periodicidad, el intervalo de integración [0, 1] se puede sustituir por cualquier otro
de longitud 1. A este respecto, en muchos contextos se muestra conveniente introducir el
toro 𝕋 = ℝ/{𝑥 → 𝑥 + 1} obtenido al enrollar la recta real alrededor de la circunferencia
unidad (un toro unidimensional es una circunferencia), o equivalentemente pegar los
extremos del intervalo [0,1]. En este sentido, 𝑓: 𝕋 → ℂ representa una función 1-periodica
y 𝑎𝑛 viene dado por la integral 𝕋.
7.1. EJERCICIO APLICATIVO
1. Determinar las componentes armónicos de:
Solución
𝜔𝑇 = 2𝑎 + 2𝑏 = 𝑃
2𝜋
𝑇 = 2𝑎 + 2𝑏 → 𝑎 + 𝑏 = 𝜋
𝑇
∞ ∞
𝑓(𝑥) = 𝑎0 + ∑ 𝑎𝑛 . 𝑐𝑜𝑠(𝑛𝑥) + ∑ 𝑏𝑛 . 𝑠𝑖𝑛(𝑛𝑥)
𝑛=1 𝑛=1
1 𝑃 2 𝑃
𝑎0 = ∫ 𝑓(𝑥)𝑑𝑥 𝑎𝑛 = ∫ 𝑓(𝑥)cos(𝑛𝑥)𝑑𝑥 𝑏𝑛
𝑃 0 𝑃 0
𝑃
2
= ∫ 𝑓(𝑥)sin(𝑛𝑥)𝑑𝑥
𝑃 0
2𝑎+2𝑏
2
𝑏𝑛 = ∫ 𝑓(𝑥)sin(𝑛𝑥)𝑑𝑥
2𝑎 + 2𝑏 0
𝑏 3𝑏
+𝑎 2𝑎+
1 2 2
= [∫ 𝐸𝑠𝑖𝑛(𝑛𝑥)𝑑𝑥 + ∫ −𝐸𝑠𝑖𝑛(𝑛𝑥)𝑑𝑥]
𝑎 + 𝑏 𝑏/2 𝑎+
3𝑏
2
𝑏 3𝑏
+𝑎 2𝑎+
𝐸 −1 2 −1 2
𝑏𝑛 = [[ cos(𝑛𝑥)] − [ cos(𝑛𝑥)] 3𝑏 ]
𝑎+𝑏 𝑛 𝑏/2 𝑛 𝑎+
2
𝑏 3𝑏
𝐸 2
+𝑎 2𝑎+
𝑏𝑛 = − [[cos(𝑛𝑥)]𝑏/2 − [cos(𝑛𝑥)] 3𝑏2 ]
(𝑎 + 𝑏)𝑛 𝑎+
2
−𝐸 𝑏 𝑏 3𝑏 3𝑏
𝑏𝑛 = [cosn ( + 𝑎) − 𝑐𝑜𝑠𝑛 − 𝑐𝑜𝑠𝑛 (2𝑎 + ) + 𝑐𝑜𝑠𝑛 (𝑎 + )]
𝜋𝑛 2 2 2 2
3𝑏 𝑏 𝑏
→ 𝑐𝑜𝑠𝑛 (2𝑎 + ) = 𝑐𝑜𝑠𝑛 [(𝑎 + 𝑏) + ( + 𝑎)] = 𝑐𝑜𝑠 [𝑛𝜋 + 𝑛 ( + 𝑎)]
2 2 2
3𝑏 𝑏 𝑏
→ 𝑐𝑜𝑠𝑛 (𝑎 + ) = 𝑐𝑜𝑠𝑛 [(𝑎 + 𝑏) + ( )] = 𝑐𝑜𝑠 [𝑛𝜋 + 𝑛 ( )]
2 2 2
Si n es par
−𝐸 𝑏 𝑏 𝑏 𝑏
𝑏𝑛 = [cosn ( + 𝑎) − 𝑐𝑜𝑠𝑛 − 𝑐𝑜𝑠𝑛 ( + 𝑎) + 𝑐𝑜𝑠𝑛 ( )]
𝜋𝑛 2 2 2 2
Si n es impar
−𝐸 𝑏 𝑏 𝑏 𝑏
𝑏𝑛 = [cosn ( + 𝑎) − 𝑐𝑜𝑠𝑛 + 𝑐𝑜𝑠𝑛 ( + 𝑎) − 𝑐𝑜𝑠𝑛 ( )]
𝜋𝑛 2 2 2 2
−2𝐸 𝑏 𝑏
𝑏𝑛 = [cosn ( + 𝑎) − 𝑐𝑜𝑠𝑛 ] … … … (1)
𝜋𝑛 2 2
𝑏 𝑎+𝑏 𝑎 𝑛𝜋 𝑎
→ cosn ( + 𝑎) = cosn ( + ) = cos ( + 𝑛 )
2 2 2 2 2
Llegamos a la conclusión que n es impar entonces:
𝑛𝜋 𝑎 𝑎
cos ( + 𝑛 ) = −𝑠𝑖𝑛 (𝑛 )
2 2 2
Luego al reemplazar en (1) obtenemos
−2𝐸 𝑎 𝑏
𝑏𝑛 = [−𝑠𝑖𝑛 (𝑛 ) − 𝑐𝑜𝑠𝑛 ]
𝜋𝑛 2 2
Por ultimo reemplazamos 𝑏𝑛 en la serie de Fourier:
∞ ∞
𝑓(𝑥) = 𝑎0 + ∑ 𝑎𝑛 . 𝑐𝑜𝑠(𝑛𝑥) + ∑ 𝑏𝑛 . 𝑠𝑖𝑛(𝑛𝑥)
𝑛=1 𝑛=1
∞
2𝐸 𝑎 𝑏
𝑓(𝜔𝑡) = ∑ [𝑠𝑖𝑛 (𝑛 ) + 𝑐𝑜𝑠𝑛 ] . 𝑠𝑖𝑛(𝑛𝜔𝑡)
𝜋𝑛 2 2
𝑛=1