Reporte2 Tercerparcial

Instituto Politécnico Nacional
Unidad Profesional Interdisciplinaria de Biotecnología
Métodos numéricos
Reporte No. 2 tercer parcial ‘Aproximación a la solución de

ecuaciones diferenciales’
Profesores:
Diana Bueno Hernández
Jesús Granados Hernández
Integrantes:
Castillo Hernández Itzel - Ing. Biotecnológica

Fragoso Jiménez Daniela Jazmín - Ing. Biotecnológica
Hernández Castro Andrea Daniela - Ing. Biotecnológica
Honorato García Carolina - Ing. Biotecnológica
Nava Vargas Andrea Sarahi - Ing. Ambiental
Fecha: 12/11/2019
INTRODUCCIÓN:
Son muchas las fuentes de ecuaciones diferenciales, algunas de ellas son los
problemas físicos. Donde, por ejemplo, es necesario conocer la posición de una
partícula sumergida en un medio viscoso en un tiempo determinado, donde dicha
posición dependen tanto de la velocidad inicial de la partícula como de la fricción
producida por el medio, que depende de la velocidad de la misma. Así mismo, son
muy pocos los modelos que llevan a soluciones algebráicas exáctas. Es necesario
recurrir a métodos que permitan obtener una solución aproximada, donde existirá
un error, que dependerá tanto del algoritmo utilizado como de parámetros
computacionales.
La interpretación aproximada de algunos fenómenos naturales ha llevado a
introducir en ciertos tipos de ecuaciones diferenciales cambios en la variable
temporal llamados retardos, lo cual hace que dichas ecuaciones y sus soluciones
tengan un comportamiento más acorde con la realidad. Estas ecuaciones,
denominadas ecuaciones diferenciales con retardo requieren de métodos complejos
para su solución y en la mayoría de casos, tan sólo se logra realizar una
aproximación numérica.
CLASE:
Toma de XB Excreción del XB
C C
Para determinar el modelo de unatexperimentación en la que sólo tienes los datos

t
obtenidos, se pueden seguir los siguientes pasos:
1) Modelo empírico lineal
𝑑𝐶
= 𝑎𝑡 + 𝑏
𝑑𝑡
2) Hallar la derivada
Hacia adelante
𝑓(𝑥𝑖 + 1) − 𝑓(𝑥𝑖)
𝐹 ′ (𝑥𝑖) =
(𝑥𝑖 + 1) − 𝑥𝑖
h
x
xi-1
i
xi+1
h
En donde:
(xi+1) – xi=h
xi - (xi-1)= h
(xi-1) – xi= -h
Centrada
𝑓(𝑥𝑖 + 1) − 𝑓(𝑥𝑖 − 1)
𝐹 ′ (𝑥𝑖) =
(𝑥𝑖 + 1) − (𝑥𝑖 − 1)
Hacia atrás
𝑓(𝑥𝑖 − 1) − 𝑓(𝑥𝑖)
𝐹 ′ (𝑥𝑖) =
(𝑥𝑖 − 1) − 𝑥𝑖
3) Cambio de variable
𝑑𝐶
= 𝑌; Y= at+b si t=X entonces: Y=aX+b
𝑑𝑡
4) Ajuste por mínimos cuadrados y encontrar a y b

