Integracin numrica Regla de Simpson 13 simple El

  • Slides: 13
Download presentation
Integración numérica Regla de Simpson 1/3 simple El método consiste en usar polinomios de

Integración numérica Regla de Simpson 1/3 simple El método consiste en usar polinomios de grado superior para aproximar la curva de la función y tomar las integrales bajo tales polinomios. Polinomio de lagrange P 2(x) de segundo grado como una aproximación de f(x) Fórmulas

Integración numérica Regla de Simpson 1/3 Simple Utilizar la regla de Simpson 1/3 simple

Integración numérica Regla de Simpson 1/3 Simple Utilizar la regla de Simpson 1/3 simple para aproximar la integral Fórmulas

Integración numérica Regla de Simpson 1/3 compuesta En el caso de que el intervalo

Integración numérica Regla de Simpson 1/3 compuesta En el caso de que el intervalo [a, b] no sea lo suficientemente pequeño, el error al calcular la integral puede ser muy grande. Para ello, se recurre a la fórmula compuesta de Simpson 1/3. Fórmulas

Integración numérica Regla de Simpson 1/3 compuesta Utilizar la regla de Simpson 1/3 compuesta

Integración numérica Regla de Simpson 1/3 compuesta Utilizar la regla de Simpson 1/3 compuesta con n=4 Subintervalos para aproximar la integral Hojas de cálculo

ECUACIONES DIFERENCIALES MÉTODO DE EULER Una de las técnicas más simples para aproximar soluciones

ECUACIONES DIFERENCIALES MÉTODO DE EULER Una de las técnicas más simples para aproximar soluciones de una ecuación diferencial es el método de Euler, o de las rectas tangentes. La idea del método de Euler es muy sencilla y está basada en el significado geométrico de la derivada de una función en un punto dado. Condición inicial Fórmulas y(0) = x 0 Xn+1=xn+h yn+1 = yn + f(xn, yn)h

MÉTODO DE EULER (Ejercicio) Dada la siguiente ecuación diferencial 2 xy = dy/dx con

MÉTODO DE EULER (Ejercicio) Dada la siguiente ecuación diferencial 2 xy = dy/dx con la condición inicial: y(0) = 1 (a) Resuélvala analíticamente. (b) Utilice el método de Euler en el intervalo de x = 0 a 1 con tamaño de paso 0. 1 (c) Aproxime y(0. 5) (d) Grafique y compare los resultados X =x +h n+1 n yn+1 = yn + f(xn, yn)h (a) Resuélvala analíticamente. y = ex 2 (b) Utilice el método de Euler en el intervalo de x = 0 a 1 con tamaño de paso 0. 1 Condición inicial: x = 0 → y = 1 , (x 0 , y 0 ) = ( 0 ; 1 ) n=0 x 1 = x 0 + h = 0 + 0, 1 = 0, 1 y 1 = y 0 + f (x 0 , y 0 ). h= 1+ 2(0)(1) ⋅ 0, 1 = 1 n=1 (x 1 , y 1 ) = ( 0, 1 ; 1 ) x 2 = x 1 + h = 0, 1 + 01 = 0, 2 y 2 = y 1 + f (x 1 , y 1 ). h= 1+ 2(0, 1)(1) ⋅ 0, 1 = 1, 02 i 0 1 2 3 4 5 6 7 8 9 10 xi 0 0, 1 0, 2 0, 3 0, 4 0, 5 0, 6 0, 7 0, 8 0, 9 1 yi 1 1 1, 02 1, 0608 1, 124448 1, 21440384 1, 33584422 1, 49614553 1, 70560591 1, 97850285 2, 33463336

MÉTODO DE EULER (Ejercicio) Dada la siguiente ecuación diferencial 2 xy = dy/dx con

MÉTODO DE EULER (Ejercicio) Dada la siguiente ecuación diferencial 2 xy = dy/dx con la condición inicial: y(0) = 1 (a) Resuélvala analíticamente. (b) Utilice el método de Euler en el intervalo de x = 0 a 1 con tamaño de paso 0. 1 (c) Aproxime y(0. 5) (d) Grafique y compare los resultados X =x +h n+1 n yn+1 = yn + f(xn, yn)h (c) Aproxime y(0. 5) El error relativo porcentual que se cometió al aplicar la formula de Euler esta dado por: Er =5. 423% (d) Grafique y compare los resultados y 3 2. 5 2 1. 5 1 0. 5 0 0 0. 1 0. 2 0. 3 0. 4 0. 5 0. 6 0. 7 0. 8 0. 9 Solución numérica Solución analítica 1 xi

