Ejemplo Frobenius

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

UNIVERSIDAD PEDAGÓGICA NACIONAL

FRANCISCO MORAZÁN
CENTRO UNIVERSITARIO DE EDUCACIÓN A DISTANCIA
ESPACIO DE ECUACIONES DIFERENCIALES
DE LOS DOCUMENTOS DE MAURICIO MONCADA

EJERCICIOS VARIOS SERIES DE POTENCIAS

Ejemplo 1:

 k! 
Sea la serie 2
k0
k!
k
x k , donde a k     k x k  . Utilice el criterio de la razón para
2 
determinar el radio de convergencia de la serie.


 k! 
Sea la serie
k0
2
k!
k
2
 
x k , donde a k   k x k  . Utilice el criterio de la razón para

determinar el radio de convergencia de la serie.


k
 
Teorema: Si lim a k existe, entonces a
k0
k
es convergente en un radio R.

 L  1  a k converge

ak  1 ak  1  L  1  a k diverge
Demostración: Criterio de la razón: lim ,si lim 
 L  1  a k no se puede
k ak k ak
 concluir

 k  1 ! x k 1
 k  1 ! x k
x
k
lim 2 k 1  lim
2 2
k k! k k k! x k
x
2k 2k

2 k  k  1 ! x k x
 lim
k
2 k 2 k! x k

 lim
 k  1 ! x , como  k  1 !   k  1 k !,resulta :
k 2 k!

 lim
 k  1 k ! x
k
2 k!
x
 lim k  1 , siendo que k  entonces k  1  0, luego:
2 k
x
 lim k  1, evaluando se tiene:
2 k
0  1, si x  0

        , si x  0
x
2   , si x  0


k!
La serie 2 k
x k es divergente para todo x   0 , es decir, únicamente converge
k0

cuando x  0.

Ejemplo 2:
Resolver por series de potencias  x  1 y   2  x  y  y  0, y  0   2, y  0    1

  
y   c n x n , y   nc n x n  1 , y   n  n  1 c n x n  2
n0 n 1 n2

 x  1 y   2  x  y  y  0
    
  n  n  1 c n x n  2  1   n  n  1 c n x n  2  2  nc n x n  1   nc n x n  1  1   c n x n  0
n2 n2 n 1 n 1 n0
    
  n  n  1 c n x n  1   n  n  1 c n x n  2  2  nc n x n  1   nc n x n   c n x n  0
n2 n2 n 1 n 1 n0

    
  k  k  1 c k  1x k    k  2  k  1 c k  2 x k  2   k  1 c k  1x k   k c k x k   c k x k  0
k 1 k 0 k 0 k 1 k 0

Para k  0, donde se requiera, se tiene:


  
  k  k  1 c k  1x k  2 c 2    k  2  k  1 c k  2 x k  2 c 1  2   k  1 c k  1x k 
k 1 k 1 k 1
 
  k c kx k  c 0   c kx k  0
k 1 k 1

 2 c 2  2 c 1  c 0    k  k  1 c k  1   k  2  k  1 c k  2  2  k  1 c k  1  k c k  c k  x k  0
k 1

2c1  c0
i) 2 c 2  2 c 1  c 0  0  c 2 
2

ii ) k  k  1 c k  1   k  2  k  1 c k  2  2  k  1 c k  1  k c k  c k  c k  2 
2  kc k 1
ck
, k  1, 2, 3,...
k2
Evaluar para k:
2c1  c0
c 2  c1  c1 c0
k  1 c3   2  
3 3 6
2c1  c0
c2 2 c0  2c1
k  2  c4   
4 4 8
 c0  2c1   c0  c  2c c0  3 c 0  6 c 1  4c 0
      0 1

c4  c3  8 6
k  3  c5        8 6  24
5 5 5 5
6c1  c0

120
 6 c 1  c 0   c 0  2 c 1   12 c  2c  15c  30 c
2 
 120    1 0 0 1
 2 c 5  c4    8   120 18 c 1  17 c 0
k  4  c6   
6 6 6 720

Solución para los valores arbitrarios: a) c 0  1, c 1  0 b) c 0  0, c 1  1


