M étodos Numéricos

1. Métodos Numéricos Juan Manuel Rodríguez Prieto I.M., M.Sc., Ph.D.

2. Integración numérica

3. Integración numérica • Objetivo: aproximar el valor de la integral • Limitaciones de la integración analítica • Las expresiones analíticas de son desconocidas • tiene una integral analítica complicada o desconocida

4. Integración numéricade Newton-Cotes • Métodos que remplazan una función complicada o datos tabulados por un polinomio de aproximación que es fácil de integral • Donde es un polinomio de la forma

5. Integración numéricade Newton-Cotes Aproximación de una integral mediante el área bajo una línea recta.

6. Integración numéricade Newton-Cotes Aproximación de una integral mediante el área bajo una parabola.

7. La regla del trapecio La regla del trapecio es la primera de las formulas de integración de Newton-Cotes. Corresponde al caso donde el polinomio es de primer grado (línea recta). La recta que pasa por los puntos y esta dada por: El área bajo la línea recta f1(x) es una aproximación de la integral de f(x) entre los límites a y b

8. La regla del trapecio Geométricamente, la regla del trapecio equivale a aproximar el área bajo la curva f(x), como el área del trapecio que se forma al unir los puntos

9. Error de la regla del trapecio Cuando usamos la integral bajo un segmento de línea recta para aproximar la integral bajo una curva, se tiene un error que puede ser importante. Una estimación al error de truncamiento para una sola aplicación de la regla del trapecio es Donde esta en algún lugar del intervalo de a a b

10. Ejemplo de la regla del trapecio Aproxime la integral de la curva Entre a=0 y b=0.8 usando la regla del trapecio

11. Ejemplo de la regla del trapecio Aproxime la integral de la curva

12. Ejemplo de la regla del trapecio Aproxime la integral de la curva El área bajo la curva es mucho mayor que el área debajo de la aproximación lineal, de acuerdo a la gráfica mostrada en la diapositiva anterior. ¿Cuanto es el error debido al aproximar la integral del polinomio de grado usando la regla del trapecio?

13. Ejemplo de la regla del trapecio El error porcentual es de ¿Cómo se puede disminuir el error? Se puede dividir el intervalo de interés, en más intervalos. En otras palabras aplicar varias veces la regla del trapecio en el intervalo de interes

14. La regla del trapecio de aplicación múltiple Una forma de mejorar la precisión de la regla del trapecio consiste en dividir el intervalo de integración de a a b en varios segmentos, y aplicar el método a cada uno de ellos. La área asociada a cada uno de los intervalos se suman después para obtener la integral en todo el intervalo. Las ecuaciones resultantes se llaman fórmulas de integración, de aplicación múltiple o compuesta. Vamos a dividir el intervalo de interés en n segmentos del mismo ancho, es decir tendremos n+1 puntos igualmente espaciados.

15. La regla del trapecio de aplicación múltiple Aplicando la regla del trapecio a cada una de las integrales

16. La regla del trapecio de aplicación múltiple Simplificando, se obtiene

17. Erro de la regla del trapecio de aplicación múltiple Simplificando, se obtiene Así, si se duplica el número de segmentos, el error de truncamiento se divide entre 4

18. Ejemplo de la regla del trapecio Aproxime la integral de la curva Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple Usaremos inicialmente n=2, lo cual da un h=0.4

19. Ejemplo de la regla del trapecio Aproxime la integral de la curva Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple Usaremos inicialmente n=2, lo cual da un h=0.4

20. Ejemplo de la regla del trapecio Aproxime la integral de la curva Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple Usaremos inicialmente n=2, lo cual da un h=0.4 Utilizando dos divisiones del intervalo el error ha disminuido de un 89.5% a un 34.9%

21. Ejemplo de la regla del trapecio Aproxime la integral de la curva Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple Usaremos n=3, lo cual da un h=0.2667

22. Ejemplo de la regla del trapecio Aproxime la integral de la curva Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple Usaremos inicialmente n=3, lo cual da un h=0.2667

