Definición y concepto

La propagación de incertidumbre se define como el efecto que ejercen las variables de incertidumbre sobre la incertidumbre resultante de una función matemática dependiente de ellas. Este concepto es fundamental en estadística y en la medición experimental, donde las variables suelen representar valores obtenidos a través de mediciones sujetas a limitaciones intrínsecas del proceso de medición. Cuando estas variables se combinan dentro de una función, sus respectivas incertidumbres se propagan, afectando la precisión del valor final calculado.

Representación de la incertidumbre

La incertidumbre puede definirse comúnmente mediante el error absoluto o el error relativo, aunque con mayor frecuencia se cuantifica utilizando la desviación estándar. El error absoluto indica la magnitud del desvío en las mismas unidades que la variable medida, mientras que el error relativo expresa esta desviación como una fracción o porcentaje del valor medido, facilitando la comparación entre magnitudes de diferente orden. La desviación estándar ofrece una medida estadística de la dispersión de los datos alrededor de la media, proporcionando una estimación robusta de la incertidumbre en conjuntos de mediciones repetidas.

En la práctica experimental, un valor medido y su asociada incertidumbre se expresan típicamente en la forma x ± Δx, donde x representa el valor más probable (generalmente la media) y Δx representa la magnitud de la incertidumbre. Esta notación comunica no solo el resultado de la medición, sino también el rango dentro del cual se espera que se encuentre el valor verdadero con un cierto nivel de confianza.

Intervalo de confianza y distribución normal

El concepto de intervalo de confianza está estrechamente ligado a la interpretación de la desviación estándar. Para una distribución normal (gaussiana), aproximadamente el 68% de las mediciones caen dentro de una desviación estándar de la media. Esto significa que si se repite la medición muchas veces bajo las mismas condiciones, se espera que el valor verdadero se encuentre dentro del intervalo x ± σ (donde σ es la desviación estándar) en el 68% de los casos. Este porcentaje aumenta a aproximadamente el 95% para dos desviaciones estándar y al 99.7% para tres desviaciones estándar, proporcionando una base cuantitativa para evaluar la fiabilidad de los resultados experimentales.

¿Cómo se calcula la propagación en combinaciones lineales?

El cálculo de la propagación de incertidumbre en combinaciones lineales de n variables se fundamenta en el álgebra de matrices, permitiendo cuantificar cómo las incertidumbres individuales y sus interrelaciones afectan al resultado final de la función. Este enfoque es particularmente útil cuando la función matemática que relaciona las variables de entrada puede expresarse como una suma ponderada o una combinación lineal directa.

Matriz de covarianza y cálculo de la incertidumbre

Para describir la incertidumbre de un conjunto de variables, se utiliza la matriz de covarianza, denotada como Mx. Esta matriz cuadrada de dimensión n x n contiene las varianzas de cada variable en la diagonal principal y las covarianzas entre pares de variables en las posiciones fuera de la diagonal. La estructura de esta matriz captura tanto la magnitud del error absoluto o relativo de cada medición como la correlación estadística entre ellas.

Estructura de la matriz de covarianza Mx para n variables
Varianza de x1 Covarianza (x1, x2) ... Covarianza (x1, xn)
Covarianza (x2, x1) Varianza de x2 ... Covarianza (x2, xn)
... ... Varianza de xi ...
Covarianza (xn, x1) Covarianza (xn, x2) ... Varianza de xn

La matriz de covarianza de la función resultante, Mf, se calcula mediante operaciones matriciales que involucran la matriz de covarianza de las variables de entrada y los coeficientes de la combinación lineal. Específicamente, si la función depende linealmente de las variables con coeficientes representados por el vector o matriz a, la relación se expresa como Mf = a^T Mx a. Esta fórmula permite transformar las incertidumbres de las variables originales al espacio de la función objetivo, preservando la información sobre las correlaciones.

Simplificación para variables no correlacionadas

