Funciones

Ortogonales y 

Series de Fourier
 Contenidos
• 12.1 Funciones Ortogonales
• 12.2 Series de Fourier
• 12.3 Series de Fourier de Cosenos y Senos
• 12.4 Series d eFourier Complejas
• 12.5 Problema de Sturm-Liouville
• 12.6 Series de Bessel y Legendre
 12.1 Funciones Ortogonales
DEFINICIÓN 12.1
Productos Interiores de Funciones
El producto interior de dos funciones f1 y f2 en un
intervalo [a, b] es el número

b
( f1, f 2 ) = ∫ f1 ( x) f 2 ( x)dx
a

DEFINICIÓN 12.2
Funciones Ortogonales
Dos funciones f1 y f2 son ortogonales en un
intervalo [a, b] si

b
( f1 , f 2 ) = ∫ f1 ( x) f 2 ( x) dx = 0
a
 Ejemplo
• Las funciones f1(x) = x2, f2(x) = x3 son
ortogonales en el intervalo [−1, 1] puesto que

1
1 2 3 1 6
( f1 , f 2 ) = ∫ x x dx = x =0
−1 6 −1
• Verificar si las funciones f1(x) = x3, f2(x) = x2 + 1
son ortogonales en el intervalo [−1, 1]
 DEFINICIÓN 12.3
Conjunto Ortogonal
Se dice que un conjunto de funciones de valores reales
{φ0(x), φ1(x), φ2(x), …} es ortogonal en un intervalo [a, b]
si
b
( φm , φn ) = ∫ φ m ( x), φ n ( x) dx
a
= 0, m ≠ n (2)

 Conjuntos Ortonormales

• La expresión (u, u) = ||u||2 se llama norma

cuadrada. Por tanto podemso definir la norma
cuadrada de una función como



2 b 2 b 2 (3)

φn ( x) = ∫ φn dx , φn ( x) = ∫a φn ( x)dx

 a
Si {φn(x)} es un conjunto ortogonal en [a, b] con
la propiedad de que ||φn(x)|| = 1 para todo n,
entonces se llama conjunto ortonormal en [a, b].
 Ejemplo 1
Demuestre que el conjunto {1, cos x, cos 2x, …}
es ortogonal en [−π, π].
Solución 

Sea φ0(x) = 1, φn(x) = cos nx, comprobamos que

π π
(φ0 , φn ) = ∫ φ0 ( x)φn ( x)dx = ∫ cos nxdx
−π −π
1 π
= sin nx −π = 0, for n ≠ 0
n
 Ejemplo 1 (2)
y π π
(φm ,φn ) = ∫ φm ( x)φn ( x)dx = ∫ cos mx ⋅ cos nxdx
−π −π
1 π
= ∫ [cos(m + n) x + cos(m − n) x]dx
2 −π
π
1 'sin( m + n) x sin( m − n) x $
= % + = 0, m ≠ n
2& m+ n m − n "# −π
 Ejemplo 2
Determine la norma de cada función del Ejemplo 1.
Solución
2 π
φ0 = 1, φ0 = ∫ dx = 2π ⇒ φ0 = 2π
−π

φn = cos nx,
2 π 2 1 π
|| φn || = ∫ cos nxdx = ∫ (1 + cos 2nx)dx = π
−π 2 −π
|| φn ||= π , n > 0
 Analogía con Vectores
• Recordando de la teoría de vectores en 3
dimensiones que

u = c 1v1 + c 2 v 2 + c 3 v 3 , (4)

tenemos


 (u, v ) ( u , v ) (u , v ) 3
(u, v n )
1 2 3
u = 2 v1 + v +
2 2
v = ∑
2 3 v (5)

2 n
|| v1 || || v 2 || || v 3 || n =1 || v n ||


Así podemos hacer una analogía entre funciones y
vectores.
 Desarrollo en Series Ortogonales
• Suponga que {φn(x)} es un conjunto ortogonal en [a,
b]. Si f(x) está definida en [a, b], escribimos primero


f ( x) = c0φ0 ( x) + c1φ1 ( x) + ! + cnφn ( x) + !?
(6)
b
Then ∫a f ( x)φm ( x) dx
b b
= c0 ∫ φ0 ( x)φm ( x) dx + c1 ∫ φ1 ( x)φm ( x) dx + !
a a
b
+ cn ∫ φn ( x)φm ( x) dx + !
a

= c0 (φ0 , φm ) + c1 (φ1 , φm ) + ! + cn (φn , φm ) + !

• Como {φn(x)} es un conjunto ortogonal en [a, b], cada
término en el lado derecho es nulo, excepto m = n. En
este caso tenemos

b b 2
∫a f ( x)φn ( x)dx = cn (φn ,φn ) = cn ∫a φn ( x)dx
b
∫ f ( x)φn ( x)dx
cn = a b , n = 0,1,2,...
2
∫a φn ( x)dx
 En otras palabras,


 ∞
f ( x) = ∑ cnφn ( x ), (7)


 n =0
b

 f ( x)φn ( x) dx
cn = a
∫ (8)

2

 || φn ( x) ||
Entonces (7) se transforma en 


 ∞
( f , φn )
f ( x) = ∑ 2 φn ( x ) (9)
n =0 || φn ( x ) ||
 DEFINICIÓN 12.4
Conjunto Ortogonal y Función Peso
Se dice que un conjunto de funciones de valores reales
{φ0(x), φ1(x), φ2(x), …} es ortogonal con respecto a una
función peso w(x) en [a, b], si 


 b
∫a w( x)φm ( x)φn ( x)dx = 0, m ≠ n
• Bajo la condición de la definición anterior, tenemos 


 b
f ( x) w( x)φn ( x) dx
c = a ∫ (10)

