Confiabilidad Operacional

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

UNIVERSIDAD TECNOLÓGICA METROPOLITANA

FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.


DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

UNIDAD 5: ECUACIONES DIFERENCIALES LINEALES DE ORDEN N.

Forma general

d n 1 y d n2 y
       a n 1 ( x)  a n ( x) y
dny dy
n 1 n2
a 0 ( x) n
a1 ( x ) a 2 ( x ) F ( x)
dx dx dx dx

si ai son constantes i , la ecuación se llama ecuación diferencial lineal con coeficientes


constante. En caso contrario ecuación diferencial lineal de coeficientes variables.

Ejemplos
3y’’’+ 2y’’–4y’+8y = 3e-x-5x2 ec. diferencial lineal con coeficientes constantes, de orden 3
xy’’ – 2y’ + x2 y = x3 – 2 ec. diferencial lineal con coeficiente variables, de orden 2.

2 3dy 2 d2y
Usando los operadores : D, D , D , ….., donde Dy = , Dy = , …...
dx dx 2
Anotaremos

ao Dny + a1 Dn-1y + ---+ an-1 Dy + any = F o bien (ao Dn + a1 Dn-1 + --- + an-1 D + an) y = F

Si I (D) = ao Dn + a1 Dn-1 + --- + an , anotaremos I (D) y = F

Propiedades de los operadores


1) Dn (u+v) = Dn (u) + Dn (v) , u,v funciones diferenciables
2) Dn (D u) = D Dn (u) , D = constante
3) I (D) (u + v) = I (D)(u) + I (D)(v)
4) I (D) (D u) = D I (D) u, D = constante
? Dn y I (D) son operadores lineales

Sea I (D) y = F(x) , F(x) z 0 ecuación diferencial general de orden n.


Definición :

Se puede asociar con la ecuación I (D) y = 0 llamada ecuación complementaria o ecuación


homogénea de la ecuación diferencial lineal general.

Si y = u (x) es cualquier solución de la ecuación I (D) y = F(x) e y = v(x) es cualquier


Teorema fundamental I

solución de la ecuación I (D) y = 0, entonces y = u(x) + v(x) es solución de la ecuación I (D)


y = F(x).
Si v(x) tiene n constantes arbitrarias, donde n es el orden de la ecuación, entonces es la solución

Si u(x) no tiene constantes arbitrarias, es una solución particular de I (D) y = F(x).


general de la ecuación complementaria.

73
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Definición : Solución complementaria de la ecuación I (D) y = F(x)


Corresponde a la solución general de la ecuación complementaria y la denotaremos por yc.

Definición : Solución particular de la ecuación I (D) y = F(x)


Es una solución seleccionada de la ec. dada (sin constantes arbitrarias) y la denotaremos por yp

La solución general de la ecuación diferencial I (D) y = F(x) viene dada por


Teorema fundamental II
y = yc + yp

Si ai (x) y F(x) son funciones continuas en >a, b@ y ao (x) z 0 y po, …., pn, son constantes
Teorema de existencia – unicidad

Entonces ! solución y(x), tal que >ao(x) Dn + a1(x) Dn-1 + …. + an(x)@ y = F(x), cumple las
dadas.

condiciones y(c) = po , y’(c) = p1, ….., y(n-1) (c) = pn, c  >a, b@

Nota: El teorema proporciona sólo condiciones suficientes, es decir, aún si las condiciones
enunciadas no se satisfacen todas, pueden existir soluciones únicas.

Solución general de la ecuación complementaria de una ecuación diferencial de coeficientes


constantes
Considerar una ecuación auxiliar, que se obtiene al reemplazar el operador D por m en la
ecuación diferencial, se obtiene un polinomio en m y se calculan sus raíces. Se pueden presentar
tres casos.

Caso 1: Raíces reales distintas


Si m1, m2,...., mk son raíces reales distintas de la ecuación auxiliar en m, entonces
y1 e m1x , y 2 e m2 x ......... y n e mn x son soluciones linealmente independientes de la ecuación
diferencial complementaria I (D) y = 0, y también y = c1 e m1x  c 2 e m2 x  !  c k e mk x es solución
de la ecuación complementaria.

Si la ecuación diferencial es de orden n y se obtiene n raíces reales distintas m1, m2, ..., mn,
entonces y = c1 e m1x + ......... + cn e mn x es la solución general de la ecuación complementaria, es
decir, es la solución complementaria.

Ejemplos: Resolver las siguientes ecuaciones diferenciales:

