stringtranslate.com

Función gamma incompleta

La función gamma superior incompleta para algunos valores de s: 0 (azul), 1 (rojo), 2 (verde), 3 (naranja), 4 (púrpura).
Gráfico de la función gamma incompleta regularizada Q(2,z) en el plano complejo desde -2-2i hasta 2+2i con colores creados con la función ComplexPlot3D de Mathematica 13.1
Gráfico de la función gamma incompleta regularizada Q(2,z) en el plano complejo desde -2-2i hasta 2+2i con colores creados con la función ComplexPlot3D de Mathematica 13.1

En matemáticas , las funciones gamma incompletas superior e inferior son tipos de funciones especiales que surgen como soluciones a diversos problemas matemáticos, como ciertas integrales .

Sus respectivos nombres se derivan de sus definiciones integrales, que se definen de manera similar a la función gamma , pero con límites integrales diferentes o "incompletos". La función gamma se define como una integral de cero a infinito. Esto contrasta con la función gamma incompleta inferior, que se define como una integral de cero a un límite superior variable. De manera similar, la función gamma incompleta superior se define como una integral de un límite inferior variable a infinito.

Definición

La función gamma incompleta superior se define como: mientras que la función gamma incompleta inferior se define como: En ambos casos, s es un parámetro complejo, tal que la parte real de s es positiva.

Propiedades

Por integración por partes encontramos las relaciones de recurrencia y Dado que la función gamma ordinaria se define como tenemos y

Continuación de valores complejos

La función gamma incompleta inferior y la función gamma incompleta superior, como se definieron anteriormente para s y x positivos reales , se pueden desarrollar en funciones holomorfas , con respecto tanto a x como a s , definidas para casi todas las combinaciones de x y s complejos . [1] El análisis complejo muestra cómo las propiedades de las funciones gamma incompletas reales se extienden a sus contrapartes holomorfas.

Función gamma incompleta inferior

Extensión holomorfa

La aplicación repetida de la relación de recurrencia para la función gamma incompleta inferior conduce a la expansión en serie de potencias : [2] Dado el rápido crecimiento en valor absoluto de Γ( z + k ) cuando k → ∞ , y el hecho de que el recíproco de Γ( z ) es una función entera , los coeficientes en la suma más a la derecha están bien definidos, y localmente la suma converge uniformemente para todos los complejos s y x . Por un teorema de Weierstrass , [3] la función límite, a veces denotada como , [4] es entera con respecto tanto a z (para s fijo ) como a s (para z fijo ), [1] y, por lo tanto, holomorfa en C × C por el teorema de Hartog . [5] Por lo tanto, la siguiente descomposición [1] extiende la función gamma incompleta inferior real como una función holomorfa , tanto conjunta como separadamente en z y s . De las propiedades de y de la función Γ se deduce que los dos primeros factores capturan las singularidades de (en z = 0 o s un entero no positivo), mientras que el último factor contribuye a sus ceros.

Multivaloridad

El logaritmo complejo log z = log | z | + i arg z se determina únicamente hasta un múltiplo de 2 πi , lo que lo hace multivaluado . Las funciones que involucran el logaritmo complejo generalmente heredan esta propiedad. Entre ellas se encuentran la potencia compleja , y, dado que z s aparece en su descomposición, también la función γ .

La indeterminación de funciones multivaluadas introduce complicaciones, ya que se debe indicar cómo seleccionar un valor. Las estrategias para manejar esto son:

El siguiente conjunto de reglas se puede utilizar para interpretar correctamente las fórmulas de esta sección. A menos que se indique lo contrario, se supone lo siguiente:

Sectores

Los sectores en C que tienen su vértice en z = 0 a menudo resultan ser dominios apropiados para expresiones complejas. Un sector D consiste en todos los complejos z que satisfacen z ≠ 0 y αδ < arg z < α + δ con algún α y 0 < δπ . A menudo, α puede elegirse arbitrariamente y no se especifica entonces. Si δ no se da, se supone que es π , y el sector es de hecho todo el plano C , con la excepción de una semilínea que se origina en z = 0 y apunta en la dirección de α , que generalmente sirve como un corte de rama . Nota: En muchas aplicaciones y textos, α se toma silenciosamente como 0, lo que centra el sector alrededor del eje real positivo.

Sucursales