23. Ejemplo de la regla del trapecio Aproxime la integral de la curva Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple Usaremos inicialmente n=3, lo cual da un h=0.2667 Utilizando dos divisiones del intervalo el error ha disminuido de un 89.5% a un 34.9%

24. Ejemplo de la regla del trapecio múltiple en Matlab n = 5; a = 0; b = 0.8; x = linspace(a,b,n+1); x2 = x.*x; x3 = x2.*x; x4 = x3.*x; x5 = x4.*x; y = 0.2+25*x-200*x2+675*x3-900*x4+400*x5; integral = y(1); for i = 2:n integral = integral + 2*y(i); end integral = integral + y(n+1); h = (b-a)/n; integral = integral*h/2; integral

25. Ejemplo de la regla del trapecio múltiple en Matlab for n= 2:15 a = 0; b = 0.8; x = linspace(a,b,n+1); x2 = x.*x; x3 = x2.*x; x4 = x3.*x; x5 = x4.*x; y = 0.2+25*x-200*x2+675*x3-900*x4+400*x5; integral = y(1); for i = 2:n integral = integral + 2*y(i); end integral = integral + y(n+1); h = (b-a)/n; integral = integral*h/2; resultado(n-1,:)= [n,h, integral]; end Resultados obtenidos para diferentes valores de n A medida que n incrementa, el valor de la integral que se obtiene usando la regla del trapecio múltiple se aproxima a la solución analítica

26. Regla de Simpson 1/3 Además de aplicar la regla del trapecio con una segmentación fina, otra forma de obtener una estimación más exacta de una integral consiste en usar polinomios de grados superior para unir los puntos. Por ejemplo, otro punto entre la mitad entre f(a) y f(b), los tres puntos se pueden unir con una parábola. Si hay dos puntos igualmente espaciados entre f(a) y f(b), los cuatro puntos se pueden unir con mediante un polinomio de tercer grado. Las formulas que resultan de tomar las integrales bajo esos polinomios se conocen como reglas de Simpson

27. Regla de Simpson 1/3

28. Regla de Simpson 1/3Error Así, la regla de Simpson 1/3 es más exacta que la regla del trapecio. Además, es más exacta de lo esperado, porque en lugar de ser el error proporcional a la tercera derivada, el error es proporcional a la cuarta derivada. En otras palabras, da resultados exactos para polinomios cúbicos aun cuando se obtenga un parábola.

29. Aplicación de la regla de Simpson 1/3 Aproxime el valor de la integral de la siguiente función Usando la regla de Simpson 1/3

30. Aplicación de la regla de Simpson 1/3Error Aproxime el valor de la integral de la siguiente función Es aproximadamente 5 veces más precisa que una sola aplicación de la regla del trapecio

31. Regla de Simpson 1/3 de aplicación múltiple La regla de Simpson se mejora al dividir el intervalo de integración en varios segmentos de un mismo tamaño. Se debe utilizar un número par de segmentos para implementar el método Aproxime el valor de la integral de la siguiente función Usando la regla de Simpson 1/3 de aplicación múltiple

32. Regla de Simpson3/8 De manera similar a la obtención de la regla del trapecio y Simpson 1/3, es posible ajustar un polinomio de Lagrange de tercer grado a cuatro puntos e integrarlo La regla 3/8 es más exacta que la regla de 1/3. Por lo general, se prefiere la regla de Simpson 1/3, ya que alcanza una exactitud de tercer orden con tres puntos en lugar de los cuatro puntos requeridos en la versión 3/8. No obstante, la regla de 3/8 es útil cuando el número de segmentos es impar.

33. Regla de Simpson3/8 Aproxime el valor de la integral de la siguiente función Usando la regla de Simpson 3/8

34. Tarea Evalué la integral siguiente: De forma analítica Con una aplicación de la regla de trapecio Con aplicación múltiple de la regla de trapecio n=3 Con una aplicación de la regla de Simpson 1/3 Con una aplicación de la regla de Simpson 3/8