1) y’’ – 3y’ + 2y = 0
Solución:
m2 – 3m + 2 = 0
(m – 1) (m – 2) = 0 Ÿ m = 1 › m = 2 Ÿ y = ex , y = e2x sea solución
y = c1 ex + c2 e2x es la solución general.

74
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

2) y’’’ – 6y’’ + 11 y’ – 6y = 0
Solución:
m3 – 6m2 + 11 m – 6 = 0
(m – 1) (m – 2) (m – 3) = 0 Ÿ m = 1, 2, 3.
y = c1 ex + c2 e2x + c3 e3x es la solución general.

Caso 2: Raíces reales iguales


Si la ecuación auxiliar tiene mk como raíz con multiplicidad k, entonces
y1 e m1x , y 2 xe m2 x ......... y k x k 1e mk x son soluciones linealmente independientes de la
ecuación diferencial complementaria I (D) y = 0, y también y = (c1 + c2 x + .... + ck xk-1) e mk x es
solución de la ecuación diferencial.

Ejemplos: Resolver las siguientes ecuaciones diferenciales:

1) (D2 – 6D + 9)y = 0
Solución:
m2 – 6 m + 9 = 0, luego (m – 3)2 = 0, es decir m = 3 con multiplicidad 2.
? y = (c1 + c2 x ) e3x solución general.

2) (D3 – 6D2 + 12 D - 8)y = 0


Solución:
m3 – 6m2 + 12 m - 8 = 0, luego (m – 2)3 = 0, es decir m = 2, con multiplicidad 3.
y = (c1 + c2 x + c3 x2) e2x es la solución general.

3) (D6 – 6D5 + 12 D4 – 6D3 – 9D2 + 12D – 4) y = 0


Solución:
m6 – 6m5 + 12m4 – 6m3 – 9m2 + 12m – 4 = 0
(m – 1)3 (m – 2)2 (m + 1) = 0
? m = 1 con multiplicidad 3
m = 2 con multiplicidad 2
m = - 1 con multiplicidad 1

? y = (c1 + c2 x + c3 x2 ) ex + ( c4 + c5 x ) e2x + (c6 e-x) es la solución general.

75
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Caso 3: Raíces complejas


Si la ecuación auxiliar presenta raíces complejas conjugadas, de la forma: m1 = a +ib, m2 = a – ib,
entonces una solución de la ecuación diferencial es

y = eax >c1 sen (b x) + c2 cos (b x)@

pues e(a + ib)x = eax ˜ eibx = eax (cos b x + i sen bx)

1    ....
x x2 x3
ex
1! 2! 3!

Nota: e ibx 1     .... 1  ibx    ....


ibx (ibx) 2 (ibx) 3 b 2 x 2 ib 3 x 3
1! 2! 3! 2! 3!
ª b2 x2 º ª º
e ibx «1   ....»  i «bx   ....» cos(bx)  isen(bx)
b3 x3
¬ 2! ¼ ¬ 3! ¼
ax ax
y1 = e cos(bx) y2 = e sen(bx) son soluciones independientes

Ejemplos

1) (D2 + 2D + 5) y = 0
m2 + 2m + 5 = 0
m=-1r 2i
y = e-x >c1 sen 2x + c2 cos 2x@

2) (D4 – 5D2 + 12 D + 28) y = 0


m4 – 5m2 + 12 m + 28 = 0
raíces : -2 , -2, 2 r 3 i
y = (c1 + c2x) e-2x + e2x (c3 sen 3 x + c4 cos 3 x)

3) (D2 – 2D + 5)2 y = 0
(m2 – 2 m + 5)2 = 0
m2 – 2m + 5 = 0 Ÿ m = 1 r 2i con multiplicidad 2.
? y = ex (c1 cos 2x + c2 sen 2x) + x ex (c3 cos 2x + c4 sen 2x).

76
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Ejercicios propuestos:

1. Encuentre las soluciones generales de cada una de las siguientes ecuaciones:


a) y’’ + 4y’ – 5y = 0 b) (4D2 – 25)y = 0 c) y’’ = 4y
d) 2y’’’ – 5y’’ + 2y’ = 0 e) I’’(t) – 4I’(t) + 2I(t) = 0 f) (D3 +2D2 – 5D – 6)y = 0

2. Encuentre las soluciones que satisfagan las condiciones dadas.


a) y’’ – y = 0, y(0) = 2, y’(0) = -3 b) (D2 –3D +2)y = 0, y(0) = -1, y’(0) = 0
c) (D3 –16D)y = 0, y(0) = 0, y’(0) = 0, y’’(0) = 16

3. Encuentre las soluciones generales de cada una de las siguientes ecuaciones:


