En general, los métodos numéricos no constituyen pruebas. Si todo lo que tenemos es una función desconocida de caja negra $f:\mathbb{R} \to \mathbb{R}$ del que no sabemos nada, y todo lo que podemos hacer es calcular su valor en (finitamente) muchos puntos, entonces simplemente no podemos demostrar que $f$ es positivo.
Sin embargo, en casos concretos, podríamos tener argumentos basados en métodos numéricos que son válidos. Típicamente, tendríamos que hacer una aproximación numérica, y luego probar, utilizando métodos no numéricos que nuestra aproximación numérica es lo suficientemente precisa para que se cumpla el teorema. Por lo tanto, la forma en que los métodos numéricos pueden ayudarnos a demostrar un enunciado es muy específica.
Tomemos, por ejemplo, el siguiente problema: demostrar que $f(x) = x^2 + 1 > 0 \ \forall x \in [-1, 1]$ .
Prueba inválida: Calculamos $f(x)$ en $10^{10^{10000}}$ puntos aleatorios y se utilizó la interpolación lineal entre ellos. Aquí hay un gráfico. Podemos ver que $f(x)$ es siempre positivo.
Prueba válida 1: Calculamos $f(x)$ en los puntos tres puntos: $f(-1) = 2$ , $f(0) = 1$ y $f(1)=2$ . Sea $g(x)$ sea la interpolación lineal de los puntos $(-1, 2)$ , $(0, 1)$ y $(1, 2)$ . $g$ alcanza su mínimo en $g(0) = 1$ . Desde $f^{\prime \prime} = 2$ podemos calcular un límite de error en nuestra interpolación (véase https://ccrma.stanford.edu/~jos/resample/Linear_Interpolation_Error_Bound.html ): $|f(x) - g(x)| \leq \frac{2}{8}$ . Por lo tanto, podemos concluir que $f(x) \geq \frac{3}{4} > 0$ .
Nota: A menudo, si tenemos que recurrir a métodos numéricos, será igual de difícil calcular las derivadas. Sin embargo, no necesitamos las derivadas reales, sólo necesitamos un límite superior. Cuanto mejor sea el límite, menos puntos necesitaremos para evaluar $f(x)$ en. Además, es suficiente con la primera derivada, pero tener la segunda también podría reducir el número de puntos necesarios.
Prueba válida 2: Sabemos que $f(x)$ es convexo. Usamos un método numérico para calcular su mínimo. Encontramos que $\min f(x) \approx 1.0000000075$ . También tenemos un (verdadero, no numérico ) en nuestra aproximación: $|1.0000000075 - \min f(x)| < 0.001$ . Por lo tanto, $f(x) > 1.0000000075 - 0.001 > 0$ .
Por último, no importa realmente si las pruebas analíticas existen o no. La validez de cualquier prueba sólo está determinada por esa prueba y no por otras.
De hecho, se ha demostrado que no todas las afirmaciones verdaderas pueden ser probadas. Pero eso no es motivo para reducir nuestros niveles de rigor.
16 votos
¿Cómo es? probado por métodos numéricos? Existen pruebas informáticas que utilizan métodos numéricos fiables basados en la aritmética de intervalos. Véase scicomp.stackexchange.com/a/1028/640
0 votos
@lhf He editado la pregunta, por favor, échale un vistazo.
12 votos
En realidad, no. Podrías llamarlo conjetura, y citar el cálculo como prueba de que es cierto.
0 votos
@TheoBendit Gracias por tu comentario. En realidad, en la demostración del enunciado muestro que la única forma de demostrarlo son las técnicas numéricas, ya que no se puede demostrar analíticamente.
3 votos
@sane Esto podría ser algo que tendría que ver los detalles para comentar adecuadamente, pero en general, si quieres llamarlo teorema, necesitarás tener una prueba adjunta. Es posible que hayas muestreado 100.000 puntos, pero ¿por qué ese número? ¿Podría cambiar tu opinión si muestreases 100 000 0000 puntos? ¿Por qué no ha tomado una muestra de, por ejemplo, 10 puntos? ¿Qué tiene esta propiedad que le hace estar absolutamente seguro de que 100.000 puntos le dirán definitivamente si una función la tiene o no, pero el muestreo de 10 no? Una prueba no debe dejar lugar a dudas, cuando se verifica.
2 votos
¿Cómo sabe con seguridad que lo que ha observado en sus datos simulados no es una casualidad?
0 votos
@TheoBendit Así que tengo que demostrar que la función $x(\gamma)>0$ sobre el dominio $\gamma \in[0,1]$ . tomo $m=100 000$ puntos $\gamma^1,...\gamma^m$ Por lo tanto, construye la función basada en $m$ puntos mediante técnicas de interpolación. Aunque creo que lo especificado $m$ observaciones son suficientes, sin embargo puedo generar 100 000 000 puntos de $\gamma \in [0,1]$ si cambiará algo.
0 votos
@littleO Estoy seguro porque controlo el proceso de generación de datos. Es decir, manualmente estimo los puntos $\gamma_i$ y especificar el error de estimación (tolerancia).
4 votos
@sane No cambiará las cosas. Aportará más pruebas, pero ningún número (finito) de puntos constituirá una prueba en sí mismo. No tienes forma de saber si la función se agacha por debajo de cero en otros puntos, tal vez en un valor mucho mayor que los valores probados, o tal vez en un punto complicado entre los puntos que has probado. La mayoría de los números no pueden ser expresados por la notación, ¡y mucho menos por los ordenadores! Es decir, la mayoría de los números no tienen ninguna posibilidad de ser probados numéricamente. ¿Cómo sabes que la función no se vuelve negativa ahí?
0 votos
¿Hay algún axioma que haga rigurosa una prueba numérica?
1 votos
Cuando dices que no se puede demostrar, ¿quieres decir que todavía no está en nuestras manos, o que está demostrado que no se puede demostrar?
0 votos
@dmtri Me refiero a que está demostrado que no se puede demostrar.
13 votos
@sane ¿Cómo sabes que no se puede demostrar?
0 votos
@NoahSchweber Obtengo los datos $x_i$ 's de $n$ polinomio de grado, por lo que no puedo resolver para $x_i$ 's para $n>4$ . En su lugar, utilizo métodos numéricos, para estimar $x_i$ 's.
0 votos
¿Así que tienes la función disponible explícitamente y es un polinomio? Entonces puedes usar definitivamente métodos numéricos como la aritmética de intervalos mencionada por @lhf, que te dan una prueba válida. También podrías intentar algo como la programación SOS (suma de cuadrados). Si publicas el polinomio, tal vez podamos decir algo más.
3 votos
En el mejor de los casos, se trata de una prueba empírica, posiblemente de un ajuste excesivo, y en el peor de los casos, de una selección selectiva. Y "teorema" tiene un significado muy específico en matemáticas. Ningún teorema se puede demostrar de la manera que usted piensa.
0 votos
@user21820 Entonces, ¿cómo debo llamar a mis hallazgos?
0 votos
Exactamente como dije. Si entiendes todo las trampas de este tipo de pruebas, entonces tal vez pueda llamarse evidencia empírica. Si no lo hace, lo más probable es que sólo pueda llamarlo "ajuste excesivo" o "selección".
0 votos
@user21820 En realidad, reconozco todas las trampas en este tipo de pruebas, pero por otro lado estoy bastante seguro de que aquí no podemos tener ninguna trampa, ya que como he mencionado el dominio es $[0,1]$ y recojo $100 000$ puntos de ese intervalo. Sin embargo, estoy seguro, pero no es una prueba.
0 votos
@sane Así que para que tengas razón, el dominio de $\gamma$ es $[0, 1]$ pero ¿cuál es el dominio de $x$ y $y$ ? Si se realiza un muestreo, también hay que especificarlo. ¿También están acotados?
0 votos
@sane Vale, y quieres demostrar que $x > 0$ para todos $\gamma \in [0, 1]$ y todos $y \in [0, 1]$ ¿correcto?
2 votos
@sane, el argumento que le das a Noah Schweber no es ciertamente una prueba de que tu afirmación no pueda ser demostrada. Hay muchos resultados (teoremas) relativos a objetos que no pueden (¿todavía?) ser calculados explícitamente.
0 votos
Por favor, actualiza tu pregunta con la ecuación real en lugar de ponerlas en diferentes comentarios, porque si no es difícil de seguir.
1 votos
Hay formas que no son "métodos numéricos" que pueden demostrar que un polinomio es positivo en cualquier parte del dominio $[0,1],$ incluso para un polinomio de grado muy alto. No hay que suponer que no existe una prueba analítica sólo porque el grado sea mayor que $4.$ Al fin y al cabo, no es necesario saber dónde están los ceros, siempre que se sepa que no están en $[0,1].$
2 votos
No puedo creer que nadie lo haya posteado aún, pero no puedo verlo aquí así que: math.stackexchange.com/questions/514/
0 votos
Si pudiéramos, la hipótesis de Riemann ya sería un teorema.
1 votos
@sane: Tu comentario da a entender que, efectivamente, desconoces los escollos. Te reto a que demuestres numéricamente que $x^{x+1/\ln(x)/e}-1>0$ para todos $x[0,1]$ .
0 votos
(convertido de una respuesta a un comentario); tal vez no sea directamente relevante para su pregunta precisa, pero un resultado famoso cuya prueba se basa en cálculos numéricos (rigurosos) es la prueba de Warwick Tucker de la existencia del atractor de Lorenz, por ejemplo ensta-bretagne.fr/jaulin/swim13_slides_tucker.pdf