stringtranslate.com

función gamma

En matemáticas , la función gamma (representada por Γ , la letra mayúscula gamma del alfabeto griego ) es una extensión comúnmente utilizada de la función factorial para números complejos . La función gamma está definida para todos los números complejos excepto los enteros no positivos. Para cada entero positivo n ,

Derivada por Daniel Bernoulli , para números complejos con parte real positiva, la función gamma se define mediante una integral impropia convergente :

La función gamma entonces se define como la continuación analítica de esta función integral a una función meromorfa que es holomorfa en todo el plano complejo excepto cero y los enteros negativos, donde la función tiene polos simples . [ se necesita aclaración ]

La función gamma no tiene ceros, por lo que la función gamma recíproca 1/Γ( z ) es una función completa . De hecho, la función gamma corresponde a la transformada de Mellin de la función exponencial negativa:

Existen otras extensiones de la función factorial, pero la función gamma es la más popular y útil. Es un componente de varias funciones de distribución de probabilidad y, como tal, es aplicable en los campos de la probabilidad y la estadística , así como en la combinatoria .

Motivación

interpola la función factorial a valores no enteros.

La función gamma puede verse como una solución al problema de interpolación de encontrar una curva suave que conecte los puntos de la secuencia factorial: para todos los valores enteros positivos de . La fórmula simple para el factorial, x ! = 1 × 2 × ⋯ × x solo es válido cuando x es un entero positivo, y ninguna función elemental tiene esta propiedad, pero una buena solución es la función gamma . [1]

La función gamma no sólo es fluida sino también analítica (excepto en los números enteros no positivos) y puede definirse de varias maneras explícitas. Sin embargo, no es la única función analítica que extiende el factorial, ya que se puede agregar cualquier función analítica que sea cero en los números enteros positivos, como para un número entero . [1] Esta función se conoce como función pseudogamma , siendo la más famosa la función de Hadamard . [2]

La función gamma, Γ( z ) en azul, representada junto con Γ( z ) + sin(π z ) en verde. Observe la intersección en números enteros positivos. Ambas son extensiones válidas de los factoriales a una función meromórfica en el plano complejo.

Un requisito más restrictivo es la ecuación funcional que interpola el factorial desplazado  : [3] [4]

Pero esto todavía no da una solución única, ya que permite la multiplicación por cualquier función periódica con y , como por ejemplo . Una forma de resolver la ambigüedad es el teorema de Bohr-Mollerup , que muestra que la única función de interpolación sobre los reales positivos es logarítmicamente convexa (superconvexa), [5] lo que significa que es convexa (donde está el logaritmo natural ). [6]

Definición

Definición principal

La notación se debe a Legendre . [1] Si la parte real del número complejo  z es estrictamente positiva ( ), entonces la integral converge absolutamente , y se conoce como integral de Euler de segunda especie . (La integral de Euler de primer tipo es la función beta . [1] ) Usando la integración por partes , se ve que:

Valor absoluto (vertical) y argumento (color) de la función gamma en el plano complejo
Valor absoluto (vertical) y argumento (color) de la función gamma en el plano complejo

Reconociendo que como

Podemos calcular :

Así podemos demostrar que para cualquier entero positivo n por inducción . Específicamente, el caso base es ese y el paso de inducción es que

La identidad se puede usar (o, dando el mismo resultado, se puede usar la continuación analítica ) para extender de manera única la formulación integral de a una función meromórfica definida para todos los números complejos z , excepto los enteros menores o iguales a cero. [1] Es esta versión extendida la que comúnmente se conoce como función gamma. [1]

Definiciones alternativas

Hay muchas definiciones equivalentes.

La definición de Euler como producto infinito.

Para un número entero fijo , a medida que el número entero aumenta, tenemos que [7]

Si no es un número entero, entonces no es posible decir si esta ecuación es verdadera porque todavía no hemos definido (en esta sección) la función factorial para números no enteros. Sin embargo, obtenemos una extensión única de la función factorial a los números no enteros al insistir en que esta ecuación continúe siendo válida cuando el número entero arbitrario se reemplaza por un número complejo arbitrario .

Multiplicar ambos lados por da este producto infinito , que se debe a Euler, [8] converge para todos los números complejos excepto los enteros no positivos, que fallan debido a una división por cero. Intuitivamente, esta fórmula indica que es aproximadamente el resultado de calcular un número entero grande , multiplicar por para aproximar y usar la relación hacia atrás para obtener una aproximación ; y además que esta aproximación se vuelve exacta a medida que aumenta hasta el infinito.

El producto infinito del recíproco es una función entera , convergente para todo número complejo z .

Definición de Weierstrass

La definición de la función gamma debida a Weierstrass también es válida para todos los números complejos  z excepto los enteros no positivos: donde está la constante de Euler-Mascheroni . [1] Este es el producto de Hadamard en una forma reescrita. No se puede exagerar la utilidad de esta definición, ya que aparece en una determinada identidad que involucra a pi. [ cita necesaria ]

Propiedades

General

Además de la propiedad fundamental discutida anteriormente: otras ecuaciones funcionales importantes para la función gamma son la fórmula de reflexión de Euler que implica y la fórmula de duplicación de Legendre.

La fórmula de duplicación es un caso especial del teorema de la multiplicación (ver  [9] Ec. 5.5.6):

Una propiedad simple pero útil, que puede verse en la definición del límite, es:

En particular, con z = a + bi , este producto es

Si la parte real es un número entero o un semientero, esto se puede expresar de forma finita en forma cerrada :

Quizás el valor más conocido de la función gamma en un argumento no entero es el que se puede encontrar estableciendo las fórmulas de reflexión o duplicación, usando la relación con la función beta que se proporciona a continuación con , o simplemente haciendo la sustitución en el Definición integral de la función gamma, lo que da como resultado una integral gaussiana . En general, para valores enteros no negativos de tenemos: donde el doble factorial . Consulte Valores particulares de la función gamma para conocer los valores calculados.

Podría resultar tentador generalizar el resultado de que al buscar una fórmula para otros valores individuales donde sea racional, especialmente porque según el teorema digamma de Gauss , es posible hacerlo para la función digamma estrechamente relacionada en cada valor racional. Sin embargo, no se sabe que estos números sean expresables por sí mismos en términos de funciones elementales. Se ha demostrado que es un número trascendental y algebraicamente independiente de cualquier número entero y de cada una de las fracciones . [10] En general, al calcular los valores de la función gamma, debemos conformarnos con aproximaciones numéricas.

