Variacion de Parametros y Transformada de Laplace

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 14

VARIACIÓN DE PARAMETROS.

1- Resolver por el método de variación de parámetros la ecuación diferencial

𝒚′′ + 𝟒𝒚 = 𝟐𝑪𝒐𝒔 𝟐𝒙

𝑃𝑟𝑖𝑚𝑒𝑟𝑜 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑚𝑜𝑠 𝑙𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 ℎ𝑜𝑚𝑜𝑔𝑒𝑛𝑒𝑎:

𝑦 ′′ + 4𝑦 = 0 − −→ 𝑚2 + 4 = 0 𝐷𝑒 𝑠𝑢𝑠 𝑟𝑎í𝑐𝑒𝑠 𝑒𝑠𝑡á 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑑𝑜 𝑝𝑜𝑟:

𝑚 = √−4 = ±2𝑖 ∴ 𝑚1 = 2𝑖 𝑚2 = −2𝑖 𝑚 = 𝛼 ± 𝛽𝑖

𝑦ℎ = 𝑒 𝛼𝑥 (𝐴𝑐𝑜𝑠 𝛽𝑥 + 𝐵𝑠𝑖𝑛 𝛽𝑥) = 𝑒 0(𝑥) (𝐴𝑐𝑜𝑠 2𝑥 + 𝐵𝑠𝑖𝑛2𝑥) = 𝐴𝑐𝑜𝑠 2𝑥 + 𝐵𝑠𝑖𝑛2𝑥

𝑦ℎ = 𝐴𝑐𝑜𝑠 2𝑥 + 𝐵𝑠𝑖𝑛2𝑥 ∴ 𝑦1 = cos 2𝑥 𝑦2 = sin 2𝑥


𝒚𝟏 𝒚𝟐
𝐿𝑜 𝑠𝑖𝑔𝑢𝑖𝑒𝑛𝑡𝑒 𝑒𝑠 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑟 𝑒𝑙 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒𝑙 𝑊𝑟𝑜𝑛𝑠𝑘𝑖𝑎𝑛𝑜 𝑤(𝑦1 , 𝑦2 ) = | |
𝒚´𝟏 𝒚´𝟐
𝑦1 = cos 2𝑥 𝑦2 = sin 2𝑥 𝑦1′ = −2sin 2𝑥 𝑦2′ = 2cos 2𝑥
cos 2𝑥 sin 2𝑥
𝑤(𝑦1 , 𝑦2 ) = | | = 2 cos2 2𝑥 − (−2sin2 2𝑥) = 2(cos 2 2𝑥 + sin2 2𝑥) = 2
−2sin 2𝑥 2cos 2𝑥
𝑤=2

𝐸𝑙 𝑠𝑖𝑔𝑢𝑖𝑒𝑛𝑡𝑒 𝑝𝑎𝑠𝑜 𝑒𝑠 𝑖𝑑𝑒𝑛𝑡𝑖𝑓𝑖𝑐𝑎𝑟 𝑙𝑎 𝑓𝑢𝑛𝑐𝑖ó𝑛 𝑟(𝑥) = 2 cos 2𝑥


𝑂𝑏𝑡𝑒𝑛𝑖𝑒𝑛𝑑𝑜 𝑒𝑠𝑜𝑠 𝑑𝑎𝑡𝑜𝑠, 𝑝𝑟𝑜𝑐𝑒𝑑𝑒𝑚𝑜𝑠 𝑎 𝑠𝑢𝑠𝑡𝑖𝑡𝑢𝑖𝑟 𝑙𝑎𝑠 𝑓𝑢𝑛𝑐𝑖𝑜𝑛𝑒𝑠 𝑦1 , 𝑦2 , 𝑟(𝑥) 𝑦 𝑤(𝑥) 𝑒𝑛 𝑙𝑎 𝑟𝑒𝑙𝑎𝑐𝑖ó𝑛 𝑝𝑎𝑟𝑎
𝑦2 𝑟(𝑥) 𝑦 𝑟(𝑥)
𝑝𝑎𝑟𝑎 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑟 𝑙𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑝𝑎𝑟𝑡𝑖𝑐𝑢𝑙𝑎𝑟 𝑦𝑝 = −𝑦1 ∫ 𝑑𝑥 + 𝑦2 ∫ 1 𝑑𝑥
𝑤 𝑤
sin 2𝑥 (2 cos 2𝑥) cos 2𝑥 (2 cos 2𝑥)
𝑦𝑝 = − cos 2𝑥 ∫ 𝑑𝑥 + sin 2𝑥 ∫ 𝑑𝑥 =
2 2
cos 2𝑥
𝑦𝑝 = − ∫ sin 2𝑥 2 cos 2𝑥 𝑑𝑥 + sin 2𝑥 ∫ cos 2𝑥 ∗ cos 2𝑥 𝑑𝑥 =
2
cos 2𝑥 (sin 4𝑥)
𝑦𝑝 = − + sin 2𝑥 ∫ cos 2 2𝑥 𝑑𝑥 =
2
cos 2𝑥 (sin 4𝑥 ) 1 1
𝑦𝑝 = − + sin 2𝑥 ( 𝑥 + 𝑠𝑖𝑛4𝑥) =
2 2 8
cos 2𝑥 (sin 4𝑥) xsin 2𝑥 sin 2𝑥 sin2 4𝑥
𝑦𝑝 = − + +
2 2 4
2- Resolver la ecuación 𝒚′′ − 𝒚′ = 𝒆𝒙 𝒄𝒐𝒔 𝒙

𝑃𝑟𝑖𝑚𝑒𝑟𝑜 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑚𝑜𝑠 𝑙𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 ℎ𝑜𝑚𝑜𝑔𝑒𝑛𝑒𝑎:

𝑦 ′′ − 𝑦 ′ = 0 − −−→ 𝑚2 − 𝑚 = 0

−𝑏 ± √𝑏 2 − 4𝑎𝑐
𝑚1,2 =
2𝑎
−(−1) ± √(−1)2 − 4(1)(0) 1 ± 1
𝑚1,2 = =
2(1) 2
1+1 1−1
𝑚1 = =1 𝑚2 = = 0 𝑐𝑜𝑚𝑜 𝑚1 ≠ 𝑚2
2 2
𝑦ℎ = 𝑐1 𝑒 𝑥 + 𝑐2 𝑒 (0)𝑥 = 𝑐1 𝑒 𝑥 + 𝑐2 − −−→ 𝑦1 = 𝑒 𝑥 𝑦2 = 1
𝒚𝟏 𝒚𝟐
𝐿𝑜 𝑠𝑖𝑔𝑢𝑖𝑒𝑛𝑡𝑒 𝑒𝑠 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑟 𝑒𝑙 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒𝑙 𝑊𝑟𝑜𝑛𝑠𝑘𝑖𝑎𝑛𝑜 𝑤(𝑦1 , 𝑦2 ) = | |
𝒚´𝟏 𝒚´𝟐
𝑦1 = 𝑒 𝑥 𝑦2 = 1 𝑦1′ = 𝑒 𝑥 𝑦2′ = 0
𝑥
𝑤(𝑦1 , 𝑦2 ) = |𝑒𝑥 1| = 0 − 𝑒𝑥 = −𝑒𝑥
𝑒 0
𝑤 = −𝑒 𝑥
𝐸𝑙 𝑠𝑖𝑔𝑢𝑖𝑒𝑛𝑡𝑒 𝑝𝑎𝑠𝑜 𝑒𝑠 𝑖𝑑𝑒𝑛𝑡𝑖𝑓𝑖𝑐𝑎𝑟 𝑙𝑎 𝑓𝑢𝑛𝑐𝑖ó𝑛 𝑟(𝑥) = 𝑒 𝑥 𝑐𝑜𝑠 𝑥
𝑂𝑏𝑡𝑒𝑛𝑖𝑒𝑛𝑑𝑜 𝑒𝑠𝑜𝑠 𝑑𝑎𝑡𝑜𝑠, 𝑝𝑟𝑜𝑐𝑒𝑑𝑒𝑚𝑜𝑠 𝑎 𝑠𝑢𝑠𝑡𝑖𝑡𝑢𝑖𝑟 𝑙𝑎𝑠 𝑓𝑢𝑛𝑐𝑖𝑜𝑛𝑒𝑠 𝑦1 , 𝑦2 , 𝑟(𝑥) 𝑦 𝑤(𝑥) 𝑒𝑛 𝑙𝑎 𝑟𝑒𝑙𝑎𝑐𝑖ó𝑛 𝑝𝑎𝑟𝑎
𝑦2 𝑟(𝑥) 𝑦 𝑟(𝑥)
𝑝𝑎𝑟𝑎 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑟 𝑙𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑝𝑎𝑟𝑡𝑖𝑐𝑢𝑙𝑎𝑟 𝑦𝑝 = −𝑦1 ∫ 𝑑𝑥 + 𝑦2 ∫ 1 𝑑𝑥
𝑤 𝑤
𝑒 𝑥 𝑐𝑜𝑠 𝑥 𝑒𝑥 (𝑒 𝑥 𝑐𝑜𝑠 𝑥)
𝑦𝑝 = −𝑒𝑥 ∫ 𝑑𝑥 + ∫ 𝑑𝑥 =
−𝑒 𝑥 −𝑒 𝑥

𝑦𝑝 = 𝑒𝑥 ∫ cos 𝑥 𝑑𝑥 − ∫ 𝑒 𝑥 𝑐𝑜𝑠 𝑥 𝑑𝑥 = 𝑒𝑥 (sin 𝑥) − ∫ 𝑒 𝑥 𝑐𝑜𝑠 𝑥 𝑑𝑥 =

∫ 𝑒 𝑥 𝑐𝑜𝑠 𝑥 𝑑𝑥 𝑆𝑖 𝑏𝑢𝑠𝑐𝑎𝑚𝑜𝑠 𝑒𝑛 𝑛𝑢𝑒𝑠𝑡𝑟𝑎 𝑡𝑎𝑏𝑙𝑎 𝑑𝑒 𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑙𝑒𝑠 𝑛𝑜𝑠 𝑑𝑎𝑚𝑜𝑠 𝑐𝑢𝑒𝑛𝑡𝑎 𝑞𝑢𝑒 𝑙𝑎 𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑙 𝑞𝑢𝑒𝑑𝑎 𝑎𝑠í:

1
∫ 𝑒 𝑥 𝑐𝑜𝑠 𝑥 𝑑𝑥 = 𝑒 𝑥 (𝑐𝑜𝑠𝑥 + 𝑠𝑖𝑛𝑥) 𝑠𝑢𝑠𝑡𝑖𝑡𝑢𝑦𝑒𝑛𝑑𝑜 𝑒𝑛 𝑙𝑎 𝑦𝑝
2
1 1
𝑦𝑝 = 𝑒𝑥 (sin 𝑥) − ( 𝑒 𝑥 (𝑐𝑜𝑠𝑥 + 𝑠𝑖𝑛𝑥)) = 𝑒𝑥 (sin 𝑥) − 𝑒 𝑥 (𝑐𝑜𝑠𝑥 + 𝑠𝑖𝑛𝑥) =
2 2

𝑒 𝑥 𝑐𝑜𝑠𝑥 𝑒 𝑥 𝑠𝑖𝑛𝑥 𝑐𝑜𝑠𝑥 𝑠𝑖𝑛𝑥 𝑒𝑥


