Confiabilidad Operacional
Confiabilidad Operacional
Confiabilidad Operacional
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
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
73
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA
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.
Nota: El teorema proporciona sólo condiciones suficientes, es decir, aún si las condiciones
enunciadas no se satisfacen todas, pueden existir soluciones únicas.
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.
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.
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.
75
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA
1 ....
x x2 x3
ex
1! 2! 3!
Ejemplos
1) (D2 + 2D + 5) y = 0
m2 + 2m + 5 = 0
m=-1r 2i
y = e-x >c1 sen 2x + c2 cos 2x@
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:
77
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA
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
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
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
Ejemplo:
Resolver y’’ + 4y = 4 e2x
Respuesta:
Sea y = a e2x una posible solución
y’ = 2a e2x
y’’ = 4a e2x
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
30 72
Resolviendo el sistema se tiene: a = , b=
30 72
169 169
? yp = sen 3x + cos 3x solución particular
169 169
Ejemplo:
Resolver (D2 + 4D + 9)y = x2 + 3x.
Respuesta:
Sea y = ax2 + bx + c
y’ = 2ax + b
y’’ = 2a.
1 19 94
Luego: a=,b ,c
9 81 729
? yp = x 2 x
1 19 94
solución particular
9 81 729
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
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
? yp =
8 6 3
sen 2x - cos 2x + 3x – 4 + ex
25 25 4
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
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)
30 72
yp = c3 sen 3x + c4 cos 3x
Reemplazando yp en la ecuación dada.: c3 = , c4 =
169 169
? 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:
83
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA
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
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 :
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
85
UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA
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.
? 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’
(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’ y1’ + v2’ y2’ = F(x) v1’ y1 + v2’ y2 = 0
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
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.
tal que
c1’ y1 + .................. + cn’ yn = 0
³
Vk ( x) F ( x)
? ck ( x)
Vk ( x ) F ( x )
Entonces c k' ( x) dx
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 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
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
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
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
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
91