1 Integración Numérica Ultima actualización: 13/04/20174 Teoría+1 Prácticas+2 Laboratorio MATLAB Tmas. básicos de Cálculo – Desarrollos de Taylor – Sistemas lineales –

2 Descripción del problemaIntegración Numérica Descripción del problema Evaluación de la integral de una función en un intervalo a partir de los valores de la función Función muy compleja (difícil de calcular su integral) No se conoce la función sino sólo sus valores en algunos puntos No existe la primitiva de la función Fórmulas de integración numérica Obtención de los valores de los coeficientes Determinación del error cometido  Descripción  Objetivos  Temario  Bibliografía

3 Entender cómo se aplica el algoritmo de integración de Romberg. Integración Numérica Objetivos Comprender la generación las fórmulas de integración numérica de Newton-Cotes. Calcular las formulas exactas de Newton-Cotes para polinomios de primer, segundo y tercer grado. Entender cómo se aplica el algoritmo de integración de Romberg. Entender la diferencia fundamental entre las fórmulas de Newton-Cotes y la cuadratura gausiana. Comprender el funcionamiento de la integración adaptativa  Descripción  Objetivos  Temario  Bibliografía

4 Teorema del valor medio para integralesIntegración Numérica Introducción Teorema del valor medio para integrales Sean f(x), g(x) C[a, b] donde g(x) es integrable y no cambia de signo en [a,b], entonces existe un punto c en (a,b) tal que  Descripción  Objetivos Temario Introducción NC Reglas simples NC Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía Definición Una fórmula (de integración numérica) es de orden n cuando se anula el error con polinomios de grado menor o igual a n

5 Obtención de las fórmulasIntegración Numérica Obtención de las fórmulas Fórmulas de Newton-Cotes Integrando los polinomios de interpolación (Habitualmente con puntos equiespaciados) Fórmulas de cuadratura (gausiana) Utilizando las raíces de polinomios ortogonales como puntos de interpolación (Dependen del intervalo y de la función) Método de coeficientes indeterminados Plantean un sistema de ecuaciones, dependiendo del orden de la fórmula, cuya resolución obtiene los coeficientes  Descripción  Objetivos Temario Introducción NC Reglas simples NC Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

6 Ejemplo: Mét. coeficientes indeterminados (I)Integración Numérica Ejemplo: Mét. coeficientes indeterminados (I) Obtener los coeficientes de la siguiente fórmula para que tenga el mayor orden posible  Descripción  Objetivos Temario Introducción NC Reglas simples NC Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía Planteamiento: debe ser exacta para polinomios del mayor orden, luego comenzamos en orden creciente

7 Ejemplo: Mét. coeficientes indeterminados (II)Integración Numérica Ejemplo: Mét. coeficientes indeterminados (II) Resolución del sistema  Descripción  Objetivos Temario Introducción NC Reglas simples NC Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía Orden de la fórmula Fórmula de orden 3

8 Reglas de Newton Cotes SimplesIntegración Numérica Reglas de Newton Cotes Simples Son formulas de tipo interpolatorio que se definen por el número de puntos del polinomio interpolante Las fórmulas genéricas se toman para puntos equiespaciados, aunque su método de obtención es válido para cualesquiera puntos Motivación  Descripción  Objetivos Temario Introducción NC Reglas simples NC Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

9 Ejemplos de reglas simplesIntegración Numérica Ejemplos de reglas simples Interpolante de grado 0 Valor en el extremo del intervalo Valor en el punto medio Interpolante de grado 1  Descripción  Objetivos Temario Introducción NC Reglas simples - Cerradas - Abiertas - Ejemplos NC Reglas Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

10 Ejemplos de reglas simplesIntegración Numérica Ejemplos de reglas simples Interpolante de grado 2 Combinamos trapecio y punto medio  Descripción  Objetivos Temario Introducción NC Reglas simples - Cerradas - Abiertas NC Reglas Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

11 Integración Numérica Reglas cerradas El pol. de interpolación SI incluye los extremos del intervalo  Descripción  Objetivos Temario Introducción NC Reglas simples - Cerradas - Abiertas NC Reglas Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía Tma: Dada una fórmula cerrada de n+1 puntos el error es Fórmulas cerradas usuales n=1 Regla del trapecio n=2 Regla de Simpson n=3 n=4

