Ejemplo de propagacion de errores

Se quiere determinar el error en el calculo del area y un lado de un terreno triangular para el que se midieron:
A b a C
18.22 15.31 522535
18.29 15.36 522545
18.24 15.24 522522
18.2 15.38 522530
18.27 15.35
C B 18.31

Primera depuracion:
2
b bi b (b ib ) Zi
18.22 -0.035 0.001225 -0.827259054
18.29 0.035 0.001225 0.827259054
18.24 -0.015 0.000225 -0.354539594
18.2 -0.055 0.003025 -1.299978513
18.27 0.015 0.000225 0.354539594
18.31 0.055 0.003025 1.299978513
Suma 0 0.00895
Promedio 18.255
Numero datos 6
Desviacion est 0.042308392

El valor mas probable del lado "b"

s
x+ t LimiteSuperior
2 n

s
xt LimiteInferior
2 n
{ } { } {}
 s
x+ t LimiteSuperior
2 n

s
xt LimiteInferior
2 n
{ } { } {}

Los valores para la distribucion segn el inervalo de confianza son:

Con los datos que quedaron despues de la depuracion establecemos el intervalo de confianza , que sera mayor si
asumimos un mayor cubrimiento poblacional , para este caso adoptamos un 95%

2
b bi b (b ib )
18.22 -0.035 0.001225
18.29 0.035 0.001225
18.24 -0.015 0.000225
18.27 0.015 0.000225
Suma 0 0.0029
Promedio 18.255
Numero datos 4
Desviacion est 0.031091264

Recordemos que el rango donde se encontrara la media poblacional esta dado por .

s
x+ t LimiteSuperior
2 n

s
xt LimiteInferior
2 n
{ } { } {}

El t para el 90 % es con v= 4-1=3 (grados de libertad): 2.353

2
Evaluando

0. 03109126
18 .255+ 2. 353 LimiteSuperior
4
0 . 03109126
18 .2552 . 353 LimiteInferior
4
{ } { } {}
Limite Superior 18.29157887 error 0.0365788715
Limite Inferior 18.21842113

Redondeamos segn la exactitud de la mediciony el valor mas probable del lado "b" es:

Limite Superior 18.29 error 0.04

Limite Inferior 18.22

El valor mas probable del lado "a"

Primera depuracion:
2
a ai a ( a ia ) Zi
15.31 -0.018 0.000324 -0.425447513
15.36 0.032 0.001024 0.756351135
15.24 -0.088 0.007744 -2.079965621
15.38 0.052 0.002704 1.229070594
15.35 0.022 0.000484 0.519991405
Suma 3.552714E-15 0.01228
Promedio 15.328
Numero datos 5
 Desviacion est 0.055407581

2
a ai a ( a ia )
15.31 -0.03 0.0009
15.36 0.02 0.0004
15.35 0.01 0.0001
Suma -5.32907E-15 0.0014
Promedio 15.34
Numero datos 3
Desviacion est 0.026457513

El t para el 90 % es con v= 3-1=2(grados de libertad): 2.92

2
Evaluando

0 . 02645751
15 .34 +2 . 920 LimiteSuperior
3
0. 02645751
15 .342 . 920 LimiteInferior
3
{ } { } {}
Limite Superior 15.38460374 error 0.0446037368
Limite Inferior 15.29539626

Redondeamos segn la exactitud de la mediciony el valor mas probable del lado "b" es:

Limite Superior 15.38 error 0.04

Limite Inferior 15.30

alor mas probable del angulo "C"

Primera depuracion:
C C i C ( C iC )
2 Zi

35 2 4 0.207763008
45 12 144 1.24657805
22 -11 121 -1.142696546
30 -3 9 -0.311644512
Suma 0 278
Promedio 33
Numero datos 4
Desviacion est 9.626352719

2
C C i C ( C iC )
35 2.5 6.25
30 -2.5 6.25
Suma 0 12.5
Promedio 32.5
Numero datos 2
Desviacion est 3.535533906

El t para el 90 % es con v= 2-1=1(grados de libertad): 6.314

2
Evaluando

3 .53553391
32. 5+6 . 314 LimiteSuperior
2
3 . 53553391
32. 56 . 314 LimiteInferior
2
{ } { } {}
Limite Superior 48.285 error 15.785
Limite Inferior 16.715

Redondeamos segn la exactitud de la mediciony el valor mas probable del lado "b" es:

Limite Superior 48 error 16

Limite Inferior 17

En sintesis los valores mas probables son:

b= 18.26 0.04
a= 15.34 0.04
C= 522533 16

La propagacion de errores para el area y el lado "c" estara dada por:

f ( x )2=G x2GT
donde
 f ( x )2=Errores en la( s) Variable( s) dependientes
G= Matriz de las derivadas parciales de la(s) var iables dependientes respecto
a las independientes u observaciones
x 2=Matriz de los errores en las observaciones o var iables dependientes

A) El primer procedimiento consiste en determinar las variables dependientes

en este caso el rea y el lado "c"
1
A = abSen
2
c= a2 +b 2 2 abCos C

La matriz "G"
A A A
, ,
b a C
c c c
, ,
b a C

ri gh



1 1 1
bSenC , a SenC , abCosC
2 2 2
a bCos C b aCos C ab Sen C
, ,
a2 + b 2 2 a bC os C a2 + b 2 2 a bCo s C a2 + b2 2 ab Cos C

ri gh



[ ]
G=

Evaluando con los valores dados:

A A A
, ,
a b C
c c c
, ,
a b C

