Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 19 vistas

    Entrega2 Resolucion

    Cargado por

    Ruben
    1) El método de Euler para este problema es decreciente. El error local está acotado por 2h y el error global por (e^2 - 1)h. 2) Para que el error al aproximar y(1) sea menor que 10^-2, el paso h debe ser menor que 0,001565176.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd

    Entrega2 Resolucion

    Cargado por

    Ruben
    19 vistas3 páginas
    1) El método de Euler para este problema es decreciente. El error local está acotado por 2h y el error global por (e^2 - 1)h. 2) Para que el error al aproximar y(1) sea menor que 10^-2, el paso h debe ser menor que 0,001565176.

    Descripción original:

    Calculo Numerico

    Título original

    Entrega2_Resolucion

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    1) El método de Euler para este problema es decreciente. El error local está acotado por 2h y el error global por (e^2 - 1)h. 2) Para que el error al aproximar y(1) sea menor que 10^-2, el paso h debe ser menor que 0,001565176.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    19 vistas3 páginas

    Entrega2 Resolucion

    Cargado por

    Ruben
    1) El método de Euler para este problema es decreciente. El error local está acotado por 2h y el error global por (e^2 - 1)h. 2) Para que el error al aproximar y(1) sea menor que 10^-2, el paso h debe ser menor que 0,001565176.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    de 3

    Elementos de Cálculo Numérico - Cálculo Numérico

    Segundo Cuatrimestre de 2020


    Entrega n◦ 2

     2
    y 0 (t) = −e−t (y(t))2
    1. Dado el problema: .
    y(0) = 1

    a) Expresar la iteración del método de Euler correspondiente a este problema y


    mostrar que es decreciente.
    b) Mostrar que y es decreciente y, si se sabe que es siempre positiva, estimar el error
    de truncado local para t ∈ [0, 1].
    c) Si se sabe que las soluciones numérica y exacta son siempre positivas, hallar una
    cota del paso h para que el error al aproximar y(1) sea menor que 10−2 .
    Elementos de Cálculo Numérico - Cálculo Numérico
    Segundo Cuatrimestre de 2020
    Entrega n◦ 2 - Resolución del ejercicio

    1a) La iteración del método de Euler con paso h es:


    ( 2
    yi+1 = yi − he−ti yi2 , para 0 ≤ i ≤ N − 1,
    y0 = 1.
    2
    Como he−ti yi2 ≥ 0 para todo 0 ≤ i ≤ N − 1, tenemos que yi+1 ≤ yi para todo
    0 ≤ i ≤ N − 1, y por lo tanto es decreciente.

    2
    1b) Primero notemos que −e−t (y(t))2 < 0 para todo t. Por lo que y 0 < 0 para todo t e y es
    decreciente. Ahora, el error de truncado local para el método de Euler es τi = h2 y 00 (ξi ),
    2
    para ξi ∈ [ti , ti+1 ], 0 ≤ i ≤ N − 1. Como y 00 (t) = dtd f (t, y(t)) = dtd (−e−t (y(t))2 ) =
    2 2
    2e−t (y(t))2 [t + e−t y(t)], tenemos:

    h 00 2 2
    τ= y (ξ) = he−ξ (y(ξ))2 [ξ + e−ξ y(ξ)].
    2
    Para acotar esta expresión (en módulo) debemos considerar que ξ ∈ [0, 1], y que
    2
    0 < y(ξ) ≤ y(0) = 1 y e−ξ ≤ e0 = 1, para todo ξ ∈ [0, 1]. Luego:

    |τ | ≤ h[1 + 1] = 2h.

    De esta forma, el valor máximo que puede tomar |τ |, τM AX , se puede acotar por
    τM AX ≤ 2h.

    1c) Para encontrar un valor del paso h para que el error cometido al aproximar y(1) sea
    menor que 10−2 , debemos trabajar con el error global del método. Esto es, |eN | =
    |y(tf ) − yN |, donde tf es el tiempo final y N = (tf − t0 )/h para t0 el tiempo inicial.
    Sabemos que
    τM AX K(tf −t0 )
    |eN | ≤ (e − 1),
    K
    donde, para el método de Euler, K es tal que es independiente de t y de h y vale que
    |f (t, y) − f (t, z)| ≤ K|y − z|.
    Veamos:
    2 2 2
    |f (t, y) − f (t, z)| = | − e−t y 2 − (−e−t z 2 )| = e−t |y 2 − z 2 |.
    2
    Aquı́, como t ∈ [0, 1], tenemos e−t ≤ 1. Para acotar |y 2 − z 2 | podemos usar el Teorema
    del Valor Medio, para ϕ(x) = x2 y nos queda |y 2 − z 2 | = |2ξ||y − z|:

    |f (t, y) − f (t, z)| ≤ 2|ξ||y − z|,


    para ξ entre y y z. Si miramos la demostración del Teorema 8.10 del apunte de la
    materia, notaremos que los roles de y y z son los de las soluciones exacta y numérica,
    respectivamente. Como sabemos que cumplen 0 ≤ y(t) ≤ y(0) = 1 para todo t ∈ [0, 1]
    y 0 ≤ yi ≤ y0 = 1 para todo 0 ≤ i ≤ N , tendremos ξ ∈ [0, 1], y por lo tanto podemos
    tomar K = 2.
    También podrı́amos haber acotado |y 2 −z 2 | = |y +z||y −z| ≤ (|y|+|z|)|y −z| ≤ 2|y −z|.
    De esta forma, nos queda:
    τM AX K(1−0) 2h 2
    |eN | ≤ (e − 1) ≤ (e − 1) = (e2 − 1)h.
    K 2

    Ahora busquemos un valor de h para que este error sea menor que 10−2 :

    1
    (e2 − 1)h < 10−2 ⇔ h < ∼ 0,001565176.
    100(e2 − 1)
    Por lo tanto, cualquier h < 0,001565176 sirve para tal propósito.

    También podría gustarte