MAPIQ

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

MAPIQ 2

Solución TAREA #3
1) Encuentra los eigenvalores y eigenfunciones para las siguientes ecuaciones con las
condiciones frontera que se marcan
b) Identifica la función de peso 𝑝(𝑥 ) = 0
𝐿
c) Cual sería la condición de ortogonalidad ∫0 𝑝(𝑥 ) 𝑦𝑛 (𝑥 )𝑦𝑚 (𝑥) = 0 ∴ 𝑛 ≠ 𝑚

Respuesta:
𝒊) 𝒚′′ − 𝝀𝒚 = 𝟎; 𝒚(𝟎) + 𝒚′ (𝟎) = 𝟎, 𝒚(𝟐) = 𝟎
Sustituyendo por: 𝑦 = 𝑒 𝑚𝑥

𝑒 𝑚𝑥 (𝑚2 + 𝜆) = 0 ∴ 𝑚1,2 = ±√−𝜆


Sustituyendo por: 𝜆 = −𝛼 2 (Solución imaginaria)

𝑚1,2 = ±�−𝛼 2 = ±𝛼𝑖


𝑦 = 𝐶1 cos 𝛼𝑥 + 𝐶2 𝑠𝑒𝑛 𝛼𝑥

Sustituyendo la primera condición frontera:


𝐶1 cos 0 + 𝐶2 𝑠𝑒𝑛 0 − 𝛼𝐶1 𝑠𝑒𝑛 0 + 𝛼𝐶2 𝑐𝑜𝑠 0 = 0
𝐶1 + 𝛼𝐶2 = 0 ∴ 𝐶1 = −𝛼𝐶2
Sustituyendo la segunda condición frontera:
𝐶1 cos 2𝛼 + 𝐶2 𝑠𝑒𝑛 2𝛼 = 0 ∴ −𝛼𝐶2 cos 2𝛼 + 𝐶2 𝑠𝑒𝑛 2𝛼 = 0
𝐶2 = 0 𝑁𝑜 𝑒𝑠 𝑝𝑜𝑠𝑖𝑏𝑙𝑒
𝐶2 [− 𝛼cos 2𝛼 + 𝑠𝑒𝑛 2𝛼 ] = 0 �
− 𝛼cos 2𝛼 + 𝑠𝑒𝑛 2𝛼 = 0 𝑆𝑖 𝑒𝑠 𝑝𝑜𝑠𝑖𝑏𝑙𝑒
Cambiando la variable 𝛼 𝑝𝑜𝑟 𝑥 y graficando
-x*cos(2*x)+sin(2*x)
20

18

16

14

12

Valores "y"
10

0
0 2* 3* 4* 5* 6*
Valores "x"

Los cruces de la gráfica con cero son aproximadamente:


𝛼1 ~2.13; 𝛼2 ~3.79; 𝛼3 ~5.40; 𝛼4 ~6.99; 𝛼5 ~8.58 … ..
Por lo tanto los eigenvalores son:
𝜆𝑛 = (𝛼𝑛 )2 ; 𝑑𝑒 𝑙𝑎 𝑔𝑟á𝑓𝑖𝑐𝑎: (− 𝛼cos 2𝛼 + 𝑠𝑒𝑛 2𝛼 = 0)
Por lo tanto las eigenfunciones son:
𝑦 = 𝐶2 (−𝛼 cos 𝛼𝑥 + 𝑠𝑒𝑛 𝛼𝑥)
Haciendo un conjunto de soluciones

𝑦 = � 𝐶𝑛 (−𝛼𝑛 cos 𝛼𝑛 𝑥 + 𝑠𝑒𝑛 𝛼𝑛 𝑥)


𝑛=1

𝑦𝑛 = −𝛼𝑛 cos 𝛼𝑛 𝑥 + 𝑠𝑒𝑛 𝛼𝑛 𝑥


b) La autoadjunta
𝑦 ′′ − 𝜆𝑦 = 0
𝑑 ′
(𝑦 ) − 𝜆𝑦 = 0
𝑑𝑥
𝑟 (𝑥 ) = 1
𝑞 (𝑥 ) = 0
𝑝(𝑥 ) = −1
c) La condición de ortogonalidad:
1
� −1 [−𝛼𝑛 cos 𝛼𝑛 𝑥 + 𝑠𝑒𝑛 𝛼𝑛 𝑥 ][−𝛼𝑚 cos 𝛼𝑚 𝑥 + 𝑠𝑒𝑛 𝛼𝑚 𝑥 ] = 0 ∴ 𝑛 ≠ 𝑚
0
𝒊𝒊) 𝒙𝟐 𝒚′′ + 𝒙𝒚′ − 𝝀𝒚 = 𝟎; 𝒚(𝟏) = 𝟎, 𝒚(𝟐𝝅) = 𝟎