1 1 1 1 17
a) c 2   c3   c4  c5  c6  
2 6 8 120 720
y1  x   1 
1 2 1 3 1 4 1 5 17 6
x  x  x  x  x  ...
2 6 8 120 720

b) Siguiendo el mismo procedimiento anterior, resulta:

y1  x   x  x 2  x 4 
1 1 5 1 6
x  x  ...
4 20 40
c) Solución general:
 
y  x  C 1  1  x 2  x 3  x 4 
1 1 1 1 5 17 6
x  x  ...  
 2 6 8 120 720 
 1 1 5 1 6 
C2  x  x 2  x 4  x  x  ... 
 4 20 40 
Resolver el problema de valor inicial: y  0   2, y  0    1
 
y  x   C 1   x  x 2  x 3  x 4 
1 1 1 17 5
x  ...  
 2 2 24 120 
 1 3 
 C 2  1  2 x  x 3  x 4  x 5  ... 
 4 20 
Sistema:
C 1  1  C 2  0   2  C 1  2

C 1  0   C 2  1   1  C 2   1
Un miembro particular de la familia de biparamétrica es:
 
y  x  2 1  x 2  x 3  x 4 
1 1 1 1 5 17 6
x  x  ... 
 2 6 8 120 720 
 1 1 5 1 6 
x  x2  x4  x  x  ... 
 4 20 40 

EL MÉTODO DE FROBENIUS
Resumen de la teoría
Definición: Sea la ecuación lineal x 2 y  x  xP  x   y   x 2Q  x   y  0, donde
n n
p  x   xP  x  , q  x   x 2Q  x  , además P  x    a k x k , Q  x    b k x k
k 0 k 0

son analíticos en x  0. A la ecuación r  r  a 0r  b 0  0 se le llama


2

 
ecuación indicial y a las raíces r1 , r2 se les llaman raíces indiciales.
De acuerdo a las soluciones se presentan los siguientes casos:
1) Si r1  r2 y además r1  r2  , entonces existen dos soluciones linealmente
 
independientes: y 1  x    c n x , y 2  x    b nx
n  r1 n  r2

n0 n0

2) Si r1  r2 y además r1  r2  , entonces existen dos soluciones linealmente


 
independientes: y 1  x    c n x , con c 0  0, y 2  x   c y 1  x  ln x   b n x
n  r1 n  r2
, con b 0  0
n0 n0

3) Si r1  r2 , entonces existen dos soluciones linealmente independientes:


 
y 1  x    c nx , con c 0  0, y 2  x   y 1  x  ln x   b n x
n  r1 n  r2

n0 n0

La función y sus derivadas:


  
y   c n x n  r , y   c n  n  r  x n  r  1 , y   c n  n  r  n  r  1 x n  r  2
n0 n0 n0

 2
Ejemplo 3: Resolver x 2 y   x   y  0
 9

 2 
 2
x 2 y   x   y  0  x 2  c n  n  r  n  r  1 x n  r  2   x   y  0
 9 n0  9
 
  c n  n  r  n  r  1 x n  r  2 2   c n x n  r  1 
n0 n0

2
 
9 n0
c nx n  r  0
 
   n  r  n  r  1 c n x n  r   c n x n  r  1 
n0 n0

2
 
9 n0
c n x nr  0

Factorizar x r , luego homologar las sumatorias y x k :


 
2 
   n  r  n  r  1 c n x n  r   c n x n  1 r   c n x nr  0
n0 n0 9 n0
 
2  
    n  r  n  r  1 c n x n x r   c n x n  1x r   c n x n x r   0
n  0 n0 9 n0 
 
 
   
    n  r  n  r  1 c n x n   c n x n  1   c n x n  x r  0
2
n  0 n0 9 n0 
 kn k  n1 kn 
 n  k 1 
 
2  
    k  r  k  r  1 c k x k   c k  1 x k  1  1   c k x k  x r  0
k 0 k 1 9 k 0 
  
2   
  r  r  1 c 0    k  r  k  r  1 c k   c k  1  c 0   c k  x k  x r  0
2
 k 1 k 1 9 9 k  1  
  2    2  
   r  r  1   c 0      k  r  k  r  1   c k  c k  1  x k   x r  0
  9  k 1  9 
  
