Ecuaciones diferenciales de orden superior - Canek - UAM

Ecuaciones diferenciales de orden superior - Canek - UAM Ecuaciones diferenciales de orden superior - Canek - UAM

03.07.2013 Views

2 Ecuaciones diferenciales ordinarias donde u1 & u2 son funciones de x desconocidas. Para determinarlas, buscaremos qué condiciones deben cumplir para que la función propuesta yp sea solución de la ED (4.1). Lo primero que necesitamos hacer es calcular las derivadas y 0 00 p , yp para utilizarlas en (4.1) y ver cuáles propiedades deben cumplir u1, u2: Antes de obtener y 00 p , suponemos que yp D u1 1 C u2 2 ) y 0 p D u 0 1 1 C u1 0 1 C u 0 2 2 C u2 0 2 : (4.3) u 0 1 1 C u 0 2 2 D 0: Esto se hace con la finalidad de que en la expresión de y 00 p no aparezcan u 00 1 00 & u2 , ya que la inclusión de estas segundas derivadas en y 00 p haría mucho más compleja la obtención de las funciones u1 & u2. Luego entonces, usando la condición anterior en (4.3): por lo cual, al volver a derivar, obtenemos: y 0 p D u1 0 1 C u2 0 2 ; y 00 p D u 0 1 0 00 1 C u1 1 C u 0 2 0 00 2 C u2 2 : Ahora bien, yp es solución de la ED lineal y 00 C p.x/y 0 C q.x/y D g.x/ sólo si se cumple que y 00 p C p.x/y 0 p C q.x/yp D g.x/: Esta condición, al usar las expresiones obtenidas de yp, y 0 p Al factorizar u1, u2: Œu 0 1 0 00 1 C u1 y 00 p u 0 1 0 1 C u1Œ 1 C u 0 2 0 2 C u2 00 2 C p.x/Œu1 0 1 C u2 0 š2 0 yp 00 1 C p.x/ 0 1 C q.x/ 1 C u 0 2 0 2 C u2Œ Pero por ser 1 & 2 soluciones de la ED homogénea se cumple: 00 1 C p.x/ Entonces debe cumplirse que 0 1 C q.x/ 1 D 0 & u 0 1 0 1 C u 0 2 0 2 D g.x/: 00 & yp nos da: Hemos obtenido un sistema de dos ecuaciones con dos incógnitas: u 0 1 1 C u 0 2 2 D 0I u 0 1 0 1 C u 0 2 0 2 D g.x/ : El determinante s del sistema es s D 1 2 0 0 1 2 š yp C q.x/Œu1 1 C u2 2 D g.x/: 00 2 C p.x/ 0 2 C q.x/ 2 D g.x/: 00 2 C p.x/ 0 2 C q.x/ 2 D 0: D W. 1; 2/: Y debido a que el wronskiano W. 1; 2/.x/ ¤ 0, entonces s ¤ 0, por lo que el sistema de ecuaciones tiene una única solución. Dicha solución única la calculamos usando la regla de Cramer: u 0 1 D 0 2 0 g.x/ 2 D s g.x/ 2 W. 1; 2/ ) u 0 g.x/ 2.x/ 1 .x/ D W. 1.x/; 2.x// : u 0 2 D 1 0 0 1 g.x/ D s g.x/ 1 W. 1; 2/ ) u 0 g.x/ 1.x/ 2 .x/ D W. 1.x/; 2.x// : (4.4)

4.7 Variación de parámetros 3 De donde obtenemos u1 & u2 mediante integración u1 D g.x/ 2.x/ W Œ 1.x/; 2.x/ dx & u2 D g.x/ 1.x/ W Œ 1.x/; 2.x/ dx: Sustituyendo u1.x/ & u2.x/ en yp.x/ se tiene que la solución particular de la lineal es yp.x/ D 1.x/ g.x/ 2.x/ W Œ 1.x/; 2.x/ dx C 2.x/ Finalmente, podemos escribir la solución general de la ED lineal como con la yp.x/ obtenida previamente. y.x/ D yp.x/ C .x/I y.x/ D yp.x/ C Œc1 1.x/ C c2 2.x/; g.x/ 1.x/ dx: (4.5) W Œ 1.x/; 2.x/ Ejemplo 4.7.1 Utilizando el método de variación de parámetros, calcular una solución particular y escribir la solución general de la ED lineal x 2 y 00 4xy 0 C 6y D 1 x ; considerando que y1 D x 2 & y2 D x 3 forman un conjunto fundamental de soluciones para la ED homogénea asociada x 2 y 00 4xy 0 C 6y D 0: H Sea yp.x/ D u1.x/y1.x/ C u2.x/y2.x/ una solución particular. Entonces: yp D u1x 2 C u2x 3 ) y 0 p D u 0 1 x2 C 2u1x C u 0 2 x3 C 3u2x 2 : Suponiendo que se tiene que Sustituyendo en la ED normalizada se obtiene: Entonces u 0 1 & u 0 2 u 0 1 x2 C u 0 2 x3 D 0; (4.6) y 0 p D 2u1x C 3u2x 2 & y 00 p D 2u 0 1 x C 2u1 C 3u 0 2 x2 C 6u2x: y 00 4 x y 0 C 6 1 y D ; x2 x3 .2u 0 1x C 2u1 C 3u 0 2x2 C 6u2x/ 4 x .2u1x C 3u2x 2 / C 6 x2 .u1x 2 C u2x 3 / D 1 x3 ) ) 2xu 0 1 C u1.2 8 C 6/ C 3x 2 u 0 2 C u2.6x 12x C 6x/ D 1 x3 ) ) 2xu 0 1 C 3x2u 0 1 2 D : x3 (4.7) El determinante del sistema es satisfacen el sistema formado por las ecuaciones (4.6) y (4.7) x 2 u 0 1 C x3 u 0 2 D 0I 2xu 0 1 C 3x2 u 0 2 D x 3 : W D x2 x3 2x 3x2 D 3x4 2x 4 D x 4 :

