Ecuaciones Diferenciales - APLICACIONES
Ecuaciones Diferenciales - APLICACIONES
DIFERENCIALES
LINEALES
6
DE PRIMER ORDEN
M É T O D O DE SO L U C IÓ N
y p ( x ) y - q(x) (6.1)
I( x ) = e ]* x)i' {6.2)
/ (x ) y + p (x )I (x )y = l(x )q (x ) (6.3)
^ = /* ( x )
dx
integrar a m b o s la d o s d e esta ú ltim a e c u a c ió n c o n resp ecto a jc, y lu e g o reso lv e r para y la e c u a c ió n resu ltan te.
R E D U C C IÓ N D E E C U A C IO N E S DE B E R N O U L L I
/ + P ( x ) y = q ( x ) y" (6.4)
Z = y1“" (6.5)
42
www.FreeLibros.me
P roblem as resuelto s 43
PR O BLEM A S R ESUELTOS
La ecuación diferencial tiene la forma de la ecuación (6.1), con p(x) = - 3 y q(x) = 6, y es lineal. Aquí
/ p ( x )d x = f —’i d x = - 3 x
I ( x ) = e ¡i*x)dx = e - 3x ( 1)
Multiplicando la ecuación por el factor de integración definido por (1) del problema 6.1, obtenemos
f — (ye~3x)dx = f 6e~3xdx
J dx J
ye~3x = - 2 e ~ 3x + c
y = ce3x - 2
La ecuación diferencial tiene la forma de la ecuación (6.1), con p(x) = -2 x y q(x) = x, y es lineal. Aquí
fp(x)dx = f ( - 2 x ) d x = - x 2
I ( x ) = e ^ p<-X)ix =■«-** ( 1)
Multiplicando la ecuación diferencial por el factor de integración definido por (1) del problema 6.3, obtenemos
Integrando ambos lados de esta última ecuación con respecto a x, encontramos que
d .
f — (ye~x> )dx = f xe~x’dx
J dx J
-x ’ 1 s ,
ye = - —e +c
y = c e -----1
. 2
La ecuación diferencial tiene la forma de la ecuación (6.1), con p ( x ) — 4 /x y q(x) = x 4, y es lineal. Aquí
J p ( x ) d x = J '^ á x = 41n|x| = ln x 4
www.FreeLibros.me
44 C a p ít u l o 6 E c u a c io n e s d i f e r e n c ia l e s l in e a l e s d e p r i m e r orden
Multiplicando la ecuación diferencial por el factor de integración definido por (7) del problema 6.5, obtenemos
4 1 9 , u* C 1 j
yx = —x + c ob ten y = —r + —x
9 x4 9
6 .7 . R esu elv a / + y = se n x.
Aquí p(x) = 1; por lo tanto 7(x) = e^ldx = ex. Multiplicando la ecuación diferencial por /(x), obtenemos
Integrando ambos lados de la ecuación con respecto a x (para integrar el lado derecho, usamos dos veces integración por
partes), encontramos que
y = ce -* + —
1 1
s e n x ----- cosx
2 2
Aquí, p (x ) = —5 y 7(x) = e ^ ~ !>dx = e~,x . Multiplicando la ecuación diferencial por 7(x), obtenemos
dz
6 .1 0 . R e s u e lv a xz = - x .
dx
Ésta es una ecuación diferencial para la función desconocida z(x). Tiene la forma de la ecuación (6.1) con y reem
plazado por z y p (x) = q( x) = —x. El factor de integración es
www.FreeLibros.me
P ro blem as resuelto s 45
ze~x‘^2 = e~x't2 + c
La solución para esta ecuación diferencial está dada en el problema 6.10 como
z(x )= l+ c e * ¡2
- 4 = z(0) = 1 + ce0 = 1 + c
o bien c = - 5 . D e este modo
z ( x ) = l - 5 e x‘/2
x~2z ' - 2 x ~ * z = —x 2
3
o bien — (x _ 2 z) = —x 2
dx 3
-i 2 ,
X 2Z = — X 3 + C
9
2
de donde *(*) = cx 2 + - x s
9
Ésta es una ecuación diferencial lineal para la función desconocida Q(t). Tiene la forma de la ecuación (ó. 1) con y
reemplazada por Q, x reemplazada por t, p ( t ) = 2/(10 + 2r) y q(t) = 4. El factor de integración es
www.FreeLibros.me
46 C a p ít u l o 6 E c u a c io n e s d i f e r e n c ia l e s l in e a l e s d e p r i m e r o r d e n
(1 0 + 2 » ) ^ + 2 2 = 40 + 8»
dt
. 40»+ 4 »2 + c ,
de donde C(») = — --------------- (t> -5 )
10 + 2 »
dO 2
6.14. R esu elv a e l problem a d e valor in icia l —— q--------------Q = 4; 0 ( 2 ) = 100.
dt 10 + 2»
* > - * £ 5 “ •> -»
Aplicando directamente la condición inicial, tenemos
100 = g ( 2 ) = 4 0 (2 )+ _4(4) + c
10 + 2 (2)
l ( t ) = eSk,u = e b
eh — + keh T = \0 0 k e b
dt
Teb — 100eh + c
Esta ecuación no es lineal. Sin embargo, es una ecuación diferencial de Bem oulli que tiene la forma de la ecuación
(6.4) con p(x) = q(x) - x y n = 2. Hacemos la sustitución sugerida por (6.5), específicamente, z = y 1 -2 = y-1 , de lo que
sigue
www.FreeLibros.me
P roblem as resuelto s 47
Z , x x
— r- + — = - r o bien z —xz = - x
Zl Z z2
Esta última ecuación es lineal. Su solución se encuentra en e l problema 6 .1 0 com o z = ce* '2 + 1 . Entonces, la solución
de la ecuación diferencial original es
1 1
Z ce*'12 + 1
6.17. R e s u e lv a y ' —^ y = x * y ^ 3 .
Esta última ecuación e s lineal. Su solución se encuentra en el problema 6.12 com o z = ex 2 + | x 5. D ado que z = y 2^ , la
solución del problema original está implícitamente dada por y 2/3 = ex 2 + | x 5, o explícitamente por y - ± { c x 2 + | x 5) 3^2.
6.18. D em u e str e q u e e l fa cto r d e in te g ra ció n h a lla d o e n e l p ro b lem a 6 .1 e s ta m b ién un fa cto r d e in teg ra ció n tal c o m o
s e le d e fin e e n e l c a p ítu lo 5 , e c u a c ió n (5 .7 ).
4
dx
= 3, + 6
d y = ( 3 y + 6)dx
o bien (3 y + 6 )á x + ( - l ) 4 y = 0 (7)
Estableciendo M ( x , y ) = 3ye~3x + 6 e ~ 3* y N ( x , y ) = - e ~ 3x
dM , _ 3x dN
tenem os = 3e = ——
ay ax
de lo cual concluim os que ( 2 ) es una ecuación diferencial exacta.
6.19. ■ E n cu en tre la fo rm a g e n e r a l d e la s o lu c ió n d e la e c u a c ió n ( 6 .7 ).
www.FreeLibros.me
48 C a p ít u l o 6 E c u a c io n e s d if e r e n c ia l e s l in e a l e s d e p r im e r o r d e n
se sigue que a partir de la regla de derivación del producto tenemos que el lado izquierdo de ( / ) se iguala a —
De este modo, ( / ) se puede volver a escribir así ^
í\.jg/t-MAy] = e /r(-)A ,W (2 )
f í \ M ')iz¥ =f eI PWdXcl(x)dx
y= e/o W *
PROBLEMAS ADICIONALES
En los problem as del 6.20 al 6.49, resuelva las ecuaciones diferenciales dadas.
6.20. 6.21.
tdx + 5> = ° ?dx - 5> = °
II
£
- n ,
J
dN
6.44. * - ± z=o 6.45. — = kN , (k = una constante)
* 21 dt
www.FreeLibros.me
P r o b l e m a s a d ic io n a l e s
www.FreeLibros.me
A p l ic a c io n e s
DE LAS ECUACIONES
DIFERENCIALES
7
DE PRIMER ORDEN
FR O BLEM A S D E C R E C IM IE N T O Y D EC A IM IENTO
S u p o n g a m o s q u e N (t) d en o ta la cantidad d e su stan cia (o p o b la ció n ) q u e está en crec im ien to o b ien en d ecaim ien to.
S i a su m im o s q u e d N /d t, la razón d e ca m b io en e l tiem p o de esta can tidad de su stan cia, e s p rop orcion al a la cantidad
d e su stan cia p resen te, en to n c e s d N /d t = kN , o bien
-d » - k N = 0 (7 .1 )
dt
— + k T = kTm (7 .2 )
dt m
d o n d e k e s una co n sta n te de p rop orcion alid ad p o s itiv a . U n a v e z q u e k se e s c o g e p o sitiv a se req u iere e l sig n o m en o s
e n la le y d e N e w to n para h acer q u e d T /d t se a n egativa e n un p ro ceso d e en friam ien to, cu an d o T e s m ayor q u e Tm. y
p o sitiv a en un p r o c e so d e calen tam ien to, cu an d o T e s m en or q u e Tm (v éa n se p rob lem as 7 .8 -7 .1 0 ). •
PR O BLEM A S DE C AÍDA DE C U ER PO S
50
www.FreeLibros.me
P ro blem as d e c a íd a d e c u e r p o s 51
(7.3)
Para el problem a que n os ocupa, existen d os fuerzas que actúan sobre e l cuerpo: 1) la fuerza debida a la gravedad
dada por e l p eso w del cuerpo, que se igu ala a m g, y 2 ) la fuerza debida a la resisten cia d el aire dada por - k v , donde
k > 0 e s una con stan te de proporcionalidad. S e n ecesita e l sig n o m en os porque esta fuerza se op on e a la velocidad;
e s decir, actúa en la dirección ascen dente, o negativa (véase figura 7 -1). La fuerza neta F sobre el cuerpo es, por lo
tanto, F = m g - k v . S u stituyendo este resultado en (7 .3 ), ob tenem os
, dv
m g — kv = m —
dt
o bien (7 .4 )
dv
— = g (7 .5 )
(7 .6 )
V mg
V0 gal
Suelo /g al/m in
Dirección x positiva
www.FreeLibros.me
52 C a p ít u l o 7 A p l ic a c io n e s d e l a s e c u a c io n e s d if e r e n c ia l e s o e p r im e r o r d e n
A d verten cia : Las ecu a cion es (7 .4 ), (7 .5 ) y (7 .6 ) só lo son válidas si se satisfacen las con d icion es dadas. Estas
ecu a cio n es no son válidas si, por ejem plo, la resistencia del aire no e s proporcional a la velocidad sin o al cuadrado
d e la velocidad, o si la dirección ascendente se tom a com o positiva. (V éan se los problem as 7 .1 4 y 7 .1 5 .)
PROBLEMAS DE DISOLUCIÓN
C on sid érese un tanque qué inicialm ente con tien e Vogal de salmuera que con tien e a Ib de sal. Otra solu ción de sal
m uera, que contiene b Ib de sal por galón, se vierte en el tanque a una tasa o ritmo de e gal/m in en tanto que, fim u l-
táneam ente, la solu ción bien agitada abandona el tanque a un ritmo d e /g a l/m in (figura 7 -2). El problem a e s en con
trar la cantidad de sal en e l tanque en cualquier tiem po t.
A q u í, Q denota la cantidad (en libras) de sal que se encuentra en e l tanque en cualquier tiem po t. La razón o tasa
de cam bio de Q respecto al tiem po, d Q /d t, se iguala al ritmo al cual la sal ingresa al tanque m en os el ritmo al cual la
sal deja e l tanque. La sal entra al tanque a un ritmo de b e lb/m in. Para determ inar e l ritm o al cual la sal abandona el
tanque, primero ca lcu lam os e l volu m en de salm uera en e l tanque en un tiem po t determ inado, que es el volum en
inicial V0 m ás el volum en de salm uera e l agregado m enos e l volum en de salm uera f t extraído. A sí, e l volum en de
salm uera en cualquier tiem po es
V0 + e t - f t (7.7)
L a concentración de sal en el tanque en un m om ento dado es Q /(V0 + e t - f t) , de lo que se desprende que la sal sale
d el tanque a una tasa de
/ lb/m in
V0 + et - f t
dQ
D e este m odo, = b e —f
dt V0 + e t - f t
dQ f
o bien Q = be (7.8)
dt ' V0 + ( e - f ) t
(V éan se problem as 7 .1 6 -7 .1 8 .)
CIRCUITOS ELÉCTRICOS
L a ecuación básica que gobierna la cantidad de corriente I (en am perios) en un circuito RL sim p le (figura 7-3 ) co n
sistente en una resistencia R (en oh m ios), un inductor L (en henrios) y una fuerza electrom otriz (abreviado fem ) E
(en v oltios) es
7d t +7/
L =7
L (79)
Para un circuito RC con sistente en una resistencia, una capacitancia C (en faradios), una fem , y sin inductancia (figu
ra 7 -4 ), la ecuación que gobierna la cantidad de carga eléctrica q (en cu lom b ios) sobre el capacitor es
^■ + — q = ~ (7 .1 0 )
dt RC R
L a relación entre q e l e s
/ = — (7-11)
dt
www.FreeLibros.me
P roblem as resuelto s 53
I
F igu ra 7-4
T R A Y E C T O R IA S O R T O G O N A L E S
F (x ,y ,c ) = 0 ( 7 .7 2 )
d o n d e c in d ic a e l p arám etro. E l p ro b lem a c o n sis te en en con trar otra fa m ilia d e cu rvas d e un p arám etro, lla m a d a s las
tr a y e c to r ia s o r to g o n a le s d e la fa m ilia ( 7 .7 2 ) y d ad as a n a lítica m en te por
G (x , y , k ) = 0 (7 .7 5 )
= /( * .> ) { 7 .1 4 )
dx
dy _ 1
( 7 .7 5 )
d x ~ ~ f ( x ,y )
( V é a n s e p ro b lem a s 7 .2 3 - 7 .2 5 .)
Para m u c h a s fa m ilia s d e cu rv a s, n o s e p u ed e re so lv e r e x p líc ita m e n te para d y ld x y ob ten er u n a e c u a c ió n d ife ren cia l
d e la fo rm a (7 .7 4 ). N o s e co n sid e r a n ta le s cu rv a s en e s t e lib ro.
P R O B L E M A S R E SU E L T O S
Aquí, N (t) indica e l balance en la cuenta en cualquier tiem po t. Inicialm ernte, N (0) = 20 000. El balance de la
. cuenta crece por m edio de los pagos de intereses, que son proporcionales a la cantidad de dinero en la cuenta. La con s
tante de proporcionalidad es la tasa de interés. En este caso, k - 0.05 y la ecuación (7.7) se convierte en
— -0 .0 5 N = 0
dt
Esta ecuación diferencial es tanto lineal com o separable. Su solución es
N ( t) = ce 0 05' (7)
www.FreeLibros.me
54 C a p ít u l o 7 A p l ic a c i o n e s d e l a s e c u a c i o n e s d if e r e n c ia l e s d e p r i m e r o r d e n
20 0 00 = ce 0 05(0) = c
b) Buscamos el tiempo t en el que e l balance sea N (t) = $40 000. Sustituyendo estos valores en (2) y resolviendo para
t, obtenemos
40 000 = 20 000e° 051
2 = e 0 05'
ln|2| = 0.05r
7 .2 . U n a p ersona d ep o sita $5 0 0 0 e n una cu en ta que acum ula interés co m p u e sto de m anera con tinu a. A su m ien d o
q u e n o hay ex tra ccio n es n i d e p ó sito s a d icio n a les, ¿cuánto habrá en la cu en ta d esp u és de sie te añ os si la tasa
d e in terés es d e l 8 .5 por cien to co n sta n te durante lo s prim eros cuatro añ os y d e l 9 .2 5 p or cie n to con stan te lo s
tres años sig u ien tes?
Aquí, N(t) denota el balance de la cuenta en un tiempo t. Inicialmente, N(0) = 5 000. Para los primeros cuatro años,
k = 0.085 y la ecuación (7.7) se convierte en
— - 0 .0 8 5 V = 0
dt
Su solución es
N (t) = ceo m ' (0 < r < 4 ) (/)
5 0 0 0 = ce 0 085(0) = c
y (7) se convierte en
N ( r ) = 5 0 0 0 e ° 085' (0 < f < 4 ) ( 2)
Esta cantidad también representa el balance para el com ienzo del periodo de los últimos tres años.
En los tres últimos años, la tasa de interés es de 9.25 por ciento y (7.7) se convierte en
d/V
- 0 .0 9 2 5 IV = 0 (4 < r < 7 )
dt
Su solución es
7V(r) = cea m 5 t (4 < í< 7 ) (3)
y (3) se convierte en
www.FreeLibros.me
P ro blem as resuelto s 55
7.3. ¿ Q u é tasa d e interés con stan te se requiere si e l d ep ó sito in icia l c o lo c a d o en una cu en ta q u e acum ula interés
co m p u esto co n tin u a m en te d eb e d u p licar su valor en se is añ os?
dN
- — kN = 0
dt
que tiene com o solución
N (t) = ce1“ (I)
N o nos dan una cantidad por el depósito inicial, de modo que lo indicamos com o iV0. En t = 0, N(0) = N0, que cuando se
sustituye en (1) da
N0 = ce*(0> = c
y (1) se convierte en
Buscamos el valor de k para el cual N = 2N0 cuando r = 6 . Sustituyendo estos valores en (2) y resolviendo para k, encon
tramos que
2N 0 = N 0ekm
e6k = 2
6* = ln | 2 |
* = i l n |2 | = 0 .1 l5 5
6
7.4. S e sa b e qu e un cu ltiv o d e b acterias crec e a un ritm o p rop orcion al a la can tidad p resen te. D esp u é s de una hora,
se ob servan en e l c u ltiv o 1 0 0 0 co lo n ia s; y lu e g o d e tres horas, 3 0 0 0 c o lo n ia s. E ncuentre a.)una ex p resió n
para e l n ú m ero ap ro xim ad o d e c o lo n ia s d e b acterias p resen tes en e l c u ltiv o en un tiem p o t y b) e l núm ero
ap roxim ado de c o lo n ia s d e bacterias q u e había origin alm en te en e l cu ltiv o .
a) Aquí, N(t) indica el número de colonias de bacterias en el cultivo en un tiempo t. D e (6.1), dN /dt - kN = 0, que es
tanto lineal com o separable, su solución es
N (t) = ceh W
En f = 1, N = 1 000; de aquí
1000 = ce* (2 )
En / = 4, N = 3 000; de aquí
3000 = ce4* (3 )
www.FreeLibros.me
56 C a p ít u l o 7 A p l ic a c io n e s d e l a s e c u a c io n e s d i f e r e n c ia l e s d e p r i m e r o r d e n
7 .5 . S e sabe que la p o b lación de cierto p aís se increm enta a un ritm o proporcional al núm ero de person as que viven
en él. S i d esp u és de d o s años la p o b la ció n se ha d uplicado, y si lu eg o d e tres añ os la p ob lación e s de 2 0 0 0 0 ,
estim e e l núm ero de personas que v ivían in icialm en te en el país.
Si N indica el número de personas que viven en el país en un tiempo f, y N0 denota el número de personas que vivían
inicialmente en él, entonces, de (7.1),
— -iW = 0
dt
cuya solución es
N = ceh (1)
N = N 0 e fa (2)
7 .6 . S e sabe que cierto m aterial radiactivo d ecae a un ritm o proporcional a la cantidad presente. S i in icialm en te
hay p resen tes 5 0 m iligram os d e m aterial y d esp u és de d os horas se ob serva q u e e l m aterial ha perdido 10 por
cien to de su m asa origin al, en cuentre á ) una exp resión para la m asa d e m aterial rem anente en un tiem p o t,
b ) la m asa d el m aterial lu e g o de cuatro horas y c ) el tiem p o en que e l m aterial ha d eca íd o h asta la m itad d e
su m asa in icial.
— ~kN = 0
dt
Ésta es una ecuación diferencial separable y lineal; su solución es
N = ceb ( /)
En t = 0, se dijo que N = 50. Por lo tanto, a partir de (1), 50 = ce*(0), o bien c = 50. D e este modo,
N = 50eb ( 2)
En t = 2 ha decaído 10 por ciento de la masa original de 50 mg, o sea 5 mg. Por esto, en r = 2, N = 50 - 5 = 45.
Sustituyendo estos valores en (2), y resolviendo para k, encontramos que
Sustituyendo este valor en (2), obtenemos la cantidad de masa presente en un tiempo cualquiera t como
N = 50e"° °53' W
donde t se mide en horas.
www.FreeLibros.me
P ro blem as resuelto s 57
c) Necesitamos t cuando N = 50/2 = 25. Sustituyendo N = 25 en (5) y resolviendo para t, encontramos que
Sea N(t) lo que indique el número de ratones con la enfermedad en el tiempo t. Se dijo que N(0) = 5, de lo que se
desprende que 500 - N(t) es el número de ratones sin la enfermedad en el tiempo t. La teoría predice que
dN
— = kN (5 0 0 - N ) ( 1)
dt
donde k es una constante de proporcionalidad. Esta ecuación es diferente a (7.1) porque la tasa de cambio ya no es más
proporcional al número justo de ratones que tiene la enfermedad. La ecuación (1) tiene la forma diferencial
dN
- —k d t = 0 ( 2)
N (5 0 0 -N )
que es separable. Utilizando la descom posición en fracciones parciales, tenemos
1 1 /5 0 0 1 /5 0 0
-+ -
Af(500 —W) N 500 - N
1 í 1 • 1
d N -k d t = 0
50oU + 500 —IVj
Su solución es
o bien ^ ( l n |A / |- l n |5 0 0 - / V |) - f r = c
N . _ e S00(c4-*<) (2)
500 —N
— 2 — = c1e 500b
500- N 1
- = c ,e 500‘« » = <
495
de m odo que C\ = 1/99 y (4) se convierte en
N 1
(5)
500 —IV 99
www.FreeLibros.me
58 C a p ít u l o 7 A p l ic a c io n e s d e l a s e c u a c io n e s d if e r e n c ia l e s d e p r i m e r o r d e n
Podríamos resolver (5) para N, pero esto no es necesario. Buscamos un valor de t donde N = 250, la mitad de la población.
Sustituyendo N = 250 en (5) y resolviendo para t, obtenemos
l = — eS00h
. 99
99 = e 500b
ln 99 = 500it
o bien / = 0.00919//: unidades de tiempo. Sin tener información adicional, no podemos obtener un valor numérico para la
constante de proporcionalidad k o bien ser más definitivos con respecto a t.
7 .8 . S e c o lo c a una barra de m etal a 100° F en un cuarto a temperatura con stan te de 0° F. S i d esp u és d e 2 0 m inu tos,
la temperatura de la barra e s de 5 0 ° F, en cuentre a ) e l tiem p o que tom ará para q u e la barra alcan ce la tem p e
ratura de 2 5 ° F y tí) la tem peratura d e la barra lu eg o de 10 m inu tos.
Use la ecuación (7.2) con Tm = 0; el medio aquí es el cuarto que se está m antenendo a una temperatura de 0° F. De
este modo, tenemos
— + kT = 0
dt
cuya solución es
T = ce ( /)
Puesto que T = 100 en t = 0 (la temperatura de la barra es inicialmente 100° F), de ( i) tenemos que 100 = ce- *(0) o bien
100 = c. Sustituyendo este valor en ( /), obtenemos
7 = 100e-b (2)
En t = 20, se dijo que T = 50; por lo tanto, de (2),
7 = lOOe'0 035' ( 3)
b) Requerimos T cuando t = 10. Sustituyendo f = 10 en (3) y luego resolviendo para T, encontramos que
Debería notarse que dado que la ley de Newton sólo es válida para pequeñas diferencias de temperatura, los
cálculos anteriores representan solamente una primera aproximación a la situación física.
Usando (7.2) con Tm = 100 (el medio circundante es el aire exterior), tenemos
— + kT = 100k
dt
www.FreeLibros.me
P r o b lem a s resuelto s 59
Esta ecuación diferencial es lineal. Su solución está dada en el problema 6.15 com o
r = c e - * '+ 1 0 0 (i)
Dado que T = 50 cuando t = 0, se sigue de (1) que 50 = ce~ t(0> + 100, o bien c = - 5 0 . Sustituyendo este valor en (7),
obtenemos
En t = 5, sabemos que 7 = 60; por eso, de (2), 60 = -5 Q e~ sk + 100. R esolviendo para k, obtenemos
—1 40 —1
-4 0 = -5 0 e o bien k = — ln — = - ( - 0 . 2 2 3 ) = 0.045
5 50 5
Sustituyendo este valor en (2), obtenemos la temperatura del cuerpo en cualquier tiempo t com o
- 0 .0 4 5 f = l n i o bien f = 15.4m in
2
De (7.2),
dT
— + * 7 = 30*
dt
Resolviendo, obtenemos
T = ce~ b + 30 ( 1)
www.FreeLibros.me
60 C a p ít u l o 7 A p l ic a c io n e s d e l a s e c u a c io n e s d i f e r e n c ia l e s d e p r im e r o r d e n
Sustituyendo estos valores en (1), tenemos que la temperatura del cuerpo en cualquier momento t es
T = —6 0 e -0069' + 30 (4)
r 0 = - 6 0 e <'° ° 69KO) + 30 = - 6 0 + 30 = - 3 0 ° F
7 .1 1 . U n cuerpo que tien e c in c o u nidades técn icas de m asa se deja caer d esd e una altura de 100 p ie s c o n velocid ad
cero. A su m ien d o que n o hay resisten cia d e l aire, encuentre a ) una exp resión para la v elo cid a d d el cu erp o en
u n determ inado tiem p o t, b ) una exp resión para la p o sició n d el cuerpo en cualquier tiem p o t y c ) e l tiem p o
requerido para llegar al su elo.
a) Escoja el sistema de coordenadas com o en la figura 7-5. Luego, dado que no hay resistencia del aire, se aplica (7.5):
dvldt = g. Esta ecuación diferencial es lineal o, en forma diferencial, separable; su solución es v = gt + c. Cuando
I = 0, v = 0 (inicialmente el cuerpo tiene velocidad 0); de aquí que 0 = g(0) + c, o bien c = 0. D e este modo, v = gt o
bien, asumiendo g = 32 pies/seg2,
v = 32í U)
b) Recuerde que la velocidad es la razón de cambio en el tiempo del desplazamiento, designado aquí por x, Por eso,
v = dx/dt, y (1) se convierte en dx/dt = 32r. Esta ecuación diferencial es tanto lineal como separable; su solución es
z = 16t2 + c 1 (2 )
Pero en t = 0, x = 0 (véase figura 7-5). D e este modo, 0 = (16)(0)2 + q , o bien q = 0. Sustituyendo este valor en
( 2 ), tenemos
jr = 16r2 (5)
7 .1 2 . U n a b o la de acero que p esa 2 Ib se d eja caer d esd e una altura d e 3 0 0 0 p ies sin v e lo c id a d . M ien tras cae,
encuentra una resisten cia d el aire num éricam en te igu al a v/ 8 (en libras), d on d e v in d ica la velo cid a d d e la b ola
(en p ies por se g u n d o ). E ncu en tre a ) la v e lo c id a d lím ite para la b o la y b ) el tiem p o req u erid o para que la
b o la im pacte en e l su elo.
Ubique el sistema de coordenadas com o en la figura 7-5 con el suelo situado ahora en x = 3 000. Aquí, w = 2 Ib y
k = 1/8. Asumiendo que la gravedad g es 32 p ie s/seg12, tenemos de la fórmula w = mg que 2 = m(32) o que la masa de
la bola es m = 1/16 unidades técnicas de masa. La ecuación (7.4) se convierte en
www.FreeLibros.me
P ro blem as resuelto s 61
0 = ce~2m + 16e = c + 16
v(í) = —16e~2' + 16 (2 )
a) D e (7) o bien (2), vem os que conforme r —>oo, v —»16 de modo que la velocidad límite es 16 pies/seg2.
b) Para hallar el tiempo quetarda la bola en impactar en el suelo (x = 3 000), necesitamos una expresión para la posición
de la bola en cualquier tiempo t. Dado que v = dx/dt, (2) se puede volver a escribir como
— = - 1 6 e - 2' + 1 6
dt
Integrando ambos lados de esta última ecuación directamente con respecto a /, tenemos
donde ci indica una constante de integración. En t = 0, x = 0. Sustituyendo estos valores en (3), obtenemos
0 = 8 e _2(0) + 1 6 (0 ) + c, = 8 + o,
La bola golpea el suelo cuando x(t) = 3 000. Sustituyendo este valor en (4), tenemos
Aunque (5) no se puede resolver explícitamente para r, podemos aproximar la solución por medio de la prueba de
ensayo y error, sustituyendo diferentes valores de t en (5) hasta que encontremos una solución del grado de exactitud
que necesitamos. D e manera alternativa, vem os que para cualquier valor grande de r, el término exponencial negati
vo puede ser despreciable. Una buena aproximación consiste en colocar 2 1 = 376 o bien t = 188 seg. Para este valor
de í, el exponente es esencialmente 0 .
U n cu erp o q u e p esa 6 4 Ib se d eja caer d esd e una altura d e 1 0 0 p ies c o n u na velo cid a d in icia l d e 10 p ies/seg .
A su m a que la resisten cia d el aire e s p rop orcion al a la v elo cid a d d el cuerpo. S i se sabe que la v elo cid a d lím i
te e s 128 p ie s/s e g , en cu en tre a ) una ex p resió n para la v elo cid a d d e l cu erp o e n cu alq u ier tiem p o t, b ) una
ex p resió n para la p o sic ió n d el cuerpo en cualq u ier tiem p o t.
a) Localice el sistema de coordenadas com o en la figura 7-5. Aquí, w = 64 Ib. Dado que w = mg, tenemos que mg = 64,
o bien m = 2 unidades técnicas de masa. Dado que v, = 128 pies/seg, de (7.6) tenemos que 128 = 64/ k , o bien k =
Sustituyendo estos valores en (6.4), obtenemos la ecuación diferencial
dv 1
— + - v = 32
. dt 4
que tiene la solución
En / = 0, sabemos que v = 10. Sustituyendo estos valores en (7), tenemos 10 = ce0 + 128, o bien c = - 1 1 8 . La
velocidad en cualquier tiempo t está dada por
v = -1 1 8 e-'''4 + 1 2 8 (2 )
www.FreeLibros.me
62 C a p ít u l o 7 A p l ic a c io n e s d e l a s e c u a c io n e s d if e r e n c ia l e s d e p r im e r o r d e n
— = - 1 1 8 « - ,''4 + 128
di
Esta última ecuación, en forma diferencial, es separable; su solución es
7 .1 4 . U n cuerpo de m asa m se arroja verticalm ente al aire co n una velocid ad inicial v0. S i el cuerpo experim enta
una resistencia d el aire proporcional a su velocid ad , encuentre a ) la ecu ación de m ovim ien to en e l sistem a de
coorde.tadas d e la figura 7 -6 , b ) una exp resión para la velocid ad d el cuerpo en cualquier tiem p o I y c) el
tiem p o en e l cual el cuerpo alcanza su altura máxima.
| Dirección x positiva
a) En este sistema de coordenadas, la ecuación (7.4) tal vez no sea la ecuación de movimiento. Para deducir la ecuación
adecuada, vemos que hay dos fuerzas sobre el cuerpo; 1 ) la fuerza debida a la gravedad, dada por mg y 2 ) la fuerza
debida a la resistencia del aire, dada por kv, que impide la velocidad del cuerpo. Dado que estas dos fuerzas actúan
en la dirección descendente o negativa, la fuerza neta sobre el cuerpo es -m g -kv. Utilizando (7.3) y reagrupando,
obtenemos
dv k ,
— + —v = -g U)
di m
como la ecuación de movimiento.
v = | v 0 + M | e-<V",) » _ M (2)
www.FreeLibros.me
P ro blem as r esu elto s 63
c) El cuerpo alcanza su máxima altura cuando v = 0. Por esto, necesitam os t cuando v = 0. Sustituyendo v = 0 en (2) y
resolviendo para t, hallamos
-(*/«)■ 1 .
1+ -
mg
1
—(* / m )t = In
1+ M
mg
< k l mg
La fuerza debida a la resistencia del aire es -k v 1, de modo que la segunda ley del movimiento de Newton se con
vierte en
dv , 2 . . „dv , ,
m — = mg — kv o b len 2— = 6 4 -á v
dt dt
2
dv —d t = 0 (/)
64 —kv2
que es separable. Por m edio de fracciones parciales,
1
_ _ _ _2_ _ _ ■ 2 s o| o
64 — ( 8 - V * v X 8 + V jtv)_ 8 -V ifcv 8 + V*v
---------------------^
—¡ = - \ d v - d t — Q
8 (,8 -\ík v S + yfkv)
7 f í + í- p - j d v - f d t= 'c
\% -s lk v 8 + vifevJ J
o bien -^ H 8 —^ v| + - ^ ! n |8 + yfkv\ —t — c
8 + V¡tv
ln = &\[kt + iy fk c
8 - sfkv
8 + yfkv
o bien = c ,e 8Æ (c, = ± e s Æ )
8 - Æ
8 + Vfcv g /j, 8 , . rr
= o b le n v = - 5= ta n h 4 V/cí
8 -V á v Jk
Obsérvese que sin información adicional, no podemos obtener un valor numérico para la constante k.
www.FreeLibros.me
64 C a p ít u l o 7 A p l ic a c io n e s d e l a s e c u a c io n e s d if e r e n c ia l e s d e p r im e r o r d e n
7 .1 6 . U n tanque con tien e inicialm ente 100 gal d e una solu ción de salm uera c o n 20 Ib de sal. En t = 0 se vierte agua
fresca al tanque a un ritm o de 5 gal/m in, m ientras que la solu ción bien agitada deja el tanque al m ism o ritmo.
Encuentre la cantidad de sal en e l tanque en cualquier tiem po t.
^ - + — Q= 0
dt 20
La solución de esta ecuación lineal es
Q = ce 0 /)
En í = 0, sabemos que Q = a = 20. Sustituyendo estos valores en (1) encontramos que c = 20, de modo que ( /) se puede
volver a escribir como Q = 20e~"i0. Obsérvese que conforme / —» °° , Q —> 0 tal como debiera de ser, puesto que sólo se
está agregando agua fresca.
7 .1 7 . U n tanque con tien e in icialm en te 100 gal d e una solu ción de salm uera c o n 1 Ib de sal. E n t = 0 se vierte otra
so lu ció n de salm uera que con tien e 1 Ib de sal por galón a un ritmo de 3 gal/m in , en tanto que la salm uera bien
agitada abandona el tanque al m ism o ritmo. Encuentre a ) la cantidad de sal en el tanque en cualquier tiem po r,
y b ) el tiem po en el cual la m ezcla en el tanque con tien e 2 Ib de sal.
— + 0.03@ = 3
dt
La solución a esta ecuación diferencial es
Q=
= ^
ce —003. + 100 ( /)
En t = 0, Q = a = 1. Sustituyendo estos valores en ( /), encontramos que 1 = ce0 + 1 0 0 , o bien c = —99. Entonces
( /) se puede volver a escribir como
7 .1 8 . U n tanque con tien e in icialm en te 10 gal d e agua fresca. En t = 0 se vierte al tanque una so lu ció n de salm uera
que co n tien e 1 Ib de sal por galón a un ritm o de 4 gal/m in , m ientras q u e la m ezcla bien agitada abandona el
tanque a un ritm o de 2 gal/m in. Encuentre a ) la cantidad de tiem p o requerido para que ocurra e l derrame o
desborde, y b) la cantidad de sal en e l tanque al m om en to d el derrame.
Ésta es una ecuación lineal; su solución está dada en el problema 6.13 como
40r + 4r2 + c
e= - U)
10 + 2/
www.FreeLibros.me
P ro b lem a s resueltos 65
En t = 0, Q = a = 0. Sustituyendo estos valores en (7), encontramos que c = 0. Necesitamos Q al momento del derra
me, que de la parte (a) es t = 20. De este modo,
40(201 + 4(20)* = 4 g lb
10 + 2(20)
7 .1 9 . U n circuito R L tiene una fe m (fuerza electrom otriz) d e 5 v o ltio s, una resisten cia d e 5 0 o h m io s, una inductan-
c ia de 1 henrio, y ninguna corriente in icial. Encuentre la corriente en e l circuito en cualquier tiem p o t.
— + 5 0 /= 5
dt
La ecuación es lineal; su solución es
/ = - - L e - * > '+ - L V)
10 10
1 .a cantidad - '¡oe en ( /) se llama corriente transitoria, pues esta cantidad va hasta cero ("se extingue”) con
forme t —» °°. La cantidad -¡q en ( /) se llama corriente de estado estacionario. Conforme t —» °°, la corriente / se aproxi
ma al valor de la corriente de estado estacionario.
7 .2 0 . U n circuito R L tien e una fe m (en v o ltio s) dada por 3 sen 2r, una resisten cia de 10 o h m io s, una inductancia
de 0 .5 h en rio s, y una corrien te in icia l d e 6 am perios. E ncu en tre la corriente en e l circu ito para cualquier
tiem p o t.
— + 20/ = 6 sen 2 r
dt
Esta ecuación es lineal, con solución (véase capítulo 6 )
Llevando a cabo las operaciones de integración (la segunda integral requiere dos veces la aplicación de operaciones de
integración por partes), obtenemos
20 2 3
6 = ce_JO<0) H sen 2 ( 0 ) eos 2 (0 ) o bien 6 = c --------
101 101 101
de donde c = 609/101. La corriente en cualquier tiempo t es
, 609 _20, 30 3
/ = e A sen 2 f co s 2 r
101 101 101
Como en el problema 7.18, la corriente es la suma de una corriente transitoria, dada por (609 / 101)e-20', y una corrien
te de estado estacionario,
3 0 - 3
— sen 2 r -------- eos 2 r
101 101
www.FreeLibros.me
66 C a p it u l o 7 A p l ic a c io n e s d e l a s e c u a c io n e s d if e r e n c ia l e s d e p r im e r o r d e n
7 .2 1 . R eescriba la corriente de estado estacion ario d el problem a 7 .2 0 en la form a A sen (2t — <¡>). E l ángulo <p es
llam ado á n g u lo d e fa s e .
30 3
/ , = — sen 2 1 ------- c o s 2 í = A cosó sen 2/ ~ Asen b eos 2;
' 101 101
30 3
De este modo. A eos ó = — y A sen t¡>= — , de lo que se desprende que
101 101
A señó
y tan^ = -------- - J -L 1 /Í Ü L L J -
A cos 0 U o ij/ U oiJ io
En consecuencia, 1, tiene la forma requerida si
909 3 1
l= J y = ¡ y O = arelan— = 0.0997 radianes
\(1 0 1 )2 VÍ01 r 10
7 .2 2 . U n circuito RC tiene una fem (en v oltios) dada por 4 0 0 eo s 2r, una resistencia de 100 oh m ios y una cap aci
tancia de 10 " 2 faradios. In icialm en te no hay carga en el capacitor. Encuentre la corriente en el circuito en
cualquier tiem po t.
Primero hallamos la carga q y luego usamos (.7.11) para obtener la corriente. Aquí, E = 400 cos'2f, R = 100 y
C = 10"2; por lo tanto, (7.10) se convierte en
da
— + q = 4 cos2f
dt
La ecuación es lineal, y su solución es (se requieren dos operaciones de integración por partes)
q —ce - 1 + j s e n 2 í + jC O s 2 r
En r = 0, q = 0; de aquí,
4 8 4
de este modo q = — e " ' + - s e n 2 fH— co s 2 r
* 5 5 5
, dq 4 _ 16 „ 8 _
/ = - = - ? ' + — co s 2 r — sen 2 r
dt 5 5 5
La familia, que está dada por (7.12) con F(x, y, c) = jc2 + y2 - c2, consiste de círculos con centros en el origen y
radios c. Derivando implícitamente la ecuación dada con respecto a x. obtenemos
2x + 2yy'= 0 ob ien — = - —
dx y
www.FreeLibros.me
P ro blem as resuelto s 67
£( y _ y
dx x
y = kx ( 1)
La familia, que está dada por (7.12) con F(x, y , c ) = y - ex2, consiste de parábolas simétricas alrededor del eje y con
vértices en el origen. Derivando la ecuación dada con respecto a x, obtenemos dxJdy = 2ex. Para eliminar c, observamos,
de la ecuación dada, que c = y / x 2-, por lo tanto, dy/dx = 2y/x. Aquí f ( x , y ) = 2y Ix , así (7.15) se convierte en
La solución de esta ecuación separable es %x2 + y 2 = k. Estas trayectorias ortogonales son elipses. Algunos miembros
de esta familia, junto con ciertos miembros de la familia original de parábolas, se muestran en la figura 7-8. Obsérvese
■ que cada elipse intersecta a cada parábola en ángulos rectos.
Aquí, F (x, y , c ) = x 2 + y 2 - ex. Derivando implícitamente la ecuación dada con respecto a x, obtenemos,
2 x+ 2 y— = c
dx
www.FreeLibros.me