5) Sustituir en el modelo original
Ejemplo:
Administramos por vía intravenosa 1.5 gramos de un determinado fármaco, es decir,
1500 mg. Las cantidades (mg) y concentración (mg/L) son las siguientes en los
tiempos que se indican:
Tiempo (h) Concentración
(mg/L)
4 150
6 100
8 80
10 50
12 30
14 20
16 10
18 7
20 6
22 5
24 4
El modelo empírico es una ecuación diferencial que expresa la concentración
remanente en el organismo (cinética de orden 1)
𝑑𝐶
= 𝛽𝐶 + 𝛼
𝑑𝑡
a) Determine las constantes α y β
t=[4 6 8 10 12 141 6 18 20 22 24] o t=4:2:24
C=[150 100 80 50 30 20 10 7 6 5 4 ]
Donde esté X en las ecuaciones de derivadas hacia adelante, atrás y centrada,
pondrás los valores de t y dónde esté f se pondrán las concentraciones.
I=1
Y(1)=C(i+1) – C(¡) / ((t(i+1) - t(i) %hacia adelante
For i=2:(length(t-1)) %centrado
Y(i)=(C(i+1) – C(i-1) / (t(i+1) – t(i-1))
End
I=length(t) %atrás
Y(i)= (C(i-1)-C(i)) / ((t(i-1) – t(i))
𝑑𝐶
%modelo= 𝑑𝑡 = 𝑌 = 𝛽 + 𝛼𝐶
Encontrar α y β por medio de matriz inversa

B=C.^0
B(2,:)=C
Sol=inv(B*B’)*(B*Y’)
ALFA=Sol(1)
BETA=Sol(2)
Plot(C,Y,’*r’)
Syms C1
F=Sol(1,1)+Sol(2,1)*C1
T1=4:0.1:24
Y1=subs(f,t1)
Hold on
Plot(t1,Y1)
ALFA=-0.8398
BETA=-0.1656
%MODELO Y=-0.8398 – 0.1656*C
APROCIMACIÓN DE ECUACIONES DIFERENCIALES
𝑑𝐶
= −0.8398 – 0.1656 ∗ 𝐶
𝑑𝑡
Syms tie
C2=vpa(dsolve(‘Dc=-0.1656 * C -0.8398,’C(4)=150’,tie))
C2=300.8 + Exp(-0.1656*t)--5.071 %solución de la ecuación diferencial
Figure
Plot(t,C,’*r’)
Hold on
t2=4:0.1:24
Yana=subs(c2,t2)
MÉTODO DE EULER
𝑑𝐶
= −0.8398 – 0.1656 ∗ 𝐶
𝑑𝑡
Si sólo se confieren un por de coordenadas (x,y), es posible determinar la ecuación
de una curva
Si tienes Y1 y X1
𝑑𝑌
𝑌= ∗ (𝑥 − 𝑥1) ∗ 𝑌1
𝑑𝑋
Y=f(x,y)(x.-x1) + Y1 %ecuación de la recta
Y2= F(x2-x1)+Y1 %evaluación en x1
Es una aproximación con rectas y puntos…
Y=f(x2,y2)(x-x2)+Y2
Y3=f(x2,y2)(x3,x2)+Y2
Y=f(x3,y3)(x-x3)+Y3
Y4=F(x3,y3)(x4-x3)+Y3
En general Yn=f(Xn-1, Yn-1)(Xn-Xn-1)+Yn-1
Aproximación a la solución de la recta fe la ecuación diferencial
CONCLUSIONES:
Al tener sólo datos de experimentación es importante lograr predecir los fenómenos
que se presentan en esta, por lo mismo, es importante determinar la ecuación del
modelo de experimentación esto puede hacerse con ecuaciones diferenciales.
REFERENCIAS:
http://repositorio.utp.edu.co/dspace/bitstream/handle/11059/7192/51535T866.
pdf?sequence=1&isAllowed=y
http://www.dmae.upct.es/~jose/metodos/numer4.pdf

Reporte2 Tercerparcial

Cargado por

Copyright:

Formatos disponibles

Reporte2 Tercerparcial

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Reporte2 Tercerparcial

Cargado por

Copyright:

Formatos disponibles

Instituto Politécnico Nacional

Unidad Profesional Interdisciplinaria de Biotecnología

Reporte No. 2 tercer parcial ‘Aproximación a la solución de

Diana Bueno Hernández

Jesús Granados Hernández

Castillo Hernández Itzel - Ing. Biotecnológica

Hernández Castro Andrea Daniela - Ing. Biotecnológica

Honorato García Carolina - Ing. Biotecnológica

Nava Vargas Andrea Sarahi - Ing. Ambiental

Para determinar el modelo de unatexperimentación en la que sólo tienes los datos

4) Ajuste por mínimos cuadrados y encontrar a y b

Encontrar α y β por medio de matriz inversa

También podría gustarte