UNIVERSIDAD NACIONAL DE COLOMBIA (SEDE BOGOTÁ) 1

Cálculo de su Incertidumbre de Medición

Juan Diego Ocampo, Juan Manuel Marı́n
Universidad Nacional de Colombia (Sede Bogotá), Facultad de Ingenierı́a
Departamento de Ingenierı́a Eléctrica y Electrónica, Instrumentación y Medidas
[email protected], [email protected]

Abstract—En este informe se presenta el diseño de un sensor IV. DATASHEETS

de temperatura mediante el uso de un Termistor NTC-3950 junto
al cálculos detallados de la incertidumbre asociada a su medida.
Ver anexos:
Se presentan también los resultados esperados en las mediciones • http://www.bapihvac.com/wp/content/uploads/2010/11-
para un rango de temperaturas entre 17o C y 17o C junto con las Thermistor/100K.pdf
respectivas gráficas de VoutT vs. T y las Curvas S para los rangos • https://www.murata.com/media/webrenewal/support-
de temperatura 17o C-17o C y 0o C-100o C.
library/catalog/products/thermistor/ntc/r44e.ashx
Index Terms—Temperatura, incertidumbre, termistor, instru- • http://www.farnell.com/datasheets/33552.pdf
mento, medición, acondicionamiento, escala de temperatura ISO-
IEC17025, constante térmica, divisor de tensión, resistencia de
carga, curva S. V. C IRCUITOS ESQUEM ÁTICOS

I. I NTRODUCCI ÓN

El propósito de este experimento es diseñar un sensor de

temperatura utilizandp un termistor y aprender a realizar los
cálculos de las incertidumbres que tiene el instrumento de
medición que diseñemos mediante herramientas estadı́sticas y
otras no estadı́sticas.

Para hallar la incertidumbre tipo A de las mediciones

tomadas se utilizan métodos estadı́stico que consisten en
hallar las desviaciones estándar de cada serie de mediciones.
Las incertidumbres tipo B son aquellas que no están Fig. 1: Circuito esquemático del Sensor de Temperatura
determinadas por patrones estadı́sticos, como puede ser la diseñado en donde RT = 100kΩ a T = 25o C.
resolución del instrumento.

II. M ATERIALES VI. P ROCEDIMIENTO

Cálculo de la incertidumbre en la medición
• Resistores de precisión de 20kΩ ± 1%, 2, 2kΩ ± 1% y
470Ω ± 1% respectivamente. El circuito de alimentación del sensor es un circuito divisor
• Fuente de tensión DC SIGLENT SPD3303. de tensión, de forma que sea posible modificar los valores
• Termistor de 100kΩ NTC-3950. de RL y de la fuente Vi para poder lograr cierto grado
• Protoboard. de linealidad en la medida del sensor dentro del rango de
• Multı́metro True RMS Fluke 177. medición, que será de (17 o C − 37 o C).

Sabemos que el circuito consta de una resistencia de

III. M ARCO T E ÓRICO carga RL , una resistencia asociada al termistor RT y una
fuente de alimentación Vi y que RT viene dada por la
Incertidumbre siguiente ecuación.
En metrologı́a la incertidumbre es un parámetro que carac- 1 1
RT = Ro eB( T − T0 ) [Ω]
teriza el intervalo de valores dentro del cual se espera que esté
el valor de la cantidad que se mide. Pueden ser incertidumbres Si se asume una temperatura ambiente de 17 o C, la cual
de tipo A (aquella quw pueden ser evaluada a partir de se encuentra dentro del rango de temperatura a medir, se
datos estadı́sticos, que pueden caracterizarse por desviaciones esperarı́a medir una tensión de salida presente en la resistencia
estándar) o de tipo B (aquellas que no corresponden a la RL de aproximadamente 0.6635[V ], teniendo en cuenta que
clasificación tipo A) [1]. el circuito con una tensión Vi = 5V , que B = 4250 o K,
UNIVERSIDAD NACIONAL DE COLOMBIA (SEDE BOGOTÁ) 2

T0 = 298.15 o K y que RT0 = 100kΩ. Ahora:



∂Vo 5 · (470 + 20000 + 2200)
Puesto que el circuito que alimenta el sensor es sencillamente θB =
∂B
=
1 1



B· −
un divisor de tensión basa con despejar la variable T de (470 + 20000 + 2200) + 10000e
273.15+17 298.15 
B

la ecuación de la resistencia del termistor en función de la 

0.007194 · (1.00327)B 

temperatura (RT ). La ecuación del divisor de tensión es la θB =
2  = −5.3221 · 10−5
B
siguiente: (1.00327) + 0.441112  B=4250
Vf · RL
Vo =

RL + RT Vf · (470 + 20000 + 2200)

∂Vo 
θsrc = =

