Contreras Roque Adal Eduardo 201935255

SERIES DE FOURIER DE COSENOS Y DE SENOS

Introducción

El esfuerzo que se invierte en la evaluación de las integrales definidas que calculas los coeficientes
𝑎0 , 𝑎𝑛 , 𝑏𝑛 al desarrollar una función f en una serie de Fourier se reduce significativamente f f es
una función par o impar. Recuerde que se dice que una función f es

Par si f(-x) =f(x) e impar si f(-x)= -f(x)

En un intervalo simétrico tal como (-p,p), la gráfica de una función par tiene simetría respecto el
eje y, mientras que la de una función impar tiene simetría respecto al origen. FUNCIONES PAR E
IMPAR Es muy probable que el origen de los terminos par e impar sea consecuencia del hecho de
que las graficas de funciones polinomiales de potencias pares de x son simétricas respecto al eje
y, mientras que las graficas de de potencias impares de x son simétricas respecto al origen. Por
ejemplo,

F(x)= 𝑥 2 es par, ya que f(-x) = (𝑥 2 )= 𝑥 2 = f(x)

F(x) =𝑥 3 es impar, ya que f(-x) = (−𝑥)3 = -𝑥 3 = -f(x)

PROPIEDADES

El teorema siguiente lista algunas propiedades de las funciones

pares e impares.

TEOREMA 11.3.1 Propiedades de funciones pares/ impares

Contreras Roque Adal Eduardo 201935255

a) El producto de dos funciones pares es par.

b) El producto de dos funciones impares es par.

c) El producto de una función impar y una función par es impar.

d) La suma (diferencia) de dos funciones pares es par.

e) La suma (diferencia) de dos funciones impares es impar.

𝑎 𝑎
f) Si f es par, entonces ∫−𝑎 𝑓(𝑥)𝑑𝑥 =2 ∫0 𝑓(𝑥)𝑑𝑥
𝑎
g) Si f es impar, entonces a a ∫−𝑎 𝑓(𝑥)𝑑𝑥 =0

SERIES DE COSENOS Y DE SENOS

Si f es una función par en (-p,p), entonces, en vista de las propiedades anteriores, los coeficientes
(9),(10), (11) se convierten en

1 𝑝 𝑝
𝑎0 =𝑝 ∫−𝑝 𝑓(𝑥)𝑑𝑥=2/p ∫0 𝑓(𝑥)𝑑𝑥

1 𝑝 𝑛𝜋 𝑝 𝑛𝜋
𝑎𝑛 =𝑝 ∫−𝑝 𝑓(𝑥)𝑐𝑜𝑠 𝑝
𝑥𝑑𝑥 =2/p ∫0 𝑓(𝑥)𝑐𝑜𝑠 𝑝 𝑥𝑑𝑥

1 𝑝 𝑛𝜋
𝑏𝑛 =𝑝 ∫−𝑝 𝑓(𝑥)𝑠𝑒𝑛 𝑝
𝑥𝑑𝑥=0

De la misma manera, cuando f es impar en el intervalo (-p, p)

1 𝑝 𝑛𝜋
𝑎𝑛 = 0, 𝑛 = 0,1,2, 𝑏𝑛 =𝑝 ∫−𝑝 𝑓(𝑥)𝑠𝑒𝑛 𝑝
𝑥𝑑𝑥

Resumiremos los siguientes resultados en la siguiente definición

DEFINICIÓN 11.3.1 Series de Fourier de cosenos y de senos

Contreras Roque Adal Eduardo 201935255

i) La serie de Fourier de una función f par en el intervalo (-p, p) es la serie de cosenos

𝑎0 𝑛𝜋
F(x) =
2
+ ∑∞
𝑛=1 𝑎𝑛 𝑜𝑠 𝑝
𝑥

2 𝑝
𝑎0 =𝑝 ∫0 𝑓(𝑥) 𝑑𝑥

2 𝑝 𝜋𝑛
𝑎𝑛 = 𝑓(𝑥) cos
𝑝 ∫0 𝑝

ii) La serie de Fourier de una función f impar en el intervalo (-p, p) es la serie de senos
𝜋𝑛
F(x) = ∑∞
𝑛=1 𝑏𝑛 sen 𝑥,
𝑝

Donde

2 𝑝 𝜋𝑛
𝑏𝑛 = ∫0 𝑓(𝑥) sen 𝑥 𝑑𝑥,
𝑝 𝑝

