A partir de los datos de la tabla, se requiere determinar la densidad para 30
valores en el rango de temperatura 105≤T≤315, utilizar una interpolación lineal y
grafique los datos obtenidos.
CODIGO EN MATLAD
%PRograma que calculela densidad para 30 valores en un rango determinado de
%temperatura, con nterpolacion y graficas
%Variables
%d=densidad
%T=temperatura
%dtinterpo=interpolacion de temperatura y de densidad
clear
clc
d = [929, 902, 860];
T = [105, 205, 315];
Tinterpo = linspace(105, 315, 30);
dinterpo = interp1(T, d, Tinterpo, 'linear');
figure;
plot(T, d, 'm-o');
hold on;
plot(Tinterpo, dinterpo, 'o--');
xlabel('Temperatura (°C)');
ylabel('Densidad (kg/dm^3)');
title('Interpolación Lineal de Datos');
legend('show');
grid on;
disp('Valores interpolados:');
disp([Tinterpo', dinterpo']);
• Para determinar la composición de una mezcla a partir de un análisis
espectrofotométrico se obtuvieron las siguientes ecuaciones:
11𝐶1 + 118𝐶2 + 9𝐶3 = 26
27𝐶1 + 27𝐶2 + 85𝐶3 = 39
98𝐶1 + 9𝐶2 + 2𝐶3 = 35
Verifique que el análisis fue realizado correctamente obteniendo el determinante
del sistema y en caso de que tenga solución obtenga las concentraciones.
CODIGO EN MATLAD
%Programa para verificar el sistema si hay solucion, y obtener las
%concentraciones
%VARIABLES
%A=matriz de las ecuacion
%B=matriz de las ecuaciones
%C1,2,3=Concentraciones de las ecuaciones
clear
clc
A = [11, 118, 9; 27, 27, 85; 98, 9, 2];
B = [26; 39; 35];
detA = det(A);
if detA ~= 0
C1 = det([B A(:,2:3)]) / detA;
C2 = det([A(:,1) B A(:,3)]) / detA;
C3 = det([A(:,1:2) B]) / detA;
disp('Concentraciones:');
disp(['C1 = ' num2str(C1)]);
disp(['C2 = ' num2str(C2)]);
disp(['C3 = ' num2str(C3)]);
else
disp('El sistema no tiene solución.');
end
• Los datos de presión-temperatura-volumen específico para el etano se muestran
en la tabla siguiente, donde la Temperatura (T) está en °C, la presión (P) en
atmosferas y el volumen específico en moles/litro. a) Calcular el volumen
específico en moles/litro para una presión de 5 atmosferas y una temperatura de
225 °C. b) Calcular el volumen especifico en moles/litro para una temperatura de
200 °C y una presión de 7.5 atmosferas.
CODIGO EN MATLAD
% Programa que calcula el volumen especifico en moles/L dependeindo de la
% atmosfera y temperaturas
%Variables
%T=Temperaturas
%P=presiones
%V=volumen especifico
clear
clc
T = [75, 150, 200, 250];
P = [1, 2, 4, 6, 8];
V = [
24.95, 31.89, 36.44, 40.87;
43.80, 59.31, 69.38, 79.16;
68.89, 106.06, 130.18, 153.59;
85.95, 151.38, 194.53, 237.38;
104.38, 207.66, 276.76, 345.38];
Tinter = 225;
Pinter = 5;
Vatm= bilinear_interpolation(T, P, V, Tinter, Pinter);
disp(['El volumen específico es: ' num2str(Vatm) ' moles/litro']);
Tinter2 = 200;
Pinter2= 7.5;
Vatm2= bilinear_interpolation(T, P, V, Tinter2, Pinter2);
disp(['la atmosferas es: ' num2str(Vatm2) ' ;moles/litro']);
function result = bilinear_interpolation(X, Y, Z, x, y)
[X, Y] = meshgrid(X, Y);
result = interp2(X, Y, Z, x, y);
end
• Dada la siguiente expresión: 𝑒^𝑥/4𝜋 tan(0.75𝑥) – 2 a) Obtener la derivada. b)
Calcular el área bajo la curva en el intervalo de 0.5 a 1.85
CODIGO EN MATLAD
%Programa para obtener la derivada y calcular el area bajo la curva dentro
%d eun intervalo de 0.5 a 1.85
%VARIABLES
%x=variable declarada syms
%expr=la expresion dada
%dexpr=derivada de la expresion
%limi,lims=limite inferior y superior
%ac=area bajo la curva
clear
clc
syms x;
expr= exp(x/(4*pi)) * tan(0.75*x) - 2;
dexpr = diff(expr, x);
disp(['La derivada de la expresión es: ' char(dexpr)]);
expr1= @(x) exp(x/(4*pi)).* tan(0.75*x) - 2;
limi= 0.5;
lims= 1.85;
ac= integral(expr1, limi, lims);
disp(['El área bajo la curva es: ' num2str(ac)]);
• Resuelva los siguientes sistemas de ecuaciones:
0.5𝑠𝑒𝑛(𝑥1𝑥2 ) − 𝑥2/4𝜋 − 0.5𝑥1 = 0
0.920423(𝑒^𝑥1 − 𝑒) + 8.65256𝑥2 − 2𝑒^𝑥1 = 0
50𝑎 − 20𝑏 + 10𝑐 = 25
−10𝑎 + 35𝑏 − 20𝑐 = 30
10𝑎 + 20𝑏 + 48𝑐 = 26
CODIGO EN MATLAD
% Programa para resolver las siguientes ecuaciones dadas
%Variables
%EC1=primer ecuación
%Ec2=segunda ecuación
%Ec3=tercera ecuación
%Ec4=cuarta ecuación
%Ec5=quinta ecuación
clear
clc
syms x1 x2 e a b c;
ec1 = 0.5 * sin(x1 * x2) - x2/(4*pi) - 0.5*x1 == 0;
ec2 = 0.920423 * (exp(x1) - exp(1)) + 8.65256*x2 - 2*exp(x1) == 0;
ec3 = 50*a - 20*b + 10*c == 25;
ec4 = -10*a + 35*b - 20*c == 30;
ec5 = 10*a + 20*b + 48*c == 26;
res= solve([ec1, ec2, ec3, ec4, ec5], [x1, x2, a, b, c]);
fprintf('La solución del sistema es: ');
disp(['x1 = ' char(res.x1)]);
disp(['x2 = ' char(res.x2)]);
disp(['a = ' char(res.a)]);
disp(['b = ' char(res.b)]);
disp(['c = ' char(res.c)])
