Notación polaca inversa
Notación polaca inversa también conocida como notación Łukasiewicz inversa, notación posfija polaca o simplemente La notación posfija es una notación... (leer más)
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 diferencia. 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 la diferencia, comúnmente denotado Δ Δ {displaystyle Delta } es el operador que mapea una función f a la función Δ Δ [f]{displaystyle Delta [f]} definidas por
Una ecuación en diferencias es una ecuación funcional que involucra el operador de diferencias finitas de la misma manera que una ecuación diferencial involucra derivadas. Hay muchas similitudes entre las ecuaciones en diferencias y las ecuaciones diferenciales, especialmente en los métodos de resolución. Ciertas relaciones de recurrencia se pueden escribir como ecuaciones en diferencias reemplazando la notación de iteración con diferencias finitas.
En el análisis numérico, las diferencias finitas se usan mucho para aproximar derivadas, y el término "diferencia finita" se utiliza a menudo como una abreviatura de "aproximación de derivadas en diferencias finitas". Las aproximaciones de diferencias finitas son cocientes de diferencias finitas en la terminología empleada anteriormente.
Brook Taylor introdujo las diferencias finitas en 1715 y también se han estudiado como objetos matemáticos autónomos abstractos en obras de George Boole (1860), L. M. Milne-Thomson (1933) y Károly Jordan 1592) y el trabajo de otros, incluido Isaac Newton. El cálculo formal de diferencias finitas puede verse como una alternativa al cálculo de infinitesimales.
(1939). Las diferencias finitas tienen su origen en uno de los algoritmos de Jost Bürgi (c.Comúnmente se consideran tres tipos básicos: diferencias finitas hacia delante, hacia atrás y central.
A diferencia de futuro, denotado Δ Δ h[f],{displaystyle Delta _{h}[f],} de una función f es una función definida
Dependiendo de la aplicación, el espaciado h puede ser variable o constante. Cuando se omite, h se toma como 1; es decir,
Una diferencia hacia atrás usa los valores de función en x y x − h, en lugar de los valores en x + h y x:
Finalmente, la diferencia central viene dada por
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 acercarse a cero, entonces el lado derecho de la la ecuación anterior se escribiría
Por lo tanto, la diferencia directa dividida por h se aproxima a la derivada cuando h es pequeño. 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 vale 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 tres veces diferenciable,
Sin embargo, el problema principal con el método de la diferencia central es que las funciones oscilantes pueden producir derivada cero. Si f (nh) = 1 para n impar, y f (nh) = 2 para n incluso, 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 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).
De manera análoga, se pueden obtener aproximaciones en diferencias finitas a derivadas de orden superior y operadores diferenciales. Por ejemplo, usando 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 en diferencia central de la segunda derivada de f:
Del mismo modo, podemos aplicar otras fórmulas de diferenciación de forma recursiva.
De forma más general, las diferencias de nésimo orden adelante, atrás y central están dadas por, respectivamente,
o para h = 1,
Estas ecuaciones usan coeficientes binomiales después del signo de suma que se muestra como (n
i). Cada fila del triángulo de Pascal proporciona el coeficiente para cada valor de i.
Tenga en cuenta que la diferencia central, para impar n, tendrá h multiplicado por números no enteros. Esto suele ser un problema porque equivale a cambiar el intervalo de discretización. El problema se puede remediar tomando el promedio de δn[ f ](x − h/2) y δn[ f ](x + h/2).
Las diferencias directas aplicadas a una secuencia a veces se denominan transformadas binomiales de la secuencia y tienen varias propiedades combinatorias interesantes. Las diferencias directas pueden evaluarse mediante 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 diferencia directa puede ser extremadamente difícil de evaluar numéricamente, porque los coeficientes binomiales crecen rápidamente para grandes n.
La relación de estas diferencias de orden superior con las respectivas derivadas es sencilla,
Las diferencias de orden superior también se pueden usar 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 plazo de pedido h2. 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 en cualquier punto mezclando diferencias hacia adelante, hacia atrás y centrales.
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 l.o.t.:
P()x)=axn+bxn− − 1+l.o.t.{displaystyle P(x)=ax^{n}+bx^{n-1}+l.o.t.}
Después de n diferencias por pares, se puede lograr el siguiente resultado, donde h ≠ 0 es un número real que marca la diferencia aritmética:
Δ Δ hn[P]()x)=ahnn!{displaystyle Delta _{n}[P](x)=ah^{n}n!
Solo queda el coeficiente del término de mayor orden. Como este resultado es constante con respecto a x, cualquier otra diferencia por pares tendrá el valor 0.
Sea Q(x) un polinomio de grado 1:
Δ Δ h[Q]()x)=Q()x+h)− − Q()x)=[a()x+h)+b]− − [ax+b]=ah=ah11!{displaystyle Delta _{h}[Q](x)=Q(x+h)-Q(x)=[a(x+h)+b] - [ax+b]=ah=ah^{1}1}
Esto lo prueba para el caso base.
Sea R(x) un polinomio de grado m-1 donde m ≥ 2 y el coeficiente del término de mayor orden es a ≠ 0. Suponiendo que lo siguiente sea cierto para todos los polinomios de grado m-1:
Δ Δ hm− − 1[R]()x)=ahm− − 1()m− − 1)!{displaystyle Delta _{h}{m-1}[R](x)=ah^{m-1}(m-1)}
Sea S(x) un polinomio de grado m. Con una diferencia por pares:
Δ Δ h[S]()x)=[a()x+h)m+b()x+h)m− − 1+l.o.t.]− − [axm+bxm− − 1+l.o.t.]=ahmxm− − 1+l.o.t.=T()x){displaystyle Delta _{h}[S](x)=[a(x+h)^{m}+b(x+h)^{m-1}+l.o.t.] - [ax^{m}+bx^{m-1}+l.o.t.]=ahmx^{m-1}+l.o.t=T(x)}
Como ahm ≠ 0, esto da como resultado un polinomio T(x) de grado m-1, con ahm como el coeficiente del término de mayor orden. Dada la suposición 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:
Δ Δ hm− − 1[T]()x)=ahm⋅ ⋅ hm− − 1()m− − 1)!=ahmm!{displaystyle Delta _{h}{m-1}[T](x)=ahmcdot h^{m-1}(m-1)=ah^{m}m!}
Esto completa la prueba.
Esta identidad se puede usar para encontrar el polinomio de menor grado que intercepta un número de puntos (x, y) donde la diferencia en el eje x de un punto a la siguiente es una constante h ≠ 0. Por ejemplo, dados los siguientes puntos:
x | Sí. |
---|---|
1 | 4 |
4 | 109 |
7 | 772 |
10 | 2641 |
13 | 6364 |
Podemos usar una tabla de diferencias, donde todas las celdas a la derecha de la primera y, la siguiente relación con las celdas en la columna inmediatamente a la izquierda existe para una celda (a+1, b+1), con la celda superior izquierda en la coordenada (0, 0):
()a+1,b+1)=()a,b)− − ()a,b+1){displaystyle (a+1,b+1)=(a,b)-(a,b+1)}
Para encontrar el primer término, se puede usar la siguiente tabla:
x | Sí. | Δy | Δ2Sí. | Δ3Sí. |
---|---|---|---|---|
1 | 4 | |||
4 | 109 | 105 | ||
7 | 772 | 663 | 558 | |
10 | 2641 | 1869 | 1206 | 648 |
13 | 6364 | 3723 | 1854 | 648 |
Esto llega a una constante 648. La diferencia aritmética es h=3, como se estableció anteriormente. Dada la cantidad de diferencias por pares necesarias para alcanzar la constante, se puede suponer que se trata de un polinomio de grado 3. Por lo tanto, usando la identidad anterior:
648=a⋅ ⋅ 33⋅ ⋅ 3!=a⋅ ⋅ 27⋅ ⋅ 6=a⋅ ⋅ 162{displaystyle 648=acdot 3}cdot 3!=acdot 27cdot 6=acdot 162}
Resolviendo para a, se puede encontrar que tiene el valor 4. Por lo tanto, el primer término del polinomio es 4x3.
Luego, restando el primer término, que reduce el grado del polinomio, y volviendo a encontrar la diferencia finita:
x | Sí. | Δy | Δ2Sí. |
---|---|---|---|
1 | 4 a 4 1)3 = 4 - 4 = 0 | ||
4 | 109 - 4(4)3 = 109 - 256 = -147 | -147 | |
7 | 772 - 4(7)3 = 772 - 1372 = -600 | -453 | -306 |
10 | 2641 - 4(10)3 = 2641 - 4000 = -1359 | -759 | -306 |
13 | 6364 - 4(13)3 = 6364 - 8788 = -2424 | -1065 | -306 |
Aquí, la constante se logra después de solo 2 diferencias por pares, por lo tanto, el siguiente resultado:
− − 306=a⋅ ⋅ 32⋅ ⋅ 2!=a⋅ ⋅ 18{displaystyle -306=acdot 3^{2}cdot 2!=acdot 18}
Resolviendo a, que es -17, el segundo término del polinomio es -17x2.
Pasando al siguiente término, restando el segundo término:
x | Sí. | Δy |
---|---|---|
1 | 0 - (-17 1)2) = 0 + 17 = 17 | |
4 | -147 - (-17(4)2) = -147 + 272 = 125 | 108 |
7 | -600 - (-17(7)2) = -600 + 833 = 233 | 108 |
10 | -1359 - (-17(10)2) = -1359 + 1700 = 341 | 108 |
13 | -2424 - (-17(13)2) = -2424 + 2873 = 449 | 108 |
Por lo tanto, la constante se logra después de solo 1 diferencia por pares:
108=a⋅ ⋅ 31⋅ ⋅ 1!=a⋅ ⋅ 3{displaystyle 108=acdot 3^{1}cdot 1!=acdot 3}
Se puede encontrar que a = 36 y por lo tanto el tercer término del polinomio es 36x. Restando el tercer término:
x | Sí. |
---|---|
1 | 17 - 36(1) = 17 - 36 = -19 |
4 | 125 - 36(4) = 125 - 144 = -19 |
7 | 233 - 36(7) = 233 - 252 = -19 |
10 | 341 - 36(10) = 341 - 360 = -19 |
13 | 449 – 36(13) = 449 – 468 = –19 |
Sin ninguna diferencia 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 grado más bajo que intercepta todos los puntos en la primera tabla:
4x3− − 17x2+36x− − 19{displaystyle 4x^{3}-17x^{2}+36x-19}
Usando á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 derivado. 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.
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 de diferencia finita para plantillas no estándar (e incluso no enteras) dada una plantilla arbitraria y un orden derivado deseado.
Una aplicación importante de las diferencias finitas es el análisis numérico, especialmente en las ecuaciones diferenciales numéricas, cuyo objetivo es la solución numérica de las 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 disciplinas de ingeniería y ciencia computacional, como ingeniería térmica, mecánica de fluidos, etc.
La serie de Newton consta de los términos de la ecuación de diferencia directa de Newton, llamada así por Isaac Newton; en esencia, es la fórmula de interpolación de Newton, publicada por primera vez en su Principia Mathematica en 1687, es decir, el análogo discreto de la expansión continua de Taylor,
f()x)=.. k=0JUEGO JUEGO Δ Δ k[f]()a)k!()x− − a)k=.. k=0JUEGO JUEGO ()x− − ak)Δ Δ k[f]()a),{displaystyle f(x)=sum _{k=0}{infty }{frac [Delta ^{k](f] {k!},(x-a)_{k}=sum ¿Qué? } {binom {x-a}{k}, Delta.
que tiene para cualquier función polinomio f y para muchas (pero no todas) funciones analíticas. (No se mantiene cuando f es el tipo exponencial π π {displaystyle pi}. Esto se ve fácilmente, ya que la función sine desaparece en múltiples enteros π π {displaystyle pi}; la correspondiente serie Newton es idénticamente cero, ya que todas las diferencias finitas son cero en este caso. Sin embargo, claramente, la función sine 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, se asumen pasos unitarios para los cambios en los valores de x, h = 1 de la generalización a continuación.
Observe la correspondencia formal de este resultado con el teorema de Taylor. Históricamente, esto, así como la identidad Chu-Vandermonde,
(siguiendo de ello, y correspondientes al teorema del binomio), se incluyen en las observaciones que maduraron al 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 (consulte la transformación de Holstein-Primakoff), funciones de operadores bosónicos o estadísticas de conteo discretas.
Para ilustrar cómo se puede usar 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, primero calculando 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,
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 asumir que f es simplemente continua.
El teorema de Carlson proporciona las condiciones necesarias y suficientes para que una serie de Newton sea única, si existe. Sin embargo, una serie de Newton, en general, 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 un poco más general y con nodos equidistantes, la fórmula dice
La diferencia directa se puede considerar como un operador, llamado operador de diferencia, que asigna la función f a Δh[ f ]. Este operador equivale a
donde Th es el operador de desplazamiento con paso h, definido por Th[ f ](x) = f (x + h), y estilo I es el operador de identidad.
La diferencia finita de órdenes superiores se puede definir de forma recursiva como Δ n
h ≡ Δh(Δn − 1
h). Otra definición equivalente es Δn
h = [Th − I]n.
El operador diferencia Δh es un operador lineal, como tal satisface Δh[αf + βg](x) = α Δh[ f ](x) + β Δh[g](x).
También cumple una regla especial de Leibniz indicada anteriormente, Δh(f (x)g(x)) = (Δhf (x)) g(x+h) + f (x) (Δ hg(x)). Declaraciones similares valen para las diferencias hacia atrás y central.
Aplicando formalmente la serie de Taylor con respecto a h, se obtiene la fórmula
donde D denota el operador derivado 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ñas. Por lo tanto, Th = ehD, e invirtiendo formalmente los rendimientos exponenciales
Esta fórmula se cumple 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 de la derivada. Por ejemplo, conservar los dos primeros términos de la serie produce la aproximación de segundo orden a f ′(x) mencionado al final de la sección Diferencias de orden superior.
Las fórmulas análogas para los operadores de diferencia hacia atrás y central 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 (h → 0 límites),
[Δ Δ hh,xTh− − 1]=[D,x]=I.{displaystyle left[{frac {Delta [D,x]=I.
Un gran número de relaciones diferenciales formales de cálculo estándar que implican
funciones f (x) por lo tanto se asignan sistemáticamente a análogos de diferencia finita umbral que involucran f (xT−1
h).
Por ejemplo, el análogo umbral de un monomio xn es una generalización de lo anterior factorial descendente (símbolo k de Pochhammer),
para 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
Como en el límite continuo, la función propia de Δh /h también resulta ser exponencial,
y, por lo tanto, las sumas de Fourier de las funciones continuas se asignan fácilmente a las sumas de Fourier umbral fielmente, es decir, involucran los mismos coeficientes de Fourier que multiplican estos exponenciales de base umbral. Este exponencial umbral equivale así 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,
y así sucesivamente. Las ecuaciones en diferencias a menudo se pueden resolver con técnicas muy similares a las que se usan para resolver ecuaciones diferenciales.
El operador inverso del operador de diferencia directa, entonces la integral umbral, es el operador de suma indefinida o antidiferencia.
Al igual que las reglas para encontrar la derivada, tenemos:
Todas las reglas anteriores se aplican igualmente bien a cualquier operador de diferencia, incluido ∇ como a Δ.
Ver referencias.
Se pueden considerar diferencias finitas en más de una variable. Son análogas a las derivadas parciales en varias variables.
Algunas aproximaciones de derivadas parciales son:
Alternativamente, para aplicaciones en las que el cálculo de f es el paso más costoso, y tanto la primera como la segunda derivada deben ser calculado, una fórmula más eficiente para el último caso es
dado 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).
Notación polaca inversa también conocida como notación Łukasiewicz inversa, notación posfija polaca o simplemente La notación posfija es una notación... (leer más)
Exponenciación es una operación matemática, escrita como bn< /span>, que involucra dos números, la base b y el exponente o poder n, y se pronuncia como... (leer más)
En matemáticas, se dice que dos conjuntos son disjuntos si no tienen ningún elemento en común. De manera equivalente, dos conjuntos disjuntos son conjuntos... (leer más)