n

 || φn ( x) ||2
2 b 2 (11)
|| φn ( x) || = ∫ w( x)φn ( x) dx
a
 Conjuntos Completos
• Un conjunto ortogonal es completo si la única
función ortogonal continua a cada miembro
del conjunto es función nula.
 12.2 Series de Fourier
• Una Serie Trigonométrica

Podemos demostrar que el conjunto

&
 π 2π 3π π 2π 3π #
1
% , cos x , cos x , cos x , … , sin x , sin x , sin x , …
$
p p p p p p " (1)

!


es orthogonal en [−p, p]. Así una función f
definida en [−p, p] puede escribirse como


 a0 ∞ ( nπ nπ %

 f ( x) = + ∑ & an cos x + bn sin x#
2 n=1 ' p p $
(2)
 • Ahora calculamos los coeficientes.

p 
 a0 p ∞
* p nπ p nπ '
∫− p f ( x)dx = 2 ∫− pdx + ∑ ( an ∫− p cos
n =1 ) p
x dx + bn ∫ sin xdx % (3)

−p p &


Como cos(nπx/p) y sin(nπx/p) son ortogonales a 1
en este intervalo, entonces (3) se transforma en 


 p a0 p a0
p

 ∫− p f ( x)dx = ∫− p dx = x = pa0
2 2 −p

• Así tenemos
 1 p

 a0 = ∫ f ( x)dx
p −p
(4)
• Además, 


 p mπ
∫ − p f 
( x) cos p xdx
= a0 ∫ p cos mπ xdx + (5)

2 −p p

∞
* p mπ mπ p mπ nπ '

 ∑ ( an ∫− p cos p x cos p x dx + bn ∫− p cos p x sin p x dx %
n =1 ) &


por ortogonalidad tenemos

p mπ
∫− p cos p xdx = 0, m > 0
p mπ nπ
∫− p cos p x sin p xdx = 0
• y
 p mπ nπ # 0, m ≠ n

 ∫− p cos p x cos p xdx = "! p, m = n


Así (5) se reduce a


 p nπ

 ∫ f ( x ) cos xdx = an p
−p p
y por tanto 


 1 p nπ
an = ∫− p f ( x) cos x dx (6)
p p
• Finalmente, si multiplicamos (2) por sin(mπx/
p) y usamos


 p mπ
∫ sin xdx = 0, m > 0

 −p p

 p mπ nπ
∫ sin x sin xdx = 0

 −p p p
y 
 p mπ nπ # 0, m ≠ n

 ∫− p sin p x sin p xdx = "! p, m = n


obtenemos que



= % 0 , m ≠ n (7)
p sin m π x sin nπ x dx
∫− p p p
$
# p, m = n
 DEFINICIÓN 12.5
Series de Fourier
La serie de Fourier de una función f definida en el
intervalo (−p, p) se determina mediante



a0 ∞ nπ nπ
f ( x) = + ∑ (an cos
x + bn sin
x) (8)

2 n=1 p p
donde 

1 p
a 0 = ∫ f ( x ) dx (9)


 p −p

a = f ( x ) cos n π x dx
1 p
(10)



n
p ∫−p p
1 p n π (11)

bn = ∫ f ( x) sin x dx
p −p p
 Ejemplo 1
#0, −π < x < 0
Desarrolle f ( x) = " (12)

!π − x, 0≤ x <π
en una serie de Fourier.
Solución

La gráfica de f se muestra en la Fig 12.1 con p = π.
1 π
a0 = ∫ f ( x) dx =
π −π
1
π
[∫ 0dx + ∫ (π − x) dx]
0
−π 0
π

2 π
1' x $ π
= %π x − " =
π& 2 #0 2
 Ejemplo 1 (2)
1 π
an = ∫ f ( x) cos nx dx
π −π
[
1 0
π −π
π
= ∫ 0dx + ∫ (π − x) cos nx dx
0
]
π
1 ( sin nx 1 π %
= &(π − x) + ∫ sin nx dx #
π' n 0 n 0 $
π
1 cos nx
=− ←cos nπ = (-1)n
nπ n 0

− cos nπ + 1 1 − (−1) n
= 2
=
nπ n 2π
 Ejemplo 1 (3)
De (11) tenemos


 1 π 1
bn = ∫ (π − x) sin nxdx =

 π 0 n


Po tanto


 π ∞
( 1 − (−1) n
1 %
f ( x) = + ∑ ' 2
cos nx + sin nx $ (13)

4 n=1 & n π n #
Fig 12.1
 TEOREMA 12.1
Condiciones de Convergencia

Sean f y f’ continuas por partes en el intervalo (−p, p); esto es,

sean f y f’ continuas excepto en un número finito de puntos en
el intervalo y discontinuidades finitas sólo en estos puntos.
Entonces al serie de Fourier de f en el intervalo converge a f(x)
en un punto de continuidad. En un punto de discontinuidad, la
serie de Fourier converge al promedio


 f ( x + ) + f ( x −)

 2
donde f(x+) y f(x－) denotan el limite de f en x por la derecha y
por la izquierda, respectivamente.
 Ejemplo 2
• La función f en el Ejemplo 1, es continua en
(−π, π) excepto en x = 0. Así que la serie (13)
converge a 


 f (0 + ) + f (0 −) π + 0 π

 = =
2 2 2
en x = 0.
 Extensión Periódica
• Fig 12.2 es la extensión periódica de la función
f del Ejemplo 1. Así que la discontinuidad en x
= 0, ±2π, ±4π, … converge a


 f (0 + ) + f (0 −) π
=

 2 2
y en x = ±π, ±3π, … converge a

f (π +) + f (0−)
=0
2
Fig 12.2

Ch12_33
 Secuencia de Sumas Parciales
• Secuencia de Sumas Parciales

Para (13), escribimos las sums parciales como


 π π 2

 S1 = , S 2 = + cos x + sin x,
4 4 π

 π 2 1

 S3 = + cos x + sin x + sin 2 x
4 π 2


Fig 12.3.
 Fig 12.3


















 12.3 Series de Fourier de Coseno y Seno

• Funciones Pares e Impares

– par si f(−x) = f(x) 


– impar si f(−x) = −f(x) 


Fig 12.4 Función Par
Fig 12.5 Función Impar
 TEOREMA 12.2 Propiedades de Funciones Pares
e Impares

(a) El producto de dos funciones pares es par.

(b) El producto de dos funciones impares es impar.
(c) El producto de una función par y uan función impar
es impar.
(d) La suma (diferencia) de dos funciones pares es par.
(e) La suma (diferencia) de dos funciones impares es
impar. a a
∫ f ( x)dx = 2 ∫ f ( x)dx
(f) Si f es par, entonces −a 0
a
(g) Si f es impar, entonces ∫−a f ( x)dx = 0
 Series de Cosenos y Senos
• Si f es par en (−p, p) entonces


 1 p 2 p

 a0 = p ∫− p f ( x)dx = p ∫0 f ( x)dx


1 p nπ 2 p nπ

 an = ∫ f ( x) cos xdx = ∫ f ( x) cos xdx

 p −p p p 0 p
1 p nπ
bn = ∫ f ( x) sin xdx = 0
p −p p
• De manera similar, si f es impar en (−p, p) entonces
2 p nπ
an = 0, n = 0,1,2,... bn = ∫ f ( x) sin xdx
p 0 p
 DEFINICIÓN 12.6
Series de Fourier de Cosenos y Senos

(i) La serie de Fourier de una función par f en el


intervalo (−p, p) es la serie de cosenos



a0 ∞ nπ
f ( x) = + ∑ an cos
x (1)

2 n=1 p
donde


 2 p
a0 = ∫ f ( x) dx
p 0 (2)


 2 p nπ
a n = ∫ f ( x) cos x dx (3)
p 0 p
 (continuación)
DEFINICIÓN 12.6
Series de Fourier de Cosenos y Senos

(ii) La serie de Fourier de una función impar f en el


intervalo (−p, p) es la serie de senos 


 ∞
nπ
f ( x) = ∑ bn sin x (4)

i =1 p
donde


 2 p nπ
bn = ∫
p 0 f ( x ) sin x
p dx (5)
 Ejemplo 1
Desarrolle f(x) = x, −2 < x < 2 en una serie de Fourier.
Solución

Estudio de la Fig 12.6, muestra que es una función
par en (−2, 2) y p = 2. 


 2 nπ 4(−1) n+1

 bn = ∫ x sin x dx =
Así
 0 2 nπ


4 ∞
(−1) n+1 nπ (6)

f ( x) = ∑ sin x
πn =1 n 2
Fig 12.7 es la extensión periódica de la función del
Ejemplo 1.
Fig 12.6
Fig 12.7
 Ejemplo 2
• La función
 #− 1, − π < x < 0
f ( x) ="

 ! 1, 0≤ x <π


representada en la Fig 12.8 es impar en (−π, π) con p
= π.

De (5), 


 2 π 2 1 − (−1) n


 bn = ∫ (1) sin nx dx =
π 0 π n


y por tanto 


 2 1 − (−1)
∞ n

f ( x ) = π ∑ n sin nx (7)
n =1
Fig 12.8
 Fenómeno de Gibbs
• Fig 12.9 muestra las sumas parciales de (7). 

Podemos ver qeu la gráfica tiene picos
pronunciados cerca de las discontinuidades.
Este “exceso” SN no se alisa sino que
permanece constante aun cuando el valor de N
sea grande. Este comportamiento se conoce
como el fenómeno de Gibbs.
Fig 12.9
 Desarrollos en Semiintervalos
• Si una función f está definida sólo para 0 < x < L,
podemos suministrar una función arbitraria para −L < x
< 0.
• Si y = f(x) está definida para 0 < x < L,
– Reflejar al gráfica respecto al eje y en −L < x < 0; la
función hora es par. Fig 12.10.
– Reflejar la gráfica por el origen sobre −L < x < 0; la
función ahora es impar. Fig 12.11.
– definir f en −L < x < 0 mediante f(x) = f(x + L). Fig
12.12.
Fig 12.10
Fig 12.11
Fig 12.12
 Ejemplo 3
Desarrolle f(x) = x2, 0 < x < L, (a) en una serie de
cosenos, (b) en una serie de senos (c) en una serie de
Fourier.
Solución 

La gráfica está representada en la Fig 12.13.

 Ejemplo 3 (2)
(a) 
 2 L 2 2 2
a0 = ∫ x dx = L ,
L 0 3

 2 n

 2 L 2 nπ 4 L (−1)
an = ∫ x cos x dx = 2 2

 L 0 L nπ


Entonces 


 2 2 ∞ n
L 4L (−1) nπ
f ( x) = + 2 ∑ n 2 cos L x (8)
3 π n =1
 Ejemplo 3 (3)
(b)
 2 n +1 2

 bn = 2 ∫ x 2 sin nπ x dx = 2 L (−1) 4L
L
+ 3 3 [(−1) n − 1]
L 0 L nπ nπ




De ahí que


 2 ∞
2L ( (−1) n +1
2 % nπ
n
f ( x) = π ∑ ' n + n 3π 2 [( −1) − 1]$ sin L x (9)
n =1 & #
 Ejemplo 3 (4)
(c) Con p = L/2, nπ/p = 2nπ/L, tenemos 


 2 L 2 2 2 2 L 2 2 nπ L2
a
 0 = ∫ x dx = L , an = ∫ x cos xdx = 2 2

 L 0 3 L 0 L n π
2
2 2 n L
b
 n = ∫ x 2 sin
L π −
xdx =

 L 0 L nπ
Por tanto


 2
L L 2 ∞
( 1 2nπ 1 2nπ %
f ( x) = + ∑ ' 2 cos x − sin x$
3 π n= 1 & n π L n L # (10)

La gráfica de esta extensión se muestra en la Fig
12.14.
Fig 12.14
 Fuerza Impulsora Periódica
• Considere el siguiente sistema físico


 2
m x + k x = f (t )
d (11)

2

 dt
donde 


 ∞
x p (t ) = ∑ Bn sin
nπ
t
n=1 p (12)



es un desarrollo en serie de senos en un
semiintervalo.
 Ejemplo 4
• Recurriendo a (11), m = 1/16 de slug, k = 4 lb/pie, la
fuerza f(t) con período 2 se muestra en la Fig 12.15.
Aunque f(t) actúa en el sistema para t > 0, podemos
ampliar al gráfica con período 2 al eje t negativo para
obtener una función impar. Con p = 1, de (5) obtenemos




 1 2(−1) n+1
bn = 2 ∫ πt sin nπtdt =
De (11) obtenemos
0 n


2 ∞ n+1 (13)

1d x 2(−1)
2
+ 4x = ∑ sin nπ t
16 dt n =1 n
 Ejemplo 4 (1)
Para hallar la solución particular xp(t), sustituimos
(12) en (13). Así 


 1 2 2 2(−1) n +1


 (− n π + 4) Bn =
16 n

 n +1

 32(−1)
Bn = 2 2

 n(64 − n π )
Por tanto


 ∞
32(−1) n+1
x p (t ) = ∑ 2 2
sin nπt (14)
n =1 n(64 − n π )
 12.4 Series de Fourier Complejas
• Formula de Euler 

eix = cos x + i sin x

e-ix = cos x − i sin x (1)
 Series de Fourier Complejas
• De (1), tenemos 


 eix + e −ix eix − e −ix
cos x = 2 , sin
x=
2i (2)



Usando (2) para remplazar cos(nπx/p) y sin(nπx/p), se
tiene



a
 0 ∞ ) einπ x/p + e −inπ x/p einπ x/p + e −inπ x/p &
+ ∑ 'an + bn $

2 n=1 ( 2 2 %

 a ∞
1 1

= + ∑ (an − ibn )e
0 ) inπ x/p
+ ( a + ib ) e −inπ x/p &
n n

 2 n=1 '( 2 2 $%
∞
inπ x/p
∞

−inπ x/p (3)
= c0 + ∑ cne + ∑ c−ne
n =1 n =1
donde c0 = a0/2, cn = (an − ibn)/2, c-n = (an + ibn)/2.
Donde la función f es real, cn y c-n son números
complejos conjugados.

Tenemos


 c0 =
1 1 p
f ( x) dx
2 p
∫−p
(4)
 1
cn = (an − ibn )

 2

 = 1 ., 1 p f ( x) cos nπ x dx − i 1 p f ( x) sin nπ x dx +)

 2 - p ∫−p p p ∫−p p *

 1 p ( nπ nπ %
= ∫ f ( x)&cos x − i sin x # dx

 2 p −p ' p p $

 1 p
= ∫ f ( x)e −inπ x / p dx

 2 −p
(5)
 1
c−n = (an + ibn )

 2

 1. 1 p
= , ∫ f ( x) cos
nπ 1 p
x dx + i ∫ f ( x) sin
nπ +
x dx )

 2- p −p p p −p p *

 1 p ( nπ nπ %
= ∫ f ( x) &cos x + i sin # dx

 2 p −p ' p p$

 1 p inπ x / p
= ∫ f ( x )e dx

 2 −p
(6)
 DEFINICIÓN 12.7
Series de Fourier Complejas
Las Series de Fourier Complejas de función f definida
en un intervalo (−p, p) están dadas por


 ∞
inπ x / p
f ( x ) = ∑ cne (7)


 n = −∞

donde 
 1 p −inπ x / p
nc = f
∫− p ( x ) e dx , n = 0 , ± 1, ± 2, … (8)
2p
• Si f satisface la hopótesis del Teorema 12.1,
una serie d eFourier compleja converge a f(x)
en un punto de continuidad y al promedio 


 f ( x + ) + f ( x −)

 2


en un punto de discontinuidad.


 Ejemplo 1
Desarrolle f(x) = e-x, −π < x <π, en una serie de
Fourier compleja.
Solución

con p = π, (8) se obtiene

1 π − x −inx 1 π −( in +1) x
cn = ∫−πe e dx = ∫−πe dx
2π 2π
1 −( in +1)π ( in +1)π
= [e −e ]
2π (in + 1)
 Ejemplo 1 (2)
Empleando la fórmula de Euler


 e −( in +1)π −π n −π
= e (cos nπ − i sin nπ ) = (−1) e

 ( in +1)π π n π
e = e (cos nπ + i sin nπ ) = (−1) e


De ahí se tiene que 




 n ( eπ
− e −π
) n sinh π 1 − in
cn = (− 1) 2(in + 1)π = (−1 ) π n 2 + 1 (9)
 Ejemplo 1 (3)
Entonces la serie de Fourier compleja es


 sinh π ∞ n 1 − in inx
f ( x) = ∑ (−1) 2 e
π n=−∞ n + 1 (10)



La serie (10) converge al desarrollo de período
2π de f.
 Frecuencia Fundamental
• El período fundamental es T = 2p y por tanto p
= T/2. 

La serie de Fourier se transforma en 


 a0 ∞
+ ∑ (an cos nωx + bn sin nωx)

 2 n=1
∞ (11)

inω x

 ∑n c e

 n = −∞

donde ω = 2π/T se llama frecuencai angular

fundamental.
 Espectro de Frecuencias
• Si f es periódica y tiene período fundamental T,
el conjunto de puntos (nω, |cn|) se llama
espectro de frecuencias de f.
 Ejemplo 2
• En el Ejemplo 1, ω = 1, por lo cual nω ecibe
valores de 0, ±1, ±2, … 

Usando , vemos de (9) que

α + i β = α 2
+ β 2


 sinh π 1
| cn | =

 π 2
n +1
Fig 12.17.
 Fig 12.17
 Ejemplo 3
• Halle el espectro de la onda mostrada en Fig12.18. La
onda es la extensión periódica de la función f:
% 0, − 12 < x < − 14
"
f ( x) = $1, − 14 < x < 14
" 0, 1 < x < 12
# 4
 Ejemplo 3 (2)
Solución

Aquí T = 1 = 2p so p = ½. Como f es 0 en (−½, −¼) y (¼,
½), (8) se transforma en
12 2inπx 14 2inπx
cn = ∫ f ( x )e dx = ∫ (1)e dx
−1 2 −1 4

e 2inπx 1 4 1 einπ / 2 − e −inπ / 2

= =
2inπ − 1 4 nπ 2i
1 nπ
cn = sin
nπ 2
 Ejemplo 3 (3)
Es fácil de comprobar que


 14 1

 c0 = ∫ dx =
−1 4 2


Fig 12.19 ilustra el espectro de frecuencias de f.
Fig 12.19
 12.5 Problema de Sturm-Liouville
• Valores propios y funciones propias

Recuerde el ejemplo Ejemplo 2, Sec 3.9



y!! + λy = 0 , y (0) = 0 , y ( L) = 0 (1)


Esta ecuación posee soluciones no triviales sólo

cuando λ toma valores λn = n2π2/L2, n = 1, 2, 3,…
llamados valores propios. Las soluciones no triviales
correspondientes y = c2 sin(nπx/L) o simplemente y =
sin(nπx/L) se llaman funciones propias.
 Ejemplo 1
• Se deja como ejercicio demostrar que los tres
casos posibles: λ = 0, λ = −α2 < 0, λ = α2 > 0, (α >
0), que los valores propios y las funciones
propias para



y!! + λy = 0 , y !(0) = 0 , y!( L) = 0 (2)



son respectivamente λn = αn2 = n2π2/L2, n = 0, 1,
2, …y y = c1 cos(nπx/L), c1 ≠ 0.
 Problema Regular de Sturm-Liouville
• Sean p, q, r y rʹ funciones de valores reales
continuas en [a, b], y sea r(x) > 0 y p(x) > 0 para todo
x en el intervalo. Entonces se dice que


 d
Resolver dx [r ( x) y !] + (q ( x) + λ p ( x)) y = 0 (3)



Sujeta a 
 α1 y (a ) + β1 y!(a ) = 0
(4)

α 2 y (a ) + β 2 y !(a ) = 0 (5)

es un problema regular de Sturm-Liouville. Los
coeficientes en (4), (5) se suponen reales e
independientes de λ.
 TEOREMA12.3 Propiedades del Problema Regular
de Strum-Liouville
(a) Existe un número infinito de valores propios reales

que se pueden arreglar en orden creciente
λ1 < λ2 < λ3 < … < λn < … tal que λn → ∞ cuando n → ∞.
(b) Para cada vlor propio hay sólo uan función propia 

(excepto para multiplos constantes no nulos).
(c) Las funciones propias que corresponden a diferentes
valores propios son linealmente independientes.
(d) El conunto d efunciones propias que corresponden
al conjunto de valores propios es ortogonal con
respecto a la función pesop(x) en el intervalo [a, b].
Ch12_83
Demostración de(d)

Sean ym e yn be funciones propias correspondientes a
valores propios λm y λn. Entonces


 d
[r ( x) ym! ] + (q ( x) + λm p ( x)) ym = 0
dx (6)


 d
[r ( x) yn! ] + (q ( x) + λn p ( x)) yn = 0
dx (7)



De (6)⋅yn − (7)⋅ym tenemos



d d
(λm − λn ) p ( x) ym yn = ym [r ( x) yn ']− yn [r ( x) ym ']
dx dx
Integrando la ecuación anterior de a a b, se tiene


 b
(λm − λn ) ∫ p ( x) ym yn dx
a (8)


 = r (b)[ ym (b) yn! (b) − yn (b) ym! (b)]

 − r (a )[ ym (a ) yn! (a ) − yn (a ) ym! (a )]

Como todas las soluciones deben satisfacer las

condiciones de frontera (4) y (5), de (4) tenemos

A1 ym (a ) + B1 ym ' (a ) = 0
A1 yn (a ) + B1 yn ' (a ) = 0
Para que A1 y B1 no nulas ambas, satisfagan el sistema,
el determinante de los coeficientes debe valer cero




 ym ( a ) yn ' ( a ) − yn ( a ) ym ' ( a ) = 0
De manera similar de (5), tenemos




 ym (b) yn ' (b) − yn (b) ym ' (b) = 0
Así el lado derecho de (8) vale cero.

De ahí tenemos la relación de ortogonalidad



b p ( x) y ( x) y ( x) dx =
0, λ ≠ λ (9)
∫a m n m n
 Ejemplo 2
Resolver

y! ! + λy =
0, y ( 0) = 0, y (1) + y!(1) = 0 (10)
Solución 

Se debería verificar que para λ = 0 y λ < 0, (10) sólo
posee la solución trivial. Para λ = α2 > 0, α > 0, la
solución general es y = c1 cos αx + c2 sin αx. Ahora la
condición y(0) = 0 implica c1 = 0, así que y = c2 sin αx. La
segunda condición y(1) + yʹ(1) = 0 implica c2 sin α + c2α
cosα. = 0.
 Ejemplo 2 (2)
Escogiendo c2 ≠ 0, tenemos

tan
α = −α (11)

De Fig 12.20, vemos que hay infinitas soluciones para α
> 0. Es fácil obtener los valores de α > 0. Así que los
valores propios son λn = αn2, n = 1, 2, 3, …

y las funciones propias correspondientes son

yn = sin αnx. 

 Fig 12.20








 Problema Singular de Sturm-Liouville
• Existen varias condiciones para (3)
– r(a) = y se especifica una condición de frontera del
tipo provisto en (5) en x = b; (12)
– r(b) = 0 y se especifica una condición de frontera
del tipo provisto en (4) en x = a. (13)
– r(a) = r(b) = 0 y no se especifica ninguna condición
de frontera en x = a ni en x = b; (14)
– r(a) = r(b) y las condiciones de frontera y(a) = y(b),
y’(a) = y’(b). (15)
 Observaciones:
• La ecuación (3) que satisface (12) y (13) es un
problema singular de valores en la frontera.

La ecuación (3) que satisface (15) es un problema
periódico de valores en la frontera.



• Al suponer que las soluciones de (3) están
acotadas en [a, b], de (8) se tiene
– Si r(a) = 0, entonces la relación de ortogonalidad
(9) se cumple sin ninguna condición en la frontera
en x = a; (16)
– Si r(b) = 0 , entonces la relación de ortogonalidad
(9) se cumple sin ninguna condición en la frontera
en x = b; (17)
– Si r(a) = r(b) = 0, entonces la relación de
ortogonalidad (9) se cumple sin ninguna condición
en la frontera en x = a ni en x = b; (18)
– Si r(a) = r(b), entonces la relación de ortogonalidad
(9) se cumple con las condiciones de frontera
periódicas y(a) = y(b), y’(a) = y’(b). (19)
 Forma Autoconjunta
• En realidad (3) es al misma que

r ( x) y!! + r !( x) y! + (q( x) + λp( x )) y = 0 (20)

Así podemos escribir la ecuación diferencial de
Legendre (1 ) y "−2 xy'+ n(n + 1 ) y = 0 como

− x 2



d 2
[(1 − ) y! ] + n(n + 1) y = 0
x (21)

dx

Aquí hallamos que el coeficiente de yʹ es al

derivada del coeficiente de yʺ.
• Además, si los coeficientes son continuos y a(x) ≠ 0 para
todo x en un intervalo, entonces cualquier eduación
diferencial de segundo orden 

a ( x) y!! + b( x) y! + (c( x) + λd ( x )) y = 0 (22)

swe puede reformular en la llamada forma autoadjunta
(3).
• Para entender el hecho anterior, empezamos desde

a1(x)yʹ + a0(x)y = 0

Sea P = a0/a1, µ = exp(∫ Pdx), µʹ = Pµ, entonces

yʹ + Py = 0, µyʹ + µPy = 0, 

Así d(µy)/dx = 0.

 • Ahora para (22), sea Y = yʹ, el factor de integración e∫
[b(x)/a(x)] dx. 


En este caso (22) se transforma en 




∫ b ( x ) / a ( x ) dx b( x) ∫ b ( x ) / a ( x ) dx d & ∫ b ( x ) / a ( x ) dx #
e 
 Y '+ e Y + ... = e Y + ...
a( x) dx $% !"


En resumen, (22) puede transformarse en
∫ ( b / a ) dx b( x) ∫ (b / a ) dx
e y"+ e y'
a( x)
& c( x) ∫ (b / a ) dx d ( x) ∫ (b / a ) dx #
+$ e +λ e !y = 0 (23)
% a( x) a( x) "
 • Además, (23) es la misma que (3)
d - ∫ (b / a ) dx * ' ∫ (b / a ) dx d ( x) ∫ (b / a ) dx $
+ e y '( + %% e +λ e "" y = 0
dx , ) & a( x) #
∫ ( b / a ) dx ∫ ( b / a ) dx d ( x) ∫ (b / a ) dx
donde r ( x) = e , q( x) = e , p( x) = λ e
a( x)
 Ejemplo 3
• En la Sec 5.3, vimos que la solución general de al
ecuación diferencial paramétrica de Bessel 


 x 2y"+ xy '+(α 2 x 2 − n 2 ) y = 0, n = 0, 1, 2, ...

 es y = c1 J n (αx) + c2Yn (αx)


Dividiendo la ecuación de Bessel entre x2 y
multiplicando la ecuación resultante por el factor de
integración e∫ [(1/x)] dx = eln x = x, tenemos
2 2
n d n
xy"+ y '+(α 2 x − ) y = 0, or [xy']+ (α 2 x − ) y = 0
x dx x
donde r = x, q = −n 2 / x, p = x, λ = α 2
 Ejemplo 3 (2)
Ahora r(0) = 0, y de las dos soluciones Jn(αx) y Yn(αx)
sólo Jn(αx) está acotada en x = 0. De (16), el conjunto
{Jn(αix)}, i = 1, 2, 3, …, es ortogonal con respecto a la
función peso p(x) = x en [0, b]. 

Así 
 b
∫0 xJ n (λi x) J n (λ j x) dx = 0, λi ≠ λ j , (24)


Siempre queαi y por consiguiente los valores propios λi

= αi2 se definen por medio de un acondición límite en x
= b del tipo provisto en (5):

A2Jn(λb) + B2αJnʹ(αb) = 0 (25)
 Ejemplo 4
• De (21), identificamos q(x) = 0, p(x) = 1 y λ = n(n + 1).
Recuerde de la Sec 5.3 que cuando n = 0, 1, 2, …, la
ED de Legendre posee soluciones polinomiales Pn(x).
Observamos que r(−1) = r(1) = 0 junto con el hecho de
que Pn(x) son las únicas soluciones de (21) que están
acotadas en [−1, 1], para concluir que el conjunto
{Pn(x)}, n = 0, 1, 2, …, es ortogonal con respecto a la
función peso p(x) = 1 en [−1, 1]. Así
1
∫ Pm ( x) Pn ( x)dx = 0, m ≠ n
-1
 12.6 Series de Bessel y Legendre
• Series de Fourier-Bessel 

Hemos demostrado que{Jn(αix)}, i = 1, 2, 3, …es
ortogonal con respecto a p(x) = x en [0, b] cuando αi
esté definida por medio de
 !
A2 J n (αb) + B2α J n (αb) = 0
(1)

Esta serie ortogonal, o serie de Fourierde generalizada,
el desarrollo de una función f definida en (0, b) en
términos de este conjunto ortogonal es

∞
f ( x) = ∑ ci J n (α i x) (2)

i =1
b
donde
 ∫ x J n (α i x) f ( x) dx
0
ci = 2 (3)
J n (α i x)
 La norma cuadrada Jn(αix) se define mediante


 2 b 2
J n (α i x ) =
0
∫ xJ n (α i x ) dx
(4)



Esta serie (2) se llama series de Fourier-Bessel.
 Relaciones de Recurrencia Diferenciales

• Recordando (20) y (21) da la Sec 5.3, tenemos

las relaciones de recurrencia diferenciales
como


 d n n
dx [ x J n ( x)]
= x J n−1 ( x) (5)


 d −n −n
[ x
dx n J ( x )] = − x J n+1 ( x), (6)
 Norma Cuadrada
• El valor de (4) depende de λi = αi2. Si y = Jn(αx)
tenemos que


 d ' 2 n2 $

 [xy']+ %% α x − "" y = 0
dx & x#


Al multiplicar por 2xy’, se tiene
d 2 d
[xy'] + (α x − n ) [ y ] = 0
2 2 2 2
dx dx
• Integrandopor partes [0, b], se obtiene


 2 b 2 b

 (
2 2 2 2
2α ∫ xy dx = [xy '] + (α x − n ) y
0
2
) 0


Como y = Jn(αx), el límite inferior es 0 para n > 0, porque
Jn(0) = 0. Para n = 0, 

en x = 0. Así 


 2 b 2
2
∫0 n (αx ) dx
α xJ (7)


 2 2 2 2 2 2 2
= α b [ J n (αb)] + (α b − n )[ J n (αb)] ,
#
donde yʹ = αJnʹ(αx).
• Ahora se consideran tres casos de (1).
– Caso I: Si se elige A2 = 1 y B2 = 0, entonces (1) es

J n (α
b) = 0 (8)

Hay un número infinito de raíces positivas xi = αib
de (8) (see Fig 5.3), que definen αi = xi/b. 

Los valores propios son positivos y λi = αi2 = (xi/b)2. 

De las raíces negativas de (8) no resulta ningún
nuevo valor propio puesto que Jn(−x) = (−1)nJn(x).
 El número 0 no es un valor propio de para ningún
n puesto que Jn(0) = 0, n= 1, 2, 3, … y J0(0) = 1. 

Cuando (6) se escribe como xJnʹ(x) = nJn(x) –
xJn+1(x), de (7) y (8) se deduce

2

 2 b 2
|| J nα i ( x) || = J n+1 (α ib).
2 (9)


– Caso II: Si se elige A2 = h ≥ 0 y B2 = b, entonces (1)
es

hJ n ( α b ) + α b J ! (αb) = 0.

n (10)

Hya un número infinito de raíces positvas xi = αib
para n = 1, 2, 3, …. Como antes, λi = αi2 = (xi/b)2. 

λ = 0 no es un valor propio para n = 1, 2, 3, ….
Sustituyendo αibJnʹ(αib) = – hJn(αib) en (7), se tiene


 2 2 2 2
2 αi b − n + h 2
|| J (
n i
α x ) || = 2 n (α ib). (11)
J
2α i
– Caso III: Si h = 0 y n = 0 en (10), αi se definen da
las raíces 

J 0! (α b) = 0
(12)

Aunque (12) es sólo un caso especial de (10), es la
única solución para la cual λ = 0 es un valor propio.
Para n = 0, el resultado en (6) implica que J0ʹ
(αb) = 0 es equivalente a J1(αb) = 0.
 Como x1 = α1b = 0 es una raíz de la última ecuación
y puesto que J0(0) = 1 no es trivial, deducimos de
λ1 = α12 = (x1/b)2 que λ1 es un valor propio. Pero
no podemso utilizar (11) cuando α1 = 0, h = 0, n
= 0, y n = 0. Sin embargo de (4) tenemos


 2 b b 2

|| 1 || = ∫ x dx = (13)

0 2
Para αi > 0 podemos usar (11) con h = 0 y n = 0:


 b 2
2 2
|| J 0 (α i x) || = J 0 (α ib ). (14)
2
 DEFINICIÓN 12.8
Serie de Fourier-Bessel

La serie de Fourier-Bessel de una función f definida en

el intervalo (0, b) se expresa mediante
(i) ∞
f ( x) = ∑ ci J n (α i x)
i =1 (15)
2 b
ci = 2 2 ∫ x J n (α i x) f ( x) dx
b J n+1 ( α ib) 0 (16)
donde αi se definen mediante Jn(αb) = 0.
 (continuación)
DEFINICIÓN 12.8
Serie de Fourier-Bessel

(ii) ∞

f ( x) = ∑ c i J n (α i x) (17)
i =1

2α i2 b
ci = 2 2 2 2 2 ∫ x J n (α i x) f ( x) dx (18)
(α i b − n + h ) J n (α ib) 0

donde αi se definen mediante hJn(αb) + αbJ’n(αb) = 0.

 (continuación)
DEFINICIÓN 12.8
Serie de Fourier-Bessel

(iii) ∞
f ( x) = c1 + ∑ ci J 0 (α i x) (19)
i =2

2 b 2 b
c 1 = 2 ∫ x f ( x ) dx, c i = 2 2 ∫ x J 0 ( α i x) f ( x) dx(20)
b 0 b J 0 (α ib) 0

donde the αi se definen mediante J’0(αb) = 0.

 TEOREMA 12.4
Condiciones para Convergencia

Si f y f’ son continuas por partes en el intervalo abierto

(0, b), entonces un desarrollo de Fourier-Bessel de f
converge a f(x) en algún punto donde f es continua y
al promedio [f(x+) + f(x-)] / 2 en algún punto donde
f es discontinua.
 Ejemplo 1
Desarrolle f(x) = x, 0 < x < 3, enn una serie de
Fourier-Bessel , usando función de Bessel de
orden uno que satisfacen la condición límite
J1(3α) = 0.
Solución

Empleamos (15) donde ci se expresan mediante
(16) con b = 3:
 2 3 2
ci = 2 2 ∫ x J1 (α i x)dx
3 J2 (3α ) 0
i
 Ejemplo 1 (2)
Sea t = αi x, dx = dt/αi, x2 = t2/αi2, y use (5) en la
forma d[t2J2(t)]/dt = t2J1(t):

2 3α d
2 i 2
ci = 3 2 ∫ [t J 2 (t )]dt =
9α i J 2 (3α )
i
0 dt α i J 2 (3α i )
Por tanto el desarrollo es
∞
2
f ( x ) = 2∑ J1 (α i x)
i =1 α i J 2 (3α i )
 Ejemplo 2
• Si las αi del Ejemplo 1 se definen mediante

J1(3α) + αJ1ʹ(3α) = 0, lo único que cambia en el
desarrollo es el valor de la norma cuadrada. Como
3J1(3α) + 3αJ1ʹ(3α) = 0 que concuerda con (10)
cuando h = 3, b = 3 y n = 1. Así (18) y (17) producen a
su vez c = 18α i J 2 (3α i )
i 2 2
(9α i + 8) J1 (3α i )
∞
α i J 2 (3α i )
f ( x) = 18∑ 2 2
J1 (α i x)
i =1 (9α i + 8) J1 (3α i )
 DEFINICIÓN 12.9
Serie de Fourier-Legendre

La serie de Fourier-Legendre de una función f definida

en el intervalo (－1, 1) se expresa mediante
(i) ∞

f ( x) = ∑ cn Pn ( x) (21)
n =0
2n + 1 1
cn = ∫ f ( x) Pn ( x) dx
2 −1 (22)
donde αi se definen mediante Jn(αb) = 0.
 TEOREMA 12.5
Condiciones de Convergencia

Si f y f’ son continuas por partes en el intervalo abierto

(－1, 1), entonces un desarrollo en serie de
Fourier-Legendre (21) converge a f(x) en algún punto
donde f es continua y al promedio [f(x+) + f(x-) / 2 en un
punto donde f es discontinua.
 Ejemplo 3
Escriba los cuatro primeros términos no nulos del
desarrollo de Fourier-Legendre de


 %0, − 1 < x < 0
f ( x) = $
#1, 0 ≤ x < 1
Solución 

De la página 269 y (22):

1 1 1 1 1
c0 = ∫ f ( x) P0 ( x)dx = ∫ 1 ⋅ 1dx =
2 −1 2 0 2
3 1 3 1 3
c1 = ∫ f ( x) P1 ( x)dx = ∫ 1 ⋅ xdx =
2 −1 2 0 4
 Ejemplo 3 (2)
5 1 5 1 1 2
c2 = ∫ f ( x) P2 ( x)dx = ∫ 1 ⋅ (3 x − 1)dx = 0
2 −1 2 0 2
7 1 7 1 1 3 7
c3 = ∫ f ( x) P3 ( x)dx = ∫ 1⋅ (5 x − 3x)dx = −
2 −1 2 0 2 16
9 1 9 1 1
c4 = ∫ f ( x) P4 ( x)dx = ∫ 1⋅ (35 x 4 − 30 x 2 + 3)dx = 0
2 −1 2 0 8
11 1 11 1 1 5 3 11
c5 = ∫ f ( x) P5 ( x)dx = ∫ 1⋅ (63x − 70 x + 15 x)dx =
2 −1 2 0 8 32
De ahí que
1 3 7 11
f ( x) = P0 ( x) + P1 ( x) − P3 ( x) + P5 ( x) + ...
2 4 16 32
Véase Fig 12.22.
 Fig 12.22












 Otra Forma de la Serie
• Si se establece x = cos θ, x = 1 implica que θ = 0, x
= −1 implica que θ = π. Como dx = −sin θ dθ, (21) y
(22) se transforma, respectivamente, en


 ∞
F (θ ) = ∑ c n Pn (cos
θ) (23)

n =0


2n + 1 π
cn = ∫ 0 F (θ ) Pn (cos θ ) sin θ dθ , (24)

2


donde f(cos θ) se ha remplazado por F(θ).