12 2 puntos (Trapecio) 3 puntos (Simpson) Ejemplo (I)  DescripciónIntegración Numérica Ejemplo (I) 2 puntos (Trapecio) 3 puntos (Simpson)  Descripción  Objetivos Temario Introducción NC Reglas simples - Cerradas - Abiertas NC Reglas Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

13 4 puntos (Simpson 3/8) 5 puntos Ejemplo (II)  Descripción  ObjetivosIntegración Numérica Ejemplo (II) 4 puntos (Simpson 3/8) 5 puntos  Descripción  Objetivos Temario Introducción NC Reglas simples - Cerradas - Abiertas NC Reglas Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

14 Integración Numérica Reglas abiertas El pol. de interpolación NO incluye los extremos del intervalo  Descripción  Objetivos Temario Introducción NC Reglas simples - Cerradas - Abiertas NC Reglas Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía Tma: Dada una fórmula abierta de n+1 puntos el error es Fórmulas abiertas usuales n=0 Regla del punto medio n=1 n=2 n=3

15 1 punto 2 puntos Ejemplo (I)  Descripción  Objetivos TemarioIntegración Numérica Ejemplo (I) 1 punto 2 puntos  Descripción  Objetivos Temario Introducción NC Reglas simples - Cerradas - Abiertas NC Reglas Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

16 3 puntos 4 puntos Ejemplo (II)  Descripción  Objetivos TemarioIntegración Numérica Ejemplo (II) 3 puntos 4 puntos  Descripción  Objetivos Temario Introducción NC Reglas simples - Cerradas - Abiertas NC Reglas Compuestas Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

17 Integración Numérica Reglas Compuestas Subdivisión en subintervalos y aplicación de la regla simple a cada subintervalo Reglas más habituales Regla del trapecio compuesta Regla de Simpson compuesta  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas - Trapecio y Simpson - Control del error - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

18 Ejemplo: trapecio compuestaIntegración Numérica Ejemplo: trapecio compuesta  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas - Trapecio y Simpson - Control del error - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía Regla 2 intervalos (3 puntos) 3 intervalos (4 puntos) 4 intervalos (5 puntos) Cota y error real 2 intervalos 3 intervalos 4 intervalos

19 Ejemplo: simpson compuestaIntegración Numérica Ejemplo: simpson compuesta  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas - Trapecio y Simpson - Control del error - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía Regla 2 intervalos (5 puntos) 3 intervalos (7 puntos) 4 intervalos (9 puntos) Cota y error real 2 int. 3 int. 4 int.

20 Control de error y ejemploIntegración Numérica Control de error y ejemplo Las reglas compuestas permiten controlar el error variando el número de subintervalos (modifican el valor de h). Cotas de error Trapecio Simpson Ejemplo Obtener el numero de subintervalos para poder asegurar que el cálculo de la integral indicada mediante las reglas compuestas tenga un error no superior a una millonésima  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas - Trapecio y Simpson - Control del error - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

21 Ejemplo: simpson compuestaIntegración Numérica Ejemplo: simpson compuesta  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas - Trapecio y Simpson - Control del error - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía Regla 2 intervalos (5 puntos) 3 intervalos (7 puntos) 4 intervalos (9 puntos) Cota y error real 2 int. 3 int. 4 int.

22 Introducción y definiciónIntegración Numérica Introducción y definición No fija los puntos. Valores óptimos para obtener mayor grado Ejemplo  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas Cuadratura Gausiana - Introducción - Polin. Ortogonales - Ejemplos Int. Romberg Int. Adaptativa  Bibliografía Polinomios ortogonales Base de polinomios Definidos en un intervalo [a,b] Ortogonales para una determinada función w(x), denominada peso

23 Propiedades Relación de recurrencia  Descripción  Objetivos TemarioIntegración Numérica Propiedades Relación de recurrencia Orden de las fórmulas y error  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas Cuadratura Gausiana - Introducción - Polin. Ortogonales - Ejemplos Int. Romberg Int. Adaptativa  Bibliografía

