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.
Propagació dels errors aleatoris
Descripció del bloc
Introducció «
La Fig. 1 es presentà en el tema sobre la propagació de la incertesa per funcions de n variables. El valor d'una magnitud es calculava amb una funció de quatre variables els valors de les quals s'havien determinat experimentalment i tenien determinades incerteses. Les quatre variables eren independents i la probabilitat que el valor real de totes elles fes augmentar el valor de la magnitud calculada havia de ser més petita que la probabilitat que almenys alguna els fes minvar.
Figura 1. Aquesta figura es presentà en un tema anterior a l'apartat «Exemple amb la mesura de la relació u = e/m».
Quan es dóna un resultat en la forma x = x0 ± δx s'ha de saber quin tipus d'error és δx perquè la interpretació del resultat canvia. En alguns casos, δx pot ser una cota absoluta dels valors de x. En tal cas, s'està segur que el valor de x està dins l'interval [x0 – δx, x0 + δx]. Però amb errors aleatoris, δx és habitualment l'error estàndard sx (tal vegada corregit pel factor de Student). En tal cas, només s'està segur que el resultat estarà dins l'interval [x0 – δx, x0 + δx] amb una probabilitat del 68%.
Fer una mesura d'una magnitud en presència de diverses fonts d'errors aleatoris és com mirar el valor que pren a l'atzar una variable aleatòria gaussiana. La mitjana i la variància de la funció densitat de probabilitat d'aquesta variable són desconegudes, però es poden estimar amb la mitjana i l'error estàndard d'un conjunt de mesures de la magnitud considerada.
Amb errors aleatoris, escriure x1 = m1 ± s1 és com dir que hi ha una funció densitat de probabilitat gaussiana Gμ,σ(x) de mitjana μ = m1 i desviació estàndard s1 que es pot usar per determinar la probabilitat que el valor real de x estigui dins un determinat interval.
La mateixa interpretació es pot fer per a una altra magnitud x2 = m2 ± s2. Llavors resulta immediat plantejar la qüestió de quina funció densitat de probabilitat tendran les magnituds calculades a partir de x1 i x2. Les operacions bàsiques considerades explícitament en aquest bloc de temes seran multiplicar una magnitud per una constant (c x1), sumar dues magnituds (x1 + x2) i dividir-les (x1 / x2). En general, s'usaran funcions d'una variable u(x1) o de dues, v(x1, x2).
Els casos c x1 i x1 + x2 tenen una transcendència més gran de la que pot semblar atesa la simplicitat de les expressions. Si la funció densitat de probabilitat de les variables són gaussianes, Gμ1,s1(x) i Gμ2,s2(x), la funció densitat de probabilitat dels valors c x1 també és gaussiana, amb mitjana c μ1 i variància c σ1. La funció densitat de probabilitat dels valors x1 + x2 es dóna a la Taula 1. Els resultats d'aquests dos casos són el punt de partida per analitzar valors calculats amb funcions d'una, dues o més variables.
Quan les incerteses són petites, les funcions es poden aproximar per sèries de Taylor de primer ordre. Així, els valors de les funcions es poden calcular amb expressions dels tipus
u(x1, x2) = c0 + c1 x1 + c2 x2,
on c0 és una constant (que no canvia la variància) i c1 i c2 són les derivades parcials de la funció avaluades a x1 = μ1 i x2 = μ2. La multiplicació per una constant o la suma amb altres variables gaussianes no canvia el fet que la funció densitat de probabilitat dels valors obtinguts sigui gaussiana. La Taula 1, recull les expressions de les mitjanes i la variàncies que es demostraran en els temes següents d'aquest bloc. En el darrer tema, la funció densitat de probabilitat dels valors x1 / x2, que es pot calcular de manera exacta de manera no massa difícil, es compararà amb la gaussiana calculada de manera aproximada desenvolupant x1 / x2 en sèrie de Taylor a primer ordre.
variable | mitjana | variància |
Taula 1. Relació de les mitjanes i variàncies de les variables indicades amb les mitjanes i variàncies de les variables x1 i x2. El subíndex μ a les derivades indica que s'han d'avaluar a x1 = μ1, x2 = μ2.