stringtranslate.com

Aproximación de Stirling

Comparación de la aproximación de Stirling con el factorial

En matemáticas , la aproximación de Stirling (o fórmula de Stirling ) es una aproximación asintótica para los factoriales . Es una buena aproximación, que conduce a resultados precisos incluso para valores pequeños de . Recibe su nombre de James Stirling , aunque Abraham de Moivre fue el primero en formular un resultado relacionado pero menos preciso . [1] [2] [3]

Una forma de expresar la aproximación implica el logaritmo del factorial: donde la notación O mayúscula significa que, para todos los valores suficientemente grandes de , la diferencia entre y será como máximo proporcional al logaritmo de . En aplicaciones informáticas como el límite inferior del peor caso para la clasificación por comparación , es conveniente utilizar en su lugar el logaritmo binario , dando la forma equivalente El término de error en cualquiera de las bases se puede expresar de forma más precisa como , correspondiente a una fórmula aproximada para el factorial en sí, Aquí el signo significa que las dos cantidades son asintóticas , es decir, que su relación tiende a 1 cuando tiende a infinito. La siguiente versión del límite se cumple para todos los n ≥ 1 {\displaystyle n\geq 1} , en lugar de solo de forma asintótica:

Derivación

En términos generales, la versión más simple de la fórmula de Stirling se puede obtener rápidamente aproximando la suma con una integral :

La fórmula completa, junto con estimaciones precisas de su error, se puede obtener de la siguiente manera. En lugar de aproximar , se considera su logaritmo natural , ya que se trata de una función que varía lentamente :

El lado derecho de esta ecuación menos es la aproximación por la regla del trapezoide de la integral

y el error en esta aproximación viene dado por la fórmula de Euler-Maclaurin :

donde es un número de Bernoulli y R m , n es el término restante en la fórmula de Euler-Maclaurin. Tome límites para encontrar que

Denotemos este límite como . Debido a que el resto R m , n en la fórmula de Euler-Maclaurin satisface

donde se utiliza la notación big-O , la combinación de las ecuaciones anteriores produce la fórmula de aproximación en su forma logarítmica:

Tomando la exponencial de ambos lados y eligiendo cualquier entero positivo , se obtiene una fórmula que involucra una cantidad desconocida . Para m = 1 , la fórmula es

La cantidad se puede hallar tomando el límite en ambos lados a medida que tiende a infinito y utilizando el producto de Wallis , que demuestra que . Por lo tanto, se obtiene la fórmula de Stirling:

Derivaciones alternativas

Una fórmula alternativa para utilizar la función gamma es (como se puede ver por la integración repetida por partes). Reescribiendo y cambiando las variables x = ny , se obtiene Aplicando el método de Laplace se tiene que recupera la fórmula de Stirling:

Órdenes superiores

De hecho, también se pueden obtener correcciones adicionales utilizando el método de Laplace. Del resultado anterior, sabemos que , por lo que "quitamos" este término dominante y luego realizamos dos cambios de variables para obtener: Para verificar esto: .

Ahora la función es unimodal, con valor máximo cero. Localmente alrededor de cero, parece , por lo que podemos realizar el método de Laplace. Para extender el método de Laplace a órdenes superiores, realizamos otro cambio de variables mediante . Esta ecuación no se puede resolver en forma cerrada, pero se puede resolver mediante expansión en serie, lo que nos da . Ahora, volvamos a la ecuación para obtener. Observe cómo no necesitamos encontrar realmente , ya que se cancela con la integral. Se pueden lograr órdenes superiores calculando más términos en , que se pueden obtener mediante programación. [nota 1]

De esta manera obtenemos la fórmula de Stirling de dos órdenes:

Versión analítica compleja

Una versión de análisis complejo de este método [4] es considerar como un coeficiente de Taylor de la función exponencial , calculado mediante la fórmula integral de Cauchy como

Esta integral de línea puede entonces aproximarse utilizando el método del punto de silla con una elección apropiada del radio del contorno . La porción dominante de la integral cerca del punto de silla se aproxima entonces mediante una integral real y el método de Laplace, mientras que la porción restante de la integral puede acotarse por encima para dar un término de error.

Velocidad de convergencia y estimaciones de error

El error relativo en una serie de Stirling truncada vs. , para 0 a 5 términos. Los puntos de quiebre en las curvas representan puntos donde la serie truncada coincide con Γ( n + 1) .

La fórmula de Stirling es de hecho la primera aproximación a la siguiente serie (ahora llamada serie de Stirling ): [5]

G. Nemes dio una fórmula explícita para los coeficientes de esta serie. [6] Otros términos se enumeran en la Enciclopedia en línea de secuencias de enteros como A001163 y A001164. El primer gráfico de esta sección muestra el error relativo vs. , para 1 a través de los 5 términos enumerados anteriormente. (Bender y Orszag [7] p. 218) da la fórmula asintótica para los coeficientes: que muestra que crece superexponencialmente y que, por la prueba de la razón, el radio de convergencia es cero.

El error relativo en una serie Stirling truncada en función del número de términos utilizados

