´ MATEMATICAS AVANZADAS ACTIVIDADES CLASE No. 3 Luis Miguel P´erez Pertuz Nelson Forero Salcedo February 27, 2016 Part
Views 121 Downloads 2 File size 319KB
´ MATEMATICAS AVANZADAS ACTIVIDADES CLASE No. 3 Luis Miguel P´erez Pertuz Nelson Forero Salcedo February 27, 2016
Parte A : Ordinary Differential Equations (ODEs) Chapter 1: First-Order ODEs 1. SET 1.5. Problema 9, 28. 2. SET 1.5. Modeling. Escoja uno cualquiera de los problemas desde el 32 hasta el 40. Realice lo siguiente: (a) Encuentre la ecuaci´on diferencial que lo modela. (b) Encuentre una soluci´on general anal´ıtica. (c) Encuentre una soluci´on particular anal´ıtica. 3. Para los tres ejercicios anteriores use Matlab o cualquier otra herramienta para: (a) Graficar un campo de direcciones. (b) Encontrar una soluci´on particular en cualquier punto. (c) Evaluar la soluci´ on particular en cualquier punto. (d) Graficar la soluci´ on particular.
1
1. Problem Set 1.5 3-13 Soluci´ on general. Problemas de valores iniciales. Encuentre la soluci´ on general. Si se da una condici´on inicial, hallar la soluci´on particular correspondiente y graficala o bosquejala. 9. y 0 + ysin(x) = ecos(x) , y(0) = −2.5 Soluci´ on: y 0 + ysin(x) = ecos(x) Esta es una ecuaci´ on lineal de primer orden de la forma: dy dx
+ P (x)y = f (x)
En el libro de ecuaciones diferenciales del autor Dennis Zill s´eptima edici´on en la pagina 53 capitulo 2.3. Se define la soluci´on a las ecuaciones diferenciales de esta forma. Seg´ un este libro la soluci´on de esta ecuaci´on est´a determinado por: R R R R y = ce− P (x)dx + e− P (x)dx e P (x)dx f (x)dx. Por lo tanto tenemos: p(x) = sen(x) f (x) =Recos(x) R − P (x)dx ==> − sen(x)dx = cos(x) reemplazando en la soluci´ on tenemos: R R R R P (x)dx y = ce− P (x)dx + e− P (x)dx e f (x)dx. R cos(x) e dx. y = cecos(x) + ecos(x) e−cos(x) R y = cecos(x) + ecos(x) e( 0)dx. y = cecos(x) + ecos(x) x. y = ecos(x) (c + x). Aplicando condiciones iniciales: y(0) = −2.5 y = ecos(x) (c + x). −2.5 = ecos(0) (c + 0). C=
−2.5 e
Soluci´ on Particular:
2
y = ecos(x) (x −
−2.5 e ).
Figure 1: Gr´ afica de la soluci´on particular y = ecos(x) (x −
3
−2.5 e ).
22-28 EDO’s no lineares Usando un m´etodo de esta secci´on o variables separables, encuentre una soluci´on general. Si se da una condici´on inicial, encuentre la soluci´on particular y bosqueja o graf´ıcala. 28. 2xyy 0 + (x − 1)y 2 = x2 ex ; (set y 2 = z) Soluci´ on:
2xyy 0 + (x − 1)y 2 = x2 ex y2 = Z 0 y 0 = 2z√z √ 0 2x z 2z√z + (x − 1)z = x2 ex xz 0 + (x − 1)z = x2 ex z0 +
(x−1 x z
= xex
Esta es una ecuaci´ on lineal de primer orden de la forma: dy dx
y = ce−
R
P (x)dx
+ P (x)y = f (x)
+ e−
R
P (x)dx
R R e
P (x)dx
f (x)dx.
utilizando la soluci´ on de la forma de Bernoulli. p(x) = (x−1 x f (x) = xex R y = ce−
p(x)dx = R
(x−1 x dx
R
(x−1 2x dx
+ e−
R
= (x − ln(x))
(x−1 x dx
z = ce−(x−ln(x)) + e−(x−ln(x)) z = ce−x+ln(x)) + e−x+ln(x))
R R e
z = ce−x+ln(x)) + e−x+ln(x))
R
e(x−ln(x)) xex dx e(2x−ln(x)) xdx
z = ce−x+ln(x)) + e−x+ln(x)) 4
xex dx
e(x−ln(x)) xex dx
R
R
(x−1 x dx
R
e(2x) xdx eln(x)
z = ce−x+ln(x)) + e−x+ln(x))
R
e(2x) dx (2x)
z = ce−x+ln(x)) + e−x+ln(x)) e 2 (2x)
z = cxe−x + xe−x e 2 z = xe−x (c +
e(2x) 2 ) (2x)
z = xe−x ( c+e2
(2x)
z = xe−x ( c+e2
) )
(2x)
y 2 = xe−x ( c+e2 y=
q
)
(2x)
xe−x ( c+e2
)
teniendo la soluci´ on general hacemos y(0)=0 para hallar la constante C. Haciendo y = 0 y y x = 0 C = −e2x C = −e2(0) = -1 C = -1 q (2x) ) y = xe−x ( −1+e 2
5
2. Problem Set 1.5 31-40 Modelado. Aplicaciones adicionales. 33.Inyecci´ on de drogas encuentre y resuelva el modelo para la inyecci´on de droga en el torrente sangu´ıneo si, en t=0, una cantidad constante A (g/min) es inyectado y la droga es simult´aneamente removida a una rata proporcional a la cantidad de droga presente en un tiempo t. Soluci´ on: consideremos y= a la cantidad de droga que entra en el organismo. por lo tanto la raz´ on de cambio de esta queda definida por la ecuaci´on diferencial:
y 0 = A − ky como condici´ on inicial tenemos que en un tiempo t=0 =¿ y(0)=0.
y 0 = A − ky (a) y0 A−ky dy dt
A−ky
R
=1 =1
dy A−ky
R = dt u = −ky + A du du = −kdy → −k = dy 1 −k
R
du u
1 −k ln(u)
=
R
dt
=t+c
ln(u) = −kt + c u = e−kt+c −ky + A = e−kt+c −ky + A = ce−kt y=
ce−kt −A −k
6
y=
Ace−kt −A −k
A(1−ce−kt ) k
y=
(b)
Evaluando condici´ on inicial: y(0)=0 0=
A(1−ce−k(0) ) k
0=
A(1−c) k
0=1−c c=1 La soluci´ on particular es: y=
A(1−e−kt ) k
7
(c)
3. Problem Set 1.5 Gr´ aficas de campos de direcciones y soluciones particulares
Figure 2: Gr´ afica del campo de direcciones con soluciones particulares de la ED: y 0 + ysin(x) = ecos(x)
8
Figure 3: Gr´ afica del campo direccional con soluciones particulares de la ED: 2xyy 0 + (x − 1)y 2 = x2 ex
9
Figure 4: Gr´ afica del campo de direcciones y soluciones particulares de ED: y 0 = A − ky donde A = 10(g/min) y k = 1 con valores de C distintos
10