En particular, existe un logaritmo holomorfo y de un solo valor en cualquier sector D que tenga su parte imaginaria limitada al rango ( αδ , α + δ ) . Con base en un logaritmo tan restringido, z s y las funciones gamma incompletas a su vez colapsan en funciones holomorfas y de un solo valor en D (o C × D ), llamadas ramas de sus contrapartes de múltiples valores en D . Agregar un múltiplo de 2 π a α produce un conjunto diferente de ramas correlacionadas en el mismo conjunto D . Sin embargo, en cualquier contexto dado aquí, α se supone fijo y todas las ramas involucradas están asociadas a él. Si | α | < δ , las ramas se llaman principales , porque son iguales a sus análogos reales en el eje real positivo. Nota: En muchas aplicaciones y textos, las fórmulas solo son válidas para ramas principales.

Relación entre ramas

Los valores de diferentes ramas tanto de la función de potencia compleja como de la función gamma incompleta inferior se pueden derivar entre sí mediante la multiplicación de , [1] para k, un entero adecuado.

Comportamiento cerca del punto de ramificación

La descomposición anterior muestra además que γ se comporta cerca de z = 0 asintóticamente como:

Para x , y y s reales positivos , x y /y → 0 , cuando ( x , y ) → (0, s ) . Esto parece justificar establecer γ ( s , 0) = 0 para s real > 0 . Sin embargo, las cosas son algo diferentes en el ámbito complejo. Solo si (a) la parte real de s es positiva, y (b) los valores u v se toman de solo un conjunto finito de ramas, se garantiza que convergen a cero cuando ( u , v ) → (0, s ) , y lo mismo ocurre con γ ( u , v ) . En una sola rama de γ ( b ) se cumple naturalmente, por lo que γ ( s , 0) = 0 para s con parte real positiva es un límite continuo . También tenga en cuenta que dicha continuación no es de ninguna manera analítica .

Relaciones algebraicas

Todas las relaciones algebraicas y ecuaciones diferenciales observadas por la función γ ( s , z ) real se cumplen también para su contraparte holomorfa. Esto es una consecuencia del teorema de identidad, que establece que las ecuaciones entre funciones holomorfas válidas en un intervalo real se cumplen en todas partes. En particular, la relación de recurrencia [2] y ∂γ ( s , z )/ ∂z = z s −1 e z [2] se conservan en las ramas correspondientes.

Representación integral

La última relación nos dice que, para s fijo , γ es una primitiva o antiderivada de la función holomorfa z s −1 e z . En consecuencia, para cualquier complejo u , v ≠ 0 , se cumple, siempre que el camino de integración esté completamente contenido en el dominio de una rama del integrando. Si, además, la parte real de s es positiva, entonces se aplica el límite γ ( s , u ) → 0 para u → 0 , llegando finalmente a la definición integral compleja de γ [1]

Aquí es válido cualquier camino de integración que contenga sólo 0 en su inicio, o de lo contrario restringido al dominio de una rama del integrando, por ejemplo, la línea recta que une 0 y z .

Límite paraz → +∞
Valores reales

Dada la representación integral de una rama principal de γ , la siguiente ecuación se cumple para todos los números reales positivos s , x : [7]

scomplejo

Este resultado se extiende a los complejos s . Supongamos primero que 1 ≤ Re( s ) ≤ 2 y 1 < a < b . Entonces donde [8] se ha utilizado en el medio. Dado que la integral final se vuelve arbitrariamente pequeña si solo a es lo suficientemente grande, γ ( s , x ) converge uniformemente para x → ∞ en la franja 1 ≤ Re(s) ≤ 2 hacia una función holomorfa, [3] que debe ser Γ(s) debido al teorema de identidad. Tomando el límite en la relación de recurrencia γ ( s , x ) = ( s − 1) γ ( s − 1, x ) − x s − 1 e x y notando que lim x n e x = 0 para x → ∞ y todo n , se muestra que γ ( s , x ) converge también fuera de la franja hacia una función que obedece la relación de recurrencia de la función Γ. Se sigue que para todo complejo s no es un entero no positivo, x es real y γ es principal.

Convergencia sectorial

Ahora sea u del sector | arg z | < δ < π /2 con algún δ fijo ( α = 0 ), γ sea la rama principal en este sector, y observe

Como se muestra arriba, la primera diferencia puede hacerse arbitrariamente pequeña, si | u | es suficientemente grande. La segunda diferencia permite la siguiente estimación: donde hicimos uso de la representación integral de γ y la fórmula sobre | z s | anterior. Si integramos a lo largo del arco con radio R = | u | alrededor de 0 que conecta u y | u | , entonces la última integral es donde M = δ (cos δ ) −Re s e Im es una constante independiente de u o R . Nuevamente refiriéndonos al comportamiento de x n e x para x grande , vemos que la última expresión se acerca a 0 a medida que R aumenta hacia . En total ahora tenemos: si s no es un entero no negativo, 0 < ε < π /2 es arbitrariamente pequeño, pero fijo, y γ denota la rama principal en este dominio.

Descripción general

