La PizarraMatemática.

Máster y doctorado

Investigación computacional reproducible

Un ordenador puede explorar, aproximar o certificar. Esas tres cosas no prometen lo mismo. Empezaremos con una multiplicación exacta, veremos cómo el redondeo puede borrar una diferencia y aprenderemos qué necesita una cota para permitir una conclusión rigurosa.

1. Declara si buscas un patrón, una aproximación o una prueba

Explorar ejemplos puede sugerir una conjetura. Aproximar estima un valor con una precisión que hay que justificar. Certificar produce evidencia que un procedimiento separado puede comprobar para una afirmación concreta.

Una lista larga de casos favorables no prueba sola una afirmación infinita. Un certificado pequeño puede bastar para una afirmación finita, aunque encontrarlo haya exigido una búsqueda grande.

Un certificado de composición

Un programa propone que 221 no es primo.

  1. Evidencia

    Da los factores 13 y 17.

  2. Multiplicar

    13·17=221.
  3. Condición

    Ambos enteros son mayores que uno y menores que 221.

Resultado

La multiplicación certifica que 221 es compuesto.

Lo comprobamos

No hace falta confiar en la búsqueda que encontró esos factores.

Una tabla no es el mismo tipo de evidencia

Se prueba una propiedad en los primeros diez mil enteros.

  1. Cobertura

    Se conocen esos casos concretos, si se evaluaron correctamente.

  2. Lo que falta

    Quedan enteros positivos fuera de la lista.

  3. Conclusión

    Para universalizar hace falta un argumento que cubra también esos casos.

Resultado

El registro debe decir qué dominio se comprobó.

Lo comprobamos

No se etiqueta como teorema universal solo por aumentar el número de pruebas.

2. Exactitud aritmética no implica exactitud de los datos o del modelo

Los enteros de precisión arbitraria y los racionales permiten representar ciertas cantidades exactamente. Eso evita una tolerancia de coma flotante en una identidad racional. El algoritmo todavía puede implementar otra fórmula o recibir una entrada incorrecta.

Una medida escrita como 1,2 puede tener incertidumbre. Representar el decimal como 6/5 exactamente no demuestra que la magnitud física sea exactamente 6/5. Se distinguen valor escrito, cantidad medida y regla matemática.

Una identidad racional exacta

Sumamos 1/10+1/5.

  1. Denominador

    1/5=2/10.
  2. Suma

    3/10.

  3. Comprobar

    No hay redondeo si se guardan numerador y denominador como enteros exactos.

Resultado

La respuesta racional es exactamente 3/10.

Lo comprobamos

Convertir luego a un formato binario finito puede introducir aproximación.

Convertir una medida no elimina su error

Una longitud se informa como 1,2±0,05.

  1. Dato escrito

    El centro 1,2 es 6/5.

  2. Información

    La longitud declarada está entre 1,15 y 1,25.

  3. Conclusión

    El racional 6/5 representa el centro, no una prueba de longitud exacta.

Resultado

Aritmética exacta y conocimiento exacto son cuestiones distintas.

Lo comprobamos

Se propaga la incertidumbre de entrada cuando forma parte del problema.

3. El redondeo puede borrar información antes de la última operación

La coma flotante usa un conjunto finito de números representables. Una operación puede redondearse y una resta de cantidades cercanas puede dejar expuestos errores previos. Cambiar una fórmula algebraicamente equivalente puede mejorar estabilidad.

El condicionamiento describe sensibilidad del problema a los datos; la estabilidad, cómo el método maneja errores. Más dígitos o dos programas que coinciden no son por sí solos una cota rigurosa. También se revisan desbordamiento, subdesbordamiento y operaciones fuera de dominio.

Una unidad que desaparece

En aritmética binaria64 habitual, sea x=10¹⁶.

  1. Suma

    x+1 se redondea al mismo número que x.

  2. Expresión directa

    √(x+1)−√x resulta cero en esa evaluación.

  3. Valor matemático

    Para x real positivo, la diferencia exacta es positiva.

Resultado

El cero mostrado no demuestra que el valor exacto sea cero.

Lo comprobamos

Se identifica el formato y el redondeo del experimento; no se afirma esto de toda aritmética.

