ICF415

Sistemas continuos

Dr. Fabián Torres.

Departamento de Ciencias Fı́sicas

Universidad de La Frontera.

1/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 1 / 28
 Índice

1 Introducción

2 Ecuación de ondas
Soluciones en una cuerda

3 Serie de Fourier
Marco teórico

2/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 2 / 28
 Introducción

Introducción

En la clase anterior vimos como el análisis del movimiento en una cuerda

conducen a una ecuación diferencial en derivadas parciales llamada
Ecuación de ondas.
∂2f 1 ∂2f
2
= 2 2 (1)
∂x v ∂t
Hoy continuaremos con el análisis de sistemas continuos e introduciremos
una nueva serie que es de mucha utilidad en el análisis de señales...

3/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 3 / 28
 Ecuación de ondas Soluciones en una cuerda

Índice

1 Introducción

2 Ecuación de ondas
Soluciones en una cuerda

3 Serie de Fourier
Marco teórico

4/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 4 / 28
 Ecuación de ondas Soluciones en una cuerda

Introducción

Como vimos en la clase anterior, los desplazamientos transversales y(x, t) de

una cuerda se describen mediante la ecuación de onda
∂2y 1 ∂2y
2
= 2 2 (2)
∂x v ∂t
en donde la velocidad de propagación de las ondas en la cuerda está dada por
s
T
v= (3)
µ

con T la tensión en la cuerda, y µ la densidad lineal de masa de la cuerda.

5/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 5 / 28
 Ecuación de ondas Soluciones en una cuerda

Para solucionar esta ecuación en derivadas parciales recurrimos a un

método conocido como separación de variables.
Como la función y(x, t) depende de una parte espacial (x) y una parte
temporal (t), entonces para implementar el método de separación de
variables, suponemos que y se puede escribir como el producto de una
función que depende solamente de x y de una función que depende
solamente de t.
Supongamos por ejemplo que y(x, t) = f (x) · cos(ωt), en donde f (x) es la
función que depende exclusivamente de x y la función de t es cos(ωt)

6/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 6 / 28
 Ecuación de ondas Soluciones en una cuerda

Derivando estas expresiones con respecto a x y a t se tiene

∂y df (x) ∂2y d2 f (x)

= cos(ωt), 2
= cos(ωt) (4)
∂x dx ∂x dx2
∂y ∂2y
= −f (x)ω sin(ωt), = −f (x)ω 2 cos(ωt) (5)
∂t ∂t2
Reemplazando estas expresiones en la ecuación diferencial se tiene

∂2y 1 ∂2y
= (6)
∂x2 v 2 ∂t2
2
d f (x) ω2
2
cos(ωt) = − 2 f (x) cos(ωt) (7)
dx v
d2 f (x) ω2
= − 2 f (x) (8)
dx2 v

7/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 7 / 28
 Ecuación de ondas Soluciones en una cuerda

Reordenando la expresión anterior, obtenemos una ecuación de oscilador

armónico
d2 f (x) ω 2
+ 2 f (x) = 0 (9)
dx2 v
solo que en lugar de depender armónicamente del tiempo, depende
armónicamente de la coordenada x.
Como esta ecuación es la de un oscilador armónico, entonces sabemos que
la solución para f (x) es

f (x) = A cos(ω ′ x + φ) (10)

con ω ′ = ωv
Como la cuerda que estamos considerando está anclada por ambos lados a
muros, entonces las condiciones de bordes son

f (x = 0) = f (x = L) = 0 (11)

es decir, en los extremos la cuerda no se mueve nunca.

8/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 8 / 28
 Ecuación de ondas Soluciones en una cuerda

Utilizando la primera condición de borde (f (x = 0) = 0), se tiene entonces que

π
cos(φ) = 0, ∴ φ= (12)
2
Como φ = π2 entonces podemos expresar f (x) = A sin(ω ′ x)
Para la segunda condición de borde (f (x = L) = 0) se tiene que

sin(ω ′ L) = 0, ∴ ω ′ L = nπ, n = 1, 2, 3, ... (13)

Si en lugar de usar la frecuencia angular, usamos la frecuencia de oscilación

ω
ν = 2π , entonces vemos que las frecuencias de oscilación de la cuerda son

ω 2✚
π νL
ω′ L = L= = n✚
π (14)
v v s
nv n T
ν= = (15)
2L 2L µ

de modo que las frecuencias son discretas y numerables a través de n. Usaremos

entonces n como un ı́ndice en las cantidades que dependen de el.
9/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 9 / 28
 Ecuación de ondas Soluciones en una cuerda

Longitud de onda

Podemos aprovechar de introducir el concepto de Longitud de Onda asociada

al modo n como
2L
λn = , n = 1, 2, 3, ... (16)
n
Con esta expresión, las frecuencias de oscilación se pueden expresar como
nv v ω 2π
νn = = , v = λn ν n , ó ω ′ = = (17)
2L λn v λn

De modo que la expresión para f (x) se puede escribir como

 
2π  nπ 
f (x) = An sin x = An sin x (18)
λn L