Las derivadas de la función gamma se describen en términos de la función poligammaψ (0) ( z ) : Para un entero positivo  m, la derivada de la función gamma se puede calcular de la siguiente manera:

Función gamma en el plano complejo con colores mostrando su argumento
Colores que muestran el argumento de la función gamma en el plano complejo desde −2 − 2 i hasta 6 + 2 i

donde H(m) es el número armónico m y γ es la constante de Euler-Mascheroni .

Para la derivada enésima de la función gamma es: (Esto se puede derivar diferenciando la forma integral de la función gamma con respecto a y usando la técnica de diferenciación bajo el signo integral ).

Usando la identidad donde es la función zeta de Riemann y es el -ésimo polinomio de Bell , tenemos en particular la expansión en serie de Laurent de la función gamma [11]

Desigualdades

Cuando se restringe a los números reales positivos, la función gamma es una función estrictamente logarítmicamente convexa . Esta propiedad podrá expresarse de cualquiera de las tres formas equivalentes siguientes:

La última de estas afirmaciones es, esencialmente por definición, la misma que la afirmación de que , donde está la función poligamma de orden 1. Para demostrar la convexidad logarítmica de la función gamma, basta observar que tiene una representación en serie que, por positivo real x , consta únicamente de términos positivos.

La convexidad logarítmica y la desigualdad de Jensen juntas implican, para cualquier número real positivo y ,

También existen límites en las proporciones de funciones gamma. La más conocida es la desigualdad de Gautschi , que dice que para cualquier número real positivo x y cualquier s ∈ (0, 1) ,

La fórmula de Stirling.

Representación de la función gamma en el plano complejo. Cada punto está coloreado según el argumento de . También se muestra el gráfico de contorno del módulo .
Gráfico tridimensional del valor absoluto de la función gamma compleja

El comportamiento de para una variable real positiva creciente viene dado por la fórmula de Stirling donde el símbolo significa convergencia asintótica: la relación de los dos lados converge a 1 en el límite . [1] Este crecimiento es más rápido que exponencial, , para cualquier valor fijo de .

Otro límite útil para aproximaciones asintóticas es:

Al escribir el término de error como un producto infinito, se puede utilizar la fórmula de Stirling para definir la función gamma: [12]

Residuos

El comportamiento de los no positivos es más complejo. La integral de Euler no converge para , pero la función que define en el semiplano complejo positivo tiene una continuación analítica única en el semiplano negativo. Una forma de encontrar esa continuación analítica es utilizar la integral de Euler para argumentos positivos y extender el dominio a números negativos mediante la aplicación repetida de la fórmula de recurrencia, [1] eligiendo los que sean positivos. El producto en el denominador es cero cuando es igual a cualquiera de los números enteros . Así, la función gamma debe estar indefinida en esos puntos para evitar la división por cero ; es una función meromórfica con polos simples en los números enteros no positivos. [1]

Para una función de una variable compleja , en un polo simple , el residuo de viene dado por:

Para el polo simple reescribimos la fórmula de recurrencia como: El numerador en es y el denominador Entonces los residuos de la función gamma en esos puntos son: [13] La función gamma es distinta de cero en todas partes a lo largo de la línea real, aunque se acerca arbitrariamente a cero cuando z → −∞ . De hecho, no existe ningún número complejo para el cual y, por tanto, la función gamma recíproca es una función completa , con ceros en . [1]

Mínimos y máximos

En la línea real, la función gamma tiene un mínimo local en z min+1.46163 21449 68362 34126 [14] donde alcanza el valor Γ( z min ) ≈ +0.88560 31944 10888 70027 . [15] La función gamma se eleva a ambos lados de este mínimo. La solución a Γ( z − 0,5) = Γ( z + 0,5) es z = +1,5 y el valor común es Γ(1) = Γ(2) = +1 . La solución positiva a Γ( z − 1) = Γ( z + 1) es z = φ ≈ +1.618 , la proporción áurea , y el valor común es Γ( φ − 1) = Γ( φ + 1) = φ . ≈ +1,44922 96022 69896 60037 . [dieciséis]

La función gamma debe alternar el signo entre sus polos en los números enteros no positivos porque el producto en la recurrencia directa contiene un número impar de factores negativos si el número de polos entre y es impar, y un número par si el número de polos es par. . [13] Los valores en los extremos locales de la función gamma a lo largo del eje real entre los números enteros no positivos son:

Γ( −0,50408 30082 64455 40925... [17] ) = −3,54464 36111 55005 08912... ,
Γ( −1,57349 84731 62390 45877... [18] ) = 2,30240 72583 39680 13582... ,
Γ( −2,61072 08684 44144 65000... [19] ) = −0,88813 63584 01241 92009... ,
Γ( −3,63529 33664 36901 09783... [20] ) = 0,24512 75398 34366 25043... ,
Γ( −4,65323 77617 43142 44171... [21] ) = −0,05277 96395 87319 40076... , etc.

Representaciones integrales

Hay muchas fórmulas, además de la integral de Euler de segundo tipo, que expresan la función gamma como una integral. Por ejemplo, cuando la parte real de z es positiva, [22] y [23] donde las tres integrales se derivan respectivamente de las sustituciones , [24] y [25] en la segunda integral de Euler. La última integral en particular deja clara la conexión entre la función gamma en argumentos semienteros y la integral gaussiana : si dejamos, obtenemos .

La primera fórmula integral de Binet para la función gamma establece que, cuando la parte real de z es positiva, entonces: [26] La integral del lado derecho puede interpretarse como una transformada de Laplace . Eso es,

La segunda fórmula integral de Binet establece que, nuevamente, cuando la parte real de z es positiva, entonces: [27]

Sea C un contorno de Hankel , es decir, un camino que comienza y termina en el punto de la esfera de Riemann , cuyo vector unitario tangente converge a −1 al inicio del camino y a 1 al final, que tiene el devanado número 1 alrededor 0 , y que no cruza [0, ∞) . Fije una rama de tomando una rama cortada a lo largo de [0, ∞) y tomando como real cuando t está en el eje real negativo. Supongamos que z no es un número entero. Entonces la fórmula de Hankel para la función gamma es: [28] donde se interpreta como . La fórmula de reflexión conduce a la expresión estrechamente relacionada que vuelve a ser válida siempre que z no sea un número entero.