∂Vf 4250· 1 − 1
Al hacer uso de estas ecuaciones y reconociendo que para el (470 + 20000 + 2200) + 10000e
273.15+17 298.15 
Vo
montaje en protoboard se utilizaron 3 resistores de precisión
(de 20kΩ ± 1%, 2, 2kΩ ± 1% y 470Ω ± 1% respectivamente) θsrc = 0.1327
como resistencia de carga, se tiene que la temperatura viene
dada de la siguiente forma. θdisplay = θsrc = 0.1327
V
1 (R1L + R2L + R3L ) Vfo − (R1L + R2L + R3L ) 1 −1 θmultimeter = 1
T =( ln( )+ )
B R0 T0
Por consiguiente solo resta calcular la incertidumbre total de la
Una vez conocida la ecuación que relaciona la tensión de medición, teniendo en cuenta que esta debe estar multiplicada
salida del circuito con la temperatura se procede a realizar por 2, puesto que se busca hallar la incertidumbre al 95%.
el cálculo de incertidumbre. Para ello se deben conocer en q
primer lugar las incertidumbres asociadas a los componentes UB = 2· 3(URL θRL )2 + (UB θB )2 + (Us θs )2 + (Ud θd )2 + (Um θm )2
del montaje y de los instrumentos utilizados.
UR1L (20kΩ) = ±200Ω UB = ±0.018976V

UR2L (2.2kΩ) = ±22Ω VRLout = 0.6635 ± 0.018976 V 95%

UR3L (470Ω) = ±4.7Ω Urelativa = ±2.5%

UB = ±127.5o K
Error
Usource = ±0.004V Para calcular el error presente en la medición (ocasionado
por la disipación de potencia térmica asociada al termistor) se
Udisplay = ±0.01V
debe calcular en primera instancia la potencia disipada por la
Umultimeter = ±0.00018V resistencia de carga.
Ahora se procede a calcular el factor de sensibilidad asociado Vo2
a cada una fuente de incertidumbre tipo B. P = = 29.2668µW
RL

P


5 · RL1 + 2200 + 470
θR =
∂Vo
=
 =
L1 ∂RL1
4250· 1 − 1

 K
273.15+17 298.15 
RL1 + 2200 + 470 + 10000e
RL1 Donde K es la constante de disipación térmica, que para el
 termistor NTC-3950 es igual a 2.7mW/ o C. Por lo tanto, se
∂Vo 0.04751
= 3.2177 · 10−5

θRL1 = = 
2 tiene que el error asociado a la medición es:
∂RL1 (RL1 + 2670.01) RL1 =20000

  = 0.01084o C

5 · RL2 + 20000 + 470
∂Vo

θR = =

L2 ∂RL2
4250· 1
273.15+17
− 1
298.15 
VII. C ONCLUSIONES
RL2 + 20000 + 470 + 10000e
RL2
 • La curva Vo vs. T no es completamente lineal en ningún
∂Vo 0.04751 punto, sin embargo, al cambiar valores de tension y
= 3.2177 · 10−5

θRL2 = = 
∂RL2 (RL2 + 20470) 2 resistencia RL , es posible obtener secciones de la curva
RL2 =2200
 muy semejantes a una lı́nea recta.
∂Vo
5· RL3 + 20000 + 2200


 • A partir de la ecuación caracterı́stica de un termistor es
θR
L3
=
∂RL3
=

4250· 1 − 1

 posible calcular la temperatura registrada por el mismo
273.15+17 298.15 
RL3 + 20000 + 2200 + 10000e
RL3
en base a la tensión de salida en el circuito VT (la del
 termistor).
∂Vo 0.04751 Al cambiar los valores de RL y Vs es posible jugar con
= 3.2177 · 10−5

θRL3 = =  •
2
∂RL3 (RL3 + 22200) RL3 =2200 la inclinación y la pronunciación de la curva S, todo
Como podemos ver: ello a favor de los intereses del diseñador del circuito de
linealización, donde se espera obtener una región lo más
θRL1 = θRL3 = θRL3 = 3.2177 · 10−5 lineal posible en un rango determinado de temperatura.
UNIVERSIDAD NACIONAL DE COLOMBIA (SEDE BOGOTÁ) 3

• Es de observar que la curva S obtenida posee un punto de

inflexión, donde el rango de temperaturas aledañas posee
un comportamiento prácticamente lineal, sin embargo se
garantizó que el rango de temperaturas a medir (17 −
37 o C) fuera lo más lineal posible, para ello se jugó con
la pronunciación e inclinación de la curva.

VIII. R ESULTADOS E SPERADOS

Fig. 4: Curva S para el rango de temperatura 17o C-37o C

Fig. 2: Tensión de salida del termistor vs. Temperatura

Fig. 5: Curva VoutT vs. T para el rango de temperatura 0o C-
100o C

Fig. 3: Curva S para el rango de temperatura 0o C-100o C

Referencias
• [1] ”Medidas Indirectas”, Propagación de In-
certidumbres, 2017. [Online]. Available: http:- Fig. 6: Curva VoutT vs. T para el rango de temperatura 17o C-
www2.fisica.unlp.edu.armaterias-FEI-teoriasTM-clase-209-
20exp4-20medida-20indirecta.pdf. [Accessed: 15- Sep- 2017].
37o C