Analisis Iv
Analisis Iv
Analisis Iv
ALUNMOS (GRUPO 8)
AYACUCHO – PERU
2020
EJERCICIOS 1.3
18)Un barril de s pies de diámetro con peso de w lb está flotando en el
agua, como lo muestra la figura 134a). Después de un hundimiento
inicial, el barril exhibe un movimiento oscilatorio hacia arriba y hacia
abajo a lo largo de una línea vertical. Mediante la figura 1.34b).
determine una ecuación diferencial para el desplazamiento vertical y(t)
si se considera que el origen está sobre el eje vertical en la superficie
del agua cuando el barril está en reposo.
Datos
Diámetro= s
Radio= s/2
Peso= w
Densidad del agua = 62.4
Solución
Dado que el barril en la figura está sumergido “y” pies adicionales por
debajo de su posición de equilibrio, entonces el número de pies cúbicos
en la porción sumergida adicional es el volumen del cilindro circular.
𝑉 = 𝜋 ∗ 𝑅2 ∗ 𝑦
𝑠 2
𝑉 = 𝜋∗( ) ∗𝑦
2
Por el principio de Arquímedes
𝐹𝑢𝑒𝑟𝑧𝑎 𝑎𝑠𝑐𝑒𝑛𝑑𝑒𝑛𝑡𝑒 𝑑𝑒𝑙 𝑎𝑔𝑢𝑎 = 𝑃𝑒𝑠𝑜 𝑑𝑒𝑙 𝑎𝑔𝑢𝑎 𝑑𝑒𝑠𝑝𝑙𝑎𝑧𝑎𝑑𝑎 𝑠𝑜𝑏𝑟𝑒 𝑒𝑙 𝑏𝑎𝑟𝑟𝑖𝑙
= (62.4) ∗ 𝑉
𝑠 2
= 62.4 ∗ (2) ∗ 𝑦
= 15.6𝜋𝑠 2 𝑦
De la segunda ley de Newton se deduce que:
𝑤 𝑑2𝑦
= −15.6𝜋𝑠 2 𝑦
𝑔 𝑑𝑡 2
𝑑 2 𝑦 15.6𝜋𝑠 2 𝑔
+ 𝑦=0
𝑑𝑡 2 𝑤
Solución
𝑑2𝑥
𝐹 = 𝑚𝑎 = 𝑚 2 = −𝑘 (𝑠 + 𝑥 ) + 𝑚. 𝑔
𝑑𝑡
𝑑2𝑥
𝑚 = −𝑘𝑥 + 𝑚𝑔 − 𝑘𝑠
𝑑𝑡 2
𝑚𝑔 = 𝑘𝑠
𝑑2𝑥
𝑚 = −𝑘𝑥
𝑑𝑡 2
Gráfica en MATHEMATICA del ejercicio
EJERCICIO 2.3
45. Suponga que P(x) es continua en algún intervalo I y que a es
número incluido en I. ¿Qué se puede decir acerca de la solución del
problema de valor inicial 𝒚′ + 𝑷(𝒙) 𝒚 = 𝟎, 𝒚(𝒂) = 𝟎?
Solución
𝑦 ′ + 𝑃 (𝑥 )𝑦 = 0
𝑑𝑦
+ 𝑃 (𝑥 )𝑦 = 0
𝑑𝑥
Tenemos una ecuación de forma lineal
∴𝑦=0
46. Series de decaimiento radiactivo. El siguiente sistema de
ecuaciones diferenciales de encuentra en el estudio del decaimiento
de un tipo especial de series de elementos radiactivos.
𝒅𝒙
= −𝝀𝟏 𝒙,
𝒅𝒕
𝒅𝒚
= 𝝀𝟏 𝒙 − 𝝀𝟐 𝒚
𝒅𝒕
Donde 𝝀𝟏 y 𝝀𝟐 son constante. Analice como resolver este sistema
sujeto a 𝒙(𝟎) = 𝒙𝟎 , 𝒚(𝟎) = 𝒚𝟎. Ponga en práctica sus ideas.
Solución
Resolvemos la primera ecuación:
𝒅𝒙
= −𝝀𝟏 𝒙
𝒅𝒕
𝒅𝒙
∫ = − ∫ 𝝀𝟏 𝒅𝒕
𝒙
𝐥𝐧 𝒙 = −𝝀𝟏 𝒕
𝒙 = 𝒄𝟏 𝒆−𝝀𝟏 𝒕
De 𝒙(𝟎) = 𝒙𝟎, 𝒄𝟏 = 𝒙𝟎
Resolvemos la segunda ecuación:
𝒅𝒚
= 𝝀𝟏 𝒙 − 𝝀𝟐 𝒚
𝒅𝒕
𝒅𝒚
= 𝒙𝟎 𝝀𝟏 𝒆−𝝀𝟏 𝒕 − 𝝀𝟐 𝒚
𝒅𝒕
𝒅𝒚
+ 𝝀𝟐 𝒚 = 𝒙𝟎 𝝀𝟏 𝒆−𝝀𝟏 𝒕
𝒅𝒕
Encontramos que el factor de integración es: 𝒆𝝀𝟐 𝒕
𝒅 𝝀𝒕
[𝒆 𝟐 𝒚] = 𝒙𝟎 𝝀𝟏 𝒆−𝝀𝟏 𝒕 𝒆𝝀𝟐 𝒕
𝒅𝒕
𝒅
[𝒆𝝀𝟐 𝒕 𝒚]=𝒙𝟎 𝝀𝟏 𝒆(𝝀𝟐 −𝝀𝟏)𝒕
𝒅𝒕
𝒙 𝝀
𝒆𝝀𝟐 𝒕 𝒚 = 𝝀 𝟎−𝝀𝟏 𝒆(𝝀𝟐 −𝝀𝟏 )𝒕 + 𝒄𝟐
𝟐 𝟏
𝒙𝟎 𝝀𝟏
𝑦= 𝒆−𝝀𝟏 𝒕 + 𝒄𝟐 𝒆−𝝀𝟐 𝒕
𝝀𝟐 − 𝝀𝟏
De 𝒚(𝟎) = 𝒚𝟎
𝒄𝟐 = (𝒚𝟎 𝝀𝟐 − 𝒚𝟎 𝝀𝟏 − 𝒙𝟎 𝝀𝟏)/(𝝀𝟐 − 𝝀𝟏 )
La solución es:
𝒙 𝝀 𝒚𝟎 𝝀𝟐−𝒚𝟎𝝀𝟏 −𝒙𝟎 𝝀𝟏
𝑦 = 𝝀 𝟎−𝝀𝟏 𝒆−𝝀𝟏 𝒕 +
𝟐 𝟏 𝝀𝟐 −𝝀𝟏
EJERCICIO 2.7
37. Gota de lluvia evaporándose A medida que una gota de lluvia cae se
evapora, pero mientras eso sucede conserva su forma esférica. Si
suponemos adicionalmente que la velocidad de evaporación de la gota
de lluvia es proporcional a su área superficial y la resistencia del aire es
insignificante, entonces un modelo de la velocidad v(t) de la gota será.
𝑑𝑣 3(𝑘 ⁄𝜌)
+ =𝑔
𝑑𝑡 (𝑘 ⁄𝜌)𝑡 + 𝑟𝑜
Aquí 𝜌 es la densidad del agua, 𝑟𝑜 es el radio de la gota de lluvia cuando
t = 0, k < 0 es la constante de proporcionalidad y la dirección hacia
abajo se toma como positiva.
𝑘
𝑑𝑣 3( ) 𝑑𝑣 3𝑎
𝜌
+ =𝑔 → + 𝑣=𝑔
𝑑𝑡 (𝑘 ) 𝑡 + 𝑟 𝑑𝑡 𝑎𝑡 + 𝑟𝑜
𝜌 𝑜
3𝑎
𝑃 (𝑡 ) =
𝑎𝑡 + 𝑟𝑜
3𝑎
∫𝑎𝑡+𝑟 𝑑𝑡
𝑒 𝑜 = 𝑒 3ln(𝑎𝑡+𝑟𝑜 ) = (𝑎𝑡 + 𝑟𝑜 )3
3𝑎
3∫ 𝑑𝑡 → 𝑢 = 𝑎𝑡 + 𝑟𝑜 , 𝑑𝑢 = 𝑎𝑑𝑡
𝑎𝑡 + 𝑟𝑜
V(t)
𝑑𝑣
(𝑎𝑡 + 𝑟𝑜 )3 + 3𝑎(𝑎𝑡 + 𝑟𝑜 )2 𝑣 = (𝑎𝑡 + 𝑟𝑜 )3 𝑔
𝑑𝑡
𝑑
∫ [(𝑎𝑡 + 𝑟𝑜 )3 𝑣 ]𝑑𝑡 = ∫(𝑎𝑡 + 𝑟𝑜 )3 𝑔𝑑𝑡
𝑑𝑡
t(s)
𝑑𝑢
∫(𝑎𝑡 + 𝑟𝑜 )3 𝑔𝑑𝑡 → 𝑢 = 𝑎𝑡 + 𝑟𝑜 , 𝑑𝑢 = 𝑎𝑑𝑡 → = 𝑑𝑡
𝑎
𝑔
(𝑎𝑡 + 𝑟𝑜 )3 𝑣 = 4𝑎 (𝑎𝑡 + 𝑟𝑜 )4 +c
Despejando v tenemos
𝑔
𝑣= (𝑎𝑡 + 𝑟𝑜 ) + 𝐶(𝑎𝑡 + 𝑟𝑜 )−3
4𝑎
𝑔𝑟𝑜 𝑔(𝑟𝑜 )4
𝑣(0) = 0 . 0 = + 𝐶(𝑟𝑜 )−3 , 𝐶=−
4𝑎 4𝑎
𝑔 𝑔(𝑟𝑜 )4 𝑘
𝑣= (𝑎𝑡 + 𝑟𝑜 ) − , 𝑎=
4𝑎 4𝑎(𝑎𝑡 + 𝑟𝑜 )3 𝜌
𝛒𝐠 𝐤𝐭 𝛒𝐠(𝐫𝐨 )𝟒
𝐯(𝐭) = ( + 𝐫𝐨 ) − 𝟑
𝟒𝐤 𝛒 𝐤𝐭
𝟒𝐤 ( 𝛒 + 𝐫𝐨 )
t(s)
𝛒𝐠 𝐤𝐭 𝛒𝐠(𝐫𝐨 )𝟒
𝐯(𝐭) = ( + 𝐫𝐨 ) − 𝟑
𝟒𝐤 𝛒 𝐤𝐭
𝟒𝐤 ( 𝛒 + 𝐫𝐨 )
r(t) = 0 ¿ ?
k
r(t) = ( ) t + ro
ρ
k
0.007 = ( ) 10 + 0.01
ρ
k 0.007 − 0.01
= = −0.0003
ρ 10
r(t) = 0
0 = −0.0003t + 0.01
−0.01
t = −0.0003 = 33.333 segundos.
P(t)
t
1. Hacemos grafico para una población PO = 100 𝑝𝑒𝑟𝑠𝑜𝑛𝑎𝑠 y un intervalo
de 𝑘 ∈ [0,10]
t
2. Hacemos grafico para una población PO = 500 𝑝𝑒𝑟𝑠𝑜𝑛𝑎𝑠 y un intervalo
de 𝑘 ∈ [1,10]
→ P(t) = 500 ∗ eKsen(t) , 𝑘 ∈ [0,10]
Manipulate[Plot[100*E^(k*Sin[t]), {t, 0, 10*Pi}], {k, 0,10}]
P
900
800
700
600
500
400
300
t
5 10 15 20 25 30
EJERCICIO 3.3
𝜋
53) Verifique si 𝑦 = sin ℎ𝑥 − 2 cos (𝑥 + 6 ) es una solución particular de
𝑦 (4) − 𝑦 = 0. Haga coincidir esta solución particular con la solución
general de la ecuación diferencial.
Solución
𝑒 𝑥 − 𝑒 −𝑥 𝜋 𝜋
𝑦= − 2 [cosx (cos ) − sin 𝑥 (sin )]
2 6 6
𝑒 𝑥 −𝑒 −𝑥 √3 1
𝑦= 2
− 2 [cosx ( 2 ) − sin 𝑥 (2)]
1 1
𝑦 = 𝑒 𝑥 − 𝑒 −𝑥 − √3 cos 𝑥 + sin 𝑥
2 2
𝑦" + 𝜆𝑦 = 0
( )=0
𝑃𝐹𝑉 = { 𝑦 0𝜋
𝑦( ) = 0
2
Solución
Ecuación auxiliar:
𝑟 2 +∝= 0 , considerando tres casos donde:
𝜆 = 0, 𝜆 =∝2 > 0, , 𝜆 = −∝2 < 0
CASO I: ∝= 0 → 𝑦" = 0
→ 𝑦 = 𝑐1 + 𝑐2 𝑥
Para : 𝑦(0) = 0 → 𝑐1 = 0
𝜋
𝑦 ( ) = 0 → 𝑐2 = 0
2
∴el problema posee una sola solución trivial
CASO II: 𝜆 = −∝2 < 0 → 𝑦" −∝2 𝑦 = 0
Ecuación auxiliar:
𝑟 2 −∝= 0 → 𝑟 = ±∝.
→ 𝑦 = 𝑐1 𝑒 ∝𝑥 + 𝑐2 𝑒 −∝𝑥 𝑜
𝑦 = 𝑐1 cos ℎ(∝ 𝑥) + 𝑐2 sin ℎ(∝ 𝑥)
Para: : 𝑦(0) = 0 → 𝑐1 = 0
→ 𝑦 = 𝑐2 sin ℎ( ∝ 𝑥).
Para:
𝜋
𝑦 ( ) = 0 → 𝑐2 = 0
2
𝜆 = 𝛼 2 > 0 ⟶ 𝑌 ′′ + 𝛼 2 𝑌 = 0
Ec. Auxiliar: 𝑟 2 + 𝛼 2 = 0 → 𝑟 = ±𝛼𝑖
𝑦 = 𝑐1 cos 𝛼𝑥 + 𝑐2 sin 𝛼𝑥
Para 𝑦(0) = 0 → 𝑐1 = 0
𝑦 = 𝑐2 sin 𝛼𝑥
𝜋
Para 𝑦 (2 ) = 0, 𝜆 = 𝛼 2 = 2𝑘
→ 𝑦 = 𝑐2 sin 2𝑘 𝑥
∴ 𝐸𝑙 𝑝𝑟𝑜𝑏𝑙𝑒𝑚𝑎 𝑡𝑒𝑛𝑑𝑟𝑎 𝑠𝑜𝑙𝑢𝑐𝑖𝑜𝑛𝑒𝑠 𝑛𝑜 𝑡𝑟𝑖𝑣𝑖𝑎𝑙𝑒𝑠
EJERCICIO 3.3
6)Una fuerza de 400 newton alarga 2 metros un resorte. Una masa de 50
kilogramos se une al extremo del resorte y se libera inicialmente desde
la posición de equilibrio con una velocidad ascendente de 10 m/s.
encuentre la ecuación de movimiento.
Sabemos que:
m.x``(t) + k.x(t)=0
pero; F=k.x → 400=k.2
K=200
m=50
entonces: 50.x``(t) + 200.x(t)=0
hallamos la solución general:
m2 + 4=0 → m=±2𝑖
entonces la solución general sería:
x(t)= c1.cos(2t) + c2.sen(2t)
derivando tendríamos:
x`(t)= -2.c1.sen(2t) + 2c2.cos(2t)
50. x``(t) + 200. x(t) = 0
{ 𝑥 (0) = 0
𝑥`(0) = −10
→ reemplazamos valores iniciales:
0 = c1.cos(0) + c2.sen(0)
0 = c1
-10 = -2.c1.sen(0) + 2c2.cos(0)
-5 = c2
→la ecuación del movimiento sería:
x(t)= -5.sen(2t)
Gráfica en MATHEMATICA del ejercicio.
𝑑𝑥
37) 𝑑𝑡 + 4𝑥 = −5𝑠𝑒𝑛2𝑡 + 3𝑐𝑜𝑠2𝑡 , x(0) = -1 , x`(0) = 1
𝑧1 𝑧2 𝑧2 ∗ f(t)∗dt 𝑧1 ∗ f(t)∗dt
W = |𝑧` 𝑧`2 | ; u1 = − ∫ u2 = ∫
1 𝑤 𝑤
5sen2t∗sen2t∗dt 5sen2t∗cos2t∗dt
u1 = ∫ u2 = − ∫
2 2
5 (1−𝑐𝑜𝑠4𝑡)∗𝑑𝑡 −5 𝑠𝑒𝑛4𝑡∗𝑑𝑡
u1 = 2 ∫ 2
u2 = 2
∫ 2
5 𝑠𝑒𝑛4𝑡 5
u1 = 4 (𝑡 − ) u2 = * cos4t
4 16
𝑧1 𝑧2 𝑧2 ∗ f(t)∗dt 𝑧1 ∗ f(t)∗dt
W = |𝑧` 𝑧`2 | ; u1 = − ∫ u2 = ∫
1 𝑤 𝑤
sen2t∗3cos2t∗dt 3cos2t∗cos2t∗dt
u1 = − ∫ u2 = ∫
2 2
−3 (𝑠𝑒𝑛4𝑡)∗𝑑𝑡 𝑐𝑜𝑠2𝑡∗𝑐𝑜𝑠2𝑡∗𝑑𝑡
u1 = ∫ u2 = 3 ∫
2 2 2
3 3 𝑠𝑒𝑛4𝑡
u1 = 16
(𝑐𝑜𝑠4𝑡) u2 = 4*(t + 4
)
Entonces la ecuación particular obtenida para z2 sería:
Zp = u1.cos(2t) + u2.sen(2t)
3 3 𝑠𝑒𝑛4𝑡
Zp = (𝑐𝑜𝑠4𝑡).cos(2t) + *(t + ).sen(2t)
16 4 4
Derivando obtenemos:
1
X`(t) = 16
[20 𝑐𝑜𝑠2𝑡 − 40𝑡. 𝑠𝑒𝑛2𝑡 − 20𝑐𝑜𝑠4𝑡. 𝑐𝑜𝑠2𝑡 + 10𝑠𝑒𝑛4𝑡. 𝑠𝑒𝑛2𝑡 −
1
12𝑠𝑒𝑛4𝑡. 𝑐𝑜𝑠2𝑡 − 6𝑐𝑜𝑠4𝑡. 𝑠𝑒𝑛2𝑡 − 32. 𝑐1. 𝑠𝑒𝑛2𝑡] + 16
[12𝑠𝑒𝑛2𝑡 + 24𝑡. 𝑐𝑜𝑠2𝑡 +
12𝑐𝑜𝑠4𝑡. 𝑠𝑒𝑛2𝑡 + 6𝑠𝑒𝑛4𝑡. 𝑐𝑜𝑠2𝑡 − 20𝑠𝑒𝑛4𝑡. 𝑠𝑒𝑛2𝑡 + 10𝑐𝑜𝑠4𝑡. 𝑐𝑜𝑠2𝑡 + 32. 𝑐2 . 𝑐𝑜𝑠2𝑡]
Reemplazando los valores iniciales obtenemos:
-1 = (0 – 0 + 3/16 + c1)
→ c1 = -19/16
1 = 1/16∗ [20 − 20] + 1/16∗ [10 + 32 ∗ 𝑐2 ]
1 = 1/16∗ [10 + 32 ∗ 𝑐2 ]
3/16 = c2
Finalmente encontramos la solución particular:
−19 3 5 𝑠𝑒𝑛4𝑡 5
X(t) = ∗ cos(2t) + ∗sen(2t) + (𝑡 − ).cos(2t) +( * cos4t).sen(2t) +
16 16 4 4 16
3 3 𝑠𝑒𝑛4𝑡
(𝑐𝑜𝑠4𝑡).cos(2t) + *(t + ).sen(2t)
16 4 4
Gráfica en MATHEMATICA del ejercicio.
EJERCICIO 5.1.1
34) 𝑦 " + 𝑒 𝑥 𝑦 , + 𝑦 = 0
Solución
∞
𝑛=2
1 1
+ (1 + 𝑥 + 𝑥 2 + 𝑥 3 + ⋯ ) (𝑐1 + 2𝑐2 + 3𝑐3 𝑥 2 + 4𝑐4 𝑥 3 + ⋯ )
2 6
∞
− ∑ 𝑐𝑛 𝑥 𝑛
𝑛=0
Igualando a cero :
1 1
2𝑐2 + 𝑐1 − 𝑐0 = 0 → 𝑐2 = 𝑐0 − 𝑐1
2 2
1
6𝑐3 + 2𝑐2 = 0 → 𝑐3 = − 𝑐2
3
1 1 1 1
12𝑐4 + 3𝑐3 + 𝑐1 = 0 → 𝑐4 = − 𝑐3 + 𝑐2 − 𝑐1
2 4 12 24
Para: 𝑐0 = 1 y 𝑐1 = 0 , tenemos:
1 1 1
𝑐2 = − 2 , 𝑐3 = 6 , 𝑐4 = − 24
Plot[𝑥 − 1⁄2 ∗ 𝑥^2 + 1⁄6 ∗ 𝑥^3 − 1⁄24 ∗ 𝑥^4, {𝑥, −∞, +∞}]
35)sin resolver realmente la ecuación diferencial cos 𝑥 𝑦 " + 𝑦 , + 5𝑦 = 0,
encuentre un límite menor para el radio de convergencia de las
soluciones en forma de series de potencias en torno a x=0.en torno a
x=1.
Solución
𝜋
Los puntos singulares de cos 𝑥 𝑦 " + 𝑦 , + 5𝑦 = 0 son múltiplos enteros de 2
𝜋 𝜋
La distancia de 0 para cualquier ± 2 es 2 .
𝜋
El punto singular más cercano a 1 es
2
𝜋
La distancia más cercana de 1 al punto singular es entonces − 1.
2