a) (D2 – 4D + 4)y = 0 b) 16y’’ – 8y’ + y = 0 c) 4I’’(t) – 12I’(t) + 9I(t) = 0
d) (D6 – 4D4)y = 0 e) (D4 – 2D3 + D2)y = 0 f) 4y(IV) – 20y’’ + 25y = 0

4. Encuentre las soluciones que satisfagan las condiciones dadas:


a) (D2 – 2D + 1)y = 0; y(0) = 1, y’(0) = -2
b) (D3– D2)y = 0; y(0) = 1, y’(0) = y’’(0) = 0
c) 2 16  64s; s 0, 4 donde t 0
d 2s ds ds
dt dt dt

5. Encuentre la solución general de cada una de las siguientes ecuaciones:


c) 4 2 9s
2 d 2s
a) y’’ + 4y = 0 b) (D + 4D + 5)y = 0
dt
d) 4y’’ – 8y’ + 7y = 0 e) y(IV) = -16y’’ f) (D3 + D2 - 2)y = 0

6. Encuentre las soluciones que satisfagan las condiciones dadas.


a) (D2 + 1)y = 0; y(0)=4, y’(0)=0 b)U’’(t) = -16U(t); U(0)=0, U’(0) = 4
c) I’’(t) + 2I’(t) + 5I(t) = 0; I(0) = 2, I’(0) = 0

77
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Independencia lineal y Wronskianos

Definición: Llamaremos wronskiano de y1 e y2 al determinante:


y y2
W (y1, y2) = 1 = y1 y’2 – y2 y’1
y1 ' y 2 '

Teorema 1
Sean y1, y2 dos soluciones de la ecuación diferencial
(ao D2 + a1 D + a2) y = 0
ao z 0, a1, a2 funciones continuas de x en algún intervalo. Entonces
 ³ a1o ( x ) dx
a ( x)
y1 y2
W (y1, y2) = ce (identidad de Abel)
y '1 y'2

Teorema 2
Sean y1 , y2 soluciones de la ecuación diferencial
(ao D2 + a1 D + a2) y = 0 en algún intervalo I. Entonces
a) y1, y2 son linealmente dependientes (l.d.) en I œ W (y1, y2) = 0 en I.
b) y1, y2 son linealmente independiente (l.i.) en I œ W (y1, y2) z 0 en I.

Teorema 3
Sea y1 una solución de la ecuación diferencial (ao D2 + a1 D + a2) y = 0.
Entonces una solución l. i con y1 está dada por
 ³
³
a1
dx
ao
e
y2 = y1 dx
y12

Nota
Este teorema permite determinar una solución y2 conocida otra solución y1.

Ejemplo
Encontrar la solución general de x2 y’’ + xy’ – y = 0, si se sabe que y = x es una solución.
Respuesta
ao = x2 , a1 = x, y1 = x
³
³ x³ 2 x ³ 2 dx x ³ x dx
1

e  ln x x 1 x 2
dx

3
x
e
2
y2 = x 2
dx x
x x x 2x

2
1
? y = c1 x + c 2 ˜ š w(y1 , y2) = z 0, x
x
1
1 x
x 1 x
x2

78
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Teorema 4
Si y1 , y2 son soluciones l.i. de la ecuación diferencial (aoD2 + a1 D + a2)y = 0 e yp es una
solución particular de (a0 D2 + a1 D + a2) y = F(x), entonces y = c1 y1 + c2 y2 es una solución de
la ecuación complementaria e y = c1 y1 + c2 y2 + yp es una solución de la ecuación general.

Nota
Los teoremas anteriores se pueden generalizar al caso de n funciones y ecuación diferencial de
orden n.
(ao Dn + a1 Dn-1 + ….. + an ) y = 0 (*)

y1 ! yn
y '1 ! y'n
1) w (y1 ,....., yn) =
# #
( n 1)
y1 ! y n( n 1)
 ³ a1o dx
a

2) w (y1,....., yn) = c e identidad de Abel.

3) y1 ,...., yn son l.d. en I œ w = 0 en I.


y1 ,...., yn son l.i. en I œ w z 0 en I.
4) Si y1, ...., yn son soluciones l.i. de la ecuación diferencial ( ) e yp es una solución particular de
la ecuación general (ao Dn + .... + an) y = F(x) entonces y = c1 y1 + .... + cn yn es una solución de
la ecuación diferencial ( ) e y = c1 y1 + ..... + cn yn + yp es una solución de la ecuación
diferencial general.

Ejercicio:

w = 7 z 0 l.i.
Determine si las siguientes funciones son linealmente independientes
a) x + 2 , 2x – 3
2 2

