Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 43 vistas

    Primer Examen - Programacion Chavez Quispe Alan Ruben

    Cargado por

    Ruth Ticona
    Este documento contiene tres tareas relacionadas con análisis numérico. La primera involucra calcular números de Fibonacci usando el triángulo de Pascal y sumando las diagonales. La segunda resuelve la ecuación de estado de Redlich-Kwong para metano. La tercera desarrolla un programa en MATLAB para leer datos termodinámicos de una base de datos, calcular el calor específico en función de la temperatura y graficarlo, así como calcular el calor necesario para calentar un compuesto.

    Copyright:

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

    Primer Examen - Programacion Chavez Quispe Alan Ruben

    Cargado por

    Ruth Ticona
    43 vistas9 páginas
    Este documento contiene tres tareas relacionadas con análisis numérico. La primera involucra calcular números de Fibonacci usando el triángulo de Pascal y sumando las diagonales. La segunda resuelve la ecuación de estado de Redlich-Kwong para metano. La tercera desarrolla un programa en MATLAB para leer datos termodinámicos de una base de datos, calcular el calor específico en función de la temperatura y graficarlo, así como calcular el calor necesario para calentar un compuesto.

    Descripción original:

    Examen

    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 contiene tres tareas relacionadas con análisis numérico. La primera involucra calcular números de Fibonacci usando el triángulo de Pascal y sumando las diagonales. La segunda resuelve la ecuación de estado de Redlich-Kwong para metano. La tercera desarrolla un programa en MATLAB para leer datos termodinámicos de una base de datos, calcular el calor específico en función de la temperatura y graficarlo, así como calcular el calor necesario para calentar un compuesto.

    Copyright:

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

    Primer Examen - Programacion Chavez Quispe Alan Ruben

    Cargado por

    Ruth Ticona
    Este documento contiene tres tareas relacionadas con análisis numérico. La primera involucra calcular números de Fibonacci usando el triángulo de Pascal y sumando las diagonales. La segunda resuelve la ecuación de estado de Redlich-Kwong para metano. La tercera desarrolla un programa en MATLAB para leer datos termodinámicos de una base de datos, calcular el calor específico en función de la temperatura y graficarlo, así como calcular el calor necesario para calentar un compuesto.

    Copyright:

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

    UNIVERSIDAD MAYOR DE SAN ANDRES

    FACULTAD DE INGENIERIA

    Tarea Recuperatoria Primer Examen de Programación y Análisis


    Numérico
    DOCENTE:
    Ing. M.Sc. Roberto Parra Zeballos.
    ESTUDIANTE:
    Universitario Alan Ruben Chávez Quispe.
    CARRERA:
    Ingeniería Química.
    MATERIA:
    Programación y Análisis Numérico.
    RU:
    1789105
    1.- A partir del triángulo de Pascal, encuentre n números de la serie de Fibonacci.

    Algoritmo:
    1. Dado el valor del exponente “n” en el binomio de newton construir el triángulo de
    pascal
    2. Llenar la primera columna y diagonal de “1” en la matriz “t”
    3. Para i = 3 hasta n + 1
    a. Para j = 2 hasta i
    i. t(i,j) = t(i-1,j-1) + t(i-1,j)
    4. Sumar las diagonales para formar la serie de Fibonacci guardando las sumas en “f(i)”
    5. Mostrar el vector “f” serie de Fibonacci

    Diagrama de flujo:
    Programa en Matlab.
    Sub – programa T(n+1) x (n+1)
    function T(n+1)x(n+1)
    N = input('Ingresar el valor de N = ');
    for i = 1:n+1:i+1
    t(i,1)=1
    t(i,i)=1
    end
    for i = 3:n+1:i+1
    for j = 2:i:i+1
    t(i,j)= t(i-1,j-i)+t(i-1,j)
    end
    end
    end

    Programa principal
    clc, clear
    function T(n+1)x(n+1);%El primer cálculo de la función
    for i = 1:n+1:i+1
    Sum=0
    for j = 1:i:j+1
    Sum= Sum+ T(j,i-j+1)
    end
    F(i)=Sum
    end
    fprintf('Nuevo vector de la serie de Fibonacci (f):
    %5.10f',f)
    end
    2.- Para gases en condiciones reales, un modelo adecuado es la ecuación de estado de
    Redlich-Kwong que está dada por:

    Algoritmo:
    1. Crear el Sub - programa FUNCTION tomando en cuenta que para el metano: Pc=4580 kPa y

    Tc= 191K (R=0.518).

    2. A partir de la tolerancia “tol”, aproximación inicial de v “v” y primera evaluación de la función en


    “f= F(v)”

    3. Mientras que cumpla que I f I >=tol

    a) Intentar con otro valor de “v” incrementándolo con un paso, sea: v= v+ 10^ (-12)
    b) Evaluar la función en “f= F(v)”
    𝑉 3
    4. Calcular la masa: 𝑚 = 𝑣
    →𝑚=𝑣

    5. Mostrar resultados “m” y “v”.


    Diagrama de flujo:
    Sub – programa FUNCTION

    Programa principal
    Programa en Matlab.
    Sub – programa FUNCTION
    function [f] = F(v)
    %ENTRADAS: v
    %SALIDA: f
    R = 0.518 ;
    Tc = 191 ;
    Pc = 4580 ;
    a = 0.427*((R^2) *(Tc^2.5))/Pc ;
    b = 0.0866*(R*Tc)/Pc ;
    f = 65000 - (0.518*223)/ (v - b) + (a)/ (v*(v + b) *sqrt (223)) ;
    end

    Programa principal
    clc, clear
    tol = 10^(-4) ;
    v = 0.002740; %1ra aproximación
    f = F(v) ; %El primer cálculo de la función
    while abs(f) >= tol
    v = v + 10^(-12) ;
    f = F(v) ;
    end
    fprintf('volumen especifico es: %5.10f
    [m^37kg]\n',v) m = 3/v ;
    fprintf('masa de metano es: %5.5f [kg] \n',m)
    3. El calor específico a presión constante (Cp) para una sustancia química, se calcula según la
    siguiente ecuación, donde A, B, C son constantes características de cada compuesto:

    Usando la base de datos de PROPIEDADES.xlsx (Ud. debe construir esta base de datos en
    EXCEL, por lo menos de 50 compuestos), realizar un programa en MATLAB capaz de ejecutar
    la siguientes tareas: i) Leer el nombre de un compuesto por teclado, ii) Para este compuesto
    calcular el Cp en función de la temperatura (el rango deberá ser definido por teclado), iii)
    Realizar una gráfica, que muestre el comportamiento del Cp con la temperatura en el rango
    especificado, iv) Calcule el calor necesario para calentar un compuesto líquido desde la
    temperatura de 5 °C hasta su temperatura de ebullición, debe leer por teclado el nombre del
    compuesto y debe calcular la temperatura de ebullición a la presión local de 495 mm Hg.
    PROPIEDADES.xlsx es una base de datos que tiene en la primera hoja (sheet1) de nombre “Pv”,
    datos para calcular la Pv a diferentes temperaturas con la ecuación de Antoine, la segunda hoja
    (sheet2) es Cp que son los datos para calcular el calor específico a presión constante. El número
    de compuestos debe ser por lo menos 700 registros que tienen, por ejemplo, 5 campos: 1)
    Nro, 2) Nombre Compuesto, 3) A, 4) B, 5) C.

    Algoritmo:
    1. Leer el nombre del compuesto.

    2. Exportar la Base de datos PROPIEDADES..xlsx asignando “n” al número de registros y

    “m” al número de campos.

    3. Realizar la búsqueda del compuesto desde i=1 hasta i=n

    4. Una vez identificado guardar el número del registro en la variable “N”.

    5. Asignar a las variables “A”,“B” y “C” los valores para el compuesto encontrado.

    6. Leer el rango de temperatura a graficar “ti” y “ts” y una temperatura “t”.

    7. Calculamos los vectores a graficar en el eje horizontal “temperatura” y en el eje vertical

    “Capacidad calorífica”, a partir de la ecuación de Cp = A + BT + CT^2

    8. Mostramos el gráfico y resultados.


    Diagrama de flujo:

    Base de satos de EXCEL


    Programa en Matlab.
    clc,clear
    comp = input('Ingresar el nombre del compuesto = ') ; P =
    readtable('PROPIEDADES.xlsx','Range','A1:E11')
    ;%Exportar datos
    [n,m] = size(P) ; %Dimensiones de la tabla

    %Encuentra la posicion del compuesto


    for i = 1:n
    if 1 == strcmpi(comp,P.Nombre(i))
    N = i ;
    end
    end

    A = P.A(N) ;
    B = P.B(N) ;
    C = P.C(N) ;
    %ENTRADA DE DATOS
    ti = input('Ingresar la temperatura inferior ºC = ') ;
    ts = input('Ingresar la temperatura superior ºC = ') ;
    t = input('Ingresar la temperatura ºC = ') ;
    Cp = (A + B*t + C*t^2)*1000 ; %[J/mol K]

    %GRAFICANDO
    x = ti:1:ts ; % temperatura
    y = (A + B*x + C*x.^2)*1000 ; % Cp
    hold on
    plot(x,y)
    plot(t,Cp,'ro')
    xlabel('Temperatura [ºC]')
    ylabel('El calor específico a presión constante [J/mol K]')
    title('Cp en funcion de la temperatura')
    hold off

    %IMPRESIONES
    disp(strcat('El compuesto es: ',{comp}))
    fprintf('A la temperatura %5.2f[ºC] el Cp = %5.2f[J/mol
    K] \n',t,Cp)

    También podría gustarte