2 <strong>Ecuaciones</strong> <strong>diferenciales</strong> ordinarias<br />

don<strong>de</strong> u1 & u2 son funciones <strong>de</strong> x <strong>de</strong>sconocidas.<br />

Para <strong>de</strong>terminarlas, buscaremos qué condiciones <strong>de</strong>ben cumplir para que la función propuesta yp sea solución<br />

<strong>de</strong> la ED (4.1). Lo primero que necesitamos hacer es calcular las <strong>de</strong>rivadas y 0 00<br />

p , yp para utilizarlas en<br />

(4.1) y ver cuáles propieda<strong>de</strong>s <strong>de</strong>ben cumplir u1, u2:<br />

Antes <strong>de</strong> obtener y 00<br />

p , suponemos que<br />

yp D u1 1 C u2 2 ) y 0 p D u 0 1 1 C u1 0 1 C u 0 2 2 C u2 0 2 : (4.3)<br />

u 0 1 1 C u 0 2 2 D 0:<br />

Esto se hace con la finalidad <strong>de</strong> que en la expresión <strong>de</strong> y 00<br />

p<br />

no aparezcan u 00<br />

1<br />

00 & u2 , ya que la inclusión <strong>de</strong><br />

estas segundas <strong>de</strong>rivadas en y 00<br />

p haría mucho más compleja la obtención <strong>de</strong> las funciones u1 & u2.<br />

Luego entonces, usando la condición anterior en (4.3):<br />

por lo cual, al volver a <strong>de</strong>rivar, obtenemos:<br />

y 0 p D u1 0 1 C u2 0 2 ;<br />

y 00<br />

p D u 0 1 0 00<br />

1 C u1 1 C u 0 2 0 00<br />

2 C u2 2 :<br />

Ahora bien, yp es solución <strong>de</strong> la ED lineal y 00 C p.x/y 0 C q.x/y D g.x/ sólo si se cumple que<br />

y 00<br />

p C p.x/y 0 p C q.x/yp D g.x/:<br />

Esta condición, al usar las expresiones obtenidas <strong>de</strong> yp, y 0 p<br />

Al factorizar u1, u2:<br />

Œu 0 1 0 00<br />

1 C u1<br />

y 00<br />

p<br />

u 0 1 0 1<br />

C u1Œ<br />

1 C u 0 2 0 2<br />

C u2 00<br />

2<br />

C p.x/Œu1 0 1 C u2 0 š2 0<br />

yp 00<br />

1 C p.x/ 0 1 C q.x/ 1 C u 0 2 0 2 C u2Œ<br />

Pero por ser 1 & 2 soluciones <strong>de</strong> la ED homogénea se cumple:<br />

00<br />

1 C p.x/<br />

Entonces <strong>de</strong>be cumplirse que<br />

0<br />

1 C q.x/ 1 D 0 &<br />

u 0 1 0 1 C u 0 2 0 2 D g.x/:<br />

00 & yp nos da:<br />

Hemos obtenido un sistema <strong>de</strong> dos ecuaciones con dos incógnitas:<br />

<br />

u 0 1 1 C u 0 2 2 D 0I<br />

u 0 1 0 1 C u 0 2 0 2 D g.x/ :<br />

El <strong>de</strong>terminante s <strong>de</strong>l sistema es<br />

<br />

<br />

s D <br />

<br />

1 2<br />

0 0<br />

1 2<br />

š yp<br />

C q.x/Œu1 1 C u2 2<br />

D g.x/:<br />

00<br />

2 C p.x/<br />

0<br />

2 C q.x/ 2 D g.x/:<br />

00<br />

2 C p.x/ 0 2 C q.x/ 2 D 0:<br />

<br />

<br />

<br />

D W. 1; 2/:<br />

Y <strong>de</strong>bido a que el wronskiano W. 1; 2/.x/ ¤ 0, entonces s ¤ 0, por lo que el sistema <strong>de</strong> ecuaciones tiene<br />

una única solución. Dicha solución única la calculamos usando la regla <strong>de</strong> Cramer:<br />

u 0 1 D<br />

<br />

<br />

<br />

0 <br />

2<br />

0 g.x/ <br />

2<br />

D<br />

s<br />

g.x/ 2<br />

W. 1; 2/ ) u 0 g.x/ 2.x/<br />

1 .x/ D<br />

W. 1.x/; 2.x// :<br />

u 0 2 D<br />

<br />

<br />

<br />

1 0<br />

0<br />

1 g.x/<br />

<br />

<br />

<br />

<br />

D<br />

s<br />

g.x/ 1<br />

W. 1; 2/ ) u 0 g.x/ 1.x/<br />

2 .x/ D<br />

W. 1.x/; 2.x// :<br />

(4.4)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!