Representación de fracción continua

La función gamma también se puede representar mediante la suma de dos fracciones continuas : [29] [30] donde .

Expansión en serie de Fourier

El logaritmo de la función gamma tiene la siguiente expansión en serie de Fourier , que durante mucho tiempo se atribuyó a Ernst Kummer , quien la derivó en 1847. [31] [32] Sin embargo, Iaroslav Blagouchine descubrió que Carl Johan Malmsten derivó por primera vez esta serie en 1842. [33] [34]

La fórmula de Raabe

En 1840, Joseph Ludwig Raabe demostró que , en particular, si entonces

Este último se puede derivar tomando el logaritmo en la fórmula de multiplicación anterior, que da una expresión para la suma de Riemann del integrando. Tomando el límite para da la fórmula.

función pi

Una notación alternativa que fue introducida originalmente por Gauss es la función -, que, en términos de la función gamma, es tal que para todo número entero no negativo .

Usando la función pi, la fórmula de reflexión toma la forma donde sinc es la función sinc normalizada , mientras que el teorema de la multiplicación toma la forma

A veces también encontramos que es una función completa , definida para cada número complejo, al igual que la función gamma recíproca . Eso es completo implica que no tiene polos, por lo que no tiene ceros .

El volumen de un n -elipsoide con radios r 1 ,…, r n se puede expresar como

Relación con otras funciones

Valores particulares

Incluyendo hasta los primeros 20 dígitos después del punto decimal, algunos valores particulares de la función gamma son: (Estos números se pueden encontrar en la OEIS . [35] [36] [37] [38] [39] [40] Los Los valores presentados aquí están truncados en lugar de redondeados). La función gamma de valores complejos no está definida para números enteros no positivos, pero en estos casos el valor se puede definir en la esfera de Riemann como . La función gamma recíproca está bien definida y es analítica en estos valores (y en todo el plano complejo ):

Función log-gamma

La función analítica logΓ( z )

Debido a que las funciones gamma y factorial crecen tan rápidamente para argumentos moderadamente grandes, muchos entornos informáticos incluyen una función que devuelve el logaritmo natural de la función gamma (a menudo se le da el nombre lgammaen lngammaentornos de programación o gammalnen hojas de cálculo); esto crece mucho más lentamente y, para cálculos combinatorios, permite sumar y restar registros en lugar de multiplicar y dividir valores muy grandes. A menudo se define como [41]

La función digamma , que es la derivada de esta función, también se ve comúnmente. En el contexto de aplicaciones técnicas y físicas, por ejemplo, con la propagación de ondas, la ecuación funcional

Función gamma logarítmica en el plano complejo desde −2 − 2i hasta 2 + 2i con colores
Función gamma logarítmica en el plano complejo desde −2 − 2i hasta 2 + 2i con colores

se utiliza a menudo ya que permite determinar valores de función en una franja de ancho 1 en z de la franja vecina. En particular, comenzando con una buena aproximación para  az con una parte real grande, se puede ir paso a paso hasta la  z deseada . Siguiendo una indicación de Carl Friedrich Gauss , Rocktaeschel (1922) propuso para logΓ( z ) una aproximación para Re( z ) grande :

Esto se puede utilizar para aproximar con precisión logΓ( z ) para z con un Re( z ) más pequeño mediante (PEBöhmer, 1939)

Se puede obtener una aproximación más precisa utilizando más términos de las expansiones asintóticas de logΓ( z ) y Γ( z ) , que se basan en la aproximación de Stirling.

como | z | → ∞ en constante | arg( z ) | < π . (Ver secuencias A001163 y A001164 en OEIS ).

En una presentación más "natural":

como | z | → ∞ en constante | arg( z ) | < π . (Ver secuencias A046968 y A046969 en OEIS ).

Los coeficientes de los términos con k > 1 de z 1− k en la última expansión son simplemente donde los B k son los números de Bernoulli .

La función gamma también tiene una Serie de Stirling (derivada por Charles Hermite en 1900) igual a [42]

Propiedades

El teorema de Bohr-Mollerup establece que entre todas las funciones que extienden las funciones factoriales a los números reales positivos, sólo la función gamma es log-convexa , es decir, su logaritmo natural es convexo en el eje real positivo. Otra caracterización viene dada por el teorema de Wielandt .

La función gamma es la única función que satisface simultáneamente

  1. ,
  2. para todos los números complejos excepto los enteros no positivos, y,
  3. para entero n , para todos los números complejos . [1]

En cierto sentido, la función log-gamma es la forma más natural; aclara algunos atributos intrínsecos de la función. Un ejemplo sorprendente es la serie de Taylor de logΓ alrededor de 1: donde ζ ( k ) denota la función zeta de Riemann en k .

Entonces, usando la siguiente propiedad: podemos encontrar una representación integral para la función log-gamma: o, estableciendo z = 1 para obtener una integral para γ , podemos reemplazar el término γ con su integral e incorporarlo a la fórmula anterior, Llegar:

También existen fórmulas especiales para el logaritmo de la función gamma para z racional . Por ejemplo, si y son números enteros con y entonces [43] Esta fórmula se utiliza a veces para cálculos numéricos, ya que el integrando disminuye muy rápidamente.

Integración sobre log-gamma

La integral se puede expresar en términos de la función G de Barnes [44] [45] (ver función G de Barnes para una prueba): donde Re( z ) > −1 .

También se puede escribir en términos de la función zeta de Hurwitz : [46] [47]

Cuando se sigue eso y esto también es una consecuencia de la fórmula de Raabe . O. Espinosa y V. Moll derivaron una fórmula similar para la integral del cuadrado de : [48] donde es .

DH Bailey y sus coautores [49] evaluaron cuándo en términos de la función zeta de Tornheim-Witten y sus derivadas.

Además, también se sabe que [50]

Aproximaciones

Comparación gamma (línea azul) con el factorial (puntos azules) y la aproximación de Stirling (línea roja)