En muchos casos prácticos de medición experimental, las variables pueden considerarse estadísticamente independientes, lo que implica que sus errores no están correlacionados. En esta situación, las covarianzas entre pares de variables distintas se aproximan a cero, simplificando significativamente la matriz de covarianza Mx, que se vuelve esencialmente diagonal. Esta simplificación elimina los términos de covarianza de las fórmulas de propagación, reduciendo la complejidad del cálculo de la incertidumbre total. La desviación estándar de la función puede estimarse entonces sumando las varianzas ponderadas de cada variable individual, facilitando la interpretación de los resultados en estadística y análisis de datos experimentales.

Tratamiento de funciones no lineales

El tratamiento de funciones no lineales representa una etapa fundamental en la propagación de errores, ya que la mayoría de las relaciones entre variables en la práctica experimental no siguen una dependencia estrictamente lineal. Para abordar esta complejidad, se emplea la técnica de linealización mediante el desarrollo en serie de Taylor de primer orden. Este método permite aproximar una función arbitraria alrededor de los valores medios de las variables independientes, transformando el problema en uno manejable mediante las reglas de combinación de varianzas.

Linealización mediante la serie de Taylor

La expansión de la serie de Taylor permite expresar una función f(x1, x2,..., xn) en términos de sus derivadas parciales evaluadas en el punto medio. Al truncar la serie en el primer orden, se desprecian los términos de segundo orden y superiores, lo que introduce una aproximación válida siempre que las incertidumbres de las variables sean pequeñas en comparación con la curvatura de la función. Esta aproximación es la base de la mayoría de las fórmulas de propagación de incertidumbre utilizadas en física y química experimental.

En esta expansión, la función se descompone en una parte constante y términos lineales. Es crucial notar que el término constante, que corresponde al valor de la función evaluada en los promedios de las variables, no contribuye a la incertidumbre total. Esto se debe a que la incertidumbre mide la dispersión o variabilidad de las mediciones, y una constante no varía. Por lo tanto, solo los términos que dependen de las variables aleatorias afectan al resultado final de la propagación de errores.

El papel de las derivadas parciales

En el caso de funciones lineales simples, los coeficientes que multiplican a las variables son constantes fijas. Sin embargo, en la linealización de funciones no lineales, estos coeficientes son reemplazados por las derivadas parciales de la función con respecto a cada variable independiente. Estas derivadas actúan como factores de escala que indican cuánto cambia el resultado de la función ante una pequeña variación en cada una de las variables de entrada.

La magnitud de cada derivada parcial determina el peso que tiene la incertidumbre de esa variable específica en la incertidumbre total de la función. Si una variable tiene una derivada parcial grande, pequeñas variaciones en su medición tendrán un impacto significativo en el resultado final. Por el contrario, si la derivada parcial es pequeña, esa variable contribuye menos a la incertidumbre global, incluso si su error de medición es relativamente alto.

Esta sustitución de coeficientes por derivadas parciales permite generalizar las fórmulas de propagación de errores a casi cualquier función matemática continua y diferenciable. El resultado es una expresión donde la varianza de la función se calcula como una combinación lineal de las varianzas y covarianzas de las variables originales, ponderadas por el cuadrado de sus respectivas derivadas parciales. Este enfoque proporciona una herramienta poderosa y versátil para cuantificar la precisión de mediciones indirectas en ciencias experimentales.

Ejercicios resueltos

La aplicación práctica de la propagación de incertidumbre requiere el uso sistemático de las fórmulas derivadas de la linealización mediante series de Taylor. A continuación, se presentan ejemplos resueltos que ilustran el procedimiento para funciones lineales y no lineales, así como la tabla de referencia para operaciones básicas.

Ejemplo 1: Propagación en una función no lineal general

Considere una función matemática f(a,b) dependiente de dos variables aleatorias a y b. Para estimar la varianza de f, se utiliza la expansión en serie de Taylor de primer orden alrededor de los valores medios de las variables. La fórmula general para la varianza σf2​ incluye las derivadas parciales y la covarianza entre las variables:

σ f 2 ≈ ( ∂ f ∂ a ) 2 σ a 2 + ( ∂ f ∂ b ) 2 σ b 2 + 2 ∂ f ∂ a ∂ f ∂ b ( cov ( a, b ) )

Esta expresión demuestra cómo la incertidumbre de las variables individuales y su correlación afectan el resultado final. El término de covarianza es crucial cuando las variables no son independientes.

Ejemplo 2: Caso particular de la multiplicación

