Modelagem de uma coluna de destilação no Matlab

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

%====QUESTÃO 2 - COLUNA DE DESTILAÇÃO====
%CÁLCULO DO BOLHA P UTILIZANDO O MODELO DE WILSON
%ENTRAR COM T e Xi PARA CALCULAR P e Yi
clear all, clc
%==================DADOS=================
A = xlsread('valores_de_x');
X1 = A(:,1);
X2 = 1-X1;
T = 50;
TK = T + 273.15; %K
R = 1.9872; %cal/mol.K
 
%Coeficientes de Wilson
L12=0.209163999;
L21=0.82284181;
 
%=========CONSTANTES DE ANTOINE=========
%ETANOL (1) E ÁGUA (2)
A1 = 18.5242; B1 = 3578.91; C1 = -50.5;
A2 = 18.3036; B2 = 3816.44; C2 = -46.13;
 
%====CÁLCULO DA PRESSÃO DE SATURAÇÃO====
P1sat = exp(A1-(B1/(TK+C1)));
P2sat = exp(A2-(B2/(TK+C2)));
fprintf('P1sat(etanol): %f\n', P1sat)
fprintf('P2sat(agua): %f\n', P2sat)
 
%=CÁLCULO DOS COEFICIENTES DE ATIVIDADE=
gama1 = exp(-log(X1+(X2.*L12))+(X2.*((L12./(X1+(X2.*L12)))-(L21./(X2+X1.*L21)))));
gama2 = exp(-log(X2+(X1.*L21))-(X1.*((L12./(X1+(X2.*L12)))-(L21./(X2+X1.*L21)))));
%===========CÁCULO DA PRESSÃO===========
P = (gama1.*X1.*P1sat)+(gama2.*X2.*P2sat);
fprintf ('P:%f\n', P)
%===CÁLCULO DAS COMPOSIÇÕES DE VAPOR====
Y1 = (gama1.*X1.*P1sat)./P;
%Y2 = 1-y1;
fprintf ('Y1:%f\n', Y1)