Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 38 vistas

    Prueba Escrita - Unidad 5-SOL

    Cargado por

    Wilson Choque
    Este documento presenta tres problemas relacionados con ecuaciones diferenciales ordinarias. El primer problema involucra resolver analíticamente y numéricamente una ecuación diferencial de valor inicial, y calcular el error de truncamiento. El segundo problema implica resolver analíticamente y numéricamente otra ecuación diferencial de valor inicial utilizando el método de Heun, y calcular el error absoluto exacto. El tercer problema pide resolver una ecuación diferencial mediante el método de Runge-Kutta de cuarto orden en MatLab, y elaborar un programa para graficar las soluc

    Copyright:

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

    Prueba Escrita - Unidad 5-SOL

    Cargado por

    Wilson Choque
    38 vistas6 páginas
    Este documento presenta tres problemas relacionados con ecuaciones diferenciales ordinarias. El primer problema involucra resolver analíticamente y numéricamente una ecuación diferencial de valor inicial, y calcular el error de truncamiento. El segundo problema implica resolver analíticamente y numéricamente otra ecuación diferencial de valor inicial utilizando el método de Heun, y calcular el error absoluto exacto. El tercer problema pide resolver una ecuación diferencial mediante el método de Runge-Kutta de cuarto orden en MatLab, y elaborar un programa para graficar las soluc

    Título original

    Prueba escrita - Unidad 5-SOL

    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 tres problemas relacionados con ecuaciones diferenciales ordinarias. El primer problema involucra resolver analíticamente y numéricamente una ecuación diferencial de valor inicial, y calcular el error de truncamiento. El segundo problema implica resolver analíticamente y numéricamente otra ecuación diferencial de valor inicial utilizando el método de Heun, y calcular el error absoluto exacto. El tercer problema pide resolver una ecuación diferencial mediante el método de Runge-Kutta de cuarto orden en MatLab, y elaborar un programa para graficar las soluc

    Copyright:

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

    Prueba Escrita - Unidad 5-SOL

    Cargado por

    Wilson Choque
    Este documento presenta tres problemas relacionados con ecuaciones diferenciales ordinarias. El primer problema involucra resolver analíticamente y numéricamente una ecuación diferencial de valor inicial, y calcular el error de truncamiento. El segundo problema implica resolver analíticamente y numéricamente otra ecuación diferencial de valor inicial utilizando el método de Heun, y calcular el error absoluto exacto. El tercer problema pide resolver una ecuación diferencial mediante el método de Runge-Kutta de cuarto orden en MatLab, y elaborar un programa para graficar las soluc

    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 DEL VALLE

    Prueba escrita

    Estudiante: Carrera:
    Asignatura: Métodos numéricos Grupo: B

    SedeCentral
    Gestión: II/2021 Fecha: 16/11/2021
    SubsSedes:LPZ SUC TDD
    Unidades temáticas a evaluar: UNIDAD 5: Ecuaciones diferenciales ordinarias
    Nombre y Firma Docente(s): Lucy Copa Gerónimo
    1. Resuelva el siguiente problema de valor inicial en el intervalo de x  1 a 4, donde
    y1  1
    dy
     4  6x 2  8x 3
    dx
    a) (5 puntos) Analíticamente (Cálculo exacto)
    b) (10 puntos) Método de Euler (Con 6 interaciones)
    c) (15 puntos) Calcular el error de truncamiento local.

    Solución:
    a)
    dy
     4  6x 2  8x 3
    dx
     dy  4  6x 2  8x 3 dx
    y  4x  2x 3  2x 4  C
    Reemplazando la condición inicial y1  1
    y  41  21 3  21 4  C  1  C  1
    y  4x  2x 3  2x 4  1
    b)
    h  0. 5
    n xi y i  y i1  hfx i1 , y i1  fx, y  4  6x 2  8x 3
    0 1 1 Condición inicial  10
    1 1. 5 1  0. 510  4 36. 5
    2 2 4  0. 536. 5  22. 25 84
    3 2. 5 22. 25 0. 584  64. 25 158. 5
    4 3  64. 25  0. 5158. 5  143. 5 266
    5 3. 5 143. 5  0. 5266  276. 5 412. 5
    6 4 276. 5  0. 5412. 5  482. 75
    f  x i , y i  2 f  x i , y i  3 f 3 x i , y i  4
    Et  h  h  h
    2! 3! 4!

    c)
    f  x i , y i   12x  24x 2
    f  x i , y i   12  48x
    f 3 x i , y i   48

    El error de truncamiento del término de segundo orden:


    f  x i , y i  2 121  241 2
    E t,2  h  0. 5 2  4. 5
    2! 2!
    Para el término de tercer orden:
    f  x i , y i  3 12  481
    E t,3  h  0. 5 3  1. 25
    3! 3!
    Para el término de cuarto orden:
    f 3 x i , y i  4
    E t,4  h  48 0. 5 4  0. 125
    4! 4!
    E t  E t,2  E t,3  E t,4
    E t  4. 5  1. 25  0. 125  5. 875
    2. Dada la ecuación diferencial
    dy
    x 2  1  xy  0, y5  1, x  5, 8
    dx
    a) (5 puntos) Resolver analíticamente.
    b) (20 puntos) Resolver utilizando el método de Heun h  1
    c) (10 puntos) Calcular el error absoluto exacto.

    Solución:
    a)
    dy ln y  ln C
    x 2  1  xy  0
    dx x2  1
    dy
    x 2  1  xy
    dx
    y  C
    2
    x 1
    dy x
    y  x 2  1 dx
    Si y5  1

    1  C  C  24
    ln y  1 lnx 2  1  ln C 2
    2 5 1

    dy xy
    b)  2
    dx x 1
    x0  5 ; y0  1
    Iteración 1: i  0

    x1  x0  h  5  1  6
    Predictor: y 01  y 0  hfx 0 , y 0   1  10. 208 33  0. 791 67
    fx 0 , y 0   fx 1 , y 01 
    Corrector: y 1  y 0  h
    2
     1  1 0. 208 33  0. 135 71  0. 827 98
    2

    Iteración 2: i  1
    x2  x1  h  6  1  7
    Predictor: y 02  y 1  hfx 1 , y 1   0. 827 98  10. 141 94  0. 686 04
    fx 1 , y 1   fx 2 , y 02 
    Corrector: y 2  y 1  h
    2
     0. 827 98  1 0. 141 94  0. 100 05  0. 706 99
    2

    Iteración 3: i  2
    x3  x2  h  7  1  8
    Predictor: y 03  y 2  hfx 2 , y 2   0. 706 99  10. 103 1  0. 603 89

    fx 2 , y 2   fx 3 , y 03 
    Corrector: y 3  y 2  h
    2
     0. 706 99  1 0. 103 1  0. 076684  0. 617 10
    2

    c)
    Iteración y aprox y exacto error absoluto exacto
    1 0. 827 98 0. 828 08 0. 000 1
    2 0. 706 99 0. 707 11 0. 000 12
    3 0. 617 10 0. 617 21 0. 000 11

    3. Dada la ecuación diferencial a valor inicial,


    dy 2x
     y0  1, x  0, 0. 5
    dx 3y 2 x 2  1
    a) (20 puntos) Resolver la ecuación diferencial mediante el método de Runge Kutta de
    cuarto orden con tamaño de paso   0.5 mediante la programación en MatLab
    b) (15 puntos) Elaborar un programa que realice el cálculo de la solución exacta, la
    solución aproximada por el método de Euler y también por el método de Runge Kutta de
    cuarto orden, que visualice la gráfica de las tres soluciones en una misma figura
    (combinar los programas realizados en clase).

    Solución:
    x 0  0 , y 0  1, x f  0. 5 , h  0. 5
    k 1  fx 0 , y 0   0
    k 2  f x 0  h , y 0  h k 1  0. 1569
    2 2
    k 3  f x 0  , y 0  h k 2  0. 1452
    h
    2 2
    k 4  fx 1 , y 0  hk 3   0. 2318

    k 1  2k 2  k 3   k 4
    y1  y0  h
    6

    0  20. 1569  0. 1452  0. 2318


    y 1  1  0. 5  1. 069 7
    6

    b)
    PROGRAMA:
    %runge Kutta de cuarto orden
    %Datos: una función f(x,y) continua con fy también continua
    % sobre un rectángulo R, un punto (x0,y0 ) en el interior
    % de R, y el número de iteraciones adicionales n, y el
    % tamaño de paso h
    %syms x y y(x)
    f  input(’f(x,y)  ’);
    x0  input(’x0  ’);
    y0  input(’y0  ’);
    h  input(’h  ’);
    n  input(’iter  ’);
    %Calculo con Euler
    X1  zeros(n1,1);
    Y1  zeros(n1,1);
    X1(1)  x0;
    Y1(1)  y0;
    for i  1 : n
    X1(i1)  X1(i)  h;
    Y1(i1)  Y1(i)  h * subs(f, {x,y}, {X1(i),Y1(i)});
    end
    %Calculo con Runge Kutta de cuarto orden
    X  zeros(n1,1);
    Y  zeros(n1,1);
    X(1)  x0;
    Y(1)  y0;
    hm  h / 2;
    hs  h / 6;
    for i  1 : n
    X(i1)  X(i)  h;
    k1  double(subs(f, {x, y} , {X(i), Y(i)})); %k1f(x0,y0)
    k2  double(subs(f, {x, y}, {X(i)  hm, Y(i)  hm * k1})); %k2f(x0h/2,y0h/2*k1)
    k3  double(subs(f, {x, y}, {X(i)  hm, Y(i)  hm * k2})); %k3f(x0h/2,y0h/2*k2)
    k4  double(subs(f, {x, y}, {X(i1), Y(i)  h * k3})); %k4f(x0h,y0h*k3)
    Y(i1)  double(Y(i)  hs * (k1  2 * (k2  k3)  k4)); %y1y0h/6*(k12(k2k3)k4)
    end
    %Solucion verdadera
    yt(x)  dsolve(diff(y)  f, y(x0)  y0);
    Yt  double(yt(X));
    %Graficas
    close all
    plot(X,Y,’r’)
    hold on
    plot(X1,Y1,’g’)
    XX  X(1):0.01:X(n1);
    plot(XX,double(yt(XX)),’b’)
    legend(’Runge Kutta 4’,’Euler’ ,’Verdadera’)
    grid on
    %Mostrar resultados
    fprintf(’Solucion verdadera y(x)  %s \n’, char(simplify(yt)))
    disp(’X(i)  ’)
    fprintf(’%.7f\n’, X)
    disp(’Solucion verdadera, Y(X)  ’)
    fprintf(’%.7f\n’, Yt)
    disp(’Solucion aproximada, Yrk4(X) ’)
    fprintf(’%.7f\n’, Y)
    disp(’Solucion aproximada, Yeuler(X) ’)
    fprintf(’%.7f\n’, Y1)

    SALIDA:
    f(x,y)  2*x/((3*y^2*x^2)(3*y^2))
    x0  0
    y0  1
    h  0.5
    iter  1
    Solucion verdadera y(x)  (log(x^2  1)  1)^(1/3)
    X(i) 
    0.0000000
    0.5000000
    Solucion verdadera, Y(X) 
    1.0000000
    1.0694467
    Solucion aproximada, Yrk4(X) 
    1.0000000
    1.0696666
    Solucion aproximada, Yeuler(X) 
    1.0000000
    1.0000000

    También podría gustarte