Para la función de multiplicación f=ab, se calculan las derivadas parciales: ∂a∂f​=b y ∂b∂f​=a. Sustituyendo estas derivadas en la fórmula general de varianza y asumiendo que las variables son no correlacionadas (covarianza cero), se obtiene:

σ f 2 = b 2 σ a 2 + a 2 σ b 2

Para expresar el error relativo, se divide toda la ecuación por f2=a2b2:

σ f f = ( σ a a ) 2 + ( σ b b ) 2

Esto indica que la varianza relativa de la función es la suma de las varianzas relativas de las variables individuales.

Tabla de fórmulas de propagación para variables A y B

Función f(A, B) Varianza σf2​ (No correlacionadas) Varianza σf2​ (Con correlación)
Suma: A+B σA2​+σB2​ σA2​+σB2​+2cov(A,B)
Resta: A−B σA2​+σB2​ σA2​+σB2​−2cov(A,B)
Multiplicación: A⋅B B2σA2​+A2σB2​ B2σA2​+A2σB2​+2ABcov(A,B)
División: A/B B21​σA2​+B4A2​σB2​ B21​σA2​+B4A2​σB2​−B32A​cov(A,B)

Aplicaciones prácticas en medición

La aplicación práctica de la propagación de incertidumbre es fundamental en la medición experimental, donde las variables de entrada poseen incertidumbre debido a las limitaciones de medición. Este efecto se propaga a la combinación de variables en la función matemática resultante, determinando la precisión final del resultado. La incertidumbre se define comúnmente por el error absoluto o relativo, y más frecuentemente por la desviación estándar, lo que permite cuantificar el rango de confianza de la medida.

Cálculo de la resistencia eléctrica

Un ejemplo clásico es la determinación de la resistencia eléctrica utilizando la ley de Ohm, expresada como R = V/I, donde V representa el voltaje e I la corriente. Al medir estas magnitudes, cada una presenta una incertidumbre asociada, denotada como ΔV y ΔI respectivamente. Para calcular la incertidumbre propagada en la resistencia, ΔR, se aplican las reglas de propagación de errores para funciones no lineales, que requieren linealización mediante series de Taylor de primer orden.

En este caso, la incertidumbre relativa de la resistencia depende de las incertidumbres relativas del voltaje y la corriente. Si las variables son no correlacionadas, las fórmulas de propagación se simplifican significativamente, eliminando los términos de covarianza. Esto implica que la varianza de la resistencia puede estimarse sumando las contribuciones individuales de las incertidumbres del voltaje y la corriente, ponderadas por las derivadas parciales de la función R con respecto a cada variable.

Funciones trigonométricas y derivadas parciales

Otro caso relevante es el cálculo de la incertidumbre en funciones trigonométricas, como la función tangente inversa. Al utilizar derivadas parciales para linealizar la función, se puede estimar cómo las incertidumbres en las variables de entrada afectan el resultado final. Este método introduce posibles sesgos en la estimación del error, especialmente cuando las funciones son altamente no lineales o las incertidumbres son grandes. La precisión de la estimación depende de la calidad de la linealización y de la independencia de las variables medidas.

La matriz de covarianza de las funciones se calcula mediante operaciones matriciales con la matriz de covarianza de las variables, lo que permite una evaluación más completa de las interacciones entre las incertidumbres. Este enfoque es esencial en mediciones complejas donde múltiples variables influyen en el resultado final, asegurando que la incertidumbre total refleje adecuadamente las limitaciones experimentales.

¿Qué limitaciones tienen estos métodos?

Sesgo en funciones no lineales

Los métodos estándar de propagación de errores dependen críticamente de la linealización de las funciones mediante series de Taylor de primer orden. Esta aproximación introduce sesgos sistemáticos en la estimación de la incertidumbre cuando las funciones subyacentes presentan no linealidades significativas. El error de truncación en la expansión de Taylor significa que los términos de orden superior, al ser ignorados, afectan la precisión de la varianza resultante. Este sesgo depende intrínsecamente de la naturaleza funcional específica; por ejemplo, en transformaciones logarítmicas como log(x), la asimetría introducida por la no linealidad puede distorsionar la distribución de la variable resultante si la incertidumbre relativa de x es grande. Por tanto, la suposición de normalidad en la variable derivada puede volverse cuestionable, requiriendo correcciones o métodos más robustos para mantener la validez estadística.

