Proyecto Listado Txy Xy

Post on 17-Dec-2015

212 views 0 download

description

Matlab

Transcript of Proyecto Listado Txy Xy

% Calcular las temperaturas de ebullicin de los componentes puros y determinar cul es el compuesto ms voltilclear all%Ecuacin de Antoine para Benceno y Hexano Ab= 6.89745; Bb=1206.350; Cb=220.237;Ah= 6.87773; Bh=1171.530; Ch=224.366;P=input('Introduzca la presin a una atmsfera o a la presin del laboratorio de la en mm de Hg')Tbb=(Bb./(Ab-(log10(P))))-Cb;Tbh=(Bh./(Ah-(log10(P))))-Ch;%Generar un vector de temperaturas desde la temperatura de ebullicin del ms voltil hasta la del menos voltil.%Para cada temperatura del vector generado calcular la presin de vapor de cada componente, la fraccin mol en la fase lquida (x) y la fraccin mol en la fase vapor (y), stas ltimas del compuesto ms voltil, segn sea el casocondicion1=TbbTbhif condicion1T=Tbb:0.1:TbhPb= 10.^(Ab-(Bb./(T+Cb)));Ph= 10.^(Ah-(Bh./(T+Cb)));x=(P-Ph)./(Pb-Ph)y=(x.*Pb)./Pelseif condicion2T=Tbh:0.1:TbbPb=10.^(Ab-(Bb./(T+Cb)));Ph=10.^(Ah-(Bh./(T+Ch)));x=(P-Pb)./(Ph-Pb)y=(x.*Ph)./Pend%Graficar: Temperatura vs. x, y (diagrama de fases) plot(x,T,y,T)% Graficar y vs. x (curva de equilibrio lquido-vapor)plot(x,y)