Ecuaciones Diferenciales II

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

Matemática III

Instituto de Ingeniería y Agronomía




3. Ecuaciones diferenciales ordinarias lineales de segundo orden


En forma general una EDO lineal de segundo orden tiene la forma:

y 00 (x) + a1 (x)y 0 (x) + a2 (x)y(x) = b(x) (9)

Y su ecuación homogenea asociada:

y 00 (x) + a1 (x)y 0 (x) + a2 (x)y(x) = 0 (10)

Nuestro objetivo es encontrar soluciones a estas ecuaciones, previo a lo cual analizaremos


tres propiedades algebraicas que nos permitirán desarrollar métodos de resolución.

Propiedad 0: y(x) = 0 es solución de toda ecuación ordinaria lineal de segundo orden


homogénea, es decir:

y(x) = 0 es solución de: y 00 (x) + a1 (x)y 0 (x) + a2 (x)y(x) = 0


Es fácil vericarlo por simple inspección. Si y(x) = 0 entonces y 0 (x) = y 00 (x) = 0 y la ecuación
resulta ser idénticamente nula. A y(x) = 0 se la conoce como solución trivial.

Propiedad 1: y1 (x) e y2 (x) dos soluciones de y 00 (x) + a1 (x)y 0 (x) + a2 (x)y(x) = 0.


Sean
00
Si y(x) = C1 y1 (x) + C2 y2 (x) entonces y(x) también es solución de y (x) + a1 (x)y (x) +
0

a2 (x)y(x) = 0.
Es decir, si encuentro dos soluciones particulares de mi EDO homogenea, cualquier combina-
ción lineal de ellas también será solución. Veamoslo: Si y1 (x) e y2 (x) son soluciones luego:

C1 [y100 (x) + a1 (x)y10 (x) + a2 (x)y1 (x)] = 0




C2 [y200 (x) + a1 (x)y20 (x) + a2 (x)y2 (x)] = 0


Sumando miembro a miembro:

(C1 y100 (x) + C2 y200 (x)) +a1 (x) (C1 y10 (x) + C2 y20 (x)) +a2 (x) (C1 y1 (x) + C2 y2 (x)) = 0
| {z } | {z } | {z }
y 00 (x) y 0 (x) y(x)
00 0
⇒ y (x) + a1 y (x) + a2 y(x) = 0
Esta propiedad se conoce como superposición y juega un rol muy importante en el estu-
dio de diversos sistemas y fenómenos físicos.

Las propiedades 0 y 1, nos serán utiles para trabajar la siguiente idea5 .


5 No es una idea sino un resultado formal de la teoría de las ecuaciones diferenciales, pero no será demos-
trada aquí.

14
Solución general de una EDO lineal homogénea Consideremos el conjunto de solucio-
00
nes de y (x) + a1 (x)y 0 (x) + a2 (x)y(x) = 0, al que llamaremos S

El cero del conjunto será la función y(x) = 0 que pertenece al conjunto de acuerdo a la
propiedad 1. Por otro lado, para cualquier par de funcionesy1 e y2 que pertenezcan al conjunto,
la combinación lineal de ellas también formará parte del conjunto, de acuerdo a la propiedad
2. Esto es: 
0∈S
C1 y1 + C2 y2 ∈ S
Para cualesquiera C1 y C2 en R (o C).
Un conjunto que cumple con estas dos propiedades constituye un Espacio Vectorial . La
6

dimensión de nuestro espacio vectorial S (el conjunto de todas las soluciones de nuestra ecua-
ción homogénea) es igual al orden de la ecuación, es decir:

El conjunto de soluciones de una EDO lineal homogénea forma un Espacio vectorial de


dimensión 2.

Por lo tanto, el conjunto S admite una base. Esto es un subconjunto de dos elementos
(su dimensión) que permiten generar cualquier solución que pertenezca a S y que además son
linealmente independientes entre sí. Es decir, si:

β = {u1 (x), u2 (x)} es base del conjunto S entonces


Si
para todo y(x) ∈ S podemos escribir y(x) = C1 u1 (x) + C2 u2 (x)
y ademásu1 (x) y u2 (x) son linealmente independientes.

Esto facilita mucho nuestro trabajo para hallar soluciones.

Siempre que podamos encontrar dos soluciones linealmente independientes podremos com-
binarlas linealmente para obtener la solución general de nuestra ecuación de segundo orden
lineal y homogenea.