Correlación en parámetros de regresión

En el contexto del análisis de datos experimentales, particularmente en la regresión lineal, es fundamental reconocer que los errores asociados a los parámetros estimados rara vez son independientes. Específicamente, la pendiente y la intercepción de una línea de regresión suelen presentar una correlación significativa. Ignorar esta covarianza al calcular la incertidumbre de predicciones futuras o al combinar estos parámetros en nuevas funciones matemáticas conduce a una subestimación o sobreestimación de la varianza total. La matriz de covarianza de los parámetros debe incorporarse explícitamente en las operaciones matriciales de propagación para capturar adecuadamente la estructura de la incertidumbre conjunta, asegurando que los intervalos de confianza reflejen la verdadera dispersión de los datos.

Casos patológicos y alternativas numéricas

Existen casos donde las suposiciones clásicas de la propagación de errores se rompen drásticamente. Un ejemplo crítico es la inversión de una variable aleatoria B que sigue una distribución normal estándar, resultando en la variable 1/B. Esta transformación genera una distribución de Cauchy, la cual carece de una varianza finita definible bajo las condiciones estándar. En tal escenario, aplicar fórmulas de propagación lineales lleva a resultados matemáticamente inconsistentes o a una divergencia de la varianza. Para abordar estas limitaciones y otros casos complejos de no linealidad, se emplean alternativas numéricas como la simulación Monte Carlo, que muestrea directamente la distribución conjunta de las variables, o la transformación de Geary-Hinkley, que ajusta la distribución de la función resultante para mejorar la aproximación normal. Estas herramientas permiten manejar distribuciones no gaussianas y dependencias complejas que escapan al alcance del análisis analítico de primer orden.

Contexto histórico y referencias

El desarrollo de la propagación de errores está intrínsecamente ligado a la evolución de la estadística experimental y la necesidad de cuantificar la precisión en las mediciones científicas. Este concepto surge de la observación de que las variables utilizadas en funciones matemáticas poseen incertidumbre inherente, derivada de las limitaciones propias de los instrumentos y métodos de medición. Cuando estas variables se combinan, sus respectivas incertidumbres no desaparecen, sino que se propagan a través de la función, afectando el resultado final. Este principio es fundamental para interpretar datos en disciplinas como la física, la química y la ingeniería, donde la distinción entre el valor medido y el valor verdadero depende directamente del manejo adecuado del error.

Fuentes académicas y manuales de laboratorio

La enseñanza y la aplicación práctica de estos conceptos se han consolidado a través de diversas publicaciones académicas que buscan estandarizar el análisis de datos experimentales. Entre los recursos educativos destacados, el Manual de Laboratorio de Mecánica de la Case Western Reserve University ofrece un enfoque práctico para comprender cómo las incertidumbres se transmiten en contextos mecánicos específicos. Este tipo de manuales es esencial para estudiantes e investigadores que requieren aplicar las fórmulas de propagación en entornos experimentales reales, proporcionando ejemplos concretos de cómo calcular el error absoluto y relativo en mediciones físicas.

En el ámbito de la investigación contemporánea, trabajos como el realizado por Mathieu Rouaud en 2013 sobre probabilidad y estadística aportan perspectivas actualizadas sobre el tratamiento de la incertidumbre. Estos estudios refuerzan la importancia de considerar no solo las desviaciones estándar individuales, sino también las relaciones entre variables, como la covarianza, para obtener estimaciones más precisas. La integración de estos enfoques teóricos y prácticos permite a los científicos evaluar la fiabilidad de sus resultados con mayor rigor, asegurando que las conclusiones derivadas de los datos experimentales estén respaldadas por un análisis estadístico sólido y bien fundamentado.

Referencias

  1. «Propagación de errores» en Wikipedia en español
  2. Propagation of Error - Wolfram MathWorld
  3. Guide to the Expression of Uncertainty in Measurement (GUM)
  4. Error Propagation - HyperPhysics (Georgia State University)
  5. Uncertainty in Measurement - NIST (National Institute of Standards and Technology)