Los valores complejos de la función gamma se pueden aproximar utilizando la aproximación de Stirling o la aproximación de Lanczos . Esto es preciso en el sentido de que la relación entre la aproximación y el valor verdadero se acerca a 1 en el límite como | z | va al infinito.

La función gamma se puede calcular con una precisión fija aplicando la integración por partes a la integral de Euler. Para cualquier número positivo  x, la función gamma se puede escribir

Cuando Re( z ) ∈ [1,2] y , el valor absoluto de la última integral es menor que . Al elegir un valor lo suficientemente grande , esta última expresión se puede hacer más pequeña que cualquier valor deseado . Por lo tanto, la función gamma se puede evaluar con precisión con la serie anterior.

EA Karatsuba construyó un algoritmo rápido para calcular la función gamma de Euler para cualquier argumento algebraico (incluido el racional). [51] [52] [53]

Para argumentos que son múltiplos enteros de 1/24 , la función gamma también se puede evaluar rápidamente utilizando iteraciones de medias aritméticas-geométricas (ver valores particulares de la función gamma ). [54]

Implementaciones prácticas

A diferencia de muchas otras funciones, como la Distribución Normal , no se encuentra fácilmente una implementación obvia, rápida y precisa que sea fácil de implementar para la Función Gamma . Por lo tanto, vale la pena investigar posibles soluciones. En el caso de que la velocidad sea más importante que la precisión, las tablas publicadas se pueden encontrar fácilmente mediante una búsqueda en Internet, como la Biblioteca Wiley en línea. Estas tablas se pueden utilizar con interpolación lineal . Se puede obtener una mayor precisión con el uso de la interpolación cúbica a costa de una mayor sobrecarga computacional. Dado que las tablas generalmente se publican para valores de argumento entre 1 y 2, la propiedad se puede usar para traducir rápida y fácilmente todos los valores reales al rango , de modo que solo sea necesario usar valores tabulados entre 1 y 2. [55]

Si las tablas de interpolación no son deseables, entonces la aproximación de Lanczos mencionada anteriormente funciona bien para obtener una precisión de 1 a 2 dígitos para valores pequeños y comúnmente utilizados de z. Si la aproximación de Lanczos no es lo suficientemente precisa, se puede utilizar la fórmula de Stirling para la función gamma .

Aplicaciones

Un autor describe la función gamma como "Posiblemente, la función especial más común, o la menos 'especial' de ellas. Las otras funciones trascendentales [...] se llaman 'especiales' porque es posible evitar algunas de ellas al mantenerse alejado de muchas temas matemáticos especializados. Por otro lado, la función gamma Γ( z ) es la más difícil de evitar." [56]

Problemas de integración

La función gamma encuentra aplicación en campos tan diversos como la física cuántica , la astrofísica y la dinámica de fluidos . [57] La ​​distribución gamma , que se formula en términos de la función gamma, se utiliza en estadística para modelar una amplia gama de procesos; por ejemplo, el tiempo entre ocurrencias de terremotos. [58]

La razón principal de la utilidad de la función gamma en tales contextos es la prevalencia de expresiones del tipo que describen procesos que decaen exponencialmente en el tiempo o el espacio. En ocasiones, las integrales de tales expresiones se pueden resolver en términos de la función gamma cuando no existe una solución elemental. Por ejemplo, si f es una función potencia y g es una función lineal, un simple cambio de variables da la evaluación

El hecho de que la integración se realice a lo largo de toda la línea real positiva podría significar que la función gamma describe la acumulación de un proceso dependiente del tiempo que continúa indefinidamente, o que el valor podría ser el total de una distribución en un espacio infinito.

Por supuesto, con frecuencia es útil tomar límites de integración distintos de 0 y para describir la acumulación de un proceso finito, en cuyo caso la función gamma ordinaria ya no es una solución; la solución entonces se llama función gamma incompleta . (La función gamma ordinaria, obtenida integrando toda la línea real positiva, a veces se denomina función gamma completa para contrastar).

Una categoría importante de funciones que decaen exponencialmente es la de las funciones gaussianas y sus integrales, como la función de error . Existen muchas interrelaciones entre estas funciones y la función gamma; en particular, el factor obtenido al evaluar es el "mismo" que se encuentra en el factor de normalización de la función de error y la distribución normal .

Las integrales que hemos comentado hasta ahora involucran funciones trascendentales , pero la función gamma también surge de integrales de funciones puramente algebraicas. En particular, las longitudes de arco de las elipses y de la lemniscata , que son curvas definidas por ecuaciones algebraicas, vienen dadas por integrales elípticas que en casos especiales pueden evaluarse en términos de la función gamma. La función gamma también se puede utilizar para calcular el "volumen" y el "área" de hiperesferas de n dimensiones .

Calcular productos

La capacidad de la función gamma para generalizar productos factoriales conduce inmediatamente a aplicaciones en muchas áreas de las matemáticas; en combinatoria , y por extensión en áreas como la teoría de la probabilidad y el cálculo de series de potencias . Muchas expresiones que involucran productos de números enteros sucesivos se pueden escribir como alguna combinación de factoriales, siendo quizás el ejemplo más importante el del coeficiente binomial . Por ejemplo, para cualquier número complejo z y n , con | z | < 1 , podemos escribir cuál se parece mucho al coeficiente binomial cuando n es un número entero no negativo,

El ejemplo de los coeficientes binomiales motiva por qué las propiedades de la función gamma cuando se extiende a números negativos son naturales. Un coeficiente binomial indica el número de formas de elegir k elementos de un conjunto de n elementos; Si k > n , por supuesto no hay formas. Si k > n , ( nk )! es el factorial de un entero negativo y, por tanto, infinito si utilizamos la definición de factoriales de la función gamma: dividir por infinito da el valor esperado de 0.

Podemos reemplazar el factorial por una función gamma para extender cualquier fórmula de este tipo a los números complejos. Generalmente, esto funciona para cualquier producto en el que cada factor sea una función racional de la variable índice, factorizando la función racional en expresiones lineales. Si P y Q son polinomios mónicos de grado myn con raíces respectivas p 1 , …, p m y q 1 , …, q n , tenemos