𝜋𝑛
El término sen es 0 en x =-p, x= 0 y x= p, entonces la serie de senos (4) converge a 0 en esos
𝑝
puntos sin importar si f esta definida en estos puntos

EJEMPLO 1 Desarrollo en una serie de senos

Desarrolle f (x) =x, -2 < x < 2 en una serie de Fourier

solución EL examen d ela figura muestra que la función es impar en el intervalo (-2,2) así que
desarrollamos f en una serie de senos, Identificando 2p=4 tenemos p=2 Por lo que la ecuación
después de integrar

2 𝑝 𝜋𝑛 4(−1)𝑛+1
𝑏𝑛 = ∫0 𝑓(𝑥) sen 𝑥 𝑑𝑥 =
𝑝 𝑝 𝑛𝜋

4 (−1)𝑛+1 𝜋𝑛
Por tanto f(x)= ∑∞ sen x
𝜋 𝑛=1 𝑛 𝑝

La función del ejemplo 1 satisface las condiciones del teorema 11.2.1. Por tanto la serie (6)
converge a la función en el intervalo (
2, 2) y la extensión periódica
Contreras Roque Adal Eduardo 201935255

Ejemplo 2

1 Las funciones sin  y cos son funciones periódicas de periodo 2. Las funciones
constantes son funciones periódicas de cualquier periodo (en el sentido de la definición).
Ejemplo de funciones que no son periódicas son , 2,  y ln .

FENÓMENO DE GIBBS

En la figura con un SAX hemos trazado las graficas de 𝑆1 (𝑥), 𝑆3 (𝑥), 𝑆15 (𝑥) de las sumas parciales
de los términos distintos de cero de la expresión. Como se muestra en la figura tiene picos
notables cerca de las discontinuidades en x=0, x=π , x=π etcétera. Este “exceso ” de las sumas
parciales 𝑆𝑁 , respecto a los valores de la respecto a los valores de la función cerca de un punto de
discontinuidad no se empareja, sino que permanece bastante constante, aunque el valor de N sea
muy grande. A este comportamiento de una serie de Fourier cerca de un punto en el que f es
discontinua se le llama fenómeno de Gibbs.

La extensión periódica de f en el ejemplo 2, sobre todo el eje x , es una función serpenteante

DESARROLLOS EN SEMIINTERVALOS

En el análisis anterior hemos sobreentendido que una función f esta definida en un

intercalo con el origen en su punto medio, es decir(-p,p) Sin embargo, en muchos casos
nos interesa representar una función f que esta definida solo por 0<x<L con una serie
trigonométrica. Esto se puede hacer de muchas formas distintas dando una definición
arbitraria de f(x) para -L<x<0. Por brevedad consideraremos los tres casos mas
importantes. SI y =f(x) esta definida en el intervalo (0,L ) entonces

i) Reflejar la grafica de f respecto al eje en (-L); la función ahora es par en (-L,L )

ii) Reflejar la grafica de f respecto al origen (-L , 0 ); la función ahora es impara en
(-L,L )
iii) Definir f en (-L,0 ) con y=f(x+L)
Contreras Roque Adal Eduardo 201935255

EJEMPLO 3 Desarrollo en tres series

Desarrolle f(x) = 𝑥 2 , 0<x< L

a) En una serie de cosenos b) en una serie de senos c) en una serie de Fourier

Solución En la figura se presenta la grafica de esta función

a) Tenemos
2 𝑙 2 2 𝐿 4 𝐿2 (−1)𝑛
𝑎0 = ∫0 𝑋 2 dx= 𝐿3 , 𝑎𝑛 = ∫0 𝑋 2 xdx =
𝐿 3 𝐿 𝑛2 𝜋 2

donde hemos integrado por partes dos veces en la evaluación de 𝑎𝑛

𝐿2 4 𝐿2 ∞ (−1)𝑛 𝑛𝜋
Por tanto f(x) = 3
+ 𝜋2
∑𝑛=1 𝑛2 cos 𝐿
x

b) En este caso debemos nuevamente integrar por partes dos veces:

2 𝐿 𝑛𝜋 2𝐿2 (−1)𝑛+1 4 𝐿2
𝑏𝑛 = ∫0 𝑋 2 x sen dx = + [(−1)𝑛 − 1]
𝐿 𝐿 𝑛𝜋 𝑛3 𝜋 3

