stringtranslate.com

Diferencia finita

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 ba , 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 .

El operador de diferencia , comúnmente denominado como , es el operador que asigna una función f a la función definida por Una ecuación de diferencia es una ecuación funcional que involucra al operador de diferencia finita de la misma manera que una ecuación diferencial involucra derivadas . Existen muchas similitudes entre las ecuaciones de diferencia y las ecuaciones diferenciales, especialmente en los métodos de resolución. Ciertas relaciones de recurrencia se pueden escribir como ecuaciones de diferencia reemplazando la notación de iteración con diferencias finitas.

En el análisis numérico , las diferencias finitas se utilizan ampliamente para aproximar derivadas, y el término "diferencia finita" se utiliza a menudo como abreviatura de "aproximación de diferencias finitas de derivadas". [1] [2] [3] Las aproximaciones de diferencias finitas son cocientes de diferencias finitas en la terminología empleada anteriormente.

Las diferencias finitas fueron introducidas por Brook Taylor en 1715 y también han sido estudiadas como objetos matemáticos abstractos y autónomos en obras de George Boole (1860), LM Milne-Thomson (1933) y Károly Jordan  [de] (1939). Las diferencias finitas tienen su origen en uno de los algoritmos de Jost Bürgi ( c.  1592 ) y en el trabajo de otros, incluido Isaac Newton . El cálculo formal de diferencias finitas puede considerarse una alternativa al cálculo de infinitesimales . [4]

Tipos básicos

Los tres tipos de diferencias finitas. La diferencia central respecto de x proporciona la mejor aproximación de la derivada de la función en x .

Se consideran comúnmente tres tipos básicos: diferencias finitas hacia adelante , hacia atrás y centrales . [1] [2] [3]

Adiferencia hacia adelante , denotadapor unafunciónfes una función definida como

Dependiendo de la aplicación, el espaciamiento h puede ser variable o constante. Cuando se omite, h se toma como 1; es decir,

ALa diferencia hacia atrás utiliza los valores de la función enxy x h , en lugar de los valores en x + hx:

Finalmente, elLa diferencia central viene dada por

Relación con los derivados

La diferencia finita se utiliza a menudo como una aproximación de la derivada, normalmente 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 acercarse a cero, entonces el lado derecho de la ecuación anterior se escribiría

Por lo tanto, la diferencia hacia delante dividida por h aproxima la derivada cuando h es pequeña. El error en esta aproximación se puede derivar del teorema de Taylor . Suponiendo que f es dos veces diferenciable, tenemos

La misma fórmula se aplica a la diferencia hacia atrás:

Sin embargo, la diferencia central (también llamada centrada) produce una aproximación más precisa. Si f es tres veces diferenciable,

Sin embargo, el problema principal [ cita requerida ] con el método de diferencia central es que las funciones oscilantes pueden dar como resultado cero derivada. 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. Véase también Derivada simétrica .