es:

Función gamma superior incompleta

En cuanto a la función gamma incompleta superior , una extensión holomorfa , con respecto a z o s , está dada por [1] en los puntos ( s , z ) , donde existe el lado derecho. Dado que es multivaluado, lo mismo se aplica a , pero una restricción a valores principales solo produce la rama principal univaluada de .

Cuando s es un entero no positivo en la ecuación anterior, ninguna parte de la diferencia está definida y un proceso de limitación , desarrollado aquí para s → 0 , completa los valores faltantes. El análisis complejo garantiza la holomorficidad , porque resulta estar acotado en un entorno de ese límite para un z fijo .

Para determinar el límite, es útil la serie de potencias de en z = 0. Al reemplazar por su serie de potencias en la definición integral de , se obtiene (supongamos x , s reales positivos por ahora): o [4] que, como representación en serie de la función completa, converge para todo complejo x (y todo complejo s que no sea un entero no positivo).

Al eliminarse la restricción a los valores reales, la serie permite la expansión:

Cuando s → 0 : [9] ( aquí es la constante de Euler-Mascheroni ), por lo tanto, es la función límite de la función gamma incompleta superior cuando s → 0 , también conocida como integral exponencial . [10]

Por medio de la relación de recurrencia, los valores de para números enteros positivos n pueden derivarse de este resultado, [11] por lo que la función gamma incompleta superior demuestra existir y ser holomorfa, con respecto tanto a z como a s , para todos los s y z ≠ 0 .

es:

Valores especiales

Aquí, es la integral exponencial , es la integral exponencial generalizada , es la función de error , y es la función de error complementaria , .

Comportamiento asintótico

Fórmulas de evaluación

La función gamma inferior se puede evaluar utilizando la expansión de la serie de potencias: [15] donde es el símbolo de Pochhammer .

Una expansión alternativa es donde M es la función hipergeométrica confluente de Kummer .

Conexión con la función hipergeométrica confluente de Kummer

Cuando la parte real de z es positiva, donde tiene un radio de convergencia infinito.

Nuevamente con funciones hipergeométricas confluentes y empleando la identidad de Kummer,

Para el cálculo real de valores numéricos, la fracción continua de Gauss proporciona una expansión útil:

Esta fracción continua converge para todo complejo z , siempre que s no sea un entero negativo.

La función gamma superior tiene la fracción continua [16] y [ cita requerida ]

Teorema de multiplicación

El siguiente teorema de multiplicación es válido:

Implementación de software

Las funciones gamma incompletas están disponibles en varios de los sistemas de álgebra computacional .

Sin embargo, incluso si no están disponibles directamente, los valores de funciones incompletas se pueden calcular utilizando funciones que se incluyen comúnmente en hojas de cálculo (y paquetes de álgebra computacional). En Excel , por ejemplo, se pueden calcular utilizando la función gamma combinada con la función de distribución gamma .

Estos se desprenden de la definición de la función de distribución acumulativa de la distribución gamma .

En Python , la biblioteca Scipy proporciona implementaciones de funciones gamma incompletas en scipy.special, sin embargo, no admite valores negativos para el primer argumento. La función gammaincde la biblioteca mpmath admite todos los argumentos complejos.

Funciones gamma regularizadas y variables aleatorias de Poisson

Dos funciones relacionadas son las funciones gamma regularizadas: es la función de distribución acumulativa para variables aleatorias gamma con parámetro de forma y parámetro de escala 1.

Cuando es un número entero, es la función de distribución acumulativa para las variables aleatorias de Poisson : Si es una variable aleatoria entonces

Esta fórmula se puede derivar mediante integración repetida por partes.

En el contexto de la distribución de conteo estable , el parámetro puede considerarse como el inverso del parámetro de estabilidad de Lévy : donde es una distribución de conteo estable estándar de forma .

y se implementan como [17] y [18] en scipy .gammaincgammaincc

Derivados

Usando la representación integral anterior, la derivada de la función gamma incompleta superior con respecto a x es La derivada con respecto a su primer argumento está dada por [19] y la segunda derivada por donde la función es un caso especial de la función G de Meijer Este caso especial particular tiene propiedades de cierre interno propias porque puede usarse para expresar todas las derivadas sucesivas. En general, donde es la permutación definida por el símbolo de Pochhammer : Todas estas derivadas pueden generarse en sucesión a partir de: y Esta función puede calcularse a partir de su representación en serie válida para , con el entendimiento de que s no es un entero negativo o cero. En tal caso, uno debe usar un límite. Los resultados para pueden obtenerse por continuación analítica . Algunos casos especiales de esta función pueden simplificarse. Por ejemplo, , , donde es la integral exponencial . Estas derivadas y la función proporcionan soluciones exactas a varias integrales mediante la diferenciación repetida de la definición integral de la función gamma incompleta superior. [20] [21] Por ejemplo, esta fórmula se puede inflar o generalizar aún más a una enorme clase de transformadas de Laplace y transformadas de Mellin . Cuando se combina con un sistema de álgebra computacional , la explotación de funciones especiales proporciona un método poderoso para resolver integrales definidas, en particular aquellas encontradas en aplicaciones prácticas de ingeniería (ver Integración simbólica para más detalles).