w = 12 z 0 l.i.
b) sen x , cos x , 2 w = 0 l.d.
c) 1, x, x2 , x3.

79
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Solución particular de la ecuación diferencial I (D)y = f(x).

Se aplica a la ecuación diferencial I(D)y = F(x), donde F(x) contiene un polinomio, términos de
Método de los coeficientes indeterminados

la forma sen(rx), cos(rx), erx , donde r es constante o combinaciones de sumas y productos de


estos.

Caso a : si F(x) contiene la forma erx


Se trata de buscar una función y cuya derivada n-ésima produzca F(x). Luego yp = a erx

Ejemplo:
Resolver y’’ + 4y = 4 e2x
Respuesta:
Sea y = a e2x una posible solución
y’ = 2a e2x
y’’ = 4a e2x

reemplazando en la ecuación, se tiene


4a e2x + 4a e2x = 4 e2x,
luego 8a e2x = 4 e2x,
1
:. a=
2
?
1 2x
yp = e es la solución particular.
2

La ecuación complementaria y’’ + 4y = 0 tiene solución general yc = c1 cos 2x + c2 sen 2x

En efecto: m2 + 4 = 0 Ÿ m = r 2i raíces complejas conjugadas.

?la solución general de la ecuación dada es y = yc + yp = c1 cos 2x + c2 sen 2x +


e2x
2

Caso b : si F(x) contiene la forma sen(rx) o cos(rx)


Se trata de buscar una función y cuya derivada n-ésima produzca sen(rx) o cos(rx).
Luego
yp = a sen(rx) + b cos(rx)
______________________________________________________________________________

Ejemplo:
Resolver (D2 + 4D + 4) y = 6 sen(3x)
Respuesta:
Sea yp = a sen(3x) + b cos(3x)
yp’ = 3a cos(3x) - 3b sen(3x)
yp’’ = - 9a sen(3x) - 9b cos(3x).

80
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Reemplazando en la ecuación dada, se tiene

- 9a sen(3x) – 9b cos(3x) + 4 (3a cos(3x) – 3b sen (3x)) + 4 (a sen(3x) + b cos(3x)) = 6 sen(3x)

? (- 5a – 12b) = 6 š (12a – 5b) = 0


(-5a – 12b) sen(3x) + (12a – 5b) cos(3x) = 6 sen(3x).

 30  72
Resolviendo el sistema se tiene: a = , b=
 30  72
169 169
? yp = sen 3x + cos 3x solución particular
169 169

La Ecuación complementaria (D2 + 4D + 4)y = 0 tiene solución general yc = (c1 + c2 x) e-2x

En efecto: m2 + 4m + 4 = 0 Ÿ (m + 2)2 = 0 Ÿ m = -2 con multiplicidad 2.

? la solución general de la ecuación dada es :


 30  72
y = yc + yp = (c1 + c2 x) e-2x + sen 3x + cos 3x
169 169

Caso c: si F(x) contiene la forma de un polinomio de grado n.


Se trata de buscar una función y, cuya derivada n-ésima produzca F(x).
Luego yp = an xn + ….. + ao

Ejemplo:
Resolver (D2 + 4D + 9)y = x2 + 3x.
Respuesta:
Sea y = ax2 + bx + c
y’ = 2ax + b
y’’ = 2a.

Reemplazando en la ecuación dada, se tiene: 2a + 8ax + 4b + 9ax2 + 9bx + 9c = x2 + 3x.


Resolviendo el sistema se tiene:
9a = 1
8a + 9b = 3
2a + 4b + 9c = 0

1 19  94
Luego: a=,b ,c
9 81 729
? yp = x 2  x 
1 19 94
solución particular
9 81 729

La ecuación complementaria (D2 + 4D + 9) y = 0 tiene solución general


yc = e-2x (c1 cos 5 x + c2 sen 5 x)

81
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

En efecto:
m2 + 4m + 9 = 0 Ÿ m = - 2 r 5 i raíces complejas conjugadas

? la solución general de la ecuación dada es :