𝑦𝑝 = 𝑒𝑥 (sin 𝑥) − − = 𝑒 𝑥 (𝑠𝑖𝑛𝑥 − − )= (𝑠𝑖𝑛𝑥 − cos 𝑥)
2 2 2 2 2
𝑒𝑥
𝑦𝑝 = (𝑠𝑖𝑛𝑥 − cos 𝑥)
2

3- Resolver mediante variación de parámetros la ecuación diferencial

𝒚′′ − 𝟖𝒚′ + 𝟏𝟔𝒚 = 𝒙𝒆𝟒𝒙

𝑃𝑟𝑖𝑚𝑒𝑟𝑜 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑚𝑜𝑠 𝑙𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 ℎ𝑜𝑚𝑜𝑔𝑒𝑛𝑒𝑎:

𝑚2 − 8𝑚 + 16 = 0 − −→ (𝑚 − 4)2 = 0 ∴ 𝑚1 = 𝑚2 = 4
𝑦ℎ = 𝑐1 𝑒 𝑥 + 𝑐2 𝑥𝑒 𝑥 = 𝑐1 𝑒 4𝑥 + 𝑐2 𝑥𝑒 4𝑥
𝑦1 = 𝑒 4𝑥 𝑦2 = 𝑥𝑒 4𝑥
𝑦1 𝑦2
𝐿𝑜 𝑠𝑖𝑔𝑢𝑖𝑒𝑛𝑡𝑒 𝑒𝑠 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑟 𝑒𝑙 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒𝑙 𝑊𝑟𝑜𝑛𝑠𝑘𝑖𝑎𝑛𝑜 𝑤(𝑦1 , 𝑦2 ) = | |
𝑦´1 𝑦´2
𝑦1 = 𝑒 4𝑥 𝑦2 = 𝑥𝑒 4𝑥 𝑦1′ = 4𝑒 4𝑥 𝑦2′ = 4𝑥𝑒 4𝑥 + 𝑒 4𝑥

𝑤(𝑦1 , 𝑦2 ) = | 𝑒4𝑥 𝑥𝑒4𝑥 |=


4𝑒 4𝑥 4𝑥𝑒4𝑥 + 𝑒4𝑥
= 4𝑥𝑒 8𝑥 + 𝑒 8𝑥 − 4𝑥𝑒 8𝑥 = 𝑒 8𝑥 ∴ 𝑤 = 𝑒 8𝑥
𝐸𝑙 𝑠𝑖𝑔𝑢𝑖𝑒𝑛𝑡𝑒 𝑝𝑎𝑠𝑜 𝑒𝑠 𝑖𝑑𝑒𝑛𝑡𝑖𝑓𝑖𝑐𝑎𝑟 𝑙𝑎 𝑓𝑢𝑛𝑐𝑖ó𝑛 𝑟(𝑥) = 𝑥𝑒 4𝑥
𝑂𝑏𝑡𝑒𝑛𝑖𝑒𝑛𝑑𝑜 𝑒𝑠𝑜𝑠 𝑑𝑎𝑡𝑜𝑠, 𝑝𝑟𝑜𝑐𝑒𝑑𝑒𝑚𝑜𝑠 𝑎 𝑠𝑢𝑠𝑡𝑖𝑡𝑢𝑖𝑟 𝑙𝑎𝑠 𝑓𝑢𝑛𝑐𝑖𝑜𝑛𝑒𝑠 𝑦1 , 𝑦2 , 𝑟(𝑥) 𝑦 𝑤(𝑥) 𝑒𝑛 𝑙𝑎 𝑟𝑒𝑙𝑎𝑐𝑖ó𝑛 𝑝𝑎𝑟𝑎
𝑦2 𝑟(𝑥) 𝑦 𝑟(𝑥)
𝑝𝑎𝑟𝑎 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑟 𝑙𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑝𝑎𝑟𝑡𝑖𝑐𝑢𝑙𝑎𝑟 𝑦𝑝 = −𝑦1 ∫ 𝑑𝑥 + 𝑦2 ∫ 1 𝑑𝑥
𝑤 𝑤

𝑥𝑒4𝑥 (𝑥𝑒 4𝑥 ) 𝑒4𝑥 (𝑥𝑒 4𝑥 ) 𝑥2 𝑒8𝑥 𝑥𝑒 8𝑥


𝑦𝑝 = −𝑒4𝑥 ∫ 𝑑𝑥 + 𝑥𝑒4𝑥 ∫ 𝑑𝑥 = −𝑒4𝑥 ∫ 𝑑𝑥 + 𝑥𝑒4𝑥 ∫ 8𝑥 𝑑𝑥 =
𝑒8𝑥 𝑒8𝑥 𝑒8𝑥 𝑒
𝑒 4𝑥 𝑥 3 𝑥𝑒 4𝑥 𝑥 2 𝑥3 𝑥3 𝑒 4𝑥 𝑥 3
𝑦𝑝 = −𝑒4𝑥 ∫ 𝑥 2 𝑑𝑥 + 𝑥𝑒4𝑥 ∫ 𝑥 𝑑𝑥 = − + = 𝑒 4𝑥 ( − ) =
3 2 2 3 6

𝑥 3 𝑒 4𝑥
𝑦𝑝 =
6
TRANSFORMADA DE LAPLACE

4- Resuelva la ecuación diferencial mediante la transformada de Laplace:

𝑦 ′′ − 2𝑦 ′ − 3𝑦 = 𝑒 𝑡 𝑦(0) = 2 𝑦 ′ (0) = 4