2𝐿2 (−1)𝑛+1 2 𝑛𝜋
Por tanto f(X) = [∑∞
𝑛=1 + ] (−1)𝑛 − 1] sen 𝑋
𝜋 𝑛 𝑛3 𝜋3 𝐿

c) Con p= L/2 , 1/p=2/L y nπ/p=2nπ/L ,tenemos

2 𝑙 2 2 𝐿 2𝑛𝜋 𝐿2
𝑎0 = ∫0 𝑋 2 dx= 𝐿2 , 𝑎𝑛 = ∫0 𝑋 2 cos xdx =
𝐿 3 𝐿 𝐿 𝑛2 𝜋 2

2 𝐿 𝑛𝜋 𝐿2
Y 𝑏𝑛 = 𝐿 ∫0 𝑋 2 x sen 𝐿
dx = − 𝑛𝜋

Por tanto

𝐿2 𝐿2 ∞ 1 2𝑛𝜋 1 2𝑛𝜋
f(x) = 3
+ ∑
𝜋 𝑛=1 𝑛2 𝜋
cos 𝐿
x- 𝑛sen 𝐿 𝑥]
Contreras Roque Adal Eduardo 201935255

Las series convergen hacia ;a extensión periódica par de periodo 2L de f, la extensión periódica
impar de periodo 2L de f y la extensión periódica de periodo L de f, respectivamente. En la figura
se presentan las graficas de esas extensiones periódicas.

FUERZA IMPULSORA PERIÓDICA Algunas veces las series de Fourier son útiles para determinar una
solución particular de la ecuación diferencial que describe un sistema físico en el que la entrada o
fuerza impulsora f (t) es periódica. En el siguiente ejemplo encontraremos una solución particular
de la ecuación diferencial.

𝑑2𝑥
m + kx =f(t)
𝑑𝑡 2

representando primero f por el desarrollo en serie de senos en un semiintervalo y después

suponiendo una solución particular de la forma

𝑛𝜋
𝑋𝑝 (𝑡) = ∑∞
𝑛=1 𝐵𝑛 sen 𝑡
𝑝

EJEMPLO 4 Solución particular de una ED

1
Un sistema resorte-masa no amortiguado en el que la masa es m = slug y la constante del
16
resorte es k 4 lb/pie, es impulsado por una fuerza externa f (t) de periodo 2 como se muestra en la
figura. Aunque la fuerza f(t) actúa sobre el mismo sistema cuando t>0, observe que si se extiendo
la grafica de la función hacia la parte negativa del eje t para que su periodo sea 2, obtenemos una
función impar. En términos prácticos esto significa que solo necesitamos encontrar el desarrollo
en una serie de senos en un semiintervalo de f(t)= π t, 0 < t<1 . Con p=1 utilizando la ecuación e Commented [ACR1]:
integrando por partes se tiene que

1 2(−1)𝑛+1
𝑏𝑛 = 2 ∫0 𝑡 π sen nπ t dt = 𝑛

De la ecuación (11) la ecuación diferencial de movimiento es

1 𝑑2𝑥 2 (−1)𝑛+1
16 𝑑𝑡 2
+ 4x = ∑∞
𝑛=1 𝑛
sen nt
Contreras Roque Adal Eduardo 201935255

Para encontrar una solución particula

𝑥𝑝 (𝑡) 𝑑𝑒 𝑙𝑎 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 , 𝑠𝑢𝑠𝑡𝑖𝑡𝑢𝑖𝑚𝑜𝑠 𝑒𝑛 𝑙𝑎 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 𝑒 𝑖𝑔𝑢𝑎𝑙𝑎𝑚𝑜𝑠 𝑙𝑜𝑠 𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒𝑠 𝑑𝑒 𝑠𝑒𝑛 𝑛𝜋t .
Así obtenemos

1 2 (−1)𝑛+1 32 (−1)𝑛+1
(- 𝑛2 𝜋 2 + 4)𝐵𝑛 = o 𝐵𝑛 =
16 𝑛 𝑛(64− 𝑛2 𝜋2 )

Por tanto

32 (−1)𝑛+1
𝑋𝑝 (𝑡) = ∑∞
𝑛=1 𝑛(64− 𝑛2 𝜋2 )
sen nπt

Observe que en la solución no hay entero n ≥1 para el cual el denominador de 𝐵𝑛 , que es 64 -