Reescribir antes de calcular

Para x≥0, racionalizamos la diferencia.

  1. Identidad

    \sqrt{x+1}-\sqrt x=1/(\sqrt{x+1}+\sqrt x).
  2. Mismo ejemplo

    La expresión de la derecha da aproximadamente 5·10⁻⁹ en binaria64.

  3. Alcance

    Es una evaluación más estable de ese caso, todavía aproximada.

Resultado

Una identidad permite comprobar el mecanismo del fallo.

Lo comprobamos

Dar una cifra plausible no sustituye un intervalo de error demostrado.

4. Un intervalo riguroso conserva inclusión en cada operación

Un intervalo [a,b] representa todos los valores entre sus extremos. Las operaciones deben contener cualquier resultado de los valores permitidos. Para sumas [a,b]+[c,d]=[a+c,b+d]; para productos se toman mínimo y máximo de los cuatro productos de extremos.

Con extremos de coma flotante se necesita redondeo exterior o un control equivalente. Una librería de bolas puede acumular errores en un radio, pero hay que usar sus garantías y formatos correctamente. Dibujar un margen alrededor de un decimal no lo certifica.

Propagar un producto con ambos signos

x∈[−1,2], y∈[3,4].

  1. Extremos

    −1·3=−3, −1·4=−4, 2·3=6, 2·4=8.
  2. Incluir

    Todo xy está en [−4,8].

  3. Comprobar

    Los extremos −4 y 8 sí se alcanzan con entradas permitidas.

Resultado

Se revisan las cuatro combinaciones.

Lo comprobamos

Multiplicar solo mínimos entre sí y máximos entre sí fallaría aquí.

Repetir una variable puede ensanchar

x∈[0,1]; queremos x−x.

  1. Intervalos independientes

    [0,1]−[0,1]=[−1,1].
  2. Dependencia real

    Las dos x son el mismo valor, así que x−x=0.

  3. Interpretar

    [−1,1] contiene cero y es válido, pero poco preciso.

Resultado

Inclusión no significa cota óptima.

Lo comprobamos

Dividir por un intervalo que contiene cero tampoco se resuelve usando solo sus extremos sin tratar la singularidad.

5. La conclusión debe valer para todo el intervalo certificado

Si x∈[a,b] está demostrado y a>0, entonces x>0.
Si a=0 solo se garantiza x≥0.

Si el intervalo incluye negativos y positivos, no permite decidir signo. El punto medio no reemplaza a los extremos.

Los extremos deciden lo que permite afirmar una cota
Una inclusión permite concluirZoom de una cota de √21,41421,4143Todos sus valores son positivos.Otro dato: [−0,1;0,2]−0,100,2Contiene ambos signos y cero.El centro no basta para decidir.
Dos intervalos con escalas independientes identificadas. Arriba se amplía el intervalo racional de raíz de dos entre uno coma cuatro uno cuatro dos y uno coma cuatro uno cuatro tres, cuyos extremos se justifican al elevar al cuadrado exactamente. Todos sus valores son positivos. Abajo otro dato, no la raíz de dos, pertenece al intervalo menos cero coma uno a cero coma dos; la recta marca cero en su posición proporcional, dentro del intervalo, por lo que no decide signo. Los segmentos representan conjuntos completos, no solo un centro con confianza asignada. La inclusión debe proceder de un cálculo riguroso antes de extraer las conclusiones.

También podemos certificar un número irracional mediante desigualdades exactas. Para acotar √2 entre números positivos basta comprobar sus cuadrados a ambos lados de 2, porque el cuadrado es creciente en los no negativos.

Un intervalo de √2 comprobado con enteros

Proponemos a=14142/10000 y b=14143/10000.

  1. Cuadrados

    14142²=199996164<200000000; 14143²=200024449>200000000.
  2. Dividir

    a²<2<b², y a,b son positivos.
  3. Concluir

    1,4142<√2<1,4143; la anchura es 1/10000.
Resultado

Los extremos se justifican exactamente, sin confiar en dígitos impresos de √2.

Lo comprobamos

Esta cota no certifica automáticamente más cifras de las que permite.

Tres intervalos, tres conclusiones

