Cuaderno de Trabajo Ecuaciones Diferenciales Ordinarias
Cuaderno de Trabajo Ecuaciones Diferenciales Ordinarias
Cuaderno de Trabajo Ecuaciones Diferenciales Ordinarias
Ordinarias
Cuaderno de Ejercicios
Recopilación
ANEXO B. .............................................................................................................................................. 88
ANEXO C ............................................................................................................................................... 89
ANEXO D. .............................................................................................................................................. 91
Ecuación Diferencial y su Clasificación
Definición 1.1. Si una ecuación contiene las derivadas o las diferenciales de una o más variables
dependientes con respecto a una o más variables independientes, se dice que es una ecuación
diferencial (ED).
Clasificación de las ED
Las ecuaciones diferenciales se pueden clasificar según tres características: tipo, orden y
linealidad.
Definición 1.2. Si la ecuación contiene derivadas ordinarias de una o más variables dependientes
con respecto a una sola variable independiente entonces la ecuación se dice que es una ecuación
diferencial ordinaria (EDO).
dy
Ejemplo 1 3 + 4y =
5
dx
du dv
Ejemplo 3 u +u =
x
dx dx
Definición 1.3. Si la ecuación contiene derivadas parciales de una o más variables dependientes
con respecto a una o más variables independientes, se dice que es una ecuación en derivadas
parciales. (EDP)
∂u ∂u
Ejemplo 4 = −
∂y ∂x
∂ 2u
Ejemplo 5 = y−x
∂x∂y
Definición 1.4. El orden de una ecuación diferencial lo determina el orden de la más alta derivada
presente en ella.
dy y
Ejemplo 7 = , La ED es de orden 1
dx x
dny d n −1 y dy
an ( x) n
+ an −1 ( x ) n −1
+ .....a1 ( x) + a0 ( x) y =
g ( x)
dx dx dx
NOTA
i.- La variable dependiente y todas sus derivadas tienen exponente sólo pueden tener exponente 1.
d3y d2y dy
Ejemplo 8 x2 3
+ cos x 2
+ senx + x 2 y =
ex ED lineal de orden 3
dx dx dx
d2y dy
Ejemplo 9 y2 2
+ y + xy =
x ED no lineal ya que la variable dependiente, depende
dx dx
de ella misma
Clasifique las siguientes ecuaciones diferenciales de acuerdo a su tipo, orden, grado y linealidad.
− 0.1(1 − y 2 ) + 9 y = 0
7. d2y dy
dx 2 dx
∂3x ∂3x
8. + 2 y = y + 2c
∂c 3 ∂y 3
d 2ρ 4 dρ
2
9. = ρ +
dθ dθ
2
Clasifique las siguientes ecuaciones diferenciales de acuerdo a su tipo, orden, grado y linealidad.
3. yy ′ + 2 y = 1 + x 2
4. x 2 dy + ( y − xy − xe x )dx = 0
5. x 3 y ( 4 ) − x 2 y ′′ + 4 xy ′ − 3 y = 0
d2y
6. + 9 y = sen y
dx 2
2
dy d2y
7. = 1 + 2
dx dx
d 2r k
8. 2
=− 2
dt r
9. ( sen x) y ′′′ − (cos x) y ′ = 2
10. (1 − y 2 )dx + xdy = 0
Definición 1.6. Una función f, definida en algún intervalo I, es solución de una ecuación
diferencial en dicho intervalo, si al sustituirla en la ED la reduce a una identidad.
Las soluciones de las ecuaciones diferenciales pueden ser explícitas o implícitas. Una ED tiene,
generalmente, un número infinito de soluciones o más bien una familia n-paramétrica de
soluciones. El número de parámetros, n, depende del orden de la ED. Cuando se dan valores
específicos a los parámetros arbitrarios, es decir, cuando se asignan valores numéricos a los
parámetros, se obtiene una solución particular de la ED. En algunas ocasiones se tiene una solución
que no pertenece a la familia n-paramétrica, a tales soluciones se les llama singulares.
Derivación Implícita
dy dy cy
=1 ln(cy ) +
dx dx cy
dy dy
=1 ln(cy ) +
dx dx
dy
=1 (ln(cy ) + 1)
dx
dy 1
=
dx (ln(cy ) + 1)
y ln(cy ) + y
=
ln(cy ) + 1
y (ln(cy ) + 1)
= = y
ln(cy ) + 1
y=y
d2y
Ejemplo 11 φ (x ) = sen x + x es solución de
2
2
+ y = x2 + 2
dx
′( x) cos x + 2 x
φ=
φ ′( x) =
− s enx + 2
− s enx + 2 + senx + x 2 = x 2 + 2
x2 + 2 = x2 + 2
Definición 1.7 Una ED acompañada de unas condiciones iniciales se le llama un problema de valor
inicial (PVI.), con frecuencia es importante saber si un problema de valor inicial tiene solución y
también deseamos saber si esta solución es única, aunque no podamos conseguir explícitamente la
solución, el siguiente
Ejercicios
1
1. y = ce x de y ' − y = 0 2. y = 2e − 2 x + e x de y ' + 2 y = e x
3
64 x sen x
3. y = 8 ln x + c de y ' = 4. y = de xy ' + y = cos x
x3 3x
x = cos t
5. y = 8e x + xe x de y ' ' − 2 y ' + y = 0 6. de yy ' + 4 x = 0
y = 2sen t
1 3
7. y − = 0 de y ' − y tan x = 0 8. y = − de y' = 3 y 2
cos x 3x + 2
II. Muestre que φ(x) es una solución explícita de la ecuación diferencial dada.
1. φ ( x ) = 2x 3 ; 2. φ ( x ) = e x − x ; + y 2 = e 2 x + (1 − 2 x )e x + x 2 − 1
dy dy
x = 3y .
dx dx
d2y d2y
3. φ ( x ) = x 2 − x −1 ; x 2 = 2y 4. φ ( x ) = sen x + x 2 ; + y = x2 + 2
dx 2 dx 2
dy 1 dy 2 xy
1. y2 + x − 3 = 0; =− 2. y − ln y = x 2 + 1; =
dx 2y dx y − 1
3. xy 3 − xy 3 sen x = 1 ; dy e − xy − y
4. e xy + y = x − 1 ; =
dy (x cos x + sen x − 1) y dx e − xy + x
=
dx 3(x − xsenx )
5. x 2 − sen ( x + y ) = 1;
= 2 x sec( x + y ) − 1
dy
dx
Tarea 2
dy
1. 2 y ' + y = 0; y = e − x / 2 2. − 2 y = e 3 x ; y = e 3 x + 10e 2 x
dx
y ' + y = sen x;
3. 1 1 4. 2 xydx + ( x 2 + 2 y )dy = 0; x 2 y + y 2 = C1
y = sen x − cos x + 10e − x
2 2
( x 2 + y 2 )dx + ( x 2 − xy )dy = 0;
5. y
c1 ( x + y ) = xe
2 x
Variables Separables
dy g ( x)
Definición 2.1.Se dice que una ED de la forma: = es separable o de variables separables.
dx h( y )
∫=
h( y )dy ∫ g ( x)dx + c
Nota: la constante o parámetro C, a veces es conveniente escribirla de otra manera, por ejemplo,
múltiplos de constantes o logaritmos de constantes o exponenciales de constantes o si aparecen
varias constantes reunirlas en una sola constante.
dy
Ejemplo 1.- = e3 x + 2 y
dx
Solución:
dy
= e3 x + 2 y
dx
dy dy
Factorizando = e3 x e 2 y Separar variables = e3 x dx
dx e2 y
∫e dy = ∫ e3 x dx
−2 y
Integrar
e −2 y e3 x
− +c =
2 3
e −2 y e3 x
Por tanto la solución general es: + =
c
2 3
dy dy
1. =
x 3 (1 − y ) ; y (0) =
3 2. =
(1 + y 2 ) tan x ; y (0) =
1
dx dx
dy x2 dy 3 x 2 + 4 x + 2
3. = 4. = ; y (0) = −1
dx 1 − y 2 dx 2( y − 1)
dy y cos x x2
=
5. = ; y (0) 1 6. y′ =
dx 1 + 2 y 2 y (1 + x3 )
7. y′ + y 2 senx =
0 8. y′ = ( cos 2 x )( cos 2 2 y )
dy x − e − x
9. = 10. y′ =
(1 − 2 x) y 2 ; y (0) =
−1/ 6
dx y + e y
1− 2x
11. y′ = ; y (1) = −2 12. xdx + ye − x dy = 0 ; y (0) = 1
y
2x
13. y′ = ; y (0) = −2 14. y′ =
xy 3 (1 + x 2 ) −1/2 ; y (0) =
1
y + x2 y
π
15. y ' = e x cos 2 y, y (0 ) =
4
dx
1. = sen5 x 2. dx + e3 x dy =
0
dy
dy
3. ( x + 1) =
x+6 4. xy′ = 4 y
dx
dy y 3 dx x 2 y 2
5. = 6. =
dx x 2 dy 1 + x
dy
7. = e3 x + 2 y 8. (4 y + yx 2 )dy − (2 x + xy 2 )dx =
0
dx
dx
9. 4( x 2 + 1) ; x (π / 4 ) =
= 1 10. x 2 y ' =y − xy ; y (−1) =−1
dy
Definición 2.2. f ( x, y ) Es homogénea de grado n si existe un real n tal que para todo t:
f (tx, ty ) = t n ( x, y )
Ejemplo 1.
Siempre que se tenga una ED homogénea podrá ser reducida por medio de una sustitución
adecuada a una ecuación en variables separables.
( x + ye y / x )dx − xe y / x dy =
0 , con y (1) = 0
Solución.
M ( x, y )= x + ye y / x Homogénea de orden 1
xdx − x 2 eu du =
0
xdx = x 2 eu du
dx
= eu du
x
dx
∫ x
= ∫ eu du
ln =
x eu + c
ln=
x e y / x + c Solución General
Para hallar la solución particular que pasa por el punto y(1) = 0, sustituimos en la solución general
y obtenemos:
= e0/1 + c ⇒ 0 =1 + c ⇒ c =−1
ln1
ln=
x e y / x − 1 Solución Particular
Ejercicios
1. ( x + y )dx − ( x − y )dy =
0 2. xy=′ 2 x + 3 y
−y
7. y + 2 xe x ; y (1) =
xy′ = 0 8. y ( ln y − ln x − 1) dx + xdy
= 0 y ( e=
) 1
x+ y
11. y′ = 12. xdy − ydx =
ydy
x− y
Tarea 4
2xy x y
1. y′ = 2. y=′ +
x − y2
2
y x
y y dx y
3. e y / x + ; y (1) =
y′ = 0 4. =
x sen y sen + x
x x dy x
5. (x 3
+ y 3 ) dx − xy 2 dy =
0
∂f ∂f
=
Si z = f ( x, y ) entonces dz dx + dy
∂x ∂y
∂f ∂f
dz= 0= dx + dy
∂x ∂y
Definición 2.4. La forma diferencial M ( x, y )dx + N ( x, y )dy es una diferencial exacta en una
región R del plano XY si corresponde a la diferencial total de alguna función f ( x, y ) .
Si M ( x, y )dx y N ( x, y )dy son continuas y tienen derivadas parciales de primer orden continuas
en una región R del plano XY, entonces la condición necesaria y suficiente para que la forma
diferencial.
M ( x, y )dx + N ( x, y )dy
∂M ∂N
=
∂y ∂x
Método de solución:
∂M ∂N
i).- Comprobar que es exacta, es decir, verificar que, =
∂y ∂x
∂f
ii).- Suponer que M ( x, y ) = y luego integrar con respecto a x dejando a y como constante.
∂x
=
f ( x, y ) ∫ M ( x, y)dx + g ( y) (2.2)
∂f ∂
∂y ∂y ∫
= + g ′( y ) N ( x, y )
M ( x, y )dx =
Despejar
∂
∂y ∫
g ′( y ) N ( x, y ) −
= M ( x, y )dx (2.3)
∂ ∂ ∂N ∂ ∂
N ( x, y ) − ∫ M ( x, y )dx = − ∫ M ( x, y)dx =
∂x ∂y ∂x ∂x ∂y
∂N ∂ ∂ ∂N ∂
−
∂x ∂x ∂y ∫ M ( x, y )dx = − M ( x, y ) =
∂x ∂y
0
∂f
NOTA: en ii) se pudo haber comenzado por = N ( x, y )
∂y
Solución:
∂M
= 4 xy + e x
∂y ∂M ∂N
Paso i) de donde =
∂N x ∂y ∂x
= 4 xy + e
∂x
Paso ii)
f ( x= ∫ N ( x, y)dx + h=
( x) ∫ (2 x y + e x − 1)dy + h( x)
2
, y)
x 2 y 2 + ye x − y + h( x)
Paso iii)
∂f
= M
= 2 xy 2 + ye x
∂x
∂f
= 2 xy 2 + ye x + h′( x) ⇒ h′( x) = 0
∂x
Paso iv)
x 2 y 2 + ye x − y + C1 =C
x 2 y 2 + ye x − y =
C2 Solución General
Ejercicios
(
1. e y dx + xe y − 2 y dy =
0) ( ) ( )
2. 6 xy − 2 y 2 dx + 3 x 2 − 4 xy dy =
0
3 dy 3
3. 1 − + x + y = −1 4. 2 x + y 2 + 2 xyy′ =
0
y dx x
( ) ( )
5. y cos x + 2 xe y + senx + x 2 e y − 1 y′ =
0 6. ( 2 x + 3) + (2 y − 2) y′ =
0
y
7. + 6 x dx + ( ln x − 2 ) dy =
0 8. ( 2 x − y ) dx + ( 2 y − x ) dy
= 0 ; y (1)
= 3
x
Tarea 5
1. (2 y − 2 xy 3 + 4 x + 6)dx + (2 x − 3 x 2 y 2 − 1)dy
= 0 ; y (−=
1) 0
2. (2 + 6 xy )dx + (3 x 2 − 2 y ) dy =
0
3. (2 x − 5 y + 2)dx + (1 − 6 y − 5 x)dy =
0
xy
5. + 2 x dx + 1 + x=
2
dy 0 ; y=
(0) 6
1+ x
2
Sea la ED
M ( x, y )dx + N ( x, y )dy = 0
Si µ ( x, y ) es tal que
µ ( x, y ) M ( x, y )dx + µ ( x, y ) N ( x, y )dy = 0
1 2 1
Ejemplo: xdx + ydy es la diferencial de ( x + y 2 ) ya que d ( x 2 + y 2 ) =xdx + ydy
2 2
1 1 1 1 1
=µ = ;µ = ;µ = ;µ = ;µ , son factores integrantes.
y 2
x 2
xy x +y
2 2
ax + bxy + cy 2
2
∂M ∂N dµ dµ
µ − =
N =
−M
∂y ∂x dx dy
Notas
∂M ∂N
−
∂y ∂x dµ dµ
1.- Si = f ( x) , entonces µ f ( x) = , y por lo tanto f ( x)dx =
N dx µ
∂M ∂N
−
∂y ∂x
= g ( y ) , entonces µ = e ∫
g ( y )dy
2.- Similarmente, si
−M
Solución:
∂M
M ( x, y )= 2 xy 2 − 2 y ⇒ = 4 xy − 2
∂y
∂N
N ( x, y )= 3 x 2 y − 4 x ⇒ = 6 xy − 4
∂x
Luego
∂M ∂N
− =
−2 xy + 2
∂y ∂x
Por tanto
∂M ∂N
−
∂y ∂x −2 xy + 2 2(− xy + 1)
= =
−M −2 xy + 2 y 2 y (− xy + 1)
2
Luego
1
1 ∫ dy
g ( y ) = ⇒ F .I =µ ( y ) =e y =e =y
Ln y
1. 6 xydx + (4 y + 9 x 2 )dy = 0 µ ( x, y ) = y 2
1
2. − y 2 dx + ( x 2 + xy )dy = 0 µ ( x, y ) =
x2 y
4. y ( x + y + 1)dx + ( x + 2 y )dy = 0 µ ( x, y ) = e x
5. (2 y 2 + 3 x)dx + 2 xydy = 0 µ ( x, y ) = x
II. Determine un factor integrante para convertir la ecuación diferencial en exacta y resuelva.
( )
1. 3 xy + y 2 + ( x 2 + xy ) y′ =
0 2. (3 x 2 y + 2 xy + y 3 )dx + ( x 2 + y 2 )dy =
0
x
3. y′ = e 2 x + y − 1 4. dx + − seny dy =
0
y
5. ydx + (2 xy − e −2 y )dy =
0 ( )
6. e x dx + e x cot y + 2 y csc y dy =
0
( ) ( )
7. 4 x 3 / y 2 + (3 / y ) dx + 3 x / y 2 + 4 y ) dy =
0
6 x2 y dy
8. 3 x + + +3 = 0
y y x dx
( ) ( )
10. 3 y 2 − x dx + 2 y 3 − 6 xy dy = 0
1. 3 x 2 ydx + ydy =
0 2. x tan xdx − y cos xdy =0 ; y (0) =2
3. x 2 senxdx + xydy =
0 4. ( y + x + 2)dx + dy =0
( )
5. xy + y + y 2 dx + ( x + 2 y )dy =0
Definición 1.6
dy
a1 ( x) + ao ( x) y =
g ( x)
dx
dy
+ p( x) y =
Q( x)
dx
Donde:
a0 ( x) h( x )
=p( x) = y Q( x)
a1 ( x) a1 ( x)
y ′+ p ( x ) y =
Q( x)
Es:
ye ∫ ∫e
∫
p ( x ) dx p ( x ) dx
= Q( x)dx + c
dy
1. x − 4y =
x6e x
dx
Solución:
4 x6e x
p( x) =
− Q( x) =
x x
4 4
ye ∫ x e ∫ x x 5e x dx + c
− dx − dx
= ∫
ye −4ln x ∫e
−4ln x
= x 5e x dx + c
yx −4 ∫x
−4 5 x
= x e dx + c
yx −4
= ∫ xe dx + c
x
yx −4 = xe x − e x + c
y = x 5e x − x 4 e x + cx 4
Ejercicios
dv dy dy
1. ( 6 − 2µ v ) + v2 =
0 2. + 2y = 0 3. 3 + 12 y = 4
dµ dx dx
dy dy
4. + y = e3x 5. = y + ex 6. y′ + 3 x 2 y = x 2
dx dx
dx
7. x 2 y′ + xy = 1 8. xdy = ( xsenx − y )dx 9. = x+ y
dy
dy dy
1. + 5 y = 20, y (0) = 2 2. + y cot x = 2 cos x
dx dx
dy dy
3. cos x + y sen x = 1 4. (1 + e x ) + ex y = 0
dx dx
5. (1 + x) y′ − xy = x + x 2
dy
Definición 2.7n Una E.D. de la forma Q( x) y n con n ≠ 0 y n ≠ 1 , se llama una
+ P( x) y =
dx
E.D. de Bernoulli. Obsérvese que es una E.D. no Lineal
La sustitución w = y1− n convierte la E.D. de Bernoulli en una E.D. Lineal en W de primer orden.
dw
+ (1 − n) p ( x) w = (1 − n)Q( x)
dx
dy
Ejemplo: xy (1 + xy 2 ) =
1
dx
Solución:
dy 1 dx
= ⇒ = xy (1 + xy 2 )
dx xy (1 + xy )
2
dy
dx dx
= xy + x 2 y 3 ⇒ − xy = x 2 y 3 (1)
dy dy
La ecuación anterior tiene la forma de E:D de Bernoulli con variable dependiente x, y Con n=2.
dx dw
= − w−2
dy dy
Sustituir en (1)
dw
− w−2 − yw−1 =
y 3 w−2
dy
Multiplicamos por − w2
dw
+ yw =
− y 3 La E.D anterior es lineal en w (resolver por el Método de Lineales)
dy
2
y
1 −
=− y 2 + 2 − e 2
x
dy dy 1
1. x + y = x2 y2 2. x +y= 2
dx dx y
dy dy
3. = y ( xy 3 − 1) 4. x − (1 + x) y = xy 2
dx dx
dy dy
5. x 2 + y 2 = xy 6. 3(1 + x 2 ) = 2 xy ( y 3 − 1)
dx dx
dy dy y x
7. − y = ex y2 8. 2 = − ; y (1) =
1
dx dx x y 2
Tarea 8
3x 2 x
1. y′ = 2. y′ =
x3 + y + 1 x y + y3
2
cos x
3. xy′ + y =x4 y3 4. xy 2 y′ + y 3 =
x
3/2
dx 2 x
=
5. − x y= 2
; y (1) 1 6. (1 − 2 xy 2 )dy
= y 3 dx ; y (0)
= 1
dy y y
Modelos de Población
dP(t )
= kP (t ) P (t0 ) = p0
dt
Donde:
k Constante de Proporcionalidad
Ejemplo:
Un cultivo de bacterias tiene una cantidad de 200 individuos inicialmente. Después de una hora la
cantidad de bacterias creció habiendo 3/2 de la población inicial. Si el crecimiento bacteriano
responde a un modelo de crecimiento proporcional, ¿cuánto tiempo tardará la población en
triplicarse?
Solución
Condiciones
dP(t )
∫ P(t )
= k ∫ dt ⇒ ln( P (t )) =
kt + c ⇒ P(t ) =
ce kt
k ≈ 0.41
P(t ) = ce0.41t
P(t ) = 600
200e0.41t = 600
t = 2.67 hr
Ley de enfriamiento
La ley de enfriamiento de Newton dice que en un cuerpo que se esta enfriando, la rapidez con que
la temperatura T(t) cambia es proporcional a la diferencia entre la temperatura del cuerpo y la
temperatura constante T0 del medio que lo rodea.
dT (t )
= k (T (t ) − Tm ) T (t0 ) = T0
dt
Donde:
k Constante de Proporcionalidad
Se piensa hacer cubos de hielo y se vierte agua en un molde para este propósito. Si el agua se
encuentra a temperatura ambiente (aprox. 20ºC), y después de 10 minutos el agua del molde está a -
3ºC, sabiendo que el congelador enfría a una temperatura constante de -15ºC.
Solución
Condiciones
dT (t )
= k (T − (−15)) {T (0) = 20
dt
dT (t )
∫ T (t ) + 15 = k ∫ dt ⇒ ln(T (t ) + 15) =kt + c ⇒ T (t ) =ce kt − 15
T (10) = −3
T (=
t ) ce kt − 15 ⇒ −= 3 35e k (10) − 15
1 12
=k ln ≈ −0.107
10 35
=T (t ) ce −0.107 t − 15
T (t ) = 0
1.- 35e −0.107 t − 15 =
0
t = 7.91min utos
T (t ) = −13
2.- 35e −0.107 t − 15 =
−13
t = 26.74 min utos
Vaciado de Tanques
Consideraciones
1) Cantidad de agua que sale por el orificio = Cantidad de agua que desciende en el orificio.
2) El volumen que bajó en el cilindro es V = A∆h (con signo negativo por ser decrecimiento)
3) El volumen que sale por el orificio es V = B∆m , donde ∆m es la distancia que recorre el
agua durante ∆t segundos si el chorro saliera horizontalmente.
dm
4) v = es la velocidad instantánea de la caída del líquido.
dt
5) Tomaremos v = 2 gh en condiciones ideales (masa del agua=su energía cinética),
suponiendo que no hay pérdidas.
∆h ∆m
−A =
B
∆t ∆t
Cuando ∆t → 0
dh dm
Tenemos − A =
B
dt dt
dh
Por lo tanto − A =
Bv
dt
dh
Y −A =
B 2 gh
dt
Que ya es la ecuación diferencial del proceso, con la condición inicial de que h = h0 , cuando t=0.
dh B
= − 2 gdt
h A
Se obtiene
h = h0 , cuando t = 0
2 h0 = c
B
∴2 h =
− 2 gt + 2 h0
A
−0.0004
=
Entonces 2 h 2 gt + 2 6
2
2 h=
−0.0008854t + 4.8989795
h=
(−0.0004427t + 2.4494898) 2
Para calcular el tiempo necesario que se necesita para vaciar el recipiente, tomamos h=0
4.8989795
=
Entonces: t = 5533.07 seg
= ⇒ t 92.22
= min 1.53horas
0.0008854
Ejercicios
1. Un cultivo tiene una cantidad inicial N o de bacterias. Cuando t=1h, la cantidad medida de
3
bacterias es N o . Si la razón de reproducción es proporcional a la cantidad de bacterias
2
presentes, calcule el tiempo necesario para triplicar la cantidad inicial de los microorganismos.
4. La población de una comunidad crece con una tasa proporcional a la población en cualquier
momento. Su población inicial es 500 y aumenta el 15% en 10 años. ¿Cuál será la población
pasados 30 años?
´
5. En un trozo de madera quemada o carbón vegetal se determinó que el 85.5% de su C-14 se
había desintegrado. Con la información del ejercicio 2 determine la edad aproximada de la
madera. Éstos son precisamente los datos que usaron los arqueólogos para fechar los murales
prehistóricos de una caverna en Lascaux, Francia.
6. Un termómetro se saca de un recinto donde la temperatura del aire es 70°F y se lleva al exterior,
donde la temperatura es de 10ºF. Pasando ½ minuto el termómetro indica 50ºF. ¿Cuál es la
lectura cuando t=1 minuto?, ¿Cuánto tiempo se necesita para que el termómetro llegue a 15ºF?
7. Se aplica una fuerza electromotriz de 100 volts a un circuito en serie RC, donde la resistencia es
8. Una ecuación diferencial que describe la velocidad v de una masa m en caída sujeta a una
resistencia del aire proporcional a la velocidad instantánea es
dv
m = mg − kv
dt
En que k es una constante de proporcionalidad positiva.
a) Resuelva la ecuación sujeta a la condición inicial v(0)=v0
b) Calcule la velocidad límite (o terminal) de la masa.
c) Si la distancia s se relaciona con la velocidad por medio de ds/dt=v, deduzca una ecuación
explicita para s, si también se sabe que s(0)=s0.
1. Era el mediodía en un frío día de diciembre en Tampa: 16°C. el detective Taylos llegó a la
escena del crimen para hallar al sargento sobre el cadáver. El sargento dijo que había varios
sospechosos. El detective Taylor sacó un termómetro y midió la temperatura del cuerpo:
34.5°C. Luego salió a comer. Al regresar, a la 1:00PM, halló que la temperatura del cuerpo era
de 33.7°C. ¿en qué momento ocurrió el asesinato¡. [Sugerencia: la temperatura normal del
cuerpo es de 37°C].
2. El Sr. Pérez coloca al mismo tiempo dos tazas de café en la mesa del desayunador. De
inmediato vierte crema en su taza, con una jarra que estaba desde hace mucho en esa mesa. Lee
el diario durante cinco minutos y toma su primer sorbo. Llega la Sra. Pérez cinco minutos
después de que las tazas fueron colocadas en la mesa, vierte crema en la suya y toma un sorbo.
Suponga que la pareja agrega exactamente la misma cantidad de crema, ¿quién y porqué toma
su café más caliente?. Base su aseveración en ecuaciones matemáticas.
Fundamentos
Definición 3.1
Teorema 3.1
Definición 3.2
c1 f 1 ( x) + c 2 f 2 ( x) + 2 + c n f n ( x) = 0
Supóngase que cada una de las funciones f 1 ( x), f 2 ( x),2 , f n ( x) posee n − 1 derivadas al
menos. El determinante
f1 f2 2 fn
f1′ f 2′ 2 f n′
W ( f1 , f 2 ,2, f n ) =
( n −1) ( n −1) ( n −1)
f1 f2 2 fn
W ( y1 , y 2 , 2 , y n ) ≠ 0
y = c1 y1 ( x) + c 2 y 2 ( x) + 2 + c n y n ( x)
y = c1 y1 ( x) + c 2 y 2 ( x) + 2 + c n y n ( x) + y p
y p = y p1 ( x) + y p 2 ( x) + + y pk ( x)
Ejercicios
I. En los problemas siguientes compruebe si las funciones respectivas son linealmente independientes
o dependientes en (−∞, ∞) .
1. f 1 ( x ) = x, f 2 ( x ) = x 2 , f 3 ( x ) = 4 x − 3 x 2
2. f1 ( x) = 5, f 2 ( x) = cos 2 x, f 3 ( x) = sen 2 x
3. f1 ( x) = 0, f 2 ( x) = x, f 3 ( x) = e x
= =
4. f1 ( x) cos 2 x, f 2 ( x) sen 2 x
5. f1 ( x) = x, f 2 ( x) = x − 1, f 3 ( x) = x + 3
6. f 1 ( x ) = 1 + x, f 2 ( x ) = x, f 3 ( x ) = x 2
7.=
f1 ( x) e=
x
, f 2 ( x) e − x
II. En los problemas siguientes compruebe que las funciones dadas forman un conjunto fundamental de
soluciones de la ecuación diferencial en el intervalo indicado. Forme la solución general.
1. y ′′ − y ′ − 12 y = 0; e −3 x , e 4 x , (−∞, ∞)
2. y ′′ − 4 y = 0; cosh 2 x, senh 2 x, (−∞, ∞)
3. y ′′ − 2 y ′ + 5 y 0;
= e x cos 2 x, e x sen 2 x, (−∞, ∞)
4. x y ′′ − 6 xy ′ + 12 y 0; x , x ,
= (0, ∞)
2 3 4
Tarea 10
Halle el Wronskiano de las siguientes funciones y forme la solución General.
1. y′′ − 9 y 0; e3 x , e −3 x , (−∞, ∞)
2. x3 y′′′ − 2 xy′ + 4 y 0; x 2 , x 2 ln x,
= (−∞, ∞)
3. y′′ − 2 y′ + 5 y 0;
= e x cos 2 x, e x sen 2 x, (−∞, ∞)
4. 4 y′′ − 4 y′ + y 0;
= e x /2 , xe x /2 , (−∞, ∞)
5. y′′ − 2 y′ + 5 y 0;
= e x cos 2 x, e x sen 2 x, (−∞, ∞)
Second-Order Equations
y′′ + a1 y′ + a0 y =
0 (5.1)
m 2 + a1m + a0 =
0 (5.2)
nth-Order Equations
y ( n ) + an −1 y n −1 + ....... + a1 y′ + a0 y =
0 (5.3)
=
with constant coefficients aj ( j 0,1,..., n − 1) is
Caution!
Characteristic equations are only defined for linear homogeneous differential equations with
constant coefficients.
(m − m1 )(m − m2 ) =
0
Note!
The roots of the characteristic polynomial determine the solution of the differential equation.
Case 1. m1 and m2 both real and distinct. Two linearly independent solutions are e m1x and e m2 x ,
and the general solution is
=y c1e m 1 x + c2 e m 2 x
Case 2. m1 = m2 . Two linearly independent solutions are e m1x and xe m1x , and the general solution
is
=y c1e m 1 x + c2 xe m 1 x
Case 3. m1= α + i β , a complex number. Since a1 and a0 in 5.1 and 5.2 are assumed real, the
roots of 5.2 must appear in conjugate pairs; thus, the other root is m2= α − i β . Two linearly
independent solutions are e(α +iβ ) x , and e(α −iβ ) x , and the general complex solution is
Warning: The above solutions are not valid if the differential equation is not linear or does not
have constant coefficients. Consider, for example, the equation y′′ − x 2 y =
0 . The roots of the
characteristic equation are m = x and m = − x , but the solution is not:
Solved Problems
Since the roots m1 = 2 and m2 = −1 are real and distinct, so the general solution is
=y c1e 2 x + c2 e − x
=
The roots m1 m=
2 4 are real and equal, so the general solution is
=y c1e 4 x + c2 xe 4 x
The roots m1 =−2 + i 3 and m1 =−2 − i 3 are complex, so the general solution is
1. 2 y ′′ − 5 y ′ − 3 y = 0 2. y ′′ − 10 y ′ + 25 y = 0 3. y′′ + y =
0
4. 4 y′′ − 4 y′ + y′ =
0 5. y ′′ − y ′ − 6 y = 0 6. y ′′ − 3 y ′ + 2 y = 0
d2y dy d2y dy
7. 2
+ 8 + 16 y = 0 8. 2
− 10 + 25 y = 0 9. y′′ − 2 y′ + 3 y =
0
dx dx dx dx
d4y d2y
15. y (5) + 5 y (4) − 2 y′′′ − 10 y′′ + y′ + 5 y =
0 16. 16 + 24 + 9y =
0
dx 4 dx 2
Tarea 11
Resuelva cada ecuación diferencial sujeta a las condiciones iniciales.
1. 2 y′′ + 6 y′ + 2 y = 0
2. y′′ − 2 y′ + 5 y = 0
3. y′′′ − 3 y′′ + 9 y′ + 13 y =
0
4. y′′ + 6 y′ + 5=y 0, y (0)= 0, y′(0)
= 3
5. 2 y′′ − 2 y′ + y =0, y ( 0 ) =−1, y′ ( 0 ) =0
Bibliografía Complementaria del Tema:
Sea L(D)y = f(x) una E.D. lineal, no homogénea, de coeficientes constantes y de orden n; si f(x)
tiene una de las siguientes formas:
a) f(x) = k, k constante
b) f(x) = polinomio en x
Método de Solución
Ejemplo:
y′′ + 4 y′ − 2 y = 2 x 2 − 3 x + 6
y′′ + 4 y′ − 2 y =
0; Función _ Complementaria
2 x 2 − 3 x + 6; Función _ Particular
y′′ + 4 y′ − 2 y =
0
Resolver Por Formula General
m 2 + 4m − 2 =0
A = −1
5
B = −5 / 2 Por lo Que y p =− x2 − x −9
2
C = −9
5
∴ y c1e − (2+
= 6)x
+ c2 e( −2+ 6)x
− x2 − x − 9
2
Ejercicios
1. y ′′ + 4 y ′ − 2 y = 2 x 2 − 3 x + 6 2. y ′′ − y ′ + y = 2 sen 3 x
3. y ′′ − 5 y ′ + 4 y = 8e x 4. y ′′ − 2 y ′ − 3 y = 4 x − 5 + 6 xe 2 x
5. y′′′ − 4 y′ =x + 3cos x + e −2 x
6. y′′′ + 4 y′ = x ; y (0) = y′(0) = 0, y′′(0) = 1
1. y ′′ + 3 y ′ − 2 y = 6 2. y′′ − 10 y′ + 25 y = 30 x + 3
1
3. y ′′ + y ′ + y = x 2 − 2 x 4. y′′ + y′ −= (0) 0, y′=
2 y 2t ; y= (0) 1
4
Sea a2 ( x) y′′ + a1 ( x) y′ + a0 ( x) y =
h( x )
y′′ + p ( x) y′ + g ( x) y =
f ( x)
Donde:
a1 ( x) a0 ( x) h( x )
p( x) = g ( x) = f ( x) =
a2 ( x) a2 ( x) a2 ( x)
y1′′+ p ( x) y1′ + g ( x) y1 =
0
y2′′ + p ( x) y2′ + g ( x) y2 =
0
=
yh c1 y1 + c2 y2
yp =
u1 ( x) y1 ( x) + u2 ( x) y2 ( x) =
u1 y1 + u2 y2
y hallemos u1 ,u2
Luego
u1′ y1 + u2′ y2 =
0
Luego,
y′′ + p ( x) y′ + g ( x) y =
f ( x)
Luego,
En resumen;
y1u1′ + y2u2′ =
0
y1′u1′ + y2′u2′ =
f ( x)
0 y2 y1 0
f ( x) y2′ y f ( x) y1′ f ( x) y1 f ( x)
u1′ = = − 2 =u2′ =
y1 y2 W ( y1 , y2 ) y1 y2 W ( y1 , y2 )
y1′ y2′ y1′ y2′
y′′ + p ( x) y′ + g ( x) y =
f ( x)
Asociada: y′′ + p ( x) y′ + g ( x) y =
0
2.- Hallamos W ( y1 , y2 )
y2 f ( x ) y1 f ( x)
3.- Hallamos u '1 = − u '2 =
W ( y1 , y2 ) W ( y1 , y2 )
Ejemplo:
y′′ + 3 y′ + 2 y =
sen(e x )
Solución:
y′′ + 3 y′ + 2 y =
0
m 2 + 3m + 2 =0
1.-
=yh C1e −2 x + C2 e − x
−2 x
=y1 e= , y2 e − x
e −2 x e− x
=
2.- W ( y1 , y2 ) = −2 x −2 x
e −3 x
−2e −e
y f ( x) e − x sen(e x )
u1′ =
− 2 =
− =
−e 2 x sen(e x )
W ( y1 , y2 ) e −3 x
3.-
′ y1 f ( x) e −2 x sen(e x )
u2 == − =
e x sen(e x )
W ( y1 , y2 ) e −3 x
u2 = ∫ u2′
= ∫ e sen(e )dx = − cos(e x )
x x
=
y p u1 y1 + u2 y2
=e x cos(e x ) − sen(e x ) e −2 x − e − x cos(e x )
5.- = −e −2 x sen(e x )
=
y yh + y p
= C1e −2 x + C2 e − x − e −2 x sen(e x )
Ejercicios
1
1. y ′′ + 4 y ′ + 4 y = ( x + 1)e 2 x 2. y ′′ + 9 y = csc 3 x
4
3. y′′ −=
y sec3 x − sec x 4. y ′′ + y = sec x
5. y ′′ + y = tan x 6. y ′′ + y = sec x tan x
7. y ′′ + y = cos 2 x 8. y′′ + 2 y′ + y =e − x ln x
ex
11. y′′ − 2 y′ + y = 12. y′′ + 4 y′ + 5 y =
e −2 x sec x
x
3. y′′ + 4 y′ =
tan 2 x 4. y′′ − 3 y′ + 2 y =
e x senx
x
5
5. y′′ − y′ + y =2
e cos x
4
=
Con la sustitución z ln=
x o x e x convertimos la ecuación anterior que es de coeficientes
variables en una de coeficientes constantes.
dz 1 dy dy dz
z =ln x ⇒ = ⇒ =
dx x dx dz dx
dy 1
=
dz x
dy dy
⇒x =
dx dz
d dy d dy
x =
dx dx dx dz
d 2 y dy d dy dz d 2 y 1
x= + =
dx 2 dx dz dz dx dz 2 x
d2y dy d 2 y
x2 + x =
dx 2 dx dz 2
2
2 d y d 2 y dy
luego x 2
= 2
− = Dz2 y − Dz y = Dz ( Dz − 1) y
dx dz dz
Similarmente
d3y
x3 = Dz ( Dz − 1)( Dz − 2) y
dx3
dny
y en general , x n = Dz ( Dz − 1)( Dz − 2).....( Dz − (n − 1)) y
dx n
Donde z = ln x
Solución:
z= ln x ⇒ Dz ( Dz − 1) y + Dz y + y= sec z
Dz2 y − Dz y + Dz y + y =
sec x
Dz2 y + =
y sec x ⇒ ( Dz2 + 1) =
y sec x
Ec. característica
m 2 + 1 =0 ⇒ m =±i
∴=
yh c1 cos z + c2 senz
Ejercicios
I. Resuelva
d2y dy d2y dy
1. x2 2
− 2x − 4 y = 0 2. 4x 2 2
+ 8x + y = 0
dx dx dx dx
d2y dy d3y 2
2 d y dy
3. x
2
2
+ 3x + 3y = 0 4. x
3
3
+ 5x 2
+ 7x + 8y = 0
dx dx dx dx dx
5. x 2 y ′′ − 3 xy ′ + 3 y = 2 x 4 e x 6. x 2 y ′′ + xy ′ + 4 y = 0
7. 25 x 2 y ′′ + 25 xy ′ + y = 0 8. x 2 y ′′ + 5 xy ′ + 4 y = 0
2 2 x
9. x3 y′′ − xy′ + y =x ln 3 x 10. y′′ − y′ + 2 y =
x x 1+ x
Tarea 14
5. x 2 y′′ + xy′ + 4 y =
sen(ln x)
Supongamos que tenemos un resorte dispuesto en forma vertical, con el extremo superior fijado a
una superficie horizontal y el otro extremo libre al cual se le fija una carga de masa m, la fuerza de
recuperación del resorte esta dada por la Ley de Hook:
F = ks
Donde:
s: elongación de resorte
Deduzcamos la E.D. que rige el movimiento de la masa sujeta al resorte ver las siguientes figuras.
d 2x
m 2 =− F + mg =−k ( x + s ) + mg
dt
− F + mg =−kx − ks + mg
d 2x d 2x
m =− kx ⇒ m + kx =0 ; si dividimos la ecuación por m
dt 2 dt 2
d 2x k
+ x=
0
dt 2 m
k d 2x
Llamemos =w2 ⇒ 2 + w2 x =0 E.D de Segundo Orden con Coeficientes Cts.
m dt
Ecuación Característica:
m 2 + w2 =⇒
0 m =− w2 =± wi
Solución General
=
x(t ) C1 cos wt + C2 s enwt
Condiciones Iniciales:
x(0) = x0
Si x0 = 0 : el cuerpo está en la PE
x′(0) = v0
C1
tan φ = Donde φ es el ángulo de fase
C2
C1 C
senφ = , cos φ = 2
A A
A ( C1 cos wt + C2 s enwt ) C C
=x(t ) = A 1 cos wt + 2 s enwt
A A A
2π
T= Periodo de vibraciones Libres
w
1
f = Frecuencia
T
Movimiento Amortiguado
Cuando el cuerpo sujeto al resorte se mueve en un medio que produce fricción sobre el cuerpo,
entonces decimos que el movimiento se efectúa con amortiguamiento Por la segunda ley de
Newton tenemos:
d 2x dx
m 2
=−k ( x + s ) + mg − β
dt dt
d 2x dx
m 2 +β + kx =
0 Dividiendo por m
dt dt
d 2x dx
2
+ 2λ + w2 x =
0
dt dt
Donde:
k β
= w 2 y 2λ =
m m
Ecuación Característica
p 2 + 2λ p + w 2 =
0
−2λ ± 4λ 2 − 4 w2 −2λ ± 2 λ 2 − w2
p1,2 = = =−λ ± λ 2 − w2
2 2
= (
x(t ) Ae − λt sen w2 − λ 2 t + φ )
Por Lo tanto lim x(t ) = 0
t →∞
Ahora supongamos que el cuerpo sujeto al resorte esta en un medio que produce fricción y
también hay una fuerza exterior F(t) que actúa sobre el cuerpo. Por la segunda ley d e Newton:
d 2x dx
m 2 = −k ( x + s ) + mg − β + F (t )
dt dt
d 2x dx
m 2 = − kx − ks + mg − β + F (t )
dt dt
d 2x dx
m 2
+β + kx =
F (t ) EC. De segundo Orden de Coef. Constantes no Homogéneas
dt dt
d 2x
1.- + 16 x =
0
dt 2
Condiciones Iniciales
dx
x(0) = 10 , t =0 =0
dt
Solución
Una formulación equivalente del problema es: se tira hacia debajo de un cuerpo que pende de un
resorte hasta que esté 10 unidades bajo la posición de equilibrio y lego se le retiene hasta t = 0 ; se
le suelta a continuación de manera que parta de un estado de reposo. Aplicando las condiciones
iniciales a la solución:
=
x(t ) C1 cos 4t + C2 sen 4t
= 10
Resulta x(0) = C1.1 + C2 .0
=
x(t ) 10 cos 4t + C2 sen4t
x′(t ) =
−40s en 4t + 4C2 cos 4t
dx
= 0= 4C2 .1
t =0
dt
x(t ) = 10 cos 4t
La solución muestra claramente que una vez que el sistema se pone en movimiento, permanece en
tal estado, con la masa desplazándose alternadamente 10 unidades hacia cada lado de la posición de
equilibrio x=0, lo cual se muestra en la siguiente figura:
=A C12 + C2 2
A= 102 + 02= 10
C2
cos φ =
A
C2 10
φ
cos= = = 1
A 10
2π 2π π
=
T = = Periodo de vibraciones Libres
w 4 2
1 2
f= = Frecuencia
T π
Datos
6 1
=
6 pu lg = pie Elongación que sufre el resorte con la masa aplicada
12 2
8 2
=
8 pu lg = pie Distancia que recorre el la masa bajo la posición de equilibrio
12 3
FR = ks ;
1 lb
2lb = k pie lo que implica que k = 4 , por lo que nos lleva a las siguientes ecuaciones:
2 pie
d 2x 1 d 2x
m 2 + kx =0 ⇒ + 4 x =0 si se multiplica por 16 se tiene
dt 16 dt 2
d 2x
+ 64 x =
0
dt 2
2 dx 4
x(0) = , t =0 = − pie / s
3 dt 3
Ec. Auxiliar
m 2 + 64 =
0
=
x(t ) C1 cos8t + C2 s en8t y aplicando las condiciones iniciales;
2 2
x(0)= = C1 cos 0 + C2 sen0 Por lo que C1 =
3 3
2
=
x(t ) cos8t + C2 s en8t Derivando y aplicando las condiciones iniciales;
3
4 16 1
x′(t ) =− =− sen0 + 8C2 cos 0 por lo que C1 = −
3 3 6
2 1
=
x(t ) cos8t − s en8t
3 6
=A C12 + C2 2
2 2
2 1
=
A +− =
4 / 9 + 1/ 36 ≈ 0.687
3 6
C2
cos φ =
A
−1/ 6
cos φ =
17 / 6
C1 2/3
φ
s en= = =
A 17 / 36
C1 2/3
φ
tan= =
C2 17 / 6
2π 2π π
=
T = = Periodo de vibraciones Libres
w 8 4
1 4
f= = Frecuencia
T π
d 2x dx
1.- 2
+ 5 + 4x =
0
dt dt
Condiciones Iniciales
Solución
El problema puede ser interpretado como una representación del movimiento sobreamortiguado de
una masa sujeta a un resorte. La masa parte desde una posición que se encuentra 1 unidad debajo
de la posición de equilibrio, con una velocidad dirigida hacia debajo de 1 m/s.
5 − t 2 −4t
x=
(t ) e − e
3 3
Para trazar la gráfica de x(t), encontramos el valor de t para el cual la primera derivada (velocidad)
es cero, si derivamos
5 − t 2 −4t
x=
(t ) e − e
3 3
5 8
x′(t ) =
− e − t + e −4t
3 3
8 1 8
e 3t = ⇒= t ln= 0.157 , con el estudio de la derivada en un entorno del punto, se deduce
5 3 5
que x(0.157) = 1.069m es un máximo. En otras palabras alcanza un desplazamiento máximo de
1.069m debajo de la posición de equilibrio.
Para saber si la gráfica cruza el eje t, esto es saber si la masa pasa por la posición de equilibrio. Se
hace
2 1 2
x(t ) = 0 , o sea e3t = , tiene la solución físicamente irrelevante t = ln = −0.305
5 3 5
2.- Un cuerpo que pesa 8 lb estira un resorte 2pie. Suponiendo que una fuerza de amortiguación
numérica igual a dos veces la velocidad instantánea actúa sobre el sistema y que el peso se suelta
desde la posición de equilibrio con una velocidad dirigida hacia arriba de 3 pie/s, determine la
ecuación del movimiento.
Datos
FR = ks ;
lb
8lb = k ( 2 pie ) lo que implica que k = 4 , por lo que nos lleva a las siguientes ecuaciones:
pie
d 2x dx 1 d 2x dx
m 2
+ β + kx =0 ⇒ 2
+ 2 + 4 x =0 si se multiplica por 4 se tiene
dt dt 4 dt dt
d 2x dx
2
+ 8 + 16 x =
0
dt dt
dx
x(0) = 0 , t =0 = −3 pie / s
dt
Ec. Auxiliar
m 2 + 8m + 16 =0 ⇒ (m + 4) 2 =0
0 = C1e0 + C2 0e0 ⇒ C1 = 0
1 d 2x dx
2
+ 1.2 + 2 x =
5cos 4t
5 dt dt
1 dx
x(0) = , t =0 =0
2 dt
Solución
El sistema tiene una masa de m = 1/ 5kg , sujeta a un resorte que tiene una constante de restitución
de k = 2 N / m . La masa se suelta, a partir del reposo, desde un punto que está ½ unidad (metros)
bajo la posición de equilibrio. El sistema tiene un amortiguador β = 1.2 y es expulsado por una
fuerza externa de f (t ) = 5cos 4t .
1 d 2x dx d 2x dx
2
+ 1.2 + 2 x =0 ⇒ 2
+ 6 + 10 x =0
5 dt dt dt dt
=
x p (t ) A cos 4t + Bsen 4t
x′p (t ) =
−4 Asen 4t + 4 B cos 4t
x′′p (t ) =
−16 A cos 4t − 16 B s en 4t
25 50
) e −3t (C1 cos t + C2 s ent ) −
x(t= cos 4t + sen 4t , aplicando condiciones iniciales se tiene:
102 51
Ejercicios:
2. Se suspende una bola de Acero que pesa 128 lb de un resorte, por lo cual el resorte se
alarga 2 pies de su longitud natural. Se pone en movimiento la bola sin velocidad inicial
desplazándola 6 pulgadas sobre la posición de equilibrio. Asumiendo que no hay
resistencia del aire, hallar la posición de la bola para t = π /128 seg , la frecuencia natural y
el período.
Tarea 15
Resuelva los siguientes problemas
1. Se suspende una masa de ¼ de slug de un resorte que tiene una constante de 1 lb/pie ,
la masa se pone en movimiento por un desplazamiento inicial de 2 pies en la dirección
hacia abajo y dándole una velocidad inicial de 2 pies/seg, en la dirección hacia arriba.
Halle el movimiento resultante de la masa, si la fuerza debida a la resistencia del aire es
•
−1x lb .
Definición Básica:
d d d
= [α f ( x) + β g ( x) ] = α f ( x) + β f ( x)
dx dx dx
b b
∫ [abab
f ( x) + g ( x) ]dx = ∫ f ( x)dx + ∫ g ( x)dx
a a
Siempre que cada integral exista. Por lo tanto la integración definida es una operación lineal.
Específicamente se estará interesado en una integral impropia que trasforma una función f (t ) en
una función del parámetro s.
DEFINICIÓN 7.1
∞ b
Ejemplo:
Calcular:
L {1}
∞ b
L {1} ∫= lim ∫ e − st dt
− st
= e (1)dt
b →∞
0 0
−e − st b
−e − sb + 1
lim = lim
b →∞ s b →∞ s
0
1
= siempre que s>0.
s
∞
NOTA: El uso del símbolo de límite se vuelve algo tedioso, por eso se adoptará la notación para
0
b
indicar lim(?)
b→∞ 0
Por ejemplo;
L {1} = ∫ e− st dt
0
−e − st ∞
1
=
s 0 s
Solo si s>0.
− st
Donde se da por entendido que el límite superior, cuando t → ∞ , e → 0 para s > 0
∫ e [a f (t ) + β g (t )]=
dt a ∫ e f (t )dt + β ∫ e − st g (t )dt
− st − st
0 a a
L [α f (t ) + β g (t )] =
α L { f (t )} + β L { g (t )}
= α F (s) + β G (s)
Ejercicios
I. Obtener la transformada de Laplace mediante la definición
1. L {1} 2. L {t} 3. L {t
2
} 4. L {e −3t }
−2t
5. L {sen 2t} 6. L {3t − 5 sen 2t} 7. L {te } 8. L {e
at
cos kt}
2 −2 t
9. L {t e } 10. L {tsenkt}
1. L {(t
2
+ 2t + 1)e −3t } 2. L {(t − 2) 2 } 3. L {(et + e − t ) 2 }
1 −2 t
7. L {Sen 3+ Cos 3t} 8. L {e Cos 4t} 9. L {t
5
+ cos 5t}
3
( t 2 + 2 )2
10. L
11. L
(
t + et )
2
13. L {tCosh9t}
e
2t 2t
e
14. L
(
et + e − t ) cos 2t
2
15. L {e 5t
+ cos 4t }
−t
e
Tarea 16
3 at
1. L {t } 2. L {cos kt} 3. L {e }
1. L {e −3t + t 3 − 2} 2. L {t
3 −2 t
e } 3. L {e −2t sen3t}
−2 t
5. L {3sen 4t + e
t
4. L {e cosh 2t} }
L
−1
{F ( s)} = f (t )
Nota:
1, si t ≥ 0 y t ≠ 1, t ≠ 2
=f (t ) =
3, si t 1
−3, si t =
2
−1
Para funciones continuas L es un operador lineal.
L
−1
[α F ( s) + β G ( s)] =
α L −1{ F ( s )} + β L
−1
{G ( s)}
−1
7s −1 A
−1 B C
L = L + +
( s − 3)( s + 2 )( s − 1) ( s − 3) ( s + 2 ) ( s − 1)
−1 1 −1 1 −1 1
=AL + B L +C L
s − 3 s + 2 s − 1
7s −1 A B C
= + +
( s − 3)( s + 2 )( s − 1) ( s − 3) ( s + 2 ) ( s − 1)
−1
7s −1 −2 t
= 2e − e − e
3t t
L
( s − 3)( s + 2 )( s − 1)
−1
s + 1 −1
A B C D E
L 2 3
=L 2+ + + +
s ( s + 2 ) s s ( s + 2 )3 ( s + 2 )2 s + 2
−1 1 −1 1 −1
1 −1
1 −1 1
=AL 2 + B L +C L 3
+D L 2
+E L
s s ( s + 2 ) ( s + 2 ) s + 2
t 2 e −2t te −2t
=At + B(1) + C +D + Ee −2t
2! 1!
s +1 A B C D E
= + + + +
s ( s + 2) s ( s + 2) ( s + 2) s + 2
2 3 2 3 2
s
1 1 1 1
A ==, B − , C=
− , D=
0, E=
8 16 4 16
−1
s + 1 1 1 1 t 2 e −2t 1
L 2 3
= t − − + Ee −2t
s ( s + 2 ) 8 16 4 2! 16
1 1
2 3
-1 -1
3. L 4. L
s ( s + 1) s + 4 s + 3s
2
s+3 s +1
2 2 2
-1 -1
5. L 6. L
( s + 4) s + 6 s + 25
s −5 2s − 1
2 2
-1 -1
7. L 8. L
s + 6 s + 13 s + 4 s + 29
3s + 1 s2
4 3
-1 -1
9. L 10. L
( s + 1) ( s + 2)
s−2 2s − 1
2 2 2
-1 -1
11. L 12. L
s ( s + 4) s ( s + 1)
7 3
3 2 2
-1 -1
13. L 14. L
s ( s − 1) s ( s − 9)
s+4
2 2
-1
15. L
s ( s + 16)
Tarea 17
Calcular la Trasformada Inversa de las siguientes funciones
s 2s + 4
2
-1 -1
1. L 2. L
s + 2s − 3 ( s − 2)( s + 4 s + 3)
2
s 2 + 6s + 9 1
2 2
-1 -1
3. L 4. L
( s − 1)( s − 2)( s + 4) s ( s + 4)
1
2
-1
5. L
s + 6 s + 13
Bibliografía Complementaria del Tema:
L {f n
(t )} (=
s ) s n F ( s ) − s n −1 f (0) − s n − 2 f ′(0) − .... − sf n − 2 (0) − f n −1 (0) en donde
F ( s ) = L { f (t )}
Demostración:
∞
L { f ′(t )} ( s) = ∫0 e − st f ′(t )dt ,
∞ − st
∞
= e − st
f (t ) 0
+ s∫ e f (t )dt
0
− f (0) + s L { f (t )} s
=
= sF ( s ) − f (0)
Demostración
∞
L { f ′′(t )} ( s) = ∫0 e − st f ′′(t )dt ,
∞ − st
∞
= e − st f ′(t ) 0
+ s∫ e f ′(t )dt u = e − st du = − se − st dt
0
∞
= ′(t ) dv
v f= ∫0
f ′′(t )dt
∞ − st
− f ′(0) + s ∫ e
= f ′(t )dt
0
− f ′(0) + s ( sF ( s ) − f (0) )
=
= s 2 F ( s ) − sf (0) − f ′(0)
∞
L { f ′′′(t )} ( s) = ∫0 e − st f ′′′(t )dt ,
∞ − st
∞
= e − st f ′′(t ) 0
+ s∫ e f ′′(t )dt u = e − st du = − se − st dt
0
∞
= ′′(t ) dv
v f= ∫0
f ′′′(t )dt
∞ − st
− f ′′(0) + s ∫ e
= f ′′(t )dt
0
Ejemplos:
L {− sent} =
s L ( cos t ) − f (0)
ζ {− sent} + 1
L {− sent} =
s L ( cos t ) − 1 ⇒= L ( cos t ) =
s
1 1 1 1 1
=
1 1
+ L {− sent} = + − 2 = −
s s s s s + 1 s s ( s 2 + 1)
1 1 1 s2 + 1 −1 1 s2 s
=1 − 2 = = =
s s + 1 s s2 + 1 s s2 + 1 s2 + 1
dy
=1.- + 3 y 13sen=
2t y (0) 6
dt
2.- y′′ − 3 y′ +=
2 y e −4t y= ′(0) 5
(0) 1, y=
dy
=3.- − y 2 cos
= 5t y (0) 0
dt
4.- x′′ + 16
= x cos 4t (0) 1, x′=
x= (0) 1
5.- y′′ − 4 y′ +=
4 y t 3e 2 t = ′(0) 0
y (0) 0, y=
7.- y′′ − 6 y′ =
+ 9y t e = ′(0) 6
y (0) 2, y=
2 3t
10.- y′′=
− y′ e cos t = y′(0) 0
y (0) 0, =
t
Tarea 18
Resuelva las siguientes E.D con la Trasformada de Laplace
1.- y′′ − 6 y′ +=
9y t (0) 0, y′=
y= (0) 1
2.- ′′ + y sen 2t
y= = y′(0) 1
y (0) 2, =
L {e at
f (t )} = L { f (t )} ( s − a)
= F ( s − a)
Donde:
F (s) = L { f (t )}
Demostración:
∞ ∞
L {eat f (t )} = ∫ e e f (t )dt = ∫ e
− st at − ( s −a )t
f (t )dt
0 0
∞ ∞
∫=
e e f (t )dt ∫ e
− st at − ( s −a )t
= f (t )dt
0 0
L {e at
t )} F ( s − a )
f (=
L {e at
f (t )} = L { f (t )}s→s −a
Ejemplo:
Calcular L {e t } =
5t 3
L {e t } = L {t 3}s→s −5
5t 3
1.- L te { }=10t
{
2.- L t e
3 −2t
}= {
3.- L t (e + e )
t 2t 2
}=
{
4.- L e sen3t =
t
} {
5.- L e senh3t =
5t
} {
6.- L te
−2t
}=
{
7.- L e (t + 3) =
t
} {
8.- L e
−2 t
sen 4t} = {
9.- L e cosh 5t =
4t
}
10.- L e { −2 t
cos 2t} =
Tarea 19
1.- L t e{ 3 5t
}= {
2.- L cos 4te
−2 t
}= {
3.- L t e
10 −7t
}=
{
4.- L e (t − 1)
2t 2
}= {
5.- L cos 5t (1 − e + 3e
t −4 t
)} =
Si L f (t ) = L
−1
{F ( s)}
=L
−1
{F ( s − a)}
=L
−1
{F ( s) s→s −a }=e at
f (t )
Ejemplos:
−1 s
Calcular L 2 =
s + 6 s + 11
Solución:
Completando el trinomio:
−1 s
L =
( s + 3) 2
+ 2
Completando la trasformada
−1 s +3−3
L =
( s + 3) + 2
2
s+3 −1 1
L
−1
−3 L =
( s + 3) + 2 ( s + 3) + 2
2 2
Resolviendo:
3 2 −3t
e −3t cos 2t − e sen 2t =
s −1 1 2
L
−1
2 −3 L 2 =
s + 2 s → s +3 s + 2 s → s +3 3 2
e −3t cos 2t − sen 2t
2
−1
1 1
1.- L 3
= 2.- L
−1
2 =
( s + 2 ) s − 6 s + 10
s −1
s
3.- L
−1
2 = 4.- L 2
=
s + 4s + 5 ( s + 1)
−1
2 s − 1 s / 2 + 5 / 3
5.- L 2 3
= 6.- L
−1
2 =
s ( s + 1) s + 4s + 6
−1 1 −1 1
7.- L 3
= 8.- L 4
=
( s + 2) ( s − 1)
−1 1 −1 1
9.- L 2 = 10.- L 2 =
s + 2s + 5 s − 2s + 3
Tarea 20
−1 s −1 1 1
1.- L 2 = 2.- L + 2 =
s + 6 s + 11 ( s − 1) s + 2 s − 8
3
−1 2s + 5 −1 1
3.- L 2
= 4.- L 2 =
( s − 3) s − 4s + 5
−1 3!
5.- L 4
=
( s − 2)
Definición 7.2
Nótese que definimos u (t − a ) solamente en el eje t no negativo ya que esto basta para estudiar la
trasformada de Laplace. En un sentido más amplio, u (t − a ) =
0 para t < a
Ejemplo:
Teorema 7.6
Si a > 0 , entonces
L { f (t − a )u (t − a )} =
e L− as
{ f (t )}
Demostración:
∞
L { f (t − a)u (t − a)=
} ∫ e− st f (t − a)u (t − a)dt
0
a ∞
∫e f (t − a )u (t − a )dt + ∫ e − st f (t − a )u (t − a )dt
− st
=
0 a
↑ ↑
cero para 0 ≤ t < a uno para t ≥ a
∫e
− st
= f (t − a )dt
a
L { f (t − a)u (t − a)} =
∫e
− s (u + a )
f (u )du
0
∞
= e − as ∫ e − su f (u )du
0
= e − as L { f (t )}
Ejemplo:
L {(t − 2) u (t − 2)} =
3
L
3
e L {t }
{(t − 2) u(t − 2)} = −2 s 3
Ejercicios: use el segundo teorema de traslación para calcular las siguientes trasformadas.
f (t − a )u (t − a ) =L
−1
{
e − as F ( s ) }
Ejemplo:
e− s
F (s) =
s ( s + 1)
−s 1 −1 − s 1
=
−1
L e − L e
s s + 1
f (t ) = u (t − 1) − u (t − 1)e − (t −1)
Ejercicios: use el segundo teorema de traslación de forma inversa para calcular las siguientes
trasformadas.
−1 e −π s /2 −1 e −2 s
1.- L 2 = 2.- L 3 =
s + 9 s
−1 e −π s −1 e− s
3.- L 2 = 4.- L =
s + 1 s ( s − 1)
Tarea 21
−1 4e − s −1 4e −3 s
1.- L 2 = 2.- L 2 =
s s
−1 1 − e −2 s 0, t <π /4
3.- L 2 = 4.- g (t ) =
s cos(t − π / 4), t ≥ π / 4
5.- f (t )= {u (t − 3) + u (t − 2)}=
Bibliografía Complementaria del Tema:
A veces se encuentra un factor integrante por inspección, después de agrupar convenientemente los términos
de una ecuación, al reconocer cierto grupo de términos como formando parte de una diferencial exacta, por
ejemplo:
1
ECUACIONES DIFERENCIALES. Frank Ayres. Editorial Mc Graw Hill. Pag. 25
g(x) Forma de yp
1(una constante) A
5x + 7 Ax + B
3x 2 − 2 Ax 2 + Bx + C
x3 − x +1 Ax 3 + Bx 2 + Cx + E
e5 x Ae 5 x
(9 x − 2)e 5 x ( Ax + B)e 5 x
x 2e5 x ( Ax 2 + Bx + C )e 5 x
2
ECUACIONES DIFERENCIALES CON PROBLEMAS DE VALORES EN LA FRONTERA. Dennis Zill.
Editorial International Thomson Editores. Pag. 147
f (t ) L { f (t )} = F ( s)
1 1
s
t 1
s2
tn n!
s n +1
senkt k
s + k2
2
cos kt s
s + k2
2
sen 2 kt 2k 2
s ( s 2 + 4k 2 )
cos 2 kt s 2 + 2k 2
s ( s 2 + 4k 2 )
e at 1
s−a
senhkt k
s − k2
2
cosh kt s
s − k2
2
3
ECUACIONES DIFERENCIALES CON PROBLEMAS DE VALORES EN LA FRONTERA. Dennis Zill.
Editorial International Thomson Editores. Pag. AP-24
t n e at n!
( s − a ) n +1
e at senkt k
( s − a)2 + k 2
e at cos kt s−a
( s − a)2 + k 2
tsenkt 2ks
(s + k 2 )2
2
t cos kt s2 − k 2
(s 2 + k 2 )2
I. Identifique las variables dependientes e independientes donde sea posible, luego clasifique las
siguientes ecuaciones diferenciales de acuerdo a su tipo, orden, grado y linealidad.
Var. Var.
Ecuación Diferencial Tipo Orden Grado Linealidad
Dep. Indep.
d2y dy
1. 2
+2 +x=0
dx dx
∂ 2u ∂ 2u
2. a2 =
∂x 2 ∂t 2
dy 2
4. y 1 + = C
dx
∂3x ∂3x
5. + 2 y = y + 2c
∂c 3 ∂y 3
II. Verifique si las siguientes funciones son solución de la correspondiente ecuación diferencial.
1
1. y− =0 de y ' − y tan x = 0
cos x
3
2. y=− de y' = 3 y 2
3x + 2
3. y = 1 + c 1 − x 2 de (1 − x )y' + xy = x
2
dy (x cos x + sen x − 1) y
4. xy 3 − xy 3 sen x = 1 ; =
dx 3(x − xsenx )
cos 2 x
1. y ' = , y (π ) = −1
y2
π
(
2. 3e x tan y dx + 2 − e x sec 2 y dy = 0, ) y (0 ) =
4
dy y 2 − 1
3. = , y (2 ) = 2
dx x 2 − 1
y (− 1) = −1
dy
4. x 2 = y − xy,
dx
y y y
5. x cos y′ = π /2
y cos − xsen ; y (1) =
x x x
y y y
6. xy cos + x 2 sen y′= y 2 cos ; y (1)= π / 2
x x x
(
7. 3 xy + x
2 3
) y ' =3 y 3
+ x 2 y ; y (1) =2
8. x e ( y/ x
− 1) y ' = e y / x ( y − x ) ; y (1)= 0
( )
9. x 2 + y 2 dx + 2 xydy = 0
( ) (
10. x x 2 + y 2 − y dx + y x 2 + y 2 − x dy )
11. (tan x − sen x sen y )dx + cos x cos ydy = 0
( ) (
12. 4 x 3 y 3 − 2 xy dx + 3 x 4 y 2 − x 2 dy = 0 )
(
13. 6 xy dx + 4 y + 9 x 2 dy = 0 )
2
14. cos x dx + 1 + sen x dy = 0
y
(
15. 10 − 6 y + e −3 x dx − 2dy = 0 )
Ecuaciones Diferenciales Ordinarias Página 92
( ) ( )
16. y 2 + xy 3 dx + 5 y 2 − xy + y 3sen y dy = 0
2.- Según la Ley de Newton de enfriamiento, la velocidad a la que se enfría una sustancia al aire
libre es proporcional a la diferencia entre la temperatura de la sustancia y la del aire. Si la
temperatura del aire es 30ºC y la sustancia se enfría de 100ºC a 70ºC en 15 minutos. ¿Cuándo será
40ºC la temperatura de la sustancia?
3.- Un termómetro que está en el interior de una habitación se lleva al exterior, en donde la
temperatura del aire es de 5ºF. Después de 1 minuto el termómetro marca 55ºF y después de 5
minutos marca 30ºF. ¿Cuál es la temperatura inicial de la habitación?
4.- Si una barra pequeña metálica, cuya temperatura inicial es de 20°C, se deja caer en un recipiente
de agua hirviendo, ¿Cuánto tiempo tardará en alcanzar los 90°C si se sabe que la temperatura
aumento 2°C en un segundo?, ¿Cuánto tiempo tardará en alcanzar los 98°C?
6.- Una fría mañana de noviembre (aproximadamente 2ºC de temperatura), en un oscuro y frío
callejón, un vagabundo encuentra el cadáver de un hombre asesinado. La policía y el forense llegan
a la escena del crimen a las 7:00 hr. de la mañana, comprobando que la temperatura del cadáver es
de 31ºC. Una hora más tarde su temperatura ha descendido a 27ºC. Un testigo afirma que sólo vio
salir del callejón a dos personas: el propio vagabundo, aproximadamente una hora antes de que
llegara la policía; y un vecino del inmueble contiguo, alrededor de las cinco de la mañana.
Suponiendo que la temperatura media de una persona es de 36ºC, ¿quién es el asesino?
1.- y′′ − 9 y =
=
0 ; y1 e=
3x
, y2 e −3 x
1.- ( x + 1) y′ − ny= e x ( x + 1)
n +1
dz
2.- ( 6 − 2 wz ) + z2 =
0
dw
( ) ( )
2
3.- x x 2 + 1 y′ + =
y x x2 + 1
1.- y′′ − 10 y′ + 25 y =
0
2.- y′′ + y′ + y =0
3.- 2 y′′ − 5 y′ − 3 y =
0
5.- y′′ + 16 y =
0 ; y (0) = 2, y′(0) = −2
1.- y′′ =
senx + sec 2 x 1, y′(0) =
y (0) = 1
3.- y′′ =
cos x − csc 2 x 1, y′(π / 2) =
y (π / 2) = 1
5.- y′′ + 25 y =
20 sen5 x
1.- y′′ + 3 y′ + 2 y =
sen(e x )
2.- y′′ + 2 y′ + y =e − x ln x
3.- y′′′ + y′ =
sec x
5.- y′′ + 3 y′ + 2 y = 1/ (1 + e x )
1.- Un peso de 8 lb se coloca en el extremo inferior de un resorte helicoidal que está suspendido del
techo. El peso se encuentra en reposo en su posición de equilibrio, por esta razón está estirado 6
pulg. Después, el peso se empuja 3 pulg hacia debajo de su posición de equilibrio y se suelta en t=0
con una velocidad inicial de 1 pie/seg, dirigida hacia abajo. Despreciando la resistencia del medio y
suponiendo que ninguna fuerza externa está presente, determine la amplitud, periodo, frecuencia y
la Ecuación de movimiento.
2.- Un peso de 32 libras está unido al extremo inferior de un resorte que está suspendido del techo.
El peso se encuentra en reposo en su posición de equilibrio y el resorte está estirado 2 pies.
Después, el peso se desplaza 6 pulg hacia debajo de su posición de equilibrio y se suelta en t=0.
Ninguna fuerza externa actúa; sin embargo la resistencia del medio en libras es numéricamente
igual a 4(dx/dt), donde dx/dt es la velocidad instantánea en pies por segundo. Determine el
movimiento resultante del peso en el resorte.
3.- Un peso de 16 lb está unido al extremo inferior de un resorte que está suspendido del techo, la
constante del resorte es 10 lb/pie. El peso se encuentra en reposo en su posición de equilibrio. Al
empezar en t=0 se aplica una fuerza externa dada por f(t)= 5cos2t. Determine el movimiento
resultante si la fuerza de amortiguamiento en libras es numéricamente igual a 2(dx/dt), donde dx/dt
es la velocidad instantánea en pie/seg.
4.- Un cuerpo que pesa 2 lb estira un resorte 6pulg. Dicho cuerpo se suelta en t = 0 desde un punto
que está 8 pulg bajo la posición de equilibrio, con una velocidad dirigida hacia arriba de
4 / 3 pie / seg . Determinar la función x(t ) que describe el movimiento libre resultante.
5.- Un cuerpo que pesa 8 lb estira un resorte 2pies. Suponiendo que una fuerza de amortiguación
numérica igual a dos veces la velocidad instantánea actúa sobre el sistema y que el peso se suelta
desde la posición de equilibrio con una velocidad dirigida hacia arriba de 3 pie / seg . Determinar
la función x(t ) que describe el movimiento resultante.
{ }=
1.- ζ t
n
2.- ζ {senkt} = 3.- ζ {cos kt} =
−1
4.- ζ
−1
{ }
e at = 5.- ζ {t cos kt} =
1.- y′′ − 8 y=
′ − 9 y 10 =
y (0) 0, y′(0) 4
=
2.- y′′ − 6 y=
′ + 8 y 2e 3 t =
y (0) 0, y′(0) 2
=
7.- y′′=
− 4y 2 =
y (0) 0, y′(0) 0
=
8.- y′′ − ( 5 / 2 )=
y′ + y 0 = y′(0) 1/ 2
y (0) 1, =
9.- y′′ − 2 y=
′ − 3y 0 =
y (0) 3, y′(0) 2
=
′′ + y 2 cos t
10.- y= = y′(0) 0
y (0) 2, =
Calcule I C
Calcule I L , i ( L)(0) = 10 A
1. Un circuito RCL tiene R=10 ohmios, C= 10-2 faradios, L=1/2 henrios, y un voltaje aplicado
E=12 voltios. Asumiendo que no hay corriente inicial y no hay carga inicial en t=0 cuando el
voltaje se aplica por primera vez, halle la corriente resultante en el sistema.
2. Un circuito RCL tiene R=180 ohmios, C= 1/280 faradios, L=20 henrios, y un voltaje aplicado
E=10sent voltios. Suponiendo que no hay carga inicial en el condensador sino una corriente inicial
de 1 amperio para t=0 cuando se aplica por primera vez el voltaje, halle la carga resultante en el
condensador.
3. Un circuito RCL tiene R=5 ohmios, C= C= 10-2 faradios, L=1/8 henrios, y un voltaje aplicado
E=sent voltios., halle la corriente en condiciones estables en el circuito. Sugerencia (las
condiciones iniciales no se necesitan).