𝑛2 𝜋 2 , sea cero. En general, si existe un valor de n, digamos N, para el cual Nπ/p=w , donde w =
√𝑘/𝑚 entonces el estado del sistema que describe la ecuación ) es un estado de resonancia pura.
Es decir, tenemos resonancia pura si el desarrollo de la función f(t) de la fuerza impulsora en serie
de Fourier contiene un término sen (Nπ/L )t (o cos (Nπ/L)t ) que tenga la misma frecuencia que la
de las vibraciones libres.

Por supuesto, si la extensión de la fuerza impulsora f con periodo 2p sobre el eje negativo de t da
como resultado una función par, entonces desarrollamos f en una serie de cosenos.
Contreras Roque Adal Eduardo 201935255

11.4 PROBLEMA DE STURM-LIOUVILLE

INTRODUCCIÓN En esta sección estudiaremos algunos tipos especiales de problemas con valores
en la frontera en los que la ecuación diferencial ordinaria en el problema contiene un parámetro Ȝ.
Los valores de Ȝ para los que el PVF tiene soluciones no triviales llamados eigenvalores y las
soluciones correspondientes se llaman eigenfunciones. Los problemas con valores en la frontera
de esta clase son especialmente importantes en los capítulos 12 y 13. En esta sección también
vemos que existe una conexión entre los conjuntos ortogonales y las eigenfunciones de un
problema con valores en la frontera.

REPASO DE LAS ED Por conveniencia, repasaremos aquí algunas EDO y sus soluciones generales
que se presentarán con frecuencia en las secciones y capítulos siguientes. El símbolo Į representa
una constante.

Ecuaciones coeficientes constantes Soluciones generales