Resolver:
 2  2 2
 r  r  1  9  c 0  0   r  r  9  c 0  0
   


 9r 2  9r  2 c0  
0
9
0

  3r  2  3r  1 c 0  0, tomando el valor arbitrario c 0  1


2 1 2 1 2 1 1
 r1 
 r2  , como  y    , entonces la
3 3 3 3 3 3 3
ecuación tiene dos soluciones linealmente independientes:
2 c k 1 c k 1 3 c k 1
a) Para r1  , c  1 y ck    , k  1, 2, 3, 4,...
3 0  2  1 2 k  3 k  1
 k  r  k  r  1  92  k  3  k  3   9
  
3c0 3
k  1  c1  

1 3  1  1  4
3
3 
3c1
  
4 9
k  2  c2 

2 3 2  1 14 
56
 9 
3 
3c2
  
56 9
k  3  c3 

3 3  3  1 30 560
 9 
3
3c3 560 
 
27
k  4  c4  

4 3 4  1 52  29120

Luego y 1  x   1 
3 9 2 9 3 27
x x  x  x 4  ...
4 56 560 29120
1 3 c k 1
b) Análogamente, se resuelve para r1  , c 0  1 y c k  , k  1, 2, 3, 4,...
3 k  3 k  1

y 2  x  1  x  ...  comprobarlo 
3 9 2 9 3 27 4
x x  x 
2 20 160 7040

Solución general: y  x   C 1x 1 y 1  x   C 2 x 2 y 2  x 
r r

2
 
y  x   C 1x 3  1  x 
3 9 2 9 3 27
x  x  x 4  ...  
 4 56 560 29120 
1
 3 9 2 9 3 27 4 
C 2 x 3  1  x  x  x  x  ... 
 2 20 160 7040 

Ejemplo del método de Frobenius cuando ciertos valores de las constantes no están definidos

3
Ejemplo 4: y  y  2 y  0
x
Paso 1: Escribir la ecuación diferencial en términos de sumatorias y simplificar si es posible.

3 3  
y  y  2 y    n  r  n  r  1 c n x n  r  2    n  r  c n x n  r  1  2 c n x n  r
x n0 x n0 n0
  
   n  r  n  r  1 c n x n  r  2  3  n  r  c n  x n  r  1   2 c n x n  r
1
n0 n0 x  n0
  
   n  r  n  r  1 c n x n  r  2  3  n  r  c n x n  r  2  2 c n x n  r
n0 n0 n0
r
Paso 2: Dejar cada término de la ecuación con el factor x y las sumatorias iniciando con el
mismo valor de k . Además factorizar la ecuación para x r .
  

  n  r  n  r  1 c nx n  r  2  3  n  r  c nx n  r  2  2 c nx n  r  0
n0 n0 n0
 

  n  r  n  r  1  3  n  r c
n0
n x n  r  2  2  c nx n  r  0
n0
 

  n  r  n  r  1  3  c
n0
n
x n  r  2  2 c n x n  r  0
n0

 
 
   r
   n  r  n  r  2   c n x  2  c n x  x  0
n2 n

n  0 n0 
 k n2 k 0 
 n k  2 
   
      k  2
 
k22
k  2  r k  2  r  2 c x  2 c kx k  x r  0
k 2 k 0 
    r
   k  r  2  k  r  4 c k  2 x  2  c k x  x  0, luego
k k

k   2 k0 
 
k   2     2  r  2   2  r  4   c  2  2 x  2    k  r  2  k  r  4 c k  2 x k
 k  1

 
2  c k x k  x r  0, luego
k0 
 
k   1  r  r  2  c 0 x  2    1  r  2   1  r  4   c  1  2 x  1    k  r  2  k  r  4 c k  2 x k
 k0

 
2  c k x k  x r  0
k0 
   
 r  r  2  c 0 x 2
 
 r  1  r  3  c
 1 x 1
  
 k  r  2  k  r  4  c
 k2 x k
 2  c kx k  x r  0
 k0 k0 
   r
r  r  2  c 0 x   r  1 r  3 c 1x     k  r  2  k  r  4   c k  2  2c k  x  x  0
2 1 k

 k0 
r  r  2  c 0  0,  r  1 r  3 c 1  0    k  r  2  k  r  4   c k  2  2c k  0