Como n → ∞ , el error en la serie truncada es asintóticamente igual al primer término omitido. Este es un ejemplo de una expansión asintótica . No es una serie convergente ; para cualquier valor particular de solo hay una cierta cantidad de términos de la serie que mejoran la precisión, después de lo cual la precisión empeora. Esto se muestra en el siguiente gráfico, que muestra el error relativo en función del número de términos en la serie, para un mayor número de términos. Más precisamente, sea S ( n , t ) la serie de Stirling para términos evaluados en  . Los gráficos muestran que, cuando es pequeño, es esencialmente el error relativo.

Escribiendo la serie de Stirling en la forma se sabe que el error al truncar la serie es siempre de signo opuesto y como máximo de la misma magnitud que el primer término omitido.

Otros límites, debidos a Robbins, [8] válidos para todos los números enteros positivos son Este límite superior corresponde a detener la serie anterior para después del término. El límite inferior es más débil que el obtenido al detener la serie después del término. Una versión más flexible de este límite es que para todos los .

Fórmula de Stirling para la función gamma

Para todos los números enteros positivos, donde Γ denota la función gamma .

Sin embargo, la función gamma, a diferencia del factorial, está definida de manera más amplia para todos los números complejos que no sean enteros no positivos; no obstante, la fórmula de Stirling aún puede aplicarse. Si Re( z ) > 0 , entonces

La integración repetida por partes da

donde es el n ° de Bernoulli (nótese que el límite de la suma como no es convergente, por lo que esta fórmula es solo una expansión asintótica ). La fórmula es válida para valores suficientemente grandes en valor absoluto, cuando | arg( z ) | < π − ε , donde ε es positivo, con un término de error de O ( z −2 N + 1 ) . La aproximación correspondiente ahora puede escribirse:

donde la expansión es idéntica a la de la serie de Stirling anterior para , excepto que se reemplaza con z − 1 . [9]

Otra aplicación de esta expansión asintótica es para el argumento complejo z con Re( z ) constante . Véase, por ejemplo, la fórmula de Stirling aplicada en Im( z ) = t de la función theta de Riemann-Siegel en la línea recta 1/4 + eso .

Límites de error

Para cualquier entero positivo , se introduce la siguiente notación: y

Entonces [10] [11]

Para obtener más información y otros límites de error, consulte los artículos citados.

Una versión convergente de la fórmula de Stirling

Thomas Bayes demostró, en una carta a John Canton publicada por la Royal Society en 1763, que la fórmula de Stirling no daba una serie convergente . [12] Obtener una versión convergente de la fórmula de Stirling implica evaluar la fórmula de Binet :

Una forma de hacer esto es mediante una serie convergente de factoriales ascendentes invertidos . Si entonces donde donde s ( nk ) denota los números de Stirling de primera especie . De esto se obtiene una versión de la serie de Stirling que converge cuando Re( x ) > 0 . La fórmula de Stirling también puede darse en forma convergente como [13] donde

Versiones aptas para calculadoras

La aproximación y su forma equivalente se pueden obtener reordenando la fórmula extendida de Stirling y observando una coincidencia entre la serie de potencias resultante y la expansión de la serie de Taylor de la función seno hiperbólico . Esta aproximación es válida para más de 8 dígitos decimales para z con una parte real mayor que 8. Robert H. Windschitl la sugirió en 2002 para calcular la función gamma con una precisión razonable en calculadoras con memoria de programa o de registro limitada. [14]

Gergő Nemes propuso en 2007 una aproximación que da el mismo número de dígitos exactos que la aproximación Windschitl pero es mucho más simple: [15] o equivalentemente,

Una aproximación alternativa para la función gamma establecida por Srinivasa Ramanujan en el cuaderno perdido de Ramanujan [16] es para x ≥ 0 . La aproximación equivalente para ln n ! tiene un error asintótico de 1/1400 n.º 3 y viene dado por

La aproximación se puede hacer más precisa dando pares de límites superior e inferior; una de esas desigualdades es [17] [18] [19] [20]

Historia

La fórmula fue descubierta por primera vez por Abraham de Moivre [2] en la forma

De Moivre dio una expresión aproximada en números racionales para el logaritmo natural de la constante. La contribución de Stirling consistió en demostrar que la constante es precisamente . [3]

Véase también