Y’+ ay =0 Y=𝑐1 𝑒 −𝑎𝑥
Y’’+𝑎2 𝑦 = 0 , a>0 Y=𝑐1 𝑐𝑜𝑠 𝑎𝑥 + 𝑐2 𝑠𝑒𝑛 𝑎𝑥
Y’’-𝑎2 𝑦 = 0 , a>0 𝑦 = 𝑐1 𝑒 −𝑎𝑥 + 𝑐2 𝑒 𝑎𝑥
{𝑦 = 𝑐 𝑐𝑜𝑠ℎ 𝑎𝑥 + 𝑐 𝑠𝑒𝑛ℎ 𝑎𝑥
1 2

Ecuación de Cauchy-Euler Soluciones generales, x 0

𝑥 2 𝑦′′ + 𝑥𝑦′ − 𝑎2 𝑦 = 0, 𝑎 ≥ 0 𝑦 = 𝑐1 𝑥 −𝑎 + 𝑐2 𝑥 𝑎 , 𝑎 > 0
{ 𝑦 = 𝑐 + 𝑐 𝑙𝑛 𝑥, 𝑎 = 0
1 2

Ecuación paramétrica de Bessel (v 0) Solución general, x 0

x𝑦′′ + 𝑦′ − 𝑎2 𝑥𝑦 = 0, 𝑎 ≥ 0 Y= 𝑐1 𝐽0 (𝑎𝑥) + 𝑐2 𝑌0 (𝑎𝑥)

Ecuación de Legendre (n=0,1,2,…. ) Las soluciones particulares son polinomios

(1 - 𝑥 2 )𝑦 𝑙𝑙 - 2xy’ +n(n+1)y =0 Y= 𝑃0 (𝑥) = 1,
Y= 𝑃1 (𝑥) = 𝑥,
1
Y= 𝑃2 (𝑥) = (3𝑥 2 − 1), . . . . ..
2

Considerando las dos formas de la solución general de y’’- 𝑎2 𝑦 = 0 en el ejemplo 1 haremos uso
inmediatamente de la siguiente regla informal así como en análisis futuros:
Contreras Roque Adal Eduardo 201935255

Utilice la forma exponencial y= 𝑐1 𝑒 −𝑎𝑥 +𝑐2 𝑒 𝑎𝑥 cuando el dominio de x es un intervalo infinito o

semiinfinito; utilice la forma hiperbólica y = 𝑐1 𝑐𝑜𝑠ℎ 𝑎𝑥 + 𝑐2 𝑠𝑒𝑛ℎ 𝑎𝑥 cuando el dominio de x es un
intervalo finito.

11.4 PROBLEMA DE STURM-LIOUVILLE

REPASO DE LAS ED Por conveniencia, repasaremos aquí algunas EDO y sus soluciones
generales que se presentarán con frecuencia en las secciones y capítulos siguientes. El
símbolo Į representa una constante.

Ecuaciones con coeficientes constantes Soluciones generales

Y'+ αy = 0 y =c1 e-α x

y'’ + α2y=0, α>0 y=c
1 cosh αx +c2 senh αx
y'’ - α2y=0, α>0 {y= c1 e-α x c2 eα x, o
{y = c1 cosh αx +c2 senh αx

Ecuación de Cauchy-Euler Soluciones generales, x>0

x2y''+xy’- α2y=0, a ≥0 {y= c1 x-α x c2 xα, α>0

{y = c1 + c2 ln x, α=0

ecuación paramétrica de Blessel (v = 0) Solución general, x > cero

x2y''+y’- α2xy=0, a ≥0 y =c1J0(αx) + c2Y0(αx)

Ecuación de Legendre (n=0,1,...) Las soluciones particulares son polinomios

(1-x2)y’’ - 2xy’ + n(n+1)y =0, y =P0(x) =1,

y =P(x) =x,
y = P2(x) = ½(3x2-1),...

Considerando las dos formas de la solución general de y2-α2y=0, en el ejemplo 1 haremos

uso inmediatamente de la siguiente regla informal así como en análisis futuros:
Contreras Roque Adal Eduardo 201935255

Utilice la forma exponencial y = c1e-αx c2 e-αx cuando el dominio de x es un intervalo infinito

o semi infinito; utilice la forma hiperbólica y =c1 cosh αx +c2 senh αx cuando el dominio de
x es un intervalo infinito

EIGENVALORES Y EIGENFUNCIONES Las funciones ortogonales surgen al resolver

ecuaciones diferenciales. Además, se puede generar un conjunto ortogonal de funciones
al resolver un problema con valores en la frontera con dos puntos que implique una
ecuación diferencial de segundo orden lineal que tenga un parámetro λ. En el ejemplo 2
de la sección 5.2, vimos que el problema con valores en la frontera

Y'’ + λy = 0, y(0) =0, y(L) =0, (1)

tiene soluciones no triviales sólo cuando el parámetro λ toma los valores λn =n2 π2 / L2 , n
1, 2, 3, . . . , llamados eigenvalores. Las correspondientes soluciones no triviales yn c2
sen(Qʌ[L) o simplemente yn sen(Qʌ[L) se llaman eigenfunciones del problema. Por
ejemplo, para el problema con valores en la frontera (1),
El -2 No es un eigenvalor
PVF: y’’ -2y =0, y(0) =0, y(L) = 0
Solución trivial y = 0 <-- nunca es una eigenfunción

El -9 π2 / L2 <-- es un eigenvalor (n =3)

PVF: y’’ -9 π2 / L2 y =0, y(0) =0, y(L) = 0
Solución no trivial y3 = sen(3πx / L) es una eigenfunción

Para nuestros fines en este capítulo es importante reconocer que el conjunto {sen(nπx/L},
n =1,2,3... es el conjunto ortogonal de las funciones en el intervalo [0,L] que se usa como
base para la serie de Fourier de senos.

EJEMPLO 1 Eigenvalores y eigenfunciones

Considere el problema con valores en la frontera
y'’ + λy = 0, y’(0) =0, y’(L) =0, (2)
Contreras Roque Adal Eduardo 201935255

Como en el ejemplo hay tres posibles casos para el parámetro λ: cero, negativo o positio;
esto es, λ =0, λ = -α2<0 y λ = α2 > 0 , donde α >0. La solucion de las ED

y'' = 0, λ = 0 (3)
y'’ - α2y=0, λ = -α2 (4)
y'’ + α2y=0, λ = -α2 (5)
Son, respectivamente
y =c1 +c2x, (6)
y=c1 cosh αx +c2 senh αx (7)
y=c1 cos αx +c2 sen αx (8)
Cuando las condiciones en la frontera, y’(0) =0, y’(L)=0 se aplican a cada una de estas
soluciones, de la ecuación (6) se obtiene y =c1, de la ecuación (7) se obtiene y = 0 y de la
ecuación (8) se obtiene y = c1 cos αx suponiendo que α= nπ/L, n =1,2,3,… Puesto que
y =c1 satisface que la ED en (3) y las ecuaciónes de frontera para cualquier
elección distina de cero, concluimos que λ=0 es un eigenvalor. Por lo que los eigenvalores
y las correspondientes eigenfunciones del problema son λ0=0, y0=c1,c1 ≠ 0
y λn = α2n= n2 π2/L2, n = 1,2,…, yn = c1cos (nπx/L) si hacemos que n = 0 . El
conjunto {cos(nπx/L)}, n = 0,1,2,3,,…, es ortogonal en el intervalo [0,L]

PROBLEMA REGULAR DE STURM-LIOUVILLE Los problemas (1) y (2) son casos especiales
de un problema importante con valores en la frontera de dos puntos. Sean p, q,
r’ funciones de valor real continuas en un interalo [a,b] y sean r(x) > 0 y p(x) > 0 para todo
x en el intervalo. Entonces

D/dx [r(x)y’] + (q(x) + λp(x))y = 0 (9)

A1y(a) + B1y’(a) = 0 (10)
A1y(b) + B1y’(b) = 0 (11)
Se dice que es un problema regular de Sutrum-Liouville. Los coeficientes en las
condiciones de frontera (10)(11) se suponen reales e independientes de λ. Además, A1 y
B1 no son iguales a cero y A2 y B2 no son iguales a cero. Los problemas con valores en la
frontera en (1) y (2) son problemas regulares de Strum-Liouville. De (1) podemos
identificar r(x) = 1, q(x) = 0 y p(x) =1 en la ecuación diferencial (9) en la condición frontera
(10) identificamos a = 0, A1=1, B1=0 y en (11), b =L, A2=1, B2= 0 De (2) las
identificaciones serán a=0, A1=0, B1=1 en (10), b =L, A2=0, B2= 1 en (11)
Contreras Roque Adal Eduardo 201935255

La ecuación diferenciao (9) es lineal y homogénea. Las condiciones de frontera en (10) y

(11), ambas una combinación lineal de y y y’ son iguales a cero en un punto y son
también homogéneas. Una condición de frontera tal como A2 y(b) + B2y’(b) = C2,
donde C2 es una constante diferente de cero, es no homogénea. Un problema con valores
en la frontera que consiste en una ecuación diferencial lineal homogénea y de condiciones
en la frontera homogéneas es, por supuesto, llamado un PVF homogéneo; de otra
manera, es no homogéneo. Las condiciones en la frontera (10) y (11) se
llaman separadas porque cada condición implica sólo un punto en la frontera.
Puesto que un problema regular de Sturm-Liouville es un PVF homogéneo, tiene siempre
la solución trivial y = 0. Sin embargo, esta solución no es de interés para nosotros. Como
en el ejemplo 1, al resolver uno de estos problemas tratamos de buscar números λ
(eigenvalores) y soluciones no triviales y que dependan de λ (eigenfunciones).

PROPIEDADES El teorema 11.4.1 es una lista de las propiedades más importantes del
problema regular de Strum-Liouville. Sólo demostramos la última propiedad.

Teorema 11.4.1 Propiedades del problema regular de Strum-Liouville

a) Existe un número infinito de eigenvalores reales que se pueden ordeenar en forma
creciente, λ1 < λ2< λ3<…tal que λn --> ∞ conforme n --> ∞
b) Para cada eigenvalor existe sólo una eigenfunción (excepto los múltiplos diferentes de
cero).
c) Las eigenfunciones que corresponden a diferentes eigenvalores son linealmente
independientes.
d) El conjunto de eigenfunciones que corresponde al conjunto de los eigenvalores es
ortogonal respecto a la función de peso p(x) en el intervalo [a,b]

EJEMPLO 2 Un problema regular de Sturm-Liouville

Resuelva el problema con valores en la frontera
y'’ + λy =0. y(0) =0, y(1) + y’(1) =0 (16)
SOLUCIÓN Procedemos exactamente como en el ejemplo 1 considerando tres casos en los
que el parámetro λ podría ser cero, negativo o positivo: λ = 0, λ=- α2 < 0, y λ - α2 > 0
donde - α>0. Las soluciones de la ED para estos valores se muestran en las ecuaciones (3)
a (5). Para los casos λ = 0, λ - α2 < 0 encontramos que los PVF en (16) sólo tienen la
solución trivial y = 0. Para λ - α2 >0 la solución general de la ecuación diferencial es
y = c1 cos αx +c2 sen αx. Ahora la condición y(0)= 0 implica que en esta solución c1 = 0, así
nos quedamos con y c2 sen αx. La segunda condición y(1)+y’(1)=0 se satisface si
Contreras Roque Adal Eduardo 201935255

c2 sen α +c2 α cos α = 0

En vista del requisito que c2 ≠ 0, la última ecuación se puede escribir como
Tan α = - α (17)
Si por un momento consideramos en (17) que tan x = -x, enconces en la figura 11.4.1 se
muestra la factibilidad de que exista un número infinito de raíces, en partiular, las
coordenadas x de los puntos donde la gráfica de y 0 –x intersecta el número infinito de
ramas de gráfica de y = tan x. Los eigenvalores del PVF (16) son entonces λn = α2n,
donde αn, n 1, 2, 3, . . . son las raíces positivas consecutivas α1, α2 α3... de (17). Con ayuda
de un SAC se muestra con facilidad que redondeando a cuatro
decimales, α1,= 2.0288 α2 = 4.9132 α3.= 7.9787..y α4 =11.0855 y que las soluciones
correspondientes son y1= sen 2.0288x, y2= sen 4.9132x, y3= sen 7.9787x y
y4= sen 11.0855x. En general, las eigenfunciones del problema son {sen αn x}, n = 1, 2, 3, .
..
Identificando r(x)= 1, q(x) =0, p(x) =1, A1=1, B1=0 , A2=1, B2= 1, vemos que la ecuación (16)
es un problema regular de Sturm-Liouville. Concluimos que {sen αn x}, n = 1, 2, 3, . . . es
un cojunto ortogonal respecto a la función de peso p(x) = 1 en el intervalo [0,1]
En algunos casos se puede demostrar la ortogonalidad de las soluciones de (9) sin
necesidad de especificar una condición en la frontera en x = a y en x = b

PROBLEMA SINGULAR DE STURM-LIOUVILLE Existen otras condiciones importantes bajo

las que buscamos las soluciones no triviales de la ecuación diferencial (9):
• r (a) = 0, y una condición de frontera del tipo dado en (11) está dada como x
= b; (18)
• r (b) = 0, y una condición de frontera del tipo dado en (10) está dada como
x = a; (19)
• r (a) = r (b) = 0, y no hay condición de frontera dada en x = a o en
x = b; (20)
• r (a) = r (b), y las condiciones de frontera y(a) = y(b),
y’(a) = y’(b). (21)

La ecuación diferencial (9) junto con una de las condiciones (18) a (20), se dice que es un
problema singular con valores en la frontera. La ecuación (9) con las condiciones dadas en
(21) se dice que es un problema con valores en la frontera periódico (las condiciones de
frontera también se llaman periódicas). Observe que si decimos que r(a) = 0, entonces x =
a puede ser un punto singular de la ecuación diferencial y por tanto, una solución de (9)
puede crecer sin límite conforme x → a. Sin embargo, vemos de (14) que si r(a) = 0, no se
necesita condición de frontera en x = a para demostrar la ortogonalidad de
Contreras Roque Adal Eduardo 201935255

las eigenfunciones suponiendo que estas soluciones estén limitadas en ese punto. Este
último requisito asegura la existencia de las integrales que intervien. Suponiendo que las
soluciones de (9) estén acotadas en un intervalo cerrado [a, b], podemos ver del examen
de la ecuación (14) que

• si r(a) =0, entonces la relación de ortogonalidad (15) es válida, sin ninguna

condición dada en la frontera en x a; (22)
• si r(b) = 0, entonces la relación de ortogonalidad (15) es válida sin ninguna
condición dada en la frontera en x = b; (23)
• si r(a) r(b) = 0, entonces la relación de ortogonalidad (15) es válida sin ninguna
condición dada en la frontera en x =a o en x = b; (24)
• si r(a) = r(b), entonces la relación de ortogonalidad (15) es válida con las
condiciones en la frontera y(a)= y(b), y’(a) y’(b). (25)
Observe que un problema de Strum-Liouville es singular cuando el intervalo que se
considera es infinito. Véase los problemas 9 y 10 de los ejercicios 11.4
FORMA AUTOADJUNTA Realizando la derivación que se indica en (9), vemos que la
ecuación diferencial es igual a
r(x)y’’ +r’(x)y’ + (q(x)+ λp(x))y =0 (26)
El examen de la ecuación (26) para conducir a creer que el coeficiente dado de y’ es la
derivada del coeficiente de y’’, y que existen pocas ecuaciones diferenciales que tengan la
forma de la ecuación (9).Por lo contrario, si los coeficientes son continuos y a(x)≠ 0 para
toda x en algún intervalo, entonces cualquier ecuación diferencial de segundo orden
a(x)y’’ +b(x)y’ + (c(x)+ λd(x))y =0 (27)
se puede escribir en la así llamada forma autoadjunta (9). Para esto básicamente
reescribimos una ecuación homogénea lineal de primer orden a1(x)y’ + a0(x)y =0 en la
forma d/dx[μy] =0 dividiendo la ecuación entre a1 (x) y después multiplicando por el
factor integrante μ =
p(x)dx, donde, se supone que no hay factores comunes, P(x) = a0(x)/a1(x).Así que primero
, dividamos la ecuación (27) por a(x). Los primeros dos términos son Y’ + b(x)a(x)Y+...,
donde enfatizamos que hemos escrito Y= y’.Segundo, multiplicamos esta ecuación por el
factor integrante eb(x)a(x)dx, donde a(x) y b(x) se supone que no tienen factores en
común:

eb(x)a(x)dxY'+b(x)a(x)eb(x)a(x)dxY+...=dd(x)[eb(x)a(x)dxY]+...dd(x)[eb(x)a(x)dxy']+...
En resumen, dividiendo la ecuación (27) entre a(x) y después multiplicando, obtenemos
ebadxy''+b(x)a(x)ebadxy'+(c(x)a(x)ebadx+ d(x)a(x)ebadx)y = 0 (28)
Contreras Roque Adal Eduardo 201935255

Ciertamente no es necesario escribir una ecuación diferencial de segundo orden (27) en la

forma autoadjunta (9) para resolver ED. Para nuestros fines usaremos la forma dada en la
ecuación (9) para determinar la función de peso p(x) que se necesita en la relación de
ortogonalidad (15). Los dos ejemplos siguientes ilustran relaciones de ortogonalidad para
funciones de Bessel y para polinomios de Legendre.

EJEMPLO 3 Ecuación paramétrica de Bessel

En la sección 6.4 vimos que la solución general de la ecuación paramétrica de Bessel de

orden n es x2y''+x'y+(2x2-n2)y=0, donde n es un entero fijio no negativo y es un
parámetro positivo. La solución general de esta ecuación es y = c1J1(x)+c2Yn(x)Después de
dividir la ecuación paramétrica de Bessel entre el primer coeficiente x2 y multiplicando
por la ecuación resultante por el factor integrante e1xdx=elnx=x, x>0obtenemos

xy’’ +y’+ (2x-n2x)y=0 o ddx[xy']+(2x-n2x)y =0

Comparando este último resultado con la forma autoadjunta (9), hacemos las
idenficiaciónes r(x), y p(x). Ahora con r(0)=0 y de las dos soluciones el conjunto {Jn((ix}, i =
1,2,3...es ortogonal respecto a la función de peso p(x) = x en un intervalo [0,b]. La relación
de ortogonalidad es

0bxjn(i x)Jn(j x)dx =0 i j(29)

suponemos que i y por lo tanto los eigenvalores i = i2,i = 1,2,3,...,se definen por medio de
una condición en la frontera de x = b del tipo de dato en la ecuación (11)

A2Jn(b)+B2Jn'(b)=0 (30)

EJEMPLO 4 Ecuación de Legendre

La ecuación diferencial de Legendre(1-x2)y''- 2xy + n(n +1)y= 0 es exactamente de la
forma dada en la ecuación (26) con r(x) 1 – x2y r’(x)= 2x. Por lo que la forma autoadjunta
(9) es inmediata,

ddx[(1-x2)y'] + n(n+1)y =0 (31)

Contreras Roque Adal Eduardo 201935255

De la ecuación (31) podemos además identificar q(x)=0, =n(n+1)y p(x)= 0. Recuerde de la

sección 6.3 que cuando n= 0, 1, 2, . . . la ED de Legendre tiene soluciones polinomiales Pn
(x). Ahora se puede expresar la observación de que r(-1)= r(1) 0 junto con el hecho de que
los polinomios de Legendre Pn (x) que son las únicas soluciones de (31) que tienen límite
en el intervalo cerrado[1, 1] por lo que se concluye de la ecuación (24) que el conjunto {Pn
(x)}, n 0, 1, 2, . . . es ortogonal respecto a la función de peso p(x) =1 en [-1, 1]. La relación
de ortogonalidad es

-11Pm(x)Pn(x)dx =0, mn
Contreras Roque Adal Eduardo 201935255

Ejercicios

1.
Contreras Roque Adal Eduardo 201935255

2.
Contreras Roque Adal Eduardo 201935255

3.
Contreras Roque Adal Eduardo 201935255

4.
Contreras Roque Adal Eduardo 201935255

5.