x  x
1 2 19 94
y = yc + yp = e-2x (c1 cos ( 5 x) + c2 sen ( 5 x) +
9 81 729

Ejercicio
Resolver (D2 + 2D + 1) y = 2 cos(2x) + 3x + 2 + 3 ex
Respuesta:
Sea yp = a sen (2x) + b cos( 2x) + cx + d + f ex
y’ = 2a cos (2x) – 2b sen (2x )+ c + f ex
y’’ = 4a sen( 2x) – 4 b cos( 2x) + f ex

Reemplazando en la ecuación dada.

(-3a - 4b) sen 2x + (4a - 3b) cos 2x + cx + d + 2c + 4 f ex = 2 cos 2x + 3x + 2 + 3 ex


? - 3a - 4b = 0
4a - 3b = 2
c=3 Ÿ a= , b  , c = 3, d = - 4, f =
8 6 3
25 25 4

? yp =
8 6 3
sen 2x - cos 2x + 3x – 4 + ex
25 25 4

La ecuación complementaria (D2 + 2D + 1)y = 0 tiene solución general yc = (c1 + c2 x ) e-x

m2 + 2m + 1 = 0 Ÿ (m + 1)2 = 0 Ÿ m = -1 con multiplicidad 2.


En efecto:

? la solución general de la ecuación dada es :


8 6 3
y = yc + yp = (c1 + c2 x) e-x + sen 2x - cos 2x + 3x – 4 + ex
25 25 4

El método aniquilador (Justificación al método de los coeficientes indeterminados)

Dada la ecuación diferencial I (D) y = F(x) se trata de buscar un operador tal que aplicado a la
ecuación dada, se obtenga una ecuación con el lado derecho igual a cero. El operador requerido
se llama operador de aniquilación o aniquilador.

Ejemplos
1) Resolver (D2 + 4) y = 4 e2x

Respuesta: (D2 + 4) y = 4 e2x /D


D (D2 + 4) y = 8 e2x

82
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

? 2 (D2 + 4) y = 8 e2x
D (D2 + 4) y = 8 e2x

Aniquilador: (D – 2)

(D – 2) (D2 + 4)y = 0 (m – 2) (m2 + 4) = 0 ecuación auxiliar


m = 2, m = r 2i
Restando:

y = c1 cos 2x + c2 sen 2x + c3 e2x solución general de

yc = c1 cos 2x + c2 sen 2x , yp = c3 e2x


1
Reemplazando yp en la ecuación dada, se tiene c3 =
2

2) Resolver (D2 + 4D + 4) y = 6 sen 3x


Respuesta:
Aniquilador: (D2 + 9)
(D2 + 9) (D2 + 4D + 4) y = (D2 + 9) (6 sen 3x)
(D2 + 9) (D2 + 4D + 4) y = 0
(m2 + 9) (m2 + 4m) = 0 Ÿ m = r 3i , m = -2 con multiplicidad 2

y = (c1 + c2 x) e-2x + c3 sen 3x + c4 cos 3x

 30  72
yp = c3 sen 3x + c4 cos 3x
Reemplazando yp en la ecuación dada.: c3 = , c4 =
169 169

3) Resolver (D2 + 2D + 1) y = 2 cos 2x + 3 x + 2 + 3 ex


Respuesta
para 2 cos 2x usamos D2 + 4
para 3x + 2 usamos D2
para 3ex usamos D - 1

? D2 (D – 1) (D2 + 4) (D2 + 2D + 1) y = 0
? y = (c1 + c2 x) e-x + c3 sen 2x + c4 cos 2x + c5 ex + c6 x + c7
yp = c3 sen 2x + c4 cos 2x + c5 ex + c6 x + c7

6
Reemplazando yp en la ecuación dada, se tiene
8 3
c3 = c4 = c5 = c6 = 3 c7 = - 4.
25 25 4

Ejercicios propuestos:

Resolver las siguientes ecuaciones diferenciales

1) (D3 + 4D) y = ex + sen x

83
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

2) s’’(t) – 3s’ (t) + 2s (t) = 8t2 + 12 e-t , s(o) = 0, s’(o) = 2


3) y’’ + y = 6 cos2x , y(0) = 0 y (S/2) = 0
4) y’’ – 3y’ + 2y = 4 sen3 3x
5) Demuestre el uso del método de aniquilación para
a) y’’ + y = 2 e3x
b) (D2 + 2D + 1) y = 4 sen 2x
c) (D2 – 4) y = 8x2

Excepciones en el método de los coeficientes indeterminados

Ejemplo:
Resolver (D2 + 3D + 2)y = 4 e-2x
a) Usando el método de los coeficientes indeterminados
b) Usando aniquiladores
Respuesta

ecuación auxiliar m2 + 3 m + 2 = 0 Ÿ m = - 1 , m = - 2
a) Usando el método de los coeficientes indeterminados

? yc = c1 e-ex + c2 e-x
sea y = a e-2x
y’ = - 2a e-2x
y’’ = + 4 a e-2x
reemplazando
4a e-2x + 3 (-2a e-2x ) + 2ae-2x = 4e-2x
4a e-2x - 6a e-2x + 2ae-2x = 4e-2x
0 = 4e-2x contradicción
Mejor considerar y = ax e , si y = a e-2x está en yc
-2x