ℒ{𝑦 ′′ − 2𝑦 ′ − 3𝑦 = 𝑒 𝑡 } − −−→ ℒ{𝑦 ′′ } − ℒ{2𝑦 ′ } − ℒ{3𝑦} = ℒ{𝑒 𝑡 }


1
𝑠 2 ℒ{𝑦} − 𝑠𝑦(0) − 𝑦 ′(0) − 2(ℒ {𝑦} − 𝑦(0)) − 3ℒ{𝑦} =
𝑠−1
1
𝑠 2 ℒ{𝑦} − 𝑠(2) − 4 − 2𝑠ℒ{𝑦} + 2(2) − 3ℒ{𝑦} =
𝑠−1
1
𝑠 2 ℒ{𝑦} − 2𝑠 − 4 − 2𝑠ℒ{𝑦} + 4 − 3ℒ{𝑦} =
𝑠−1
1 2𝑠2 − 2𝑠 + 1
ℒ{𝑦}(𝑠 2 − 2𝑠 − 3) = + 2𝑠 − −→ ℒ{𝑦}(𝑠 2 − 2𝑠 − 3) =
𝑠−1 𝑠−1
2𝑠2 − 2𝑠 + 1 2𝑠2 − 2𝑠 + 1
ℒ{𝑦} = − −→ ℒ −1 (ℒ {𝑦} = )
𝑠 − 1(𝑠 − 3)(𝑠 + 1) 𝑠 − 1(𝑠 − 3)(𝑠 + 1)
2𝑠2 − 2𝑠 + 1 𝐴 𝐵 𝐶
𝑦 = ℒ −1 ( ) ∴ 𝑦 = ℒ −1 ( + + )
𝑠 − 1(𝑠 − 3)(𝑠 + 1) 𝑆−1 𝑆−3 𝑆+1
𝐴 𝐵 𝐶
𝑦 = ℒ−1 { −1
}+ℒ {
−1
}+ℒ { }
𝑆−1 𝑆−3 𝑆+1

𝑝𝑎𝑟𝑎 𝑒𝑛𝑐𝑜𝑛𝑡𝑟𝑎𝑟 𝑙𝑜𝑠 𝑣𝑎𝑙𝑜𝑟𝑒𝑠 𝑜𝑐𝑢𝑝𝑎𝑟𝑒𝑚𝑜𝑠 𝑓𝑟𝑎𝑐𝑐𝑖𝑜𝑛𝑒𝑠 𝑝𝑎𝑟𝑐𝑖𝑎𝑙𝑒𝑠:


2𝑠2 − 2𝑠 + 1 𝐴 𝐵 𝐶
= + +
𝑠 − 1(𝑠 − 3)(𝑠 + 1) 𝑆 − 1 𝑆 − 3 𝑆 + 1
2𝑠2 − 2𝑠 + 1 𝐴(𝑠 − 3)(𝑠 + 1) + 𝐵(𝑠 − 1)(𝑠 + 1)𝐶(𝑠 − 1)(𝑠 − 3)
=
𝑠 − 1(𝑠 − 3)(𝑠 + 1) 𝑠 − 1(𝑠 − 3)(𝑠 + 1)
2𝑠2 − 2𝑠 + 1 = 𝐴𝑠2 − 2𝐴𝑠 − 3𝐴 + 𝐵𝑠2 − 𝐵 + 𝐶𝑠2 − 4𝐶𝑠 + 3𝐶