Los autores para quienes las diferencias finitas significan aproximaciones de 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 de diferencias finitas para derivadas de orden superior y operadores diferenciales. Por ejemplo, utilizando la fórmula de diferencia central anterior para f  ′( x + yo/2 ) ​​y f  ′( xyo/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

De manera similar podemos aplicar otras fórmulas de diferenciación de manera recursiva.

Segundo orden hacia adelante
Segundo orden hacia atrás

De manera más general, las diferencias hacia adelante, hacia atrás y centrales de orden n se dan, respectivamente,

Adelante
Hacia atrás
Central

Estas ecuaciones utilizan coeficientes binomiales después del signo de suma que se muestra como (yo
)
. Cada fila deltriángulo de Pascalproporciona el coeficiente para cada valor dei.

Nótese que la diferencia central, para n impar , tendrá h multiplicada por números no enteros. Esto suele ser un problema porque equivale a cambiar el intervalo de discretización. El problema se puede solucionar sustituyendo el promedio de y

Las diferencias hacia delante aplicadas a una secuencia se denominan a veces transformada binomial de la secuencia y tienen varias propiedades combinatorias interesantes. Las diferencias hacia delante se pueden evaluar utilizando la integral de Nörlund–Rice . La representación integral para este tipo de series es interesante, porque la integral a menudo se puede evaluar utilizando técnicas de expansión asintótica o de punto de silla ; por el contrario, la serie de diferencias hacia delante puede ser extremadamente difícil de evaluar numéricamente, porque los coeficientes binomiales crecen rápidamente para valores grandes de n .

La relación de estas diferencias de orden superior con las respectivas derivadas es sencilla:

También se pueden utilizar diferencias de orden superior para construir mejores aproximaciones. Como se mencionó anteriormente, la diferencia de primer orden aproxima 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 demostrar desarrollando la expresión anterior en series de Taylor o utilizando el cálculo de diferencias finitas, que se explica a continuación.

Si es necesario, la diferencia finita se puede centrar en cualquier punto mezclando diferencias hacia adelante, hacia atrás y centrales.

Polinomios

Para un polinomio dado de grado n ≥ 1 , expresado en la función P ( x ) , con números reales a ≠ 0 y b y términos de orden inferior (si los hay) marcados como lot :

Después de n diferencias por pares, se puede obtener el siguiente resultado, donde h ≠ 0 es un número real que marca la diferencia aritmética: [5]

Solo queda el coeficiente del término de orden más alto. Como este resultado es constante con respecto a x , cualquier diferencia por pares adicional tendrá el valor 0 .

Prueba inductiva

Caso base

Sea Q ( x ) un polinomio de grado 1 :

Esto lo demuestra para el caso base.

Paso inductivo

Sea R ( x ) un polinomio de grado m  − 1 donde m ≥ 2 y el coeficiente del término de mayor orden sea a ≠ 0 . Suponiendo que lo siguiente es cierto para todos los polinomios de grado m  − 1 :

Sea S ( x ) un polinomio de grado m . Con una diferencia por pares:

Como ahm ≠ 0 , esto da como resultado un polinomio T ( x ) de grado m − 1 , con ahm como el coeficiente del término de orden más alto. Dado el supuesto anterior y m − 1 diferencias por pares (lo que da como resultado un total de m diferencias por pares para S ( x ) ), se puede encontrar que:

Con esto finaliza la prueba.

Solicitud

Esta identidad se puede utilizar para encontrar el polinomio de grado más bajo que intercepta una cantidad de puntos ( x , y ) donde la diferencia en el eje x de un punto al siguiente es una constante h ≠ 0 . Por ejemplo, dados los siguientes puntos:

Podemos utilizar una tabla de diferencias, donde para todas las celdas a la derecha de la primera y , existe la siguiente relación con las celdas en la columna inmediatamente a la izquierda para una celda ( a + 1, b + 1) , con la celda superior izquierda en la coordenada (0, 0) :

Para encontrar el primer término se puede utilizar la siguiente tabla:

Esto da como resultado una constante de 648. La diferencia aritmética es h = 3 , como se estableció anteriormente. Dado el número de diferencias por pares necesarias para alcanzar la constante, se puede suponer que se trata de un polinomio de grado 3. Por lo tanto, utilizando la identidad anterior:

Resolviendo para a , se puede encontrar que tiene el valor 4 . Por lo tanto, el primer término del polinomio es 4 x 3 .

Luego, restando el primer término, que baja el grado del polinomio, y encontrando nuevamente la diferencia finita:

Aquí, la constante se alcanza después de solo dos diferencias por pares, por lo que se obtiene el siguiente resultado:

Resolviendo para a , que es −17 , el segundo término del polinomio es −17 x 2 .

Pasando al siguiente término, restando el segundo término:

Por lo tanto, la constante se alcanza después de una sola diferencia por pares:

Se puede encontrar que a = 36 y por lo tanto el tercer término del polinomio es 36 x . Restando el tercer término:

Sin diferencias por pares, se encuentra que el cuarto y último término del polinomio es la constante −19 . Por lo tanto, se encuentra el polinomio de menor grado que intercepta todos los puntos de la primera tabla:

Granos de tamaño arbitrario

Utilizando el álgebra lineal se pueden construir aproximaciones de 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 de derivada. Esto implica resolver un sistema lineal de modo 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. Tales fórmulas se pueden representar gráficamente en una cuadrícula hexagonal o en forma de diamante. [6] 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, uno debe muestrear cada vez menos puntos en un lado. [7] Se pueden construir aproximaciones de diferencias finitas para plantillas no estándar (e incluso no enteras) dada una plantilla arbitraria y un orden de derivada deseado. [8]

Propiedades

En ecuaciones diferenciales

Una aplicación importante de las diferencias finitas es en el análisis numérico , especialmente en las ecuaciones diferenciales numéricas , que tienen como objetivo 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 disciplinas de ingeniería, como la ingeniería térmica , la mecánica de fluidos , etc.

Serie de Newton

La serie de Newton consiste en los términos de la ecuación diferencial directa de Newton , llamada así por Isaac Newton ; en esencia, es la fórmula de interpolación de Gregory-Newton [9] (llamada así por Isaac Newton y James Gregory ), publicada por primera vez en sus Principia Mathematica en 1687, [10] [11] es decir, el análogo discreto de la expansión continua de Taylor,

lo cual es válido para cualquier función polinómica f y para muchas (pero no todas) funciones analíticas . (No es válido 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, claramente, 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 siguiente.

Obsérvese la correspondencia formal de este resultado con el teorema de Taylor . Históricamente, este , así como la identidad de Chu-Vandermonde (que se deriva de él y corresponde al teorema del binomio ), se incluyen en las observaciones que maduraron hasta convertirse en el sistema de cálculo umbral .

Las expansiones de la serie de Newton pueden ser superiores a las expansiones de la serie de Taylor cuando se aplican a cantidades discretas como espines cuánticos (ver transformación de Holstein-Primakoff ), funciones de operador bosónico o estadísticas de conteo discretas. [12]

Para ilustrar cómo se puede utilizar la fórmula de Newton en la práctica real, considere los primeros términos de la duplicación de la secuencia de Fibonacci f = 2, 2, 4, ... Se puede encontrar un polinomio que reproduzca estos valores, calculando primero una tabla de diferencias y luego sustituyendo las diferencias que corresponden a x 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 , [13]

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 polinómica puede debilitarse hasta llegar a la suposición de que f es simplemente continua.

El teorema de Carlson establece las condiciones necesarias y suficientes para que una serie de Newton sea única, si es que existe. Sin embargo, en general, una serie de Newton no existe.

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 progresivas adecuadamente escaladas.

En una forma comprimida y ligeramente más general y con nodos equidistantes la fórmula se lee

Cálculo de diferencias finitas

La diferencia hacia adelante puede considerarse como un operador , llamado operador de diferencia , que asigna la función f a Δ h [  f  ] . [14] [15] Este operador equivale a donde  T h es el operador de desplazamiento con paso h , definido por T h [  f  ]( x ) = f  ( x + h )  , e  I  es el operador identidad .

La diferencia finita de órdenes superiores se puede definir de manera recursiva como Δen
≡ Δ hn - 1
hora
)
 .
Otra definición equivalente es Δen
≡ [ 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 :

Se aplican reglas de Leibniz similares para las diferencias centrales y retrógradas.

Aplicando formalmente la serie de Taylor con respecto a h , se obtiene la ecuación del operador donde D denota el operador de derivada continua convencional, que asigna 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 ; en el caso especial de que la serie de derivadas termine (cuando la función sobre la que se opera es un polinomio finito ) la expresión es exacta, para todos los tamaños de paso finitos, h . Por lo tanto, T h = e h D  , e invirtiendo formalmente el exponencial se obtiene Esta fórmula se cumple en el sentido de que ambos operadores dan el mismo resultado cuando se aplican a un polinomio.

Incluso en el caso de 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, si se conservan los dos primeros términos de la serie se obtiene 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 umbrales con sus análogos continuos ( límites h → 0 ).

Un gran número de relaciones diferenciales formales del cálculo estándar que involucran funciones   f  ( x ) se   asignan sistemáticamente a análogos de diferencias finitas umbrales que involucran f  ( x T-1
hora
)
 .

Por ejemplo, el análogo umbral de un monomio x n es una generalización del factorial descendente anterior ( símbolo k de Pochhammer ), de modo que 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

Al igual que en el límite continuo , la función propia de   Δh/yo   también resulta ser exponencial,

y por lo tanto, las sumas de Fourier de funciones continuas se asignan fácil y fielmente a sumas de Fourier umbrales , es decir, involucran los mismos coeficientes de Fourier que multiplican estos exponenciales de base umbrales. [16] Por lo tanto, este exponencial umbral 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 seno cardinal , y así sucesivamente. [17] Las ecuaciones de diferencias a menudo se pueden resolver con técnicas muy similares a las utilizadas para resolver ecuaciones diferenciales .

El operador inverso del operador de diferencia hacia adelante, y por lo tanto 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:

Todas las reglas anteriores se aplican igualmente a cualquier operador de diferencia como Δ , incluidos δ y  .

Ver referencias. [18] [19] [20] [21]

Generalizaciones

Diferencias finitas multivariadas

Las diferencias finitas pueden considerarse en más de una variable. Son análogas a las derivadas parciales en varias variables.

Algunas aproximaciones de derivadas parciales son:

Como alternativa, para aplicaciones en las que el cálculo de f es el paso más costoso y deben calcularse tanto la primera como la segunda derivada, una fórmula más eficiente para el último caso es ya que los únicos valores a calcular que no son ya necesarios para las cuatro ecuaciones anteriores son f  ( x + h , y + k ) y f  ( xh , yk ) .

Véase también

Referencias

  1. ^ abc Paul Wilmott; Sam Howison; Jeff Dewynne (1995). Las matemáticas de los derivados financieros: una introducción para estudiantes . Cambridge University Press. pág. 137. ISBN 978-0-521-49789-3.
  2. ^ abc Peter Olver (2013). Introducción a las ecuaciones diferenciales parciales . Springer Science & Business Media. pág. 182. ISBN 978-3-319-02099-0.
  3. ^ abc M Hanif Chaudhry (2007). Flujo en canal abierto . Springer. pág. 369. ISBN 978-0-387-68648-6.
  4. ^ Jordán, op. cit., p. 1 y Milne-Thomson, p. xxi. Milne-Thomson, Louis Melville (2000): El cálculo de las diferencias finitas (Chelsea Pub Co, 2000) ISBN 978-0821821077 
  5. ^ "Diferencias finitas de polinomios". 13 de febrero de 2018.
  6. ^ Fraser, Duncan C. (1 de enero de 1909). "Sobre la delineación gráfica de fórmulas de interpolación". Journal of the Institute of Actuaries . 43 (2): 235–241. doi :10.1017/S002026810002494X . Consultado el 17 de abril de 2017 .
  7. ^ notas
  8. ^ Calculadora de coeficientes de diferencias finitas
  9. ^ Burkard Polster /Mathologer (2021). "¿Por qué no enseñan el cálculo de Newton sobre 'qué viene después'?" en YouTube
  10. ^ Newton, Isaac, (1687). Principia, Libro III, Lema V, Caso 1
  11. ^ Iaroslav V. Blagouchine (2018). "Tres notas sobre las representaciones de Ser y Hasse para las funciones zeta" (PDF) . Integers (Electronic Journal of Combinatorial Number Theory) . 18A : 1–45. arXiv : 1606.02044 .
  12. ^ König, Jürgen; Hucht, Fred (2021). "Expansión de funciones de operador bosónico en serie de Newton". SciPost Physics . 10 (1): 007. arXiv : 2008.11139 . Bibcode :2021ScPP...10....7K. doi : 10.21468/SciPostPhys.10.1.007 . S2CID  221293056.
  13. ^ Richtmeyer, D. y Morton, KW, (1967). Métodos diferenciales para problemas de valor inicial , 2.ª ed., Wiley, Nueva York.
  14. ^ Boole, George (1872). Tratado sobre el cálculo de diferencias finitas (2.ª ed.). Macmillan and Company – vía Internet Archive .También, una edición reimpresa de Dover, 1960.
  15. ^ Jordan, Charles (1965) [1939]. Cálculo de diferencias finitas. Chelsea Publishing. ISBN 978-0-8284-0033-6– a través de Google Books.
  16. ^ Zachos, C. (2008). "Deformaciones umbrales en el espacio-tiempo discreto". Revista Internacional de Física Moderna A . 23 (13): 200–214. arXiv : 0710.2306 . Código Bibliográfico :2008IJMPA..23.2005Z. doi :10.1142/S0217751X08040548. S2CID  16797959.
  17. ^ Curtright, TL; Zachos, CK (2013). "Vademécum umbral". Frontiers in Physics . 1 : 15. arXiv : 1304.0429 . Bibcode :2013FrP.....1...15C. doi : 10.3389/fphy.2013.00015 . S2CID  14106142.
  18. ^ Levy, H.; Lessman, F. (1992). Ecuaciones en diferencias finitas . Dover. ISBN 0-486-67260-3.
  19. ^ Ames, WF (1977). Métodos numéricos para ecuaciones diferenciales parciales . Nueva York, NY: Academic Press. Sección 1.6. ISBN 0-12-056760-1.
  20. ^ Hildebrand, FB (1968). Ecuaciones y simulaciones de diferencias finitas . Englewood Cliffs, NJ: Prentice-Hall. Sección 2.2.
  21. ^ Flajolet, Philippe; Sedgewick, Robert (1995). "Transformadas de Mellin y asintóticas: diferencias finitas e integrales de Rice" (PDF) . Ciencias de la Computación Teórica . 144 (1–2): 101–124. doi :10.1016/0304-3975(94)00281-M.

Enlaces externos