(D2 + 3D + 2) (ax e-2x ) = - a e-2x = 4 e-2x Ÿ a = - 4


? y = c1 e-2x + c2 e-x - 4x e-2x

b) Por el método de aniquilación


(D2 + 3D + 2) y = 4 e-2x / (D + 2)
(D + 2) (D2 + 3D + 2) y = - 8 e-2x + 8 e-2x = 0

Ÿ y = c1 e-2x + c2 e-x + c3 x e-2x


(m + 2) (m + 2) (m + 1) = 0
, c3 = -4

Observación
Si un término de yp está en yc, se debe multiplicar el término de yp por una potencia mínima de x
para que esto no ocurra.

Ejercicios
1) Resolver (D2 + 4) y = 6 sen 2x + 3x2

ecuación auxiliar m2 + 4 = 0 Ÿ m = r 2i
Respuesta :

? yc = c1 cos (2x) + c2 sen (2x), f(x) = 6 sen(2x) + 3x2


como sen (2x) aparece en yc.
sea yp = x (a sen (2x) + b cos (2x) ) + cx2 + dx + f.

84
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

y’p = (a sen (2x) + b cos (2x) ) + x (2a cos (2x) - 2b sen (2x)) + 2cx + d
y’’p = 2a cos(2x) - 2b sen(2x) + 2a cos(2x) - 2b sen(2x) + x(-4a sen(2x) - 4b cos(2x)) + 2c
= 4a cos (2x) - 4b sen (2x) + x (-4a sen (2x) - 4b cos (2x)) + 2c
al sustituir en la ecuación diferencial queda:
4a cos(2x) - 4b sen(2x) + x(-4a sen(2x) - 4b cos(2x)) + 2c + 4x(a sen(2x) + b cos(2x)) + cx2 + dx
+ f) = 6 sen 2x + 3 x2

4a cos 2x – 4b sen 2x + 4cx2 + 4dx + 4f + 2c = 6 sen 2x + 3 x2


? 4a = 0 , - 4b = 6, 4c = 3 , 4d = 0 , 4f + 2c = 0
Ÿ a=0, b=-
3 3 3
, c= , d=0 f= -
2 4 8
? yp = -
3 3 2 3
x cos 2x + x -
2 4 8
? y = c1 sen 2x + c2 cos 2x - x cos 2x +
3 3 2 3
x -
2 4 8

2) Resolver (D3 – 3D2 + 3D – 1) y = 2 ex


Respuesta:
ecuación auxiliar (m3 – 3 m2 + 3m – 1) = 0
(m – 1)3 = 0
yc = (c1 + c2 x + c3 x2 ) ex
como ex está en yc y x ex está en yc y x2 ex está en yc , sea yp = a x3 ex
? (D3 – 3D2 + 3D – 1) (a x3 ex) = 6a ex = 2 ex Ÿ a =
1
3
? y = (c1 + c2 x + c3 x2 ) ex +
1 3 x
x e.
3

85
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Método de variación de parámetro

Dada la ecuación diferencial


y’’ + P(x) y’ + Q(x) y = F(x)

sabemos que tiene como solución general


y (x) = yc + yp

donde yc = c1 y1(x) + c2 y2(x) e yp es una solución particular

El método consiste en reemplazar las constantes c1 y c2 por las funciones v1 (x) y v2 (x), tal
que y = v1 y1 + v2 y2 sea una solución particular de la ecuación dada.

derivando : y’ = v1’ y1 + v1 y1’ + v2’ y2 + v2 y2’

ordenando : y’ = ( v1 y1’ + v2 y2’) + (v1’ y1 + v2’ y2)

? y’ = v1 y1’ + v2 y2’
para simplificar, hacemos v1’ y1 + v2’ y2 = 0

derivando y’
queda y’’ = v1 y1’’ + v1’ y1’ + v2 y2’’ + v2’ y2’

Reemplazando en la ecuación dada, se tiene

(v1 y1’’ + v1’ y1’ + v2 y2’’ + v2’ y2’) + P(x) (v1 y1’ + v2 y2’) + Q(x) (v1 y1 + v2 y2) = F(x).

Ordenado queda

? v1 ˜ 0 + v2 ˜ 0 + v1’ y1’ + v2’ y2’ = F(x)


v1 (y1’’ + P(x) y1’ + Q(x)y1) + v2 (y2’’ + P(x) y2’ + Q(x)y2) + v1’ y1’ + v2’ y2’ = F(x)