Paso 3: Resolver para r de r  r  2   0 y despejar para c k  2 .


a) Si r  r  2   0, c 0  1 entonces r1  0  r 2   2
b) Si  r  1 r  3   0, con c 1  0
c) Si   k  r  2  k  r  4   c k  2  2c k  0 entonces
2
c k2  c k , k  0,1,2,3,
 k  r  2  k  r  4 
Paso 4: Calcular c k para los valores de k  0,1,2,3, .
1 2
Para r1  , c 0  1, c 1  0 se tiene c k  2  c , k  0,1,2,3,
2  k  2  k  4  k
2 1 1
k  0 c 2  c0  c2  c0 
 0  2  0  4  4 4
2 2
k  1 c 3  c1  c3  c1  0
 1  2 1  4  15

c 2  c 4    
2 1 1 1 1
k  2 c 4  c2  c4 
 2  2  2  4  12 12  4  48
2 2
k  3 c 5  c 3 c 5  c 3  c 5  0
 3  2  3  4  35
1 1
Luego y 1  x   1  x 2  x 4 
4 48
2
Para r2   2, c 0  1, c 1  0 se tiene c k  2  c k , k  0,1,2,3
k  k  2
2 2
k  0 c2  c 0  c 2  c 0 = no definido
0 0  2 0
2 2
k  1 c 3  c1  c 3  c1  0
11  2  3
2 1
k  2 c4  c 2  c 4  c 2  c 4  no definido
2 2  2 4
2 2
k  3 c5  c 3 c 5  c 3  c 5  0
3 3  2 15
Puesto que para r2   2 la relación de recurrencia tiene valores no definidos entonces, se
aplica la formula y 2  x  :
 e   P  x  dx
Paso 5: Calcular la segunda solución con la formula y 2  x   y 1  x   dx
  y  x  2
 1 
Escribiendo la ecuación en forma estándar se tiene:
3 3
y  y  2 y  0  P  x  
x x
Luego sustituyendo en la fórmula:

 
3
 dx
 e x
y 2  x  y 1  x  2
dx
 1  1 x 2  1 x 4  
  4 48 

  3
1
 dx
 e x 
 y 1  x  dx
 1 2 1 4  1 2 1 4 
  1  x  x   1  x  x  
 4 48  4 48 
  3ln x

 y 1  x  e
dx
 1 2 5 4 1 6 1
 1  x  x  x  x 8

2 48 96 2304
  3
e ln x
 y 1  x  dx , x  
 1 2 5 4 1 6 1
 1  2 x  48 x  96 x  2304 x
8


x 3
 y 1  x  dx
 1 1 x2  5 x4  1 x6  1 x8

2 48 96 2304

 y 1  x 
1
dx
 3 1 2 5 4 1 6 1 8
 x 1  x  x  x  x
 2 48 96 2304 

 y 1  x 
1
dx , por división larga resulta:
 3 1 5 5 7 1 9 1
x  x  x  x  11
x
2 48 96 2304

1 5 5 7 1 9
1  0x 2  0x 4  0x 6 x3  x  x  x
2 48 96
1 2 5 4 1 6 1 1 7
1 x  x  x 3
  x
2 48 96 x 2 x 48
1 5 1
 x 2  x 4  x 6  0x 8
2 48 96
1 2 1 4 5 6 1 8
 x  x  x  x
2 4 96 192
7 4 1 6 1 8
x  x  x  0 x 10
48 24 192
7 4 7 6 35 8 7
x  x  x  x 10
48 96 2304 4608
1 23 8 7
 x6  x  x 10 
32 2304 4608
Luego:
 1 
y 2  x  y 1  x   3 
1 7
 x  dx
 x 2 x 48 
 1 
 y 1  x   3 
1 7
 x   dx
 x 2 x 48 
 1 
 y 1  x    2  ln x  x 2  
1 7
 2x 2 96 
 1 
 y 1  x  ln x  y 1  x    2  x 2 
1 7

2  2x 96 
Por tanto la solución general es y  C 3 y 3  x   C 4 y 4  x 

  1  1 
  2 C 2 y 1  x  ln x  C 2 y 1  x    2  96 x 
1 1 4 7 2
y  C1 1  x 2  x  
 4 48   2x 
