27.08.2013 Views

Propagacion de errores. Practica Nro1

Propagacion de errores. Practica Nro1

Propagacion de errores. Practica Nro1

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Propagación <strong>de</strong> <strong>errores</strong><br />

La formula más general <strong>de</strong> propagación <strong>de</strong> <strong>errores</strong> para una función z(x,y,t) que <strong>de</strong>pen<strong>de</strong><br />

<strong>de</strong> las variables x, y y t; las cuáles tienen un correspondiente error Δx, Δy, Δt es:<br />

∂z ∂z ∂z<br />

Δz≤ . Δ x+ . Δ y+<br />

. Δt<br />

(1)<br />

∂x ∂y ∂t<br />

Es <strong>de</strong>cir, el error va a ser, a lo sumo, como el valor <strong>de</strong>l término <strong>de</strong> la <strong>de</strong>recha. En este<br />

caso se supusieron tres variables x, y, t pero podría haber más o menos. Las barras<br />

indican valor absoluto, es <strong>de</strong>cir, la contribución al error es siempre positiva.<br />

Las <strong>de</strong>rivadas parciales indican que se <strong>de</strong>be <strong>de</strong>rivar respecto a una variable <strong>de</strong>jando el<br />

resto como si fueran constantes (<strong>de</strong>rivar respecto a x suponiendo y y t constantes, y así<br />

sucesivamente).<br />

Ejemplo<br />

Sea la función q(x,y)=x 2 ·y-x·y 2 , don<strong>de</strong> x e y son dos variables <strong>de</strong>terminadas<br />

experimentalmente, con su correspondiente error.<br />

Supongamos que se tiene: x=3.0 ±0.1 y=2.0±0.1.<br />

¿Cuál es el error <strong>de</strong> q?<br />

∂q ∂q<br />

2 2<br />

Δ q= . Δ x+ . Δ y = 2 xy− y . Δ x+ x −2xy . Δy<br />

∂x ∂y<br />

Ahora se <strong>de</strong>ben evaluar los valores<br />

2 2<br />

Δ q = 2 . 3 . 2 − 2 . 0.1+ 3 − 2 . 3 . 3 . 0.1 = 8 . 0.1 + 3 . 0.1=1.1<br />

Si consi<strong>de</strong>ramos una sola cifra significativa entonces Δq=1<br />

Calculando el valor <strong>de</strong> q se obtiene q = 18 – 12 = 6 por lo tanto se <strong>de</strong>bería reportar:<br />

q = 6 ± 1<br />

Un ejemplo más complejo.<br />

Supongamos que se <strong>de</strong>be <strong>de</strong>spejar v0 <strong>de</strong> las dos ecuaciones <strong>de</strong> movimiento,<br />

x = x + v θ t (2)<br />

f<br />

0 0 cos<br />

1<br />

yfy v sen t g t<br />

2<br />

2<br />

= 0 + 0 θ − (3)<br />

tal como se requería en la práctica <strong>de</strong> laboratorio. Entonces una manera <strong>de</strong> hacerlo es<br />

<strong>de</strong>spejar t <strong>de</strong> la ecuación en x:<br />

x f − x0<br />

t = (4)<br />

v θ<br />

0 cos<br />

e introducirlo en la ecuación correspondiente en el eje y:<br />

2 2<br />

x f −x0 1 ⎛ xf −x0 ⎞ 1 ⎛ xf −x0<br />

⎞<br />

yf= y0 + v0senθ − g⎜ ⎟ = y0 + tgθ( xf−x0) − g⎜<br />

v0cosθ 2 ⎝v0cosθ ⎠ 2 ⎝v0cosθ ⎠ ⎟ (5)


Ahora, <strong>de</strong>spejando v0<br />

v<br />

0<br />

=<br />

g<br />

xf−x0 2 cos θ − y + y + tgθ( x −x<br />

)<br />

f 0 f 0<br />

(6)<br />

Ahora bien, para calcular el error <strong>de</strong> v0, se <strong>de</strong>bería aplicar la ecuación (1); es <strong>de</strong>cir,<br />

hacer la propagación <strong>de</strong> <strong>errores</strong> <strong>de</strong> todas las variables que tienen error. En este caso<br />

vamos a consi<strong>de</strong>rar que las variables que aportan al error son xf, y0 y θ. Por simplicidad<br />

vamos a asumir x0 e yf iguales a cero, <strong>de</strong> modo que<br />

v<br />

0<br />

=<br />

Entonces:<br />

g x f 1<br />

2 cosθ<br />

y + x ⋅ tgθ<br />

0<br />

0 0 0<br />

0<br />

∂xf f<br />

∂y0 0<br />

∂<br />

f<br />

(7)<br />

∂v ∂v ∂v<br />

Δv ≤ . Δ x + . Δ y + . Δθ<br />

θ<br />