Si tenemos una forma de calcular numéricamente la función gamma, es muy sencillo calcular los valores numéricos de dichos productos. El número de funciones gamma en el lado derecho depende sólo del grado de los polinomios, por lo que no importa si ba es igual a 5 o 10 5 . Al tomar los límites apropiados, también se puede hacer que la ecuación se cumpla incluso cuando el producto de la izquierda contenga ceros o polos.

Al tomar límites, ciertos productos racionales con infinitos factores también pueden evaluarse en términos de la función gamma. Debido al teorema de factorización de Weierstrass , las funciones analíticas se pueden escribir como productos infinitos y, en ocasiones, se pueden representar como productos finitos o cocientes de la función gamma. Ya hemos visto un ejemplo sorprendente: la fórmula de reflexión esencialmente representa la función seno como el producto de dos funciones gamma. A partir de esta fórmula, la función exponencial, así como todas las funciones trigonométricas e hiperbólicas, se pueden expresar en términos de la función gamma.

Se pueden representar más funciones aún, incluida la función hipergeométrica y sus casos especiales, mediante integrales de contorno complejas de productos y cocientes de la función gamma, llamadas integrales de Mellin-Barnes .

Teoría analítica de números

Una aplicación de la función gamma es el estudio de la función zeta de Riemann . Una propiedad fundamental de la función zeta de Riemann es su ecuación funcional :

Esto proporciona, entre otras cosas, una forma explícita para la continuación analítica de la función zeta a una función meromorfa en el plano complejo y conduce a una prueba inmediata de que la función zeta tiene infinitos ceros llamados "triviales" en la recta real. Borwein et al. llama a esta fórmula "uno de los hallazgos más bellos de las matemáticas". [59] Otro contendiente para ese título podría ser

Ambas fórmulas fueron derivadas por Bernhard Riemann en su artículo fundamental de 1859 " Ueber die Anzahl der Primzahlen unter einer gegebenen Größe " ("Sobre el número de primos menores que una magnitud dada"), uno de los hitos en el desarrollo de la teoría analítica de números . Rama de las matemáticas que estudia los números primos utilizando las herramientas del análisis matemático.

Historia

La función gamma ha despertado el interés de algunos de los matemáticos más destacados de todos los tiempos. Su historia, notablemente documentada por Philip J. Davis en un artículo que le valió el Premio Chauvenet en 1963 , refleja muchos de los principales avances ocurridos en las matemáticas desde el siglo XVIII. En palabras de Davis, "cada generación ha encontrado algo interesante que decir sobre la función gamma. Quizás la próxima generación también lo haga". [1]

Siglo XVIII: Euler y Stirling

Carta de Daniel Bernoulli a Christian Goldbach , 6 de octubre de 1729

El problema de extender el factorial a argumentos no enteros aparentemente fue considerado por primera vez por Daniel Bernoulli y Christian Goldbach en la década de 1720. En particular, en una carta de Bernoulli a Goldbach fechada el 6 de octubre de 1729, Bernoulli introdujo la representación del producto [60] que está bien definida para valores reales de x distintos de los números enteros negativos.

Leonhard Euler dio más tarde dos definiciones diferentes: la primera no era su integral sino un producto infinito bien definido para todos los números complejos n distintos de los enteros negativos, de lo que informó a Goldbach en una carta fechada el 13 de octubre de 1729. Le escribió a Goldbach. nuevamente el 8 de enero de 1730, para anunciar su descubrimiento de la representación integral que es válida cuando la parte real del número complejo n es estrictamente mayor que −1 (es decir, ). Por el cambio de variables t = −ln s , esto se convierte en la familiar integral de Euler. Euler publicó sus resultados en el artículo "De progresivoibus trascendentibus seu quarum termini generales algebraice dari nequeunt" ("Sobre progresiones trascendentales, es decir, aquellas cuyos términos generales no pueden darse algebraicamente"), presentado a la Academia de San Petersburgo el 28 de noviembre de 1729. [ 61] Euler descubrió además algunas de las propiedades funcionales importantes de la función gamma, incluida la fórmula de reflexión.

James Stirling , contemporáneo de Euler, también intentó encontrar una expresión continua para el factorial y ideó lo que hoy se conoce como fórmula de Stirling . Aunque la fórmula de Stirling da una buena estimación de n ! , además, para números no enteros, no proporciona el valor exacto. Extensiones de su fórmula que corrigen el error fueron dadas por el propio Stirling y por Jacques Philippe Marie Binet .

Siglo XIX: Gauss, Weierstrass y Legendre

De progresivoibus trascendentibus, seu quarum termini generales algebraicae dari nequeunt
La primera página del artículo de Euler.

Carl Friedrich Gauss reescribió el producto de Euler y utilizó esta fórmula para descubrir nuevas propiedades de la función gamma. Aunque Euler fue un pionero en la teoría de variables complejas, no parece haber considerado el factorial de un número complejo, como lo hizo Gauss por primera vez. [62] Gauss también demostró el teorema de multiplicación de la función gamma e investigó la conexión entre la función gamma y las integrales elípticas .

Karl Weierstrass estableció además el papel de la función gamma en el análisis complejo , a partir de otra representación del producto, donde γ es la constante de Euler-Mascheroni . Weierstrass originalmente escribió su producto como uno para 1/Γ , en cuyo caso se toman los ceros de la función en lugar de sus polos. Inspirado por este resultado, demostró lo que se conoce como teorema de factorización de Weierstrass : que cualquier función completa puede escribirse como un producto de sus ceros en el plano complejo; una generalización del teorema fundamental del álgebra .

El nombre función gamma y el símbolo Γ fueron introducidos por Adrien-Marie Legendre alrededor de 1811; Legendre también reescribió la definición integral de Euler en su forma moderna. Aunque el símbolo es una gamma griega en mayúsculas, no existe un estándar aceptado sobre si el nombre de la función debe escribirse "función gamma" o "función gamma" (algunos autores simplemente escriben " Γ -función"). La notación alternativa de "función pi" Π( z ) = z ! debida a Gauss se encuentra a veces en la literatura más antigua, pero la notación de Legendre es dominante en las obras modernas.

