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 -2013.

Distribució normal (gaussiana)

Integral de la funció densitat de probabilitat

Introducció «

Es presenta la integral de la funció densitat de probabilitat gaussiana

"ade_1.gif"

i la taula de valors d'aquesta integral. Els fulls de càlcul i altres programes informàtics tenen funcions per calcular la integral per a qualsevol valor de x0. Les taules es presenten per plantejar l'exercici de trobar la incertesa d'un valor obtingut per interpolació.

La integral en un interval finit a partir de la integral des de menys infinit a x0 «

La probabilitat que una variable aleatòria x amb funció de distribució de probabilitat gaussiana prengui un valor dins l'interval [xmin, xmàx] és

"ade_2.gif"

amb la funció

"ade_3.gif"

Si es té present que una integral dóna l'àrea sota la corba, resultarà obvi que la integral que dóna la probabilitat P es pot calcula amb la resta (Fig. 1)

"ade_4.gif"

Aquesta separació revela la utilitat de definir la funció integral Iμ,σ de la funció gaussiana:

"ade_5.gif"

"ade_6.gif"

Figura 1. L'àrea de la zona taronja dóna la probabilitat d'obtenir x entre xmin i xmàx.

Les funcions trigonomètriques d'un angle o el logaritme d'un número es consultaven en taules fins que les calculadores les feren innecessàries.

El valor de la integral Iμ,σ es pot obtenir amb un full de càlcul. La instrucció per a excel o calc és DISTR.NORM(x0; μ; σ; c) amb c = 1 (cert). La funció dóna el valor de la gaussiana en lloc de la integral si es posa c = 0.

No és l'ús de la taula el que resulta interessant aquí sinó els exercicis que es resoldran usant-la.

Relació amb la integral de la gaussiana centrada a l'origen i variància unitat «

Les integrals d'una funció gaussiana es poden relacionar amb la integral de la gaussiana centrada a l'origen i amb desviació estàndard unitat, G0,1(x), amb el canvi de variable

"ade_7.gif"

Es tendrà dx = σ dz i el límit superior de la integral passarà a ser

"ade_8.gif"

de manera que

"ade_9.gif"

En el pas (1) s'ha fet el canvi de variable i, en el pas (2), s'ha identificat la funció amb la gaussiana de mitjana 0 i variància 1. Resulta així que

"ade_10.gif"

i, en conseqüència, coneixent

"ade_11.gif"

es podran calcular les probabilitats per a tot tipus de gaussianes. Els valors d'aquesta integral per a un conjunt de valors z0 = a + b es donen a la taula presentada a l'apartat següent.

Taula dels valors de la integral de la gaussiana de mitjana 0 i variància unitat «

A la Taula 1 es presenten els valors de la integral

"ade_12.gif"

Només s'inclouen valors de z0 = a + b positius perquè la gaussiana és simètrica i està normalitzada, de manera que

"ade_13.gif"

Aquesta relació s'il·lustra a la Fig. 2 i es basa en que l'obtenció de z > z0 és el succés complementari d'obtenir z ≤ z0.

"ade_14.gif" "ade_15.gif"

Figura 2. L'àrea total sota la corba és 1, de manera que l'àrea verda de l'esquerra és igual a 1 menys l'àrea vermella de la dreta la qual s'obté amb la funció I(z).

a\b0.0.010.020.030.040.050.060.070.080.09
0.00.500000.503990.507980.511970.515950.519940.523920.527900.531880.53586
0.10.539830.543800.547760.551720.555670.559620.563560.567490.571420.57535
0.20.579260.583170.587060.590950.594830.598710.602570.606420.610260.61409
0.30.617910.621720.625520.629300.633070.636830.640580.644310.648030.65173
0.40.655420.659100.662760.666400.670030.673640.677240.680820.684390.68793
0.50.691460.694970.698470.701940.705400.708840.712260.715660.719040.72240
0.60.725750.729070.732370.735650.738910.742150.745370.748570.751750.75490
0.70.758040.761150.764240.767300.770350.773370.776370.779350.782300.78524
0.80.788140.791030.793890.796730.799550.802340.805110.807850.810570.81327
0.90.815940.818590.821210.823810.826390.828940.831470.833980.836460.83891
1.00.841340.843750.846140.848490.850830.853140.855430.857690.859930.86214
1.10.864330.866500.868640.870760.872860.874930.876980.879000.881000.88298
1.20.884930.886860.888770.890650.892510.894350.896170.897960.899730.90147
1.30.903200.904900.906580.908240.909880.911490.913090.914660.916210.91774
1.40.919240.920730.922200.923640.925070.926470.927850.929220.930560.93189
1.50.933190.934480.935740.936990.938220.939430.940620.941790.942950.94408
1.60.945200.946300.947380.948450.949500.950530.951540.952540.953520.95449
1.70.955430.956370.957280.958180.959070.959940.960800.961640.962460.96327
1.80.964070.964850.965620.966380.967120.967840.968560.969260.969950.97062
1.90.971280.971930.972570.973200.973810.974410.975000.975580.976150.97670
2.00.977250.977780.978310.978820.979320.979820.980300.980770.981240.98169
2.10.982140.982570.983000.983410.983820.984220.984610.985000.985370.98574
2.20.986100.986450.986790.987130.987450.987780.988090.988400.988700.98899
2.30.989280.989560.989830.990100.990360.990610.990860.991110.991340.99158
2.40.991800.992020.992240.992450.992660.992860.993050.993240.993430.99361
2.50.993790.993960.994130.994300.994460.994610.994770.994920.995060.99520
2.60.995340.995470.995600.995730.995850.995980.996090.996210.996320.99643
2.70.996530.996640.996740.996830.996930.997020.997110.997200.997280.99736
2.80.997440.997520.997600.997670.997740.997810.997880.997950.998010.99807
2.90.998130.998190.998250.998310.998360.998410.998460.998510.998560.99861

