DIFERENCIACION NUMERICA


DIFERENCIACIÓN NUMÉRICA
Se desarrollarán fórmulas para aproximaciones de diferencias hacia delante, hacia atrás y centradas para la primera derivada utilizando la serie truncada de Taylor.
En el mejor de los casos, estas estimaciones presentan errores de orden O(h2); es decir, sus errores fueron proporcionales al cuadrado de su tamaño de paso. Este nivel de exactitud se debe al número de términos de la serie de Taylor.
Diferencias Divididas Finitas, 1ª derivada:
Definiendo un tamaño de paso h = xi+1 – xi
- Hacia delante:
f(x_{i + 1} ) = f(x_i ) + f'(x_i )h + {{f''(x_i )} \over {2!}}h^2 + {{f'''(x_i )} \over {3!}}h^3 + \ldots
f'(x_i ) \simeq {{f(x_{i + 1} ) - f(x_i )} \over h} + O(h)
- Hacia atrás:
f(x_{i - 1} ) = f(x_i ) - f'(x_i )h + {{f''(x_i )} \over {2!}}h^2 - {{f'''(x_i )} \over {3!}}h^3 + \ldots
f'(x_i ) \simeq {{f(x_i ) - f(x_{i - 1} )} \over h} + O(h)
- Central:
f'(x_i ) \simeq {{f(x_{i + 1} ) - f(x_{i - 1} )} \over {2h}} + O(h^2 )
Diferenciación de fórmulas de alta exactitud:
Se pueden generar fórmulas de alta exactitud al incluir términos adicionales en la expansión de la serie de Taylor. Teniendo en cuenta el término de la segunda derivada:
 f(x_{i + 1} ) = f(x_i ) + f'(x_i )h + {{f''(x_i )} \over {2!}}h^2 + \ldots
Despejando para la primera derivada:
 f(x_{i + 1} ) = f(x_i ) + f'(x_i )h + {{f''(x_i )} \over {2!}}h^2 + \ldots
De la expansión de Taylor hacia delante para f(xi+2) en términos de f(xi):
f''(x_i ) = {{f(x_{i + 2} ) - 2f(x_{i + 1} ) + f(x_i )} \over {h^2 }} + O(h) 
Despejando para f’’(xi):
f'(x_i ) = {{f(x_{i + 1} ) - f(x_i )} \over h} - {{f(x_{i + 2} ) - 2f(x_{i + 1} ) + f(x_i )} \over {2h^2 }}h
Agrupando términos y reordenando:
f'(x_i ) = {{ - f(x_{i + 2} ) + 4f(x_{i + 1} ) - 3f(x_i )} \over {2h}} + O(h^2 ) 
La inclusión del término de la segunda derivada mejoró la exactitud en O(h2).
Se pueden desarrollar fórmulas similares para las fórmulas centradas y hacia atrás.
Referencias:
Este módulo fue desarrollado por Diego López, usando notas de libro:
HEATH, Michael; Scientific Computing: An introductory survey. McGraw Hill. 1997. Capítulo 8. Página 262.
Al igual que en la integral, a veces resulta muy complicado hayar la derivada de una función. Usando la serie truncada de Taylor, se pueden desarrollar fórmulas aproximadas iterativas para resolver las derivadas numéricamente usando n puntos.
Usando hasta la expresión de Taylor de la segunda derivada se puede hayar mayor presición para la derivada de la función.









Comentarios

Publicar un comentario

Entradas más populares de este blog

unidad 6

Métodos de un paso Los métodos de un paso tienen por objetivo obtener una aproximación de la solución de un problema bien planteado de valor inicial en cada punto de la malla, basándose en el resultado obtenido para el punto anterior. Se desarrollan aquí los métodos Taylor (incluyendo Euler), y de Runge Kutta. Para ver el detalle de cada uno de los métodos, hacer click en cada uno de los siguientes vínculos. Para volver a esta página, hacer click en la solapa "métodos de un paso". Consistencia, estabilidad y convergencia Hay algunas propiedades importantes de las ecuaciones en diferencias para problemas de valor inicial de EDOs de primer orden que deben considerarse antes de que se pongan en práctica los métodos numéricos. Ellas son consistencia, estabilidad y convergencia. Métodos de pasos múltiples Se considera el problema de valores iniciales (P.V.I.) 8<: y0(x) = f(x; y(x)); x 2 [a; b]; y(a) = y0 dado, el que supondremos tiene solución única, y : [a; b] ...
Leer más