Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 148 vistas

    MB536 Metodos Numericos

    Cargado por

    Anonymous 5jQw8JS
    El documento presenta un examen final de métodos numéricos que consta de 4 problemas. El primer problema involucra encontrar raíces de un sistema de ecuaciones no lineales usando métodos iterativos como Newton-Raphson y punto fijo. El segundo problema trata de ajustar datos experimentales a una función usando mínimos cuadrados. El tercer problema aproxima una integral definida usando diferentes métodos numéricos. El cuarto problema calcula la deflexión de un mástil usando métodos de Euler y Runge-Kutta.

    Copyright:

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

    MB536 Metodos Numericos

    Cargado por

    Anonymous 5jQw8JS
    148 vistas6 páginas
    El documento presenta un examen final de métodos numéricos que consta de 4 problemas. El primer problema involucra encontrar raíces de un sistema de ecuaciones no lineales usando métodos iterativos como Newton-Raphson y punto fijo. El segundo problema trata de ajustar datos experimentales a una función usando mínimos cuadrados. El tercer problema aproxima una integral definida usando diferentes métodos numéricos. El cuarto problema calcula la deflexión de un mástil usando métodos de Euler y Runge-Kutta.

    Descripción original:

    metodos numericos

    Título original

    MB536 metodos numericos

    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?

    El documento presenta un examen final de métodos numéricos que consta de 4 problemas. El primer problema involucra encontrar raíces de un sistema de ecuaciones no lineales usando métodos iterativos como Newton-Raphson y punto fijo. El segundo problema trata de ajustar datos experimentales a una función usando mínimos cuadrados. El tercer problema aproxima una integral definida usando diferentes métodos numéricos. El cuarto problema calcula la deflexión de un mástil usando métodos de Euler y Runge-Kutta.

    Copyright:

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

    MB536 Metodos Numericos

    Cargado por

    Anonymous 5jQw8JS
    El documento presenta un examen final de métodos numéricos que consta de 4 problemas. El primer problema involucra encontrar raíces de un sistema de ecuaciones no lineales usando métodos iterativos como Newton-Raphson y punto fijo. El segundo problema trata de ajustar datos experimentales a una función usando mínimos cuadrados. El tercer problema aproxima una integral definida usando diferentes métodos numéricos. El cuarto problema calcula la deflexión de un mástil usando métodos de Euler y Runge-Kutta.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    Está en la página 1de 6

    UNIVERSIDAD NACIONAL DE INGENIERIA P.A.

    2018-3
    FACULTAD DE INGENIERIA MECANICA 27/02/2019
    DACBHCC 18:00-19:50 h

    EXAMEN FINAL DE METODOS NUMERICOS (MB536)

     SOLO SE PERMITE EL USO DE UNA HOJA DE FORMULARIO Y CALCULADORA


     ESCRIBA CLARAMENTE SUS PROCEDIMIENTOS
     PROHIBIDO EL USO DE CELULARES U OTROS EQUIPOS DE COMUNICACION ELECTRONICA
     DURACION: 110 MINUTOS

    Problema 1
    y 2 + x 2 + 5 = 8x
    Sea el sistema de ecuaciones no lineales:
    xy 2 + x + 2 = 8y

    a) (1 P) Bosqueje a mano alzada y localice las raíces, indicando valores aproximados al


    entero más cercano a la raíz.
    b) (1.5 P) Determine la raíz más alejada del origen de coordenadas realizando 02
    iteraciones paso a paso usando el algoritmo de Newton-Raphson para sistemas, a
    partir del valor obtenido en a), muestre el error usando norma infinita.
    c) (1.5 P) Encuentre una fórmula del algoritmo de punto fijo convergente aplicando el
    criterio de convergencia para la raíz buscada en b).
    d) (1 P) Realice tres iteraciones de la fórmula encontrada en c) a partir del valor inicial
    dado en a) y muestre el error (el cual debe ser decreciente).

    Problema 2
    Un móvil se mueve verticalmente hacia arriba de acuerdo a la siguiente data:

    t (seg) 0 0.2 0.4 0.6 0.8 1


    y (m) 1 1.4 1.6 1.8 2.1 2.3

    Se desea ajustar estos datos a la función 𝑦 = 𝐿𝑛(𝑎𝑡 + 𝑏𝑒 𝑡 )


    (a) (2 P) Hallar a y b óptimos para un ajuste por mínimos cuadrados
    (b) (1 P) Determine el factor de regresión y comente sus resultados
    (c) (1 P) Estime la velocidad en t=0.4 s, usando la fórmula de diferenciación central de 3
    puntos de la tabla y compárela con la derivada de la función de ajuste. ¿Cuál considera
    más precisa?
    (d) (1 P) Escriba un programa en MATLAB para resolver la parte (b)

    Problema 3
    Sea la integral con límite infinito:
    ∞
    ∫ 𝑒 −𝑥 √𝑥 𝑑𝑥
    1
    a) (1.5 P) Aproxime la integral reemplazando el limite infinito por un valor b=11
    aplique el Método de Simpson 1/3, con h=1.
    b) (1.5 P) Haga el cambio de variable x=1/t y aproxime la nueva integral con
    limites finitos usando la regla del rectángulo dividiendo el intervalo en 8
    particiones iguales.
    c) (1.5 P) Resuelva la integral dada en b) aplicando la Cuadratura de Gauss con
    N=3.
    d) (0.5 P) Evalué el error para cada caso el valor exacto es 0.4916517 y comente
    sus resultados.
    UNIVERSIDAD NACIONAL DE INGENIERIA P.A. 2018-3
    FACULTAD DE INGENIERIA MECANICA 27/02/2019
    DACBHCC 18:00-19:50 h

    Problema 4
    La ecuación siguiente se utiliza para modelar la deflexión del mástil de un bote sujeto a la
    fuerza del viento:
    𝑑2 𝑦 𝑓(𝑧)
    = (𝐿 − 𝑧)2
    𝑑𝑧 2 2𝐸𝐼
    Donde E=1.25x108, es el módulo de Elasticidad, I=0.05, es el momento de inercia, L=30, es
    la longitud del mástil además se sabe que la fuerza del viento varía con la altitud de acuerdo
    a la relación:
    200𝑧 −𝑧/15
    𝑓(𝑧) = 𝑒
    5+𝑧
    Calcule la deflexión para z=1.5 si para z=0, y=0 y dy/dz=0:
    (a) (2 P) Usando Euler con un paso de h=0.5
    (b) (2 P) Usando Runge-Kutta 2 con un paso de h=0.5
    (c) (1 P) Escriba un programa en Matlab para la solución de b)

    Ing. Robert Castro S.


    UNIVERSIDAD NACIONAL DE INGENIERIA P.A. 2018-3
    FACULTAD DE INGENIERIA MECANICA 27/02/2019
    DACBHCC 18:00-19:50 h

    Solucionario

    Problema 1
    a)

    Raíces aproximadas (1,0) y (2,3)


    b) Aplicando Newton-Raphson para la mayor raíz:
    𝒙𝟎 = 𝟐 𝒚𝟎 = 𝟑
    𝟐𝒙𝟎 − 𝟖 𝟐𝒚𝟎 ∆𝒙 𝒙𝟐 − 𝟖𝒙𝟎 + 𝒚𝟐𝟎 + 𝟓
    [ 𝟐 ] [ ] = − [ 𝟎𝟐 ]
    𝒚𝟎 + 𝟏 𝟐𝒙𝟎 𝒚𝟎 − 𝟖 ∆𝒚 𝒙𝟎 𝒚𝟎 + 𝒙𝟎 − 𝟖𝒚𝟎 + 𝟐
    −𝟒 𝟔 ∆𝒙 −𝟐
    [ ][ ] = [ ]
    𝟏𝟎 −𝟏 ∆𝒚 𝟐
    ∆𝒙 𝟎. 𝟐𝟔𝟑𝟐
    [ ]=[ ] 𝒆𝒓𝒓 = 𝟎. 𝟐𝟔𝟑𝟐
    ∆𝒚 −𝟎. 𝟏𝟓𝟕𝟗
    𝒙𝟏 = 𝒙𝟎 + ∆𝒙 = 𝟐. 𝟐𝟔𝟑𝟐 𝒚𝟏 = 𝒚𝟎 + ∆𝒚 = 𝟐. 𝟖𝟒𝟐𝟏
    ∆𝒙 𝟎. 𝟎𝟐𝟐𝟕
    [ ]=[ ] 𝒆𝒓𝒓 = 𝟎. 𝟎𝟐𝟐𝟕
    ∆𝒚 −𝟎. 𝟎𝟎𝟐𝟕
    𝒙𝟐 = 𝒙𝟏 + ∆𝒙 = 𝟐. 𝟐𝟖𝟓𝟗 𝒚𝟐 = 𝒚𝟏 + ∆𝒚 = 𝟐. 𝟖𝟑𝟗𝟒
    c)
    𝟏
    𝒚 = (𝟖𝒙 − 𝒙𝟐 − 𝟓)𝟐 =𝑮𝟏 (𝒙, 𝒚)
    𝟖𝒚 − 𝟐
    𝒙= 𝟐 = 𝑮𝟐 (𝒙, 𝒚)
    𝒚 +𝟏
    𝟎. 𝟕𝟓𝟓𝟗 𝟎
    𝑱(𝒙𝟎 , 𝒚𝟎 ) = [ ]
    𝟎 −𝟎. 𝟓𝟐𝟎𝟎
    ‖𝑱(𝒙𝟎 , 𝒚𝟎 )‖ = 𝟎. 𝟕𝟓𝟓𝟗<1
    ∞
    Por lo tanto es convergente
    d)
    Realizando iteraciones:
    𝟏
    𝒚 = (𝟖𝒙 − 𝒙𝟐 − 𝟓)𝟐
    𝟖𝒚 − 𝟐
    𝒙= 𝟐
    𝒚 +𝟏
    X Y ERR
    2 3 -----
    2.2 2.6548 0.3542
    2.3958 2.7857 0.1958
    2.3157 2.9028 0.1171
    UNIVERSIDAD NACIONAL DE INGENIERIA P.A. 2018-3
    FACULTAD DE INGENIERIA MECANICA 27/02/2019
    DACBHCC 18:00-19:50 h

    Problema 2
    a)

    i 1 2 3 4 5 6
    ti (seg) 0 0.2 0.4 0.6 0.8 1
    yi (m) 1 1.4 1.6 1.8 2.1 2.3
    Yi=eyi 2.7183 4.0552 4.9530 6.0496 8.1662 9.9742
    ̂𝑖 2.7557 3.8387 5.0568 6.4399 8.0245 9.8552
    𝑌
    𝑦 = 𝐿𝑛(𝑎𝑡 + 𝑏𝑒 𝑡 )
    𝑌 = 𝑒 𝑦 = 𝑎𝑡 + 𝑏𝑒 𝑡
    Aplicando la ecuación normal:
    𝒕𝟏 𝒆𝒕𝟏
    𝒕𝟐
    𝒆𝒚𝟏
    𝒕𝟐 𝒆 𝒆𝒚𝟐
    𝒕𝟑 𝒆𝒕𝟑 𝒂 𝒆𝒚𝟑
    𝒕𝟒 [𝒃] = 𝒆𝒚𝟒
    𝒕𝟒 𝒆
    𝒕𝟓 𝒆𝒕𝟓 𝒆𝒚𝟓
    [𝒆𝒚𝟔 ]
    [𝒕𝟔 𝒆𝒕𝟔 ]
    𝑴𝒄 = 𝑵
    𝑴 𝑴𝒄 = 𝑴𝑻 𝑵
    𝑻

    𝟐. 𝟐 𝟔. 𝟒𝟑𝟑 𝒂 𝟐𝟐. 𝟗𝟐𝟗𝟐 𝒂 𝟐. 𝟑𝟔𝟒𝟓


    [ ][ ] = [ ] [ ]=[ ]
    𝟔. 𝟒𝟑𝟑 𝟐𝟎. 𝟑𝟕𝟗𝟔 𝒃 𝟕𝟏. 𝟑𝟕𝟎𝟑 𝒃 𝟐. 𝟕𝟓𝟓𝟕
    b) Calculo del factor de regresión
    ̂𝑖 = 𝑎𝑡𝑖 + 𝑏𝑒 𝑡𝑖
    𝑌
    6
    𝑌𝑖
    𝑌𝑚 = ∑ = 5.9861
    6
    𝑖=1
    ∑ ̂𝑖 − 𝑌𝑚)𝟐
    (𝑌
    𝑹𝟐 = = 𝟎. 𝟗𝟕𝟓𝟑
    ∑ (𝑌𝑖 − 𝑌𝑚)𝟐
    c)
    𝒚(𝟎. 𝟔) − 𝒚(𝟎. 𝟐)
    𝒗(𝟎. 𝟒) = =𝟏
    𝟐𝒙𝟎. 𝟐

    𝑎 + 𝑏𝑒 𝑡
    𝑣(𝑡) = => 𝑣(0.4) = 1.2806
    𝑎𝑡 + 𝑏𝑒 𝑡
    La primera fórmula es más precisa debido a que usa datos de la tabla directamente.
    d)
    x=0:0.2:1, y=[ 1 1.4 1.6 1.8 2.1 2.3]
    Y=exp(y)'
    P=[x' exp(x)']
    M=P'*P
    N=P'*Y
    c=M\N
    a=c(1), b=c(2)
    Ys=a*x+b*exp(x)
    Ym=mean(Y)
    R2=sum((Ys-Ym).^2)/sum((Y-Ym).^2)
    UNIVERSIDAD NACIONAL DE INGENIERIA P.A. 2018-3
    FACULTAD DE INGENIERIA MECANICA 27/02/2019
    DACBHCC 18:00-19:50 h

    Problema 3

    a) Simpson 1/3
    11
    ∫ 𝑒 −𝑥 √𝑥 𝑑𝑥
    1
    h=1
    I1=h/3*(f(1)+4*f(2)+2*f(3)+ 4*f(4)+2*f(5)+ 4*f(6)+2*f(7)+ 4*f(8)+2*f(9)+ 4*f(10)+f(11))
    I1=0.50561813

    b) Regla del rectangulo:


    1
    c) ∫0 𝑒 −1/𝑡 √1/𝑡 /𝑡 2 𝑑𝑡
    h=1/8
    I2=2*h*(F(1/8)+ F(3/8)+ F(5/8)+ F(7/8))=0.49166517

    c) Cuadratura de Gauss:
    𝑠+1
    𝑡= 𝑑𝑡 = 𝑑𝑠/2
    2
    1 −2
    2 𝑑𝑠
    ∫ 𝑒 𝑠+1 √ /(𝑠 + 1)2
    0 𝑠+1 2
    5 3 8 5 3
    𝐼3 = 9 ∗ 𝐺 (−√5) + 9 ∗ 𝐺(0) + 9 ∗ 𝐺 (√5) =0.47073977

    d) Errores:
    Ie=0.50728223
    Ea=0.00166410
    Eb=0.01561706
    Ec=0.03654246

    Se observa que la regla de Simpson resultó más precisa.


    UNIVERSIDAD NACIONAL DE INGENIERIA P.A. 2018-3
    FACULTAD DE INGENIERIA MECANICA 27/02/2019
    DACBHCC 18:00-19:50 h

    Problema 4
    a) Euler
    16
    𝑑2 𝑦 106 𝑧 −𝑧/15
    = 𝑒 (30 − 𝑧)2 = 𝑓(𝑧)
    𝑑𝑧 2 5 + 𝑧
    𝑑𝑦
    =𝑣
    𝑑𝑧
    16
    𝑑𝑣 106 𝑧 −𝑧/15
    = 𝑒 (30 − 𝑧)2
    𝑑𝑧 5 + 𝑧
    𝑧0 = 𝟎 𝑦0 = 𝟎 𝑣0 = 𝟎 𝒉 = 𝟎. 𝟓
    𝑧𝑛+1 = 𝑧𝑛+1 +h

    𝑦𝑛+1 = 𝑦𝑛 + ℎ ∗ 𝑣𝑛
    𝑣𝑛+1 = 𝑣𝑛 + ℎ ∗ 𝑓(𝑧𝑛 )

    i Z y v=y’
    0 0 0 0

    1 0.5 0 0

    2 1.0 0 0.00061216

    3 1.5 0.00030608 0.00166118

    b) Runge-Kutta 2
    𝑧𝑛+1 = 𝑧𝑛+1 + ℎ
    𝑘1 = ℎ ∗ 𝑣𝑛
    𝑙1 = ℎ ∗ 𝑓(𝑧𝑛 )
    𝑘2 = ℎ ∗ (𝑣𝑛 + 𝑙1 )
    𝑙1 = ℎ ∗ 𝑓(𝑧𝑛 + ℎ)
    𝑦𝑛+1 = 𝑦𝑛 + 1/2 ∗ (𝑘1 + 𝑘2 )
    𝑣𝑛+1 = 𝑣𝑛 + 1/2 ∗ (𝑙1 + 𝑙2 )
    i z y v=y’

    0 0 0 0
    1 0.5 0 0.00030608

    2 1.0 0.00030608 0.00113667

    3 1.5 0.00113667 0.00233959

    c) Programa
    % RK2
    clear all, f=inline('16/1e6*z*exp(-z/15)*(30-z)^2/(5+z)')
    z(1)=0, y(1)=0, v(1)=0, h=0.5
    for i=1:3
    k1=h*v(i),l1=h*f(z(i)), k2=h*(v(i)+l1), l2=h*f(z(i)+h)
    z(i+1)=z(i)+h, y(i+1)=y(i)+0.5*(k1+k2), v(i+1)=v(i)+0.5*(l1+l2)
    end

    También podría gustarte