Sustituyendo por: 𝑦 = 𝑥 𝑚
𝑎 = 1; 𝑏 = 1; 𝑐 = −𝜆
La fórmula de Cauchy-Euler nos dice
𝑥 𝑚 (𝑎𝑚2 + (𝑏 − 𝑎)𝑚 + 𝑐 ) = 0
𝑥 𝑚 (𝑚2 − 𝜆) = 0 ∴ 𝑚1,2 = ±√𝜆
Sustituyendo por: 𝜆 = −𝛼 2 (Solución imaginaria)

𝑚1,2 = ±�−𝛼 2 = ±𝛼𝑖


𝑦 = 𝐶1 cos(𝛼 ln 𝑥) + 𝐶2 𝑠𝑒𝑛 (𝛼 ln 𝑥)

Sustituyendo la primera condición frontera:


𝐶1 cos 𝛼(ln 1) + 𝐶2 𝑠𝑒𝑛 𝛼 (ln 1) = 0
𝐶1 cos 0 + 𝐶2 𝑠𝑒𝑛 0 = 0 ∴ 𝐶1 = 0

Sustituyendo la segunda condición frontera:


𝐶2 𝑠𝑒𝑛 𝛼 (ln 2𝜋) = 0
𝐶 =0 𝑁𝑜 𝑒𝑠 𝑝𝑜𝑠𝑖𝑏𝑙𝑒
𝐶2 [𝑠𝑒𝑛 𝛼 (ln 2𝜋)] = 0 � 2
𝑠𝑒𝑛 𝛼 (ln 2𝜋) 𝑆𝑖 𝑒𝑠 𝑝𝑜𝑠𝑖𝑏𝑙𝑒
𝑛𝜋
𝑠𝑒𝑛 𝛼 (ln 2𝜋) = 0 = 𝑠𝑒𝑛 (𝑛𝜋) ∴ 𝛼𝑛 =
ln 2𝜋

De la gráfica: 𝑠𝑒𝑛 (𝑥 ln 2𝜋) = 0


1𝜋
= 1.71
ln 2𝜋
2𝜋
= 3.41
ln 2𝜋
3𝜋
= 5.13
ln 2𝜋
Por lo tanto los eigenvalores son:
𝑛2 𝜋 2
𝜆𝑛 = (𝛼𝑛 )2 =
ln( 2𝜋)2
Por lo tanto las eigenfunciones son:
𝑛𝜋
𝑦 = 𝐶2 𝑠𝑒𝑛 (ln 𝑥)
ln 2𝜋
Haciendo un conjunto de soluciones:

𝑛𝜋
𝑦 = � 𝐶𝑛 𝑠𝑒𝑛 (ln 𝑥)
ln 2𝜋
𝑛=1
𝑛𝜋
𝑦𝑛 = 𝑠𝑒𝑛 (ln 𝑥)
ln 2𝜋
b) La autoadjunta
𝑥 2 𝑦 ′′ + 𝑥𝑦 ′ − 𝜆𝑦 = 0
Factor integrante

𝜇 (𝑥 ) = 𝑒 ∫ 𝑃(𝑥)𝑑𝑥
Dejando la ecuación de la forma general:
𝑦 ′′ + 𝑃(𝑥)𝑦 ′ + 𝑄(𝑥)𝑦 = 0
1 𝜆
𝑦 ′′ + 𝑦 ′ − 2 𝑦 = 0
𝑥 𝑥
1
𝜇 (𝑥 ) = 𝑒 ∫𝑥𝑑𝑥 = 𝑒 ln 𝑥 = 𝑥
Multiplicando la ecuación en naranja por el factor integrante:
1 𝜆
𝑥 �𝑦 ′′ + 𝑦 ′ − 2 𝑦 = 0�
𝑥 𝑥
𝜆
𝑥𝑦 ′′ + 𝑦 ′ − 𝑦 = 0
𝑥
Es lo mismo que escribir:
𝑑 𝜆
[𝑥𝑦 ′ ] − 𝑦 = 0
𝑑𝑥 𝑥
𝑟 (𝑥 ) = 𝑥
𝑞 (𝑥 ) = 0
1
𝑝 (𝑥 ) = −
𝑥
c) La condición de ortogonalidad
1
1 𝑛𝜋 𝑚𝜋
� − �𝑠𝑒𝑛 (ln 𝑥)� �𝑠𝑒𝑛 (ln 𝑥)� = 0 ∴ 𝑛 ≠ 𝑚
0 𝑥 ln 2𝜋 ln 2𝜋