Haciendo uso de esta idea, comenzaremos a trabajar en la resolución de nuestras ecuaciones.

Ya hemos visto al estudiar soluciones de la Ecuaciones diferenciales ordinarias de primer


orden que algunos métodos de resolución se basan en proponer soluciones, es decir, se propone
una función de determinada forma y luego se ajustabala propuesta para que la misma sea
una solución.
Ese es el camino que tomaremos.

3.1. Soluciones de EDO de 2do.orden homogeneas a coecientes constantes


Recordemos nuestra denición:
y 00 (x) + a1 (x)y 0 (x) + a2 (x)y(x) = 0 (11)

6 Matemática 2 - Módulo 2; Pp:48-54. IIyA; UNAJ; 2018

15
Trabajaremos con ecuaciones a coecientes constantes, es decir, a1 (x) = a1 = cte ∧ a2 (x) =
a2 = cte. A partir de ésta deniciíon la ecuación (9) resulta ahora:

y 00 (x) + a1 y 0 (x) + a2 y(x) = 0 (12)

Proponemos para la ecuación ( ??) una solución de la forma


y(x) = ceλx (13)

donde tanto c λ son constantes indeterminadas. Resolver la ecuación


como será encontrar los
valores posibles de λ 7 que hagan de nuestra ecuación una identidad.
Para ello, hallamos las derivadas de y(x):

y 0 (x) = cλeλx y 00 (x) = cλ2 eλx (14)

Reemplazamos en ( ??) obtenemos:


y 00 (x) + a1 y 0 (x) + a2 y(x) = cλ2 eλx + a1 cλeλx + a2 ceλx = 0 (15)

Sacando factor común ceλx .


cλ2 eλx 
+ a1 cλeλx + a2 ce λx
 =0
ceλx λ2 + a1 λ + a2 =0

Dado que c 6= 0 y que la función exponencial no se anula nunca, tenemos que:

λ2 + a1 λ + a2 = 0 (16)

Esta ecuación se conoce con el nombre de ecuación característica; y su solución nos dará los
posibles valores de la variable λ a determinar.
Analicemos la ecuación (16). Podemos obtener sus raíces (λ) según:
p
−a1 ± a21 − 4(1)(a2 )
λ=
2(1)
vemos que tenemos tres posibilidades. La ecuación característica puede tener dos raíces reales,
una única raíz o dos raíces complejas, de acuerdo a si el radicando toma valores positivos,
negativos o cero. Veamos cada uno de los casos:

1. a21 − 4a2 > 0 ⇒


p p
−a1 + | a21 − 4a2 | −a1 − | a21 − 4a2 |
λ1 = λ2 =
2 2
Si λ1 6= λ2 ∈ R, entonces nuestra propuesta involucra a dos soluciones

y1 (x) = eλ1 x y2 (x) = eλ2 x


que son entre sí, linealmente independientes (ver 4). De manera que cualquier combina-
ción lineal de estas soluciones será una solución a nuestro problema.
Luego podemos escribir la solución general como:

y(x) = c1 y1 (x) + c2 y2 (x) = c1 eλ1 x + c2 eλ2 x (17)

Trabajemos sobre el siguiente

7 Observe que los valores de c solo podran ser hallados si contamos con las condiciones iniciales del problema
(PVI)

16
Ejemplo: Hallar las solución general de la ecuación y 00 (x) + y 0 (x) − 6y = 0

Proponemos como solución particular y(x) = eλx , derivando y reemplazando en nuestra


ecuación diferecial, obtenemos:
 
eλx λ2 + λ − 6 = 0 ⇒ (λ − 2)(λ + 3) = 0 ⇒ λ1 = 2 ∧ λ2 = −3

reemplazando en nuestra propuesta de solución obtenemos:

y1 (x) = e2x ∧ y2 (x) = e−3x ⇒ y(x) = c1 e2x + c2 e−3x


es la solución general a nuestra ecuación diferencial. Puede vericar este resultado por
inspección.

2. a21 − 4a2 > 0 ⇒ λ1 = λ2 = − a21


Si λ1 = λ2 ∈ R, entonces nuestra propuesta resulta sólo en una solución particular:

y1 (x) = eλx
Para construir una solución general necesitaremos otra solución particular linealmente
independiente a y1 (x).