righ



1 1 1
18 . 26Sen 52 25 33 , 15 . 34Sen 52 25 33 , 18 . 2615 . 34Cos 52 25 33
2 2 2
15 . 34 18 . 26Cos 52 25 33 18 . 26 15 . 34 Cos 52 25 33 15 . 3418 . 26Sen 52 25 33
, ,
15 . 34 2 + 18 . 26 2 215 . 3418 . 26Cos 52 25 33 15 . 34 2 + 18 . 26 2 215 . 3418 . 26Cos 52 25 33 15 . 342 + 18 . 26 2 215 . 3418 . 26

righ



[ ]
G =

G= 7.234133962 6.07897097 85.37996791

0.279291092 0.5907174 14.72956742

La matriz de los errores en las observaciones o variables independientes es:

a , a a ,b a ,C
b , a b ,b b ,C
,a , b C , C

ri gh



a2
a , b a , C

b , a b 2 b , C

, a , b C
2

ri gh



[ ] [ ]
x2 =

 a , b = b , a =0 ; a , C = C , a =0 ; b , = , b =0
Por lo tanto la matriz quedara:

( 0 . 04 ) 2 0 0
0 ( 0 . 04 ) 2 0
2

0 0 ( ( 16 / 3600 )
180 )

righ



0 . 0016 0 0
0 0 . 0016 0
0 0 6 . 01713 E-09

righ



[ ] [ ]
x 2 =

0.0016 0 0
2 0 0.0016 0
x =
0 0 0.000000006

La matriz transpuesta de G

7.234133962 0.27929109
T
G = 6.078970966 0.5907174
85.37996791 14.7295674

0.011574614 0.00972635 5.137427E-07

2
G x 0.000446866 0.00094515 8.862978E-08

0.14290239 0.00898578
f ( x )2=G x2GT 0.00898578 0.000684426

La interpretacion de los resultados de la matriz

A2
A , c

c, A c
2

2
f (x) =

righ



[ ]

El valor mas probable del Area

 1
VP= abSen A
2
Donde
A= 0.3780
El Area 112.4448

Limite Superior 112.8228

Limite Inferior 112.0668

El valor mas probable del Lado "c"

VP= a2 +b2 2abCos C c

Donde
c= 0.03
El Lado "c" 15.07

Limite Superior 15.09

Limite Inferior 15.04
5 . 34Cos 52 25 33
15 . 3418 . 26Sen 52 25 33
,
3 15 . 342 + 18 . 26 2 215 . 3418 . 26Cos 52 25 33
 Problemas Propuestos
1) La formula para la frialdad producida por el viento "C" (en grados farenheit) es:

0.16 0 .16
C=35 .74 +0 . 6215T 35. 75 v +0 .4275 Tv
Donde "v" es la veloccidad del viento en millas por hora y "T" es la temperatura en grados fahrenheit.
La veloccidad del viento y la temperatura se midieron en cierto lugar obteniendose los siguientes resultados:

Determinar el valor mas probable de la frialdad con un intervalo de

confianza del 95% para las estimaciones de v y T .

Promedio (V) 23.658

Desvest 1.442

v Zi Depurados
24.8 0.79207238 24.8
24.8 0.79207238 24.8
24.8 0.79207238 24.8
24.8 0.79207238 24.8
24.8 0.79207238 24.8
24.8 0.79207238 24.8
22.8 -0.59517303 22.8
21.1 -1.77433163 21.1
25 0.93079692 25
25.1 1.00015919
21 -1.8436939 21
23.6 -0.04027487 23.6
21.7 -1.35815801 21.7
22 -1.15007119 22
22.7 -0.6645353 22.7
21.5 -1.49688255 21.5
22.1 -1.08070892 22.1
26 1.62441962
25.2 1.06952146
24 0.23717422 24
24.7 0.72271011 24.7
22 -1.15007119 22
24.5 0.58398557 24.5
25.1 1.00015919
24.2 0.37589876 24.2
22.7 -0.6645353 22.7
24.4 0.5146233 24.4
 23.2 -0.31772395 23.2
21.6 -1.42752028 21.6
24.6 0.65334784 24.6
23.8 0.09844967 23.8

Estimacin por intervalos

23.4074074
Desviacion estandar 1.367739833
n 27
1- 0.95
ta/2 2.0555294386
error de estimacion 0.5410598578
limite superior 23.95
limite inferior 22.87
En un intervalo de confianza del 95% el valor real de la medicin estara dentro de estos limites
Promedio(T) 2.808
Desvest 1.190

T Zi Depurados
5 1.843
3.4 0.498 3.4
1.9 -0.763 1.9
2 -0.679 2
2.9 0.078 2.9
4.8 1.675
3 0.162 3
2.5 -0.259 2.5
2.1 -0.595 2.1
3.1 0.246 3.1
5 1.843
2.7 -0.091 2.7
2.9 0.078 2.9
1.6 -1.015 1.6
3.1 0.246 3.1
1.1 -1.436 1.1
1.3 -1.267 1.3
1.3 -1.267 1.3
2.7 -0.091 2.7
3.3 0.414 3.3
1.6 -1.015 1.6
1.2 -1.351 1.2
3.8 0.834 3.8
4.7 1.591
3.8 0.834 3.8
2.2 -0.511 2.2
 2.056 intervalo de confianza

2.43181818
Desviacin estandar 0.8481837
n 22
1- 0.95
ta/2 2.07961384
Error de estimacin 0.37606358
Limite superior 2.06
Limite inferior 2.81
de estos limites