10/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 10 / 28
 Ecuación de ondas Soluciones en una cuerda

Finalmente, la solución para el desplazamiento transversal de la cuerda y(x, t)

se puede expresar como

y(x, t) = f (x) cos(ωt) (19)

s
 nπ  nπ T
= An sin x cos(ωn t), ωn = (20)
L L µ

De la expresión para ωn se puede ver que todas las frecuencias angulares se

pueden escribir como multiplos enteros de una frecuencia fundamental ω1
s !
π T
ωn = n = nω1 (21)
L µ

de esta forma definimos la frecuencia angular fundamental ω1 de oscilación de

la cuerda como
s
π T
ω1 = (22)
L µ

11/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 11 / 28
 Ecuación de ondas Soluciones en una cuerda

Modos Normales

Como ya se conocen las frecuencias propias de oscilación de la cuerda (equiva-

lente a los autovalores), podemos determinar los modos normales asociados a
cada una de las frecuencias angulares de oscilación .

Definimos los nodos como los puntos en

los que la amplitud siempre es cero.
Vemos que para:
n = 1 existen dos nodos.
n = 2 existen tres nodos.
n = 3 existen cuatro nodos. etc...

12/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 12 / 28
 Ecuación de ondas Soluciones en una cuerda

Condiciones Iniciales

Si bien vimos que podemos fijar condiciones de borde para la función f (x), no
hemos dicho nada sobre las condiciones iniciales asociada a la parte temporal
de la solución.
Las condiciones iniciales corresponden en este caso a la forma (o función) que
describe a la cuerda en el instante de partida del movimiento.

13/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 13 / 28
 Ecuación de ondas Soluciones en una cuerda

Como vemos la forma inicial que toma la cuerda es muy poco probable
que tome una forma sinusoidal. Sin embargo, luego de la etapa transiente
el sistema continua moviéndose como una superposición de los modos
normales de oscilación mas bajos.
Esto quiere decir que las vibraciones en modos normales pueden ocurrir
simultáneamente e independientemente uno de otro.
Esto se debe a que la ecuación diferencial de ondas es una ecuación lineal.

IMPORTANTE!!
Dado que se conocen los modos normales y1 , y2 , y3 , ... asociados a las
frecuencias ν1 , ν2 , ν3 , ..., los cuales son soluciones de la EDP lineal de ondas,
entonces y1 + y2 + y3 + ... también es solución.
Por este motivo, las funciones y1 , y2 , y3 , ... se pueden utilizar como una base
(de un espacio vectorial) para escribir cualquier solución como una
combinación lineal de estas
La expresión de cualquier función en la base mencionada se conoce como
Serie de Fourier
14/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 14 / 28
 Serie de Fourier Marco teórico

Outline

1 Introducción

2 Ecuación de ondas
Soluciones en una cuerda

3 Serie de Fourier
Marco teórico

15/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 15 / 28
 Serie de Fourier Marco teórico

Introducción

Dentro del mundo de las matemáticas, el álgebra lineal es una de las

ramas mas importantes para la fı́sica.
El estudio espacios y subespacios vectoriales, junto con el estudio de bases
de esos espacios permite el desarrollo de técnicas de aproximación
ampliamente utilizadas tanto en Fı́sica como en programación
computacional.
Ejemplo de esto son por ejemplo el uso de matrices computacionales y
algoritmos de descomposición de estas,Series de Taylor y en particular el
tópico de interés en esta ocasión, las Series de Fourier.

16/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 16 / 28
 Serie de Fourier Marco teórico

Base de la serie de Fourier

Para las series de Fourier en una dimensión, los elementos de la base son
los modos normales de oscilación de una cuerda.
Cualquier función periódica f (x) se puede escribir como una combinación
lineal de los modos de oscilación de una cuerda.
De forma mas general, los modos normales de una cuerda (sin la
dependencia temporal) se pueden expresar como las funciones cos ωn x y
sin ωn x. Una propiedad importante de estas funciones (y la base del
desarrollo de Fourier) es que son ortogonales entre si, y para distintos
valores de n.

17/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 17 / 28
 Serie de Fourier Marco teórico

Desarrollo de Fourier
Cualquier función periódica f (x) = f (x + L), de periodo L, se puede
descomponer en una suma de senos y cosenos de la forma
∞ h
X  x  x i
f (x) = a0 + an cos 2πn + bn sin 2πn (23)
n=1
L L

Dado que el desarrollo de Fourier no está acotado a funciones espaciales, si se

desea aplicar a funciones periódicas en el tiempo, entonces se puede reemplazar
x t
L por T , donde T es el periodo de la función periódica temporal.

18/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 18 / 28
 Serie de Fourier Marco teórico

La pregunta natural que surge es entonces ¿Cómo determino los coeficientes a0 ,

an y bn de la expansión de Fourier?
Debemos recurrir a las propiedades de ortogonalidad de las funciones trigo-
nométricas.

19/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 19 / 28
 Serie de Fourier Marco teórico

Ortogonalidad