Integrales indefinidas y definidas

Las siguientes integrales indefinidas se obtienen fácilmente utilizando la integración por partes (omitiendo la constante de integración en ambos casos): La función gamma incompleta inferior y superior están conectadas a través de la transformada de Fourier : Esto se deduce, por ejemplo, de la especialización adecuada de (Gradshteyn et al. 2015, §7.642).

Notas

  1. ^ abcdef "DLMF: §8.2 Definiciones y propiedades básicas ‣ Funciones gamma incompletas ‣ Capítulo 8 Funciones gamma incompletas y relacionadas". dlmf.nist.gov .
  2. ^ abc "DLMF: §8.8 Relaciones de recurrencia y derivadas ‣ Funciones gamma incompletas ‣ Capítulo 8 Funciones gamma incompletas y relacionadas". dlmf.nist.gov .
  3. ^ ab Donald E. Marshall (otoño de 2009). "Complex Analysis" (PDF) . Matemáticas 534 (material para estudiantes). Universidad de Washington. Teorema 3.9 en la pág. 56. Archivado desde el original (PDF) el 16 de mayo de 2011. Consultado el 23 de abril de 2011 .
  4. ^ ab "DLMF: §8.7 Expansiones de series ‣ Funciones gamma incompletas ‣ Capítulo 8 Gamma incompleta y funciones relacionadas". dlmf.nist.gov .
  5. ^ Paul Garrett. "Teorema de Hartog: la analiticidad separada implica articulación" (PDF) . cse.umn.edu . Consultado el 21 de diciembre de 2023 .
  6. ^ C. Teleman. "Superficies de Riemann" (PDF) . berkeley.edu . Consultado el 21 de diciembre de 2023 .
  7. ^ "DLMF: §5.2 Definiciones ‣ Propiedades ‣ Capítulo 5 Función Gamma". dlmf.nist.gov .
  8. ^ "DLMF: §4.4 Valores especiales y límites ‣ Logaritmo, exponencial, potencias ‣ Capítulo 4 Funciones elementales". dlmf.nist.gov .
  9. ^ ver última ecuación.
  10. ^ "DLMF: §8.4 Valores especiales ‣ Funciones gamma incompletas ‣ Capítulo 8 Gamma incompleta y funciones relacionadas". dlmf.nist.gov .
  11. ^ "DLMF: 8.4 Valores especiales".
  12. ^ Weisstein, Eric W. "Función gamma incompleta". MathWorld .(ecuación 2)
  13. ^ ab Bender y Orszag (1978). Métodos matemáticos avanzados para científicos e ingenieros . Springer.
  14. ^ "DLMF: §8.11 Aproximaciones y expansiones asintóticas ‣ Funciones gamma incompletas ‣ Capítulo 8 Funciones gamma incompletas y relacionadas". dlmf.nist.gov .
  15. ^ "DLMF: §8.11 Aproximaciones y expansiones asintóticas ‣ Funciones gamma incompletas ‣ Capítulo 8 Funciones gamma incompletas y relacionadas". dlmf.nist.gov .
  16. ^ Abramowitz y Stegun p. 263, 6.5.31
  17. ^ "scipy.special.gammainc — Manual de SciPy v1.11.4". docs.scipy.org .
  18. ^ "scipy.special.gammaincc — Manual de SciPy v1.11.4". docs.scipy.org .
  19. ^ KO Geddes , ML Glasser, RA Moore y TC Scott, Evaluación de clases de integrales definidas que involucran funciones elementales mediante la diferenciación de funciones especiales , AAECC (Álgebra aplicable en ingeniería, comunicación y computación), vol. 1, (1990), págs. 149-165, [1]
  20. ^ Milgram, MS (1985). "La función integro-exponencial generalizada". Matemáticas. Comp . 44 (170): 443–458. doi : 10.1090/S0025-5718-1985-0777276-4 . MR  0777276.
  21. ^ Mathar (2009). "Evaluación numérica de la integral oscilatoria sobre exp(i*pi*x)*x^(1/x) entre 1 e infinito". arXiv : 0912.3844 [math.CA]., Aplicación B

Referencias

Enlaces externos