UNIVERSIDAD AUSTRAL DE CHILE
FACULTAD DE CIENCIAS DE LA INGENIERIA
CENTRO DE DOCENCIA DE CIENCIAS BASICAS
PARA INGENIERIA.
BAIN041 ECUACIONES DIFERENCIALES PARA INGENIERIA
Prueba Parcial N 2 1
Lunes 27 de octubre de 2014
Nombre:............................................................................Grupo.........................Sala.........
Conteste en forma ordenada identificando la pregunta e tem que corresponde.
Cada soluci
on debe llevar desarrollo y respuesta.
1.- (10 pts.) ........................
Aproxime a tres decimales e indique donde aproxima.
2.- (10 pts.) ........................
Est
a permitido el uso de calculadora estandar.
3.- (20 pts.) ........................
No est
a permitido el uso de celular.
4.- (20 pts.) ........................
Tiempo: 120 minutos.
Problema 1 ESCRIBA (sin resolver) la forma de una soluci
on particular de la ecuaci
on diferencial
y (4) y = xex + sin(x) .
Soluci
on:
La ecuacion caracterstica asociada a la EDO homogenea es:
m4 1 = 0
(m2 1)(m2 + 1) = 0
cuyas soluciones son m1 = 1, m2 = 1, m3 = i y m4 = i.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Luego, como 1 es raz simple, el metodo de coeficientes indeterminados implica que la parte de xex tiene
asociada una soluci
on particular de la forma (Ax + B)xex .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Por otro lado, como i es raz simple, el metodo de coeficientes indeterminados implica que la parte de
sin(x) tiene asociada una soluci
on de la forma (C cos(x) + D sin(x))x.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
As, por superposici
on se concluye que una solucion particular tiene la forma
yp (x) = (Ax + B)ex x + (C cos(x) + D sin(x))x , x R ,
con A, B, C y D coeficientes reales por determinar.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Problema 2 Encuentre la soluci
on general de la siguiente ecuaci
on diferencial:
x2 y 00 xy 0 + y = x(x + 1) , x > 0 .
Soluci
on: Sabemos que la soluci
on general y se escribe como
y(x) = yH (x) + yp (x) , x > 0 ,
donde yH denota la soluci
on general de la EDO homogenea asociada e yp una solucion particular.
1
CAG/SJC/AMR
Soluci
on homog
enea
La EDO es tipo Cauchy-Euler. Entonces
y(x) = xm ,
y 0 (x) = mxm1
y 00 (x) = m(m 1)xm2
Reemplazando en la EDO homogenea se tiene:
m(m 1) m + 1 = 0
m2 2m + 1 = 0
(m 1)2 = 0
Las races de la ecuaci
on son m1 = m2 = 1 o m = 1 con multiplicidad 2.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Por tanto, la soluci
on general de la EDO homogenea es
yH (x) = C1 x + C2 x ln(x);
C1 , C1 R , x > 0 .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Soluci
on particular
Usando variaci
on de par
ametros, una solucion particular se escribe de la forma
yp (x) = v1 (x)y1 (x) + v2 (x)y2 (x)
Con y1 (x) = x e y2 (x) = x ln(x). Ademas, R(x) =
x
x ln(x)
W (y1 , y2 )(x) =
1 1 + ln(x)
As
W1 =
= x(1 + ln(x)) x ln(x) = x
x ln(x)
= x ln(x)
1 + ln(x)
0
x+1
x
x+1
. Entonces,
x
x+1
x
= (x + 1) ln(x)
W =x.
W1 = (x + 1) ln(x)
Luego,
Z
v1 (x) =
W1
dx =
W
(x + 1) ln(x)
dx =
x
Z
ln(x)
ln2 (x)
ln(x) +
dx = x ln(x) + x
x
2
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 punto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Por otro lado,
x
W2 =
1
x+1
x
=x x+1 =x+1
x
W2 = x + 1
As,
v2 (x) =
W2
dx =
W
x+1
dx =
x
Z
1
1+
dx = x + ln(|x|) = x + ln(x)
x
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 punto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Entonces una soluci
on particular para la EDO es
yp (x) = v1 (x)y1 (x) + v2 (x)y2 (x)
ln2 (x)
x + (x + ln(x))x ln(x)
=
x ln(x) + x
2
x(ln(x))2
+ x2 ln(x) + x(ln(x))2
= x2 ln(x) + x2
2
= x2 + 21 x(ln(x))2
.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Soluci
on general
La solucion general de la EDO es
1
y(x) = yH (x) + yp (x) = C1 x + C2 x ln(x) + x2 + x(ln(x))2 ;
2
C1 , C1 R
, x > 0
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Problema 3
Resuelva el siguiente problema de valores iniciales
( 00
y (t) + 4y 0 (t) + 4y(t) = 4 cos(2t) + U(t 1) sin(t 1) , t 0
.
y(0) = 0 , y 0 (0) = 1
Soluci
on : Usando propiedades de la transformada de Laplace y transformada de Laplace de funciones
basicas, se obtiene
L y 00 (s) + 4L y 0 (s) + 4L(y)(s) = 4L cos(2t) (s) + L (t 1) sin(t 1) (s)
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
s2 L(y)(s) 1 + 4sL(y)(s) + 4L(y)(s) =
4s
s2 +4
s
+ se2 +1
(s + 2)2 L(y)(s) =
4s
s2 +4
s
+ 1 + se2 +1
(s + 2)2 L(y)(s) =
4s+s2 +4
s2 +4
L(y)(s) =
1
s2 +4
s
+ se2 +1 ,
o equivalentemente ,
es
(s2 +1)(s+2)2
para todo s suficientemente grande.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Por otro lado, la descomposici
on en fracciones parciales asociada a
1
(s2 +1)(s+2)2
As+B
s2 +1
C
s+2
D
(s+2)2
1
(s2 +1)(s+2)2
queda
1 = (As + B)(s + 2)2 + C(1 + s2 )(s + 2) + D(1 + s2 )
As, considerando la extensi
on en C, al evaluar en s = i (o bien considerando la igualdad coeficiente a
coeficiente en las respectivas potencias) se concluye
1 = (Ai + B)(i + 2)2 3A + 4B = 0 ; 4A + 3B = 1 . Luego ,
A=
4
3
, B=
.
25
25
Para los coeficientes restantes, evaluando en s = 2 se obtiene D = 51 mientras que al evaluar en s = 1
4
se obtiene C = 25
.
Luego, usando la descomposici
on en fracciones parciales anterior y transformada de funciones conocidas
se obtiene
4s 3
+
4
1
L(y)(s) = L 12 sin(2t) (s) + es s252 +125 + 25(s+2)
+ 5(s+2)
2
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4
1
L(y)(s) = L 12 sin(2t) (s) + es 25(s4s2 +1) + 25(s32 +1) + 25(s+2)
+ 5(s+2)
2
2t
2t
= L 12 sin(2t) (s) + es L 4 cos(t)
+ 3 sin(t)
+ 4e25 + te5 (s)
25
25
de donde se concluye a partir de la segunda formula de traslacion
sin(2t)
y(t) =
+ (t 1)
2
4 cos(t 1) 3 sin(t 1) 4e2(t1) (t 1)e2(t1)
+
+
+
25
25
25
5
!
, t 0.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Problema 4 Responda UNO Y SOLO UNO de los siguientes problemas:
a) Un cuerpo de masa 1 Kg se sujeta a un resorte suspendido del techo, el cual alcanza su posici
on
de equilibrio al estirarse 1 m. Tal sistema es sometido a amortiguacion, cuya constante asociada es
7 N s/m, y a la acci
on de una fuerza externa dada por f (t) = 100 sin(10t) N .
En el instante t = 0, el cuerpo se desplaza 2 m bajo su posicion de equilibrio y se libera con velocidad
inicial nula. Si se considera g = 10 m/s2 :
i) Establezca el PVI que modela el movimiento de la masa.
Soluci
on:
Definici
on de variable
x(t) : Posici
on del cuerpo de masa m en el tiempo t.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Planteamiento PVI
Antes de plantear el PVI, necesitamos la constante del resorte
mg = ks k =
10 1
N
k = 10
1
m
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Luego, el PVI es
00
x + 7x0 + 10x = 100 sin(10t)
x(0) = 2, x0 (0) = 0
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
ii) Encuentre la ecuaci
on de movimiento.
Soluci
on:
Soluci
on homog
enea
La ecuaci
on caracterstica asociada a la EDO homogenea es:
2 + 7 + 10 = 0
Sus races son: 1 = 5 y 2 = 2
Por lo tanto, la soluci
on general de la EDO homogenea es:
xh (t) = C1 e5t + C2 e2t ; C1 , C2 R , 0 t
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Soluci
on particular
xp (t) = A cos(10t) + B sin(10t)
x0p (t) = 10A sin(10t) + 10B cos(10t)
x00p (t) = 100A cos(10t) 100B sin(10t)
Reemplazando en la EDO original, se tiene:
(90A + 70B) cos(10t) + (70A 90B) sin(10t) = 100 sin(10t)
El sistema de ecuaciones lineales es
9A + 7B = 0
,
7A 9B = 10
de donde se obtiene A =
7
9
yB= .
13
13
Luego, una soluci
on particular de la EDO es
xp (t) =
9
7
cos(10t)
sin(10t) , 0 t .
13
13
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Soluci
on general
x(t) = xh (t) + xp (t)
x(t) = C1 e5t + C2 e2t
7
9
cos(10t)
sin(10t) ; C1 , C2 R , 0 t
13
13
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Ahora debemos determinar las constantes C1 y C2 .
x0 (t) = 5C1 e5t 2C2 e2t +
7 10
9 10
sin(10t)
cos(10t)
13
13
Utilizando las condiciones iniciales, se obtiene
x(0) = 2 = C1 + C2
7
13
x0 (0) = 0 = 5C1 2C2
90
13
Resolviendo el sistema de ecuaciones lineales anterior se obtiene C1 =
52
85
= 4 y C2 = .
13
13
As, la soluci
on del PVI es
x(t) = 4e5t +
7
9
85 2t
e
cos(10t)
sin(10t) , 0 t .
13
13
13
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
iii) Cu
al es la soluci
on de estado estable o estacionario?
Soluci
on:
De la parte anterior, lm
t+
4e5t
85
+ e2t
13
= 0, de donde se concluye que el termino de
estado estable de la soluci
on es
xestable (t) =
7
9
cos(10t)
sin(10t)
13
13
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
iv) Determine la velocidad de la masa en cualquier instante de tiempo t.
Soluci
on: Derivando la soluci
on se obtiene directamente
170 2t 70
90
e
+
sin(10t)
cos(10t) , 0 t .
13
13
13
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
v(t) = x0 (t) = 20e5t
b) Un circuito RCL en serie, tiene una fueza electromotriz dada por E(t) = 100 sin(10t) V , un resistor
de 7 , un inductor de 1 H y un capacitor de 0.1 F . Si la corriente inicial y la carga inicial son nulas:
i) Establezca el PVI que modela la carga.
Soluci
on:
Definici
on de variable
q(t) : Cantidad de la carga electrica del circuito en el tiempo t.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Planteamiento PVI
Los datos son: E(t) = 100 sin(10t), L = 1, R = 7, C = 0.1, q(0) = 0 y I(0) = 0. Entonces, el
PVI es
00
q + 7q 0 + 10q = 100 sin(10t)
q 0 (0) = q(0) = 0
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
ii) Encuentre la ecuaci
on de la carga del circuito.
Soluci
on:
Soluci
on homog
enea
La ecuaci
on caracterstica asociada a la EDO homogenea es:
2 + 7 + 10 = 0
Sus races son: 1 = 5 y 2 = 2
Por lo tanto, la soluci
on de la EDO homogenea es:
qh (t) = C1 e5t + C2 e2t ; C1 , C2 R , 0 t .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Soluci
on particular
qp (t) = A cos(10t) + B sin(10t)
qp0 (t) = 10A sin(10t) + 10B cos(10t)
.
qp00 (t) = 100A cos(10t) 100B sin(10t)
Reemplazando en la EDO original, se tiene
(90A + 70B) cos(10t) + (70A 90B) sin(10t) = 100 sin(10t)
El sistema de ecuaciones lineales es
9A + 7B = 0
,
7A 9B = 10
de donde se obtiene A =
7
9
yB= .
13
13
As, una soluci
on particular de la EDO es
qp (t) =
9
7
cos(10t)
sin(10t) , 0 t .
13
13
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Soluci
on general
q(t) = qh (t) + qp (t)
q(t) = C1 e5t + C2 e2t
7
9
cos(10t)
sin(10t) ; C1 , C2 R , 0 t .
13
13
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Ahora debemos determinar las constantes C1 y C2 .
q 0 (t) = 5C1 e5t 2C2 e2t +
7 10
9 10
sin(10t)
cos(10t)
13
13
Utilizando las condiciones iniciales, se obtiene
q(0) = 0 = C1 + C2
7
13
q 0 (0) = 0 = 5C1 2C2
Resolviendo el sistema anterior se llega a C1 =
90
.
13
8
125
y C2 =
.
3
39
Por lo tanto, la soluci
on del PVI es
8
125 2t
7
9
q(t) = e5t +
e
cos(10t)
sin(10t) , 0 t .
3
39
13
13
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
iii) Cu
al es la soluci
on de estado estable o estacionario?
Soluci
on:
8
125 2t
De la parte anterior, lm e5t +
e
t+
3
39
estado estable de la soluci
on es
qestable (t) =
= 0, de donde se concluye que el termino de
7
9
cos(10t)
sin(10t) .
13
13
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
iv) Determine la corriente electrica del circuito en cualquier instante de tiempo t.
Soluci
on:
Derivando la soluci
on se obtiene directamente
90
40 5t 250 2t 70
e
e
+
sin(10t)
cos(10t) , 0 t .
3
13
13
13
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 puntos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
i(t) = q 0 (t) =