Una diferencia finita es una expresión matemática de la forma f ( x + b ) - f ( x + a ) . Si una diferencia finita se divide por b - a , se obtiene un cociente de diferencias . La aproximación de derivadas por diferencias finitas juega un papel central en los métodos de diferencias finitas para la solución numérica de ecuaciones diferenciales , especialmente problemas de valores en la frontera .
Ciertas relaciones de recurrencia se pueden escribir como ecuaciones en diferencias reemplazando la notación de iteración con diferencias finitas.
Hoy en día, el término "diferencia finita" a menudo se toma como sinónimo de aproximaciones de derivadas en diferencias finitas , especialmente en el contexto de métodos numéricos . [1] [2] [3] Las aproximaciones en diferencias finitas son cocientes en diferencias finitas en la terminología empleada anteriormente.
Las diferencias finitas fueron introducidas por Brook Taylor en 1715 y también se han estudiado como objetos matemáticos abstractos autónomos en trabajos de George Boole (1860), LM Milne-Thomson (1933) y Károly Jordan (1939). Las diferencias finitas tienen su origen en uno de los algoritmos de Jost Bürgi ( c. 1592 ) y funcionan con otros, incluido Isaac Newton . El cálculo formal de diferencias finitas puede verse como una alternativa al cálculo de infinitesimales . [4]
Tipos basicos
Normalmente se consideran tres tipos básicos: diferencias finitas hacia adelante , hacia atrás y centrales . [1] [2] [3]
Una diferencia hacia adelante es una expresión de la forma
Dependiendo de la aplicación, la separación h puede ser variable o constante. Cuando se omite, h se toma como 1: Δ [ f ] ( x ) = Δ 1 [ f ] ( x ) .
A diferencia hacia atrás utiliza los valores de la función en x y x - h , en lugar de los valores en x + h y x :
Finalmente, la diferencia central viene dada por
Relación con derivados
La diferencia finita se usa a menudo como una aproximación de la derivada, típicamente en la diferenciación numérica .
La derivada de una función f en un punto x está definida por el límite .
Si h tiene un valor fijo (distinto de cero) en lugar de aproximarse a cero, entonces el lado derecho de la ecuación anterior se escribiría
Por tanto, la diferencia directa dividida por h se aproxima a la derivada cuando h es pequeña. El error en esta aproximación se puede derivar del teorema de Taylor . Suponiendo que f es diferenciable, tenemos
La misma fórmula es válida para la diferencia hacia atrás:
Sin embargo, la diferencia central (también llamada centrada) produce una aproximación más precisa. Si f es dos veces diferenciable,
El principal problema [ cita requerida ] con el método de diferencia central, sin embargo, es que las funciones oscilantes pueden producir derivada cero. Si f ( nh ) = 1 para n impar y f ( nh ) = 2 para n par, entonces f ′ ( nh ) = 0 si se calcula con el esquema de diferencia central. Esto es particularmente problemático si el dominio de f es discreto. Ver también derivada simétrica
Los autores para quienes diferencias finitas significan aproximaciones en diferencias finitas definen las diferencias hacia adelante / hacia atrás / centrales como los cocientes dados en esta sección (en lugar de emplear las definiciones dadas en la sección anterior). [1] [2] [3]
Diferencias de orden superior
De manera análoga, se pueden obtener aproximaciones en diferencias finitas para derivadas de orden superior y operadores diferenciales. Por ejemplo, utilizando la fórmula de diferencia central anterior para f ′ ( x +h/2) y f ′ ( x - h/2) y aplicando una fórmula de diferencia central para la derivada de f ′ en x , obtenemos la aproximación de diferencia central de la segunda derivada de f :
- Central de segundo orden
Del mismo modo, podemos aplicar otras fórmulas de diferenciación de forma recursiva.
- Segundo orden hacia adelante
- Segundo orden al revés
Más generalmente, el n º orden hacia adelante, hacia atrás y centrales diferencias están dadas por, respectivamente,
- Hacia adelante
o para h = 1 ,
- Hacia atrás
- Central
Estas ecuaciones usan coeficientes binomiales después del signo de suma que se muestra como (n
yo) . Cada fila deltriángulodePascalproporciona el coeficiente para cada valor dei.
Tenga en cuenta que la diferencia central, para n impar , tendrá h multiplicado por no enteros. Esto suele ser un problema porque equivale a cambiar el intervalo de discretización. El problema se puede solucionar tomando el promedio de δ n [ f ] ( x - h/2) y δ n [ f ] ( x + h/2) .
Las diferencias hacia adelante aplicadas a una secuencia a veces se denominan transformada binomial de la secuencia y tienen varias propiedades combinatorias interesantes. Las diferencias hacia adelante se pueden evaluar utilizando la integral de Nörlund-Rice . La representación integral para estos tipos de series es interesante, porque la integral a menudo se puede evaluar usando técnicas de expansión asintótica o de punto silla ; por el contrario, la serie de diferencias progresivas puede ser extremadamente difícil de evaluar numéricamente, porque los coeficientes binomiales crecen rápidamente para n grandes .
La relación de estas diferencias de orden superior con los respectivos derivados es sencilla,
Las diferencias de orden superior también se pueden utilizar para construir mejores aproximaciones. Como se mencionó anteriormente, la diferencia de primer orden se aproxima a la derivada de primer orden hasta un término de orden h . Sin embargo, la combinación
aproxima f ′ ( x ) hasta un término de orden h 2 . Esto se puede probar expandiendo la expresión anterior en series de Taylor , o usando el cálculo de diferencias finitas, que se explica a continuación.
Si es necesario, la diferencia finita se puede centrar alrededor de cualquier punto mezclando diferencias hacia adelante, hacia atrás y centrales.
Granos de tamaño arbitrario
Usando álgebra lineal se pueden construir aproximaciones en diferencias finitas que utilizan un número arbitrario de puntos a la izquierda y un número (posiblemente diferente) de puntos a la derecha del punto de evaluación, para cualquier orden derivada. Esto implica resolver un sistema lineal tal que la expansión de Taylor de la suma de esos puntos alrededor del punto de evaluación se aproxime mejor a la expansión de Taylor de la derivada deseada. Estas fórmulas se pueden representar gráficamente en una cuadrícula hexagonal o en forma de diamante. [5]
Esto es útil para diferenciar una función en una cuadrícula, donde, a medida que uno se acerca al borde de la cuadrícula, debe muestrear cada vez menos puntos en un lado.
Los detalles se describen en estas notas .
La calculadora de coeficientes de diferencia finita construye aproximaciones en diferencias finitas para plantillas no estándar (e incluso no enteras) dadas una plantilla arbitraria y un orden derivado deseado.
Propiedades
- Para todos los k y n positivos
En ecuaciones diferenciales
Una aplicación importante de las diferencias finitas se encuentra en el análisis numérico , especialmente en las ecuaciones diferenciales numéricas , que apuntan a la solución numérica de ecuaciones diferenciales ordinarias y parciales . La idea es reemplazar las derivadas que aparecen en la ecuación diferencial por diferencias finitas que las aproximen. Los métodos resultantes se denominan métodos de diferencias finitas .
Las aplicaciones comunes del método de diferencias finitas se encuentran en las ciencias computacionales y las disciplinas de la ingeniería, como la ingeniería térmica , la mecánica de fluidos , etc.
Serie de Newton
La serie de Newton consta de los términos de la ecuación de diferencias progresivas de Newton , que lleva el nombre de Isaac Newton ; en esencia, es la fórmula de interpolación de Newton , publicada por primera vez en su Principia Mathematica en 1687, [6] es decir, el análogo discreto de la expansión continua de Taylor,
que se cumple para cualquier función polinomial f y para muchas (pero no todas) funciones analíticas (no se cumple cuando f es de tipo exponencial . Esto se ve fácilmente, ya que la función seno se desvanece en múltiplos enteros de; la serie de Newton correspondiente es idénticamente cero, ya que todas las diferencias finitas son cero en este caso. Sin embargo, está claro que la función seno no es cero). Aquí, la expresión
es el coeficiente binomial , y
es el " factorial descendente " o "factorial inferior", mientras que el producto vacío ( x ) 0 se define como 1. En este caso particular, hay una suposición de pasos unitarios para los cambios en los valores de x , h = 1 de la generalización a continuación.
Note la correspondencia formal de este resultado con el teorema de Taylor . Históricamente, esto, así como la identidad Chu-Vandermonde ,
(que siguen de él, y que corresponden al teorema del binomio ), se incluyen en las observaciones que maduraron al sistema de cálculo umbral .
Para ilustrar cómo se puede usar la fórmula de Newton en la práctica, considere los primeros términos de duplicar la secuencia de Fibonacci f = 2, 2, 4, ... Uno puede encontrar un polinomio que reproduzca estos valores, calculando primero una tabla de diferencias, y luego sustituyendo las diferencias que corresponden ax 0 (subrayado) en la fórmula de la siguiente manera,
Para el caso de pasos no uniformes en los valores de x , Newton calcula las diferencias divididas ,
la serie de productos,
y el polinomio resultante es el producto escalar , [7]
- .
En el análisis con números p -ádicos , el teorema de Mahler establece que la suposición de que f es una función polinomial puede debilitarse hasta la suposición de que f es meramente continua.
El teorema de Carlson proporciona las condiciones necesarias y suficientes para que una serie de Newton sea única, si es que existe. Sin embargo, en general, no existe una serie de Newton.
La serie de Newton, junto con la serie de Stirling y la serie de Selberg , es un caso especial de la serie de diferencias generales , todas las cuales se definen en términos de diferencias hacia adelante adecuadamente escaladas.
En una forma comprimida y ligeramente más general y en nodos equidistantes, la fórmula dice
Cálculo de diferencias finitas
La diferencia hacia adelante se puede considerar como un operador , llamado eloperador de diferencia , que asigna la funciónfaΔ h [ f ]. [8] [9] Este operador asciende a
donde T h es el operador de desplazamiento con el paso h , definido por T h [ f ] ( x ) = f ( x + h ) , e I es el operador de identidad .
La diferencia finita de órdenes superiores se puede definir de manera recursiva como Δn
h≡ Δ h (Δn - 1
hora) . Otra definición equivalente es Δn
h= [ T h - I ] n .
El operador de diferencia Δ h es un operador lineal , como tal, satisface Δ h [ αf + βg ] ( x ) = α Δ h [ f ] ( x ) + β Δ h [ g ] ( x ) .
También satisface una regla especial de Leibniz indicada anteriormente, Δ h ( f ( x ) g ( x )) = (Δ h f ( x )) g ( x + h ) + f ( x ) (Δ h g ( x )) . Declaraciones similares son válidas para las diferencias centrales y hacia atrás.
Al aplicar formalmente la serie de Taylor con respecto a h , se obtiene la fórmula
donde D denota el operador de la derivada del continuo, mapeando f a su derivada f ′ . La expansión es válida cuando ambos lados actúan sobre funciones analíticas , para h suficientemente pequeño . Por lo tanto, T h = e hD , y formalmente invirtiendo los rendimientos exponenciales
Esta fórmula es válida en el sentido de que ambos operadores dan el mismo resultado cuando se aplican a un polinomio.
Incluso para las funciones analíticas, no se garantiza que la serie de la derecha converja; puede ser una serie asintótica . Sin embargo, se puede utilizar para obtener aproximaciones más precisas para la derivada. Por ejemplo, retener los dos primeros términos de la serie produce la aproximación de segundo orden a f ′ ( x ) mencionada al final de la sección Diferencias de orden superior .
Las fórmulas análogas para los operadores de diferencia central y hacia atrás son
El cálculo de diferencias finitas está relacionado con el cálculo umbral de la combinatoria. Esta correspondencia notablemente sistemática se debe a la identidad de los conmutadores de las cantidades umbral con sus análogos continuos ( límites h → 0 ),
Por tanto, un gran número de relaciones diferenciales formales de cálculo estándar que involucran funciones f ( x ) se asignan sistemáticamente a análogos umbral de diferencias finitas que involucran f ( x−1
h) .
Por ejemplo, el análogo umbral de un monomio x n es una generalización del factorial descendente anterior ( símbolo k de Pochhammer ),
así que eso
de ahí la fórmula de interpolación de Newton anterior (haciendo coincidir los coeficientes en la expansión de una función arbitraria f ( x ) en tales símbolos), y así sucesivamente.
Por ejemplo, el seno umbral es
Como en el límite del continuo, la función propia de Δ h/h también resulta ser exponencial,
y por lo tanto, las sumas de Fourier de funciones continuas se asignan fácilmente a sumas umbral de Fourier fielmente , es decir, involucrando los mismos coeficientes de Fourier multiplicando estos exponenciales de base umbral. [10] Este umbral exponencial equivale a la función generadora exponencial de los símbolos de Pochhammer .
Así, por ejemplo, la función delta de Dirac se asigna a su correspondiente umbral, la función del seno cardinal ,
Etcétera. [11] Las ecuaciones en diferencias a menudo se pueden resolver con técnicas muy similares a las que se utilizan para resolver ecuaciones diferenciales .
El operador inverso del operador de diferencia hacia adelante, entonces la integral umbral, es el operador de suma indefinida o antidiferencia.
Reglas para el cálculo de operadores de diferencias finitas
De manera análoga a las reglas para encontrar la derivada , tenemos:
- Regla de la constante : si c es una constante , entonces
- Linealidad : si un y b son constantes ,
Todas las reglas anteriores se aplican igualmente bien a cualquier operador de diferencia, incluido ∇ como Δ .
- Regla de producto :
- Regla del cociente :
- o
- Reglas de suma :
Ver referencias. [12] [13] [14] [15]
Generalizaciones
- Una diferencia finita generalizada se define generalmente como donde μ = ( μ 0 ,…, μ N ) es su vector de coeficientes. Una diferencia infinita es una generalización adicional, donde la suma finita anterior se reemplaza por una serie infinita . Otra forma de generalización es hacer que los coeficientes μ k dependan del punto x : μ k = μ k ( x ) , considerando así la diferencia finita ponderada . También se puede hacer que el paso h dependa del punto x : h = h ( x ) . Estas generalizaciones son útiles para construir diferentes módulos de continuidad .
- La diferencia generalizada puede verse como los anillos polinomiales R [ T h ] . Conduce a álgebras de diferencias.
- El operador de diferencia se generaliza a la inversión de Möbius sobre un conjunto parcialmente ordenado .
- Como operador de convolución: a través del formalismo de álgebras de incidencia , los operadores de diferencia y otras inversiones de Möbius se pueden representar mediante convolución con una función en el poset, llamada función de Möbius μ ; para el operador de diferencia, μ es la secuencia (1, −1, 0, 0, 0,…) .
Diferencias finitas multivariadas
Las diferencias finitas se pueden considerar en más de una variable. Son análogos a las derivadas parciales en varias variables.
Algunas aproximaciones derivadas parciales son:
Alternativamente, para aplicaciones en las que el cálculo de f es el paso más costoso, y se deben calcular tanto la primera como la segunda derivada, una fórmula más eficiente para el último caso es
ya que los únicos valores para calcular que no son necesarios para las cuatro ecuaciones anteriores son f ( x + h , y + k ) y f ( x - h , y - k ) .
Ver también
- Esquema de diferenciación central
- Cálculo discreto
- Diferencias divididas
- Coeficientes de diferencias finitas
- Método de diferencias finitas
- Método de volumen finito
- Plantilla de cinco puntos
- Conjetura de Gilbreath
- Módulo de continuidad
- Integral de Nörlund – Rice
- Diferenciación numérica
- Secuencia de Sheffer
- Suma por partes
- Cálculo de escala de tiempo
- Cálculo umbral
- Esquema de diferenciación de ceñida para convección
- FDTD
- Esquema FTCS
Referencias
- ^ a b c Paul Wilmott; Sam Howison; Jeff Dewynne (1995). Las matemáticas de los derivados financieros: una introducción del estudiante . Prensa de la Universidad de Cambridge. pag. 137 . ISBN 978-0-521-49789-3.
- ^ a b c Peter Olver (2013). Introducción a las ecuaciones diferenciales parciales . Springer Science & Business Media. pag. 182. ISBN 978-3-319-02099-0.
- ^ a b c M Hanif Chaudhry (2007). Flujo de canal abierto . Saltador. pag. 369. ISBN 978-0-387-68648-6.
- ^ Jordán, op. cit., pág. 1 y Milne-Thomson, pág. xxi. Milne-Thomson, Louis Melville (2000): El cálculo de diferencias finitas (Chelsea Pub Co, 2000) ISBN 978-0821821077
- ^ Fraser, Duncan C. (1 de enero de 1909). "Sobre la delimitación gráfica de las fórmulas de interpolación" . Revista del Instituto de Actuarios . 43 (2): 235–241. doi : 10.1017 / S002026810002494X . Consultado el 17 de abril de 2017 .
- ^ Newton, Isaac, (1687). Principia , Libro III, Lema V, Caso 1
- ^ Richtmeyer, D. y Morton, KW, (1967). Métodos de diferencia para problemas de valor inicial , 2ª ed., Wiley, Nueva York.
- ↑ Boole, George , (1872). Tratado sobre el cálculo de diferencias finitas , 2ª ed., Macmillan and Company. En línea . Además, [edición de Dover 1960]
- ↑ Jordan, Charles, (1939/1965). "Cálculo de diferencias finitas", Chelsea Publishing. En línea: [1]
- ^ Zachos, C. (2008). "Deformaciones umbrales en espacio-tiempo discreto". International Journal of Modern Physics A . 23 (13): 2005-2014. arXiv : 0710.2306 . Código bibliográfico : 2008IJMPA..23.2005Z . doi : 10.1142 / S0217751X08040548 . S2CID 16797959 .
- ^ Curtright, TL; Zachos, CK (2013). "Umbral Vade Mecum". Fronteras en física . 1 : 15. arXiv : 1304.0429 . Código Bibliográfico : 2013FrP ..... 1 ... 15C . doi : 10.3389 / fphy.2013.00015 . S2CID 14106142 .
- ^ Levy, H .; Lessman, F. (1992). Ecuaciones en diferencias finitas . Dover. ISBN 0-486-67260-3.
- ^ Ames, WF, (1977). Métodos numéricos para ecuaciones diferenciales parciales , Sección 1.6. Academic Press, Nueva York. ISBN 0-12-056760-1 .
- ^ Hildebrand, FB , (1968). Ecuaciones y simulaciones en diferencias finitas , Sección 2.2, Prentice-Hall, Englewood Cliffs, Nueva Jersey.
- ^ Flajolet, Philippe; Sedgewick, Robert (1995). "Transformadas de Mellin y asintóticas: diferencias finitas e integrales de Rice" (PDF) . Informática Teórica . 144 (1-2): 101-124. doi : 10.1016 / 0304-3975 (94) 00281-M ..
- Richardson, CH (1954): Introducción al cálculo de diferencias finitas (Van Nostrand (1954) copia en línea
- Mickens, RE (1991): Ecuaciones en diferencias: teoría y aplicaciones (Chapman y Hall / CRC) ISBN 978-0442001360
enlaces externos
- "Cálculo de diferencias finitas" , Enciclopedia de Matemáticas , EMS Press , 2001 [1994]
- Tabla de fórmulas de diferencias finitas útiles generadas con Mathematica
- D. Gleich (2005), Cálculo finito: un tutorial para resolver sumas desagradables
- Segunda derivada discreta de puntos espaciados desigualmente