Paso 6: Resolver para r de  r  1 r  3   0 y despejar para c k  2 .
a) Si  r  1 r  3   0, c 1  1 entonces r 3   1  r 2   3
b) Si r  r  2   0, con c 0  0
c) Si   k  r  2  k  r  4   c k  2  2c k  0 entonces
2
c k2  c k , k  0,1,2,3,
 k  r  2  k  r  4 
Paso 7: Calcular c k para los valores de k  0,1,2,3, .
2
Para r1   1, c 0  0, c 1  1 se tiene c k  2  c k , k  0,1,2,3,
 k  1 k  3 
2 2
k  0 c 2  c0  c2 
c0  0
 0  1 0  3  3
2 1 1
k  1 c 3  c1  c3  c1 
 1  11  3  4 4
2 2
k  2 c 4  c2  c4  c2  0
 2  1 2  3  15

c 3  c 5    
2 1 1 1 1
k  3 c 5  c 3 c 5 
 3  1 3  3  12 12  4  48
1 3 1 3
Luego y 3  x   x  x  x 
4 48
2
Para r2   3, c 0  0, c 1  1 se tiene c k  2  c k , k  0,1,2,3
 k  1 k  1
2
k  0 c2  c 0  c 2   2c 0  0
 k  1 k  1
2 2
k  1 c 3  c1  c3  c  no definido
 k  1 k  1 0 1
2 2
k  2 c4  c2  c4  c  c4  0
 k  1 k  1 3 2
2 1
k  3 c5  c3 c5  c  c 5  no ndefinido
 k  1 k  1 4 3
Puesto que para r2   3 la relación de recurrencia tiene valores indefinidos entonces se aplica
la formula y 4  x 
 e   P  x  dx
Paso 7: Calcular la segunda solución con la formula y 4  x   y 3  x   dx
  y  x  2
 3 
Luego sustituyendo en la fórmula:
   dx
3

 e x
y 4  x  y 3  x 2 dx
 1
 x  x 
3 1 
x  
5
 4 48 

  3  dx
1

 e x
y 4  x  y 3  x dx

x   
1 1 5 1 1 5
 x  x3   x x3  x  
 4 48  4 48 
 e  3 ln x
 y 3  x dx
 x 2  1 x 4  5 x 6  1 x 8  1 x 10 

2 48 96 2304
 3

 e ln x
 y 3  x dx , x  
 x 2  1 x 4  5 x 6  1 x 8  1 x 10 
 2 48 96 2304
 x 3
 y 3  x  dx
 x 2  1 x 4  5 x 6  1 x 8  1 x 10 

2 48 96 2304
 1
 y 3  x dx
 3 2 1 4 5 6 1 8 1 
 x x  x  x  x  x  
10

 2 48 96 2304 
 1
 y 3  x dx
 5 1 5 1 1
x  x  x 9  x 11  x 13 
7

2 48 96 2304

1 5 1
Resolviendo por división larga 1  x 3  x 5  x 7  x 9  se tiene
2 48 96
1 7 5 9 1 11
1  0x 2  0x 4  0x 6 x5  x  x  x 
2 48 96
1 2 5 4 1 6 1 1 7
1 x  x  x 5
 3 
2 48 96 x 2x 48 x
1 5 1
 x 2  x 4  x 6  0x 8
2 48 96
1 1 5 6 1 8
 x2  x4  x  x
2 4 96 192
7 4 1 6 1 8
x  x  x  0 x 10
48 24 192
7 4 7 6 35 8 7
x  x  x  x 10
48 96 2304 4608
1 23 8 7
 x6  x  x 10 
32 2304 4608
Luego
 1 1 7 
y 4  x  y 3  x  5  3    dx
 x 2x 48 x 

 y 3  x    4   ln x  
1 1 7
 4x 
2
4x 48
7  1 1 
 y 3  x  ln x  y 3  x    4   
48  4x 4x 2 

Por tanto la solución general es y  C 3 y 3  x   C 4 y 4  x 

  7  1 
  48 C 4 y 3  x  ln x  C 4 y 3  x    4  2 
1 1 3 1
y  C3 x  x3  x  
 4 48   4x 4x 

También podría gustarte