Facultad de Ingenieria Universidad Nacio
Facultad de Ingenieria Universidad Nacio
CUADERNO DE EJERCICIOS DE
ECUACIONES DIFERENCIALES
Se invita a los estudiantes y profesores a que comuniquen a los autores las observaciones y
sugerencias que mejoren el contenido de la obra, con el fin de que se incorporen en una
futura edición definitiva.
PRÓLOGO
Conscientes de la importancia de que los alumnos cuenten con el material didáctico que
complemente su clase de teoría de la asignatura, y que los profesores que la imparten
tengan un apoyo para sus clases, se ha trabajado con entusiasmo y compromiso para que
esta obra en su versión actualizada contribuya al propósito con el que originalmente fue
creada.
Cabe señalar que aunque este trabajo fue revisado minuciosamente, existe la posibilidad de
que se hayan cometido errores involuntarios, por lo que se agradece de antemano a los
profesores y alumnos los comentarios que hagan sobre esta obra, lo que contribuirá a
mejorar la calidad de futuras ediciones.
ECUACIONES DIFERENCIALES
= k ( 4 − x) ( 1 − x) , donde k es constante
dx
a)
dt
RESOLUCIÓN
1
2) Determine los valores que pueden tomar las constantes k , m , n , t de modo que
⎛d y⎞ ⎛ d ky ⎞
la ecuación diferencial x ⎜ ⎟ + 4⎜ k ⎟
= x n sea lineal, de segundo orden
m
⎝dx⎠ ⎝d x ⎠
t
y de coeficientes constantes.
RESOLUCIÓN
Si m = 1 es lineal
Si k = 2 es de segundo orden
Si t = 0 x t = x 0 = 1 es constante y la ecuación es de coeficientes constantes
2
1.2 SOLUCIÓN DE LA ECUACIÓN DIFERENCIAL.
1) Sea la función y ( x ) = C 1 e − x + C 2 e 2 x
− − 2y = 0
d 2y dy
2
dx dx
RESOLUCIÓN
= − C 1 e − x + 2C 2 e 2 x
dy
dx
Al derivar nuevamente se obtiene
= C 1 e − x + 4C 2 e 2 x
d 2y
2
dx
Enseguida se sustituyen estas funciones en la ecuación diferencial:
(
C 1 e − x + 4 C 2 e 2 x − − C 1 e − x + 2C 2 e 2 x ) − 2( C 1e
−x
)
+ C 2e 2x = 0
esto es
3
2 = C1 e( ) + C 2 e( )
0 0
2 = C1 + C2 .............. ( 1)
la función solución
1 = − C 1 e( ) + 2C 2 e( )
0 0
1 = − C 1 + 2C 2 .............. ( 2 )
y y ''' = −
1
2
RESOLUCIÓN
La función dada debe satisfacer la ecuación diferencial, por lo que es necesario calcular
y ' , y '' , y '''
Para y ' se tiene
3⎛ 2 ⎞
y' = ⎜ ⎟ ( x) 2 = 3 ( x) 2
1 1
2⎝ 3⎠
4
Para y '' se tiene
( x) 2
−
y '' =
1
3
2
Finalmente se tiene y '''
( x) 2 = − ( x) 2
− −
y ''' = −
3 3
1 3 3
2 2 4
⎛ 2 ⎞ ⎛ − ⎞
⎜⎜ ( x ) 2 ⎟⎟ ⎜⎜ − ( x ) 2 ⎟⎟ = −
3 3
3 1
⎝ 3 ⎠ ⎝ 4 ⎠ 2
− = −
1 1
2 2
y = (x −C) 2 + C
RESOLUCIÓN
5
b) Para determinar si existe alguna solución singular, se deriva la ecuación diferencial
parcialmente respecto a y ' :
∂
⎡( y ' ) 2 + 2 y ' + 4 x − 4 y ⎤ = 0
∂ y' ⎣ ⎦
de aquí se obtiene
2 y' + 2 = 0
y' = −1
Este valor de y ' es sustituido en la ecuación diferencial
( − 1) 2 + 2 ( − 1) + 4 x − 4 y = 0
6
en a) indicando, en su caso, qué tipo de solución es:
I) y = x2
II) y = −1
III) y = x 2 + 2 x
RESOLUCIÓN
y' = 2 (x − C)
lo cual se logra derivando dicha solución
es decir,
(x − C) =
y'
2
este resultado se sustituye en la solución general, con lo que se tiene
⎛ y'⎞
y = ⎜ ⎟ − 1
2
⎝2⎠
que es la ecuación diferencial y es posible reescribir como
4 y = ( y ') − 4
2
( y ') = 4y + 4
2
solución general.
7
Para y = x 2 + 2 x , se deriva
y' = 2 x + 2
( )
al sustituir en la ecuación diferencial
(2 x + 2) 2 = 4 x 2 + 2 x + 4 ,
4 x 2 + 8 x + 4 = 4x 2 + 8 x + 4
por lo tanto sí es solución y es particular pues se obtiene de la solución general cuando
C = − 1.
RESOLUCIÓN
función y , de donde
m = = 4 − 2x
dy
dx
que representa la ecuación diferencial de la familia.
b) Para resolver este inciso se requiere la solución general de la ecuación diferencial, es decir,
la función y . Entonces al integrar
d y = (4 − 2 x) d x
∫ dy = ∫ (4 − 2 x) d x
de donde resulta
y = 4 x − x 2 + C , que es la solución general.
8
Para la condición x = 0 , se tiene y = 0 y al sustituir en la solución general resulta
0 = C
Finalmente, sustituyendo este valor en la solución general, se obtiene
y = 4 x − x 2 que representa el elemento de la familia que pasa por (0 , 0)
6) Obtenga la ecuación diferencial cuya solución general es la familia de rectas tales que su
pendiente y su abscisa al origen son iguales.
RESOLUCIÓN
y − y 0 = m(x − x 0 )
La ecuación de la recta a considerar corresponde a la de punto y pendiente,
Finalmente
y = ax − a2
que es la ecuación de la familia de rectas.
( y ') 2 − y' x + y = 0
9
7) Obtenga la ecuación diferencial cuya solución es la familia de parábolas con vértice sobre
el eje x , eje focal paralelo al eje y , con distancia del foco al vértice igual a a .
RESOLUCIÓN
4 a y ' = 2( x − h )
donde h es la única constante arbitraria; al derivar una vez se tiene
despejando x − h
x − h = 2ay ' ……………………...…………………(2)
sustituyendo (2) en (1)
4a y = ( 2a y ' ) 2
Finalmente
4 a y = 4 a 2( y ' ) 2 es la ecuación diferencial.
RESOLUCIÓN
(x − h)2 + (y − k )2 = r2
De acuerdo a la descripción del enunciado, el centro de las circunferencias tiene ordenada cero,
es decir, k = 0 ; en lo que respecta al radio, es también variable. Entonces, la ecuación que
describe a la familia es
(x − C1 )2 + y2 = C2 2
10
donde C 1 y C 2 son constantes arbitrarias.
La ecuación diferencial deseada es de segundo orden (el orden corresponde al número de
constantes arbitrarias), por lo que es necesario derivar dos veces y eliminar las constantes
C1 y C2
2( x − C 1 )
Para la primera derivada
+ 2 y y' = 0
Al simplificar
x − C 1 + y y' = 0
C 1 = x + y y'
Nuevamente se deriva
0 = 1 + y y '' + y ' y '
Finalmente
( y ') 2 + y y '' + 1 = 0
es la ecuación diferencial.
RESOLUCIÓN
11
C1 =
(y − γ sen x )
eαx
..........................................................................(3)
y ' − γ cos x
C1 =
αe α x
............................................................................(4)
y al simplificar
α y − α γ sen x = y ' − γ cos x
Finalmente, al reordenar los términos se obtiene
y ' − α y = γ cos x − α γ sen x
que es la ecuación diferencial.
= P( 1 − P )
dP
dt
RESOLUCIÓN
Para que la función dada sea solución de la ecuación diferencial debe satisfacerla, entonces es
necesario calcular la derivada de la función.
Se observa que la función es implícita por lo que es conveniente expresarla de la siguiente
manera
P = C 1et (1 − P )
P' = − C 1 e t P' + C 1 e t − C 1 e t P
P ' = C1 e t ( 1 − P − P ' )
12
pero sabemos que C 1 e t =
P
1 − P
, por lo que
P' = (1 − P − P ')
P
1 − P
P ' ( 1 − P ) = P ( 1 − P − P ')
P ' = P (1 − P )
RESOLUCIÓN
Para verificar es necesario que la función satisfaga a la ecuación diferencial; obtenemos y'
y y ''
y ' = C1 e x − C2 e − x − 4
y '' = C1 e x + C2 e − x
(C e ) − (C e )
enseguida se sustituyen en la ecuación diferencial
1
x
+ C2 e − x 1
x
+ C2 e − x − 4 x = 4 x
13
inicial y ' ( 0 ) = 0
Para obtener la solución particular se requiere derivar la solución general y aplicar la condición
Al derivar se obtiene
y 'G = C1 e x − C2 e − x − 4
2 = C1 + C2
0 = C1 − C2 − 4
RESOLUCIÓN
Inicialmente se calcula la derivada de la función que representa a la solución general; dado que
no está en forma explícita, se deriva implícitamente
2y y ' + 2C = 0
o bien
yy' + C = 0
de donde
C = − yy'
Enseguida se sustituye en la solución general
y 2 + 2 ( − y y ') x = ( − y y ') 2
efectuando operaciones
14
y 2 − 2 y y ' x = y 2 ( y ') 2
dividiendo entre y
y ( y ') 2 + 2 x y ' − y = 0
Q ( y) ⎜ ⎟ + P ( x)
⎛ dy ⎞
− y = 0
2
dy
⎝ dx ⎠ dx
Q ( y ) = y , P ( x ) = 2x
Por lo que al comparar los coeficientes se obtiene
15
1.3 PROBLEMA DEL VALOR INICIAL.
4 x ( y ') 2 + 2 x y ' − y = 0
inicial y ( 1) = 2
y obtenga una solución particular que satisfaga la condición
RESOLUCIÓN
La solución general contiene una sola constante arbitraria, por lo que es necesario derivar una
2 ( y − C)y' = C
vez la solución general, esto es
⎛ ⎞
y ' = ⎜⎜ ⎟⎟
⎝ 2 ( y − C) ⎠
C
⎛ ⎞ ⎛ ⎞
4 x ( y ') + 2 x y ' − y = 4 x ⎜⎜ ⎟⎟ + 2 x⎜⎜ ⎟⎟ − y
2
⎝ 2 ( y − C) ⎠ ⎝ 2 ( y − C) ⎠
2 C C
x C 2 + xC ( y − C ) − y ( y − C )
=
2
(y − C)
2
+ x y C + x C 2 − y( y − C)
=
(y − C) 2
2 2
xC
x yC − y( y − C)2
=
(y − C)2
x y C − y (C x )
ecuación diferencial anterior se escribe
= = 0
Cx
por lo tanto se verifica que sí es solución.
16
Para obtener la solución particular se aplica la condición inicial y ( 1) = 2 en la solución
(2 − C ) 2 = C ( 1)
C2 − 5C + 4 = 0
Esta última expresión es una ecuación de segundo grado que al resolver determina
C = 1 y C = 4
Para cada valor obtenido se tiene una solución particular, es decir
Para C = 1 (y − 1) = x
2
Para C = 4 (y − 4) 2 = 4 x
17
EJERCICIOS PROPUESTOS
x (2 x + C ) − ( x x + C x) = x x
Solución
x x = x x
Se verifica.
a) x y '' − 4 y ' − 5 y = e 3 x
∂2y ∂y ∂y
= + ky −
∂t ∂t ∂s
b) 2
donde k es una constante
c) y '' + 2 x 3 y ' − ( x − 1) y = x y 2
3
Solución
a) x y '' − 4 y ' − 5 y = e 3 x ; ordinaria, coeficientes constantes, lineal, orden 2, y
variable dependiente, x variable independiente.
∂2 y ∂y ∂y
= + ky − k = cte ; parcial, coeficientes constantes, lineal, orden
∂t ∂t ∂s
b) 2
18
c) y '' + 2 x y ' − ( x − 1) y = x y
3
3 2 ; ordinaria, coeficientes variables, no lineal,
( x + y y ') 2 =
Solución
⎛ x + y y'⎞
x 2 + y 2 − 2⎜ ⎟ ( x + y) +
⎝ x + y' ⎠ ( x + y ') 2
0
⎛x2 ⎞
uniparamétrica de soluciones y = ⎜ + C⎟
2
⎝ 4 ⎠
Determine
a) Una solución particular.
b) Una solución singular.
Solución
⎛x2⎞
a) Si C = 0 , una solución particular es y = ⎜ ⎟
2
⎝ 4 ⎠
b) Una solución singular es y = 0
Verifique que y =
2
1 − Ce x
5) , donde C es una constante arbitraria, es una familia
y ( y − 2)
=
dy
uniparamétrica de soluciones de la ecuación diferencial
dx 2
19
6) Determine los valores que pueden tomar las constantes k , m , n , t de modo que la
⎛ d ny ⎞ ⎛ d ky ⎞
+ 4⎜ = x sea de segundo orden y de
3
ecuación diferencial x ⎜
⎜ d x n ⎟⎟ ⎜ d x k ⎟⎟
m
t
⎝ ⎠ ⎝ ⎠
segundo grado.
Solución.
⎛ d ny ⎞ ⎛ d ky ⎞
+ 4⎜ = x
3
Para que la ecuación diferencial x ⎜
⎜ d x n ⎟⎟ ⎜ d x k ⎟⎟
m
t
⎝ ⎠ ⎝ ⎠
sea de segundo orden
y1 = − y2 = x + 1 , y3 = x + 3
1 2
x ,
4
son solución de la ecuación obtenida en el inciso anterior e indique, en cada caso, qué tipo
de solución es.
Solución
a) y = x y ' + ( y ') 2
y 3 no es solución
y G = C 1 + C 2 e t − 20 t
y ' ( 0 ) = 40
20
Solución
y p = − 60 + 60 e t − 20 t
x2 + (y − C) 2 = 4
(4 ) ( y ')
9) Si es la solución general de la ecuación diferencial
x2 = − x2 2
, obtenga, si existe, una solución singular de la ecuación
diferencial dada.
Solución
Existen dos soluciones singulares que son x = 2 , x = −2
− 4x = 0
dy
y
dx
21
TEMA 2
RESOLUCIÓN
y 2 dx − x yd y = x 2 y d y
( )
Agrupando términos para cada diferencial se tiene
y2 dx = x2 y + x y d y
=
dx ydy
x + x
2
y2
El paso inmediato es integrar ambos miembros
∫ =
∫
dx ydy
x + x
2
y2
∫ ∫
donde se obtiene
⎛1 1 ⎞
⎜ − ⎟ dx =
dy
⎝x x + 1⎠ y
integrando
Ln x − Ln x + 1 = Ln y + C
22
por propiedades de logaritmos se obtiene
⎛ x ⎞
Ln ⎜ ⎟ = Ln ( C y )
⎝ x + 1⎠
y aplicando la función exponencial
= Cy
x
x + 1
finalmente
y =
C ( x + 1)
x
, que es la solución general
2) (
Resuelva la ecuación diferencial y ' 1 + e x )=e x−y
utilizando el método de
separación de variables.
RESOLUCIÓN
( )
La ecuación diferencial también se puede escribir de la siguiente manera
1 + ex =
dy ex
dx ey
al separar variables se obtiene
⎛ ex ⎞
ey dy = ⎜ x ⎟
⎝1 + e ⎠
dx
∫ ∫
posteriormente se integran ambos miembros
ey dy =
e x dx
1 + ex
( )+C
de donde se obtiene
e y = Ln 1 + e x
23
(3 x − 6 y + 4) + x − 2y + 3 = 0
dy
dx
Sugerencia: utilizar el cambio de variable v = x − 2 y
RESOLUCIÓN
(3 x − 6 y + 4) + x − 2y + 3 = 0
dy
dx
al utilizar la sustitución v = x − 2 y se despeja a la variable y , asimismo se calcula
la derivada con respecto a x
x − v 1⎛ dv ⎞
y = = ⎜1 − ⎟
dy
2⎝ dx ⎠
;
2 dx
⎛ ⎛x − v⎞ ⎞⎛ 1 1 dv ⎞ ⎛ x − v⎞
⎜ 3 x − 6 ⎜ 2 ⎟ + 4 ⎟ ⎜ 2 − 2 dx ⎟ + x − 2 ⎜ 2 ⎟ + 3 = 0
⎝ ⎝ ⎠ ⎠⎝ ⎠ ⎝ ⎠
se efectúan las operaciones indicadas y se reducen términos
⎛1 1 dv ⎞
(3 x − 3 x + 3 v + 4) ⎜ − ⎟ + x − x + v + 3 = 0
⎝2 2 dx ⎠
⎛ dv ⎞
(3 v + 4) ⎜ 1 − ⎟ + 2v + 6 = 0
⎝ dx ⎠
3v + 4 − 3v − 4 + 2v + 6 = 0
dv dv
dx dx
(3 v + 4 ) dv = ( 5 v + 10 ) dx
se llega a una ecuación diferencial de variables separables
1 ( 3 v + 4)
dx =
5 ( v + 2)
dv
∫ ∫
integrando
3v + 4
dx =
1
v + 2
dv
5
para el segundo miembro
24
3v + 4
= 3 −
2
v + 2 v + 2
∫ ∫
entonces
dx = 3 −
1 2
v + 2
dv
5
x =
1
5
( 3 v − 2 Ln ( v + 2) ) + C
Sin embargo, de la sustitución v = x − 2 y se obtiene
x =
1
5
( 3 x − 6 y − 2 Ln ( x − 2 y + 2) ) + C
que es la solución general.
RESOLUCIÓN
Se tiene
= tan 2 ( x + y )
dy
dx
La forma del argumento de la función sugiere emplear la sustitución
u = x + y ...........................................(1)
lo que implica calcular du , entonces
du = dx + d y .........................................(2)
dy
ahora bien, en la ecuación diferencial es necesario el término ; lo obtenemos de (2)
dx
= − 1 .........................................(3)
dy du
dx dx
Enseguida se sustituyen (1) y (3) en la ecuación diferencial
25
− 1 = tan 2 ( u )
du
dx
Se tiene una ecuación de variables separables, por lo que
= tan 2 ( u ) + 1
du
dx
= dx
tan ( u ) + 1
du
2
= dx
sec 2( u )
du
cos 2( u ) du = dx
∫ cos ( u ) du = ∫ dx
Entonces ya es posible integrar ambos miembros
2
∫ ∫ dx
⎛1 ⎞
⎜ + cos 2 u ⎟ du =
1
⎝2 2 ⎠
de donde resulta
u + sen 2 u = x + C
1 1
2 4
pero de la sustitución efectuada se tiene
u = x + y
entonces del resultado anterior
( x + y ) + sen 2 ( x + y ) = x + C
1 1
2 2
2 ( x + y ) + 2 sen 2 ( x + y ) = 4 x + 4 C
o bien
2 ( y − x ) + 2 sen 2 ( x + y ) = C 1
agrupando términos
; C1 = 4C
26
⎛1 x⎞
(
⎜ − ⎟ dx + x cos y d y = 0
⎝y y⎠
2
)
sujeta a la condición inicial y ( 1) = π
RESOLUCIÓN
∫ ∫ y cos y d y
Luego de separar variables se integran ambos miembros de la igualdad
x − 1
dx =
x2
La integral del lado izquierdo se resuelve separando en dos nuevas integrales; la del lado
derecho se resuelve integrando por partes
∫ dx −
∫ =
∫ y cos y d y
x dx
∫ ∫x ∫ y cos y d y
x2 x2
− −2
dx =
dx
x
x −1
Ln x − + c 1 = y sen y + cos y + c 2
−1
Reordenando y agrupando constantes se obtiene
Ln x + = y sen y + cos y + C
1
x
27
que es la solución general de la ecuación diferencial. Sin embargo, se tienen condiciones
iniciales, las cuales se aplican en la solución general para obtener el valor de la constante C
y ( 1) = π ⇒ Ln ( 1) + = π sen π + cos π + C
1
(0 ) + 1 = ( 0 ) + ( − 1) + C
1
2 = C
Entonces, la solución particular que satisface la condición inicial dada es
y sen y + cos y + 2 = Ln x +
1
x
28
2.2 ECUACIONES DIFERENCIALES DE COEFICIENTES
HOMOGÉNEOS.
RESOLUCIÓN
(x )
Se tiene la ecuación
2
+ y 2 dx − x y d y = 0
Se analiza qué tipo de ecuación es; se descarta que sea de separación de variables y, por
inspección, se puede afirmar que es de coeficientes homogéneos de grado 2.
Se realiza la sustitución y = u x teniendo presente que es necesario calcular d y . Al
calcularlo se obtiene d y = u dx + x du , que al sustituirlo en la ecuación diferencial
(x )
nos lleva a
2
+ ( u x ) 2 dx − x ( u x ) ( u dx + x du )
(x )
Efectuando operaciones y simplificando
2
+ u 2 x 2 dx − x 2 u ( u dx + x du ) = 0
x 2 dx + u 2 x 2 dx − u 2 x 2 dx − u x 3 du = 0
se obtiene una ecuación diferencial de variables separables
x 2 dx − x 3 u du = 0
∫ ∫ u du = 0
Separando variables
−
dx
x
e integrando
Ln x − = C
u2
2
Finalmente se obtiene
29
Ln x − = C
y2
, que es la solución general de la ecuación diferencial.
2x 2
(x + )
2) Resuelva la ecuación diferencial de coeficientes homogéneos
y dx = y 2 − x 2 dy
RESOLUCIÓN
( ) (u dx + x du )
Enseguida se sustituye en la ecuación diferencial
u x dx = x + u2 x2 − x2
( )
donde luego de efectuar operaciones y reordenar se tiene
− u u 2 − 1 x dx = x 2 1 + u 2 − 1 du
1 + u2 − 1
− =
dx
u u2 − 1
du
x
⎛ u2 − 1 ⎞
− = ⎜ + ⎟ du
⎜u u 2 − 1 ⎟
dx 1
x
⎝ u u 2
− 1 ⎠
al simplificar
⎛ 1⎞
− = ⎜ + ⎟ du
⎜u u 2 − 1 u ⎟
dx 1
⎝ ⎠
x
Posteriormente se integra
30
−
∫ =
∫u du +
∫
dx 1 du
x u 2
− 1 u
− Ln x = ang sec u + Ln u + C
− Ln x = ang sec + Ln + C
y y
x x
que es la solución general de la ecuación diferencial.
(x + )
grado de homogeneidad.
x y + y2
F (x , y ) = − xy − +
2 x5
y
x y4
RESOLUCIÓN
F (λ x , λ y ) = λ n F ( x , y )
(λ x) (λ y ) + (λ y ) 2 (λ x) 5
entonces
F (λ x , λ y ) = λ x + (λ y ) 2 − (λ x) (λ y ) − +
(λ x) (λ y) 4
F (λ x , λ y ) = λ x + λ − xy −
(
λ2 x y + y2 ) +
λ5 x5
λx λ5 y4
2
y
⎡
(
F (λ x , λ y ) = λ ⎢ x +
⎣
y2 − x y − ) x y + y2
x
x5⎤
+ 4⎥
y ⎦
F (λ x , λ y ) = λ F ( x , y )
31
2.3 ECUACIONES DIFERENCIALES EXACTAS. FACTOR
INTEGRANTE.
+ P ( x ) y = x es µ ( x ) = e x ,
dy 2
1) Si un factor integrante de la ecuación diferencial
dx
a) P ( x )
obtenga:
b) la solución general y = y ( x )
RESOLUCIÓN
+ P ( x) y = x
dy
dx
∫ P ( x ) dx
Esta ecuación diferencial es lineal y de primer orden; un factor integrante es
µ ( x) = e ∫ y como µ ( x ) = e x
P ( x) d x
= x2
2
se tiene que por lo tanto
P ( x) = 2 x
+ 2x y = x
dy
dx
d y = ( x − 2 x y ) dx
que en forma diferencial se expresa como
(2 x y − x ) dx + d y = 0
2
multiplicando por su factor integrante e x se obtiene
e x ( 2 x y − x ) dx + e x dy = 0
2 2
( ) = ∂∂x ( e )
Esta ecuación diferencial, según se observa, es diferencial exacta pues cumple
∂
2x yex − xex
∂y
2 2
x2
32
= 2xex
2 2
2xex
∂ f ( x , y)
= M ( x , y ) = 2 x y e x − x e x …………….(1)
∂x
2 2
∂ f ( x , y)
así como
= N ( x , y ) = e x ……………………………..(2)
∂y
2
∫
Al integrar (2) resulta
f (x , y ) = e x d y + h ( x)
2
f ( x , y ) = y e x + h ( x ) ……………...…(3)
2
Por lo tanto
h '( x) = − x e x
2
h ( x) = − + C
2
ex
2
finalmente, la solución general es
− e = C
2 1 x2
yex
2
y en forma explícita
y = C e −x +
2 1
2
− x seny d y − ( 5 x − 5 x cos y ) dx = 0
2) Resuelva la ecuación diferencial
33
RESOLUCIÓN
x du − ( 5 x − 5 x u ) dx = 0
sustituyendo en la ecuación diferencial
v ( x) = e ∫
P ( x) d x
v ( x) = e ∫ = e5x
5 dx
se transforma en exacta
e 5 x du + 5 e 5 x u dx = 5 e 5 x d x
e 5 x du + 5 e 5 x u dx − 5 e 5 x d x = 0 ……………(3)
(5 e 5x
)
− 5 e 5 x u dx − e 5 x du = 0
M dx + N du = 0
∂M ∂N
= ∴ se comprueba que es exacta.
∂u ∂x
∫ 5e ∫ 5e
Integrando (3) con respecto a x se tiene
5x
u dx = 5x
dx
e5x u = e5x + C
u = 1 + C e −5x
pero u = cos y ; entonces
cos y = 1 + C e − 5 x
34
es la solución general.
⎛ 1⎞
Otra forma de resolver la ecuación (1) es multiplicarla por ⎜ − ⎟ , de donde
⎝ x⎠
sen y d y + 5 ( 1 − cos y ) dx = 0
d y + 5 dx = 0
sen y
1 − cos y
Ln ( 1 − cos y ) + 5 x = C 1
integrando resulta
1 − cos y = C 2 e − 5 x
que es la solución general en forma implícita.
(2 x y )
− sec 2 x dx + (x 2
)
+ 2 y dy = 0
RESOLUCIÓN
∂M
= 2x
∂y
N (x , y ) = x 2 + 2 y
∂N
= 2x
∂x
se cumple
35
∂M ∂N
=
∂y ∂x
f (x , y )
entonces es una ecuación diferencial exacta.
Al ser exacta se sabe que existe una función tal que
∂f ∂f
= dx +
∂x ∂y
df d y …………….………….…….(1)
d f = M ( x , y ) dx + N ( x , y ) d y ………..………(2)
donde
∫ N (x , y )dy + g (x)
De (4)
f (x , y ) =
∫
f (x , y ) = (x + 2 y )dy + g (x)
2
Luego de integrar
f (x , y ) = x 2 y + y 2 + g(x)
M ( x , y ) = 2 x y + g '( x )
igualando los segundos miembros de (3) y (5)
al sustituir al coeficiente M
2 x y − sec 2 x = 2 x y + g ' ( x )
resulta
g ' ( x ) = − sec 2 x
36
∫ sec
y al integrar
g(x) = − 2
x dx
g ( x ) = − tan x
f ( x , y ) = x 2 y + y 2 − tan x
f (x , y ) = C
Teniendo presente que la solución general es de la forma
se puede escribir
x 2 y + y 2 − tan x = C
que es la solución general en forma implícita.
RESOLUCIÓN
M ( x , y ) dx + N ( x , y ) d y = 0
verá si es o no exacta, para lo cual la llevamos a la forma
entonces se tiene
∂M ∂N
∂y ∂x
Enseguida se determinan y
∂M
= cos y
∂y
37
∂N
= sec y + tan y
∂x
∂M ∂N
≠ ∴
∂y ∂x
no es exacta
Dado que no es exacta, se procede a obtener un factor integrante adecuado para transformar
∫e
la ecuación diferencial en exacta.
∂M ∂N
−
∂y ∂x cos y − ( sec y + tan y )
g ( x) = =
N − 2 y cos y + x ( sec y + tan y )
∫e
se observa que no depende sólo de x .
µ( y) =
h ( y) d y
Ahora se propone un factor integrante , para el cual
∂N ∂M
−
∂x ∂y sen y + tan y − cos y
h( y) = =
M 1 + sen y
1 + sen y 1 + sen y
− cos y
h( y) = = −
cos y cos y
1 + sen y 1 + sen y
cos y
1 + sen y
h( y) = −
1 cos y
cos y 1 + sen y
Enseguida se procede a integrar
38
∫ ∫
⎛ 1 cos y ⎞
h( y) d y = ⎜ − ⎟ dy
⎝ + ⎠
∫ ∫
cos y 1 sen y
= sec y d y −
cos y
1 + sen y
dy
⎛ sec y + tan y ⎞
= Ln ⎜ ⎟
⎝ 1 + sen y ⎠
⎛ 1 + sen y ⎞
⎜ cos y ⎟
= Ln ⎜ ⎟
⎜ 1 + sen y ⎟
⎜ ⎟
⎝ ⎠
⎛ 1 ⎞
= Ln ⎜ ⎟
⎝ cos y ⎠
e∫
h ( y) d y
El elemento que se requiere es , entonces
⎛ 1 ⎞
e∫
h ( y) d y
Ln ⎜ ⎟
= e ⎝ cos y ⎠ =
1
cos y
Así, el factor integrante es
µ( y) =
1
cos y
{
µ ( y ) ( 1 + sen y ) dx − ⎡⎣ 2 y cos y − x ( sec y + tan y ) ⎤⎦ d y = 0 }
1
cos y
{(1 + sen y ) dx − ⎡⎣ 2 y cos y − x ( sec y + tan y ) ⎤⎦ d y = 0 }
de aquí se obtiene
⎛ 1 + sen y ⎞ ⎡ ⎤
⎜ ⎟ dx − ⎢ 2 y − ( sec y + tan y )⎥ d y = 0 …………..( 2 )
x
⎝ cos y ⎠ ⎣ cos y ⎦
Se procede a verificar si es o no exacta
1 + sen y
M (x , y ) =
cos y
39
N (x , y ) = − 2 y + ( sec y + tan y )
x
cos y
+
∂N sec y + tan y
1 sen y
= =
∂x
cos y cos y
cos y cos y
1 + sen y
=
cos y
cos y
1 + sen y
=
cos 2 y
∂M ∂N
=
∂y ∂x
Así, se cumple por lo que se verifica que ( 2 ) es exacta.
f ( x , y ) tal que
Por otro lado, se tiene que si una ecuación diferencial es exacta, entonces existe una función
d f = M (x , y )d x + N (x , y )d y
donde
∂
M (x , y ) = f (x , y )
∂x
∂
N (x , y ) = f (x , y )
∂y
f (x , y ) :
∫ M ( x , y ) d x + h( y)
Se puede proceder a obtener
f (x , y ) =
∫
1 + sen y
f (x , y ) = dx + h ( y )
cos y
⎛ 1 + sen y ⎞
f (x , y ) = x ⎜ ⎟ + h ( y ) ……………( 3 )
⎝ cos y ⎠
Por otro lado
40
∂
f (x , y ) = N (x , y )
∂y
entonces
∂ ⎛ 1 + sen y ⎞
f (x , y ) = x ⎜ ⎟ + h '( y)
∂y ⎝ cos y ⎠
2
⎛ 1 sen y ⎞
⎜ cos y cos y ⎟
= − 2 y + x⎜ + ⎟
⎜ cos y cos y ⎟
⎜ ⎟
⎝ ⎠
⎛ 1 + sen y ⎞
= − 2 y + x⎜ ⎟
⎝ cos y ⎠
2
h '( y ) = − 2 y
∫ y dy
de aquí se procede a integrar en ambos miembros, de donde
h( y) = − 2
h( y) = − y 2
f (x , y ) = C
Se recordará que la solución está dada por
0 − ( 1) 2 = C
41
− 1 = C
por lo que la solución particular que satisface la condición es
⎛ 1 + sen y ⎞
x⎜ ⎟ − y = − 1
⎝ cos y ⎠
2
(x ) ( )
5) Demuestre que
+ 3 x 3 sen y dx + x 4 cos y d y = 0
−1
no es exacta, pero que multiplicando esta ecuación por el factor x se obtiene
una ecuación exacta. Posteriormente obtenga la solución de esta última.
RESOLUCIÓN
Para que una ecuación diferencial sea exacta, es condición necesaria y suficiente que
∂ ∂
M (x , y ) = N (x , y )
∂y ∂x
por lo que inicialmente se procede a identificar las funciones M y N ,
M ( x , y ) = x + 3 x 3 sen y
N ( x , y ) = x 4 cos y
factor integrante µ ( x )
42
x − 1⎡
⎣ (x ) (
+ 3 x 3 sen y dx + x 4 cos y dy ⎤ = 0
⎦ )
(1 + ) ( )
de donde se obtiene
N ( x , y ) = x 3cos y
∂M
= 3 x 2cos y
∂y
∂N
= 3 x 2cos y
∂x
se cumple que
∂M ∂N
=
∂y ∂x
por lo tanto es exacta
∂f ∂f
f (x , y ) = dx +
∂x ∂y
dy
donde
∂f ∂f
= M (x , y ) ; = N (x , y )
∂x ∂y
f (x , y ) :
∫ N ( x , y ) d y + g ( x)
Se procede a obtener
f (x , y ) =
f (x , y ) =
∫ ( x cos y ) d y + g ( x )
3
f ( x , y ) = x 3 sen y + g ( x )
∂f
= M (x , y )
∂x
43
∂f
= 3 x 2 sen y + g ' ( x ) = 1 + 3 x 2 sen y
∂x
g '( x) = 1
de donde se obtiene
∫ dx = x
e integrando
g ( x) =
por lo que
f ( x , y ) = x 3 sen y + x
x 3 sen y + x = C
es la solución general.
⎣ y cos ( x y ) + e ⎦ dx + N ( x , y ) d y = 0
⎡ x
⎤
RESOLUCIÓN
M ( x , y ) = y cos ( x y ) + e x
44
y bajo la condición necesaria y suficiente establecida por la ecuación ( A ) se tiene
∂N
= − x y sen ( x y ) + cos ( x y ) ………….( B )
∂y
N (x , y ) =
∫ ( − x y ) sen ( x y ) dx
de aquí tenemos dos integrales, a saber
I1 =
I2 =
∫ cos ( x y ) dx
Para I 1 se tiene una integral por partes donde
u = − xy , du = − y dx
∫ sen ( x y ) dx =
− 1
y dv = sen ( x y ) dx , v = cos ( x y )
y
∫ cos ( x y ) dx
Entonces
I 1 = x cos ( x y ) −
I 1 = x cos ( x y ) − sen ( x y )
1
y
∫
Para I 2 se tiene una integral inmediata
I2 = cos ( x y ) d x = sen ( x y )
1
y
De lo anterior se tiene
N ( x , y ) = x cos ( x y ) + h ( y )
que es la función general que garantiza que la ecuación diferencial dada sea exacta.
45
( )
7) Determine si la ecuación diferencial
cos θ dr − r sen θ − e θ d θ = 0
RESOLUCIÓN
M (r , θ ) d r + N (r , θ ) d θ = 0
La ecuación diferencial es de la forma
M ( r , θ ) = cos θ
donde se identifican a las funciones M y N
(
N ( r , θ ) = r sen θ − e θ )
La condición necesaria y suficiente para que la ecuación diferencial sea exacta es que
∂M ∂N
=
∂θ ∂r
por lo que se procede a obtener las derivadas parciales
∂M
= − sen θ
∂θ
∂N
= − sen θ
∂r
Ahora bien, si es exacta, se asegura que existe una función f (r , θ ) tal que su
(r , θ ) dr + N ( r , θ ) d θ
o bien
df = M
por lo que
46
∂f ∂f
= M (r , θ) = N (r , θ )
∂r ∂θ
,
∫ M ( r , θ ) dr + h ( θ )
como se muestra enseguida
f (r , θ ) =
f (r , θ) =
∫ cos θ dr + h ( θ )
f ( r , θ ) = r cos θ + h ( θ ) ……………………( A )
∂f
= N ( r , θ ) ……………………( B )
∂θ
de donde
∂f
∂θ
=
∂
∂θ
( r cos θ + h ( θ ) )
∂f
= − r sen θ + h ' ( θ )
∂θ
entonces, de acuerdo a la ecuación ( B ) se tiene
despejando h ' ( θ )
h ' (θ ) = e θ
∫ e dθ
posteriormente se integra
h (θ ) = θ
h (θ ) = e θ
f ( r , θ ) = r cos θ + e θ
47
Finalmente, se tiene que f (r , θ ) = C es la forma de la solución general de la ecuación
48
EJERCICIOS PROPUESTOS
(2 y ) ( )
1) Resuelva la ecuación diferencial
2
+ 3 x y − 2 y + 6 x dx + x 2 + 2 x y − x d y = 0
Solución.
x 2y 2 + x3y − x2y + 2 x3 = C
⎡ π⎤
⎢⎣ 0 ,
2 ⎥⎦
Solución.
49
( ) ( )
3) Resuelva la ecuación diferencial x 2 + y 2 + 1 dx + x 2 − 2 x y d y = 0
Solución.
x − − + y = C
y2 1
x x
Solución.
y2 = k
Trazando para k = ± 2 , ± 1 , ±
1
, 0 , se obtiene
2
50
(x + ) ( )
5) Resuelva la siguiente ecuación diferencial homogénea
y 2 − x y y' = y ; y − 3 e− 4 = e− 4
Solución.
⎛ x⎞
− 2 ⎜1 − ⎟
1
y = e ⎝ y⎠
2
6) Si la ecuación propuesta no es exacta obtenga un factor integrante que dependa de una sola
(x )
de las variables y posteriormente resuelva la ecuación diferencial
+ y 2 dx − 2 x y d y = 0
Solución.
− + Ln x = C
y2
x
(x y + )
7) Resuelva la ecuación diferencial
y 2 + x 2 dx − x 2 d y = 0
y = x tan ( Ln C x )
Solución.
F (x , y )
(y ) ( )dy
soluciones
2
e 2 x + 1 dx + 2 y e 2 x − e x
= 0
51
Solución.
y2e x
− e −x
− y = C
( )
9) Resuelva la ecuación diferencial
x d y − y dx − 1 − x 2 dx = 0
Solución.
y = C x − 1 − x2
Solución.
Ln ( x + 1) + 2 y = 2
Solución.
x y2 + x2 − y3 = C
= + 2
dy y
12) Resuelva la ecuación diferencial
dx x
Solución.
= Ce
y
2 x
x
52
x d y + y dx + ( x + y ) dx = 0
13) Resuelva la ecuación diferencial
Solución.
x 3 + 3 x 2y = C
Solución.
y = Cex + 1
Solución.
x = Ce
y2
2x2
= + x + y
dy y 2 2
16) Resuelva la ecuación diferencial
dx x
Sugerencia: hacer el cambio de variable y = vx .
Solución.
⎛x2 ⎞
y = x tan ⎜ + C⎟
⎝ 2 ⎠
53
(x + ) ddxy = y
17) Resuelva la ecuación diferencial
y 2
− xy
Solución.
− 2 − xy
− Ln y = C
2
y
y
r cos θ d θ + ( r − sen θ ) dr = 0
18) Resuelva la ecuación diferencial
Solución.
sen θ + Ln r = C
1
r
( 2 x − y + 1) 2
20) Resuelva la ecuación diferencial
=
dy
dx 2x − y
Sugerencia: Realice el cambio de variable u = 2 x − y .
( )
Solución.
− Ln ( 2 x − y ) 2 + 1 − x = C
1
2
54
TEMA 3
y ( x) = C e − 1
sen 2 x
2
RESOLUCIÓN
yh = Ce = −1
sen 2 x
2 y y p
y h' = ( 2) C e
sen 2 x
cos 2 x 2
2
y h ' = C cos 2 x e
sen 2 x
2
a la no homogénea, se tiene
+ P ( x)C = 0
sen 2 x sen 2 x
C cos 2 x e 2 e 2
55
De esta ecuación se obtiene P ( x )
C ( cos 2 x ) e
P ( x) = −
sen 2 x
2
sen 2 x
Ce 2
P ( x ) = − cos 2 x
De donde resulta
y p' = 0
( 0) + P ( x ) ( − 1) = Q ( x ) ; pero P ( x ) = − cos 2 x
tiene
(− cos 2 x ) ( − 1) = Q ( x )
Q ( x ) = cos 2 x
Finalmente
− y = x 2 cos 2 x
dy
x
dx
RESOLUCIÓN
y ' + P ( x ) y = Q ( x ) ……………..(A)
variables, que puede asociarse con la expresión
56
Por lo anterior, se observa la conveniencia de normalizar la ecuación diferencial dada.
Dividiendo entre x se obtiene.
− y = x cos 2 x
dy 1
dx x
y ya es posible identificar las funciones
P ( x) = −
1
Q ( x ) = x cos 2 x
x
Debe tenerse presente que la solución general de (A) está dada por
yG = yh + y p
yh = Ce ∫
−
∫ Q ( x) e ∫
P( x) d x
= e ∫
− P( x) d x P ( x ) dx
y p
∫ P( x) dx ∫ x dx = − Ln x = Ln ⎜⎝ x ⎟⎠
Enseguida se procede a calcular las integrales anteriores requeridas:
⎛1⎞
= −
1
−
∫ P( x) dx = ∫ x dx = Ln x
1
∫ Q ( x ) e ∫ dx = ∫ x cos 2 x e
P( x) d x ( ) − Ln x dx
dx
∫
⎛1⎞
= x cos 2 x ⎜ ⎟ dx
⎝ x⎠
=
∫ cos 2 x dx =
1
sen 2 x
2
A continuación se tiene y h y y p , esto es,
yh = Ce ∫ = C e Ln x = C x
−
∫
P( x) d x
yp = C e ∫ Q ( x) e ∫
− P( x) d x P( x) d x
57
⎛1 ⎞
= ⎜ sen 2 x ⎟ e Ln x
⎝2 ⎠
y p
⎛1 ⎞
= ⎜ sen 2 x ⎟ x
⎝2 ⎠
y p
y = Cx +
x
sen 2 x
2
y ' + P ( x) y = Q ( x) ,
3) Sabiendo que la solución general de una ecuación diferencial de la forma
es
x3 − 3x
( ) ( )
y' = +
c
3 x + x 2
3 x + x 2
RESOLUCIÓN
∫
p
yh = ce ∫ = e ∫ Q ( x) e ∫
− P( x) d x − P( x) d x P ( x ) dx
donde y y p dx por lo que de la
58
donde y h es la solución de la homogénea asociada y y p es una solución particular
de la no homogénea.
y ' + P ( x ) y = 0 …………….( C )
Así, para la ecuación homogénea asociada
− c ( 2 x + 1)
( )
y h' =
3 x2 + x 2
− c ( 2 x + 1)
la cual se sustituye en la ecuación ( C )
+ P ( x)
( ) ( )
= 0
c
3 x2 + x 2
3 x2 + x
P ( x) =
(2 x + 1)
(x 2
+ x )
Por otro lado, asumiendo que y es una solución particular de
y ' + P ( x) y = Q ( x)
p
∫
y se sabe que la forma de y p está dada por
= e ∫ Q ( x) e ∫
− P( x) d x P( x) d x
y p dx
∫ ∫(
(2 x + 1)
( )
⎛ 1 ⎞
)
− P( x) dx = − dx = − Ln x 2 + x = Ln ⎜ 2 ⎟ …….( D )
x + x ⎝ x + x⎠
2
∫ ∫ ( )
2x + 1
P( x) dx = dx = Ln x 2 + x
x + x
2
……………...( E )
entonces y es
∫
p
( )
⎛ 1 ⎞
Q ( x) e
Ln ⎜ 2 ⎟ Ln x 2 + x
y p = e ⎝x + x⎠
dx
simplificando se obtiene
59
=
∫ (
Q ( x ) x 2 + x dx ………………….( F ) )
1
x + x
y p 2
Por otro lado, se tiene el valor conocido y p dado por ( B ), entonces se pueden igualar ( B )
∫ ( )
y (F)
x 3 − 3x
Q ( x ) x 2 + x dx =
( )
1
x + x
2
3 x2 + x
(x ) (x )
∫
de donde
( )
− 3x + x
Q ( x ) x + x dx =
( )
3 2
3 x2 + x
2
∫ (
Q ( x ) x 2 + x dx = ) (
x 3 − 3x )
1
3
( ) ( )
Derivando ambos miembros de la ecuación se tiene
Q ( x) x 2 + x = 3x 2 − 3
1
( )
3
= 3 x2 − 1
1
3
= x2 − 1
(x + 1)( x − 1)
Finalmente, se despeja a la función de interés
x2 − 1
Q ( x) =
(x )
=
2
+ x x ( x + 1)
x − 1
Q ( x) =
x
− 2 x = t Ln 3 t
dx
t
dt
RESOLUCIÓN
60
xG = xh + x p
∫
expresión que involucra factores integrantes:
xG = Ce ∫ + e ∫ Q (t ) e ∫
− P (t ) d t − P (t ) d t P (t ) d t
dt
x ' + P (t ) x = Q (t )
transformarla posteriormente en la forma normalizada
− 2 x = t Ln ( 3 t )
dx
t
dt
Dividiendo entre “ t “
− x = Ln ( 3 t )
dx 2
dt t
x ' + P (t ) x = Q (t )
que es una ecuación diferencial de la forma
− x = 0
dx 2
dt t
Efectuamos separación de variables
=
dx 2
x
dt t
= 2
dx dt
x t
∫ ∫
Integramos
= 2
dx dt
x t
61
Ln x + C 1 = 2 Ln t + C 2
( )
Que también se puede expresar en forma explícita; aplicamos propiedades según se indica
Ln ( x ) = Ln t 2 + Ln C
Ln ( x ) = Ln C t 2 ( ) , y al aplicar la exponencial
xh = Ct2
∫
Sabemos que una solución particular de la no homogénea tiene la forma:
= e ∫ Q (t ) e ∫
− P (t ) d t P (t ) d t
x p dt
y de la ecuación normalizada
− 2
P (t ) = ; Q ( t ) = Ln ( 3 t )
t
∫ ∫
enseguida se integran las funciones respectivas:
− P(t ) dt = 2 = 2 Ln t = Ln t 2
dt
t
⎛ 1 ⎞
e − ∫ P (t ) d t
= e
( ) = t 2 ; e ∫ P (t ) d t dt = e Ln ⎜⎝ t
Ln t 2 2 ⎟
⎠
=
1
t2
∫ ∫ ∫ ( ) dt
⎛ 1 ⎞
Q (t ) e ∫ = Ln ( 3t ) ⎜ 2 ⎟ dt = Ln ( 3 t ) t − = I
P (t ) d t
⎝t ⎠
2
t −1 − 1
du = dt = dt v = =
3 1
− 1
∫ ∫t
3t t t
− 1
I = Ln ( 3 t ) = ⋅ dt = − Ln ( 3 t ) + −2
1 1 1
dt
t t t t
62
− 1
Ln ( 3 t ) + = − Ln ( 3 t ) −
t −1
I =
1 1
t −1 t t
Entonces la solución particular es
= t 2 ⎜ − Ln ( 3 t ) − ⎟
⎛ 1 1⎞
⎝ t t⎠
x p
y' + x 2 y = x 2 1 − e
1 1 x2
1 2
x
RESOLUCIÓN
⎛ ⎞2
1
y ' + x y = x ⎜1 − e ⎟
x2
⎜ ⎟
2
⎝ ⎠
y ' + P ( x) y = Q ( x)
es una ecuación diferencial de la forma
∫ Q ( x) e ∫
enseguida
yG = C e ∫ + e ∫
− P( x) d x − P( x) d x P( x) d x
dx
donde
⎛ ⎞
P ( x) = x Q ( x ) = x ⎜1 − e
1
⎟
x2 2
⎜ ⎟
⎝ ⎠
, 2
63
∫ P( x) d x ∫ ; e ∫
−
− = − x dx = − = e
−
x2
x2 P ( x ) dx
2
∫ P( x) d x = ∫
2
x dx = ; e∫ = e
x2
x2 P ( x ) dx
2
2
∫ ∫ ∫u
⎛ ⎞2 ⎛ x ⎞
1
dx = x ⎜1 − e 2 ⎟ ⎜e 2 ⎟ dx =
2 2
x
⎜ ⎟ ⎜ ⎟
P( x) d x n
⎝ ⎠ ⎝ ⎠
Q( x) e du
u = 1 −e
x2
2
du = − e ( x ) dx = −xe
x2 x2
2 2 dx
Por lo que
∫
2⎛ ⎞2
3
u n du = − ⎜ 1 − e 2 ⎟
2
x
3⎜ ⎟
⎝ ⎠
Entonces, al sustituir en la forma de la solución general se obtiene
⎡ ⎤
⎢ 2⎛ ⎞ 2⎥
3
− −
yG = Ce + ⎢ − ⎜⎜ 1 − e ⎟ ⎥
x2 x2 x2
⎟
⎢ 3⎝ ⎠ ⎥
2 e 2 2
⎣ ⎦
Reordenando y factorizando se tiene finalmente
⎡ ⎤
⎢ 2⎛ ⎞ 2⎥
3
−
yG = e ⎢C − ⎜⎜ 1 − e ⎟ ⎥
x2 x2
⎟
⎢ ⎝ ⎠ ⎥
2 2
⎣ ⎦
3
RESOLUCIÓN
64
Se trata de una ecuación diferencial de primer orden, no homogénea y coeficiente
variable; para resolverla es necesario normalizarla por lo que se divide entre la variable
x , de donde se obtiene
y' − y = xe−
1 3x
P ( x) = − Q ( x) = x e −
1 3x
y
x
∫
La solución general de una ecuación de la forma anterior está dada por
yG = Ce ∫ + e ∫ Q ( x) e ∫
− P( x) d x − P( x) d x P( x) d x
dx
donde y h = C e ∫
−
∫
P( x) d x
es la solución de la homogénea asociada y
= e ∫ Q ( x) e ∫
− P( x) d x P( x) d x
y p dx es una solución particular de la no homogénea.
∫ ∫
Efectuamos el cálculo de las diferentes integrales, según se indica.
⎛1⎞
P( x) dx = − = − Ln x = Ln ⎜ ⎟
dx
⎝ x⎠
∫ ∫
x
− P( x) dx = = Ln x
dx
x
e ∫ = e Ln x = x
− P( x) d x
⎛1⎞
e∫
Ln ⎜ ⎟
= e ⎝ x⎠
=
P( x) d x 1
∫ Q ( x) e ∫ ∫ xe ∫e
x
dx = − 3x
dx = − 3x
P( x) d x 1
dx
x
= −
1 − 3x
e
3
Así, se obtiene
yh = C x = −
x − 3x
y y p e
3
Y la solución general es
yh = C x −
x − 3x
e
3
65
7) Obtenga la solución general de la ecuación diferencial
+ 2 y =
dy x cos 3 x
dx x + 1 x2 + 1
RESOLUCIÓN
y ' + P ( x) y = Q ( x)
forma
∫ Q ( x) e ∫
su solución general está dada por
yG = Ce ∫ + e ∫
− P( x) d x − P( x) d x P( x) d x
dx
∫
donde
yh = Ce ∫ = e ∫ Q ( x) e ∫
− P( x) d x − P( x) d x P( x) d x
; y p dx
Para resolverla es conveniente identificar a las funciones que intervienen en las integrales,
es decir,
P ( x) =
x
x + 1
2
Q ( x) =
cos 3 x
x2 + 1
∫ ∫
A continuación efectuamos las operaciones necesarias para y h y y
( )
p
− P ( x ) dx = − dx = − Ln x 2 + 1
x 1
x + 1
2
2
⎛ ⎞
= Ln ⎜ ⎟
⎜ x2 + 1⎟
1
⎝ ⎠
∫ P ( x ) dx =
∫ dx = Ln x 2 + 1 ( )
x 1
x + 1
2
2
66
= Ln x 2 + 1
entonces
⎛ ⎞
Ln ⎜ ⎟
⎜ x2 + 1⎟
yh = Ce =
1
⎝ ⎠ C
x2 + 1
∫
( )
⎛ ⎞
Ln ⎜ ⎟
⎜ x2 + 1⎟ x2 + 1
= e
1
⎝ ⎠ cos 3 x Ln
x2 + 1
y p e dx
=
∫ x 2 + 1 dx
1 cos 3 x
x + 1 x + 1
y p
∫ cos 3 x dx
2 2
=
1
x + 1
y p
2
⎛ ⎞ sen 3 x
= ⎜ ⎟
⎜ x2 + 1⎟
1
⎝ ⎠
y p
3
y = +
C sen 3 x
x2 + 1 3 x2 + 1
⎛π ⎞
Sujeta a la condición inicial y ⎜ ⎟ =
⎝4⎠
2
RESOLUCIÓN
Se tiene una ecuación diferencial de primer orden, que puede escribirse de la siguiente
67
y ' + P ( x) y = Q ( x)
la cual es una ecuación de la forma
Q ( x ) = sen x
∫
variables, cuya solución general es de la forma
y = Ce ∫ + e ∫ Q ( x) e ∫
− P ( x ) dx − P ( x ) dx P ( x ) dx
dx
donde
yh = Ce ∫
−
∫ Q ( x) e ∫
P ( x ) dx
= e ∫
− P ( x ) dx P ( x ) dx
y p dx
∫ ∫ ∫
A continuación se calculan las integrales necesarias, esto es
e ∫
Ln ( sec x )
= e = sec x
− P ( x ) dx
por lo que
y h = C sec x
∫ P ( x ) dx = − ∫ tan x d x = Ln ( cos x )
Por otro lado
e∫
Ln ( cos x )
= e = cos x
P ( x ) dx
∫ ∫
( sen x ) 2
entonces
Q ( x) e ∫ P ( x ) dx
dx = sen x cos x dx =
2
( sen x ) 2
por lo que
y p = sec x
2
y la solución general es
68
y = C sec x + sec x
( sen x ) 2
2
⎛π ⎞
Aplicando la condición inicial y ⎜ ⎟ =
⎝4⎠
2 se tiene
⎛ ⎛ π ⎞⎞
⎜ sen ⎜ ⎟ ⎟
2
⎛π ⎞ ⎛π ⎞ ⎝ 4 ⎠⎠
2 = C sec ⎜ ⎟ + sec ⎜ ⎟ ⎝
⎝4⎠ ⎝4⎠ 2
⎛ 1 ⎞
⎜ ⎟
2
2 ⎝
2⎠
2 = C 2 +
2
⎛1⎞
2 = C 2 + 2 ⎜ ⎟ ⇒ C =
3
⎝4⎠ 4
( sen x ) 2
Finalmente para la condición inicial dada
y = sec x + sec x
3
4 2
69
3.2 SOLUCIÓN DE LA ECUACIÓN DIFERENCIAL LINEAL
HOMOGÉNEA.
RESOLUCIÓN
(D )
constantes. Inicialmente la expresamos en términos del operador diferencial
2
− D − 2 y = 0
RESOLUCIÓN
(D )
constantes, que expresada en términos del operador diferencial es
2
− 6D + 9 y = 0
70
y = C 1 e 3x + C 2 xe 3x
RESOLUCIÓN
(D )
constantes; expresada en términos del operador diferencial es
2
− 3D + 4 y = 0
y = C1e x + C 2 e 2 sen
3 3
x
2 cos
7 x 7
x
2 2
RESOLUCIÓN
(D )
expresada en términos del operador diferencial es
3
− 2D 2 − D + 2 y = 0
71
λ 3 − 2λ 2 − λ + 2 = 0
Se procede a obtener sus raíces, las cuales son λ 1 = 1 , λ 2 = 2 y λ 3 = − 1
Las raíces son reales y distintas, por lo que la solución general es
y = C 1 e x + C 2 e2x + C 3 e −x
RESOLUCIÓN
(D )
en términos del operador diferencial es
3
− 3D 2 + 3D − 1 y = 0
Su ecuación característica es
λ 3 − 3λ 2 + 3λ − 1 = 0
siendo sus raíces λ 1 = λ 2 = λ 3 = 1
Se tiene el caso de raíces reales repetidas por lo que la solución general de la ecuación
diferencial es
y = C 1 e x + C 2 xe x + C 3 x 2 e x
RESOLUCIÓN
72
(D 4
)
+ 2D 2 + 1 y = 0
y su ecuación característica es
λ 4 + 2 λ 2 + 1 = 0 …………..( A )
Se determinan sus raíces reales mediante un cambio de variable adecuada que transforme la
ecuación de cuarto grado en una ecuación de segundo grado, esto es
w = λ 2 , w2 = λ 4
por lo que se tiene
w 2 + 2w + 1 = 0
(w + 1) 2 = 0
73
3.3 SOLUCIÓN DE LA ECUACIÓN DIFERENCIAL LINEAL NO
HOMOGÉNEA.
RESOLUCIÓN
= x2 − 1
d 2y
dx 2
Al integrar esta ecuación se tiene
= − x + C1
dy x3
dx 3
∫
nuevamente se integra
⎛x3 ⎞
y = ⎜ − x + C 1 ⎟ dx
⎝ 3 ⎠
de donde se obtiene
y = − + C1 x + C2
x4 x2
12 2
que constituye la solución general de la ecuación diferencial ( A )
Segundo método:
74
La ecuación dada es de segundo orden, coeficientes constantes y no homogénea y su
solución general es de la forma
yG = yh + y p
la no homogénea.
Para obtener y h se considera la ecuación
y '' = 0
que expresada en términos del operador diferencial es
D2y = 0
y cuya ecuación característica es
λ2 = 0
De esta última obtenemos sus raíces, a saber,
λ1 =λ 2 = 0
Se tienen valores característicos repetidos nulos, por lo que la solución correspondiente es
yh = C 1 + C 2 x
Q ( x) = x 2 − 1
p
Al aplicar este operador a ambos miembros de la ecuación (A ), esta última en términos del
( ) ( )
operador D se tiene
D3 D2 y = D3 x2 − 1
lleva a la expresión
D5 y = 0
75
(y h)
(y p)
Esta solución involucra la solución de la ecuación homogénea asociada y una
solución es de la forma
y h = C 1 + C 2 x …………………….………….( B )
y p = C 3 x 2 + C 4 x 3 + C 5 x 4 …………………( C )
obtener
y p ' , y p ''
y p' = 2 C 3 x + 3 C 4 x 2 + 4 C 5 x 3
y p '' = 2 C 3 + 6 C 4 x + 12 C 5 x 2
6C 4 = 0
12 C 5 = 1
que al resolver permite obtener
C3 = − , C4 = 0 , C5 =
1 1
2 12
Estos valores se sustituyen en la ecuación ( C ) con lo que se tiene
= − x +
1 2 1 4
y p x
2 12
Finalmente, la solución general es
yG = C1 + C2 x − x +
1 2 1 4
x ……………………….( D )
2 12
la cual es igual a la obtenida por el primer método.
76
yG = C2 − x2 +
1 3
x ……………………..( E )
2
Esta ecuación junto con ( D ), al aplicar las condiciones iniciales constituyen un sistema de
dos ecuaciones con dos incógnitas: C 1 y C 2 .
1 = C1
2 = C2
y = 1 + 2x − x +
1 2 1 4
x
2 12
que constituye una solución particular de la ecuación diferencial ( A ) .
ecuación ( D )
2 = − + C1 + C2
1 1
12 2
2 = − + C1 + C2
5
12
finalmente
2 + = C1 + C2
5
12
= C 1 + C 2 …………….…………( F )
29
12
Asimismo
5 = C 1 + 3C 2 x − ( 3) 2 + ( 3) 4
1 1
2 12
5 = C 1 + 3C 2 x − +
9 81
2 12
5 + − = C 1 + 3C 2
9 81
2 12
77
finalmente
= C 1 + 3 C 2 …………………………….( G )
33
12
Enseguida se resuelve el sistema formado por las ecuaciones ( F ) y ( G ), de donde se tiene
C1 = C2 =
27 1
y
12 6
Al sustituir estos valores en ( D ), se obtiene la solución
y = + x − x2 +
27 1 1 1 4
x
RESOLUCIÓN
Se tiene la ecuación
y '' + 3 y ' − 10 y = cot 3 x sen 3 x
Q ( x)
que es una ecuación diferencial no homogénea, de segundo orden y de coeficientes
constantes; asimismo se tiene una función para la cual aparentemente no existe
Q ( x ) = cos 3 x
78
yG = yh + y p
(D )
y en términos del operador diferencial es
2
+ 3 D − 10 y = 0
Q ( x ) = cos 3 x
un operador que anule a la función
el cual resulta
P ( x) = D 2 + 9
(
P ( D ) D 2 + 3 D − 10 y = P ( D ) ( cos 3 x ) )
(D 2
+ 9 ) (D 2
)
+ 3 D − 10 y = 0 ……………………..( A )
que es una ecuación diferencial de orden mayor a la original, cuya ecuación característica y
(λ )( )
raíces respectivas son
2
+ 9 λ 2 + 3 λ − 10 = 0
⇒ λ 1 = − 5 , λ 2 = 2 , λ 3 = 3i , λ 4 = − 3i
Así, la solución general de ( A ) es
y G = C 1 e − 5 x + C 2 e 2 x + C 3 cos 3 x + C 4 sen 3 x
y h = C 1 e −5x + C 2 e 2x
y p = C 3 cos 3 x + C 4 sen 3 x
79
En y p se tienen justamente los coeficientes indeterminados C3 y C4 cuyos
la ecuación no homogénea
y p ' = − 3 C 3 sen 3 x + 3 C 4 cos 3 x
− 9 C 3 − 19 C 4 = 0
cuya solución es
C3 = − ; C4 =
19 9
442 442
Entonces y p resulta
= − cos 3 x +
19 9
y p sen 3 x
442 442
por lo que la solución general es
y G = C 1 e −5x + C 2 e 2x − cos 3 x +
19 9
sen 3 x
442 442
RESOLUCIÓN
80
La ecuación diferencial es no homogénea, lineal, de segundo orden y de coeficientes
constantes; su solución general es de la forma
yG = yh + y p
(D )
y en términos del operador diferencial
2
− 4D + 8 y = 0
su ecuación característica es
λ 2 − 4λ + 8 = 0
cuyas raíces son
λ 1 = 2 + 2i y λ 2 = 2 − 2i
De acuerdo a las raíces obtenidas resulta
y h = C 1 e 2 x cos 2 x + C 2 e 2 x sen 2 x
⎡ ⎤ ⎡u ' ( x ) ⎤ ⎡ 0 ⎤
⎢ 2x ⎥ =
2 e 2 x ( cos 2 x + sen 2 x ) ⎦⎥ ⎣ v ' ( x ) ⎦ ⎣ ( )⎦
⎢ ⎥ ⎢Q x ⎥
⎣⎢ 2 e ( cos 2 x − sen 2 x )
e 2 x cos 2 x e 2 x sen 2 x
Para resolver este sistema por regla de Cramer es necesario obtener los determinantes
siguientes:
⎡ ⎤
∆1 = ⎢ 2x ⎥
⎢⎣ 2 e ( cos 2 x − sen 2 x ) 2 e 2 x ( cos 2 x + sen 2 x ) ⎥⎦
e 2 x cos 2 x e 2 x sen 2 x
81
∆ 1 = 2 e 4 x cos 2 2 x + 2 e 4 x sen 2 x cos 2 x + 2 e 4 x sen 2 2 x − 2 e 4 x cos 2 x sen 2 x
∆ 1 = 2e 4 x
⎡ 0 ⎤
∆2 = ⎢ ⎥ = − 4 e 4 x sen 2 x sec 2 x + 4 e 4 x sen 2 x csc 2 x
⎣⎢Q ( x ) 2 e ( cos 2 x + sen 2 x ) ⎦⎥
e 2 x sen 2 x
2x
∆ 2 = − 4 e 4x sen 2 x ⋅ + 4 e 4 x sen 2 x ⋅
1 1
cos 2 x sen 2 x
∆ 1 = − 4 e 4 x tan 2 x + 4 e 4 x
⎡ 0 ⎤
∆3 = ⎢ 2x ⎥ = 4 e sec 2 x cos 2 x − 4 e cos 2 x csc 2 x
⎣⎢ 2 e ( cos 2 x − sen 2 x ) Q ( x ) ⎥⎦
e 2 x cos 2 x 4x 4x
∆3 = 4 e 4 x cos 2 x − 4 e 4 x cos 2 x
1 1
cos 2 x sen 2 x
∆ 3 = 4 e 4 x − 4 e 4 x cot 2 x
∆2 − 4 e 4 x ( tan 2 x − 1)
u '( x) = = = − 2 ( tan 2 x − 1)
∆1 2e4x
∆3 4 e 4 x ( 1 − cot 2 x )
v '( x) = = = 2 ( 1 − cot 2 x )
∆1 2e4x
∫ ( tan 2 x − 1) dx
obtenidas
u ( x) = − 2
∫
= − 2 tan 2 x dx + 2 dx
∫
= − 2
∫ cos 2 x dx + 2 x
sen 2 x
= Ln ( cos 2 x ) + 2 x
v ( x) = 2
∫ (1 − cot 2 x ) dx
∫
= 2 dx − 2
∫ sen 2 x dx
cos 2 x
= 2 x − ln ( sen 2 x )
82
De lo anterior se establece la solución y
y p
+ e 2 x sen 2 x ⎡⎣ 2 x − ln ( sen 2 x ) ⎤⎦
y '' + y ' − 2 y = ( 2 + 6 x ) e x
indeterminados
RESOLUCIÓN
(D )
que en términos del operador diferencial es
2
+ D − 2 y = 0
Q ( x ) , a saber
mediante el método de los coeficientes indeterminados, es necesario identificar a la función
83
Q ( x) = (2 + 6 x) e x
= 2ex + 6 xex
P ( D ) = ( D − 1) 2
P ( D)Q ( x) = 0
(D (
− 1) 2 2 e x + 6 x e x )=0
Al aplicar P ( D ) a la ecuación no homogénea original se tiene
(
P ( D) D 2 + D − 2 y = P ( D) Q ( x))
(D ( )
− 1) 2 D 2 + D − 2 y = 0 ……………………..( A )
( )
característica correspondiente es
(λ − 1) 2 λ 2 + λ − 2 = 0
y h = C 1 e x + C 2 e −2x
y p = C3 xex + C4 x2 ex
y p ' y y p '' .
yp ' = C 3 x e x + C 3 e x + C 4 x 2 e x + 2 C 4 x e x
y p '' = C 3 x e x + C 3 e x + C 3 e x + C 4 x 2 e x + 2 C 4 x e x + 2 C 4 x e x + 2 C 4 e x
y p '' = C 3 x e x + 2 C 3 e x + C 4 x 2 e x + 4 C 4 x e x + 2 C 4 e x
84
enseguida se sustituye en la ecuación diferencial original, esto es,
C 3 x e x + 2 C 3 e x + C 4 x 2 e x + 4C 4 x e x
+ 2C 4 e x + C 3 x e x + C 3 e x + C 4 x 2 e x
+ 2C 4 x e x − 2C 3 x e x − 2C 4 x 2 e x = 2e x + 6 x e x
e x( 3 C 3 + 2 C 4 ) + 6 C 4 x e x = 2 e x + 6 x e x
6C 4 = 6
cuya solución es C 3 = 0 , C 4 = 1
y p = x2ex
RESOLUCIÓN
(y h)
la cual es una ecuación de segundo orden, no homogénea de coeficientes constantes. La
(y p) .
solución general está dada por la solución de la ecuación homogénea asociada y
85
y '' − 2 y ' + y = 0
(D )
que expresada en términos del operador diferencial es
2
− 2D + 1 y = 0
Las raíces de esta ecuación son λ 1 = λ 2 = 1 , raíces reales repetidas; para este caso la
solución es
yh = C1ex + C2 xex
que corresponde a la solución de la ecuación homogénea asociada a la ecuación homogénea
asociada a la ecuación no homogénea ( A ).
Para obtener y p es conveniente emplear el método de coeficientes indeterminados, pues
( ) ( )
miembros
P ( D ) D 2 − 2 D + 1 y = P ( D ) x 2− 1
( )
D 3 D 2 − 2 D + 1 y = D 3 x 2− 1( )
La aplicación de P ( D ) = D 3 en el segundo miembro determina que el término Q ( x )
( )
se anule, por lo que se tiene
D 3 D 2 − 2 D + 1 y = 0 …………………….( B )
86
(
λ 3 λ 2 − 2λ + 1 = 0 )
λ 3( λ − 1) 2 = 0
y sus valores característicos son
λ1 = λ2 = 1 , λ3 = λ4 = λ5 = 0
Así, la solución general de ( B ) es
y = C1ex + C2 xex + C3 + C4 x + C5 x2
yh = C1ex + C2 xex
y p = C3 + C4 x + C5 x2
yp ' = C 4 + 2 C 5 x
y p '' = 2 C 5
Sustituyendo en ( A )
2C 5 − 2C 4 − 4C 5 x + C 3 + C 4 x + C 5 x 2 = x 2 − 1
87
y ''' − 3 y '' + 9 y ' + 13 y = − 6 e − x + 3 sen 3 x
RESOLUCIÓN
(D )
La ecuación homogénea asociada correspondiente
3
− 3 D 2 + 9 D + 13 y = 0 ………………..( B )
⎩ λ 3 = 2 − 3i
por lo que la solución de ( B ) es
y h = C 1 e − x + e 2 x ( C 2 cos 3 x + C 3 sen 3 x )
( ) (D
el operador anulador es
P ( D) = D 2 + 9 + 1)
(D ) (D ( )
por lo que, al aplicarlo en la ecuación diferencial ( A ) se tiene
2
+ 9 + 1) D 3 − 3 D 2 + 9 D + 13 y = 0 ………….( C )
(λ ) ( )
cuya ecuación característica es
2
+ 9 ( λ + 1) λ 3 − 3 λ 2 + 9 λ + 13 y = 0
88
donde se identifican yh y y p , esto es
y h = C 1 e − x + C 2 e x cos 3 x + C 3 e 2 x sen 3 x
+ 6 C 4 e − x + 27 C 5 cos 3 x + 27 C 6 sen 3 x − 9 C 4 x e − x + 9 C 4 e − x
+ 13 C 6 sen 3 x = − 6 e − x + 3 sen 3 x
al reducir términos semejantes
18 C 4 e − x + 40 C 5 cos 3 x + 40 C 6 sen 3 x = − 6 e − x + 3 sen 3 x
de donde se obtienen los valores buscados
C4 = − C5 = 0 C6 =
1 3
, ,
3 40
al sustituir estos valores en la forma de y p se tiene
= − x e −x +
1 3
y p sen 3 x
3 40
Finalmente, la solución general es
+ A + B y = Q (t )
d 2y dy
2
dt dt
89
de la que se conoce una solución particular y p y
= −
3
y p t cos 2 t
4
y h = C 1 cos 2 t + C 2 sen 2 t
Obtenga
Q (t )
a) los valores de las constantes A y B
b) la función
RESOLUCIÓN
+ A + By = 0
d 2y dy
2
dt dt
(D )
y la ecuación diferencial homogénea expresada en términos del operador diferencial resulta
2
+ 4 y = 0
o bien
+ 4y = 0
d 2y
dt 2
de donde los valores constantes son A = 0 , B = 4
Otra forma para determinar estos valores consiste en obtener yh ' , yh '' y posteriormente
90
y h = C 1 cos 2 t + C 2 sen 2 t …………………………….( I )
+ 4y = 0
d 2y
dt 2
por lo que las constantes buscadas son A = 0 , B = 4 .
diferencial es de la forma
+ 4 y = Q ( t ) ……………………..( IV )
d 2y
dt 2
⎛ 3 ⎞
3 t cos 2 t + 3 sen 2 t + 4 ⎜ − t cos 2 t ⎟ = Q(t )
⎝ 4 ⎠
3 sen 2 t = Q(t )
91
8) Resuelva la siguiente ecuación diferencial y '' + y = 2 cos x utilizando el método de
RESOLUCIÓN
(D )
en términos del operador diferencial es
2
+ 1 y = 0
y su ecuación característica
λ2 + 1 = 0 ⇒ λ2 = − 1
y sus valores característicos λ 1 = i , λ 2 = − i de donde
y h = C 1 cos x + C 2 sen x
= u ( x ) cos x + v ( x ) sen x
de parámetros, se tiene
y
u ( x ) y v ( x ) satisfacen el sistema
p
v ' ( x ) = 2 cos 2 x
donde, al integrar
92
u ( x) = − 2
∫ cos x sen x dx = − sen 2 x
v ( x) = 2
∫ cos 2 x dx = x +
1
sen 2 x
2
y p ( ) ⎛1 ⎞
= − sen 2 x cos x + x + ⎜ sen 2 x ⎟ sen x
⎝2 ⎠
y p ( )
= − sen 2 x cos x + ( x + sen x cos x ) sen x
( ) (
= − sen 2 x cos x + x sen x + sen 2 x cos x )
= x sen x
C1 = 5
C2 = 2
Finalmente
y = 5 cos x + 2 sen x + x sen x
93
10) Sea la ecuación diferencial homogénea de coeficientes constantes
L ( y ) = 3 y '' + A y ' + B y = 0 , de la cual se conocen las soluciones y 1 = 4 e x / 3
y 2 = 3e 3x −
1 x/3
y e ; donde L es el operador diferencial.
RESOLUCIÓN
y2 ' = 9 e 3 x − y2 '' = 27 e 3 x −
x x
1 3 1
e , e3
9 27
Dado que y 1 y y 2 son soluciones de la ecuación diferencial, se sustituyen en ella.
Para y 1 se tiene
⎛4 x⎞ ⎛4 x⎞ ⎛ ⎞
3⎜ e 3 ⎟ + A⎜ e 3 ⎟ + B ⎜4 e 3 ⎟ = 0
x
⎜9 ⎟ ⎜3 ⎟ ⎜ ⎟
⎝ ⎠ ⎝ ⎠ ⎝ ⎠
e + A e3 + 4Be3 = 0
x x x
4 3 4
3 3
Multiplicando por 3
4 e 3 + 4 A e 3 + 12 B e 3 = 0 ………………….( 1 )
x x x
94
Para y 2 se tiene
⎛ 1 3⎞ ⎛ 1 3⎞ ⎛ 3x 1 3⎞
3 ⎜ 27 e − e ⎟ + A⎜ 9 e − e ⎟ + B⎜3 e − e ⎟ = 0
x x x
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ ⎠ ⎝ ⎠ ⎝ ⎠
3x 3x
27 9 3
81 e 3 x − e + 9 Ae 3x − e + 3 B e3x − e = 0
x x x
1 3 A 3 B 3
9 9 3
Multiplicando por 9
− e + 81 A e − Ae + 27 B e − 3Be = 0 ………..( II )
x x x
3x 3 3x 3 3x 3
729 e
Debe resolverse el sistema formado por ( I ) y ( II ) .
Resolviendo para ( A) en ( I ) resulta
A = −3B − 1
Al sustituir en ( II ) y reducir términos semejantes se obtiene
648 = 216 B ⇒ B = 3
por lo que A = − 10
Así, la ecuación diferencial homogénea resulta ser
3 y '' − 10 y ' + 3 y = 0
yG = yh + y p
(3 D )
que en términos del operador diferencial es
2
− 10 D + 3 y = 0
95
Entonces, la solución de la homogénea asociada es
yh = C1e + C2e
1
x
3x 3
Q ( x) = 9
p
P ( D) = D
un operador anulador para esta función es
( )
Y al aplicarlo en ( III ) se obtiene
D 3 D 2 − 10 D + 3 y = 0 ………….( IV )
λ ( 3 λ 2 − 10 λ + 3 ) y = 0
cuya ecuación característica es
λ1 = 3 λ2 = λ3 = 0
1
, ,
3
así, la solución general de ( IV ) es
y = C 1 e 3x + C 2 e 3 + C3
1
x
por lo que
y p = C3
Ahora bien, si y p es una solución particular de ( III ) entonces debe satisfacerla por lo
yp ' = 0
y p '' = 0
y G = C 1 e 3x + C 2 e 3 + 3
1
x
96
Para aplicar las condiciones iniciales dadas es necesario derivar a la función que representa
a la solución general, esto es
yG ' = 3 C 1 e + C2e3
1
3x 1 x
3 = C1 + C2 + 3
− 8 = 3C 1 +
1
C2
3
al resolver este sistema se obtiene
C1 = − 3 , C2 = 3
sustituyendo estos valores en la solución general resulta finalmente
y = − 3e + 3e + 3
1
x
3x 3
RESOLUCIÓN
y' − y = x e −3x
1
P ( x) = − Q ( x) = x e − 3 x
1
y
x
∫
La solución general de una ecuación de la forma anterior está dada por
yG = C e ∫ + e ∫ Q ( x) e ∫
− P( x) d x − P( x) d x P( x) d x
dx
97
donde yh = C e ∫
− P ( x ) dx
∫
es la solución de la homogénea asociada
yp = e ∫ Q ( x) e ∫
− P( x) d x P( x) d x
dx
∫ ∫
Efectuamos el cálculo de las diferentes integrales, según se indica
⎛1⎞
P( x) dx = − = − Ln x = Ln ⎜ ⎟
dx
x ⎝ x⎠
−
∫ P( x) dx =
∫ = Ln x
dx
x
e ∫ dx = e Ln x = x
− P( x) d x
⎛1⎞
e∫ dx = e =
Ln ⎜ ⎟
P( x) d x ⎝ x⎠ 1
∫ Q ( x) e ∫ ∫ xe ∫e
x
dx = −3x
dx = −3x
P ( x ) dx 1
dx
x
= −
1 −3x
e
3
Así, se obtienen
yh = C x yp = −
x −3x
y e
3
y la solución general es
yG = C x −
x −3x
e
3
12) Obtenga la ecuación diferencial lineal de coeficientes constantes de menor orden cuya
solución general es
y ( x ) = e − 3 x ( C 1 cos 2 x + C 2 sen 2 x ) + C 3 + C 4 x + 5 e 4x
RESOLUCIÓN
De la solución general
y G = C 1 e − 3 x cos 2 x + C 2 e − 3 x sen 2 x + C 3 + C 4 x + 5 e 4x
98
se tiene la solución yh y la solución y p , esto es
y h = C 1 e − 3 x cos 2 x + C 2 e − 3 x sen 2 x + C 3 + C 4 x
yp = 5 e 4x
Lo anterior significa que la ecuación buscada es no homogénea. Por otro lado, al observar
los elementos de yh , es posible determinar las raíces que están involucradas; de los
términos con funciones exponenciales y trigonométricas se tiene
λ 1 = − 3 + 2i , λ 2 = − 3 − 2i
de los términos con funciones polinomiales se tiene
λ3 = λ4 = 0
Así, la ecuación característica correspondiente a yh está dada por
⎡⎣ λ − (− 3 + 2 i ) ⎤⎦ ⎡⎣λ − (− 3 − 2 i ) ⎤⎦ λ 2 = 0
(D )
en términos del operador diferencial
4
+ 6 D 3 + 13 D 2 y = 0
o bien
y ( IV ) + 6 y ''' + 13 y '' = 0
Entonces, la ecuación diferencial buscada se puede escribir como
y ( IV ) + 6 y ''' + 13 y '' = Q ( x ) …………………………( A )
yp = 5 e 4x
y p ' = 20 e 4x
y p '' = 80 e 4x
99
y p ''' = 320 e 4x
y p ( IV ) = 1280 e 4x
+ 6 ( 320 e ) + 13 ( 80 e )
enseguida se sustituye en ( A )
1280 e 4x 4x 4x
= Q ( x)
4240 e 4x
= Q ( x)
RESOLUCIÓN
considera a la función
Q ( x) = x e 2x
100
cuyo operador anulador es
P ( D ) = ( D − 2) 2
(D − 2) 2
(D 3
)
+ 2 D 2 − 4 D − 8 y = ( D − 2) 2
Q ( x)
(D − 2) 2
(D 3
+ 2D2 − 4D − 8 y = 0 )
( )
su ecuación característica y raíces respectivas son
(λ − 2) 2 λ 3 + 2 λ 2 − 4 λ − 8 = 0 ;
λ1 = 2 , λ2 = λ3 = − 2 , λ4 = λ5 = 2
Por lo que su solución general es
yG = C1e 2x
+ C 2 e −2x + C 3 x e −2x + C 4 x e 2x + C 5 x 2 e 2x
y p = C 4 x e 2x + C 5 x 2 e 2x
yp ' = 2 C 4 x e 2x
+ C 4 e 2x + 2C 5 x 2 e 2x + 2C 5 x e 2x
y p '' = 4 C 4 x e 2x
+ 4C 4 e 2x + 4C 5 x 2 e 2x + 8C 5 x e 2x + 2C 5 e 2x
y p ''' = 8 C 4 x e 2x
+ 12 C 4 e 2 x + 8 C 5 x 2 e 2 x + 24 C 5 x e 2 x + 12 C 5 e 2 x
(16 C 4 + 16 C 5 ) e 2x
+ 32 C 5 x e 2 x = x e 2 x
C4 = − C5
resolviendo para las constantes
101
C4 = −
1
32
32 C 5 = 1 ⇒ C5 =
1
32
sustituyendo estos valores en la solución particular se tiene que
= − x e 2x +
1 1 2 2x
y p x e
32 32
Y la solución general es finalmente
(x )
14) Obtenga la solución particular de la ecuación diferencial
2
− 1 y '' − 2 x y ' + 2 y = x 2 − 1
( )
si la solución de la ecuación homogénea asociada correspondiente es
yh = C1 x + C2 1 + x2
RESOLUCIÓN
obtiene
y '' − y' + 2 y = 1
2x 2
x − 1
2
x − 1
( )
La forma de y es
= u ( x) x + v ( x) 1 + x 2
p
sistema
⎡x 1 + x 2 ⎤ ⎡u ' ( x ) ⎤ ⎡0 ⎤
⎥ = ⎢ ⎥
2 x ⎦ ⎣v '( x)⎦
⎢ ⎥⎢
⎣1 ⎣1⎦
102
Para resolver por Cramer es necesario obtener los siguientes determinantes:
⎡x 1 + x 2⎤
∆ = ⎢ ⎥ = x − 1
⎣1 2x ⎦
2
( )
⎡0 1 + x 2⎤
∆1 = ⎢ ⎥ = − x + 1
⎣1 2x ⎦
2
⎡x 0⎤
∆2 = ⎢ = x
⎣1 1 ⎥⎦
( )
de lo anterior, resulta
∆1 − x2 + 1
u '( x) = = = −1 −
(x − 1)( x + 1)
2
∆ x 2
− 1
∆2
v '( x) = =
x
∆ x − 1
2
Integrando cada función u ' ( x ) y v ' ( x ) obtenemos las funciones deseadas para yP
∫ ∫ ∫
⎛ ⎞ ⎛ ⎞
u ( x) = ⎜⎜ −1 − ⎟ dx = − dx − 2 ⎜⎜ ⎟⎟
( x − 1)( x + 1) ⎟⎠ ⎝ ( x − 1)( x + 1) ⎠
2 dx
⎝
∫
⎛ ⎞
⎜⎜ ⎟⎟ se emplea la descomposición en fracciones parciales,
( )( )
dx
⎝ − + ⎠
para obtener
x 1 x 1
esto es,
= +
(x − 1)( x + 1)
1 A B
x − 1 x + 1
1 = A ( x − 1) + B ( x − 1)
de aquí resulta A = B = −
1 1
,
2 2
∫ ∫ ∫ x + 1 = 2 Ln ( x − 1) − 2 Ln ( x + 1)
por lo que
⎛ ⎞
⎜⎜ ⎟⎟ d x = −
⎝ ( x − 1)( x + 1) ⎠
1 1 dx 1 dx 1 1
2 x − 1 2
= − ⎡ Ln ( x + 1) − Ln ( x − 1) ⎤⎦
2⎣
1
⎛ x + 1⎞
= − Ln ⎜ ⎟
1
2 ⎝ x − 1⎠
103
Regresando a u ( x ) resulta
⎛ x + 1⎞
u ( x ) = − x + Ln ⎜ ⎟
⎝ x − 1⎠
De manera similar para v ( x )
v ( x) =
∫ dx = Ln x 2 − 1 ( )
x 1
x − 1
2
2
Finalmente, una solución particular es
( )( )
⎡ ⎛ x + 1 ⎞⎤ ⎡1 ⎤
yP = ⎢ − x + Ln ⎜ ⎟ ⎥ x + ⎢ Ln x − 1 ⎥ 1 + x
⎣ ⎝ − ⎠⎦ ⎣ ⎦
2 2
x 1 2
⎛ x + 1⎞
yP = − x 2 + x Ln ⎜ ⎟ +
1
(
Ln 1 + x 2 Ln x 2 − 1 ) ( )
⎝ x − 1⎠ 2
15) {
Si 1 , x − 1
} es un conjunto fundamental de soluciones de la ecuación diferencial
x 2 y '' + 2 x y ' = 0 , obtenga la solución general de la ecuación
x 2 y '' + 2 x y ' = x − 1
RESOLUCIÓN
104
Por otro lado, la solución y p se obtiene empleando el método de variación de parámetros,
esto es
y = u ( x )( 1) + v ( x ) x − 1
⎡1 x − 1 ⎤ ⎡u ' ( x ) ⎤ ⎡ 0 ⎤
⎥ = ⎢
⎣0 − x ⎦ ⎣ v ' ( x ) ⎦ ⎣Q ( x ) ⎦
⎢ − 2⎥ ⎢ ⎥
y '' + y' = 3
2 1
x x
de esta forma se identifica a la función deseada
Q ( x) =
1
x3
Así, el sistema anterior se escribe nuevamente
⎡ 1 x − 1 ⎤ ⎡u ' ( x ) ⎤ ⎡ 0 ⎤
⎥ = ⎢ − 3⎥
⎣0 − x ⎦ ⎣ v ' ( x ) ⎦
⎢ − 2⎥ ⎢
⎣x ⎦
que se resuelve utilizando la regla de Cramer. Para resolverlo se obtiene
∆1 ∆
u '( x) = , v '( x) = 2
∆ ∆
donde
⎡1 x − 1 ⎤
∆ = ⎢ − 2⎥
= − x−
⎣0 − x ⎦
2
⎡ 0 x− ⎤
∆1 = ⎢ − 3 − 2⎥
= −x−
1
⎣x −x ⎦
4
105
⎡1 0 ⎤
∆2 = ⎢ − 3⎥
= x−
⎣0 x ⎦
3
Entonces
− x−
u '( x) = ; v '( x) =
x− 3
= x− = − x−
4
− x− − x − 2
2 1
2
∫x
Integrando
u ( x) = − 2
dx = −
1
∫
x
v ( x) = −
⎛1⎞
dx = − Ln x = Ln ⎜ ⎟
dx
x ⎝ x⎠
Enseguida se sustituye en y p , esto es
= ⎜ − ⎟ ( 1) + Ln ⎜ ⎟ ⎜ ⎟
⎛ 1⎞ ⎛1⎞⎛1⎞
⎝ x⎠ ⎝ x⎠⎝ x⎠
y p
⎡ ⎛1⎞ ⎤
= ⎢ Ln ⎜ x ⎟ − 1⎥
1
⎣ ⎝ ⎠ ⎦
y p
x
Finalmente, la solución general es
⎡ ⎛1⎞ ⎤
y = C1 + C2 x− 1 + ⎢ Ln ⎜ x ⎟ − 1⎥
1
x ⎣ ⎝ ⎠ ⎦
16) Considere una ecuación diferencial lineal no homogénea con coeficientes constantes,
RESOLUCIÓN
106
Se requiere conocer primero yh (solución de la ecuación homogénea asociada) lo cual
y h = C 1 e−x + C 2 xe−x
= u ( x) e − x + v ( x) x e − x
p
⎡ 0 ⎤
⎡ e− x xe − ⎤ ⎡u ' ( x ) ⎤ ⎢ −x ⎥
⎥ = ⎢ e
+ e ⎦ ⎣ ( )⎦
− x⎥ ⎢
⎢ ⎥
x
⎣− e − xe−
− x
⎢⎣ 1 + x 2 ⎥⎦
x
v ' x
(− )+ (x e )
⎡ e− x xe − ⎤
∆ = ⎢ − x⎥
= e− xe− + e− e− − x
x
⎣− e − xe + e ⎦
− x − x
x x x x
∆ = − xe− x
+ e− 2x
+ xe− 2x
= e− 2x
⎡ 0 ⎤
( )
xe −
⎢ ⎥ ⎛ e− x ⎞
x
∆1 = ⎢ e − x − x⎥
= − ⎜ 2 ⎟
xe−
− xe + e ⎥ ⎝1 + x ⎠
− x
x
⎢1 + x 2
⎣ ⎦
x e − 2x
∆1 = −
1 + x2
⎡ e− x ⎤
⎢ ⎥ e − 2x
0
∆2 = ⎢ − x
⎥ = 1 + x2
⎢− e
− x
2⎥
e
⎣ 1 + x ⎦
de lo anterior se tiene
x e − 2x
−
1 + x2
u '( x) = = −
x
e − 2x
1 + x2
107
u ( x) = −
∫ dx = − Ln 1 + x 2( )
x 1
1 + x 2
2
⎛ ⎞
u ( x ) = Ln ⎜ ⎟
⎜ 1 + x2 ⎟
1
⎝ ⎠
e − 2x
1 + x2
v '( x) = =
1
− 2x
1 + x2
∫
e
v ( x) = = ang tan x
dx
1 + x2
Por lo que finalmente
⎛ ⎞
= Ln ⎜ ⎟e− + ( ang tan x ) x e −
⎜ 1 + x ⎟
1 x x
⎝ ⎠
y p 2
RESOLUCIÓN
(y h),
cuyos coeficientes están determinados; lo anterior indica que la solución general involucra la
homogénea ( y ) , es decir
solución de la ecuación homogénea asociada más una solución particular de la no
yG = yh + y p
108
Para y h = C 1 e x + C 2 cos x + C 3 sen x es posible identificar las raíces que la generan,
las cuales son
λ 1 = 1 , λ 2,3 = ± i
( )
y en términos de factores se tiene
(λ − 1) λ 2 + 1 = 0
( )
ecuación diferencial en términos del operador diferencial
(D − 1) D 2 + 1 = 0
(D 3
− D2 + D − 1 = 0 )
o bien
y ''' − y '' + y ' − y = 0
Al realizar la sustitución en ( A ):
x sen x + x cos x − 3 cos x + 3 sen x + x cos x
+ 2 sen x + 2 cos x − x sen x − x sen x + cos x
Q ( x ) = 4 sen x
reduciendo términos se obtiene
109
EJERCICIOS PROPUESTOS
( )
diferencial:
( ) ( )
Solución.
y p = C 3 x + C 4 x 2 + C 5 x 3 sen 3 x + C 6 x + C 7 x 2 + C 8 x 3 cos 3 x
x y '' + y ' − y = x + x3
4
x
utilizando el método de variación de parámetros si:
yh ( x ) = C1 x 2 + C2 x − 2
Solución.
yG = C1x2 + C2x− + x Ln x + x −
2 1 2 1 4 x2
4 12 16
( )
Solución.
yG = C 1 e − + C 2e− − 2x + x2 e−
2x 3x 1 3x
2
110
4) Resuelva la ecuación diferencial
y '' + y = sec x
Solución.
y '' − 4 y ' + 4 y =
2e 2 x
x3
P ( D ) = Q ( x ) …………………………..( 1 )
6) Sea la ecuación diferencial lineal
Si { }
A = x − 2 e 3 x , 4 + 5 e 3 x , 2 − 9 x ¨ es un conjunto de algunas soluciones
Q ( x)
Solución.
a) y G = C 1 + C 2 x + C 3 e 3 x − e 4 x
b) D 2( D − 3 ) y = − 16 e 4 x
111
7) Utilice el método de coeficientes indeterminados para obtener la forma de una
(D )
solución particular de la ecuación diferencial
2
− 2 D − 3 y = x sen 2 x + x 3 e 3 x
( )
Solución.
y p = C 3 x + C 4 x 2 + C 5 x 3 + C 6 x 4 e 3 x + ( C 7 + C 9 x ) cos 2 x + ( C 8 + C 10 x ) sen 2 x
L 1 = D 2 + 1 y L 2 = ( sen x ) D
Obtenga
a) (L 1 L 2) y ;
b) (L 2 L 1) y ;
c) de acuerdo al resultado obtenido en los incisos a) y b) ¿qué puede concluir
respecto al producto de los operadores L1 y L2 ?
Solución.
a) (L 1 L 2) y = 2 ( cos x ) y '' + ( sen x ) y '''
y = C 1 e x + C 2 e−x − +
1 1
cos 2 x
2 10
Solución.
y '' − y = − cos 2 x
1 1
2 2
112
10) Obtenga la solución general de la ecuación diferencial
y' = 2 + ex − y
( )
Solución.
y = C e−x + e 2 + ex
2 −x
3
2
3
y ''' + y ' = 6 x ( 2 x + 1)
11) Determine la ecuación diferencial
Solución.
y = C 1 + C 2 x + C 3 e−x + 9 x 2 − 3 x 3 + x 4
y '' − 4 y ' + 4 y =
2e 2x
x
Solución.
y = ( C 1 + C 2 x ) e 2 x − 2 x e 2 x ( 1 − Ln x )
Solución.
y = C e sen x + 2 sen x − 2
113
14) Resuelva la ecuación diferencial
(D + 1)( D − 1) y = cos
x
2
Solución.
y = C1e + C 2 e−x −
x 4 x
cos
5 2
Solución.
AB = ( x D − 1)( D + 1) = x D 2 + x D − D − 1
B A = ( D + 1)( x D − 1) = D( x D) − D − x D − 1
Se demuestra que AB ≠ B A .
114
TEMA 4
⎡0 1 0 ⎤
C = ⎢⎢0 0 1 ⎥⎥
⎢⎣ 1 − 1 1 ⎥⎦
RESOLUCIÓN
(λ i ) 2
0 1 2
= β + β1λi + β
λi t
e 0 2 ………………….( B )
et = β 0 + β 1 + β 2 ………………………...….( C )
λ2 = i
115
eti = β 0 + β 1 i + β 2 i 2 ……………………… ( D )
eti = β 0 + β1i − β 2
β 1 = sen t
De la ecuación ( C )
β = et − β − β
β = e − sen t − β
0 1 2
t
………………….( G )
0 2
además ( F ) = ( G )
cos t + β 2 = e t − sen t − β 2
e t − sen t − cos t
β 2 =
2
( )
sustituyendo en ( G ) se tiene
⎡ β0 β β ⎤
= ⎢⎢ β2 β − β β 1 + β 2 ⎥⎥
1 2
At
e
⎢⎣ β 1 + β − β β 0 + β 1 ⎥⎦
0 2
2 1
finalmente
116
⎡ e t + cos t − sen t e t − cos t − sen t ⎤
1⎢ ⎥
2 sen t
= ⎢e t − cos t − sen t e t − cos t + sen t ⎥
2⎢ t ⎥
e At 2 cos t
⎡ 0 1 0⎤
2) Determine e At
para A = ⎢⎢ − 1 0 0 ⎥⎥
⎢⎣ 0 0 1 ⎥⎦
RESOLUCIÓN
(λ )
además, es necesario emplear la expresión
= β + β 1 λi + β
λi t 2
e 0 2 i ………………….( B )
( )
0
(1 − λ) λ 2 + 1 , λ1 = 1 , λ2 = i , λ3 = − i
et = β 0 + β 1 + β 2 ……………….……… ( C )
para λ = i
117
β 0 + β1i − β 2 = cos t + i sen t
de donde resulta
β − β = cos t
β 1 = sen t
0 2
…………………..( E )
( )
De la ecuación ( C ) y ( E ) se obtiene
β = e − sen t + cos t
1 t
( )
0
2
β = e − sen t − cos t
1 t
2
2
Ahora bien, considerando la ecuación ( A ) es conveniente expresarla en términos de β 0 ,
⎡1 0⎤ ⎡ 0 0⎤ ⎡− 1 0⎤
⎢0 ⎥ ⎢ 0 ⎥⎥ + β ⎢ 0 0 ⎥⎥
0 1 0
e At
= β ⎢ 1 0 ⎥ + β 1 ⎢− 1 0 ⎢ − 1
⎢⎣ 0 1 ⎦⎥ ⎣⎢ 0 1 ⎦⎥ ⎢⎣ 0 1 ⎦⎥
0 2
0 0 0
⎡β 0 − β 2 β ⎤
= ⎢⎢ ⎥
0
− β1 β − β
1
⎥
At
e 0
⎢⎣ β + β 1 + β 2 ⎥⎦
1 2
0 0 0
⎡ cos t 0⎤
⎢ ⎥
sen t
At
= ⎢ − sen t 0⎥
⎢ 0 e t ⎥⎦
e cos t
⎣ 0
⎡ cos t − sen t ⎤
e At = e 2t ⎢
⎣ sen t cos t ⎥⎦
Obtenga
a) las funciones β 0 y β 1 tales que e A t = β 0 I + β 1 A ,
118
RESOLUCIÓN
= β + β 1 λi
λi t
e 0
Para λ 1 = 2 + i
e(
2 + i) t
= β 0 + β 1 (2 + i)
e 2t ⋅ e t i = β 0 + 2β 1 + β 1 i ………...………...( A )
= e 2 t cos t + i e 2 t sen t
y por igualdad de números complejos
β 0 + 2 β 1 = e 2 t cos t
β 1 = e 2 t sen t
por lo que resulta
β 0 = e 2 t cos t − 2 e 2 t sen t
119
⎡β + aβ bβ1 ⎤
e At = ⎢
⎣ cβ1 β 0 + d β 1 ⎥⎦
0 1
……………………..( C )
c os t + ( sen t )( a − 2 ) = cos t
( sen t ) b = − sen t
( sen t ) c = sen t
c os t + ( sen t )( d − 2 ) = cos t
120
4.2 TRANSFORMACIÓN DE UNA ECUACIÓN DIFERENCIAL DE
ORDEN n A UN SISTEMA DE n ECUACIONES DE PRIMER
ORDEN.
y '' − 2 y ' + y = 1
1) Sea la ecuación diferencial
RESOLUCIÓN
y '' − 2 y ' + y = 1
Se tiene
y ( 0) = 1 y y '( 0) = 1
a) Para realizar la transformación se introducen tantas variables nuevas como el orden de la
ecuación diferencial; para este caso se tendrán 2 variables, a saber
y = w1 ……………………………....( 1 )
w 2' = − w 1 + 2 w 2 + 1
y en forma matricial
⎡ w 1'⎤ ⎡ 0 1⎤ ⎡ w 1 ⎤ ⎡0⎤
⎢w '⎥ = ⎢− 1 2⎥ ⎢w ⎥ + ⎢1⎥
⎣ 2⎦ ⎣ ⎦ ⎣ 2⎦ ⎣ ⎦
que es el sistema de ecuaciones diferenciales equivalente.
121
En lo que se refiera a las condiciones iniciales dadas, para obtener el vector w( 0) , se
y = w 1 , por lo que si y ( 0 ) = 1 ⇒ w 1 ( 0 ) = 1
tiene
entonces
⎡1⎤ ⎡0⎤
w( 0) = ⎢ ⎥ y b ( t ) = ⎢ ⎥ que es el vector de términos independientes
⎣1⎦ ⎣1⎦
w 1' = A w + b ( t )
b) El sistema es de la forma
∫
que es un sistema de primer orden; su solución está dada por
w = e At w (0) + e A τ b ( t − τ) d τ
t
donde
⎡ 0 1⎤
A=⎢ ⎥
⎣− 1 2⎦
y las funciones β 0 y β 1 se obtendrán considerando los valores característicos; para
determinar éstos últimos
⎡0 − λ 1 ⎤
det A − λ I = ⎢ ⎥ = ( λ − 1) 2
⎣ − 1 2 − λ⎦
entonces
λ1 = λ2 = 1
= β + β 1 λ i …………………..( 2 )
λi t
e 0
pero dado que los valores característicos son repetidos, deberá derivarse ( 2 ) con
respecto a λ ,
122
= ( β 0 + β 1 λ)
d (λ t )
dλ dλ
d
e
t eλt = β 1
y sustituyendo λ = 1 se tiene
β = t et
et = β0 + t et ⇒ β0 = e t − t e t
así, la matriz exponencial es
⎡1 0 ⎤ ⎡ 0 1⎤
e At = β ⎢0 1 ⎥ + β 1 ⎢ − 1 2 ⎥
⎣ ⎦ ⎣ ⎦
0
⎡ β β ⎤
e At = ⎢ 0
⎣− β 1 β + 2 β 1 ⎥⎦
1
⎡e t − t e t ⎤
= ⎢ ⎥
t et
⎣⎢ − t e + t e t ⎦⎥
At
e t
et
wp =
∫ e A t b ( t − τ) d τ
t
por lo que
⎡e t − t e t t e t ⎤ ⎡1⎤ ⎡e t ⎤
wh = ⎢ ⎥⎢ ⎥ = ⎢ ⎥
⎢⎣ − t e e t + t e t ⎥⎦ ⎣1⎦ ⎢⎣e t ⎥⎦
t
y para w p
⎡e τ − τ e τ τ e τ ⎤ ⎡0⎤
b ( t − τ) = ⎢ ⎥⎢ ⎥
⎣⎢ − τ e e τ + τ e τ ⎦⎥ ⎣1⎦
At
τ
e
123
∫ ∫
enseguida se procede a integrar
⎡ τeτ ⎤ ⎡I 1 ⎤
wp = e A t b ( t − τ) d τ = ⎢ τ ⎥ τ = ⎢I ⎥
t t
⎣⎢e + τ e ⎦⎥ ⎣ 2⎦
τ
d
0 0
∫
donde
I1 = τ e τ dτ = τ e τ − e τ = t et − et + 1
t t
∫
0
(e + τ e ) dτ = e
0
I2 = τ τ
+ τeτ − eτ = eτ + t et − et = t et
t t
0 0
entonces
⎡t e t − e t + 1⎤
= ⎢ ⎥
⎣⎢ ⎥⎦
wp
t et
⎡ t et + 1 ⎤
w= ⎢ t ⎥
⎣⎢e + t e ⎦⎥
t
124
4.3 SISTEMAS HOMOGÉNEOS DE PRIMER ORDEN.
x 1' = − 2 x 1 − 4 x 2
x2 ' = x1 − 2 x2
RESOLUCIÓN
e At = β 0 I + β 1 A ………………….. ( 1 )
donde
λ i , i = 1 , 2 , ... , n son los valores característicos de la matriz de coeficientes y
n es el orden de la matriz.
Los λ i se obtienen del determinante A − λI , esto es
125
⎡− 2 − λ − 4 ⎤
det A − λ I = ⎢ ⎥ = λ 2 + 4λ + 8
⎣ 1 − 2 − λ⎦
e(
− 2 + 2 i) t
= β 0 + β 1 (− 2 + 2 i)
e− 2t
e2 t i = β 0 − 2β 1 + 2 β 1 i …………………….( 3 )
Considerando la representación en forma de Euler para números complejos
e θ i = cos θ + i sen θ
e 2 t i = cos 2 t + i sen 2 t
2 β 1 = e − 2 t sen 2 t …………………….( 6 )
De la ecuación ( 6 )
β1 =
1 − 2t
e sen 2 t
2
Al sustituir β 1 en ( 5 ) y despejar β 0
⎛1 ⎞
β 0 = 2 ⎜ e − 2 t sen 2 t ⎟ + e − 2 t cos 2 t
⎝2 ⎠
β 0 = e − 2 t sen 2 t + e − 2 t cos 2 t
126
Según se observa, se han determinado los valores de β0 y β1 a partir del valor
e At = β 0 I + β 1 A
⎡1 0⎤ ⎡− 2 − 4⎤
e At = β 0 ⎢ ⎥ + β ⎢ 1 − 2⎥
⎣0 1⎦ ⎣ ⎦
1
⎡ − 2 t cos 2 t − 2 e − 2 t sen 2 t ⎤
⎡β 0 − 2 β 1 − 4β1 ⎤ ⎢ e ⎥
=⎢ ⎥ = 1 − 2t
⎣ β1 β 0 − 2 β 1⎦ ⎢ e cos 2 t ⎥
At
−
e
⎢⎣ 2 ⎥⎦
2 t
e sen 2 t
⎡ cos 2 t − 2 sen 2 t ⎤
x =e − 2t ⎢1 ⎥⎡ 1 ⎤
⎢ sen 2 t ⎢ ⎥
cos 2 t ⎥ ⎣− 1⎦
⎣2 ⎦
Finalmente
⎡ e − 2 t ( cos 2 t + 2 sen 2 t ) ⎤
⎢ ⎥
x = ⎢ − 2t ⎛ 1 ⎞⎥
e ⎜ sen 2 t − cos 2 t ⎟
⎢⎣ ⎝2 ⎠⎥⎦
x 2' = 3 x 1 + 7 x 2
utilizando el método de la matriz exponencial.
RESOLUCIÓN
127
x' = A x
y su solución está dada por
x = e A t x ( 0) ……………… ( 1 )
⎡k1 ⎤
en este caso x ( 0) = ⎢ ⎥
⎣k2 ⎦
La matriz exponencial se obtiene de
e A t = β 0 + β 1 A ………………..( 2 )
⎡1 − λ − 3 ⎤
det A − λ I = ⎢ ⎥ = λ 2 − 8 λ + 16
⎣ 3 7 − λ⎦
Así, la ecuación característica es
λ 2 − 8 λ + 16 = 0
y sus valores característicos son λ 1 = λ 2 = 4
Cada valor característico se sustituye en
= β + β 1 λ i ………………( 3 )
λi t
e 0
para λ 1 = 4
e A t = β 0 + 4 β 1 …………………….( 4 )
ahora bien, como se tienen valores característicos repetidos, la ecuación ( 3 ) debe
derivarse con respecto a λ , esto es
d λt
dλ
e ( )
=
dλ
d
( β 0 + β 1 λ)
t eλt = β 1
t e4t = β 1
y de la ecuación ( 4 )
128
β 0 = e 4 t − 4β 1
β 0 = e 4t − 4t e 4t
⎡1 0 ⎤ ⎡1 − 3⎤
e At = β 0 ⎢ ⎥ + β1⎢ ⎥
⎣0 1⎦ ⎣3 7 ⎦
⎡β + β 1 − 3β 1 ⎤
e At = ⎢ 0
⎣ 3β 1 β 0 + 7 β 1 ⎥⎦
Sustituyendo β 0 y β 1
⎡e 4 t − 3t e 4 t − 3t e 4 t ⎤
e At = ⎢ ⎥
⎢⎣ 3t e
4t
e 4 t + 3t e 4 t ⎥⎦
Finalmente, la solución general es:
x = e A t x ( 0)
⎡e 4 t − 3t e 4 t − 3t e 4 t ⎤ ⎡k 1 ⎤
x = ⎢ ⎥⎢ ⎥
⎣⎢ 3t e
4t
e 4t + 3t e 4 t ⎦⎥ ⎣k 2 ⎦
x 1 = e 4 t ⎡⎣k 1 − t ( 3 k 1 + 3 k 2 ) ⎤⎦
x 2 = e 4 t ⎡⎣k 2 + t ( 3 k 1 + 3 k 2 ) ⎤⎦
y si consideramos
C 1 = k1 , C 2 = k 2 , C 3 = 3 k 1 + 3 k 2
entonces
x 1 = e 4 t(C 1 − C 3 t )
x 2 = e 4 t(C 2 + C 3 t )
129
3) Escoger la opción que contenga la solución correcta del sistema de ecuaciones
⎧ x' = 3 x − 2 y
⎨ y determine la solución particular para t = 0 si x ( 0) = 1
⎩ y' = x + y
y y ( 0) = 2 .
⎧⎪ x ( t ) = e t ( C 1 cos 2t + C 2 sen 2t )
A) ⎨
⎪⎩ y ( t ) = e ⎡⎣( C 1 − C 2 ) cos 2t − ( C 1 − C 2 ) sen 2t ⎤⎦
t
⎧⎪ x ( t ) = C 1 e 2 t + C 2 e t
B) ⎨
⎪⎩ y ( t ) = − 3 C 1 e + 2 C 2 e
t t
⎪⎧ x ( t ) = C 1 e cost + C 1 ( C 1 − 2 C 2 ) e sent
C) ⎨
⎪⎩ y ( t ) = C 2 e cost + ( C 1 − C 2 ) e sent
2t 2t
2t 2t
RESOLUCIÓN
3−λ −2
det A − λ I = = ( 3 − λ ) (1 − λ ) + 2
1 1−λ
λ 2 − 4λ + 5 = 0 ⇒ λ1 = 2 + i , λ 2 = 2 − i
Estas raíces permiten obtener términos como e 2 t cos t y e 2 t sen t , por lo que la
solución del sistema los debe contener.
Considerando las diferentes opciones, la única posible es la del inciso ( C ) y al aplicar
las condiciones iniciales dadas se llega a
C1 = 1 , C2 = 2
por lo que la solución particular resulta
x ( t ) = e 2 t cos t − 3 e 2 t sen t
y ( t ) = 2 e 2 t cos t − e 2 t sen t
130
4) Si el sistema de ecuaciones diferenciales
x '1 = x1 + x2 + x3
x '2 = − 2 x1 − 2 x2 − 2 x3
x '3 = x1 + x2 + x3
tiene como matriz exponencial
⎡t + 1 t ⎤
= ⎢⎢ − 2 t − 2 t ⎥⎥
t
e At
− 2t + 1
⎢⎣ t t t + 1⎥⎦
β k (t ) ∀ k = 0 , 1 , 2 tal que e
λjt
= ∑ (λ β (t ))
2
k
obtener las funciones
k = 0
i k
RESOLUCIÓN
⎡0 0 0⎤
A = ⎢⎢0 0 0⎥⎥
2
⎢⎣0 0 0⎥⎦
y se obtiene la matriz
= β + β 1λ + β 2 λ 2j
λjt
e 0 j …………………….( A )
131
= ( 1 − λ ) ⎡⎣( − 2 − λ )( 1 − λ ) + 2⎤⎦
⎣ 2( 1 − λ ) + 2⎤⎦ + ⎡−
= − 1 ⎡− ⎣ 2 − ( − 2 − λ )⎤⎦
= λ3
d 2 λt
dλ
=
d ⎛ d λt ⎞
⎜
dλ ⎝ dλ
⎟ =
d2
(
β0 + β1λ + β2λ2 )
⎠ dλ
2
e e 2
dλ 2
e =
d 2 λt
dλ
d
(
t eλt = )
dλ
d
( β 1 + 2 β 2 λ ) ………..….( B )
Finalmente t 2 e λ t = 2 β 2 …………………..( C )
Enseguida se sustituye en ( A ) , ( B ) y ( C ) el valor λ = 0 y se obtiene
β0 = 1 , β1 = t y β2 =
t2
que son las funciones pedidas.
2
132
4.4 SISTEMAS NO HOMOGÉNEOS DE PRIMER ORDEN.
⎡2⎤
sujeto a las condiciones x ( 1) = ⎢ ⎥ , si la matriz exponencial correspondiente es
⎣0⎦
⎡( 1 + t ) e 2 t ⎤
=⎢ ⎥
(1 − t ) e 2 t ⎥⎦
t e 2t
⎢⎣ − t e
At
e 2t
⎡1⎤
y el vector de términos independientes es b ( t ) = ⎢ ⎥
⎣0⎦
RESOLUCIÓN
∫
orden no homogéneo, cuya solución es de la forma
x(t 0 ) + e A τ b ( t − τ) d τ
t − t0
A ( t − t0 )
x =e
xh = e ( x(t 0 )
A t − t0 )
xh = e ( x ( 1)
A t − 1)
2 t − 1 ⎡1 + ( t − 1) t −1 ⎤ ⎡2⎤
xh = e ( ) ⎢
⎣ − ( t − 1) 1 − ( t − 1) ⎦ ⎢⎣0⎥⎦
⎥
⎡ t t − 1⎤ ⎡2⎤
xh = e (
− 1)
⎢− t + 1 2 − t ⎥⎦ ⎢⎣0⎥⎦
⎣
2 t
⎡ 2t ⎤
xh = e (
− 1)
⎢2 − 2 t ⎥
⎣ ⎦
2 t
133
∫
En lo que respecta a la solución particular del sistema no homogéneo, está dada por
e Aτ b ( t − τ) d τ
t − t0
xp =
0
En esta expresión conviene trabajar el integrando de manera que tome una forma más
sencilla para integrar.
Realizando el producto
⎡( 1 + τ) e 2 τ ⎤
Aτ
b ( t − τ) d τ = ⎢ ⎥
⎣⎢ − τ e ⎦⎥
2τ
e
∫ ∫ ∫
el primer elemento se tiene
( 1 + τ) e 2 τ d τ =
t − 1 t − 1 t − 1
e 2τ d τ + τ e 2τ d τ
0 0 0
⎡1 ⎤ ⎡1 ⎤
t − 1 t − 1
= ⎢ e 2τ ⎥ + ⎢ τ e 2τ − e 2 τ ⎥
1
⎣2 ⎦ 0 ⎣2 ⎦0
( )
4
− e 0 + ⎢ ( t − 1) e (
1 2 (t − 1) ⎡1 − 1) 1 2 (t − 1) 1⎤
= − +
⎣2 4 ⎥⎦
2 t
e e
2 4
1 2 (t − 1) 1 1 2 (t − 1) 1 2 (t − 1) 1 2 (t − 1)
= − + te − − +
1
e e e
2 2 2 2 4 4
1 2 (t − 1) 1 2 (t − 1)
= − −
1
te e
2 4 4
∫
Para el siguiente elemento
⎡1 1 ⎤
t − 1 t − 1
− τ e 2τ d τ = − ⎢ τ e 2τ − e 2τ ⎥
0 ⎣2 4 ⎦0
= − ⎢ ( t − 1) e (
⎡1 − 1) 1 2 (τ − 1) 1⎤
− +
⎣2 4 ⎥⎦
2 t
e
4
⎡1 2 t 1⎤
= − ⎢ te (
− 1) 1 2 (t − 1) 1 2 (τ − 1)
− − +
⎣2 4 ⎥⎦
e e
2 4
1 2 (t − 1) 3 2 (τ − 1)
=− + +
1
te e
2 4 4
Entonces la integral
134
∫ e Aτ b ( t − τ ) d τ
t − 1
da por resultado
⎡ t 2 (t − 1) 1 2 (t − 1) 1 ⎤
⎢ 2 − −
4 ⎥
⎢ ⎥
e e
4
⎢− t e 2 (t − 1) − 3 e 2 ( t − 1) − 1 ⎥
⎢⎣ 2 4 4 ⎥⎦
Por lo que
⎡ ⎛t 1 ⎞ 2 (t − 1) 1 ⎤
⎢ ⎜2 − 4⎟e −
4 ⎥⎥
⎝ ⎠
=⎢
⎢⎛ t 3 ⎞ 2 (t − 1) 1⎥
xp
⎢⎜ − + ⎟ e − ⎥
⎣⎝ 2 4 ⎠ 4⎦
Finalmente, la forma de la solución del sistema es
x = xh + x p
1 ⎡ 3e + 2 3 e 5t − 3 ⎤
e At = ⎢ ⎥
5t
5 ⎢⎣ 2 e 5t − 2 2 e 5t + 3 ⎥⎦
RESOLUCIÓN
135
Se tiene un sistema no homogéneo de primer orden cuya representación matricial es
x ' = A x + b ( t ) , donde la matriz de coeficientes es
⎡3 3⎤
A=⎢ ⎥
⎣2 2⎦
y el vector de términos independientes es
⎡ 5⎤
b (t ) = ⎢ ⎥
⎣0⎦
Asimismo, el vector de condiciones iniciales es
⎡10⎤
x ( 0) = ⎢ ⎥
⎣5⎦
∫
forma
x = e At x ( 0) + e A τ b ( t − τ) d τ
t
∫
y una solución particular del sistema no homogéneo asociado
xp = e Aτ b ( t − τ) d τ
t
1 ⎡3 e + 2 3 e 5 t − 3 ⎤ ⎡10⎤
xh = e x ( 0) = ⎢ 5 t ⎥⎢ ⎥
5t
5 ⎣⎢ 2 e − 2 2 e 5 t + 3⎦⎥ ⎣ 5 ⎦
At
efectuando operaciones
1 ⎡45 e + 5⎤ ⎡9 e + 1⎤
xh = ⎢ =
⎥ ⎢ ⎥
5t 5t
En lo que respecta a x p , se tiene que calcular la integral de una matriz, que corresponde
136
1 ⎡3 e + 2 3 e 5 τ − 3 ⎤ ⎡5⎤
b ( t − τ) d τ = ⎢ 5 τ
5τ
⎥⎢ ⎥
Aτ
5 ⎢⎣ 2 e − 2 2 e 5 τ + 3⎥⎦ ⎣0 ⎦
e
⎡3 e 5 τ + 2⎤
= ⎢ 5τ ⎥
⎢⎣2 e − 2⎥⎦
∫ ∫
Enseguida se integra
⎡3 e 5 τ + 2⎤
e A t b ( t − τ) d τ = ⎢ 5τ ⎥ dτ
t t
0 0 ⎣⎢ 2 e − 2⎥⎦
∫
Para cada uno de los términos se tiene
(3 e 5τ
) ⎡3 ⎤
+ 2 d τ = ⎢ e 5 τ + 2 τ⎥ = e 5 t − + 2 t
t t
3 3
⎣5 ⎦0
∫
0 5 5
( ) ⎡2 ⎤
2 e 5 τ − 2 d τ = ⎢ e 5 τ − 2 τ⎥ = e 5 t − − 2 t
t t
2 2
0 ⎣5 ⎦0 5 5
Por lo que
⎡3 5t 3⎤
⎢5 e + 2t − 5⎥
= ⎢ ⎥
⎢ 2 e 5 t − 2t − 2⎥
xp
⎣⎢ 5 5 ⎥⎦
Finalmente
x = xh + x p
⎡3 5t 3⎤
⎡9 e 5 t + 1⎤ ⎢ 5 e + 2 t − 5 ⎥
x = ⎢ 5t ⎥ + ⎢ ⎥
⎣⎢6 e − 1⎥⎦ ⎢ 2 e 5 t − 2 t − 2 ⎥
⎢⎣ 5 5 ⎥⎦
⎡ 48 5 t 2⎤
⎢ 5 e + 2t + 5⎥
x = ⎢ ⎥
⎢ 32 e 5 t − 2 t − 7 ⎥
⎢⎣ 5 5 ⎥⎦
3) Resuelva el sistema
x' = x − y
y' = x + 3 y − 1
137
cuya solución está sujeta a x ( 0) = 1 , y ( 0) = 0
RESOLUCIÓN
∫
Para un sistema no homogéneo, su solución es de la forma
x = e A t x ( 0) + e A τ b ( t − τ) d τ
t
At
donde e es la matriz exponencial y se obtiene considerando las ecuaciones
e A t = β 0 I + β 1 A ………………. ( A )
= β 0 + β 1 λ i …………………( B )
λi t
e
Es necesario determinar los valores característicos correspondientes a la matriz de
coeficientes, por lo que se calcula
1−λ −1
det A − λ I = = λ 2 − 4λ + 4
1 3−λ
de donde
λ 2 − 4λ + 4 = 0
(λ − 2) 2 = 0 ⇒ λ 1 = λ 2 = 2
e 2 t = β 0 + 2 β 1 ………………………………..( C )
Para tener otra ecuación que permita obtener las funciones β 0 y β 1 , se deriva la ecuación
( B ) con respecto a λ ,
138
dλ
d
( )
eλt = β 1
t eλt = β 1
β 1 = t e 2t
β 0 = e 2t − 2t e 2t
Luego de obtener estas funciones se procede a calcular la matriz exponencial, por lo que
At
conviene desarrollar e de la siguiente manera
⎡1 0 ⎤ ⎡1 − 1⎤
e At = β1 ⎢ ⎥ + β1 ⎢ ⎥
⎣0 1⎦ ⎣1 3 ⎦
⎡β + β 1 − β1 ⎤
e At = ⎢ 0
⎣ β1 β 0 + 3 β 1 ⎥⎦
entonces
⎡e 2 t − t e 2 t − t e 2t ⎤
e At = ⎢ ⎥
⎣ te
2t
e 2t + t e 2t ⎦
⎡e 2 t − t e 2 t − t e 2 t ⎤ ⎡1⎤
x = ⎢ ⎥⎢ ⎥
⎣⎢ t e
2t
e 2 t + t e 2 t ⎦⎥ ⎣0⎦
∫
y para el sistema no homogéneo, al efectuar operaciones
⎡ x ⎤ ⎡e − t e ⎤ ⎡ τ e 2τ ⎤
=
⎢ y⎥ ⎢ ⎥ + ⎢ 2τ ⎥ dτ
2t 2t t
⎣ ⎦ ⎢⎣ t e ⎥⎦ ⎢⎣− e − τ e ⎥⎦
2t 2τ
0
⎡ τ 2τ ⎤
⎢ e − e 2τ ⎥
t
⎡ x⎤ ⎡e − te ⎤
1
⎢ y⎥ = ⎢ ⎥ + ⎢ ⎥
2t 2t
⎦⎥ ⎢− 1 e 2 τ − τ e 2 τ + 1 e 2 τ ⎥
2 4
⎣ ⎦ ⎣⎢ t e
⎣⎢ 2 ⎦⎥ 0
2t
2 4
Finalmente
139
⎡3 2t 1⎤
⎢ e − t e 2t + ⎥
⎡ x⎤
1
⎢ y⎥ = ⎢ 1 ⎥
4 2 4
⎣ ⎦ ⎢ t e 2 t − e 2 t + 1⎥
⎢⎣ 2 4 ⎥⎦
1
4
(t )
4) Sea el sistema de ecuaciones diferenciales
x '1 = 3 x 1 + x 2 + f
(t )
1
x '2 = − 3 x 2 + f
(t )
1
( )
2
x1 = e − 1 − 3t
1 3t
9
x2 = 0
RESOLUCIÓN
x' = A x + f (t )
f (t ) = x' − A x
se tiene que
( )
obtenemos el vector x '
⎡1 3 t ⎤
e − 1 − 3t ⎥
x = ⎢9
⎢ ⎥
⎣ 0 ⎦
⎡e 3 t 1⎤
⎢ − ⎥
x' = 3
⎢ 3⎥
⎢⎣ 0 ⎥⎦
Ax = ⎢
⎡1
⎥ 9
( ⎤
⎡3 1 ⎤ ⎢ e 3 t − 1 − 3 t ⎥ )
⎣0 − 3⎦ ⎢ ⎥
⎣ 0 ⎦
140
⎡1 3 t ⎤
⎢ e − − t⎥
Ax = 3
1
⎢ ⎥
3
⎣ ⎦
f (t ) :
0
entonces, al sustituir en
⎡e 3 t 1⎤ ⎡ 1 3 t ⎤
f (t ) = ⎢ 3
− ⎥ − ⎢3 e − 3 − t⎥
1
⎢ 3⎥ ⎢ ⎥
⎢⎣ 0 ⎥⎦ ⎣ 0 ⎦
⎡t ⎤
f (t ) = ⎢ ⎥
⎣0⎦
(t ) =t
de donde
f 1
f 2 (t ) =0
5) Sea el sistema
⎧⎪ x ' = x − 3 y − 9 e t
⎨
⎪⎩ y ' = 3 x + y + 9 e
t
⎡cos 3 t − sen 3 t ⎤
e At = e t ⎢ ⎥
⎣sen 3 t cos 3 t ⎦
RESOLUCIÓN
x ' = A x + b( t )
y su solución está dada por
141
x = e At x ( 0) +
∫ e ( b ( τ) d τ
− τ)
t
At
⎡x ⎤
x ( 0) = ⎢ 0 ⎥
⎣y0 ⎦
Dado que se conoce la matriz exponencial, la solución del sistema homogéneo asociado
x h puede obtenerse de la siguiente manera
⎡cos 3 t − sen 3 t ⎤ ⎡ x 0 ⎤
x h = e A t x ( 0) = e t ⎢ ⎥⎢ ⎥
⎣sen 3 t cos 3 t ⎦ ⎣ y 0 ⎦
⎡ x cos 3 t − y 0 sen 3 t ⎤
= et ⎢ 0 ⎥
⎣ x 0 sen 3 t y 0 cos 3 t ⎦
∫
Para el sistema no homogéneo, una solución particular x p se obtiene según se indica:
e ( b ( τ) d τ
− τ)
xp =
t
At
por lo que se obtienen inicialmente los elementos del integrando, esto es,
∫
⎡I 1 ⎤
e ( b ( τ) d τ = ⎢ ⎥
− τ)
t
⎢⎣I 2 ⎥⎦
A t
∫
0
donde
∫
0
Para I 1 se tiene
142
∫
⎡
I1 = − 9 e ⎢ ( cos 3t cos 3τ + sen 3 t sen 3 τ − 3 τ) dτ +
t
⎣
t
∫
0
⎤
( sen 3t cos 3τ − cos 3 t sen 3 τ) dτ⎥
t
∫
0
( sen 3t − cos 3 t ) ∫
⎤
0
⎦
⎡1
0
⎢⎣ 3 0
⎤
− ( sen 3t − cos 3 t ) cos 3τ ⎥
t
1
3 0⎥⎦
= − 3 e t ( cos 3t + sen 3 t ) sen 3t + 3 e t ( sen 3t − cos 3 t ) ( cos 3t − 1)
I 1 = 3 e t ( cos 3t − sen 3 t − 1)
∫
De manera similar para I 2 se tiene
⎡
I2 = 9e t ⎢ ( cos 3t cos 3τ + sen 3 t sen 3 τ) dτ +
t
∫
0
⎤
− ( sen 3t cos 3τ − cos 3 t sen 3 τ) dτ⎥
t
∫
0
⎡
= 9 e t ⎢( cos 3t − sen 3 t ) cos 3τ dτ +
t
∫
0
⎤
( sen 3t + cos 3 t ) ( sen 3τ) dτ⎥
t
⎦
⎡1
0
⎢⎣ 3
t
⎤
− ( sen 3t − cos 3 t ) cos 3τ ⎥
t
1
3 0⎥⎦
143
Efectuando operaciones y simplificando resulta
I 2 = 3 e t ( − 1 + sen 3t + cos 3 t )
⎡3 e t ( cos 3 t − sen 3 t − 1) ⎤
Entonces, una solución particular del sistema no homogéneo es
xp = ⎢ t ⎥
⎢⎣3 e ( cos 3 t + sen 3 t − 1) ⎥⎦
por lo que la solución es
xG = xh + x p
⎡e t ( x 0 cos 3 t − y 0 sen 3 t ) + 3 e t ( cos 3 t − sen 3 t − 1) ⎤
x = ⎢ t ⎥
⎣⎢e ( x 0 sen 3 t − y 0 cos 3 t ) + 3 e ( cos 3 t + sen 3 t − 1) ⎥⎦
t
⎡e − 2 t − t e − 2 t ⎤
⎡e A t ⎤
−1
⎣ ⎦ = ⎢ ⎥ es la matriz inversa de la matriz exponencial e de
⎣⎢ 0 ⎦⎥
At
− 2t
6) Si
e
la matriz de coeficientes del sistema de ecuaciones diferenciales
x 1 ' = 2 x 1 + x 2 + f (t )
x2' = 2 x 2 + 4 e 2t
y una solución particular del sistema es
( )
⎡ 3t + 2t 2 e 2t ⎤
x p (t ) = ⎢ ⎥ .
⎢⎣ g (t ) ⎥⎦
RESOLUCIÓN
x 1 ' = 2 x 1 + x 2 + f ( t ) ……………….( 1 )
Del sistema
x2' = 2 x 2 + 4 e 2 t ………………….( 2 )
se resuelve la ecuación ( 2 )
x 2 ' − 2 x 2 = 4 e 2t
x 2 ' + P ( t ) x 2 = Q( t )
que es no homogénea de primer orden, de la forma
144
cuya solución está dada por
x2 = x2h + x2 p
donde
x 2h = C1 e ∫
− P ( t ) dt
x 2h = C1 e ∫ = C1 e 2t
2 dt
y también
x2 p = e ∫ ∫ Q( t ) e ∫ dt
− P( t ) dt P( t ) dt
x 2 p = e 2 t ∫ 4 e 2 t e ∫ dt
− 2 dt
x 2 p = e 2t (4 t)
x 2 p = 4 t e 2t
x 1 ' − 2 x 1 = C 1 e 2 t + 4 t e 2 t + f ( t ) …………………( 4 )
(
x1 = 3t + 2t 2 e 2t )
( )
y al derivar esta última ecuación
x1 ' = 3t + 2t 2 2 e 2t + e 2t (3 + 4 t)
(
x 1 ' = 3 + 10 t − 4 t 2 e 2 t )
( ) ( )
por lo que al sustituir en (4):
x 1 ' − 2 x 1 = 3 + 10 t − 4 t 2 e 2 t − 6 t + 4 t 2 e 2 t ………………………..( 5 )
= (3 + 4 t) e 2t
Igualando ( 4 ) y ( 5 )
C 1 e 2 t + 4 t e 2 t + f (t ) = ( 3 + 4 t ) e 2 t
145
(
f (t ) = ( 3 + 4 t ) e 2 t − C 1 + 4 t e 2 t )
( )
f (t ) = 3 − C 1 e 2 t
( )
y de la ecuación ( 3 ) resulta
g (t ) = C 1 + 4 t e 2 t
7) Resuelva el sistema
x1' = x2 + 1
x2' = x1 + e t
2 ⎣⎢ e − e − t e t + e − t ⎦⎥
At
RESOLUCIÓN
x = e At x ( 0) + ∫0 e b ( t − τ) d τ
t Aτ
x = xh +x p
Del sistema se identifican
⎡0 1⎤ ⎡1⎤
A=⎢ ⎥ , b( t ) = ⎢ t ⎥
⎣1 0⎦ ⎣e ⎦
y en este caso las condiciones iniciales están dadas por elementos constantes por lo que se
tiene
146
1 ⎡e + e e t − e −t ⎤ ⎡C1 ⎤
x ( 0) = ⎢ t
−t
xh = e ⎥ ⎢ ⎥
t
2 ⎢⎣e − e − t e t + e − t ⎥⎦ ⎢⎣C2 ⎥⎦
At
1 ⎡e + e ⎤ 1 ⎡e − e ⎤
−t −t
xh = C1 ⎢ t ⎥ + 2⎢ ⎥
t t
2 ⎢⎣e − e − t ⎥⎦ 2 ⎢⎣e t + e − t ⎥⎦
C
1 ⎡e + e e τ − e −τ ⎤ ⎡ 1 ⎤
b ( t − τ) = ⎢ τ
τ −τ
⎥ ⎢ ( t − τ) ⎥
Ατ
2 ⎣⎢e − e −τ e τ + e −τ ⎦⎥ ⎣⎢e ⎦⎥
e
efectuando operaciones
1 ⎡e + e + e − e e ⎤
b ( t − τ) = ⎢ τ
τ −τ − 2τ t
⎥
Ατ
t
2 ⎣⎢e − e −τ + e t + e −2τ e t ⎦⎥
e
∫
por lo que
(e )
⎡1 ⎤
∫
⎢
τ
+ e −τ + e t − e −2t e t d τ⎥
t
e b ( t − τ) = ⎢
∫
xp = ⎥
( )
Ατ
t
⎢1 ⎥
2 0
⎢ e τ − e −τ + e t + e −2t e t d τ⎥
t
⎢⎣ 2 ⎥⎦
0
1⎡ 2t e t + e t − e − t ⎤
= ⎢ ⎥
4 ⎢⎣− 4 + 2 t e t + 3 e t + e − t ⎥⎦
finalmente resulta
x = xh + x p
1 ⎡e t + e − t ⎤ 1 ⎡e t − e − t ⎤ 1 ⎡ 2t e t + e t − e − t ⎤
x = C1 ⎢ t −t ⎥
+ C2 ⎢ t −t ⎥
+ ⎢ −t⎥
2 ⎣e − e ⎦ 2 ⎣e + e ⎦ 4 ⎣− 4 + 2t e + 3e + e ⎦
t t
147
EJERCICIOS PROPUESTOS
⎡cos t − sen t ⎤
1) Dada la matriz exponencial
e Αt = e 2t ⎢ ⎥
⎣sen t cos t ⎦
A= ⎢ ⎥
⎣c d ⎦
Solución.
β 0 = e 2 t ( cos t − 2 sen t )
β 1 = e 2 t ( sen t )
⎡2 − 1⎤
A=⎢ ⎥
⎣1 2 ⎦
x1 = x 2 + 1
x 2 = − x 1 + cot ( t )
⎛ π⎞ ⎛ π⎞
con las condiciones x 1 ⎜ ⎟ = 0 , x 2 ⎜ ⎟ = 0 .
⎝ 2⎠ ⎝ 2⎠
⎛ sen t ⎞
Solución.
x 1 = ( sen t ) Ln ⎜ ⎟
⎝ cos t + 1⎠
⎛ sen t ⎞
x 2 = ( cos t ) Ln ⎜ ⎟
⎝ cos t + 1⎠
148
Αt
3) Obtenga la matriz exponencial e , donde
⎡0 1 0⎤
A = ⎢⎢0 2 0⎥⎥
⎢⎣1 0 0⎥⎦
Solución.
⎡ e 2t −1 ⎤
⎢1 0⎥
⎢ ⎥
= ⎢0 0⎥
2
e Αt
⎢ ⎥
e 2t
⎢ e 2 t − 1 − 2t ⎥
⎢t 1⎥
⎣ 4 ⎦
x2 = 3 x1 + 2 x2
( )
Solución.
x1 = 4 t e 3t + 9 e −t + 7 e 3t
1
( )
16
x2 = 12 t e 3 t − 9 e − t + 9 e 3 t
1
16
⎡0 2⎤ ⎡ 0 0⎤
5) Dadas las matrices A = ⎢ ⎥ y B =⎢ ⎥
⎣0 0⎦ ⎣− 2 0 ⎦
( A +B ) t
a) verifique si se cumple e ⋅ e =e
At Bt
b) de acuerdo con el resultado del inciso a) ¿se puede afirmar que la expresión
e At ⋅ e B t = e (
A +B ) t
es falsa para toda matriz cuadrada. Justifique la respuesta.
149
Solución.
La igualdad no se cumple, sólo se cumple cuando las matrices A y B son conmutativas.
Solución.
x = 1 − e −t
y = x' = e − t
Solución.
⎡x' 1 ⎤ ⎡ 0 1 ⎤ ⎡x 1 ⎤ ⎡ 0 ⎤
⎢ x ' ⎥ = ⎢− 2 − 2⎥ ⎢ x ⎥ + ⎢e −t cos t ⎥
⎣ 2⎦ ⎣ ⎦ ⎣ 2⎦ ⎣ ⎦
= y + z
dx
dt
= z + x
dy
dt
= x + y
dz
dt
150
( ( ) ( ) ( ))
Solución.
x = C1 6 e − t + 3 e 2 t + C 2 − 3 e − t + 3 e 2 t + C 3 − 3 e − t + 3 e 2 t
1
a
y =
1
a
( ( ) ( ) (
C1 − 3e − t + 3e 2 t + C2 6 e − t + 3e 2 t + C3 − 3e − t + 3e 2 t ))
z =
1
a
( ( ) ( ) (
C1 − 3 e − t + 3 e 2 t + C 2 − 3 e − t + 3 e 2 t + C 3 6 e − t + 3 e 2 t ))
⎡ − 1 − 2⎤
9) Si se sabe que para A = ⎢
⎣− 2 2 ⎥⎦
, la matriz exponencial
⎡ e 3 t + 4 e −2 t − 2 e 3 t + 2 e −2 t ⎤
Αt
= ⎢ ⎥
1
⎣− 2 e + 2 e 4 e 3 t + e −2 t ⎦
−2 t
e 3t
5
(D + 1) x + 2 y = 0
calcule la solución del sistema
2 x + ( D − 2) y = t
x ( 0) = 0 , y ( 0) = 0 .
que cumpla con las condiciones
Solución.
x =− e + + t −
2 3t 1 − 2t 1 1
e
45 10 3 18
y =− e + + t −
4 3t 1 − 2t 1 5
e
45 20 6 36
151
obtenga un sistema de ecuaciones diferenciales de primer orden equivalente y expréselo
en forma matricial.
Solución.
⎡ x ' 1 ⎤ ⎡0 1 0 0⎤ ⎡x 1 ⎤ ⎡ ⎤
⎢ x ' ⎥ ⎢0 0 1 0⎥⎥ ⎢x ⎥ ⎢ 0 ⎥⎥
0
⎢ 2⎥ = ⎢ ⎢ 2⎥ + ⎢
⎢ x ' 3 ⎥ ⎢0 0 0 1⎥ ⎢x 3 ⎥ ⎢ 0 ⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎣ x ' 4 ⎦ ⎣5 −7 0 5⎦ ⎣ x 4 ⎦ ⎣− 2 e −t ⎦
152
TEMA 5
TRANSFORMADA DE LAPLACE
(
f (t ) = 1 + e 2t )
2
RESOLUCIÓN
f (t ) = 1 + 2 e 2t + e 4 t
∫
y la definición de Transformada de Laplace es:
{ f ( t )} e − s t f ( t ) dt
∞
L =
∫
0
{ f ( t )} = lim e − s t f ( t ) dt
b
L
b→∞ 0
{ f ( t )} ∫ ( )
∞
= e − s t 1 + 2 e 2 t + e 4 t dt
∫
L
{ f ( t )} = lim ( )
0
e − s t 1 + 2 e 2 t + e 4 t dt
b
L
b→∞ 0
{ f ( t )} ∫ ∫ ∫
⎡ ⎤
= lim ⎢ e − s t dt + 2 e − s t e 2 t dt + e − s t e 4 t dt ⎥
b b b
⎣ ⎦
L
b→∞ 0 0 0
{ f ( t )} ∫ ∫
⎡ 1 ⎤
= lim ⎢ − e − s t + 2 e ( ) dt + e ( ) dt ⎥
− s−2 t − s−4 t
b b b
⎣⎢ s ⎦⎥
L
b→∞ 0 0
0
153
⎡ 1 ⎤
{ f ( t )} = lim ⎢ − e − s t − e ( ) − e ( ) ⎥
b b
− s−2 t − s−4 t
( s − 2) ( s − 4)
b
2 1
b→∞ ⎢ s ⎥
L
⎣ 0 0 0⎦
( ) ( )
Evaluando en los extremos de integración y aplicando el límite a cada uno de los términos
{ f ( t )} ( )
⎡ 1 ⎤
= lim ⎢ − e − s b − 1 − e ( ) −1 − e ( ) −1 ⎥
( s − 2) ( s − 4)
2 − s−2 b 1 − s−4 b
⎣ s ⎦
L
b→∞
L { f ( t )} =
1
s
+
2
s − 2
+
1
s − 4
Así resulta
L { f ( t )} = F (s) =
1
s
+
2
s − 2
+
1
s − 4
RESOLUCIÓN
154
⎧0 0 ≤ t < 1
f (t ) = ⎨
⎩t t ≥ 1
{ f ( t )} = ∫ e − s t f ( t ) dt =
∫ e − s t ( 0 ) dt +
∫
∞ ∞
t e − s t dt
1
L
∫
0 0 1
{ f ( t )}
⎡ ⎤
= lim ⎢ t e − s t dt ⎥
b
b→∞⎣ ⎦
L
1
∫
Integrando por partes
{ f ( t )}
⎡ t ⎤
= lim ⎢ − e − s t + −s t
dt ⎥
b
1
b→∞⎣ ⎦1
L e
s s
{ f ( t )}
⎡ t ⎤
= lim ⎢ − e − s t − 2 e − s t ⎥
b
1
b→∞⎣ s ⎦1
L
s
{ f ( t )}
⎡ b ⎤
= lim ⎢ − e − s b − 2 e − s b + e − s + 2 e − s ⎥
1 1 1
b→∞⎣ ⎦
L
s s s s
{ f ( t )}
⎡ 1 −s ⎤
= ⎢ + 2 e −s ⎥
1
⎣ s ⎦
L e
s
⎡ 1 1 ⎤
⇒ F (s) = e −s ⎢ + 2⎥
⎣ s s ⎦
155
RESOLUCIÓN
{ f ( t )} = ∫ ∫ ∫
π
f (t ) d t = e − s t ( 0) d t +
∞ ∞
−s t
e − s t cos t d t
2
π
L e
{ f ( t )} ∫
0 0
2
= lim e − s t cos t dt
b
π
L
b→∞
2
∫ ∫
⎡ cos t − s t 1 ⎛ s ent − s t ⎞⎤
cos t dt = lim ⎢ − − ⎜− + cos t dt ⎟ ⎥
b
−s t −s t
b
1
π
b→∞⎣ s⎝ ⎠⎦ π
lim e e e e
b→∞ s s s
∫ ∫
2
⎡ cos t − s t ⎤
2
e− s t cos t d t ⎥
b
1 1
b→∞ ⎣ ⎦π
π
lim e
b→∞ s s s
2
2
∫
⎛ cos t − s t s ent − s t ⎞
cos t dt = lim ⎜ − + ⎟
b
−s t
b
π
b→∞⎝ ⎠π
lim e e e
b→∞ s s2
2
2
⎛ − s ( cos t ) e − s t + ( s ent ) e − s t ⎞
∫
⎜ ⎟
b
e −s t cos t dt = lim ⎜ s2 ⎟
b→∞⎜ ⎟
b
π
s2 + 1
lim
⎜⎜ ⎟⎟
b→∞
⎝ ⎠π
2
s2
2
∫
⎛ − s ( cos t ) e − s t + ( s ent ) e − s t ⎞
cos t dt = lim ⎜ ⎟⎟
b
−s t
⎜
b
π
b→∞⎝ s2 + 1 ⎠π
lim e
b→∞
2
2
⎡ ⎛π ⎞ − s
π
⎛π ⎞ − s
π
⎤
∫
⎢ − s cos b e + s enb e − s b − s cos ⎜ ⎟ e 2 + s en ⎜ ⎟ e 2 ⎥
−s b
b→∞⎢ s + 1 ⎥
π
lim
b→∞
⎢ ⎥
2
⎣ ⎦
2
156
∫
π
−
−e 2
cos t d t = 2
s
−s t
b
s + 1
lim e
b→∞ 0
{ f ( t )}
Finalmente
−1
= F (s) = 2
π
− s
s + 1
L e 2
157
5.2 LA TRANSFORMADA DE LAPLACE COMO UN OPERADOR
LINEAL.
RESOLUCIÓN
( )
⎧ cos ht ⎫ ⎧1 ⎫
L ⎨ − 3t ⎬ = L ⎨ e + e−t e ⎬
⎩e ⎭ ⎩2 ⎭
t 3t
{ e } + L { e }⎤⎦
⎧ cos ht ⎫ 1
L ⎨ − 3 t ⎬ = ⎡L
⎩e ⎭ 2⎣
4t 2t
⎡ 1 1 ⎤
= ⎢ + ⎥
1
2 ⎣s − 4 s − 2⎦
1 ⎡s − 2 + s − 4⎤
= ⎢ ⎥
2 ⎣ s 2 − 6s + 8 ⎦
1 ⎡ 2 ( s − 3) ⎤
= ⎢ ⎥
2 ⎣s 2 − 6s + 8⎦
De donde
s − 3
F (s) =
s 2
− 6s + 8
158
5.3 TEOREMA DE TRASLACIÓN EN EL DOMINIO DE “ s ”.
{e }
1) Determine
t−2
a) L sen 3 t
⎧ s + 4 ⎫
⎨ ⎬
−1
⎩ s − 2s + 5 ⎭
b) L
RESOLUCIÓN
{
f ( t ) = e t e − 2 sen 3 t }
{e }
por lo que se tiene
−2
L e t sen 3 t
{e }
De la linealidad de la transformada de Laplace resulta
−2 t
e L sen 3 t
{e }
y de acuerdo al primer teorema de traslación, expresado como sigue
L at
f (t ) = F ( s − a )
donde F ( s ) = L { f ( t )}
{e }
se obtiene
⎡ ⎤
sen 3 t = e ⎢ ⎥
⎢⎣ ( s − 1) + 9 ⎥⎦
−t −2 −2 3
L 2
⎡ −2 ⎤
= ⎢ ⎥
⎣⎢ ( s − 1)
3e
2
+ 9 ⎦⎥
159
⎪⎧ s + 4 ⎪⎫ ⎪⎧ ⎪⎫ ⎪⎧ ⎪⎫
⎨ ⎬ = L ⎨ ⎬ + L ⎨ ⎬
⎩⎪ ( s − 1) + 4 ⎭⎪ ⎩⎪ ( s − 1) + 4 ⎭⎪ ⎩⎪ ( s − 1) + 4 ⎪⎭
−1 −1 s −1 4
L 2 2 2
⎪⎧ s + 4 ⎪⎫ ⎪⎧ s − 1 + 1 ⎪⎫ ⎪⎧ ⎪⎫
= ⎨ ⎬ = L ⎨ ⎬ + 4L ⎨ ⎬
⎪⎩ ( s − 1) + 4 ⎪⎭ ⎪⎩ ( s − 1) + 4 ⎪⎭ ⎪⎩ ( s − 1) + 4 ⎪⎭
−1 −1 −1 1
L 2 2 2
⎧⎪ s − 1 ⎫⎪ ⎧⎪ ⎫⎪
= L ⎨ ⎬ + 5L ⎨ ⎬
⎪⎩ ( s − 1) + 4 ⎭⎪ ⎩⎪ ( s − 1) + 4 ⎭⎪
−1 −1 1
2 2
{ f ( t )}
⎧ ⎛ t ⎞⎫
= F ( s ) , compruebe que L ⎨ f ⎜ ⎟ ⎬ = a F ( a s )
⎩ ⎝ a ⎠⎭
2) Sea L
donde f ( t ) = t + e(
2 2 t + 1)
RESOLUCIÓN
Si f ( t ) = t 2
+ e 2t e , entonces
⎛t⎞ ⎛t⎞
⎛t⎞
2⎜ ⎟
f⎜ ⎟ = ⎜ ⎟ + ee ⎝a⎠
2 2
⎝a⎠ ⎝a⎠
⎛ 2⎞
⎜ ⎟t
= 2 + ee ⎝a⎠
t2
a
Transformando en Laplace resulta
{t }
⎧ ⎛ t ⎞⎫ ⎧⎪ ⎛⎜ 2 ⎞⎟ t ⎫⎪
L ⎨ f ⎜ ⎟⎬ = 2 L + e L ⎨e ⎝ a ⎠ ⎬
1
⎩ ⎝ a ⎠⎭ ⎪⎩ ⎭⎪
2
160
⎧ ⎛ t ⎞⎫
L ⎨ f ⎜ ⎟⎬ = 2 3 + = 2 3 +
as−2
2 e 2 e
⎩ ⎝ a ⎠⎭ a s s−
2 a s
a a
= +
2 ae
as−2
2 3 …………………...( A )
a s
{ f ( t )} = F (s)
Por otro lado
L
y al transformar f ( t ) se tiene
F (s) = L {t }2
{ }
+ eL e 2 t
= +
2 e
s 3
s−2
por lo que F ( a s ) = +
(a s)
2 e
3
as−2
a F (a s) = +
2a ae
3 3
a s as − 2
a F (a s) = +
2 ae
as − 2
2 3 ………………..( B )
a s
Al comparar (A) y (B) se comprueba que
⎧ ⎛ t ⎞⎫
L ⎨f ⎜ ⎟⎬ = a F ( a s )
⎩ ⎝ a ⎠⎭
161
5.4 TEOREMA DE TRASLACIÓN EN EL DOMINIO DE “ t ”.
⎧ t sen a t ⎫
( )
L ⎨ ⎬ =
s
⎩ 2a ⎭ s + a2
y 2 2
⎪⎩ ⎡⎣( s − 2 ) + 9 ⎤⎦
s
⎭⎪
L 2
RESOLUCIÓN
⎪⎩ ⎡⎣( s − 2 ) + 9 ⎤⎦ ⎪⎩ ⎡⎣( s − 2 ) + 9 ⎤⎦
s
⎪⎭ ⎭⎪
L 2
= L 2
⎧ ⎫
⎪ ⎪
( )
−1
⎨ 2⎬
=
s t sen a t
⎪⎩ s + a ⎪⎭
L 2 2
2a
162
⎧ ⎫
⎪ ⎪
( )
⎨ ⎬ = =
2 ( 3)
−1 s t sen 3 t 2 t e 2 t t sen 3 t
⎪⎩ s + 9 ⎪⎭
L e
( s − 2) → s
2 2 6
también
⎧ ⎫ ⎡ sen 3 t − 3 t cos 3 t ⎤ 2 t
⎪ ⎪
( )
⎨ ⎬ = ⎢ ⎥e
2 ( 3) 3
−1 1
⎪⎩ s + 9 ⎪⎭ ⎢⎣ ⎥⎦
2L 2
( s − 2) → s
2 2
= ( sen 3 t − 3 t cos 3 t )
e 2t
27
Finalmente
⎧ ⎫
⎪ ⎪
−1
⎨ 2⎬
= + ( sen 3 t − 3 t cos 3 t )
⎪⎩ ⎣⎡( s − 2 ) + 9 ⎦⎤
s e 2t t sen 3 t e 2t
⎪⎭
L 2
6 27
RESOLUCIÓN
x ( t ) = e − 3 t cos 2 ( t − 4 ) u ( t − 4 )
o bien
x (t ) = e e − 12 cos 2 ( t − 4 ) u ( t − 4 )
− 3( t − 4 )
{g (t − a ) u ( t − a )} = e − a s G ( s )
del segundo teorema de traslación
L
{ x (t ) } {e }
entonces
X (s) = L = L e − 12 cos 2 ( t − 4 ) u ( t − 4 )
− 3 (t − 4)
donde
g (t − 4) = e e − 12 cos 2 ( t − 4 )
− 3 (t − 4)
de aquí se obtiene
163
g ( t ) = e − 3 t e − 12 cos 2 t
{ }
y al transformar
G ( s ) = e − 1 2 L e − 3 t cos 2 t
164
5.5 TRANSFORMADA DE LAPLACE DE LAS FUNCIONES
IMPULSO, ESCALÓN Y RAMPA.
⎧ ; 0≤t <1
f (t ) = ⎨ 2
1
⎩t − 2 t + 2 ; 1 ≤ t
RESOLUCIÓN
f ( t ) = u ( t ) + u ( t − 1) ⎡⎣( t − 1) 2 + 1⎤⎦ − u ( t − 1)
= u ( t ) + u ( t − 1) ( t − 1) 2 + u ( t − 1) − u ( t − 1)
f ( t ) = u ( t ) + u ( t − 1) ( t − 1) 2
al simplificar se obtiene
{ f ( t )} {u ( t )} {u ( t }
de esta manera la obtención de la transformada de Laplace resulta
L = L + L − 1) ( t − 1) 2
{u ( t − a ) g ( t − a )} = e − a s G ( s ) , donde G ( s ) = L { g( t ) }
y del segundo teorema de traslación
L
{u (t }=
Interesa determinar
L − 1) ( t − 1) 2 e − s G (s)
g ( t − 1) = ( t − 1) 2 ⇒ g (t ) = t 2
{t } =
y resulta
G (s) = L 2 2
s3
165
entonces
L {u ( t − 1) g ( t − 1)} = 2e−s
s3
y para la función u ( t )
L {u ( t )} =
1
s
Finalmente
F (s) =
2e−s
+
1
s s3
Es conveniente mencionar que la transformada obtenida, se puede determinar también
empleando la definición de transformada de Laplace.
RESOLUCIÓN
f ( t ) = u ( t ) cos 2 t − u ( t − 2π ) cos 2 t
la función escalón unitario:
f ( t ) = ( cos 2 t ) u ( t ) − ⎡⎣cos 2 ( t − 2π ) ⎤⎦ u ( t − 2π )
entonces
{ f ( t )} {( cos 2 t ) u ( t )} {cos 2 ( t − 2π ) u ( t − 2π )}
transformando en Laplace y aplicando el segundo teorema de traslación
L = L − L
⎛ ⎞
= − e − 2π s ⎜ 2 ⎟
s s
s + 4
2
⎝ s + 4 ⎠
166
Finalmente
F (s) =
s e − 2π s
−
s
s2 + 4 s2 + 4
⎛ e − 2π s ⎞
F (s) = ⎜ − ⎟
s
s2 + 4 ⎝ s2 + 4⎠
1
{ f ( t )} ∫ e − s t f ( t ) dt =
∫
∞ 2π
L = e − s t cos 2 t dt
0 0
∫
Integrando dos veces por partes y simplificando
2π
2π
2 e − s t sen 2 t − s e − s t cos 2 t
−s t
cos 2 t dt =
s2 + 4
e
0
0
s e − 2π s
= −
s
s2 + 4 s2 + 4
( )
Finalmente
L { f ( t )} =
s
s + 4
2
1 − e − 2π s
RESOLUCIÓN
167
f ( t ) = r ( t − 1) − r ( t − 2 ) − r ( t − 3 ) + r ( t − 4 )
( )
s s s s
RESOLUCIÓN
⎧ 0 ; t<0 ⎧ 0 ; t<a
u (t ) = ⎨ ; u (t − a ) = ⎨
⎩ 1 ; t≥ 0 ⎩ 1 ; t≥ a
⎧ 0 ; t <0 ⎧ 0 ; t<a
r (t ) = ⎨ ; r (t − a ) = ⎨
⎩ t ; t≥ 0 ⎩ t − a ; t≥ a
168
f (t ) = 3t 2 u ( t ) − 3 t 2 u ( t − 1) + 3 r ( t − 2 ) − 3 u ( t − 3 )
− 3 r (t − 3) + 3 u ( t − 4 ) − 3 u ( t − 5 )
⎧ t ; 0 < t < 2
⎪
5) Sea f ( t ) = ⎨ 4 ; 2 < t < 4
⎪ 0 ; t > 4
⎩
RESOLUCIÓN
a) a partir de esta gráfica resulta más fácil expresar f ( t ) en términos de las funciones
r (t ) y u (t ) :
f ( t ) = r ( t ) + 2 u ( t − 2) − r ( t − 2) − 4 u ( t − 4 )
f (t ) = L { f ( t )} {r ( t )} {u ( t − 2 )} − L {r ( t − 2 )} − 4 L {u ( t − 4 )}
b)
= L + 2L
F (s) =
Finalmente
+ − 2 e −2s − e −4s
1 2 −2s 1 4
2
e
s s s s
169
5.6 DERIVADA DE LA TRANSFORMADA DE UNA FUNCIÓN.
⎪⎧ ⎛ s + 2 ⎞ ⎪⎫
⎨ Ln ⎜ ⎟⎬
−1
⎩⎪ ⎝ s − 2 ⎠ ⎭⎪
1) Calcule L
RESOLUCIÓN
−1 ⎪⎧ ⎛ s + 2 ⎞ ⎪⎫
⎨ Ln ⎜ ⎟⎬ = L
−1
{Ln ( s + 2 ) − Ln ( s − 2 )}
⎩⎪ ⎝ s − 2 ⎠ ⎭⎪
L
F (s)
Del teorema de la derivada de una transformada
{t f ( t )} = ( − 1) ; F (s) = L { f ( t )}
n
n n d
L
dsn
y si F ( s ) = Ln ( s + 2 ) − Ln ( s − 2 ) al derivar se tiene
d F (s)
= − .............................. ( A )
1 1
ds s + 2 s − 2
Ahora bien, si se antitransforma el teorema referido de la sugerencia
⎧⎪ F ( s ) ⎫⎪
⎨( − 1) ⎬ = t f (t )
−1
n
n d
⎩⎪ ⎭⎪
n
L n
ds
para n = 1 resulta
⎧ d F (s) ⎫
⎨− ⎬ = t f (t )
−1
⎩ ds ⎭
L
d F (s)
por lo que de la derivada obtenida en (A)
− = − +
1 1
ds s + 2 s − 2
Antitransformando resulta
⎧ 1 ⎫
⎨− + ⎬ = t f (t )
−1 1
⎩ s + 2 s − 2⎭
L
170
− e − 2 t + e 2 t = t f (t )
pero también
− e − 2 t + e 2 t = 2 sen h 2 t
t f ( t ) = 2 sen h 2 t
por lo que
de donde se tiene
f (t ) =
2 sen h 2 t
t
171
5.7 TEOREMA DE CONVOLUCIÓN.
⎪⎩ ( s − 1) s + 4 ⎭⎪
2
L 2
RESOLUCIÓN
{H ( s )} , donde H ( s ) =
(s − 1) s 2 + 4( )
−1 2
Se tiene L
⎡ 1 ⎤ ⎡ ⎤
( )
⎢ ⎥ ⎢ ⎥
( )
2
⎣ − ⎢
⎦ ⎣ + ⎥
⎦
2
s 1 s 4
⎩⎪ ⎣
L
⎭
2
{F ( s ) G ( s )} = f (t ) * g (t )
Del teorema de Convolución
−1
∫
L
−1
{F ( s ) G ( s )} = f (τ ) g ( t − τ ) dτ
t
L
0
Entonces se tiene
{F ( s )} − 1⎪
⎧ 1 ⎫⎪
g (t ) = L = L ⎨ ⎬ = e
⎩⎪ ( s − 1) ⎪⎭
−1 t
⎧ ⎫
f (t ) = L {G ( s )} = L − 1⎪ ⎪
( )
⎨ 2 ⎬ = sen 2 t
−1 2
⎪⎩ s − 4 ⎪⎭
g ( t − τ ) , así
172
f (τ ) = sen 2 τ
g (t − τ ) = e (
t −τ )
∫
sustituyendo estas funciones en la integral de convolución
f (t ) * g (t ) = ( sen 2 τ ) e (t − τ ) dτ
t
∫
0
f (t ) * g (t ) = e t e − τ sen 2 τ dτ
t
∫ (− )
La integral se calcula integrando dos veces por partes, de donde se obtiene
(− )
A partir de este resultado se tiene
f (t ) * g (t ) = 2 e − t cos 2 t − e − t sen 2 t + 2
et
(− )
5
f (t ) * g (t ) = 2 cos 2 t − sen 2 t + 2e
1 t
5
De donde
( 2 cos 2 t )
⎧⎪ ⎫⎪
( )
⎨ ⎬ = − + sen 2 t − 2e
−1
⎪⎩ ( s − 1) s + 4 ⎪⎭
2 1 t
L 2
5
173
5.8 RESOLUCIÓN DE ECUACIONES Y SISTEMAS DE
ECUACIONES DIFERENCIALES LINEALES.
RESOLUCIÓN
{ y ''} { y} {sen 3 t}
Transformando en Laplace ambos miembros de la ecuación diferencial se tiene
L + 4L = L
s 2 Y ( s ) − s Y (0 ) − Y '(0 ) + 4 Y ( s ) =
3
s + 9
2
s 2 Y (s) + 4 Y (s) =
3
s2 + 9
( )
esto es
Y (s) s 2 + 4 =
3
s + 9
2
Y (s) =
(s ) (s )
3
2
+ 4 2
+ 9
As + B Cs + D
Y (s) =
(s ) (s )
= +
3
2
+ 4 2
+ 9 s + 4
2
s2 + 9
de aquí resulta
174
( )
3 = ( A s + B ) s 2 + 9 + (C s + D ) s 2 + 4 ( )
3 = s 3( A + C ) + s 2 ( B + D ) + s ( 9 A + 4 C ) + 9 B + 4 D
Y (s) = 2 5
3 3
− 25
s + 4 s + 9
Al antitransformar
{Y ( s )}
⎧ 1 ⎫ ⎧ 1 ⎫
y (t ) = L −1
= −1
⎨ 2 ⎬ − L
−1
⎨ 2 ⎬
3 3
⎩s + 4⎭ ⎩s + 9⎭
L
5 5
por lo que la solución de la ecuación diferencial es
y (t ) = sen 2 t −
3 3
sen 3 t
10 15
o bien
y (t ) = sen 2 t − sen 3 t
3 1
10 5
y '' − 7 y ' + 6 y = 10 δ ( t − 2 )
2) Aplique la transformada de Laplace para resolver la ecuación diferencial
175
RESOLUCIÓN
{ y ''} − 7L { y '} + 6 L { y} = 10 L {δ ( t − 2 )}
s 2 Y (s) s Y ( 0 ) − Y ' ( 0 ) − 7 ⎡⎣ s Y ( s ) − Y ( 0 ) ⎤⎦ + 6 Y ( s )
L
− = 10 e − 2 s
s 2 Y ( s ) − 7 s Y ( s ) + 6 Y ( s ) = 10 e − 2 s
Al efectuar operaciones y aplicar las condiciones iniciales resulta
En esta última ecuación interesa obtener Y ( s ) , por lo que se factoriza y se despeja según
( )
se muestra
Y ( s ) s 2 − 7 s + 6 = 10 e− 2 s
Y (s) =
10 e− 2 s
s2 − 7s + 6
Y (s) =
10 e − 2 s
( s − 6 ) ( s − 1)
Teniendo presente que la solución de la ecuación diferencial está dada por y ( t ) , es
necesario obtener la antitransformada de la función Y ( s ) , esto es
{Y ( s )}
⎧⎪ ⎫⎪
y (t ) = L
10 e − 2 s
= L ⎨ ⎬
⎪⎩ ( s − 6 )( s − 1) ⎭⎪
−1 −1
Y ( s ) = e− 2 s
( s − 6 )( s − 1)
10
{e− a s G ( s )}
que llevaría a la aplicación del segundo teorema de traslación en forma inversa, es decir
L −1
= g (t − a ) u (t − a )
176
G (s) = = +
( s − 6 )( s − 1) s − 6 s − 1
10 A B
10 = A ( s − 1) + B ( s − 6 )
siguiendo este proceso se tiene
{G ( s )}
⎧ 1 ⎫ ⎧ 1 ⎫
g (t ) = L −1
= 2L −1
⎨ ⎬ − 2L
−1
⎨ ⎬
⎩ s − 6 ⎭ ⎩ s − 1⎭
⇒ g (t ) = 2 e 6 t − 2 e t
g ( t − a ) = g ( t − 2) = 2 e ( ) − 2 e ( )
6 t−2 t−2
y (t ) = L {Y ( s )} = L
⎪⎧ − 2 s
⎨e
⎪⎫
⎬ = g ( t − 2) u ( t − 2)
( s − 6 )( s − 1) ⎪⎭
−1 −1 10
⎪⎩
por lo que finalmente
y ( t ) = ⎡ 2 e ( ) − 2 e ( ) ⎤ u ( t − 2)
6 t−2 t−2
⎣ ⎦
y '' + y = t − ( t − 4 ) u ( t − 2 )
RESOLUCIÓN
177
y '' + y = t − ( t − 4 ) u ( t − 2 )
{ y ''} { y} {t} {( t − 4 ) u ( t − 2 )}
Enseguida se aplica la transformada de Laplace a cada uno de los términos
L + L = L − L
de donde se obtiene
⎛1 − 2 s⎞
s 2 Y ( s ) − s Y (0 ) − Y '(0 ) + Y ( s ) = − e− 2 s ⎜ ⎟
1
⎝ s ⎠
2 2
s
Al aplicar condiciones iniciales y reducir términos resulta
⎛1 − 2 s⎞
s 2 Y (s) − 1 + Y (s) = − e− 2 s ⎜ ⎟
1
⎝ s ⎠
2 2
s
( )
⎛1 − 2 s⎞
Y (s) s 2 + 1 = 1 + − e− 2 s ⎜ ⎟
1
⎝ s ⎠
2 2
s
despejando a la función de interés
s2 + 1 ⎛ 1 − 2s ⎞
Y (s) =
( ) ( )
− e− 2 s ⎜ 2 2 ⎟
+ 1 ⎜s s + 1 ⎟
⎝ ⎠
s2 s2
⎧ ⎛ 1 − 2 s ⎞ ⎫⎪
y (t ) = L {Y ( s )} = L − 1⎧ 1⎫ − 1⎪ − 2 s
( )
−1
⎨ 2⎬ − L ⎨e ⎜ ⎟⎬
⎩s ⎭ ⎪⎩ ⎜ s 2 s 2 + 1 ⎟⎪
⎝ ⎠⎭
En el segundo miembro, al antitransformar el segundo término se tiene
⎧ ⎛ 1 − 2 s ⎞ ⎫⎪
− 1⎪ − 2 s
( )
⎨e ⎜ ⎟⎬
⎪⎩ ⎜ s 2 s 2 + 1 ⎟⎪
L
⎝ ⎠⎭
{e }
y de acuerdo al segundo teorema de traslación en su forma inversa
L −1 −as
G ( s ) = g (t − a ) u (t − a )
1 − 2s Cs + D
G (s) =
( )
= + 2+
A B
s2 s2 + 1 s s s2 + 1
( )
1 − 2 s = A( s ) s 2 + 1 + B s 2 + 1 + (C s + D ) s 2 ( )
178
1 − 2 s = s 3( A + C ) + s 2( B + D ) + s ( A ) + B
de donde resulta
A + C = 0 ⇒ C = 2
B + D = 0 ⇒ D = − 1
A = − 2
B = 1
entonces
2s − 1
G (s) = − + 2+ 2
2 1
s s s + 1
g (t ) = L −1
{G ( s )} = − 2L
⎨ ⎬ + L
− 1⎧ 1 ⎫ − 1⎧ 1⎫
⎨ 2 ⎬ + 2L
− 1⎧
⎨ 2
s ⎫
⎬ − L
− 1⎧
⎨ 2
1 ⎫
⎬
⎩s ⎭ ⎩s ⎭ ⎩ s + 1⎭ ⎩ s + 1⎭
g ( t ) = − 2 + t + 2 cos t − sent
g ( t − 2 ) = − 2 + ( t − 2 ) + 2 cos ( t − 2 ) − sen ( t − 2 )
g ( t ) = t − ⎡⎣ − 2 + ( t − 2 ) + 2 cos ( t − 2 ) − sen ( t − 2 ) ⎤⎦ u ( t − 2 )
Finalmente, para la función solución
y ( t ) = t + ⎡⎣ 4 − t + sen ( t − 2 ) − 2 cos ( t − 2 ) ⎤⎦ u ( t − 2 )
y '' + 4 y ' + 3 y = 1 + δ ( t − 3 )
4) Utilice la transformada de Laplace para resolver la ecuación diferencial
RESOLUCIÓN
s 2 Y ( s ) − s Y ( 0 ) − Y ' ( 0 ) + 4 ⎡⎣ s Y ( s ) − Y ( 0 ) ⎤⎦ + 3 Y ( s ) = + e−3s
1
s
179
posteriormente se aplican condiciones iniciales y se simplifica, es decir,
( )
s
Y (s) s 2 + 4 s + 3 = + e−3s + s + 5
1
s
1 + s e−3s + s 2 + 5 s
Y ( s ) ( s + 3 ) ( s + 1) =
s
s2 + 5s + 1
Y (s) =
s e−3s
+
s ( s + 3 )( s + 1) s ( s + 3 )( s + 1)
y (t ) = L −1
{Y ( s )}
⎪⎧ s + 5 s + 1 ⎪⎫ ⎪⎧ − 3 s ⎪⎫
y (t ) = L ⎨ ⎬ + L ⎨e ⋅ ⎬
⎩⎪ s ( s + 3 )( s + 1) ⎪⎭ ( s + 3)( s + 1) ⎪⎭
−1 −1
2
1
⎩⎪
se tienen entonces dos funciones a antitransformar, a saber
⎪⎧ s + 5 s + 1 ⎪⎫
I = L ⎨ ⎬
⎪⎩ s ( s + 3 )( s + 1) ⎭⎪
−1
2
s 2 + 5 s + 1 = A ( s + 3 ) ( s + 1) + B s ( s + 1) + C s ( s + 3 )
para s = 0 : A =
1
3
para s = − 3 : B = −
5
6
para s = − 1 : C =
3
2
por lo que resulta
⎧1⎫ ⎧⎪ ⎫⎪ ⎧⎪ 1 ⎫⎪
I = ⎨ ⎬ − L ⎨ ⎬ + L ⎨ ⎬
⎪⎩ ( s + 3 ) ⎭⎪ ⎩⎪ ( s + 1) ⎭⎪
1 −1 5 −1 1 3 −1
⎩s⎭
L
3 6 2
I = − e − 3t + e −t
1 5 3
3 6 2
Por otro lado
180
⎧⎪ − 3 s ⎫⎪
II = L ⎨e ⋅ ⎬
( s + 3)( s + 1) ⎭⎪
−1 1
⎩⎪
{e }
Esta función se antitransformará considerando el segundo teorema de traslación, esto es
L −1 −as
G ( s ) = g (t − a ) u (t − a )
Entonces, la función G ( s ) es
(s + 3 )( s + 1)
1
, y por fracciones parciales
= +
(s + 3 )( s + 1)
1 A B
s + 3 s + 1
1 = A ( s + 1) + B ( s + 3 )
para s = − 1 : B =
1
2
para s = − 3 : A = −
1
⎧ 1 ⎫ ⎧ 1 ⎫
g (t ) = − −1
⎨ ⎬ + L
−1
⎨ ⎬
1 1
⎩s + 3⎭ ⎩ s + 1⎭
L
2 2
g (t ) = − + e−t
1 −3t 1
e
2 2
y de acuerdo a la expresión II , se tiene de la función e − 3 s el elemento a = 3 , por lo
g ( t − 3) = −
1 − 3 (t − 3) 1 − (t − 3)
e + e
2 2
Así, resulta
⎥⎦ u ( t − 3 )
⎡ 1 −3 t − 3) ⎤
II = ⎢ − e (
− 3) 1 −( t
+
⎣ 2
e
2
Finalmente
y (t ) = ⎥⎦ u ( t − 3 )
⎡ 1 −3 t − 3) ⎤
− e− 3 t + e− t + ⎢ e (
− 3) 1 − 3 (t
−
1 5 3
⎣ 2
e
3 6 2 2
181
5) Utilice la transformada de Laplace para resolver el sistema
x '' − y '' = t 2
x '' + y '' = 4 t
RESOLUCIÓN
{t }
Transformando en Laplace cada una de las ecuaciones del sistema resulta
{ x ''} − L { y ''} =L 2
s 2 X ( s ) − s X ( 0 ) − X ' ( 0 ) − ⎡⎣ s 2 Y ( s ) − s Y ( 0 ) − Y ' ( 0 ) ⎤⎦ = 3
2
s
s 2 X ( s ) − s X ( 0 ) − X ' ( 0 ) + ⎡⎣ s 2 Y ( s ) − s Y ( 0 ) − Y ' ( 0 ) ⎤⎦ = 2
4
s
enseguida se aplican condiciones iniciales y se reducen términos
s 2 X (s) − 8 s − s 2 Y (s) =
2
………………. (A)
s3
s 2 X (s) − 8 s − s 2 Y (s) =
4
…..………...….(B)
s2
Sumando (A) y (B) se obtiene
2 s 2 X ( s ) − 16 s = + 2
2 4
3
s s
de donde
2 + 4 s + 16 s 4
2s 2 X (s) = + + =
2 4
16 s
s3 s2 2s 3
X (s) =
2 + 4 s + 16 s 4
=
2 1 + 2s + 8s 4 ( )
2s 5 2s 5
X (s) = + 4 +
1 2 8
5
s s s
182
x (t ) = L { X ( s )}
y al antitransformar
−1
⎧ 1 ⎫ ⎧1⎫
x (t ) = L
⎧1⎫
⎨ 5 ⎬ + 2L ⎨ 4⎬ + 8L ⎨ ⎬
−1 −1 −1
⎩s ⎭ ⎩s ⎭ ⎩s⎭
x (t ) = t + t + 8 …………….( C )
1 4 2 3
s 2 Y (s) = + s 2 X (s) + 8 s
4
2
s
4 − s 4 X (s) + 8 s 3
s 2 Y (s) =
s2
4 − s 4 X (s) + 8 s 3
de donde
Y (s) =
s4
Y (s) = − X (s) +
4 8
4
s s
y (t ) = L {Y ( s )}
y al antitransformar
−1
⎧1 ⎫
y (t ) = 4 L ⎨ ⎬ − x (t )
⎧1⎫
⎨ 4⎬ + 8L
−1 −1
⎩s ⎭ ⎩s⎭
y (t ) = t + 8 − t − t3 − 8
4 3 1 4 2
3! 4! 3!
y (t ) = t −
2 3 1 4
t ……………….( D )
3! 4!
por lo que ( C ) y ( D ) constituyen la solución del sistema.
183
EJERCICIOS PROPUESTOS
⎧⎪0 ; t<5
f (t ) = ⎨ 2
⎪⎩t + e ; t ≥ 5
a)
∫
t
g ( t ) = e 3 t cos 5 t − t 2e − 2 t + senτ dτ
t
b)
0
{ f ( t )}
⎡ 2 ⎤
Solución.
= e −5s
⎢ 3 + 2 + + ⎥
10 25 es
⎣s s − 1⎦
a) L
s s
{g ( t )} s − 3
( )
= − +
( s − 3) 2 + 5 ( s + 2)
2 1
+ 1
b) L 3 2
s s
G (s) =
3
s + s
a) 3
, usando convolución
e −3s + s − 2
F (s) = 2
s − 4s − 5
b)
a) g ( t ) = 3 − 3 cos t
Solución.
⎧ e 2 t cosh 3 t ; 0 ≤ t < 3
⎪
f (t ) = ⎨ 2 t ⎛ 1 − 6
senh ( 3 t − 9 ) + cosh 3 t ⎟ ; t > 3
⎞
⎪e ⎜3e
b)
⎩ ⎝ ⎠
184
y '' + 4 y ' + 4 y = 3 e − t que satisface las condiciones y ( 0 ) = 2 ,
y '(0 ) = 1
y (t ) = 2 t e − 2 t − e − 2 t + 3 e − t
Solución.
Solución.
⎧ − 1 + 2 t + e − 2t
⎪ ; 0 ≤ t < 1
⎪
y (t ) = ⎨
( )
4
⎪ 1 + e e
− 2t
; t ≥ 1
2
⎪⎩ 4
x 1' = − x 2 + f ( t )
diferenciales
⎧t ; 0 ≤ t < 1
f (t ) = ⎨
x 1 = 1 − x 2' ⎩0 ; t ≥ 1
donde
185
⎧⎪ x 1 ( 0 ) = 0
⎨
⎪⎩ x 2 ( 0 ) = 0
con condiciones
( )
Solución.
⎧⎪ ; 0 ≤ t < 1
x 1 (t ) = 1 − e ( u ( t − 1) o bien x 1 ( t ) = ⎨
t − 1)
0
( t − 1)
⎪⎩− 1 + e ; t ≥ 1
y ' − 3 y = e 3 t y la condición y ( 1) = 0
6) Sea la ecuación diferencial
y = ( t − 1) e 3 t
Solución.
y (0 ) = − 1
x' + 8y = 2
x' − 2 x + y' = e 4 t
( )
Solución.
x (t ) = e 4 t 2 t − 4 t 2
186
8) Emplee la transformada de Laplace para resolver la ecuación diferencial
y '' + 6 y ' + 9 y = 6 t 2 e − 3 t
con y ( 0 ) = 0 , y ' ( 0 ) = 0 .
Solución.
y (t ) = L −1
{Y ( s )} =
1 4 −3t
2
t e
y (0 ) = 1 , y '(0 ) = 1
diferenciales con las condiciones dadas
y '' + y 2 = e − t ;
(D + 2) 2 y = e t ; y ( 0) = 1 , y '( 0) = 0
Solución.
y (t ) = e + e −2t + t e −2t
1 t 8 5
9 9 3
2 x ' + y '' = − e 2 t + δ ( t − 1)
y '(0 ) = 0
187
Solución.
x (t ) = − e 2 t + t e 2 t + e t − 1 senh ( t − 1) u ( t − 1)
1 1
2 2
y (t ) = − − t + e 2 t − t e 2 t − e t − 1 senh ( t − 1) u ( t − 1)
3 1 3
4 2 4
Solución.
y (t ) = sent − sen 3 t
3 1
8 8
− 4 x 1 + x 2 '' = − 4 x 2
x 2' ( 0 ) = − 1
Solución.
x 1( t ) = − sen 2 t +
2 3
sen 2 3 t
10 5
x 2( t ) = sen 2 t −
2 3
sen 2 3 t
5 10
188
Solución.
L { f ( t )} =
1 ⎡1 − e − s ⎤
⎢ ⎥
s ⎣1 + e − s ⎦
y '' − 7 y ' + 6 y = 10 δ ( t − 2 )
14) Aplique la transformada de Laplace, para resolver la ecuación diferencial
( )
Solución.
y ( t ) = 2 e ( ) − 2 e ( ) u ( t − 2)
6 t−2 t−2
( )
Solución.
y (t ) = − 1 + e − t + 2 t e − t u (t )
189
TEMA 6
U ( x , y) = f ( y) e x + g ( x)
RESOLUCIÓN
U ( x , y) = f ( y) e x + g ( x)
Se tiene a la función
∂ 2U
= f ( y ) e x + g '' ( x ) …………………………….( B )
∂x 2
∂U
= f ' ( y ) e x …………………………………...…….( C )
∂y
∂ 2U
= f ' ( y ) e x …………………………………...….( D )
∂ x∂ y
Se observa que al igualar las ecuaciones ( C ) y ( D ) se obtiene:
∂ 2U ∂U
=
∂ x∂ y ∂y
190
que es la ecuación en derivadas parciales buscada.
( )
función
u ( x , y) = f x 2 − 2 y + x 2
RESOLUCIÓN
( )
Para la función solución general
u ( x , y) = f x 2 − 2 y + x 2
( ) (2 x)
arbitraria.
∂u
= f ' x2 − 2y + 2 x …………………( A )
∂x
∂u
∂y
(
= f ' x2 − 2y ) (− 2 ) ………………..………( B )
( )
De ( A ) y ( B ) se despeja el factor común f ' x 2 − 2 y , esto es
(
f ' x2 − 2y = ) 1 ⎛ ∂u
⎜
2x ⎝∂x
⎞
− 2 x ⎟ ………………….( C )
⎠
(
f ' x2 − 2y = − ) 1 ∂u
2 ∂y
…………………………...( D )
igualando ( C ) y ( D ) resulta
1 ⎛ ∂u ⎞ 1 ∂u
⎜ − 2 x⎟ = −
2x ⎝∂x ⎠ 2 ∂y
191
∂u ∂u
+ x = 2x
∂x ∂y
5) Determine la ecuación en derivadas parciales cuya solución general es
u ( x , y) = f ( y) e + g ( y) e
e−x
−xy
+
1 − y2
xy
RESOLUCIÓN
y 2 u = y 2 ⎡⎣ f ( y ) e x y + g ( y ) e − x y ⎤⎦ + y 2
ex
1 − y2
………..( D )
y 2 ⎣⎡ f ( y ) e x y + g ( y ) e − x y ⎦⎤ = y 2 u − y 2
ex
1 − y2
……....( F )
∂2u
∂x2
= y 2
u +
ex
1 − y2
1 − y2 ( )
192
Finalmente
∂2u
= y2u + ex
∂x 2
193
6.2 EL MÉTODO DE SEPARACIÓN DE VARIABLES.
RESOLUCIÓN
u ( x , y) = F ( x) G ( y)
De acuerdo al método de separación de variables, se propone la función
x F '' G = F ( G − y G ' )
Posteriormente se reordenan términos para separar variables
G − y G'
=
x F ''
F G
Para la constante de separación indicada α = 0 se tiene
194
= α = 0 ……………………( B )
x F ''
F
G − y G'
también
= α = 0 ……………….( C )
G
Resolviendo la ecuación ( B ) resulta
= 0 ; x F '' = 0 ⇒ F '' = 0
x F ''
F
F = F ( x) = C 1 + C 2 x
que es una ecuación diferencial ordinaria homogénea de segundo orden, cuya solución es
G' = =
G dG G
;
y dy y
∫ ∫
Separando variables e integrando
=
dG dy
G y
ln G = ln y + C
G = G ( y) = C y
de donde resulta
u ( x , y) = A y + B x y
y al efectuar operaciones se tiene finalmente
195
RESOLUCIÓN
u ( x , t ) = F ( x ) G (t )
de acuerdo al método de separación de variables se propone como solución la función
G ( F '' − F ) = − F G '
F '' − F − G'
=
F G
Para la constante de separación α < 0 se considera α = − k 2 , por lo que se tiene
F '' − F
= α = − k 2
……………………( B )
F
− G'
= α = − k 2 ……………...…….( C )
G
Inicialmente se resuelve la ecuación ( B )
F '' − F
= − k 2
F
⇒ F '' − F + k 2 F = 0
( )
o bien
F '' − 1 − k 2 F = 0
196
que es una ecuación ordinaria homogénea de segundo orden, la cual también se puede
( )
escribir como
⎡D 2 − 1 − k 2 ⎤ F = 0
⎣ ⎦
( )
siendo su ecuación característica asociada
λ2 − 1 − k2 = 0
λ = ± 1 − k2
( )x (− )x
por lo que la solución de ( B ) es
F ( x) = C 1 e
1 − k 1 − k
+ C2e
2 2
G' − k 2G = 0
que es una ecuación de primer orden de variables separables; al resolver se tiene
G' = k 2G
= k 2G
dG
∫ ∫k
dt
=
dG 2
dt
G
ln G = k 2 t + C
y al aplicar la función exponencial resulta
G (t ) = C 3 e k
2
t
( ) (− )x⎤ C
Entonces, la solución es
⎡
u ( x , t ) = ⎢C 1 e
1 − k2 x 1 − k
+ C2e ⎥
2
2
⎣ ⎦
3 ek t
197
6.3 SERIE TRIGONOMÉTRICA DE FOURIER.
⎧− x − π ≤ x < 0
f ( x) = ⎨
⎩ x 0 ≤ x ≤ π
RESOLUCIÓN
La función está dada por dos reglas de correspondencia y su representación gráfica es:
Se trata de una función par, entonces los coeficientes de la Serie Trigonométrica de Fourier
son: a 0 , a n y b n = 0 y la serie en cuestión es de la forma
∑ a n cos
nπ x
f ( x) =
∞
a0 +
1
………………………( A )
2 n=1 L
∫ ∫
por la que se procede a la obtención de a 0 y a n .
nπ x
a0 = f ( x ) dx ; a n = f ( x ) cos
2 L 2 L
dx
L 0 L 0 L
En las expresiones anteriores L representa el semiperiodo, que en este caso está dado por
L = π ; se tiene así:
∫
π
π
2 ⎡x 2⎤
a0 = x dx = ⎢ ⎥ = π
π π ⎣ 2 ⎦
2
0
0
∫ ∫
Para el segundo coeficiente se tiene:
π
nπ x π
an = dx =
π π π
2 2
x cos x cos n x dx
0 0
198
∫
π
2⎛x ⎞
= ⎜ sen n x − sen n x dx ⎟
π ⎝n
1
⎠
an
n 0
π
2⎛x ⎞ 2 ⎛π 1 ⎞
= ⎜ sen n x + 2 cos n x ⎟ = ⎜ sen n π + 2 cos n π − 2 ⎟
π ⎝n π ⎝n
1 1
⎠ n ⎠
an
n 0
n
2⎛ 1 1 ⎞
an = ⎜ 2 ( − 1) − 2 ⎟ = 2 ⎡⎣( − 1) − 1⎤⎦
π ⎝n n π
2
n ⎠
n n
∑
π
f ( x) = ⎡( − 1) n − 1⎤ cos n x
∞
+
n π ⎣ ⎦
2
n=1
2
2
2) Obtenga los tres primeros términos no nulos del desarrollo en serie de Fourier de la
función
⎧ − 1 − π < x < 0
f ( x) = ⎨
⎩ 1 0 ≤ x < π
RESOLUCIÓN
La función indicada es una función impar por lo que sólo se tienen coeficientes b n ,
entonces el desarrollo de la serie trigonométrica de Fourier es de la forma
∑ b n sen
nπ x
f ( x) =
∞
n=1 L
∫
Para obtener el coeficiente b n se considera
nπ x
bn = f ( x ) sen
2 L
dx
L 0 L
y para L = π
∫
resulta
π
bn =
π
2
sen x dx
0
199
enseguida se integra y evalúa la integral
π
⎡1 − ( − 1) n ⎤
2⎡ 1 ⎤
= ⎢ − cos n x ⎥ =
π ⎣ n nπ ⎣ ⎦
2
⎦
bn
0
f ( x) = ∑ ⎡1 − ( − 1) n ⎤ sen ( n x )
∞
nπ ⎣ ⎦
2
n=1
Es necesario obtener los 3 primeros términos no nulos de la serie indicada, por lo que
dando a n distintos valores se observa que
b n = 0 , si n es par
b n ≠ 0 , si n es impar
Por lo anterior, se tiene que para la función con términos no nulos se considerará
n = 1 , 3 , 5 ; esto es,
RESOLUCIÓN
∑ a n cos
nπ x
f ( x) =
∞
a0 +
1
2 n=1 L
200
∫ ∫
nπ x
a0 = f ( x ) dx ; a n = f ( x ) cos
2 L 2 L
dx
L 0 L 0 L
∫
Al sustituir L y la función dada se tiene por un lado
⎡x 3⎤
a0 = 2 x dx = 2 ⎢ ⎥ =
1
1
2
⎣ 3 ⎦
2
0 3
0
∫
y por otro
an = 2 x 2 cos n π x dx
1
⎡ 2x ⎛ x2 2 ⎞ ⎤
= 2 ⎢ 2 2 cos n π x + ⎜ − 3 3 ⎟ sen n π x ⎥
1
⎣⎢ n π ⎝nπ n π ⎠ ⎥⎦
an
0
⎡ 2 ⎛ 1 2 ⎞ ⎤
an = 2⎢ 2 cos n π + ⎜ − 3 3 ⎟ sen n π − 0 ⎥
⎣n π ⎝nπ n π ⎠ ⎦
2
= ( − 1) n
n π
4
2 2
f ( x) = ∑ ( − 1) n cos n π
∞
+
n π
1 4
x
n=1
2 2
3
y para los cuatro primeros términos no nulos damos a n los valores 1 , 2 y 3 , pues el
primer término no nulo está dado por a 0 ; así se tiene
o bien
201
RESOLUCIÓN
Es conveniente trazar la gráfica de la función para determinar a partir de ella qué tipo de
función (par o impar) es.
Se observa que es una función par (simétrica respecto al eje de las ordenadas) por lo que el
desarrollo de la serie de Fourier estará dado por
∑ a n cos
nπ x
f ( x) =
∞
a0 +
1
2 n=1 L
∫ ∫
Los coeficientes a 0 y a n se determinan a partir de las expresiones
nπ x
a0 = f ( x ) dx ; a n = f ( x ) cos
2 L 2 L
dx
L 0 L 0 L
El periodo de la función es T = 2 , por lo que el semiperiodo es L = 1 ; entonces al
∫
sustituir en las expresiones anteriores se tiene
⎡ x2 ⎤
a0 = 2 ( − x + 1) dx = 2 ⎢ − + x⎥
1
1
0
⎣ 2 ⎦ 0
⎛ 1 ⎞
a 0 = 2⎜− + 1⎟ = 1
⎝ 2 ⎠
an = 2
∫ (− x + 1) cos n π x dx
1
∫ ∫
de aquí resultan dos integrales
⎡ ⎤
a n = 2 ⎢− x cos n π x dx + cos n π x dx ⎥
1 1
⎣ 0 0 ⎦
resolviendo
202
⎡ ⎛ cos n π x π ⎞ ⎛ sen n π x ⎞ ⎤
= 2 ⎢− ⎜ 2 2 + ⎟ + ⎜ ⎟ ⎥
1 1
⎢ ⎝ n π nπ π
x sen n x
⎠0 ⎝ ⎠ 0 ⎥⎦
an
⎣
n
⎡ − ( − 1) n + 1 ⎤
an = 2⎢ ⎥ = − 2 ⎡( − 1) n − 1⎤
n π n π ⎣ ⎦
2
⎢⎣ ⎥⎦
2 2 2
f ( x) = ∑ cos n π x
∞
+
1
2 n=1
RESOLUCIÓN
π
Se tiene una función definida en medio intervalo, entonces L =
2
El desarrollo en serie seno de Fourier está dado por
∑b
nπ x
f ( x) =
∞
sen
n =1
n
L
203
∫
nπ x
bn = f ( x ) sen
2 L
dx
L 0 L
∫
Entonces se tiene
nπ x
π
bn = cos ( 3 x ) sen
π π
2 2
dx
0
π ∫
2 2
π
bn = cos ( 3 x ) sen ( 2 n x ) dx
4 2
0
∫
1 ⎡ cos ( m + n ) u cos ( m − n ) u ⎤
argumento diferente, que se puede resolver empleando la expresión
sen m u cos n u du = − ⎢ − ⎥ + C
2⎣ m + n m − n ⎦
Así, para la integral que nos ocupa
4 ⎡ cos ( 2 n + 3 ) x cos ( 2 n − 3 ) x ⎤
π
bn = ⎢− − ⎥
2 ( 2 n + 3) 2 ( 2 n − 3) ⎦
2
π ⎣
0
4 ⎡ 2n ⎤
( )
bn = ⎢ ⎥ =
π ⎣4 n − 9⎦
8n
π 4n2 − 9
2
f ( x) = ∑ sen ( 2 n x )
( )
∞
8n
n=1 π 4n2 − 9
RESOLUCIÓN
Para la función en cuestión el desarrollo de la serie de cosenos de Fourier está dado por
204
∑ a n cos
nπ x
f ( x) =
∞
a0 +
1
2 n=1 L
∫ ∫
y los coeficientes a 0 y a n se determinan a partir de las expresiones
nπ x
a0 = f ( x ) dx y a n = f ( x ) cos
2 L 2 L
dx
L 0 L 0 L
Se tiene una función definida en medio intervalo donde L = 1 , de donde
a0 = 2
∫ e x dx = 2 ⎣⎡ e x ⎦⎤ = 2 ⎣⎡e − e 0 ⎦⎤ = 2 ( e − 1)
1 1
0 0
∫
también
⎡ e x nπ ⎛ ⎞⎤
an = 2 e cos n π x dx = 2 ⎢ 2 2 ⎜ sen n π x + cos n π
1
x ⎟⎥
1
⎣n π + 1⎝ nπ
1
⎠⎦
x
0
0
⎡ enπ ⎛ ⎞ nπ ⎛ 1 ⎞⎤
an = 2⎢ 2 2 ⎜ sen n π + cos n π ⎟ − 2 2 ⎜ ⎟⎥
⎣n π + 1⎝ nπ n π + 1⎝ n π
1
⎠ ⎠⎦
n = 1 , 2 , 3 ,…
por lo que resulta
⎡ nπ e ⎛ 1 ⎞ ⎤
an = 2⎢ 2 2 ⎜ ( − 1) n ⎟ − 2 2 ⎥
⎣n π + 1 ⎝ n π n π + 1⎦
1
⎠
simplificando
⎡ ( − 1) n e − 1 ⎤
= 2⎢ 2 2 ⎥
⎣⎢ n π + 1 ⎦⎥
an
Finalmente se obtiene
⎡ ( − 1) n e − 1 ⎤
f ( x ) = ( e − 1) + 2 ∑
∞
⎢ 2 2 ⎥ cos n π x
n=1 ⎢⎣ n π + 1 ⎥⎦
205
6.4 RESOLUCIÓN DE PROBLEMAS DE CONDICIONES INICIALES Y
DE FRONTERA.
RESOLUCIÓN
u ( x , y) = F ( x)G ( y)
como solución de la ecuación en derivadas parciales, esto es
∂2u
= F ( x ) + G '' ( y )
∂y 2
206
A continuación se resuelve cada una de las ecuaciones diferenciales ordinarias obtenidas.
Por un lado se tiene
F' = α F
y al resolver
= αF
dF
∫ ∫ αdx
dx
=
dF
F
ln F = α x + C
F ( x) = C 1 e α x
y por otro
G '' + G = α G '
G '' − α G ' + G = 0
(D )
que en términos del operador diferencial es
2
− αD + 1 G = 0
α + α − 4 α − α − 4
⇒ λ1 = , λ2 =
2 2
2 2
por lo que
⎛α + α − 4⎞ ⎛α − α − 4⎞
G ( y) = C 2 e
⎜ ⎟y ⎜ ⎟y
2 2
⎜ ⎟ ⎜ ⎟
⎝ ⎠
+ C3e ⎝ ⎠
2 2
⎜ ⎟ ⎜
⎢C 2 e
⎝ ⎠
+ C3e ⎝ ⎠
⎥
2 2
⎢ ⎥
⎣ ⎦
que al efectuar operaciones se simplifica según se indica
⎛α + α − 4⎞ ⎛α − α − 4⎞
u ( x , y) = A e
⎜ ⎟ y +α x ⎜ ⎟ y +α x
2 2
⎜ ⎟ ⎜ ⎟
⎝ ⎠
+ Be ⎝ ⎠
2 2
……………………( A )
207
Posteriormente se aplica la condición de frontera u ( 0 , y ) = 2 e y , esto es
⎛α + α − 4⎞ ⎛α − α − 4⎞
⎜ ⎟y ⎜ ⎟y
2 2
⎜ ⎟ ⎜ ⎟
⎝ ⎠
+ Be ⎝ ⎠
= 2 e y ……………………..( B )
2 2
Ae
y por igualdad de exponentes
α + α − 4
= 1
2
α − α − 4
= 1
2
2
o bien
α ± α − 4
= 1
2
α ± α 2
− 4 = 2
± α 2
− 4 = 2 − α
α − 4 = (2 − α )
2 2
α 2
− 4 = 4 − 4α + α 2
4α = 8
α =2
sustituyendo este valor en ( B )
⎛2 + 4 − 4⎞ ⎛2 − 4 − 4⎞
⎜⎜ ⎟⎟ y ⎜⎜ ⎟⎟ y
⎝ ⎠
+ Be ⎝ ⎠
= 2e y
2 2
Ae
de donde
Ae y
+ Be y
= 2e y
(A + B) e y
= 2e y
(A + B ) = 2 ……………………..( C )
y por igualdad de coeficientes
208
Al sustituir α en u ( x , y) = A e y + 2x
+ Bey + 2x
= ( A + B) e y + 2 x
(A + B) e y + 2 x = 2 e y + 2 x
u ( x , y) = 2 e y
Finalmente
+ 2x
∂ 2u ∂ 2u
= 2
∂y 2 ∂x ∂ y
RESOLUCIÓN
Se tiene la ecuación
∂2u ∂ 2u
=
∂y 2 ∂x ∂ y
2 …………………..( I )
con ciertas condiciones de frontera, las cuales llevarán a la obtención de la solución que
satisfaga dichas condiciones.
u ( x , y ) = F ( x ) G ( y ) ……………………( II )
∂2u
= F 'G'
∂x ∂ y
sustituyendo en la ecuación ( I )
F G '' = 2 F ' G '
209
separando variables resulta
=
G '' F'
2G' F
y para la constante de separación α
= α = α
G '' F'
;
2G' F
G '' − 2 α G ' = 0
cuya solución es
G ( y) = C 1 + C 2 e 2α y
y por otro
F' = α F
= αF
dF
∫ ∫ α dx
dx
=
dF
F
ln F = α x + C
F ( x) = C 3 e α x
( )
por lo que la solución completa es
u ( x , y) = C 3 e α x C 1 + C 2 e 2α y
u ( x , y) = A e α x + B e (
α x + 2 y)
u ( x , y) = A 1 e α 1 x + B 1 e
α 1 ( x + 2 y) ( x + 2 y)
+ A2 e + B2e
α2x α 2
……( III )
210
5 + e 2 y + 2 e 4 y = A 1 + B 1 e 2α 1 y + A 2 + B 2 e 2α 2 y
donde por igualdad de coeficientes se tiene
A1 + A2 = 5
B1 = 1
B2 = 2
2α 1 = 2 ⇒ α1 = 1
2α 2 = 4 ⇒ α 2 = 2
sustituyendo estos valores en ( III )
u ( x , y) = A 1 e x + e x + 2 y + A 2 e 2 x + 2 e (
2 x + 2 y)
…………..( IV )
3e x + 5 e 2x = A1 e x + e x + A2 e 2x + 2e 2x
A2 + 2 = 5 ⇒ A2 = 3
Finalmente, al sustituir en ( IV ) resulta
u ( x , y) = 2 e x + e x + 2 y + 3 e 2 x + 2 e (
2 x + 2 y)
U (0 , t ) = 0
t xx
U (L , t) = 0
⎛π x⎞ ⎛ 3π x ⎞
U ( x , 0 ) = sen ⎜ ⎟ + sen ⎜ ⎟
1
⎝ L ⎠ 2 ⎝ L ⎠
211
RESOLUCIÓN
U (0 , t ) = 0
se tienen las condiciones de frontera
U (L , t) = 0
⎛π x⎞ ⎛ 3π x ⎞
U ( x , 0 ) = sen ⎜ ⎟ + sen ⎜ ⎟
1
⎝ L ⎠ 2 ⎝ L ⎠
y se considera una constante de separación negativa, es decir se considerará de la forma
− k2 , k ∈
U ( x , t ) = F ( x ) G ( t ) …………….………..( I )
De acuerdo al método de separación de variables se tiene
∂ 2U
= F '' ( x ) G ( t )
∂x 2
enseguida se sustituye en la ecuación diferencial parcial dada:
F ( x ) G ' ( t ) = α 2 F '' ( x ) G ( t )
o bien
F G ' = α 2 F '' G
Separando variables
= 2
α G
F '' G'
F
212
para la constante de separación igual a − k 2 se tiene
= − k2 , = − k2
α G
F '' G'
2
F
Resolviendo para cada una de las ecuaciones diferenciales anteriores
= − k2 ⇒ F '' + k 2 F = 0
F ''
F
= − k2 ⇒ G' + k 2 α 2 G = 0
α G
G'
2
F ( x ) = C 1 cos k x + C 2 sen k x
se obtiene, en cada caso
G (t ) = C 3 e − k α2t
2
U (0 , t ) = F (0 ) G (t ) = 0
F (0 ) = C 1 = 0
F ( x ) = C 2 sen k x
entonces
U (L , t) = 0
Para la condición de frontera
U ( L , t ) = F ( L ) G (t ) = 0
F ( L ) = C 2 sen k L = 0
k L = ang sen ( 0 )
resulta
sen k L = 0 ⇒
k L = nπ para n = 1 , 2 ,…
nπ
k =
L
por lo que
213
nπ
F ( x ) = C 2 sen x ; n = 1 , 2 ,…
L
Sustituyendo en ( I )
⎛ ⎞
nπ ⎞⎜
⎛nπ ⎞
U ( x , t ) = ⎜ C 2 sen
⎛ −⎜ ⎟ α
⎟
2
x⎟ C 3 e ⎝ L ⎠
2
t
⎝ L ⎠⎜ ⎟
⎝ ⎠
que se puede expresar también como
∑ Cne
nπ x
⎛nπ α ⎞
U (x , t) =
∞ −⎜ ⎟ t
2
⎝ L ⎠
sen
n = 1 L
⎛π x⎞ ⎛ 3π x ⎞
Para la condición de frontera U ( x , 0 ) = sen ⎜ ⎟ + sen ⎜ ⎟
1
⎝ L ⎠ 2 ⎝ L ⎠
∑ C n sen
⎛π x ⎞ ⎛ 3π x ⎞ ∞
nπ x
sen ⎜ ⎟ + sen ⎜ ⎟ =
1
⎝ L ⎠ 2 ⎝ L ⎠ n = 1 L
n = 1 , C1 = 1 ; n = 3 , C3 =
1
2
entonces la solución es
⎛π x ⎞ ⎛ 3π x ⎞
⎛π α ⎞ ⎛ 3π α ⎞
U (x , t) = e
−⎜ ⎟ t −⎜ ⎟ t
⎟ + e
2 2
⎝ L ⎠
sen ⎜ ⎝ L ⎠
sen ⎜ ⎟
1
⎝ L ⎠ 2 ⎝ L ⎠
RESOLUCIÓN
214
Las derivadas parciales para la función propuesta como solución U ( x , y ) = F ( x ) G ( y ) ,
son
∂2u ∂u ∂ 2u
= = = F 'G'
∂x 2 ∂x ∂ x∂ y
F '' G , F ' G ,
G' − G
=
F ''
G' − G
= α ; = α
F ''
F' G
(D )
enseguida se resuelve cada una de estas ecuaciones diferenciales ordinarias:
F '' − α F ' = 0 2
− αD F = 0
λ (λ − α ) = 0
;
λ2 − αλ = 0 ⇒ ∴ λ1 = 0 , λ2 = α
F ( x) = C 1 + C 2 e α x
de donde
también
G' − G = α G = G ( α + 1)
dG
;
∫ ∫
dy
= (α + 1) d y
dG
ln G = (α + 1) y
G
por lo que
G ( y) = C 3 e ( )
1+α y
( )
así, la solución completa es
u ( x , y) = C 1 + C 2 e α x C 3 e ( )
1+α y
o bien
215
u ( x , y) = A e (
1+α) y α x + (1 + α ) y
+ Be …………………( I )
2 = Ae(
1+α) y
+ Be(
1+α) y
2 = ( A + B) e (
1+α) y
(1 + α ) y
de aquí se tiene
= 0
1+α = 0 ⇒ α = − 1
A + B = 2
A = 2 − B
sustituyendo en ( I )
u ( x , y) = (2 − B) e 0 + B e − x + 0
u ( x , y) = (2 − B) + B e − x
para la condición u ( − 1 , y ) = 0
0 = (2 − B) + B e
B (1 − e ) = 2
B = ⇒ A = 2 −
2 2
1 − e 1 − e
;
2 − 2e − 2
A = = −
2e
1 − e 1 − e
Finalmente se sustituyen los valores de A , B y α en ( I )
u ( x , y) = − e−x
2e 2
e − 1 e − 1
( )
y luego de simplificar
u ( x , y) = e − e−x
2
e − 1
216
EJERCICIOS PROPUESTOS
∂u ∂ 2U
= 2 ; 0 < x < 3 , t > 0
∂t ∂x2
( )
sabiendo que una solución completa de la ecuación en derivadas parciales es
U ( x , t ) = ( A cos k x + B sen k x ) e − 2 k
2
t
Solución.
⎡ 1 2( 2 n − 1) 2 π
(2 n − 1) π ⎤
u (x , t) = ∑⎢
∞ −
x⎥
2
π n =1⎢2 n − 1
8 t
⎥⎦
9
⎣
e sen
3
− π ≤ x ≤ π.
Solución.
∑
π
f ( x) = 1 + ⎡( − 1) n − 1⎤ cos n x
∞
+
n π ⎣ ⎦
2
n=1
2
2
∂2u
= u sujeta a u ( 0 , t ) = 2 e t + 4 e − t
∂ x ∂t
3) Resuelva la ecuación
217
4) Represente a la función f ( x ) = ( x − π ) 2 , definida en π < x < 2 π , por medio
Solución.
∑
π ⎛ 1 π ⎞
f ( x) =
∞
+ ⎜ 2 cos 2 n x − sen 2 n x ⎟
2
3 n =1 ⎝n n ⎠
π
el intervalo 0 < x <
2
Solución.
f ( x) = ∑ (16 n − 9) π (16 n 2 − 9 ) π
∞
+ cos 4 n x +
3π
4 12 16 n
sen 4 n x
n=1
2
t ∈ ( 0 , 2) .
Solución.
∑
⎛ 4 ⎞
f (t ) = cos ( n π t ) − sen ( n π t ) ⎟
∞
+ ⎜ 2
⎝n π nπ
4 4
n=1 ⎠
2
3
∂u ∂2u
= α ; 0 < x < 100 , t > 0
∂t ∂x2
2
218
u (0 , t ) = 0 , t > 0
u ( 100 , t ) = 0 , t > 0
u ( x , 0 ) = 50
con
, 0 < x < 100
Solución.
( 2 n − 1) π x e − n10000
u(x , t) = ∑
∞ π
α2t
( 2 n − 1) π
2 2
200
sen
n=1 100
Solución.
9) Obtenga los cuatro primeros términos del desarrollo en serie coseno de Fourier de la
función mostrada en la figura siguiente:
219
Solución.
⎛ 16 4 ⎞ π ⎛ 4 ⎞ ⎛ 16 4 ⎞ 3π x
f ( x) = + ⎜ − 3 + 2 ⎟ cos x + ⎜ − 2 ⎟ cos x + ⎜ + 2 ⎟
+ …
⎝ π π ⎠ ⎝ π ⎠ ⎝ 27 π 9π ⎠
5
3
cos
12 2 2
el intervalo − π ≤ x ≤ π .
Sugerencia: Verifique si la función es par o impar.
( − 1) n
Solución.
f ( x ) = − 12 A ∑
∞
sen n x
n=1 n3
220