Vayamos por parte, primero calculemos las <strong>de</strong>rivadas parciales.<br />

∂v<br />

g 1 x 1<br />

=−<br />

∂ +<br />

0<br />

f<br />

y0 2cosθy tgθx 3/2<br />

2<br />

( 0<br />

f )<br />

(9)<br />

(8)<br />

∂v0<br />

g 1 1 1 g 1 xtgθ<br />

= −<br />

∂ x f 2cosθ y + tg x 2 2cosθ<br />

y + tg x<br />

f<br />

1/2 3/2<br />

( 0 θ f ) ( 0 θ f )<br />

∂v<br />

g tgθ x 1 g 1 x x<br />

= −<br />

( y0 tgθxf ) ( y0 tgθxf )<br />

(10)<br />

0<br />

f f<br />

f<br />

1/2 3/2 2<br />

∂ θ 2 cosθ + 2 2 cosθ + cos θ (11)<br />

La primera <strong>de</strong>rivada es relativamente simple ya que y0 aparece sólo una vez en la raíz en<br />

el <strong>de</strong>nominador y se pue<strong>de</strong> buscar en tablas.<br />

En la <strong>de</strong>rivada con respecto a xf, éste aparece tanto en el numerador como en el<br />

<strong>de</strong>nominador, y por eso al <strong>de</strong>rivar aparecen dos términos.<br />

Finalmente la <strong>de</strong>rivada respecto a θ es más compleja porque aparecen dos funciones<br />

trigonométricas. La estrategia que se podría usar es suponer a la función v0 como<br />

producto <strong>de</strong> dos funciones A(θ) y B(θ), <strong>de</strong> manera que v0= A(θ).B(θ) y siendo<br />

g x f g<br />

A( θ ) = = xfsecθ(12)<br />

2cosθ 2<br />

1<br />

B(<br />

θ ) =<br />

(13)<br />

y + tgθx 0<br />

f<br />

De esta manera se pue<strong>de</strong> aplicar la <strong>de</strong>rivación en un producto <strong>de</strong> manera que<br />

∂v0<br />

= A'B+ AB'<br />

∂θ<br />

(14)


(acá la comilla significa <strong>de</strong>rivada con respecto a θ y se utilizó esta notación por<br />

simplicidad).<br />

La <strong>de</strong>rivada A’ es la <strong>de</strong>rivada <strong>de</strong> una secante que se pue<strong>de</strong> obtener <strong>de</strong> tablas.<br />

Para la <strong>de</strong>rivada <strong>de</strong> B se pue<strong>de</strong> hacer un paso adicional que es suponer B compuesta por<br />

una función C(θ). Entonces:<br />

B =<br />

siendo:<br />

1<br />

y0+ C( θ ) xf<br />

C(θ)= tgθ (16)<br />

(15)<br />

Ahora, para calcular B’ se pue<strong>de</strong> aplicar la regla <strong>de</strong> la ca<strong>de</strong>na.<br />

B'= B'( C( θ )) . C'(<br />

θ ) (17)<br />

De esta manera se llega a la ecuación (11).<br />

Finalmente el error en v0 se calcula introduciendo las ecuaciones 9, 10 y 11 en la<br />

ecuación 8:<br />

g 1 x f 1<br />

Δ v0 = − Δ y<br />

3/2 0 +<br />

2cosθ θ 2<br />

( y0+ tg xf)<br />

g 1 1 1 g 1<br />

+ −<br />

2cosθ 2 2cosθ<br />

xtgθ<br />

f<br />

1/2 3/2<br />

( y0 + tgθxf ) ( y0 + tgθxf )<br />

g tgθ xf 1 g 1 xf<br />

f<br />

1/2 3/2 2<br />

2cos 2 2cos cos<br />

x<br />

+ −<br />

Δθ<br />

θ + θ + θ<br />

( y0 tgθxf ) ( y0 tgθxf )<br />

Δ x +<br />

Para evaluar esta ecuación, es <strong>de</strong>cir evaluar los valores numéricos <strong>de</strong> y0, xf y θ, se <strong>de</strong>ben<br />

usar los valores que se midieron experimentalmente o los promedios si se realizaron<br />

muchas mediciones. Los <strong>errores</strong> Δxf, Δy0 y Δθ se <strong>de</strong>ben evaluar usando los <strong>errores</strong> <strong>de</strong><br />

los instrumentos o los <strong>errores</strong> estadísticos según el caso (siempre se usa el mayor entre<br />

error <strong>de</strong>l instrumento o error estadístico).<br />

Cabe aclarar que para Δθ <strong>de</strong>be usarse el valor en radianes, es <strong>de</strong>cir si el error era <strong>de</strong> 2º<br />

<strong>de</strong>bería hacerse la conversión:<br />

π<br />

.2 =<br />

0.0349<br />

180<br />

f<br />

(18)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!