Tarea #7
Tarea #7
Tarea #7
FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA MECÁNICA
TAREA N° 7
CURSO:
ECUACIONES DIFERENCIALES
ALUMNOS:
ESQUIVEL CASTAÑEDA HAMNER NIKOLAE
GUTIERREZ RISCO EVERTH LEONEL
PAREDES MURGA ALEXIS PAUL
PINTADO CASEUX AARON WILFREDO
QUIROZ ALVA CRISTIAN DANIEL
SOLANO HERRERA JOE
VEGA BALTODANO DIEGO ALBERTO
DOCENTE:
VARGAS VERA, LUIS FERNANDO
CICLO Y SECCIÓN:
IV CICLO – SECCIÓN A
Trujillo – Perú
2022
Sección 11.2
Ejercicio 5: Encuentre la serie de Fourier de 𝑓 en el intervalo dado.
0, −𝜋 < 𝑥 < 0
𝑓(𝑥) = { 2
𝑥 , 0≤𝑥<𝜋
Solución:
Sabemos que si tenemos la función 𝑓(𝑥) con periodo 2𝐿 y definida en el intervalo
−𝐿 < 𝑥 < 𝐿 entonces la serie de Fourier viene dada por
∞
𝑎0 𝑚𝜋𝑥 𝑚𝜋𝑥
𝑓(𝑥) = + ∑ (𝑎𝑚 cos + 𝑏𝑚 sin )
2 𝐿 𝐿
𝑚=1
Donde
1 𝐿
𝑎0 = ∫ 𝑓(𝑥)𝑑𝑥
𝐿 −𝐿
1 𝐿 𝑚𝜋𝑥
𝑎𝑚 = ∫ 𝑓(𝑥) cos 𝑑𝑥 , 𝑚 = 1, 2, …
𝐿 −𝐿 𝐿
1 𝐿 𝑚𝜋𝑥
𝑏𝑚 = ∫ 𝑓(𝑥) sin 𝑑𝑥 , 𝑚 = 1, 2, …
𝐿 −𝐿 𝐿
Ahora, evaluamos los coeficientes 𝑎0 , 𝑎𝑚 y 𝑏𝑚 para la función dada. Ahora vamos a
evaluar 𝑎0 reemplazando 𝐿 = 𝜋.
1 𝜋
𝑎0 = ∫ 𝑓(𝑥)𝑑𝑥
𝜋 −𝜋
1 0 1 𝜋
𝑎0 = ∫ 0𝑑𝑥 + ∫ 𝑥 2 𝑑𝑥
𝜋 −𝜋 𝜋 0
1 3𝜋
𝑎0 = [𝑥 ]0
3𝜋
1 3
𝑎0 = 𝜋
3𝜋
𝜋2
𝑎0 =
3
Y
𝑓(−𝜋 + ) = lim + 0 = 0
𝑥→−𝜋
𝜕 2𝑢 𝜕𝑢
𝑘 2
− ℎ𝑢 =
𝜕𝑥 𝜕𝑡
Solución:
Tendremos:
𝑢(0, 𝑡) = 0 , 𝑢(𝐿, 𝑡) = 0 𝑡>0
𝑢(𝑥, 0) = 𝑓(𝑥) , 0<𝑥<𝐿
Utilizando el método de separación de variables obtenemos:
𝑢(𝑥, 𝑡) = 𝑋(𝑥)𝑇(𝑡)
𝜕 2𝑢
= 𝑋 ′′ (𝑥)𝑇(𝑡)
𝜕𝑥 2
𝜕𝑢
= 𝑋(𝑥)𝑇 ′ (𝑡)
𝜕𝑡
𝜕𝑢
= 𝑋′(𝑥)𝑇(𝑡)
𝜕𝑡
Donde al sustituir en la ecuación inicial:
𝑘𝑋 ′′ (𝑥)𝑇(𝑡) − ℎ𝑋(𝑥)𝑇(𝑡) = 𝑋(𝑥)𝑇 ′ (𝑡)
Realizaremos una división que nos permita darle forma a esta nueva ecuación:
𝑋 ′′ (𝑥)𝑇(𝑡) − ℎ𝑋(𝑥)𝑇(𝑡) 𝑋(𝑥)𝑇 ′ (𝑡)
=
𝑘𝑇(𝑡) 𝑘𝑇(𝑡)
ℎ𝑋(𝑥) 𝑋(𝑥)𝑇 ′ (𝑡)
𝑋 ′′ (𝑥) − =
𝑘 𝑘𝑇(𝑡)
1 ℎ𝑋(𝑥) 𝑇 ′ (𝑡)
(𝑋 ′′ (𝑥) − )=
𝑋(𝑥) 𝑘 𝑘𝑇(𝑡)
𝑋 ′′ (𝑥) ℎ𝑋(𝑥) 𝑇 ′ (𝑡)
− =
𝑋(𝑥) 𝑋(𝑥)𝑘 𝑘𝑇(𝑡)
𝑋 ′′ (𝑥) ℎ 𝑇 ′ (𝑡)
− =
𝑋(𝑥) 𝑘 𝑘𝑇(𝑡)
𝑋 ′′ (𝑥) 𝑇 ′ (𝑡) ℎ
= +
𝑋(𝑥) 𝑘𝑇(𝑡) 𝑘
𝑋 ′′ (𝑥) 1 𝑇 ′ (𝑡)
= ( + ℎ)
𝑋(𝑥) 𝑘 𝑇(𝑡)
Asumimos que:
𝑋 ′′ (𝑥) 1 𝑇 ′ (𝑡)
= ( + ℎ) = −𝜆2
𝑋(𝑥) 𝑘 𝑇(𝑡)
Tendremos entonces dos términos con el mismo valor a los cuales llamaremos:
𝑋 ′′ (𝑥)
= −𝜆2 (𝑖)
𝑋(𝑥)
1 𝑇 ′ (𝑡)
( + ℎ) = −𝜆2 (𝑖𝑖)
𝑘 𝑇(𝑡)
Para 𝑖:
𝑋 ′′ (𝑥)
= −𝜆2
𝑋(𝑥)
𝑋 ′′ (𝑥) + 𝜆2 𝑋(𝑥) = 0
Considerando un 𝜆 ≠ 0, la ecuación auxiliar será:
𝑚2 + 𝜆2 = 0
𝑚2 = −𝜆2
𝑚 = ±𝑖𝜆
Tendremos la solución:
𝑋(𝑥) = 𝐶1 cos(𝜆𝑥) + 𝐶2 𝑠𝑒𝑛(𝜆𝑥)
Ahora tomando 𝑖𝑖:
1 𝑇 ′ (𝑡)
( + ℎ) = −𝜆2
𝑘 𝑇(𝑡)
𝑇 ′ (𝑡) + (𝑘𝜆2 + ℎ)𝑇(𝑡) = 0
Considerando un 𝜆 ≠ 0:
𝑇 ′ (𝑡) = −(𝑘𝜆2 + ℎ)𝑇(𝑡)
𝑇 ′ (𝑡)
= −(𝑘𝜆2 + ℎ)
𝑇(𝑡)
Integrando esta ecuación tenemos:
𝑇 ′ (𝑡)
∫ 𝑑𝑡 = ∫ −(𝑘𝜆2 + ℎ) 𝑑𝑡
𝑇(𝑡)
ln(𝑇(𝑡)) = −(𝑘𝜆2 + ℎ)𝑡 + ln(𝑐)
ln(𝑇(𝑡)) − ln(𝑐) = −(𝑘𝜆2 + ℎ)𝑡
𝑇(𝑡)
ln ( ) = −(𝑘𝜆2 + ℎ)𝑡
𝑐
Tendremos la solución:
2 +ℎ)𝑡
𝑇(𝑡) = 𝐶𝑒 −(𝑘𝜆
Si aplicamos las condiciones iniciales del problema.
Para 𝑢(0, 𝑡) = 0:
𝑢(𝑥, 𝑡) = 𝑋(𝑥)𝑇(𝑡)
𝑢(0, 𝑡) = 𝑋(0)𝑇(𝑡) = 0
𝑋(0) = 0
𝑋 ′ (0) = 0
Sustituyendo en la solución de 𝑖:
𝑋(𝑥) = 𝐶1 cos(𝜆𝑥) + 𝐶2 𝑠𝑒𝑛(𝜆𝑥)
𝑋′(𝑥) = −𝜆𝐶1 sen(𝜆𝑥) + 𝜆𝐶2 𝑐𝑜𝑠(𝜆𝑥)
𝑋′(0) = −𝜆𝐶1 sen(0) + 𝜆𝐶2 𝑐𝑜𝑠(0)
𝑋 ′(0) = 0 + 𝜆𝐶2
Ya que 𝜆 ≠ 0:
𝐶2 = 0
Tenemos entonces:
𝑋(𝑥) = 𝐶1 cos(𝜆𝑥)
Para 𝑢(𝐿, 𝑡) = 0:
𝑢(𝑥, 𝑡) = 𝑋(𝑥)𝑇(𝑡)
𝑢(𝐿, 𝑡) = 𝑋(𝐿)𝑇(𝑡) = 0
𝑋(𝐿) = 0
𝑋 ′ (𝐿) = 0
Sustituyendo en la solución de 𝑖:
𝑋(𝑥) = 𝐶1 cos(𝜆𝑥) + 𝐶2 𝑠𝑒𝑛(𝜆𝑥)
𝑋′(𝑥) = −𝜆𝐶1 sen(𝜆𝑥) + 𝜆𝐶2 𝑐𝑜𝑠(𝜆𝑥)
𝑋′(𝐿) = −𝜆𝐶1 sen(𝜆𝐿) + 𝜆𝐶2 𝑐𝑜𝑠(𝐿)
Sabiendo que 𝐶2 = 0:
𝑋′(𝐿) = −𝜆𝐶1 sen(𝜆𝐿)
−𝜆𝐶1 sen(𝜆𝐿) = 0
Para una solución no trivial, se tiene que 𝐶1 ≠ 0
Por lo tanto:
sen(𝜆𝐿) = sen(𝑛𝜋)
𝜆𝐿 = 𝑛𝜋
𝑛𝜋
𝜆=
𝐿
𝑛𝜋
Para 𝜆, 𝐶1 ≠ 0, los autovalores son 𝜆𝑛 = , donde 𝑛 = 1,2,3, …
𝐿
Donde:
2 𝐿 𝑛𝜋𝑥
𝑎𝑛 = ∫ 𝑓(𝑥)𝑠𝑒𝑛 ( ) 𝑑𝑥
𝐿 0 𝐿
Ejercicio 6: Resuelva el problema 5, si los extremos 𝑥 = 0 𝑦 𝑥 = 𝐿, se mantienen a
temperatura cero.
Solución:
Del ejercicio 5 tenemos la siguiente ecuación:
𝐾𝜕 2 𝑢 𝜕𝑢
2
− ℎ𝑢 = , 0 < 𝑥 < 𝐿, 𝑡 > 0 𝑑𝑜𝑛𝑑𝑒 ℎ 𝑒𝑠 𝑢𝑛𝑎 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
𝜕𝑥 𝜕𝑡
Como nos dicen que los extremos 𝑥 = 0 𝑦 𝑥 = 𝐿, se mantienen a temperatura cero,
entonces, tenemos las condiciones:
𝑢(0, 𝑡) = 0, 𝑢(𝐿, 𝑡) = 0, 𝑝𝑎𝑟𝑎 𝑡 > 0
Y, además:
𝑢(𝑥, 0) = 𝑓(𝑥)
Por el método de separación de variables, denotamos:
𝑢(𝑥, 𝑡) = 𝑀(𝑥)𝑁(𝑡)
Entonces tendríamos:
𝜕 2𝑢
= 𝑀′′(𝑥)𝑁(𝑡)
𝜕𝑥 2
𝜕𝑢
= 𝑀(𝑥)𝑁′(𝑡)
𝜕𝑡
Reemplazando y operando:
𝐾𝑀′′𝑁 − ℎ𝑀𝑁 = 𝑀𝑁′
𝐾𝑀′′𝑁 ℎ𝑀𝑁 𝑀𝑁′
− =
𝑀 𝑀 𝑀
𝐾𝑀′′𝑁
− ℎ𝑁 = 𝑁′
𝑀
𝐾𝑀′′𝑁
= 𝑁 ′ + ℎ𝑁
𝑀
𝑀′′ 𝑁 ′ + ℎ𝑁
=
𝑀 𝐾𝑁
𝑀′′ 𝑁 ′ ℎ𝑁
= +
𝑀 𝐾𝑁 𝐾𝑁
𝑀′′ 𝑁 ′ ℎ
= +
𝑀 𝐾𝑁 𝐾
𝑀′′ 1 𝑁′
= ( ℎ + ) = −𝜆
𝑀 𝐾 𝑁
Por lo que tendríamos:
𝑀′′ + 𝜆𝑀 = 0 y 𝑁 ′ + (ℎ + 𝐾𝜆)𝑁 = 0
Ahora, resolveremos estas ecuaciones para cuando: 𝜆 = 0, 𝜆 = −𝛼 2 < 0 𝑦 𝜆 = 𝛼 2 > 0
Para 𝜆 = 0:
𝑀′′ + (0)𝑀 = 0
𝑀′′ = 0
𝑀 = 𝑐1 𝑥 + 𝑐2
Y
𝑁 ′ + (ℎ + 𝐾(0))𝑁 = 0
𝑁 ′ + ℎ𝑁 = 0
𝑁 = 𝑐2 𝑒 −ℎ𝑡
Entonces, 𝑢(𝑥, 𝑡) sería:
𝑢(𝑥, 𝑡) = 𝑀(𝑥)𝑁(𝑡) = (𝑐1 𝑥 + 𝑐2 ) (𝑐3 𝑒 −ℎ𝑡 )
Para 𝜆 = −𝛼 2 :
𝑀′′ − 𝛼 2 𝑀 = 0
𝑀 = 𝑐4 cosh(𝛼𝑥) + 𝑐5 sinh(𝛼𝑥)
Y
𝑁 ′ + (ℎ − 𝛼 2 𝐾)𝑁 = 0
2 𝐾)𝑡
𝑁 = 𝑐6 𝑒 −(ℎ−𝛼
Entonces, 𝑢(𝑥, 𝑡) sería:
2 𝐾)𝑡
𝑢(𝑥, 𝑡) = 𝑀(𝑥)𝑁(𝑡) = (𝑐4 cosh(𝛼𝑥) + 𝑐5 sinh(𝛼𝑥)) (𝑐6 𝑒 −(ℎ−𝛼 )
Para 𝜆 = 𝛼 2 :
𝑀′′ + 𝛼 2 𝑀 = 0
𝑀 = 𝑐7 cos(𝛼𝑥) + 𝑐8 sin(𝛼𝑥)
Y
𝑁 ′ + (ℎ + 𝛼 2 𝐾)𝑁 = 0
2 𝐾)𝑡
𝑁 = 𝑐9 𝑒 −(ℎ+𝛼
Entonces, 𝑢(𝑥, 𝑡) sería:
2 𝐾)𝑡
𝑢(𝑥, 𝑡) = 𝑀(𝑥)𝑁(𝑡) = (𝑐7 cos(𝛼𝑥) + 𝑐8 sin(𝛼𝑥)) (𝑐9 𝑒 −(ℎ+𝛼 )
Analizamos por casos:
Reemplazamos las condiciones de frontera 𝑢(0, 𝑡) = 0, 𝑢(𝐿, 𝑡) = 0 en (𝑐1𝑥 +
𝑐2 ) (𝑐3 𝑒 −ℎ𝑡 )
𝑢(0, 𝑡) = (𝑐1 (0) + 𝑐2 ) (𝑐3 𝑒 −ℎ𝑡 ) = 0
𝑐2 = 0
𝑢(𝐿, 𝑡) = (𝑐1 (𝐿) + 0) (𝑐3 𝑒 −ℎ𝑡 ) = 0
𝑐1 = 0
Con esto, tendríamos la solución trivial, así que este caso no cumple.
Ahora reemplazamos las condiciones de frontera 𝑢(0, 𝑡) = 0, 𝑢(𝐿, 𝑡) = 0 en 𝑢(𝑥, 𝑡) =
2 𝐾)𝑡
𝑀(𝑥)𝑁(𝑡) = (𝑐4 cosh(𝛼𝑥) + 𝑐5 sinh(𝛼𝑥)) (𝑐6 𝑒 −(ℎ−𝛼 ):
2 𝐾)𝑡
𝑢(0, 𝑡) = (𝑐4 cosh(0) + 𝑐5 sinh(𝑥)) (𝑐6 𝑒 −(ℎ−𝛼 )=0
2 𝐾)𝑡
𝑢(0, 𝑡) = (𝑐4 (1) + 𝑐5 (0)) (𝑐6 𝑒 −(ℎ−𝛼 )=0
𝑐4 = 0
2 𝐾)𝑡
𝑢(𝐿, 𝑡) = ((0) + 𝑐5 sinh(𝛼𝐿) (𝑐6 𝑒 −(ℎ−𝛼 )=0
𝑐5 = 0 𝑜 𝑐6 = 0
Con lo que tenemos la solución trivial, por lo que este caso tampoco cumple.
Ahora reemplazamos las condiciones de frontera para (𝑐7 cos(𝛼𝑥) +
2 𝐾)𝑡
𝑐8 sin(𝛼𝑥)) (𝑐9 𝑒 −(ℎ+𝛼 ):
2 𝐾)𝑡
𝑢(0, 𝑡) = (𝑐7 cos(0) + 𝑐8 sin(0)) (𝑐9 𝑒 −(ℎ+𝛼 )=0
2 𝐾)𝑡
𝑢(0, 𝑡) = (𝑐7 ) (𝑐9 𝑒 −(ℎ+𝛼 )=0
𝑐7 = 0
2 𝐾)𝑡
𝑢(𝐿, 𝑡) = (0 cos(𝐿) + 𝑐8 sin(𝐿)) (𝑐9 𝑒 −(ℎ+𝛼 )=0
2 𝐾)𝑡
𝑢(𝐿, 𝑡) = (𝑐8 sin(𝛼𝐿)) (𝑐9 𝑒 −(ℎ+𝛼 )=0
𝑐8 = 0 𝑜 sin(𝛼𝐿) = 0
Como 𝑐8 = 0 da la solución trivial, entoncessin(𝛼𝐿) = 0.
Entonces:
sin(𝛼𝐿) = 0.
𝛼𝐿 = 𝑛𝜋
𝑛𝜋
𝛼=
𝐿
Reemplazamos y tenemos:
𝑛𝜋 𝑛𝜋 2
−(ℎ+( ) 𝐾)𝑡
𝑢(𝑥, 𝑡) = (𝑐8 sin ( 𝐿)) (𝑐9 𝑒 𝐿 )
𝐿
𝑛𝜋 𝑛𝜋 2
−(ℎ+( ) 𝐾)𝑡
𝑢(𝑥, 𝑡) = 𝑐𝑛 sin ( 𝐿) 𝑒 𝐿
𝐿
Por el principio de superposición:
𝑛𝜋 2 𝑛𝜋
−(ℎ+( ) 𝐾)𝑡
𝑢𝑛 = 𝐴𝑛 𝑒 𝐿 sin ( 𝑥)
𝐿
∞
𝑛𝜋 2 𝑛𝜋
−(ℎ+( ) 𝐾)𝑡
𝑢(𝑥, 𝑡) = ∑ 𝐴𝑛 𝑒 𝐿 sin ( 𝑥)
𝐿
𝑛=1
Ahora para 𝑢.
⇒ 𝑢(𝑥, 𝑦) = 𝑋𝑌
𝜆+√𝜆2 −4 𝜆−√𝜆2 −4
𝑢(𝑥, 𝑦) = 𝐶1 𝑒 𝜆𝑥
(𝐶2 𝑒 2 𝑦 + 𝐶3 𝑒 2 𝑦 )
𝜆+√𝜆2 −4 𝜆−√𝜆2 −4
𝐶1 𝐶2 𝑒 2 𝑦 + 𝐶1 𝐶3 𝑒 2 𝑦 = 2𝑒 𝑦
Simplificamos constantes.
𝐴1 = 𝐶1 𝐶2 𝐵1 = 𝐶1 𝐶3
𝜆+√𝜆2 −4 𝜆−√𝜆2 −4
𝐴1 𝑒 2 𝑦 + 𝐵1 𝑒 2 𝑦 = 𝑎𝑒 𝑦 + 𝑏𝑒 𝑦
Donde:
𝑎+𝑏 = 2
Igualamos términos de cada lado.
𝜆+√𝜆2 −4
𝑦
⇒ 𝐴1 𝑒 2 = 𝑎𝑒 𝑦
𝜆−√𝜆2 −4
⇒ 𝐵1 𝑒 2 𝑦 = 𝑏𝑒 𝑦
Las constantes nos quedan:
𝐴1 + 𝐵2 = 2
Para el valor de 𝜆.
𝜆+√𝜆2 −4 𝜆−√𝜆2 −4
𝑦 𝑦
𝑒 2 = 𝑒𝑦 = 𝑒 2
𝜆=2
Finalmente 𝑢(𝑥, 𝑦) nos queda:
2+√22 −4 2−√22 −4
𝑦 𝑦
𝑢(𝑥, 𝑦) = 𝑒 2𝑥 (𝐴1 𝑒 2 + 𝐵1 𝑒 2 )
𝑢(𝑥, 𝑦) = 𝑒 2𝑥 (𝐴1 𝑒 𝑦 + 𝐵1 𝑒 𝑦 )
𝒖(𝒙, 𝒚) = 𝟐𝒆𝟐𝒙+𝒚