Análisis Numérico
Profesor Adjunto
Mg. Martín Conocchiari
Ayudante
Leandro Sueldo
Derivación e Integración
Numérica
Métodos de derivación numérica
Definición:
Aproximación de la derivada en un punto
Sensible a errores de redondeo y truncamiento
Diferentes fórmulas según la cantidad de puntos (+ puntos y midpoint mejor aproximación)
H = distancia entre los puntos
Fórmulas:
2 PUNTOS:
No necesita nodos equidistantes.
Progresiva H > 0 : (f (x0 + h) - f(x0)) / h (Uso punto posterior)
Regresiva H < 0 : (f (x0 - h) - f(x0)) / - h (Uso punto anterior)
Cota de error (truncamiento / método): + |h/2| |max(f’’(ė)| (elijo la mayor, conservador)
La progresiva de un punto en la regresiva del siguiente (se conectan)
Mg. Martín Conocchiari
Métodos de derivación numérica
3 PUNTOS:
Necesita nodos equidistantes.
Endpoint Progresiva: 1 / 2h [ -3f(x0) + 4f(x0+h) – f(x0+2h) ] + [ h^2 / 3 |max f’’’(ė)| ] (usa dos posteriores)
Midpoint: 1 / 2h [ -f(x0-h) + f(x0+h) - [ h^2 / 6 |max f’’’(ė)| ] (usa anterior y siguiente, cota con signo negativo)
Endpoint: 1 / 2h [ f(x0-2h) – 4f(x0-h) + 3f(x0) ] + [ h^2 / 3 |max f’’’(ė)| ] (usa dos anteriores)
5 PUNTOS:
Necesita nodos equidistantes.
Endpoint Progresiva: 1 / 12h [ -25f(x0) + 48f(x0+h) – 36f(x0+2h) + 16f(x0+3h) – 3f(x0+4h) ] + [ h^4 / 5 |max f^v
(ė)| ] (usa cuatro posteriores)
Midpoint: 1 / 12h [ f(x0-2h) - 8f(x0-h) + 8f(x0+h) – f(x0+2h) + [ h^4 / 30 |max f^v (ė)| ] (2 anteriores y 2 post)
Endpoint: 1 / 12h [ 25f(x0) – 48f(x0-h) + 36f(x0 – 2h) – 16f(x0-3h) + 3f(x0-4h)] + [ h^4 / 5 |max f^v (ė)| ] (usa
cuatro anteriores)
Mg. Martín Conocchiari
Error de redondeo y truncamiento
Segunda derivada:
= 1 / h^2 [ f(x0-h) – 2f(x0) + f(x0+h) ] – [ h^2 / 12 |max f^iv (ė)|. (fórmula única)
Error del método (truncamiento):
Es la cota de error del método, prestar atención a los signos. Podemos determinar el sentido del error si es
monótona y mantiene su signo. Si es positiva está subestimando (la sumo en el intervalo de la derecha) y si es
negativa está sobreestimando (la resto en el intervalo de la izquierda)
Error de redondeo:
Es el error de los datos dados, puede ser error de redondeo o truncamiento del dato (no es el error del método).
No podemos determinar el sentido del error, se suma y resta en ambos lados del intervalo. No podemos presumir
que a puntos más cercanos hay menor error, ya que el error de redondeo aumenta (aunque en la mayoría de los
casos conviene usar nodos consecutivos).
Error total:
Es la suma de ambos errores, miramos este error para definir cual fórmula usar. (menor cota)
Mg. Martín Conocchiari
Error de redondeo y truncamiento
Fórmulas de error de redondeo y truncamiento:
Nos indican la cantidad de decimales de redondeo o truncamiento (K)
Truncamiento: ė = 10 x ( 10 ^ - (k+1) )
Redondeo: ė = 5 x ( 10 ^ - (k+1) ) (tiene un 0 más)
2 Puntos: ( 2 x ė ) / h (+ cota en todos)
3 Puntos Endpoint: ( 4 x ė ) / h Tips:
3 Puntos Midpoint: ė / h Siempre fijarse que los puntos elegidos
5 Puntos Endpoint: ( 32 x ė ) / 3h estén interpolando
5 Puntos Midpoint: ( 3 x ė ) / 2h Chequear la “h” si necesito que sean
equidistantes
Siempre evaluar todas las opciones de
cotas antes de elegir un método, ver el
menor con los datos disponibles
Mg. Martín Conocchiari
Métodos de integración numérica
Definición:
Aproximación de la integral en un punto. Lo utilizamos también cuando la integral es muy compleja.
No calculamos error de redondeo en integración
Fórmulas:
Regla de Trapecio:
Polinomio de primer grado (recta). No es buena aproximación para concavidad.
H /2 [f(x0) + f(x1)] – [ (h^3 / 12) x f´´(ė )]
Menor intervalo, error más chico (método estable).
Da un resultado exacto cuando la segunda derivada de la función es = 0. Grado de presición = 1.
H = distancia entre los puntos ( b – a ).
La cota se calcula como la derivada máxima en valor absoluto (conservador).
Determinamos el sentido del error si mantiene su signo
Mg. Martín Conocchiari
Métodos de integración numérica
Regla de Simpson:
Polinomio de segundo grado (curva).
H /3 [f(x0) + 4f(x1) + f(x2)] – [ (h^5 / 90) x f^iv (ė) ]
X0 = a X2 = b X1 = X0 + h h = (b-a) / 2
Error más chico que trapecio. Grado de presición = 3.
Da un resultado exacto cuando la cuarta derivada de la función es = 0.
La cota se calcula como la derivada máxima en valor absoluto (conservador).
Determinamos el sentido del error si mantiene su signo
Error real = valor real - aproximación
Fórmula cerrada de Newton Cotes: Fórmula abierta de Newton Cotes:
De esta provienen Trapecio y Simspon. De esta proviene el Punto Medio
Xi = X0 + h X0 = a Xn = b h = (b-a) / n Xi = X0 + h X0 = a + h Xn = b - h h = (b-a) / n+2
Inluye los límites del intervalo No inluye los límites del intervalo
7
Mg. Martín Conocchiari
Métodos de integración numérica
Punto Medio:
Se aproxima mediante polinomio de grado 0, calcula mediante un rectángulo.
h = (b-a) / 2 X0 = a + h Xn = b – h
2h x f(x0) + [ (h^3 / 3) x f´´ (ė) ]
Es base x altura (2h y valor de la fx en el punto medio del dominio)
En la cota valúa a y b y elijo el valor máximo
Xo y Xn no son a y b (cuidado)
Integración Compuesta:
Definición:
Se utiliza con intervalos de integración grandes, los métodos simples aumentan el error
Divide el dominio en fragmentos pequeños (n intervalos) y aplica Newton cerrado en cada uno
N subintervalos o N+1 puntos
Estabilidad respecto al error de redondeo
8
Mg. Martín Conocchiari
Métodos de integración numérica
Regla de Trapecio Compuesta:
h = (b-a) / n n = (b-a) / h Xj = a + jh
Se suman pequeños errores en cada subintervalo, menor que el método simple
Converge al valor exacto si es polinomio de grado 1 o menos, segunda derivada = 0
H / 2 [ f(a) + 2 x ∑ f(x0 + ih) + f(b) ] – [ h^2 / 12 x (b-a) x f´´(ė)]
La sumatoria va de i = 1 hasta n-1
Cota de error: [(b-a)^3 / (12 x n^2)] x f´´(ė)
Otra: (h^3 / 12) x [n x f´´(ė)]
Despejo n: Raiz cuadrada de { [(b-a)^3 / 12ė] x f´´(ė) }
(el ė suele ser dato)
Siempre evaluar las cotas de los métodos para ver cual usar si
no pide uno específico
Mg. Martín Conocchiari
Métodos de integración numérica
Regla de Simpson Compuesta:
h = (b-a) / n n = (b-a) / h Xj = a + jh
N = 2m P = 2m + 1 m = n/2 (los dividimos)
Cantidad de subintervalos siempre par (n se redondea al próximo número par)
Converge al valor exacto si es polinomio de grado 3 o menos, cuarta derivada = 0, más exacto que Trapecio.
H / 3 x [f(a) + 2 [∑ f(a + 2ih)] + 4 [∑ f(a + (2i + 1) x h)] + f(b) ] – [ (h^5 / 90) x (n/2 x f´iv (ė) ]
La primer sumatoria va de i=1 a (n/2 – 1) y la segunda sumatoria va de i=0 a (n/2 – 1)
Cota de error: (h^4 / 180) x (b-a) x f´iv (ė)
Otra: ((b-a)^5 / 180n^4) x f´iv (ė)
Despejo n: Raiz cuarta de { ((b-a)^5 / 180ė) x f´iv (ė) }
Tengo que recalcular la h y cota por el redondeo, va a ser menor.
10
Mg. Martín Conocchiari