24 Familias de polinomios ortogonalesIntegración Numérica Familias de polinomios ortogonales  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas Cuadratura Gausiana - Introducción - Polin. Ortogonales - Ejemplos Int. Romberg Int. Adaptativa  Bibliografía Nombre [a,b] w(x) Notas Legendre [-1,1] 1 Se impone que su valor en x=1 sea 1 Asociados a la ecuación diferencial Tchebishev Uso en interpolación (aproximación minimax) Hermite (-,) Laguerre [0,)

25 Polinomios de LegendreIntegración Numérica Polinomios de Legendre Polinomios Relación de recurrencia Norma Raíces y coeficientes  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas Cuadratura Gausiana - Introducción - Polin. Ortogonales - Ejemplos Int. Romberg Int. Adaptativa  Bibliografía N Raíces Coeff. 1 2 3 4 5

26 Polinomios de TchebishevIntegración Numérica Polinomios de Tchebishev Polinomios Relación de recurrencia Raíces  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas Cuadratura Gausiana - Introducción - Polin. Ortogonales - Ejemplos Int. Romberg Int. Adaptativa  Bibliografía

27 Integración de Romberg: ProcesoIntegración Numérica Integración de Romberg: Proceso Utiliza la regla del trapecio compuesta Se combina con la extrap. de Richardson Aprovecha las relaciones de la regla  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas - Trapecio - Simpson - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

28 Integración de Romberg: RelacionesIntegración Numérica Integración de Romberg: Relaciones  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas - Trapecio - Simpson - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

29 Ejemplo  Descripción  Objetivos Temario  BibliografíaIntegración Numérica Ejemplo Calcular, con un error inferior a 10-6, aplicando el método de integración de Romberg Aplicamos la fórmula de Romberg, donde h=(b-a)/n, para utilizar a continuación la extrapolación de Richardson  Descripción  Objetivos Temario Introducción NC Reglas simples NC Reglas Compuestas - Trapecio - Simpson - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía n h N0(h) N1(h) N2(h) N3(h) N4(h) 1 /2 2 /4 4 /8 8 /16 16 /32

30 Motivación Basada en la fórmula de Simpson compuesta  DescripciónIntegración Numérica Motivación Basada en la fórmula de Simpson compuesta Supone comportamiento “suave” de f(4)(x)  Descripción  Objetivos Temario Introducción Reglas simples Reglas Compuestas - Trapecio - Simpson - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía

31 Proceso Calcular  Descripción Estimar el error como e=|I-Is|Integración Numérica Proceso Calcular Estimar el error como e=|I-Is| Si e<15ε tomar Is como valor de la integral sino repetir el proceso con cada semi-intervalo, dividiendo también por dos el error admisible  Descripción  Objetivos Temario Introducción Reglas simples Reglas Compuestas - Trapecio - Simpson - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía Ejemplo: Calcular, con un error inferior a 1 centésima Representación de la función

32 Ejemplo: resultados  Descripción  Objetivos Temario  BibliografíaIntegración Numérica Ejemplo: resultados  Descripción  Objetivos Temario Introducción Reglas simples Reglas Compuestas - Trapecio - Simpson - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía [a,b] I Ii Id Is |I-Is| tol (4,20) /32 0.5715 0.3812 0.9526 1.2292 0.0100 0.4721 (4,12) /32 0.5715 0.3861 0.1947 0.5808 0.0094 0.0050 (12,20) /32 0.3812 0.8052 (12,16) /32 0.0011 0.0412 0.0025 (16,20) /32 0.0136 0.0646 0.0781 0.3053 (12,14) /32 0.0001 0.0013 (14,16) /32 0.0011 0.0555 0.0032 0.0021 (16,18) /32 0.0136 0.0715 0.0081 0.0055 0.0013 (18,20) /32 0.0646 0.0462 0.0379 0.0266 (18,19) /32 0.0217 0.0005 0.0006 (19,20) /32 0.0462 0.0448 0.0447 0.0015

33 Ejemplo: representación gráficaIntegración Numérica Ejemplo: representación gráfica  Descripción  Objetivos Temario Introducción Reglas simples Reglas Compuestas - Trapecio - Simpson - Estabilidad Cuadratura Gausiana Int. Romberg Int. Adaptativa  Bibliografía Integral: ( )+ ( )+ =