Está justificado preguntar por qué distinguimos entre el "factorial ordinario" y la función gamma usando símbolos distintos, y particularmente por qué la función gamma debería normalizarse a Γ( n + 1) = n . en lugar de simplemente usar " Γ( n ) = n ! ". Considere que la notación para exponentes, x n , se ha generalizado desde números enteros hasta números complejos x z sin ningún cambio. La motivación de Legendre para la normalización no parece ser conocida, y ha sido criticada como engorrosa por algunos (el matemático del siglo XX Cornelius Lanczos , por ejemplo, la llamó "desprovista de cualquier racionalidad" y en su lugar usaría z ! ). [63] La normalización de Legendre simplifica algunas fórmulas, pero complica otras. Desde un punto de vista moderno, la normalización de Legendre de la función gamma es la integral del carácter aditivo e x contra el carácter multiplicativo x z con respecto a la medida de Haar en el grupo de Lie R + . Así, esta normalización deja más claro que la función gamma es un análogo continuo de una suma de Gauss . [64]

Siglos XIX y XX: caracterización de la función gamma

Es algo problemático que se hayan dado tantas definiciones para la función gamma. Aunque describen la misma función, no es del todo sencillo demostrar la equivalencia. Stirling nunca demostró que su fórmula extendida corresponda exactamente a la función gamma de Euler; Charles Hermite dio por primera vez una prueba en 1900. [65] En lugar de encontrar una prueba especializada para cada fórmula, sería deseable tener un método general para identificar la función gamma.

Una forma de demostrar la equivalencia sería encontrar una ecuación diferencial que caracterice la función gamma. La mayoría de las funciones especiales en matemáticas aplicadas surgen como soluciones de ecuaciones diferenciales, cuyas soluciones son únicas. Sin embargo, la función gamma no parece satisfacer ninguna ecuación diferencial simple. Otto Hölder demostró en 1887 que la función gamma al menos no satisface ninguna ecuación diferencial algebraica al mostrar que una solución a tal ecuación no podría satisfacer la fórmula de recurrencia de la función gamma, convirtiéndola en una función trascendentalmente trascendental . Este resultado se conoce como teorema de Hölder .

Hasta 1922 no se dio una caracterización definitiva y generalmente aplicable de la función gamma. Harald Bohr y Johannes Mollerup demostraron entonces lo que se conoce como teorema de Bohr-Mollerup : que la función gamma es la única solución de la relación de recurrencia factorial que es positiva y logarítmicamente convexa para z positiva y cuyo valor en 1 es 1 (una función es logarítmicamente convexa si su logaritmo es convexo). Otra caracterización viene dada por el teorema de Wielandt .

El teorema de Bohr-Mollerup es útil porque es relativamente fácil demostrar la convexidad logarítmica de cualquiera de las diferentes fórmulas utilizadas para definir la función gamma. Yendo más allá, en lugar de definir la función gamma mediante una fórmula particular, podemos elegir las condiciones del teorema de Bohr-Mollerup como definición y luego elegir cualquier fórmula que nos guste que satisfaga las condiciones como punto de partida para estudiar la función gamma. . Este enfoque fue utilizado por el grupo Bourbaki .

Borwein y Corless [66] revisan tres siglos de trabajo sobre la función gamma.

Tablas de referencia y software.

Aunque la función gamma se puede calcular prácticamente tan fácilmente como cualquier función matemáticamente más simple con una computadora moderna (incluso con una calculadora de bolsillo programable), por supuesto, este no siempre fue así. Hasta mediados del siglo XX, los matemáticos dependían de tablas hechas a mano; en el caso de la función gamma, en particular una tabla calculada por Gauss en 1813 y otra calculada por Legendre en 1825. [67]

Un gráfico dibujado a mano del valor absoluto de la función gamma compleja, de Tablas de funciones superiores de Jahnke y Emde  [de] .

Las tablas de valores complejos de la función gamma, así como gráficos dibujados a mano, se dieron en Tablas de funciones con fórmulas y curvas de Jahnke y Emde  [de] , publicado por primera vez en Alemania en 1909. Según Michael Berry , "la publicación en J&E de un gráfico tridimensional que muestra los polos de la función gamma en el plano complejo adquirió un estatus casi icónico." [68]

De hecho, había poca necesidad práctica de algo más que valores reales de la función gamma hasta la década de 1930, cuando se descubrieron aplicaciones de la función gamma compleja en la física teórica. A medida que las computadoras electrónicas estuvieron disponibles para la producción de tablas en la década de 1950, se publicaron varias tablas extensas para la función gamma compleja para satisfacer la demanda, incluida una tabla con una precisión de 12 decimales de la Oficina Nacional de Estándares de EE. UU . [1]

Reproducción de un famoso diagrama complejo de Janhke y Emde (Tablas de funciones con fórmulas y curvas, 4ª ed., Dover, 1945) de la función gamma de −4,5 − 2,5i a 4,5 + 2,5i
Reproducción de un famoso diagrama complejo de Janhke y Emde (Tablas de funciones con fórmulas y curvas, 4ª ed., Dover, 1945) de la función gamma de −4,5 − 2,5i a 4,5 + 2,5i

Las implementaciones de punto flotante de doble precisión de la función gamma y su logaritmo ahora están disponibles en la mayoría de los software de computación científica y bibliotecas de funciones especiales, por ejemplo TK Solver , Matlab , GNU Octave y la Biblioteca Científica GNU . La función gamma también se agregó a la biblioteca estándar de C ( math.h ). Las implementaciones de precisión arbitraria están disponibles en la mayoría de los sistemas de álgebra informática , como Mathematica y Maple . PARI/GP , MPFR y MPFUN contienen implementaciones gratuitas de precisión arbitraria. En algunas calculadoras de software , por ejemplo, la Calculadora de Windows y la Calculadora de GNOME , la función factorial devuelve Γ( x + 1) cuando la entrada x es un valor no entero. [69] [70]

Ver también

