Estabilidad numérica

Estabilidad numérica

Estabilidad numérica

En el subcampo matemático del análisis numérico, la estabilidad numérica es una propiedad de los algoritmos numéricos. Describe cómo los errores en los datos de entrada se propagan a través del algoritmo. En un método estable, los errores debidos a las aproximaciones se atenúan a medida que la computación procede. En un método inestable, cualquier error en el procesamiento se magnifica conforme el cálculo procede. Métodos inestables generan rápidamente basura y son inútiles para el procesamiento numérico.

La estabilidad numérica de un método junto con el número condición (en:condition number) define cuán buen resultado podemos obtener usando métodos aproximados para calcular cierto problema matemático.

Algunas veces un sólo cálculo puede ser logrado de varias maneras, todas las cuales son algebraicamente idénticas en términos de números ideales reales o complejos, pero en la práctica producen resultados diferentes conforme(as) tienen diferentes niveles de estabilidad numérica. Una de las tareas comunes del análisis numérico es tratar de seleccionar algoritmos robustos -- esto es, que tienen una buena estabilidad numérica en un amplio intervalo(range) de situaciones. Éstos métodos están frecuentemente disponibles para usuarios de lenguajes de programación como librerías de computación matemática ( ver mathematical computing libraries). El uso apropiado de librerías (libraries) de computación matemática es usualmente muy superior a algoritmos numéricos "roll you own".

Definición

Dado un algoritmo f(x), con x los datos de entrada y ε el error en los datos de entrada, decimos que el algoritmo es numéricamente estable, es decir que el algoritmo depende continuamente de los parámetros, para el error absoluto si

x - (x + \epsilon) \simeq f(x) - f(x + \epsilon)

y numéricamente estable para el error relativo si

\frac{x - (x + \epsilon)}{x} \simeq \frac{f(x) - f(x + \epsilon)}{f(x)}

Decimos que un algoritmo es numéricamente inestable para el error absoluto si

x - (x + \epsilon) << f(x) - f(x + \epsilon)\

y numéricamente inestable para el error relativo si

\frac{x - (x + \epsilon)}{x} << \frac{f(x) - f(x + \epsilon)}{f(x)}

Notas

Cuando se calculan soluciones numéricas a ciertas ecuaciones diferenciales parciales, la estabilidad se consigue algunas veces incluyendo la difusión numérica. La difusión numérica es un término matemático que asegura que errores de redondeo y de otro tipo en los cálculos se separen(spread out) y no se sumen para causar que el cálculo se desborde "blow up".

La estabilidad numérica es la razón por la cual no se puede normalmente probar (testear) un código numérico como la simulación del clima corriéndolo hacia atrás (backward). Correr el código hacia adelante (forward) incluye métodos numéricos para asegurar que los erroes de aproximación aleatorios se vuelvan (become) menos y menos importantes mientras (as) el cálculo procede asegurando la estabilidad numérica. Correr el código hacia atrás causa que esos mecanismos magnifiquen esos errores generando resultados inútiles.

Cuando se resuleve un problema numérico con un método aproximado, dos tipos de errorres pueden ocurrir:

  • Errores de truncamiento: Sólo se puede hacer un número finito de cálculos. Ejemplos: calculando una función transcendental usando su serie de Taylor, integrando usando una suma finita de rectángulos.
  • Errores de redondeo: Ciertos números necesitan un infinito número de dígitos para ser representados (pi), cuando se redondean estos números los errores de redondeo se propagarán a través del cálculo.
Obtenido de "Estabilidad num%C3%A9rica"

Wikimedia foundation. 2010.

См. также в других словарях:

  • Estabilidad numérica — En el subcampo matemático del análisis numérico, la estabilidad numérica es una propiedad de los algoritmos numéricos. Describe cómo los errores en los datos de entrada se propagan a través del algoritmo. En un método estable, los errores debidos …   Enciclopedia Universal

  • Estabilidad — Saltar a navegación, búsqueda En ciencias,una situación es estable si se mantiene en estado estacionario, es decir, igual en el tiempo y una modificación razonablemente pequeña de las condiciones iniciales no altera significativamente el futuro… …   Wikipedia Español

  • Teoría de la estabilidad — En matemáticas, la teoría de estabilidad estudia la estabilidad de las soluciones de ecuaciones diferenciales y sistemas dinámicos, es decir, examina como difieren las soluciones bajo pequeñas modificiaciones de las condiciones iniciales. La… …   Wikipedia Español

  • Fórmula de Herón — Triángulo de lados a, b, c. En geometría, la fórmula de Herón, descubierta por Herón de Alejandría, relaciona el área de un triángulo en términos de las longitudes de sus lados a, b y c …   Wikipedia Español

  • Factorización QR — Saltar a navegación, búsqueda En álgebra lineal, la descomposición o factorización QR de una matriz es una descomposición de la misma como producto de una matriz ortogonal por una triangular superior. La descomposición QR es la base del algoritmo …   Wikipedia Español

  • Unidad Erlang — El Erlang es una unidad adimensional utilizada en telefonía como una medida estadística del volumen de tráfico. Recibe el nombre del ingeniero danés A. K. Erlang, pionero de la teoría de colas. El tráfico de un Erlang corresponde a un recurso… …   Wikipedia Español

  • Error de aproximación — Saltar a navegación, búsqueda El error de aproximación o error numérico es una medida del ajuste de la medida o cálculo de una magnitud con respecto al valor real o teórico que dicha magnitud tiene. Un aspecto importante de los errores de… …   Wikipedia Español

  • Método de flexibilidad — En ingeniería estructural, el Método de flexibilidad es el clásico método consistente en deformación para calcular fuerzas en miembros y desplazamientos en sistemas estructurales. Su versión moderna formulada en términos de la matriz de… …   Wikipedia Español

  • Socovos — Escudo …   Wikipedia Español

  • Algoritmo de Strassen — En la disciplina matemática del álgebra lineal, el algoritmo de Strassen, llamado así por Volker Strassen, es un algoritmo usado para la multiplicación de matrices. Es asintóticamente más rápido que el algoritmo de multiplicación de matrices… …   Wikipedia Español


Поделиться ссылкой на выделенное

Прямая ссылка:
Нажмите правой клавишей мыши и выберите «Копировать ссылку»