Las funciones trigonométricas forman una base en el espacio de funciones pe-

riódicas, es decir son linealmente independientes, y expanden todo el espacio.
Un punto importante a tener en mente, es que las funciones utilizadas son
simétricas con respecto a f (x) = 0, es decir
Z L Z L2
2  x  x
cos 2πn dx = sin 2πn dx = 0 (24)
−L
2
L −L 2
L

Por este motivo es necesario el término a0 .

20/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 20 / 28
 Serie de Fourier Marco teórico

Ejercicio

Calcule los valores de las siguientes integrales:

Z L
2  x  x
I0 = cos 2πn sin 2πm dx (25)
−L
2
L L
Z L
2  x  x
I1 = cos 2πn cos 2πm dx (26)
−L
2
L L
Z L
2  x  x
I2 = sin 2πn sin 2πm dx (27)
−L
2
L L

21/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 21 / 28
 Serie de Fourier Marco teórico

Solución:
Z L
2 x  x
I0 = cos 2πn sin 2πm dx = 0, ∀x, m, n (28)
−L2
L L
Z L2  L
 x  x si m = n
I1 = cos 2πn cos 2πm dx = 2 (29)
−L L L 0 si m =
6 n
2
Z L2  L
 x  x si m = n
I2 = sin 2πn sin 2πm dx = 2 (30)
−L L L 0 si m =
6 n
2

22/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 22 / 28
 Serie de Fourier Marco teórico

Usando estas integrales (propiedades de ortogonalidad entre funciones trigo-

nométricas), podemos entonces obtener los coeficientes a0 , an y bn como:

Valor de a0

∞ h
X  x  x i
f (x) = a0 + an cos 2πn + bn sin 2πn (31)
n=1
L L
L L
" ∞ h
#
✘x✘ x✘
✘
Z 2
Z 2
✘ ✘
X    i
f (x)dx = a0 + a cos ✘ 2πn sin ✘
+ bn✘ 2πn dx (32)
−L −L ✘n✘✘ L ✘ ✘ L
2 2 n=1

Por lo tanto

L
1
Z 2
a0 = f (x)dx (33)
L −L
2

23/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 23 / 28
 Serie de Fourier Marco teórico

Valor de an
2πn 2πm
Consideremos ωn = L y ωm = L , entonces
∞
X
f (x) = a0 + [an cos(ωn x) + bn sin(ωn x)]
n=1
Z L Z L
" ∞
#
2 2
✭✭m✭✭ X
f (x) cos(ωm x)dx = cos(ω
a0✭
✭ x) + an cos(ωn x) cos(ωm x)dx
−L
2 −L
2 n=1
L ∞
✭✭ ✭ ✭
Z 2
x)✭
X
+ sin(ω
bn✭ ✭ ✭ n✭ cos(ω m x) dx
−L
✭
2 n=1
Z L ∞ Z L Z L ∞
2 X 2 2 X L
f (x) cos(ωm x)dx = an cos(ωn x) cos(ωm x)dx = an δmn
−L
2 n=1 −L
2 −L
2 n=1
2
Z L
2 L
f (x) cos(ωm x)dx = am
−L
2
2

24/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 24 / 28
 Serie de Fourier Marco teórico

Ejemplo

Por lo tanto
L L  
2 2 2πn
Z 2
Z 2
an = f (x) cos(ωn x)dx = f (x) cos x dx (34)
L −L
2
L −L
2
L

25/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 25 / 28
 Serie de Fourier Marco teórico

Valor de bn
2πn 2πm
Consideremos nuevamente ωn = L y ωm = L , entonces
∞
X
f (x) = a0 + [an cos(ωn x) + bn sin(ωn x)]
n=1
" #
✭✭
L L ∞
✭✭✭
Z Z
2 2
✭✭ ✭✭ X
✭ ✭
f (x) sin(ωm x)dx = sin(ω
a0✭
✭ m x) + cos(ω
an✭ ✭ x) sin(ω x)dx
−L −L ✭ ✭ n m
2 2 n=1
Z L ∞
2 X
+ bn sin(ωn x) sin(ωm x)dx
−L
2 n=1
Z L ∞ Z L Z L ∞
2 X 2 2 X L
f (x) sin(ωm x)dx = bn sin(ωn x) sin(ωm x)dx = bn δmn
−L
2 n=1 −L
2 −L
2 n=1
2
Z L
2 L
f (x) sin(ωm x)dx = bm
−L
2
2

26/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 26 / 28
 Serie de Fourier Marco teórico

Ejemplo

Por lo tanto

L L
2 2 2π2
Z 2
Z 2
bn = f (x) sin(ωn x)dx = f (x) sin( x)dx (35)
L −L
2
L −L
2
L

27/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 27 / 28
 Serie de Fourier Marco teórico

Ejemplo

Determine la serie de Fourier de la función

 1
− 2 −1 < x < 0
f (x) = 1 (36)
2 0<x<1

28/28
F. A. Torres-Ruiz (FICA-UFRO) ICF415 Sistemas continuos 28 / 28