𝒊𝒊𝒊) 𝒚′′ − 𝒚′ + 𝝀𝒚 = 𝟎; 𝒚(𝟎) = 𝟎, 𝒚(𝟏𝟎) = 𝟎

Sustituyendo por: 𝑦 = 𝑒 𝑚𝑥
1 √1 − 4𝜆
𝑒 𝑚𝑥 (𝑚2 − 𝑚 + 𝜆) = 0 ∴ 𝑚1,2 = ±
2 2
1 − 4𝜆 = −4𝛼 2
4𝛼 2 + 1
𝜆=
4
4𝛼 2 +1
Sustituyendo por: 𝜆 =
4

1 √−4𝛼 2 1
𝑚1,2 = ± = ± 𝛼𝑖
2 2 2
1
𝑦 = 𝑒 2𝑥 [𝐶1 cos 𝛼𝑥 + 𝐶2 𝑠𝑒𝑛 𝛼𝑥 ]

Sustituyendo la primera condición frontera:


𝑒 0 [𝐶1 cos 0 + 𝐶2 𝑠𝑒𝑛 0] = 0 ∴ 𝐶1 = 0

Sustituyendo la segunda condición frontera:


1
𝑒 2𝜋 𝐶2 [𝑠𝑒𝑛 10𝛼 ] = 0
1 1
𝜋
𝑒 2𝜋 𝐶2 [𝑠𝑒𝑛 10𝛼 ] = 0 �𝐶2 = 0; 𝑒 2 𝑁𝑜 𝑒𝑠 𝑝𝑜𝑠𝑖𝑏𝑙𝑒
𝑠𝑒𝑛 10𝛼 𝑆𝑖 𝑒𝑠 𝑝𝑜𝑠𝑖𝑏𝑙𝑒
𝑛𝜋
𝑠𝑒𝑛 (10𝛼 ) = 0 = 𝑠𝑒𝑛 (𝑛𝜋) ∴ 𝛼𝑛 =
10
Por lo tanto los eigenvalores son:
𝜆𝑛 = (𝛼𝑛 )2 ; 𝜆𝑛 = 𝑛2
Por lo tanto las eigenfunciones son:
1 𝑛𝜋
𝑦 = 𝑒 2𝑥 �𝐶2 𝑠𝑒𝑛 𝑥�
10
Haciendo un conjunto de soluciones

1 𝑛𝜋
𝑦 = � 𝐶𝑛 𝑒 2𝑥 (𝑠𝑒𝑛 𝑥)
10
𝑛=1
1 𝑛𝜋
𝑦𝑛 = 𝑒 2𝑥 (𝑠𝑒𝑛 𝑥)
10
b) La autoadjunta de la ecuación
𝑦 ′′ + 𝑦 ′ + 𝜆𝑦 = 0
Factor integrante
𝜇 (𝑥 ) = 𝑒 ∫ 𝑃(𝑥)𝑑𝑥

Dejando la ecuación de la forma general:


𝑦 ′′ + 𝑃(𝑥)𝑦 ′ + 𝑄(𝑥)𝑦 = 0
𝑦 ′′ − 𝑦 ′ + 𝜆𝑦 = 0
𝜇 (𝑥 ) = 𝑒 ∫ −1𝑑𝑥 = 𝑒 −𝑥
Multiplicando la ecuación en naranja por el factor integrante:
𝑒 −𝑥 [𝑦 ′′ + 𝑦 ′ + 𝜆𝑦 = 0]
𝑒 −𝑥 𝑦 ′′ − 𝑒 −𝑥 𝑦 ′ + 𝑒 −𝑥 𝜆𝑦 = 0
Es lo mismo que escribir:
𝑑 −𝑥 ′
[𝑒 𝑦 ] + 𝑒 −𝑥 𝜆𝑦 = 0
𝑑𝑥
𝑟(𝑥 ) = 𝑒 −𝑥
𝑞 (𝑥 ) = 0
𝑝(𝑥 ) = 𝑒 −𝑥
c) La condición de ortogonalidad
1 1 1
� 𝑒 −𝑥 �𝑒 2𝑥 (𝑠𝑒𝑛 𝑛𝑥)� �𝑒 2𝑥 (𝑠𝑒𝑛 𝑚𝑥)� = 0 ∴ 𝑛 ≠ 𝑚
0

También podría gustarte