pues y1, y2 son soluciones de la ecuación diferencial complementaria.

Luego v1’y1’ + v2’y2’ = F(x)

y como agregamos la condición v1’y1 + v2’y2 = 0 se tiene el sistema

Ÿ
v1’ y1’ + v2’ y2’ = F(x) v1’ y1 + v2’ y2 = 0

v1’ y1 + v2’ y2 = 0 v1’ y1’ + v2’ y2’ = F(x)

0 y2 y1 0
F y2 ' y1 ' F
por Regla de Cramer v1’ = , v2’ =
w w

86
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

 y 2 F ( x) y1 F ( x)
de donde v1’ = , v2’ =
w( y1 , y 2 ) w( y1 , y 2 )

v1 = ³ v2 = ³ 1
 y 2 F ( x) y F ( x)
:. dx dx
w w

³ ³
 y 2 F ( x)
? yp = y1
y1 F ( x)
dx + y2 dx es la solución particular buscada.
w w

Ejemplo
Resolver y’’ + 4y’ + 4y = x e2x
Respuesta
Solución complementaria

m2 + 4m + 4 = 0 Ÿ (m + 2)2 = 0 Ÿ m = - 2 con multiplicidad 2.


y’’ + 4y + 4y = 0
?
yc = (c1 + c2 x) e-2x

Solución particular
sea y1 = e-2x , , y2 = x e-2x
-2x
y2’ = e-2x (1 – 2x)
? yp = v1 (x) y1 (x) + v2 (x) y2(x)
y1’ = -2e ,

e 2 x xe 2 x
 2e  2 x e  2 x (1  2 x)
w (y1 , y2) = = e-4x

³ e 4 x dx + x e ³ e 4 x dx
 xe 2 x xe 2 x e 2 x xe 2 x
? yp = e-2x -2x

yp = e-2x ³  x 2 e 4 x dx  xe 2 x ³ xe 4 x dx = e2x «  »
ªx 1º
¬16 32 ¼
ªx 1º
? y = yc + yp = e-2x (c1 + c2 x ) + e2x «  »
¬16 32 ¼

Ejercicios propuestos:

1) Encuentre la solución general de las siguientes ec. diferenciales con coef. constantes.
a) y’’ + 4y = tg 2x
b) y’’ – y’- 6y = e-x
c) y’’ + y = cosec x
d) y’’ + 2y’ + y = e-x ln x
e) y’’ – 2y’ – 3y = 64x e-x
f) y’’ + 2y’ + 5y = e-x sec 2x

2) Encuentre la solución general de las siguientes ecuaciones dif. con coef. variables.
a) (x2 – 1) y’’ – 2xy’ + 2y = (x2 - 1)2 , si y = x es solución de la ec. complementaria.
b) (x2 + x)y’’ + (2 – x2)y’ – (2 + x)y = x(x + 1)2, si y = ex es sol. de la ec. complementaria

87
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Observación
El método de variación de parámetros puede extenderse a ecuaciones diferenciales de orden
arbitrario.

Dada la ecuación diferencial de orden n.


y(n) + a 1 (x) y(n-1) + ............+ an(x) = F(x)

y supongamos que se conoce la solución complementaria.


yc = c1 y1(x) + ...... + cn yn (x)

se busca una solución particular de la forma.


yp = c1(x)y1(x) + ...... + cn(x) yn (x)

tal que
c1’ y1 + .................. + cn’ yn = 0

c1’ y1’ + .................. + cn’ yn’ = 0


#
#
c1‘ y 1( n  2) + ........... + cn’yn(n-2) = 0
c1’ y 1( n 1) + ...... + c n ' y n( n 1) = F(x)

sea W >y1, ....., yn@ = Wronskiano de y1, ....., yn

y sea Vk (x) = determinante obtenido de W al reemplazar su columna k por


§0·
¨ ¸
¨0¸
¨#¸
¨ ¸
¨0¸
¨1¸
© ¹

³
Vk ( x) ˜ F ( x)
? ck ( x)
Vk ( x ) F ( x )
Entonces c k' ( x) dx
w w

Luego: yp = y1(x) ³ dx  .......  y n ( x) ³ n


V1 ( x) F ( x) V ( x) F ( x)
dx es la solución particular buscada.
w w

Ejemplos
1) Hallar una solución particular yp para la ecuación 3 y’’’ + 5y’’ – 2y’ = ex

5 2
Respuesta
y ' ' ' y ' '
ex ex
y' , luego F(x) =
3 3 3 3

88
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Ecuación complementaria: 3y’’’ + 5y’’ – 2y’ = 0