Referencias

  1. ^ Dutka, Jacques (1991), "La historia temprana de la función factorial", Archivo de Historia de las Ciencias Exactas , 43 (3): 225–249, doi :10.1007/BF00389433, S2CID  122237769
  2. ^ ab Le Cam, L. (1986), "El teorema del límite central alrededor de 1935", Statistical Science , 1 (1): 78–96, doi : 10.1214/ss/1177013818 , JSTOR  2245503, MR  0833276; ver pág. 81, "El resultado, obtenido utilizando una fórmula originalmente probada por de Moivre pero ahora llamada fórmula de Stirling, aparece en su 'Doctrina de las probabilidades' de 1733".
  3. ^ ab Pearson, Karl (1924), "Nota histórica sobre el origen de la curva normal de errores", Biometrika , 16 (3/4): 402–404 [p. 403], doi :10.2307/2331714, JSTOR  2331714, Considero que el hecho de que Stirling demostrara que la constante aritmética de De Moivre era no le da derecho a afirmar el teorema, [...]
  4. ^ Flajolet, Philippe; Sedgewick, Robert (2009), Combinatoria analítica , Cambridge, Reino Unido: Cambridge University Press, pág. 555, doi :10.1017/CBO9780511801655, ISBN 978-0-521-89806-5, MR  2483235, S2CID  27509971
  5. ^ Olver, FWJ; Olde Daalhuis, AB; Lozier, DW; Schneider, BI; Boisvert, RF; Clark, CW; Miller, BR y Saunders, BV, "5.11 Propiedades de la función gamma: expansiones asintóticas", Biblioteca digital de funciones matemáticas del NIST , versión 1.0.13 del 16 de septiembre de 2016
  6. ^ Nemes, Gergő (2010), "Sobre los coeficientes de la expansión asintótica de ", Journal of Integer Sequences , 13 (6): 5
  7. ^ Bender, Carl M.; Orszag, Steven A. (2009). Métodos matemáticos avanzados para científicos e ingenieros. 1: Métodos asintóticos y teoría de perturbaciones (edición posterior). Nueva York, NY: Springer. ISBN 978-0-387-98931-0.
  8. ^ Robbins, Herbert (1955), "Una observación sobre la fórmula de Stirling", The American Mathematical Monthly , 62 (1): 26–29, doi :10.2307/2308012, JSTOR  2308012
  9. ^ Spiegel, MR (1999), Manual matemático de fórmulas y tablas , McGraw-Hill, pág. 148
  10. ^ Schäfke, FW; Sattler, A. (1990), "Restgliedabschätzungen für die Stirlingsche Reihe", Note di Matematica , 10 (supl. 2): 453–470, MR  1221957
  11. ^ G. Nemes, Límites de error y mejoras exponenciales para las expansiones asintóticas de la función gamma y su recíproca, Proc. Roy. Soc. Edinburgh Sect. A 145 (2015), 571–596.
  12. ^ Bayes, Thomas (24 de noviembre de 1763), "Una carta del difunto reverendo Sr. Thomas Bayes, FRS a John Canton, MA y FRS" (PDF) , Philosophical Transactions of the Royal Society of London, Serie I , 53 : 269, Bibcode :1763RSPT...53..269B, archivado (PDF) desde el original el 28 de enero de 2012 , consultado el 1 de marzo de 2012
  13. ^ Artin, Emil (2015). La función gamma . Dover. pág. 24.
  14. ^ Toth, VT Calculadoras programables: Calculadoras y la función gamma (2006) Archivado el 31 de diciembre de 2005 en Wayback Machine .
  15. ^ Nemes, Gergő (2010), "Nueva expansión asintótica para la función Gamma", Archiv der Mathematik , 95 (2): 161–169, doi :10.1007/s00013-010-0146-9, S2CID  121820640
  16. ^ Ramanujan, Srinivasa (14 de agosto de 1920), Cuaderno perdido y otros documentos inéditos, pág. 339 – vía Internet Archive
  17. ^ Karatsuba, Ekatherina A. (2001), "Sobre la representación asintótica de la función gamma de Euler por Ramanujan", Journal of Computational and Applied Mathematics , 135 (2): 225–240, Bibcode :2001JCoAM.135..225K, doi : 10.1016/S0377-0427(00)00586-0 , MR  1850542
  18. ^ Mortici, Cristinel (2011), "Estimación de Ramanujan para la función gamma mediante argumentos de monotonicidad", Ramanujan J. , 25 (2): 149–154, doi :10.1007/s11139-010-9265-y, S2CID  119530041
  19. ^ Mortici, Cristinel (2011), "Fórmulas asintóticas mejoradas para la función gamma", Comput. Math. Appl. , 61 (11): 3364–3369, doi :10.1016/j.camwa.2011.04.036.
  20. ^ Mortici, Cristinel (2011), "Sobre la fórmula del argumento grande de Ramanujan para la función gamma", Ramanujan J. , 26 (2): 185–192, doi :10.1007/s11139-010-9281-y, S2CID  120371952.

Lectura adicional

  1. ^ Por ejemplo, un programa en Mathematica:
    serie = tau - tau ^ 2 / 6 + tau ^ 3 / 36 + tau ^ 4 * a + tau ^ 5 * b ; (*elige el a,b correcto para que la serie sea igual a 0 en órdenes superiores*) Serie [ tau ^ 2 / 2 + 1 + t - Exp [ t ] /. t -> serie , { tau , 0 , 8 }]                       (*ahora haz la integral*) integral = Integrar [ Exp [ - x * tau ^ 2 / 2 ] * D [ serie /. a -> 0 /. b -> 0 , tau ], { tau , - Infinito , Infinito }]; Simplificar [ integral / Sqrt [ 2 * Pi ] * Sqrt [ x ]]                

Enlaces externos