Ejemplo Diferencias Finitas 1D
Deisy Chaves
[Link]@[Link]
Estrategia general para la solución de EDPs
Estrategia general para la solución de EDPs:
1. Evaluar la solución (variables independientes) en un número
finito de puntos (malla o grilla)
Estrategia general para la solución de EDPs
Estrategia general para la solución de EDPs:
1. Evaluar la solución (variables independientes) en un número
finito de puntos (malla o grilla)
2. Aproximar derivadas en termino de los valores conocidos en los
nodos/puntos de la grilla
Estrategia general para la solución de EDPs
Estrategia general para la solución de EDPs:
1. Evaluar la solución (variables independientes) en un número
finito de puntos (malla o grilla)
2. Aproximar derivadas en termino de los valores conocidos en los
nodos/puntos de la grilla
3. Reemplazar derivadas temporales por alguna de las
discretizaciones estudiadas (Derivación basada en desarrollos
de la serie de Taylor o Diferencias Finitas)
Estrategia general para la solución de EDPs
Estrategia general para la solución de EDPs:
1. Evaluar la solución (variables independientes) en un número
finito de puntos (malla o grilla)
2. Aproximar derivadas en termino de los valores conocidos en los
nodos/puntos de la grilla
3. Reemplazar derivadas temporales por alguna de las
discretizaciones estudiadas (Derivación basada en desarrollos
de la serie de Taylor o Diferencias Finitas)
4. Reducir la(s) ecuación(es) diferencial(es) parcial(es) y
condiciones de borde a un sistema de ecuaciones algebraicas
Estrategia general para la solución de EDPs
Métodos tradicionales para invertir la matriz A asociada a un
Sistema de ecuaciones (e.g. Eliminación Gausiana) requieren del
orden de O(N3) operaciones y requieren almacenar todos los
coecientes de la matriz (O(N2)).
La matriz de coeficientes asociada al sistema de ecuaciones que
da solución al EDP es:
Estrategia general para la solución de EDPs
Métodos tradicionales para invertir la matriz A asociada a un
Sistema de ecuaciones (e.g. Eliminación Gausiana) requieren del
orden de O(N3) operaciones y requieren almacenar todos los
coecientes de la matriz (O(N2)).
La matriz de coeficientes asociada al sistema de ecuaciones que
da solución al EDP es:
o sparse o con la mayoría de los elementos iguales a cero.
o Existen varios métodos para resolver este tipo de sistemas que
aprovechan la estructura de la matriz para ahorrar memoria
(solo los coecientes que no son ceros son almacenados o para
acelerar el computo).
Estrategia general para la solución de EDPs
La matriz de coeficientes asociada al sistema de ecuaciones que
da solución al EDP es:
o sparse o con la mayoría de los elementos iguales a cero.
Ejercicio
Halle la solución de la ecuación diferencial
𝑑𝑑𝑑𝑑
= sin(𝑥𝑥 2 )
𝑑𝑑𝑑𝑑
Ejercicio
Halle la solución de la ecuación diferencial
𝑑𝑑𝑑𝑑
= sin(𝑥𝑥 2 )
𝑑𝑑𝑑𝑑
Se determina la malla para discretizar y hallar la solución en un
número finito de n + 1 puntos.
Ejercicio
Halle la solución de la ecuación diferencial
𝑑𝑑𝑑𝑑
= sin(𝑥𝑥 2 )
𝑑𝑑𝑑𝑑
Se determina la malla para discretizar y hallar la solución en un
número finito de n + 1 puntos.
1
ℎ=
𝑛𝑛
Ejercicio
Reemplazar derivadas temporales por alguna de las
discretizaciones estudiadas (Derivación basada en desarrollos
de la serie de Taylor o Diferencias Finitas)
𝑑𝑑𝑑𝑑
= sin(𝑥𝑥 2 )
𝑑𝑑𝑑𝑑
Ejercicio
Reemplazar derivadas temporales por alguna de las
discretizaciones estudiadas (Derivación basada en desarrollos
de la serie de Taylor o Diferencias Finitas)
𝑑𝑑𝑑𝑑
= sin(𝑥𝑥 2 )
𝑑𝑑𝑑𝑑
Derivadas de orden superior hacia adelante
Ejercicio
Reemplazar derivadas temporales por alguna de las
discretizaciones estudiadas (Derivación basada en desarrollos
de la serie de Taylor o Diferencias Finitas)
𝑑𝑑𝑑𝑑
= sin(𝑥𝑥 2 )
𝑑𝑑𝑑𝑑
Derivadas de orden superior hacia adelante
Ejercicio
Reemplazar derivadas temporales por alguna de las
discretizaciones estudiadas (Derivación basada en desarrollos
de la serie de Taylor o Diferencias Finitas)
𝑑𝑑𝑑𝑑
= sin(𝑥𝑥 2 )
𝑑𝑑𝑑𝑑
1
𝑥𝑥1 −3𝑓𝑓 𝑥𝑥1 + 4𝑓𝑓 𝑥𝑥2 − 𝑓𝑓(𝑥𝑥3 ) = sin(𝑥𝑥1 )2
2ℎ
Ejercicio
Reemplazar derivadas temporales por alguna de las
discretizaciones estudiadas (Derivación basada en desarrollos
de la serie de Taylor o Diferencias Finitas)
𝑑𝑑𝑑𝑑
= sin(𝑥𝑥 2 )
𝑑𝑑𝑑𝑑
1
𝑥𝑥1 −3𝑓𝑓 𝑥𝑥1 + 4𝑓𝑓 𝑥𝑥2 − 𝑓𝑓(𝑥𝑥3 ) = sin(𝑥𝑥1 )2
2ℎ
Ejercicio
Reemplazar derivadas temporales por alguna de las
discretizaciones estudiadas (Derivación basada en desarrollos
de la serie de Taylor o Diferencias Finitas)
𝑑𝑑𝑑𝑑
= sin(𝑥𝑥 2 )
𝑑𝑑𝑑𝑑
1
𝑥𝑥1 −3𝑓𝑓 𝑥𝑥1 + 4𝑓𝑓 𝑥𝑥2 − 𝑓𝑓(𝑥𝑥3 ) = sin(𝑥𝑥1 )2
2ℎ
1
𝑥𝑥1 −3𝑦𝑦1 + 4𝑦𝑦2 − 𝑦𝑦3 = sin(𝑥𝑥1 )2
2ℎ
Ejercicio
Reducir la(s) ecuación(es) diferencial(es) parcial(es) y
condiciones de borde a un sistema de ecuaciones algebraicas
𝑑𝑑𝑑𝑑
= sin(𝑥𝑥 2 )
𝑑𝑑𝑑𝑑
Ejercicio
Reducir la(s) ecuación(es) diferencial(es) parcial(es) y
condiciones de borde a un sistema de ecuaciones algebraicas
Ejercicio
Reducir la(s) ecuación(es) diferencial(es) parcial(es) y
condiciones de borde a un sistema de ecuaciones algebraicas
𝐴𝐴𝑦𝑦 = 𝑏𝑏
𝑦𝑦 = 𝐴𝐴−1 𝑏𝑏