Notas

  1. ^ abcdefghijklmn Davis, PJ (1959). "La integral de Leonhard Euler: un perfil histórico de la función gamma". Mensual Matemático Estadounidense . 66 (10): 849–869. doi :10.2307/2309786. JSTOR  2309786. Archivado desde el original el 7 de noviembre de 2012 . Consultado el 3 de diciembre de 2016 .
  2. ^ "¿Está mal definida la función Gamma? O: Hadamard versus Euler: ¿quién encontró la mejor función Gamma?".
  3. ^ Beals, Richard; Wong, Roderick (2010). Funciones especiales: un texto de posgrado. Prensa de la Universidad de Cambridge. pag. 28.ISBN 978-1-139-49043-6.Extracto de la página 28
  4. ^ Ross, arcilla C. (2013). Ecuaciones diferenciales: una introducción a Mathematica (edición ilustrada). Medios de ciencia y negocios de Springer. pag. 293.ISBN 978-1-4757-3949-7.Expresión G.2 en la página 293
  5. ^ Kingman, JFC (1961). "Una propiedad de convexidad de matrices positivas". La Revista Trimestral de Matemáticas . 12 (1): 283–284. Código Bib : 1961QJMat..12..283K. doi :10.1093/qmath/12.1.283.
  6. ^ Weisstein, Eric W. "Teorema de Bohr-Mollerup". MundoMatemático .
  7. ^ Davis, Felipe. "La integral de Leonhard Euler: un perfil histórico de la función gamma" (PDF) . maa.org .
  8. ^ Bonvini, Marco (9 de octubre de 2010). "La función Gamma" (PDF) . Roma1.infn.it .
  9. ^ Askey, RA ; Roy, R. (2010), "Series Expansions", en Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), Manual de funciones matemáticas del NIST , Cambridge University Press, ISBN 978-0-521-19225-5, señor  2723248.
  10. ^ Waldschmidt, M. (2006). "Trascendencia de los períodos: el estado del arte" (PDF) . Pura aplicación. Matemáticas. Cuarto . 2 (2): 435–463. doi : 10.4310/pamq.2006.v2.n2.a3 . Archivado (PDF) desde el original el 6 de mayo de 2006.
  11. ^ "¿Cómo obtener la expansión de Laurent de la función gamma alrededor de $z=0$?". Intercambio de pilas de matemáticas . Consultado el 17 de agosto de 2022 .
  12. ^ Artín, Emil (2015). La función gamma . Dover. pag. 24.
  13. ^ ab Weisstein, Eric W. "Función gamma". MundoMatemático .
  14. ^ Sloane, Nueva Jersey (ed.). "Secuencia A030169 (Expansión decimal del número real x tal que y = Gamma(x) es un mínimo)". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  15. ^ Sloane, Nueva Jersey (ed.). "Secuencia A030171 (Expansión decimal del número real y tal que y = Gamma(x) es un mínimo)". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  16. ^ Sloane, Nueva Jersey (ed.). "Secuencia A178840 (Expansión decimal del factorial de la proporción áurea)". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  17. ^ Sloane, Nueva Jersey (ed.). "Secuencia A175472 (Expansión decimal del valor absoluto de la abscisa del máximo local de la función Gamma en el intervalo [-1,0])". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  18. ^ Sloane, Nueva Jersey (ed.). "Secuencia A175473 (Expansión decimal del valor absoluto de la abscisa del mínimo local de la función Gamma en el intervalo [-2,-1])". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  19. ^ Sloane, Nueva Jersey (ed.). "Secuencia A175474 (Expansión decimal del valor absoluto de la abscisa del máximo local de la función Gamma en el intervalo [-3,-2])". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  20. ^ Sloane, Nueva Jersey (ed.). "Secuencia A256681 (Expansión decimal de la abscisa [negada] del mínimo local de la función Gamma en el intervalo [-4,-3])". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  21. ^ Sloane, Nueva Jersey (ed.). "Secuencia A256682 (Expansión decimal de la abscisa [negada] del máximo local de la función Gamma en el intervalo [-5,-4])". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  22. ^ Gradshteyn, ES; Ryzhik, IM (2007). Tabla de integrales, series y productos (Séptima ed.). Prensa académica. pag. 893.ISBN 978-0-12-373637-6.
  23. ^ Whittaker y Watson, 12.2 ejemplo 1.
  24. ^ Detlef, Gronau. "¿Por qué la función gamma es así?" (PDF) . Imsc.uni-graz.at .
  25. ^ Pascal Sebah, Xavier Gourdon. «Introducción a la Función Gamma» (PDF) . Computación de números . Archivado desde el original (PDF) el 30 de enero de 2023 . Consultado el 30 de enero de 2023 .
  26. ^ Whittaker y Watson, 12.31.
  27. ^ Whittaker y Watson, 12.32.
  28. ^ Whittaker y Watson, 22.12.
  29. ^ "Integral exponencial E: representaciones de fracciones continuas (Fórmula 06.34.10.0005)".
  30. ^ "Integral exponencial E: representaciones de fracciones continuas (Fórmula 06.34.10.0003)".
  31. ^ Bateman, Harry; Erdélyi, Arthur (1955). Funciones trascendentales superiores . McGraw-Hill.
  32. ^ Srivastava, HM; Choi, J. (2001). Series Asociadas al Zeta y Funciones Afines . Países Bajos: Kluwer Academic.
  33. ^ Blagouchine, Iaroslav V. (2014). "Redescubrimiento de las integrales de Malmsten, su evaluación mediante métodos de integración de contornos y algunos resultados relacionados". Ramanujan J. 35 (1): 21-110. doi :10.1007/s11139-013-9528-5. S2CID  120943474.
  34. ^ Blagouchine, Iaroslav V. (2016). "Erratum y Anexo a "Redescubrimiento de las integrales de Malmsten, su evaluación mediante métodos de integración de contornos y algunos resultados relacionados"". Ramanujan J. 42 ( 3): 777–781. doi :10.1007/s11139-015-9763-z. S2CID  125198685.
  35. ^ Sloane, Nueva Jersey (ed.). "Secuencia A245886 (Expansión decimal de Gamma (-3/2), donde Gamma es la función gamma de Euler)". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  36. ^ Sloane, Nueva Jersey (ed.). "Secuencia A019707 (Expansión decimal de sqrt(Pi)/5)". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  37. ^ Sloane, Nueva Jersey (ed.). "Secuencia A002161 (Expansión decimal de la raíz cuadrada de Pi)". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  38. ^ Sloane, Nueva Jersey (ed.). "Secuencia A019704 (Expansión decimal de sqrt(Pi)/2)". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  39. ^ Sloane, Nueva Jersey (ed.). "Secuencia A245884 (Expansión decimal de Gamma (5/2), donde Gamma es la función gamma de Euler)". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  40. ^ Sloane, Nueva Jersey (ed.). "Secuencia A245885 (Expansión decimal de Gamma (7/2), donde Gamma es la función gamma de Euler)". La enciclopedia en línea de secuencias enteras . Fundación OEIS.
  41. ^ "Función de registro gamma". Wolfram MathWorld . Consultado el 3 de enero de 2019 .
  42. ^ "Integral de Leonhard Euler: un perfil histórico de la función gamma" (PDF) . Archivado (PDF) desde el original el 12 de septiembre de 2014 . Consultado el 11 de abril de 2022 .
  43. ^ Blagouchine, Iaroslav V. (2015). "Un teorema para la evaluación en forma cerrada de la primera constante de Stieltjes generalizada en argumentos racionales y algunas sumas relacionadas". Revista de teoría de números . 148 : 537–592. arXiv : 1401.3724 . doi :10.1016/j.jnt.2014.08.009.
  44. ^ Alexejewsky, WP (1894). "Über eine Classe von Funktionen, die der Gammafunktion analog sind" [Sobre una clase de funciones análogas a la función gamma]. Leipzig Weidmanncshe Buchhandluns . 46 : 268–275.
  45. ^ Barnes, EW (1899). "La teoría de la función G ". Cuarto de galón. J. Matemáticas . 31 : 264–314.
  46. ^ Adamchik, Víctor S. (1998). "Funciones poligamma de orden negativo". J. Computación. Aplica. Matemáticas . 100 (2): 191–199. doi : 10.1016/S0377-0427(98)00192-7 .
  47. ^ Gosper, RW (1997). " en funciones especiales, q -series y temas afines". Mermelada. Matemáticas. Soc . 14 .
  48. ^ Espinosa, Olivier; Moll, Víctor H. (2002). "Sobre algunas integrales que involucran la función Hurwitz Zeta: Parte 1". El diario Ramanujan . 6 (2): 159–188. doi :10.1023/A:1015706300169. S2CID  128246166.
  49. ^ Bailey, David H.; Borwein, David; Borwein, Jonathan M. (2015). "Sobre las integrales log-gamma de Euler y las funciones zeta de Tornheim-Witten". El diario Ramanujan . 36 (1–2): 43–68. doi :10.1007/s11139-012-9427-1. S2CID  7335291.
  50. ^ Amdeberhan, T.; Coffey, Mark W.; Espinosa, Olivier; Koutschan, Christoph; Maná, Dante V.; Moll, Víctor H. (2011). "Integrales de potencias de logamma". Proc. América. Matemáticas. Soc . 139 (2): 535–545. doi : 10.1090/S0002-9939-2010-10589-0 .
  51. ^ EA Karatsuba, Evaluación rápida de funciones trascendentales. Problema. inf. Trans. Vol.27, No.4, págs. 339–360 (1991).
  52. ^ EA Karatsuba, Sobre un nuevo método para la evaluación rápida de funciones trascendentales. Ruso. Matemáticas. Sobrevivir. Vol.46, n.º 2, págs. 246–247 (1991).
  53. ^ EA Karatsuba "Algoritmos rápidos y el método FEE".
  54. ^ Borwein, JM; Zucker, IJ (1992). "Evaluación rápida de la función gamma para pequeñas fracciones racionales utilizando integrales elípticas completas de primer tipo". Revista IMA de Análisis Numérico . 12 (4): 519–526. doi :10.1093/IMANUM/12.4.519.
  55. ^ Werner, Helmut; Collinge, Robert (1961). "Aproximaciones de Chebyshev a la función gamma". Matemáticas. Computación . 15 (74): 195-197. JSTOR  2004230.
  56. ^ Michon, GP "Trigonometría y funciones básicas Archivado el 9 de enero de 2010 en Wayback Machine ". Numérica . Consultado el 5 de mayo de 2007.
  57. ^ Chaudry, MA y Zubair, SM (2001). Sobre una clase de funciones gamma incompletas con aplicaciones . pag. 37
  58. ^ Arroz, JA (1995). Estadística Matemática y Análisis de Datos (Segunda Edición). pag. 52–53
  59. ^ Borwein, J.; Bailey, DH y Girgensohn, R. (2003). Experimentación en Matemáticas . AK Peters. pag. 133.ISBN 978-1-56881-136-9.
  60. ^ "Interpolando el factorial natural n! o El nacimiento de la función factorial real (1729 - 1826)".
  61. ^ El artículo de Euler fue publicado en Commentarii academiae scientiarum Petropolitanae 5, 1738, 36–57. Véase E19 - De progresivoibus trascendentibus seu quarum termini generales algebraice dari nequeunt, de The Euler Archive, que incluye una copia escaneada del artículo original.
  62. ^ Remmert, R. (2006). Temas clásicos de la teoría de funciones complejas . Traducido por Kay, LD Springer. ISBN 978-0-387-98221-2.
  63. ^ Lanczos, C. (1964). "Una aproximación de precisión de la función gamma". J. SIAM Número. Anal. Ser. B . 1 (1): 86. Código bibliográfico : 1964SJNA....1...86L. doi :10.1137/0701008.
  64. ^ Ilker Inam; Engin Büyükaşşk (2019). Notas de la Escuela Internacional de Otoño sobre Teoría Computacional de Números. Saltador. pag. 205.ISBN 978-3-030-12558-5.Extracto de la página 205
  65. ^ Knuth, DE (1997). El arte de la programación informática, volumen 1 (algoritmos fundamentales) . Addison-Wesley.
  66. ^ Borwein, Jonathan M .; Corless, Robert M. (2017). "Gamma y Factorial en la Mensual". Mensual Matemático Estadounidense . 125 (5). Asociación Matemática de América: 400–24. arXiv : 1703.05349 . Código Bib : 2017arXiv170305349B. doi :10.1080/00029890.2018.1420983. S2CID  119324101.
  67. ^ "¿Cuál es la historia de Gamma_function?". añois.com . Consultado el 5 de noviembre de 2022 .
  68. ^ Berry, M. (abril de 2001). "¿Por qué las funciones especiales son especiales?". Física hoy .
  69. ^ "microsoft/calculadora". GitHub . Consultado el 25 de diciembre de 2020 .
  70. ^ "calculadora de gnomos". GNOME.org . Consultado el 3 de marzo de 2023 .

Otras lecturas

enlaces externos