Función gamma incompleta

De Wikipedia, la enciclopedia libre
(Redirigido desde «Funcion gamma incompleta»)
Saltar a: navegación, búsqueda

En matemática, la función gamma se define como una integral definida. La función gamma incompleta se define como una integral definida del mismo integrando.

Hay dos tipos de función gamma incompleta, una para el caso en el que varía el límite inferior de integración, y otro cuando varía el límite superior. La primera se denota como \Gamma(a,x) y se define como

 \Gamma(a,x) = \int_x^{\infty} t^{a-1}\,e^{-t}\,dt .\,\!

La segunda se escribe \gamma(a,x) y se define como

 \gamma(a,x) = \int_0^x t^{a-1}\,e^{-t}\,dt .\,\!

En ambos casos, x es una variable real mayor o igual que cero, y a es una variable compleja, cuya parte real es positiva.

Propiedades[editar]

Integrando por partes se demuestra que

 \Gamma(a+1,x) = a\Gamma(a,x) + x^a e^{-x}\,
 \gamma(a+1,x) = a\gamma(a,x) - x^a e^{-x}.\,

Dado que la función gamma ordinaria se define como

 \Gamma(a) = \int_0^{\infty} t^{a-1}\,e^{-t}\,dt \,\!

tenemos que

 \gamma(a,x) + \Gamma(a,x) = \Gamma(a).\,

Además,

\Gamma(a,x)=(a-1)!e^{-x}\sum_{k=0}^{a-1}\frac{x^k}{k!} si a es un número entero. (Weisstein)
 \Gamma(a,0) = \Gamma(a)\,
 \Gamma(a) = (a-1)!\, si a es un número entero.

y

 \gamma(a,x) \rightarrow \Gamma(a) 
  \quad \mathrm{as\ } x \rightarrow \infty.  \,

También,

\Gamma(0,x) = -\mbox{Ei}(-x)\mbox{ for }x>0 \,
\Gamma\left({1 \over 2}, x\right) = \sqrt\pi\,\mbox{erfc}\left(\sqrt x\right) \,
\gamma\left({1 \over 2}, x\right) = \sqrt\pi\,\mbox{erf}\left(\sqrt x\right) \,
\Gamma(1,x) = e^{-x} \,
\gamma(1,x) = 1 - e^{-x} \,

donde Ei es la función integral exponencial, erf es la función de error, y erfc la función de error complementaria, erfc(x) = 1 − erf(x).

Funciones Gamma regularizadas[editar]

Dos funciones relacionadas son las funciones Gamma regularizadas:

P(a,x)=\frac{\gamma(a,x)}{\Gamma(a)}
Q(a,x)=\frac{\Gamma(a,x)}{\Gamma(a)}=1-P(a,x)

Derivadas[editar]

La derivada de la función gamma incompleta  \Gamma (a,x) en x es bien conocida. Es dado simplemente por el integrando de su definición completa:

 
\frac{\partial \Gamma (a,x) }{\partial x} = - x^{a-1} e^{-x}

La derivada con respecto a la parámetro a viene dada por[1]

 \frac{\partial \Gamma (a,x) }{\partial a} = \ln (x) \Gamma (a,x) + x ~T(3,a,x)

y la segunda derivada es:

 \frac{\partial^2 \Gamma (a,x) }{\partial a^2} = \ln^2 (x) \Gamma (a,x) + 2 x ~ ( \ln (x) ~ T(3,a,x) + T(4,a,x) )

donde la función "T(m,a,x)" es un caso especial de Meijer función G

 T(m,a,z) = G_{m-1, m}^{~m,~0} \left( x  \left|  \begin{array}{c} 0,0, \ldots 0 \\ -1, -1, \ldots, a-1, -1 \end{array} \right. \right) ~.

Este caso tiene propiedades internas de cierre de los suyos, dado que puede expresar todas las derivadas sucesivas. En general,

 \frac{\partial^m \Gamma (a,x) }{\partial a^m} = \ln^m (x) \Gamma (a,x) + m x ~ \sum_{i=0}^{m-1} P_i^{m-1} \ln^{m-i-1} (x) ~ T(3+i,a,x)

dóndeP_j^i es la permutación definida por el símbolo de Pochhammer:

P_j^i = \left( \begin{array}{l} i \\ j \end{array} \right) j! = \frac{i!}{(i-j)!} ~ .

Todos estos derivados pueden ser producidos a partir de:

 \frac{\partial T (m,a,x) }{\partial a} = \ln (x) ~ T(m,a,x) + (m-1) T(m+1,a,x)

y

 \frac{\partial T (m,a,x) }{\partial x} = -\frac{1}{x} (T(m-1,a,x) + T(m,a,x))

Esta función T(m,a,x) se puede calcular por su representación estándar, siempre que : |z| < 1 ,

 T(m,a,z) = - \frac{(-1)^{m-1} }{(m-2)! } \frac{d^{m-2} }{dt^{m-2} } \left. (\Gamma (a-t) z^{t-1} ) \right]_{t=0} + \sum_{i=0}^{\infty} \frac{(-1)^i z^{a-1+i}}{i! (-a-i)^{m-1} }

y siempre que el parámetro a no es un número entero negativo o cero. En este último caso, se debe utilizar un límite. Resultados de  |z| \ge 1 Se puede obtener por una extensión analítica. Algunos casos especiales de esta función se puede simplificar. Por ejemplo,

 T(2,a,x) = \frac{\Gamma(a,x)}{x}
  x ~ T(3,1,x) = E_1 (x)

donde  E_1 (x) es la función integral exponenciel. Los derivados y la función T(m,a,x) proporcionar soluciones exactas a un número de integrales por la diferenciación repetido de la definición completa de la función gamma incompleta  \Gamma (a,x) . Por ejemplo,

 
\int_{x}^{\infty} t^{a-1} \ln^m (t) ~ e^{-t} = \frac{\partial^m}{\partial a^m} \int_{x}^{\infty} t^{a-1} e^{-t} = \frac{\partial^m}{\partial a^m} \Gamma (a,x)

Esta fórmula se puede "inflar" o más generalizado a una gran clase de la transformada de Laplace o de Mellin. Una vez combinado con un sistema algebraico computacional, funcionamiento de las funciones especiales proporciona un método poderoso para resolver integrales definidas, en particular los que se enfrentan los ingenieros de aplicaciones prácticas. Este método fue inventado por el sistema Maple[2] y, más tarde imitada por Mathematica, MuPAD y otros sistemas. La función "T(m,a,x)" era conocido en el grupo de investigación de arce como una función de Scott-G.

Notas[editar]

  1. K.O Geddes, M.L. Glasser, R.A. Moore y T.C. Scott, Evaluation of Classes of Definite Integrals Involving Elementary Functions via Differentiation of Special Functions, AAECC (Applicable Algebra in Engineering, Communication and Computing), vol. 1, (1990), pp.149-165, [1]
  2. K.O. Geddes y T.C. Scott, Recipes for Classes of Definite Integrals Involving Exponentials and Logarithms, Proceedings of the 1989 Computers and Mathematics conference, (MIT 12 de junio, 1989), editado por E. Kaltofen y S.M. Watt, Springer-Verlag, New York, (1989), pp. 192-201. [2]

Enlaces externos[editar]

Referencias[editar]

  • G. Arfken y H. Weber. Mathematical Methods for Physicists. Harcourt/Academic Press, 2000. (Ver capítulo 10.)
  • W. H. Press, B. P. Flannery, S. A. Teukolsky, y W. T. Vetterling. Numerical Recipes in C. Cambridge, UK: Cambridge University Press, 1988. (Ver Sección 6.2.)