Taula 1. Àrea sota la corba de la funció de distribució gaussiana de mitjana 0 i desviació estàndard 1 des de menys infinit fins a a + b.

Valors que s'haurien de recordar «

A la Fig. 3 es mostren sis àrees els valors de les quals s'haurien de recordar.

"ade_16.gif" "ade_17.gif"

"ade_18.gif" "ade_19.gif"

"ade_20.gif" "ade_21.gif"

Figura 3. L'àrea destacada en color verd representa el percentatge indicat de l'àrea total.

Exemples «

Es consideren a continuació tres casos típics que s'han de saber resoldre.

Ús directe de la taula «

La probabilitat que un succés aleatori amb probabilitat gaussiana de mitjana 0 i variància 1 doni un valor inferior a a+ b, amb

a = 0.7,

b = 0.02,

es pot obtenir directament de la taula amb cinc decimals,

"ade_22.gif"

a\b0.0.010.020.030.040.050.060.070.080.09
0.60.725750.729070.732370.735650.738910.742150.745370.748570.751750.75490
0.70.758040.761150.764240.767300.770350.773370.776370.779350.782300.78524
0.80.788140.791030.793890.796730.799550.802340.805110.807850.810570.81327

Taula 2. Valor de la taula corresponent a a + b = 0.72 destacat en vermell.

Càlcul per interpolació «

La probabilitat que un succés aleatori amb probabilitat gaussiana de mitjana 0 i variància 1 doni un valor més petit que un valor donat amb més dos decimals, s'ha de calcular interpolant entre dos valors de la taula.

En general, si x1 < x < x2 s'usarà

"ade_23.gif"

Per exemple, I0,1(0.724) es determinarà mirant els valors de I0,1(0.72) i I0,1(0.73) i suposant que hi ha una variació lineal del valor de la integral (Fig. 3)

"ade_24.gif"

"ade_25.gif"

Figura 4. Donat el valor de la integral per a dues abscisses, el valor en un punt intermedi es pot determinar aproximadament suposant que I(x) varia linealment.  

Incertesa del valor interpolat «

Els valors de la integral I0,1(x) de la taula estan arrodonits al cinquè decimal. El valor exacte sense arrodonir és el de la taula més o menys 0.000005.

L'arrodoniment a un nombre finit de decimals introdueix incertesa. El valor exacte podria ser qualsevol que arrodonit a cinc decimals donés el que surt. Sense més informació disponible, s'ha de suposar que

"ade_26.gif"

Usant els valors 0.764235 i 0.767305, la interpolació dóna 0.765465. Però els valors extrems es poden combinar de quatre maneres i es donen quatre possibles valors de la interpolació. Els valors mínim i màxims són 0.765459 i 0.765469. El valor central és

"ade_27.gif"

Com que el màxim i el mínim disten d'aquest valor central

"ade_28.gif"

el resultat es pot expressar en la forma

"ade_29.gif"

El valor exacte amb sis decimals és 0.765467. La diferència amb el valor central és 0.000003 < 0.000005.

En general es pot cercar la incertesa sobre la integral usant les fórmules de la propagació de l'error per la funció

"ade_30.gif"

Per simplificar la notació, s'escriurà

"ade_31.gif"

Es tendrà

"ade_32.gif"

Es calcularà ara només l'error degut als dos primers termes

"ade_33.gif"

Com que x està entre x1 i x2 es poden llevar els valors absoluts i, com que les incerteses sobre I1 i I2 són iguals, només queda

"ade_34.gif"

Integral d'una gaussiana de mitjana i desviació estàndard qualsevol «

La probabilitat que en un succés aleatori amb probabilitat gaussiana de mitjana μ i variància σ s'obtengui un valor x0 s'ha de determinar amb la taula usant la relació

"ade_35.gif"

Per calcular la probabilitat d'obtenir x < 5.04 amb la funció densitat de probabilitat gaussiana de mitjana 4.8 i desviació estàndard 1/3, primer s'ha de fer el canvi de variable

"ade_36.gif"

En general el valor de z0 no sortirà a la taula i s'haurà de fer una interpolació; però, en aquest cas, es pot usar la taula directament per obtenir

"ade_37.gif"

La igualtat entre els dos primers termes està il·lustrada a la Fig. 4.

"ade_38.gif"

Figura 5. L'àrea sota una gaussiana es pot relacionar amb l'àrea sota una gaussiana centrada en l'origen i de desviació estàndard unitat fent un canvi de variable. L'àrea de la zona blava de la figura és igual a l'àrea de la zona verda.