INSTITUTO POLITÉCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS
METODOS NUMERICO.
PROYECTO.
NEWTON RAPHSON MULTIVARIABLE.
“EJERCICIO 2.3”.
INTEGRANTES:
PORTILLO OLIVERA AILEEN GRISEL.
SUSANO ARGUETA MARIA DE JESÚS.
XOLIO MEJÍA ANDREA.
PROFESORA:
GRECIA VÁZQUEZ CAMARILLO.
GRUPO: 1IM44.
PERIODO: 19/1.
FECHA DE ENTREGA: 03/10/2018.
CONTENIDO
ENUNCIADO DEL PROBLEMA ................................................................Error! Bookmark not defined.
METODO POLINOMIO DE LAGRANGE. ..................................................Error! Bookmark not defined.
ALGORITMO. ....................................................................................................................................... 4
PROGRAMA. ........................................................................................................................................ 5
EJECUCION DEL PROGRAMA. .............................................................................................................. 7
RESULTADOS Y COMPROBACION ........................................................................................................ 8
CONCLUSION. ...................................................................................................................................... 8
BIBLIOGRAFIA. ..................................................................................................................................... 8
0 2
Métodos Numéricos Aplicados a la Ingeniería
Dada la siguiente tabla calcule:
La presión de saturación (bar) @ una temperatura de 100°C.
El volumen especifico (cm3/g) @ una presión de 0.4 bar y una temperatura de 180°C.
La temperatura @ una presión de 0.8 bar y una entalpia de 3000 KJ/Kg.
Aplique polinomio de Lagrange usando para ello todos los puntos.
Complete la tabla con los valores generados.
Tabla. Propiedades de vapor sobrecalentado.
v(cm3/g), u(KJ/Kg), h(KJ(Kg) y s(KJ/Kg K).
Temp v u h s v u h s
°C.
Psat=0.06 bar (36.16 °C) Psat=0.35 bar (72.69 °C)
80 2713.2 2650.1 4625 2645.6
120 3021.9 2726.0 5163 2723.1
160 3330.2 2802.5 5696 2800.6
200 3638.3 2879.7 6228 2878.4
240 3946.2 2957.8 6758 2956.8
280 4254.0 3036.8 7787 3036.0
Psat=0.70 bar (89.95 °C) Psat=1.00 bar (99.63 °C)
80 2434 2680.0 1696 2676.2
120 2571 2719.6 1793 2716.6
160 2841 2798.2 1984 2796.2
200 3108 2876.7 2172 2875.3
240 3374 2955.5 2359 2954.5
280 3640 3035.0 2546 3034.2
0 3
ALGORITMO.
Se generan las ecuaciones que
Se escriben las 3 ecuaciones en
corresponden a la solucion del
forma de funcion con respecto a
ejercicio. Al ser 3 incognitas, son
sus 3 variables.
3 ecuaciones
Se deriban parcialmente las 3 Se asigna un valor supuesto para
funciones con respecto a sus 3 cada variable, y se sustituye en
variables, por lo tanto se las 3 funciones y en las 9
obtienen 16 derivadas parciales derivadas parciales
Se calculan los ∆ de cada
variable, por el metodo de
Se obtiene el valor del Jacobiano, crammer, siendo mi factor
con una matriz de 3x3 constante de division el
jacobiano calculado
anteriormente
Se obtienen los nuevos valores
de mis variables para la siguiente Con el nuevo valor de la variable,
iteracion con la siguiente se realiza el numero de
formula. iteraciones necesarias htasa
Vnueva=Vanterior+∆V obtener el valor buscado.
donde V es cualquier variable.
0 4
Programa.
%******************************************************
%** Interpolacion con el Metodo **
%** del Polinomio de Lagrange
clear;
clc;
%*********************************************************
%** Se ingresan los puntos de referencia para realizar la*
%** interpolacion, asi como el numero para el que se desea
%** interpolar **
%*********************************************************
fprintf('Interpolacion con el Metodo del Polinomio de Lagrange\n\n');
n=input('grado del polinolio: ');
for i=1:n+1
x(1,i)=input('dame los valores de xi:');
end
for i=1:n+1
xi(1,i)=input('dame los valores de f(xi):');
end
xi
xint=input('Numero para el que desea interpolar x: ');
0 5
%**********************************************************
%** En la siguiente seccion se realizan los calculos para *
%** determinar el valor aproximado resultado de la interpo*
%** lacion. **
%**********************************************************
fxint=0;
i=1;
while i<=n+1
L=1;
J=0;
while J<=n
if i~=J+1
L=L*(xint-x(1,J+1))/(x(1,i)-x(1,J+1));
end
J=J+1;
end
fxint=fxint+L*xi(1,i);
i=i+1;
end
fprintf('\nresultado xi: %d',fxint');
%******************************************************
%** Enseguida se grafican los puntos que se dieron al**
%** Principio del programa **
%******************************************************
0 6
plot(x,xi)
grid
xlabel('x');ylabel('y')
Ejecución del Programa.
0 7
RESULTADOS Y COMPROBACION
Conclusión.
Bibliografia.
(2018). Como programas en ti nspire cx cas - Foro de calculadoras. [online] Calculadoras.cl.
Available at: http://www.calculadoras.cl/foro/forum/calculadoras/calculadoras-texas-familia-ti-nspire-
touchpad-cx-cx-cas/9709-como-programas-en-ti-nspire-cx-cas [Accessed 23 Mar. 2018].
0 8