Ecuación auxiliar: 3 m3 + 5 m2 – 2m = 0
m (3 m2 + 5 m- 2) = 0
3 m (m + 2) (m – 1/3) = 0

? yc = c1 + c2 e-2x + c3 ex/3

y1 = 1 , y2 = e-2x , y3 = ex/3

yp = y1( x ) ³ dx  y 2 ³ 2  y3 ³ 3
V1 ( x) F ( x) V ( x) F ( x)dx V ( x) F ( x)dx
w w w

e 2 x e 2 x
x x

1 e3 0 e3
5 x 5 x
x x

0  2e  2 x  e V1 = 0  2e  2 x
e3 14 3
e3 7 3
w = e
3 9 3 3
x x

4e  2 x 4e  2 x
e3 e3
0 1
9 9

e 2 x
3
1 0 e
x x
1 0
 V3 = 0  2e  2 x 2e 2 x
3 3
e e
V2 = 0 0 0
0 4e  2 x
3 3
x 1
3
e
0 1
9
5
x

 ˜  2e  2 x ˜
³ ³ ³
7 3x ex e3 ex ex
? yp = 1 3 dx  e  2 x 3 3 dx  e 3
e x

 14 3  14 3 x  14 3
3 3 dx
5 x 5 5 x

e e e
9 9 9

e  e  e
1 x 1 x 9 x 1 x
yp = - e
2 42 14 6

2) Resolver y’’’ – 2y’’ – y’ + 2y = ex


Respuesta:
Ec. complementaria (D3 – 2D2 – D + 2) y = 0
Ec. auxiliar m3 – 2m2 – m + 2 = 0 Ÿ m = 2, 1, - 1.

89
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

? yc = c1 e2x + c2 ex + c3 e-x

sean y1 = e2x , y2 = ex , y3 = e-x

ex
1  3
e2x ex 1 1 1 1 1 1
e x
e ˜ e e 2 1 1
x x
e 0 1  3 e 2 x (3  9)  6e 2 x
3 3
2x x 2x 2x
W = 2e e e2x
4e 2 x ex e x 4 1 1 0 3 3
0 ex ex
ex
e x
2
ex
 ex
x
V1 = 0 e
ex
ee
1 ex
e2x 0 ex
e x
0  e x  (e x  2e x )
e2x
 e x
V2 = 2e 2 x 3e x
e x
2e 2 x
4e 2 x 1

e2x ex 0
e 3 x  2e 3 x  e3x
e2x ex
V3 = 2e 2 x ex 0
2e 2 x ex
4e 2 x ex 1

³  6e 2 x ³  6e 2 x ³  6e 2 x dx
 2e x x  e e
 
3e x e x 3x x
yp = e2x dx e x
dx e

 e x x ex
 e x   ˜  ex  ex  ex
2e 2 x 1 x 1
yp =
6 2 6 3 2 12

yp =  e x  e x como ex  yc , yp = - e x
1 x x
4 2 2

?y = c1 e2x + c2 ex + c3 e-x -
x x
e
2

90
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Ejercicios propuestos de ecuaciones diferenciales lineales no homogéneas

Resuelva cada una de las siguientes ecuaciones diferenciales sujetas a condiciones dadas si las
hay.
2) 2  3  2 y senx
d2y dy
1) y’’ + 3y = x2 + 1; y(0) = 0, y’(0) = 2

3) (D + 2D + 1)y = e  e 4) y' ' ' - 4y 4x  2  3e -2x


dx dx
2 x x

 2  5 I 34 cos 2t  x 8e t
d 2I dI d 4x
5) 2
6)
dt 4
7) y ' '4 y xe 2 x ; y (0) y ' (0) 8) x 2 y ' '6 y 0; y (1)
dt dt
0 2, y ' (1) 0

1 10) y ( IV )  16 y ' ' 64 cos 4 x


d3y d2y
9) 2
dx 3 dx 2
11) y’’  4 y x(1  cos x) 2r  e  2I
d 2r
dI 2
12)

13) y’’’ - 4y’’ + 4y’ = 12 e 2 x  24 x 2 14) y’’ + y = secx, y(0) = 1, y’(0) = 2


16) 4  2 2  s 100 cos 3t
d 4s d 2s
15) x2y’’ – 4xy’ + 4y = 24(x+1)
dt dt
17) 4y’’ – 4y’ + y = ln x 18) D(D – 1)(D2 – 4)y = X2 – X + ex
2

19) 4  9 2 20e .t ; I (0)


d 4I d 2I ln x
I ' (0) 0 20) x2y’’’ – xy’’ + y’ =
dt dt x

91

También podría gustarte