Si bien puede hallarse esta otra solución a partir de la que hemos obtenido (en general,
puede hallarse una segunda solución a partir de una conocida utilizando un método
conocido como reducción de orden) optaremos por proponer una segunda función, veri-
caremos que es solución y que es Linealmente Independiente con y1 (x) para, nalmente,
construir nuestra solución general.

Propondremos como solución a y2 (x) = x eλx . Para ver si efectivamente es una solución
particular vericamos reemplazando en ( ??) y obtenemos:
(2λ + xλ2 )eλx + a1 (1 + λx)eλx + a2 xeλx = 0
igualdad que nos lleva al siguiente sistema de ecuaciones:

2λ + a1 = 0
λ2 + λa1 + a2 = 0

La segunda de estas ecuaciones es, precisamente, nuestra ecuación característica (16)


que ya hemos resuelto y que tiene como solución a λ = −a1 /2. Reemplazando este valor
en la primer ecuación obtenemos una identidad lo que implica que nuestra y2 (x) es una
solución particular de ( ??).

Con ese procedimiento vericamos que la y2 (x) propuesta es solución particular de nues-
tra Ecuación diferencial. Nos resta ver que y1 (x) e y2 (x) son linealmente independientes.

Para ello consideramos el cociente (ver Anexo 4 ) y2 (x)/y1 (x) = xeλx /eλx = x 6= cons-
tante, lo que garantiza la independencia lineal de las funciones, luego construimos nuestra
solución general como combinación lineal de nuestras soluciones particulares:

y1 (x) = eλx ∧ y2 (x) = x eλx ⇒ y(x) = c1 eλx + c2 xeλx (18)

17
Ejemplo: Hallar la solución general de y 00 + 2y 0 + y = 0.

Proponemos como solución particular y(x) = eλx , derivando y reemplazando en nuestra


ecuación diferencial, obtenemos:

eλx (λ2 + λ + 1) = 0 ⇒ (λ + 1)2 = 0 de manera que λ = −1 luego:

y1 (x) = eλx ∧ y2 (x) = x eλx ⇒ y(x) = c1 eλx + c2 xeλx


es solución general de y 00 + 2y 0 + y = 0.
3. a21 − 4a2 < 0 ⇒
p p
−a1 + i| a21 − 4a2 | −a1 − i| a21 − 4a2 |
λ1 = λ2 =
2 2
p
Escribimos, por simplicidad: λ1 = −a1 /2 + i(
2
|a1 − 4a1 |) = α + iβ ∧ λ2 = −a1 /2 −
p
i( |a21 − 4a1 |) = α − iβ

La solución general, siendo que las raíces son distintas tendrá la forma:

y(x) = c1 e(α+iβ)x + c2 e(α−iβ)x


Utilizando propiedades de la funcio« exponencial y la identidad de Euler podemos rees-
cribir:
y(x) = c1 eαx ( cos(βx) + i sin(βx)) + c2 eαx ( cos(βx) − i sin(βx))
reagrupando y llamando A = c1 + c2 y B = c1 − c2 :
y(x) = A cos(βx)eαx + i B sin(βx)eαx

Un importante resultado del algebra lineal indica que dada una solución compleja de
una ecuación diferencial real, sus partes real e imaginaria son también soluciones (parti-
culares) de la ecuación (ver Propiedad 4 del Anexo 4), es decir que podemos considerar:

y1 (x) = Re(y(x)) = A cos(βx)eαx y2 (x) = Im(y(x)) = B sin(βx)eαx


Siendo además que y1 (x) e y2 (x) son linealmente independientes podemos construir
nuestra solución general real como:

y1 (x) = A cos(βx)eαx ∧ y2 (x) = B sin(βx)eαx ⇒ y(x) = A cos(βx)eαx + B sin(βx)eαx


(19)

Ejemplo: soluciones complejas. Consideremos la ED y 00 − 6y 0 + 13y = 0

Proponemos como solución particular y(x) = eλx , derivando y reemplazando en nuestra


ecuación diferencial, obtenemos:
(λ2 − 6λ + 13) = 0 ⇒ λ1 = 3 + 2i ∧ λ2 = 3 − 2i donde α=3 y β = 2, luego:

y1 (x) = e3x cos(2x) ∧ y2 (x) = e3x sin(2x) ⇒ y(x) = c1 e3x cos(2x) + c2 e3x sin(2x)
y(x) = e3x (c1 cos(2x) + c2 sin(2x))
es la solución general a nuestra ecuación diferencial.

18

También podría gustarte