Inclusiones rigurosas [0,2;0,3], [0;0,1] y [−0,1;0,2], usando coma decimal y punto y coma para separar extremos.

  1. Primero

    El extremo inferior es 0,2>0: positividad estricta.
  2. Segundo

    Puede valer cero: solo no negatividad.

  3. Tercero

    Contiene negativos, cero y positivos: signo no decidido.

Resultado

La respuesta se extrae de todos los valores permitidos.

Lo comprobamos

No se trata una comparación indeterminada como una negación demostrada.

6. Un residual mide la ecuación, no directamente la distancia al resultado

Para una raíz aproximada x̃, el residual es p(x̃). Puede ser pequeño porque p está multiplicada por un factor pequeño, incluso lejos de su raíz. Para relacionarlo con el error necesitamos una cota que use propiedades del problema.

El mismo error puede producir residuales de tamaños muy distintos
Residual pequeño ≠ error pequeñoRaíz 1Candidata 100Distancia: 99q(x)=x−1q(100)=99p(x)=10⁻¹²(x−1)p(100)=9,9·10⁻¹¹Misma raíz; mismo error 99.Hace falta relacionar ambas cosas.
Candidata cien para dos ecuaciones que tienen la misma raíz uno. La recta de arriba representa su distancia de noventa y nueve, con las etiquetas de raíz y candidata. Para q de x igual x menos uno, el residual en cien es noventa y nueve. Para p de x igual diez elevado a menos doce por x menos uno, el residual es nueve coma nueve por diez elevado a menos once. La raíz y la distancia no han cambiado al multiplicar la ecuación por un factor pequeño. La cota de error requiere dividir por la escala correspondiente en este caso lineal, y otros problemas necesitan una relación justificada propia.

En una ecuación lineal p(x)=ε(x−r), con ε≠0, la identidad exacta es |x̃−r|=|p(x̃)|/|ε|. Para funciones generales intervienen otras condiciones; raíces múltiples o escalas pequeñas pueden dificultar la inferencia.

Residual pequeño lejos de la raíz

p(x)=10⁻¹²(x−1) y x̃=100.

  1. Raíz

    La única raíz es 1.

  2. Residual

    p(100)=99·10⁻¹²=9,9·10⁻¹¹.
  3. Error

    |100−1|=99, pese al residual pequeño.
Resultado

Una tolerancia residual aislada no certifica proximidad.

Lo comprobamos

La misma raíz para q(x)=x−1 daría residual 99 en esa candidata.

Usar una relación demostrada

Para ese p se conoce |p(x̃)|≤10⁻¹⁴.

  1. Factor

    |ε|=10⁻¹².
  2. Dividir

    |x̃−1|≤10⁻¹⁴/10⁻¹²=10⁻².
  3. Cota

    Ahora sí se obtiene un error de como máximo 0,01.

Resultado

La cota usa la escala de la ecuación, no solo un decimal pequeño.

Lo comprobamos

En otro p habría que justificar la relación adecuada.

7. Una comprobación independiente puede probar otra cosa que la búsqueda

Factores se revisan multiplicando; candidatos a raíces se sustituyen y se relacionan con un conteo o una factorización completa. Volver a ejecutar el mismo algoritmo puede repetir el mismo error. La independencia útil aporta otra identidad o un verificador más sencillo.

Encontrar varias raíces no demuestra que estén todas. Tampoco observar convergencia al refinar una malla prueba por sí solo una cota. Hace falta un argumento que conecte la salida con la conclusión buscada.

Sustitución frente a cobertura completa

p(x)=x³−x.

  1. Candidata

    x=1 da residual cero, pero eso no dice que sea la única raíz.

  2. Factorizar

    p=x(x−1)(x+1).
  3. Todas

    Sobre ℝ, un producto de estos factores solo se anula en −1,0,1.

Resultado

La factorización prueba también que no faltan raíces.

Lo comprobamos

Para otros polinomios podría usarse aislamiento y conteo con hipótesis y certificados propios.

Un error de integración que sí podemos demostrar