𝐸𝑐1: 2𝑠 2 = 𝐴𝑠 2 + 𝐵𝑠 2 + 𝐶𝑠 2 − −→ 2𝑠 2 = 𝑠 2 (𝐴 + 𝐵 + 𝐶 ) − −→ 2 = 𝐴 + 𝐵 + 𝐶
{ 𝐸𝑐2: − 2𝑠 = −2𝐴𝑠 − 4𝐶𝑠 − −→ −2𝑠 = 𝑠(−2𝐴 − 4𝐶 ) − −→ 2 = 2𝐴 + 4𝐶
𝐸𝑐3: 1 = −3𝐴 − 𝐵 + 3𝐶
1 1 1 𝐴 2
|2 0 4| ∗ |𝐵 | = |2|
−3 −1 3 𝐶 1
Nos ayudaremos de un software (Octave) para facilitar el calculo de las incognitas, mediante la
formula: 𝑥 = 𝐴−1 ∗ 𝑏

𝑑𝑜𝑛𝑑𝑒 𝐴−1 𝑒𝑠 𝑙𝑎 𝑖𝑛𝑣𝑒𝑟𝑠𝑎 𝑑𝑒 𝑙𝑎 𝑚𝑎𝑡𝑟𝑖𝑧 𝑦 𝑏 𝑙𝑜𝑠 𝑟𝑒𝑠𝑢𝑙𝑡𝑎𝑑𝑜𝑠 𝑑𝑒 𝑐𝑎𝑑𝑎 𝑒𝑐𝑢𝑎𝑐𝑖ó𝑛


𝑆𝑢𝑠𝑡𝑖𝑡𝑢𝑦𝑒𝑛𝑑𝑜 𝑐𝑎𝑑𝑎 𝑣𝑎𝑙𝑜𝑟 𝑒𝑛 𝑛𝑢𝑒𝑠𝑡𝑟𝑎 𝑡𝑟𝑎𝑛𝑠𝑓𝑜𝑟𝑚𝑎𝑑𝑎, , 𝑛𝑜𝑠 𝑞𝑢𝑒𝑑𝑎:
1 13 5
−4
𝑦 = ℒ −1 ( ) + ℒ −1 ( 8 ) + ℒ −1 ( 8 )
𝑆−1 𝑆−3 𝑆+1

1 1 13 −1 1 5 1
𝑦𝑝 = − ℒ −1 ( )+ ℒ ( ) + ℒ −1 ( )
4 𝑆−1 8 𝑆−3 8 𝑆+1
1 13 3𝑡 5 −𝑡
𝑦𝑝 = − 𝑒 𝑡 + 𝑒 + 𝑒
4 8 8
5- Resuelva la siguiente ecuación diferencial por medio de la transformada de Laplace
5 1
𝑦 ′′ − 𝑦 ′ + 𝑦 = 0 𝑦(0) = 1 𝑦 ′ (0) =
2 2
5 5
ℒ (𝑦 ′′ − 𝑦 ′ + 𝑦 = 0) − −→ ℒ(𝑦 ′′ ) − ℒ(𝑦 ′ ) + ℒ(𝑦) = ℒ(0)
2 2
5 0
𝑠 2 ℒ {𝑦} − 𝑠𝑦(0) − 𝑦 ′ (0) − (𝑠ℒ{𝑦} − 𝑦(0)) + ℒ(𝑦) =
2 𝑠
1 5 5
𝑠 2 ℒ {𝑦} − 𝑠(1) − − 𝑠ℒ {𝑦} + + ℒ(𝑦) = 0
2 2 2
5 5
𝑠 2 ℒ {𝑦} − 𝑠ℒ {𝑦} + ℒ(𝑦) = 𝑠 − 2 − −→ ℒ{𝑦} (𝑠 2 − 𝑠 + 1) = 𝑠 − 2
2 2
5
Para calcular las raíces de la ecuacion: 𝑠 2 − 2 𝑠 + 1 crearemos un programa para facilitar el calculo de
sus raíces en base a la Formula General:

−𝑏 ± √𝑏 2 − 4𝑎𝑐
𝑥=
2𝑎
5
En nuestro problema el valor de 𝑎 = 1 𝑏 = −2 𝑐 = 1, a continuación se muestra una
captura de como se compone el programa al igual que el valor de las raíces:
𝐿𝑜 ú𝑛𝑖𝑐𝑜 𝑞𝑢𝑒 𝑞𝑢𝑒𝑑𝑎 𝑒𝑠 𝑑𝑒𝑠𝑝𝑒𝑗𝑎𝑟 𝑐𝑎𝑑𝑎 𝑣𝑎𝑙𝑜𝑟 𝑦 𝑡𝑒𝑛𝑒𝑚𝑜𝑠 𝑞𝑢𝑒:

5 1
𝑠 2 − 𝑠 + 1 = (𝑠 − 2) (𝑠 − ) 𝑠𝑢𝑠𝑡𝑖𝑡𝑢𝑦𝑒𝑛𝑑𝑜 𝑒𝑠𝑡𝑜 𝑒𝑛 𝑛𝑢𝑒𝑠𝑡𝑟𝑜 𝑝𝑟𝑜𝑏𝑙𝑒𝑚𝑎
2 2
1 (𝑠 − 2)
ℒ {𝑦}(𝑠 − 2) (𝑠 − ) = 𝑠 − 2 − −→ ℒ{𝑦} = 1
2 (𝑠 − 2) (𝑠 − )
2

(𝑠 − 2) 1
ℒ −1 (ℒ{𝑦} = 1
) − −→ 𝑦 = ℒ −1 ( 1
)
(𝑠 − 2) (𝑠 − ) (𝑠 − )
2 2
𝑡
𝑦𝑝 = 𝑒 2

6- Resuelva la siguiente ecuación diferencial promedio de la transformada de Laplace.

𝑦 ′′ − 4𝑦 = 2 𝑦(0) = 𝑦 ′ (0) = 0
′′
ℒ (𝑦′′ − 4𝑦 = 2) − −→ ℒ{𝑦 } − 4ℒ{𝑦} = ℒ{2}
2
𝑠 2 ℒ {𝑦} − 𝑠𝑦(0) − 𝑦 ′ (0) − 4ℒ{𝑦} =
𝑠
2 2 2
𝑠 2 ℒ {𝑦} − 4ℒ{𝑦} = − −→ ℒ {𝑦}(𝑠 2 − 4) = − −→ ℒ{𝑦} =
𝑠 𝑠 𝑠(𝑠 2 − 4)
2 1
ℒ −1 (ℒ{𝑦} = ) − −→ 𝑦 = 2ℒ −1 { }
𝑠(𝑠 2 − 4) 𝑠(𝑠 2 − 4)

Haciendo un análisis de estás formulas de transformadas de laplace, nos damos


cuenta que, la variación entre fórmula está en el signo + y -, por lo cual pordemos
decir que la transformada inversa de
1 1
ℒ −1 { 2 } − −−→ (cosh 𝑎𝑡 − 1)
𝑠(𝑠 − 𝑎2 ) 𝑎2
Aplcando esta fórmula, tenemos que
1
𝑦 = 2ℒ −1 { 2 }
𝑠(𝑠 − 4)
2 1
𝑦𝑝 = (cosh 2𝑡 − 1) ∴ 𝑦𝑝 = (cosh 2𝑡 − 1)
4 2

SERIES DE FOURIER

7- Hallar la serie de Fourier de la siguiente función.


1, −𝜋 < 𝑥 < 0
𝑓(𝑥) = {
−1, 0 < 𝑥 < 𝜋
Encontramos el primer coeficiente de Fourier

1 0 1 𝜋
𝑎0 = ∫ 1 𝑑𝑥 + ∫ −1 𝑑𝑥
2𝜋 −𝜋 2𝜋 0
1 1 1 1 1 1
𝑎0 = [𝑥]0−𝜋 − [𝑥]𝜋0 −→ [0 − (−𝜋)] − [𝜋 − 0]−→ [𝜋] − [𝜋] = 0 ∴ 𝒂𝟎 = 𝟎
2𝜋 2𝜋 2𝜋 2𝜋 2𝜋 2𝜋
Encontramos el segundo coeficiente de Fourier

1 0 1 𝜋 1 0 1 𝜋
𝑎𝑛 = ∫ cos(𝑛𝑥) 𝑑𝑥 + ∫ − cos(𝑛𝑥) 𝑑𝑥 − −−> ∫ cos(𝑛𝑥) 𝑑𝑥 − ∫ cos(𝑛𝑥) 𝑑𝑥
𝜋 −𝜋 𝜋 0 𝜋 −𝜋 𝜋 0
0 𝜋
1 1 1 1 1 1
[ sin(𝑛𝑥)] − [ sin(𝑛𝑥)] −→ [sin(0) − sin(−𝑛𝜋)] − [sin(𝑛𝜋) − sin(0)] = 0
𝜋 𝑛 −𝜋 𝜋 𝑛 0 𝑛𝜋 𝑛𝜋
∴ 𝒂𝒏 = 𝟎
Encontramos el tercer coeficiente de Fourier

1 0 1 𝜋 1 0 1 𝜋
𝑏𝑛 = ∫ sin(𝑛𝑥) 𝑑𝑥 + ∫ − sin(𝑛𝑥) 𝑑𝑥 − −→ ∫ sin(𝑛𝑥) 𝑑𝑥 − ∫ sin(𝑛𝑥) 𝑑𝑥
𝜋 −𝜋 𝜋 0 𝜋 −𝜋 𝜋 0
0 𝜋
1 1 1 1 1 1
[− cos(𝑛𝑥)] − [− cos 𝑛𝑥] − −→ − [cos(𝑛𝑥)]0−𝜋 + [cos(𝑛𝑥)]𝜋0
𝜋 𝑛 −𝜋 𝜋 𝑛 0 𝑛𝜋 𝑛𝜋
1 1 1 1
− [cos(0) − cos(−𝑛𝜋)] + [cos(𝑛𝜋) − cos(0)]−→ − [1 − (−1)𝑛 ] + [(−1)𝑛 − 1]
𝑛𝜋 𝑛𝜋 𝑛𝜋 𝑛𝜋
1 2 𝟐[(−𝟏)𝒏 − 𝟏]
[−1 + (−1)𝑛 + (−1)𝑛 − 1] − −→ ((−1)𝑛 − 1) ∴ 𝒃𝒏 =
𝑛𝜋 𝑛𝜋 𝒏𝝅
Sustituimos en la serie de Fourier los valores que acabamos de encontrar.

𝑓(𝑥) = 𝑎0 + ∑(𝑎𝑛 cos 𝑛𝑥 + 𝑏𝑛 sin 𝑛𝑥)


𝑛=1

Donde:
2[(−1)𝑛 − 1]
𝑎0 = 0; 𝑎𝑛 = 0; 𝑏𝑛 =
𝑛𝜋

2[(−1)𝑛 − 1]
𝑓(𝑥) = 0 + ∑ (0 cos 𝑛𝑥 + sin 𝑛𝑥)
𝑛𝜋
𝑛=1
∞ ∞
2[(−1)𝑛 − 1] 2 (−1)𝑛 − 1
𝑓(𝑥) = ∑ ( sin 𝑛𝑥) − −→ 𝑓(𝑥) = ∑ (sin 𝑛𝑥)
𝑛𝜋 π 𝑛
𝑛=1 𝑛=1
Evaluamos los 10 primeros términos de la serie
2 2 2 2 2
𝑓(𝑥) = (−2 sin 𝑥 − sin 3𝑥 − sin 5𝑥 − sin 7𝑥 − sin 9𝑥)
π 3 5 7 9
4 1 1 1 1 1 1
𝑓(𝑥) = − (sin 𝑥 + sin 3𝑥 + sin 5𝑥 + sin 7𝑥 + sin 9𝑥 + sin 11𝑥 + sin 13𝑥
π 3 5 7 9 11 13
1 1 1
+ sin 15𝑥 + sin 17𝑥 + sin 19𝑥)
15 17 19

Ocupando un graficador obtenemos lo siguiente:

8- Hallar la serie de Fourier de la siguiente función:


0, −𝜋 < 𝑥 < 0
𝜋
−1, 0 < 𝑥 <
𝑓(𝑥) = 2
𝜋
1, <𝑥<𝜋
{ 2
Encontramos el primer coeficiente de Fourier
𝜋
𝜋
1 0 1 2 1 𝜋 1 1 1
𝑎0 = ∫ 0 𝑑𝑥 + ∫ −1 𝑑𝑥 + ∫ 1 𝑑𝑥 − −→ [𝐶]0−𝜋 − [𝑥]02 + [𝑥]𝜋𝜋
2𝜋 −𝜋 2𝜋 0 2𝜋 𝜋 2𝜋 2𝜋 2𝜋 2
2

1 𝜋 1 𝜋 1 𝜋 1 𝜋 1 1
− [ − 0] + [𝜋 − ] − −→ − [ ]+ [ ] − −→ − + = 0 ∴ 𝒂𝟎 = 𝟎
2𝜋 2 2𝜋 2 2𝜋 2 2𝜋 2 4 4
Encontramos el segundo coeficiente de Fourier
𝜋
1 0 1 2 1 𝜋
𝑎𝑛 = ∫ 0 cos(𝑛𝑥) 𝑑𝑥 + ∫ − cos(𝑛𝑥) 𝑑𝑥 + ∫ cos(𝑛𝑥) 𝑑𝑥
𝜋 −𝜋 𝜋 0 𝜋 𝜋
2
𝜋 𝜋
𝜋
1 2 1 𝜋 1 1 2 1 1
𝑎𝑛 = − ∫ cos(𝑛𝑥) 𝑑𝑥 + ∫ cos(𝑛𝑥) 𝑑𝑥 − −→ − [ sin(𝑛𝑥)] + [ sin(𝑛𝑥)]𝜋
𝜋 0 𝜋 𝜋 𝜋 𝑛 0 𝜋 𝑛
2 2

1 𝑛𝜋 1 𝑛𝜋 1 𝑛+3 1 𝑛+3
− [sin ( ) − sin(0)] + [sin(𝑛𝜋) − sin ( )] − −→ − [(−1) 2 ] + [−(−1) 2 ]
𝑛𝜋 2 𝑛𝜋 2 𝑛𝜋 𝑛𝜋
𝑛+3 𝒏+𝟑
2 𝑛+3 2 (−1) 2 𝟐 (−𝟏) 𝟐
− [(−1) 2 ] − −→ − [ ] ∴ 𝒂𝒏 = − [ ]
𝑛𝜋 𝜋 𝑛 𝝅 𝒏

𝑛+3
𝑓(𝑥) = {(−1) 2 𝑛 = 𝑖𝑚𝑝𝑎𝑟
0 𝑛 = 𝑝𝑎𝑟
Encontramos el tercer coeficiente de Fourier
𝜋
1 0 1 2 1 𝜋
𝑏𝑛 = ∫ 0 sin(𝑛𝑥) 𝑑𝑥 − ∫ sin(𝑛𝑥) 𝑑𝑥 + ∫ sin(𝑛𝑥) 𝑑𝑥
𝜋 −𝜋 𝜋 0 𝜋 𝜋
2
𝝅
𝝅 𝝅
1 1 1 𝟐 1 1 1 1
[𝐶] − [− cos(𝑛𝑥)] + [− cos(𝑛𝑥)]𝝅 − −→ [cos(𝑛𝑥)]𝟎𝟐 − [cos(𝑛𝑥)]𝝅𝝅
𝜋 𝜋 𝑛 𝟎 𝜋 𝑛 𝑛𝜋 𝑛𝜋 𝟐
𝟐

1 𝑛𝜋 1 𝑛𝜋 1 𝑛 1 𝑛
[cos ( ) − cos(0)] − [cos(𝑛𝜋) − cos ( )] − −→ [(−1) 2 − 1] − [1 − (−1) 2 ]
𝑛𝜋 2 𝑛𝜋 2 𝑛𝜋 𝑛𝜋
𝒏
𝟐 (−𝟏)𝟐 − 𝟏
∴ 𝒃𝒏 = [ ]
𝝅 𝒏

𝑛
𝑓(𝑥) = {(−1) 2 𝑛 = 𝑝𝑎𝑟
0 𝑛 = 𝑖𝑚𝑝𝑎𝑟
Sustituimos en la serie de Fourier los valores que acabamos de encontrar.

𝑓(𝑥) = 𝑎0 + ∑(𝑎𝑛 cos 𝑛𝑥 + 𝑏𝑛 sin 𝑛𝑥)


𝑛=1

Donde:
𝑛+3 𝑛
2 (−1) 2 2 (−1) 2 − 1
𝑎0 = 0; 𝑎𝑛 = − [ ] ; 𝑏𝑛 = [ ]
𝜋 𝑛 𝜋 𝑛

𝑛+3 𝑛
𝑃𝑎𝑟𝑎 𝑎𝑛 − −→ 𝑓(𝑥) = {(−1) 2 𝑛 = 𝑖𝑚𝑝𝑎𝑟 𝑃𝑎𝑟𝑎 𝑏𝑛 − −→ 𝑓(𝑥) = {(−1) 2 𝑛 = 𝑝𝑎𝑟
0 𝑛 = 𝑝𝑎𝑟 0 𝑛 = 𝑖𝑚𝑝𝑎𝑟
∞ 𝑛+3 𝑛
2 (−1) 2 2 (−1) 2 − 1
𝑓(𝑥) = ∑ (− [ ] cos 𝑛𝑥 + [ ] sin 𝑛𝑥)
𝜋 𝑛 𝜋 𝑛
𝑛=1

∞ 𝑛+3 𝑛
2 (−1) 2 (−1) 2 − 1
𝑓(𝑥) = − ∑ ([ ] cos 𝑛𝑥 − [ ] sin 𝑛𝑥)
π 𝑛 𝑛
𝑛=1

Evaluamos los 11 primeros términos de la serie


2 1 1 1 1 1
𝑓(𝑥) = − (cos(𝑥) + sin(2𝑥) − cos 3𝑥 + cos 5𝑥 + sin 6𝑥 − cos 7𝑥 + cos 9𝑥
π 3 5 3 7 9
1 1 1 1
+ sin 10𝑥 − cos 11𝑥 + cos 13𝑥 + sin 14𝑥)
5 11 13 7
Ocupando un graficador obtenemos lo siguiente, muy parecida a la original:
9- Resolver la siguiente ecuación diferencial con series de Fourier

𝑥2
𝑓(𝑥) = +3
6

Encontramos el primer coeficiente de la serie de Fourier


𝜋
1 𝜋 𝑥2 1 𝑥3 1 𝜋 3 (−𝜋)3
𝑎0 = ∫ ( + 3) 𝑑𝑥 − −→ [ + 3𝑥] − −→ [ − + 3𝜋 − 3(−𝜋)]
2𝜋 −𝜋 6 2𝜋 18 −𝜋
2𝜋 18 18

1 𝜋3 𝜋3 1 2𝜋 3 2𝜋 𝜋 2 𝝅𝟐
[ + + 3𝜋 + 3𝜋] − −→ [ + 6𝜋] − −→ [ + 3] ∴ 𝒂𝟎 = +𝟑
2𝜋 18 18 2𝜋 18 2𝜋 18 𝟏𝟖

Calculamos el segundo coeficiente de Fourier

1 𝜋 𝑥2 1 𝜋 𝑥2 1 𝜋
𝑎𝑛 = ∫ ( + 3) cos(𝑛𝑥) 𝑑𝑥 − −→ ∫ cos(𝑛𝑥) 𝑑𝑥 + ∫ 3 cos(𝑛𝑥) 𝑑𝑥
𝜋 −𝜋 6 𝜋 −𝜋 6 𝜋 −𝜋

1 𝜋 2 3 𝜋
∫ 𝑥 cos(𝑛𝑥) 𝑑𝑥 + ∫ cos(𝑛𝑥) 𝑑𝑥
6𝜋 −𝜋 𝜋 −𝜋
𝜋 𝜋
1 2𝑥 𝑥2 2 3 1
[ cos(𝑛𝑥) + ( − 3 ) sin(𝑛𝑥)] + [ sin(𝑛𝑥)]
6𝜋 𝑛2 𝑛 𝑛 −𝜋
𝜋 𝑛 −𝜋

Evaluamos los términos por separado para hacerlo más sencillo a la vista

1- Evaluamos el primer término.

1 2𝜋 −2𝜋 𝜋2 2 (−𝜋)2 2
[ 2 cos(𝑛𝜋) − 2 cos(−𝑛𝜋) + ( − 3 ) sin(𝑛𝜋) − ( − 3 ) sin(−𝑛𝜋)]
6𝜋 𝑛 𝑛 𝑛 𝑛 𝑛 𝑛

1 2𝜋 2𝜋 2 2(−1)𝑛
[ 2 cos(𝑛𝜋) + 2 cos(−𝑛𝜋)] − −→ 2 cos(𝑛𝜋) − −→
6𝜋 𝑛 𝑛 3𝑛 3𝑛2
2- Evaluamos el segundo término.
3
[sin(𝑛𝜋) − sin(−𝑛𝜋)] = 0
𝑛𝜋
Sumamos los resultados que acabamos de obtener
2(−1)𝑛 𝟐(−𝟏)𝒏
𝑎𝑛 = + 0 ∴ 𝒂𝒏 =
3𝑛2 𝟑𝒏𝟐
Calculamos el tercer coeficiente de Fourier

1 𝜋 𝑥2 1 𝜋 𝑥2
𝑏𝑛 = ∫ ( + 3) sin(𝑛𝑥) 𝑑𝑥 − −→ ∫ ( sin(𝑛𝑥) + 3 sin(𝑛𝑥)) 𝑑𝑥
𝜋 −𝜋 6 𝜋 −𝜋 6

1 𝜋 𝑥2 1 𝜋 1 𝜋 2 3 𝜋
∫ sin(𝑛𝑥) 𝑑𝑥 + ∫ 3 sin(𝑛𝑥) 𝑑𝑥 − −→ ∫ 𝑥 sin(𝑛𝑥) 𝑑𝑥 + ∫ sin(𝑛𝑥) 𝑑𝑥
𝜋 −𝜋 6 𝜋 −𝜋 6𝜋 −𝜋 𝜋 −𝜋
𝜋 𝜋
1 2𝑥 2 𝑥2 3 1
[ 2 sin(𝑛𝑥) + ( 3 − ) cos(𝑛𝑥)] + [− cos(𝑛𝑥)]
6𝜋 𝑛 𝑛 𝑛 −𝜋
𝜋 𝑛 −𝜋

1 2 𝜋2 2 𝜋2 3
[( 3 − ) cos(𝑛𝜋) − ( 3 − ) cos(−𝑛𝜋)] − [cos(𝑛𝜋) − cos(−𝑛𝜋)]
6𝜋 𝑛 𝑛 𝑛 𝑛 𝑛𝜋

1 2 𝜋2 2 𝜋2
[( 3 − ) cos(𝑛𝜋) − ( 3 − ) cos(𝑛𝜋)] = 0 ∴ 𝒃𝒏 = 𝟎
6𝜋 𝑛 𝑛 𝑛 𝑛

Sustituimos en la serie de Fourier los valores que acabamos de encontrar.


𝑓(𝑥) = 𝑎0 + ∑(𝑎𝑛 cos 𝑛𝑥 + 𝑏𝑛 sin 𝑛𝑥)


𝑛=1

Donde:

𝜋2 2(−1)𝑛
𝑎0 = + 3; 𝑎𝑛 = ; 𝑏𝑛 = 0
18 3𝑛2

∞ ∞
𝜋2 2(−1)𝑛 𝜋2 2(−1)𝑛
𝑓(𝑥) = +3+ ∑( cos 𝑛𝑥 + 0 sin 𝑛𝑥) − −→ + 3 + ∑ ( cos 𝑛𝑥)
18 3𝑛2 18 3𝑛2
𝑛=1 𝑛=1

𝜋2 2 (−1)𝑛
𝑓(𝑥) = + 3 + ∑ ( 2 cos 𝑛𝑥)
18 3 𝑛
𝑛=1

Evaluamos los 10 primeros términos de la serie

𝜋2 2 1 1 1 1 1 1
+ 3 + (− cos 𝑥 + cos 2𝑥 − cos 3𝑥 + cos 4𝑥 − cos 5𝑥 + cos 6𝑥 − cos 7𝑥
18 3 4 9 16 25 36 49
1 1 1
+ cos 8𝑥 − cos 9𝑥 + cos 10𝑥)
64 81 100
Ocupando graficador obtenemos lo siguiente:

También podría gustarte