Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 20 vistas

    An Alisis Num Erico (ICM840-01)

    Cargado por

    Grafica
    Este documento presenta un modelo matemático para analizar el movimiento de un paracaidista en descenso. Inicialmente, se define un modelo simplificado mediante una ecuación diferencial que relaciona la aceleración, fuerzas y masa. Luego, se aproxima numéricamente la solución mediante el método de Euler, obteniendo resultados similares a la solución exacta. Finalmente, se discuten conceptos como el análisis de error, cifras significativas y criterios de tolerancia para evaluar la precisión de los métodos numéricos.

    Copyright:

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

    An Alisis Num Erico (ICM840-01)

    Cargado por

    Grafica
    20 vistas22 páginas
    Este documento presenta un modelo matemático para analizar el movimiento de un paracaidista en descenso. Inicialmente, se define un modelo simplificado mediante una ecuación diferencial que relaciona la aceleración, fuerzas y masa. Luego, se aproxima numéricamente la solución mediante el método de Euler, obteniendo resultados similares a la solución exacta. Finalmente, se discuten conceptos como el análisis de error, cifras significativas y criterios de tolerancia para evaluar la precisión de los métodos numéricos.

    Título original

    Untitled

    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?

    Este documento presenta un modelo matemático para analizar el movimiento de un paracaidista en descenso. Inicialmente, se define un modelo simplificado mediante una ecuación diferencial que relaciona la aceleración, fuerzas y masa. Luego, se aproxima numéricamente la solución mediante el método de Euler, obteniendo resultados similares a la solución exacta. Finalmente, se discuten conceptos como el análisis de error, cifras significativas y criterios de tolerancia para evaluar la precisión de los métodos numéricos.

    Copyright:

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

    An Alisis Num Erico (ICM840-01)

    Cargado por

    Grafica
    Este documento presenta un modelo matemático para analizar el movimiento de un paracaidista en descenso. Inicialmente, se define un modelo simplificado mediante una ecuación diferencial que relaciona la aceleración, fuerzas y masa. Luego, se aproxima numéricamente la solución mediante el método de Euler, obteniendo resultados similares a la solución exacta. Finalmente, se discuten conceptos como el análisis de error, cifras significativas y criterios de tolerancia para evaluar la precisión de los métodos numéricos.

    Copyright:

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

    Análisis Numérico (ICM840-01)

    Modelos Matemáticos
    Análisis de Error

    1er Trimestre 2022

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 1 / 21
    Modelos Matemáticos y Análisis de Error
    Modelación matemática en ingenierı́a

    Se define un modelo matemático como una formulación o una ecuación


    que expresa las caracterı́sticas esenciales de un sistema fı́sico o de un
    proceso en términos matemáticos.

    Por lo general, un modelo puede ser representado mediante una relación


    funcional de la forma:

    vardep = f (varindep , parametros, f uncionesf uerza )

    donde:
    vardep es una caracterı́stica que refleja el comportamiento del sistema,
    varindep son las dimensiones que determinan el comportamiento del
    sistema (por lo general, espacio y tiempo),
    parametros son el reflejo de las propiedades del sistema,
    y f uncionesf uerza son las influencias externas que actúan sobre el sistema.
    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 2 / 21
    Modelos Matemáticos
    Ejemplo

    Fuerzas que actúan sobre un paracaidista en descenso.

    Modelo simplificado :

    dv F
    a= = (1)
    dt m
    donde:
    F = FD + FU
    con:
    FD = mg
    FU = −cv

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 3 / 21
    Modelos Matemáticos
    Ejemplo

    Si el paracaidista está inicialmente en reposo (v = 0 en t = 0), entonces la


    solución analı́tica está dada por:
    mg
    v(t) = (1 − e−(c/m)t ) (2)
    c
    A la Ecuación 2 se le llama solución exacta, ya que satisface con exactitud
    la ecuación diferencial original.
    En ingenierı́a, no todas las ecuaciones diferenciales que modelan problemas
    reales pueden resolverse de forma analı́tica.
    Se emplean Métodos Numéricos que permiten aproximar la solución
    exacta.

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 4 / 21
    Modelos Matemáticos
    Ejemplo

    La razón en el cambio de la velocidad con respecto al tiempo se puede


    aproximar de la siguiente forma:
    dv ∆v v(ti+1 ) − v(ti )
    ≃ = (3)
    dt ∆t ti+1 − ti
    A esta Ecuación 3 se le llama aproximación en diferencias finitas
    dividida.
    Sustituyendo la Ecuación 3 en el modelo y reordenando, se tiene:
    c
    v(ti+1 ) = v(ti ) + [g − v(ti )](ti+1 − ti ) (4)
    m
    La ecuación diferencial se ha transformado en la Ecuación 4, algebraica
    fácil de resolver. La aproximación se ha realizado con el denominado
    método de Euler:
    V alorN uevo = V alorAnterior + P endiente · P asoT emporal
    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 5 / 21
    Modelos Matemáticos
    Ejemplo

    Utilizando el software MatLab:


    t=0:2:100; % Vector Tiempo en intérvalos de 2 segundos,
    desde 0 hasta 100
    V=zeros(1,length(t)); % Definir la matriz Velocidad, con
    una fila y la misma cantidad de columnas que el vector
    tiempo.
    for i=1:(length(t)-1) % Ciclo, desde 1 hasta la cantidad de
    elementos de la Velocidad menos uno, como no se especifica
    será de paso uno (por defecto)
    V(i+1)=V(i)+(9.81-12.5/68.1*V(i))*(t(i+1)-t(i)); % Cálculo
    de la Velocidad, elemento i+1, en cada paso del ciclo de i.
    end % Fin del ciclo
    plot(t,V,’b--o’) % Plotear (dibujar) los valores de t en
    las X y de la V en las y, con lı́neas discontinuas y
    cı́rculos
    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 6 / 21
    Modelos Matemáticos
    Ejemplo

    Al comparar gráficamente la solución exacta con la solución analı́tica, se


    tiene:

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 7 / 21
    Modelos Matemáticos
    Comentarios generales

    • El método numérico aproxima de buena forma a la solución exacta.


    • La naturaleza del método numérico (pendiente que aproxima una
    curva) genera diferencias en los resultados.
    • ¿Cómo reducir esta diferencia?
    • Usando un menor paso temporal.
    • ¿Cual es el costo?
    • Mayor cantidad de cálculos.
    • Entonces, ¿cuál es el paso de tiempo óptimo para asegurar un buen
    resultado?

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 8 / 21
    Análisis de Error
    Exactitud y Precisión

    Exactitud
    Se refiere a qué tan cercano está el valor calculado del valor verdadero.

    Precisión
    Se refiere a qué tan cercanos se encuentran, unos de otros, los valores
    calculados.

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 9 / 21
    Análisis de Error
    Cifras Significativas

    • Cuando se emplea un número para realizar un cálculo, debe haber


    seguridad de que puede usarse con confianza. El concepto de cifras o
    dı́gitos significativos ha sido desarrollado para designar formalmente la
    confiabilidad de un valor numérico.
    • Las cifras significativas (CS) de un número son aquellas que pueden
    usarse de manera confiable.
    √
    • Existen ciertas cantidades especı́ficas (π, e, 7, etc) que no se
    pueden expresar con un número finito de dı́gitos. (Ejemplo:
    π = 3,141592653589793238462643...)
    • Este tipo de números no puede ser representado con exactitud por un
    computador. Solo retendrán una cantidad determinada de CS. A la
    omisión del resto de CS se le conoce como error de redondeo.
    • Los errores de truncamiento representan la diferencia entre la solución
    exacta de un problema y su aproximación numérica.
    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 10 / 21
    Análisis de Error
    Definición de Error

    • Los errores numéricos surgen del uso de aproximaciones para


    representar operaciones y cantidades matemáticas exactas.
    • Se define el error absoluto como:

    Eabs = V alor verdadero − V alor aproximado


    • Este valor no permite visualizar el orden de magnitud del error.
    • Para tomar en cuenta las magnitudes de las cantidades que se
    evalúan, se debe normalizar el error respecto del valor verdadero. Esto
    se conoce como error relativo:

    V alor verdadero − V alor aproximado


    Erel = 100 %
    V alor verdadero

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 11 / 21
    Análisis de Error
    Definición de Error

    • En el uso de Métodos Numéricos, el valor verdadero solo se conocerá


    cuando el problema tenga solución analı́tica.
    • ¿Como estimamos entonces el error si nuestro problema no tiene
    solución analı́tica?

    Aproximacion actual − Aproximacion anterior


    Eaprox = 100 %
    Aproximacion actual
    • ¿Que criterio utilizamos para llegar a un valor de error aceptable?

    | Eaprox |< Es

    en donde Es corresponde a una tolerancia porcentual definida


    previamente.

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 12 / 21
    Análisis de Error
    Criterio de Tolerancia

    • Se recomienda relacionar el error con el número de cifras significativas


    en la aproximación.
    • Se ha demostrado (Scarborough, 1966) que si se cumple que:

    Es = (0,5 · 102−n ) %

    donde n es el numéro de cifras significativas.


    • Entonces se tiene la seguridad que el resultado es correcto en al
    menos n cifras significativas.

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 13 / 21
    Análisis de Error
    Ejemplo 1

    • La función ex puede ser aproximada mediante la Serie de MacLaurin


    de la forma:

    x2 x3 xn
    ex = 1 + x + + + ... +
    2! 3! n!
    Estime el valor de e0,5 considerando Es = 0,05 %. Considere que el
    valor verdadero es e0,5 = 1,648721.

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 14 / 21
    Análisis de Error
    Ejemplo 1

    • La función ex puede ser aproximada mediante la Serie de MacLaurin


    de la forma:

    x2 x3 xn
    ex = 1 + x + + + ... +
    2! 3! n!
    Estime el valor de e0,5 considerando Es = 0,05 %. Considere que el
    valor verdadero es e0,5 = 1,648721.
    • Solución:
    • T erminos = 6
    • Resultado = 1,648697917
    • Erel = 0,00142 % < 0,05 %
    • Eaprox = 0,0158 %

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 14 / 21
    Análisis de Error
    Ejercicio propuesto 1

    • Evalúe e−5 con el uso de dos métodos:


    • Método 1:

    x2 x3
    e−x = 1 − x + − + ...
    2! 3!
    • Método 2:

    1 1
    e−x = = x2 x3
    ex 1 + x + 2! + 3! + ...
    • Compare con el valor verdadero de 6,737947 · 10−3 . Utilice 10
    términos para evaluar cada serie. Calcule los errores.

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 15 / 21
    Análisis de Error
    Ejercicio propuesto 2

    • Determine el número de términos necesarios para aproximar cos(x) a


    8 cifras significativas con el uso de la serie de MacLaurin:

    x2 x4 x6 x8
    cos(x) = 1 −
    + − + − ...
    2! 4! 6! 8!
    • Calcule la aproximación con el empleo del valor de x = 0,3π.

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 16 / 21
    Análisis de Error
    Propagación del error

    • Si xf l está referido a la representación de punto flotante


    (computacional) del número real x
    • Debido a las limitaciones computacionales de representación, existirá
    un error de redondeo entre xf l y x. La estimación del cálculo de f (x):
    ∆f (xf l ) = |f (x) − f (xf l )|
    Usualmente x y f (x) son desconocidos, pero suponiendo que xf l es
    cercano a x se pudiera usar la 1ra expansión de la Serie de Taylor:

    f (x) − f (xf l ) ∼
    = f ′ (xf l )(x − xf l ) → ∆f (xf l ) = f ′ (xf l ) · ∆x

    Lo cual puede ser útil cuando f ′ (xf l ) y ∆x son conocidos.


    Ver material extra sobre Residuo en la Expansión de la Serie de
    Taylor.

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 17 / 21
    Análisis de Error
    Ejemplo 2

    Utilice expanciones de la Serie de Taylor con n desde 0 hasta 6 para


    aproximar f (x) = cos(x) en xi+1 = π/3 con base en el valor de f (x) y sus
    derivadas en xi = π/4
    Definiendo el paso: h = π/3 − π/4
    Aproximación de orden cero: f (x+1 ) ∼
    = f (xi )
    f (π/3) ∼
    = cos(π/4) = 0,707106781

    En este simple ejemplo se conoce el valor real exacto, situación rara en los
    problemas reales: f (π/3) = 0,5

    Error Relativo Porcentual Real:


    0,5 − 0,707106781
    εt = · 100 % = −41,4 %
    0,5
    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 18 / 21
    Análisis de Error
    Ejemplo 2

    Aproximación de primer orden: f (xi+1 ) ∼


    = f (xi ) + f ′ (xi )(xi+1 − xi )

    f (π/3) ∼
    = cos(π/4) − sen(π/4) · (π/12) = 0,521986659

    Error Relativo Porcentual Real: εt = −4,40 %


    Aproximación de segundo orden:

    f ′′ (xi )
    f (xi+1 ) ∼
    = f (xi ) + f ′ (xi )(xi+1 − xi ) + (xi+1 − xi )2
    2!
    cos(π/4)
    f (π/3) ∼
    = cos(π/4) − sen(π/4) · (π/12) + (π/12)2 = 0,497754491
    2
    Error Relativo Porcentual Real: εt = 0,449 %

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 19 / 21
    Análisis de Error
    Ejemplo 2

    Notar que la mayor mejorı́a del error es con los cuatro primeros órdenes de
    expansión de la serie.

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 20 / 21
    Análisis Numérico (ICM840-01)
    Modelos Matemáticos
    Análisis de Error

    1er Trimestre 2022

    Profesor: Angel Rodrı́guez Soto (PUCV) Análisis Numérico 1er Trimestre 2022 21 / 21

    También podría gustarte