Anàlisi de Dades Experimentals
© Antoni Amengual Colom. Departament de Física, Universitat de les Illes Balears.

Versió 1.0 publicada al setembre de 2013. DL: PM 860-2013.

Apèndix

Funció gamma

Introducció «

La funció gamma, que s'identifica amb la lletra grega d'aquest nom en majúscula, es pot definir de maneres diferents. Aquí es consideraran dues definicions: el límit d'Euler i la integral d'Euler. La primera és vàlida per a variables reals i complexes. La segona només és vàlida si la part real de la variable complexa és positiva.

Dues propietats fonamentals de la funció gamma són el valor que té per a z = 1,

Γ(1) = 1,

i la relació de recurrència

Γ(z + 1) = z Γ(z).

Aplicant aquesta darrera relació de manera repetida després de prendre z igual a un nombre sencer positiu  n fins que surti Γ(1) a la dreta, es troba

Γ(n + 1) = n!

Per aquesta igualtat, la funció gamma es considera una extensió del factorial d'un número sencer positiu als números reals i complexos.

El límit d'Euler «

El límit d'Euler és

"ade_1.gif"

Com resulta evident mirant el denominador, la funció gamma divergeix a z = 0 i a z igual a qualsevol sencer negatiu. A la Fig. 1 es mostra una representació de la funció gamma per a z igual a un número real x.

"ade_2.gif"

Figura 1. Representació de la funció gamma per a x real. Els punts taronja marquen el valor de la funció gamma per a x igual a 1, 2, 3 i 4. El mínim a la part positiva de les abscisses és Γ(1.46) = 0.886.

De manera particular i directa a partir de la definició resulta

"ade_3.gif"

Per demostrar que es compleix la relació de recurrència, s'aplica la definició a la variable z + 1:

"ade_4.gif"

El darrer terme de color verd es posa com nz per n, llavors es multiplica i divideix per z i s'agrupen termes. Així queda

"ade_5.gif"

Finalment, s'usa la relació

"ade_6.gif"

per obtenir

"ade_7.gif"

El primer terme és Γ(z), el segon val z i el producte dóna la relació de recurrència.

La integral d'Euler «

La funció gamma per a variables complexes amb la part real positiva es pot definir amb la integral d'Euler

"ade_8.gif"

De manera particular i directa a partir de la definició resulta

"ade_9.gif"

Per demostrar que es compleix la relació de recurrència, s'aplica la definició a la variable z + 1:

"ade_10.gif"

i es fa una integració per parts,

"ade_11.gif"

Per demostrar l'equivalència de la integral d'Euler en el domini d'existència amb el límit d'Euler es comença per definir

"ade_12.gif"

Es mostrarà a continuació que segons com es calculi

"ade_13.gif"

s'obté la integral o el límit d'Euler i com que el límit és únic, l'equivalència entre les dues definicions quedarà demostrada.

Primer càlcul del límit «

Per calcular el límit es té en compte que el límit de la integral és igual a la integral del límit (aquest fet s'hauria de demostrar) i la definició del número e,

"ade_14.gif"

Segon càlcul del límit «

A la integral que defineix Γ(z, n) es fa el canvi de variable

t = n u

i s'integra succesivament per parts. Amb el canvi de variable es té

"ade_15.gif"

S'integra per parts. La part integrada s'anul·la en els límit 0 i 1 i queda

"ade_16.gif"

Es repeteix la integració per parts i s'obté

"ade_17.gif"

La repetició successiva de la integració per parts n vegades dóna

"ade_18.gif"

La integral ara és immediata i s'obté finalment

"ade_19.gif"

La definició de Γ(z) era el límit per a n → ∞ d'aquesta expressió,

"ade_20.gif"

Queda així demostrat que el límit dóna Γ(z) en la forma de límit d'Euler com aquí o d'integral d'Euler com a l'apartat anterior. Per tant, les definicions són equivalents en el domini de z tal que la part real de z és positiva.

Valors de Γ(n/2) per a n sencer «

S'ha de distingir si n és parell o imparell.

n parell «

Els valors de la funció gamma per als semienters n/2 es calculen directament amb el factorial si n es parell. Escrivint n = 2k, amb k un número natural, es té

Γ(n/2) = Γ(k) = (k – 1)!

n imparell «

S'ha de partir de

"ade_21.gif"

Llavors, es pot aplicar la relació de recurrència per calcular els valors de la funció gamma per a semienters més grans. Per exemple,

"ade_22.gif"

A la Taula 1 es recullen els valors de Γ(n/2) per a n imparells des de 1 fins a 15.

n13579111315
Γ(n2)ππ23π415π8105π16945π3210395π64135135π128

Taula 1. Valors de la funció gamma per als primers valors de n / 2 amb n imparell.