Regla del trapecio uniforme para f(x)=x² en [0,1], con n subintervalos.

  1. Suma

    T_n=\frac1n[\tfrac12+\sum_{j=1}^{n-1}(j/n)^2].
  2. Identidad

    Usar la suma de cuadrados da Tₙ=1/3+1/(6n²).
  3. Comparar

    La integral exacta es 1/3; el error es exactamente 1/(6n²).

Resultado

La cota procede de una identidad, no solo de que dos mallas se parezcan.

Lo comprobamos

En n=10 el error es 1/600; cuadruplicar n lo divide por dieciséis.

8. Reproducir requiere conservar entradas, métodos y límites de la conclusión

Se guardan enunciado, datos, código, parámetros, precisión, tolerancias, versiones y procedimiento. Con aleatoriedad se registran generador, semilla y distribución del trabajo. Certificados y resultados deben enlazarse con esa configuración.

Una repetición confirma reproducibilidad de esa ejecución, no la validez del modelo ni una garantía estadística por sí sola. El informe distingue lo explorado, lo aproximado, lo certificado y lo que sigue pendiente de argumento.

Misma semilla, misma salida, misma posible equivocación

Una simulación con un generador y configuración fijados repite una secuencia.

  1. Reproducir

    Otra ejecución puede reconstruir los mismos valores.

  2. Lo que no prueba

    No demuestra independencia de todos los datos ni que la distribución modelada sea la real.

  3. Registro

    Se conserva cómo se generaron y qué hipótesis se usaron al analizarlos.

Resultado

Reproducibilidad es valiosa y tiene un alcance concreto.

Lo comprobamos

No se usa la repetición como sustituto de comprobar el modelo.

Un registro que permite revisar una cota

Resultado: √2 está entre 14142/10000 y 14143/10000.

  1. Datos

    Se guardan los dos racionales, no solo un decimal central.

  2. Certificado

    Se conservan las dos desigualdades enteras de cuadrados.

  3. Verificador

    Otra persona puede multiplicar esos enteros sin repetir la búsqueda del intervalo.

Resultado

El informe explica exactamente por qué la inclusión está demostrada.

Lo comprobamos

El formato de salida también debe conservar la inclusión al imprimir o convertir números.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Un certificado exacto

Un cálculo sugiere que 91 no es primo. Da un certificado verificable.

  1. Afirmación

    Se quiere certificar que 91 no es primo.

  2. Certificado

    Se proponen factores 7 y 13.

  3. Producto

    7·13=91 exactamente.
  4. No trivial

    1<7<91 y 1<13<91.
  5. Concluir

    91 tiene una factorización no trivial: es compuesto.

  6. Comprobar

    No se necesita tolerancia ni confianza en el algoritmo de búsqueda.

Resultado

91=7·13 certifica que es compuesto.

Lo comprobamos

No depende de una tolerancia numérica ni de repetir una búsqueda.

EJEMPLO RESUELTO 2

Una cota que permite afirmar un signo

Un procedimiento riguroso devuelve x∈[1,2; 1,3], es decir, entre 1,2 y 1,3. ¿Puede concluirse x>0?

  1. Dato

    Una inclusión rigurosa sitúa x entre 1,2 y 1,3.

  2. Extremo inferior

    1,2>0.
  3. Orden

    x≥1,2.
  4. Deducir

    x≥1,2>0.
  5. Conclusión

    Sí, x es estrictamente positivo.

  6. Comprobar

    La conclusión usa una inclusión demostrada; un centro aproximado sin cota no sería la misma evidencia.

Resultado

Sí, si la inclusión es rigurosa.

Lo comprobamos

Una salida aproximada 1,25 sin cota no sería el mismo tipo de certificado.

Los símbolos que hemos usado

Experimento computacional
Cálculo usado para explorar un patrón o una pregunta.
Certificado
Objeto verificable que justifica una conclusión mediante reglas declaradas.
Aritmética exacta
Operaciones que representan exactamente las cantidades del dominio empleado.
Intervalo riguroso
Conjunto que se garantiza que contiene el valor objetivo bajo las hipótesis del procedimiento.
Redondeo exterior
Ajuste de extremos que conserva inclusión pese al redondeo.
Residual
Incumplimiento de una ecuación por una aproximación.
Reproducibilidad
Posibilidad de reconstruir el experimento a partir de información conservada suficiente.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar