Errores en las medidas

Reglas para expresar una medida y su error

Cuando un físico mide algo debe tener gran cuidado para no producir una perturbación en el sistema que está bajo observación. Por ejemplo, cuando medimos la temperatura de un cuerpo, lo ponemos en contacto con un termómetro. Pero cuando los ponemos juntos, algo de energía o "calor" se intercambia entre el cuerpo y el termómetro, dando como resultado un pequeño cambio en la temperatura del cuerpo que deseamos medir. Así, el instrumento de medida afecta de algún modo a la cantidad que deseábamos medir

Además, todas las medidas está afectadas en algún grado por un error experimental debido a las imperfecciones inevitables del instrumento de medida, o las limitaciones impuestas por nuestros sentidos que deben de registrar la información.

Por ejemplo, al medir una cierta distancia hemos obtenido

297±2 mm.

De este modo, entendemos que la medida de dicha magnitud está en alguna parte entre 295 mm y 299 mm. En realidad, la expresión anterior no significa que se está seguro de que el valor verdadero esté entre los límites indicados, sino que hay cierta probabilidad de que esté ahí.

Una medida de una velocidad expresada de la forma

6051.78±30 m/s

es completamente ridícula, ya que la cifra de las centenas puede ser tan pequeña como 2 o tan grande como 8. Las cifras que vienen a continuación 1, 7 y 8 carecen de significado y deben de ser redondeadas. La expresión correcta es

6050±30 m/s

Una medida de 92.81 con un error de 0.3, se expresa

92.8±0.3

Con un error de 3, se expresa

93±3

Con un error de 30 se expresa

90±30

Los errores se deben dar solamente con una única cifra significativa. Únicamente, en casos excepcionales, se pueden dar una cifra y media (la segunda cifra 5 ó 0).

Expresiones incorrectas: el error se debe dar con una sola cifra significativa
24467±2928 m
23.463±0.165 cm
345.20±3.10 mm
Expresiones incorrectas: la medida y el error tienen el mismo orden de magnitud.
24467±3000 cm
43±0.06 m
345.2±3 m
Expresiones correctas
24000±3000 m
23.5±0.2 cm
345±3 m
43.00±0.06 m

Medidas directas

Un experimentador que haga la misma medida varias veces no obtendrá, en general, el mismo resultado, no sólo por causas imponderables como variaciones imprevistas de las condiciones de medida: temperatura, presión, humedad, etc., sino también, por las variaciones en las condiciones de observación del experimentador.

Si al tratar de determinar una magnitud por medida directa realizamos varias medidas con el fin de corregir los errores aleatorios, los resultados obtenidos son x1, x2, ... xn se adopta como mejor estimación del valor verdadero, el valor medio <x>, que viene dado por

<x>= x 1 + x 2 +... x n n = 1 n x i n

El valor medio, se aproximará tanto más al valor verdadero de la magnitud cuanto mayor sea el número de medidas, ya que los errores aleatorios de cada medida se va compensando unos con otros.

Sin embargo, en la práctica, no debe pasarse de un cierto número de medidas. En general, es suficiente con 10, e incluso podría bastar 4 ó 5.

Cuando la sensibilidad del método o de los aparatos utilizados es pequeña comparada con la magnitud de los errores aleatorios, puede ocurrir que la repetición de la medida nos lleve siempre al mismo resultado; en este caso, está claro que el valor medio coincidirá con el valor medido en una sola medida, y no se obtiene nada nuevo en la repetición de la medida y del cálculo del valor medio, por lo que solamente será necesario en este caso hacer una sola medida.

Los errores se producen por causas aleatorias, se toma como la mejor estimación del error, el llamado error cuadrático

Δx= 1 n ( x i <x>) 2 n(n1)

El resultado del experimento se expresa como

<x>±Δx y la unidad de medida

Es evidente, por ejemplo, tomando el caso más extremo, que si el resultado de las n medidas ha sido el mismo, el error cuadrático, de acuerdo con la formula será cero, pero eso no quiere decir que el error de la medida sea nulo. Sino, que el error instrumental es tan grande, que no permite observar diferencias entre las diferentes medidas y por tanto, el error instrumental, aquél que viene definido por la resolución del aparato de medida, será el error de la medida.

Ejemplos:

  1. Si al hacer una medida de la intensidad con un amperímetro cuya división o cifra significativa más pequeña es 0.01 A, la lectura es 0.64 A, y esta lectura es constante (no se observan variaciones al medir en diferentes instantes), tomaremos 0.64 como el valor de la medida y 0.01 A como su error. La medida se expresará así

  2. 0.64±0.01 A

  3. Supongamos que hemos medido un determinado tiempo, t, cuatro veces, y disponemos de un cronómetro que permite conocer hasta las décimas de segundo. Los resultados han sido: 6.3, 6.2, 6.4 y 6.2 s. De acuerdo a lo dicho anteriormente, tomaremos como valor medido el valor medio:

  4. <t>= 6.3+6.2+6.4+6.2 4 =6.275 s

    El error cuadrático será

    Δt= (6.36.275) 2 + (6.26.275) 2 + (6.46.275) 2 + (6.26.275) 2 43

    >> x=[6.3 6.2 6.4 6.2];
    >> media=mean(x)
    media =    6.2750
    >> error=std(x)/sqrt(length(x))
    error =    0.0479

    Este error se expresa con una sola cifra significativa, Δt=0.05 s. Pero el error cuadrático es menor que el error instrumental, que es 0.1 s, por lo que debemos tomar este último como el error de la medida y redondear en consecuencia el valor medio, por lo que el resultado final de la medida es

    t=6.3±0.1 s

  5. Consideremos un ejemplo similar al anterior, pero en que los valores obtenidos para el tiempo están más dispersos: 5.5, 5.7, 6.2 y 6.5 s. Se encuentra que el valor medio es 5.975, y el error cuadrático 0.2286737.

  6. >> x=[5.5 5.7 6.2 6.5];
    >> media=mean(x)
    media =    5.9750
    >> error=std(x)/sqrt(length(x))
    error =    0.2287

    El error cuadrático es en esta caso mayor que el error instrumental, por lo que debemos tomarlo como el error de la medida. Lo debemos redondear a 0.2 (una sola cifra significativa). Expresamos la medida finalmente como

    t=6.0±0.2 s

Medidas indirectas

En muchos casos, el valor experimental de una magnitud se obtiene, de acuerdo a una determinada expresión matemática, a partir de la medida de otras magnitudes de las que depende.

Funciones de una sola variable

  1. Se desea calcular el índice de refracción n de un vidrio midiendo el ángulo crítico θ, mediante la fórmula n=1/sinθ. Si medimos el ángulo θ es fácil calcular el índice de refracción n. Pero si conocemos el error de la medida del ángulo, necesitamos conocer el error del índice de refracción.

  2. Sea una función y=y(x). Como se aprecia en la figura, si el error Δx es pequeño. El error Δy se calcula del siguiente modo

    Δy=tanθ·Δx

    Pero tanθ es la pendiente de la recta tangente a al curva en el punto de abscisa x

    Como la pendiente puede ser positiva, si la función es creciente o negativa si la función es decreciente, en general tendremos que

    Δy=| f'(<x>) |Δx

    Sea y=cos x

    Sea x=20±3 º,

    y=cos20=0.9397

    El error Δx=3π/180=0.0524 rad

    Δy=|sin20|·0.0524=0.018

    y=0.94±0.02

  3. Supongamos que queremos medir el periodo P de un oscilador, es decir, el tiempo que tarda en efectuar una oscilación completa, disponemos de un cronómetro que aprecia las décimas de segundo, 0.1 s. Medimos el tiempo que tarda en hacer 10 oscilaciones, por ejemplo 4.6 s, dividiendo este tiempo entre 10 resulta P=0.46 s, que es el periodo "medio".

  4. P= t 10 ΔP= Δt 10

    Obtenemos para el error ΔP=0.01 s. Expresamos la medida como

    P=0.46±0.01 s

    Es evidente, que podemos aumentar indefinidamente la resolución instrumental para medir P, aumentando el número de periodos que incluimos en la medida directa de t. El límite está en nuestra paciencia y la creciente probabilidad de cometer errores cuando contamos el número de oscilaciones. Por otra parte, el oscilador no se mantiene con la misma amplitud indefinidamente, sino que se para al cabo de un cierto tiempo.

Función de varias variables

La magnitud y viene determinada por la medida de varias magnitudes p, q, r, etc., con la que está ligada por la función

y=f(p, q, r ...).

El error de la magnitud y viene dado por la siguiente expresión.

Δy= ( < f p >Δp ) 2 + ( < f q >Δq ) 2 + ( < f r >Δr ) 2 +...

Casos más frecuentes

z=x+yΔz= Δ x 2 +Δ y 2 z=xyΔz= Δ x 2 +Δ y 2 z=xy Δzz= ( Δx x ) 2 + ( Δy y ) 2 z= x y Δzz= ( Δx x ) 2 + ( Δy y ) 2

Ejemplo 1: La medida de los lados de un rectángulo son 1.53±0.06 cm, y 10.2±0.1 cm, respectivamente. Hallar el área del rectángulo y el error de la medida indirecta.

El área es z=1.53×10.2=15.606 cm2

>> x=1.53; dx=0.06;
>> y=10.2; dy=0.1;
>> z=x*y
z =   15.6060
>> dz=z*sqrt((dx/x)^2+(dy/y)^2)
dz =    0.6308

El error en el área Δz/z se obtiene aplicando la fórmula del producto de dos magnitudes.

Δz z = ( 0.06 1.53 ) 2 + ( 0.1 10.2 ) 2 =0.0404422504 Δz=( 1.5310.2 )0.0404422504=0.63083

El error con una sola cifra significativa es 0.6. La medida del área junto con el error y la unidad se escribirá correctamente

z=15.6±0.6 cm2

Ejemplo 2: Queremos calcular la aceleración de la gravedad g, midiendo el periodo P de las oscilaciones de un péndulo de longitud l

P=2π l g g=4 π 2 l P 2

La expresión del error Δg de la variable dependiente g es

Δg= ( 4 π 2 1 P 2 Δl ) 2 + ( 4 π 2 2l P 3 ΔP ) 2 =4 π 2 l P 2 ( Δl l ) 2 + ( 2ΔP P ) 2 Δg g = ( Δl l ) 2 + ( 2ΔP P ) 2

Supongamos que medimos el periodo P y la longitud l del péndulo

P=1.95±0.01 s
l
=94±1 cm

Calculamos la aceleración de la gravedad y el error

g=9.759293 m/s2
Δg=0.144

Expresamos correctamente la medida y el error de g

g=9.8±0.1 m/s2

Utilizando el Symbolic Math Toolbox de MATLAB

>> syms l P dl dP;
>> g=4*sym('pi')^2*l/P^2;
>> gg=subs(g,{P,l},{1.95,0.94})
gg =(1504*pi^2)/1521
>> double(gg)
ans =    9.7593

>> dg=sqrt((diff(g,l)*dl)^2+(diff(g,P)*dP)^2);
>> simplify(dg)
ans =4*pi^2*((P^2*dl^2 + 4*dP^2*l^2)/P^6)^(1/2)
>> dgg=subs(dg,{P,l,dP,dl},{1.95,0.94,0.01,0.01}) 
dgg =4*((1173904*pi^4)/87968594025)^(1/2)
>> double(dgg)
ans =    0.1442	 

Ejemplo 3: Queremos calcular el índice de refracción de un vidrio, midiendo el ángulo que hace el rayo incidente y el rayo refractado con la normal. Aplicamos la ley de Snell de la refracción

n= sini sinr

Cálculo del error en la medida del índice de refracción n.

Δn n = ( 1 tani Δi ) 2 + ( 1 tanr Δr ) 2

Sea i=20±1º y r=13±1º.
1º=π/180 rad

Se calcula el índice de refracción n y el error Δn

>> syms i r di dr;
>> n=sin(i)/sin(r);
>> nn=subs(n,{i,r},{20*pi/180,13*pi/180});
>> double(nn)
ans =    1.5204

>> dn=sqrt((diff(n,i)*di)^2+(diff(n,r)*dr)^2);
>> simplify(dn)
ans =((di^2*cos(i)^2*sin(r)^2 + dr^2*cos(r)^2*sin(i)^2)/sin(r)^4)^(1/2)
>> dnn=subs(dn,{i,r,di,dr},{20*pi/180,13*pi/180,pi/180,pi/180});
>> double (dnn)
ans =    0.1361 

n=1.52
Δn=0.136

Expresamos correctamente la medida y el error de n

n=1.5±0.1

Referencias

Taylor J. R. An Introduction to Error Analysis. The Study of Uncertainties in Physical Measurements. University Science Books (1982)