MÉTODO DE EULER MODIFICADO Este método se basa en la misma idea del método

MÉTODO DE EULER MODIFICADO Este método se basa en la misma idea del método anterior, pero hace un refinamiento en la aproximación, tomando un promedio entre ciertas pendientes. Esto permite obtener una mejor aproximación de la pendiente en todo el intervalo. Fórmulas Xn+1=xn+h y*n+1 = yn + f(xn, yn)h yn+1 = yn + ((f(xn, yn)+f(xn+1, y*n+1))/2)h

MÉTODO DE EULER MODIFICADO Dada la siguiente ecuación diferencial 2 xy = dy/dx con

MÉTODO DE EULER MODIFICADO Dada la siguiente ecuación diferencial 2 xy = dy/dx con la condición inicial: y(0) = 1 (a) Resuélvala analíticamente. (b) Utilice el método de Euler modificado en el intervalo de x = 0 a 1 con tamaño de paso 0. 1 (c) Aproxime y(0. 5) (d) Grafique y compare los resultados X =x +h n+1 (a) Resuélvala analíticamente. n y*n+1 = yn + f(xn, yn)h yn+1 = yn + ((f(xn, yn)+f(xn+1, y*n+1))/2)h y = ex 2 (b) Utilice el método de Euler en el intervalo de x = 0 a 1 con tamaño de paso 0. 1 n=0 Condición inicial: x = 0 → y = 1 , (x 0 , y 0 ) = ( 0 ; 1 ) i 0 1 2 3 4 5 6 7 8 9 10 xi 0 0, 1 0, 2 0, 3 0, 4 0, 5 0, 6 0, 7 0, 8 0, 9 1 yi 1 1, 040704 1, 093988 1, 173192 1, 283472 1, 432355 1, 630593 1, 893445 2, 242596 2, 709057

MÉTODO DE EULER MODIFICADO Dada la siguiente ecuación diferencial 2 xy = dy/dx con

MÉTODO DE EULER MODIFICADO Dada la siguiente ecuación diferencial 2 xy = dy/dx con la condición inicial: y(0) = 1 (a) Resuélvala analíticamente. (b) Utilice el método de Euler modificado en el intervalo de x = 0 a 1 con tamaño de paso 0. 1 (c) Aproxime y(0. 5) (d) Grafique y compare los resultados Xn+1=xn+h y*n+1 = yn + f(xn, yn)h yn+1 = yn + ((f(xn, yn)+f(xn+1, y*n+1))/2)h (c) Aproxime y(0. 5) El error relativo porcentual que se cometió al aplicar la formula de Euler modificado esta dado por: Er =0. 044% y 3 2. 5 2 1. 5 1 0. 5 0 0 0. 1 0. 2 0. 3 0. 4 0. 5 0. 6 0. 7 0. 8 0. 9 Solución numérica Solución analítica 1 xi

MÉTODO DE RUNGE –KUTTA DE CUARTO ORDEN Los Runge-Kutta no es solo un método

MÉTODO DE RUNGE –KUTTA DE CUARTO ORDEN Los Runge-Kutta no es solo un método sino una importante familia de métodos iterativos tanto implícitos como explícitios para aproximar las soluciones de ecuaciones diferenciales ordinarias. Fórmulas

MÉTODO DE RUNGE –KUTTA DE CUARTO ORDEN Continuación Fórmulas

MÉTODO DE RUNGE –KUTTA DE CUARTO ORDEN Continuación Fórmulas

MÉTODO DE RUNGE –KUTTA DE CUARTO ORDEN El error relativo porcentual que se cometió

MÉTODO DE RUNGE –KUTTA DE CUARTO ORDEN El error relativo porcentual que se cometió al aplicar la Runge-Kutta esta dado por: Er = 0. 000779% y 3 2. 5 2 1. 5 Hojas de cálculo 1 0. 5 0 0 0. 1 0. 2 0. 3 0. 4 0. 5 0. 6 0. 7 0. 8 0. 9 1 xi Solución numérica Solución analítica