solving a system of eqns
2 次查看(过去 30 天)
显示 更早的评论
I need help seeing the solutions to this problem. It is giving me a matrix answer or a huge fraction. Any help would be appreciated
%12.16
% For one of the binary systems listed in the preceding table, based on
% eq 10.5 and the Wilson eqn, make the following calculations for
% (1)Propanol/(2)water
%a. Bubl P: t=60 deg C, x1=.3
%b. Dew P: t=60 deg C, y1=.3
%c. P,T-flash t=60 deg C, P=.5(Pbubble+Pdew), z1=.3
%d. If an azeotrope exists at t=60 deg C, find Paz and x1_az and y1+az
clear
clc
%Propanol
P_sat1=exp(16.1154-(3483.67/(60+205.807)));
V1=74.05 ;
%Water
P_sat2=exp(16.3872-(3358.79/(60+212.041)));
V2=18.07;
%Given
A12=291.27;
A21=1448.01;
R=8.314;
T=60+273;
lambda12=V2/V1*exp(-A12/(R*T))
lambda21=V1/V2*exp(-A21/(R*T))
%Part A
x1=.3
x2=1-x1
gamma1=exp((x2*(lambda12/(x1+x2*lambda12)-(lambda21/(x2+x1*lambda21)))))/(x1+x2*lambda12)
gamma2=exp((-x1*(lambda12/(x1+x2*lambda12)-(lambda21/(x2+x1*lambda21)))))/(x2+x1*lambda21)
syms y1 y2 P_b
y=solve(y1*P_b==x1*gamma1*P_sat1, y2*P_b==x2*gamma2*P_sat2, y1+y2==1)
**I also tried
[y1 y2 P_b]=solve(y1*P_b==x1*gamma1*P_sat1, y2*P_b==x2*gamma2*P_sat2, y1+y2==1, y1,y2,P_b)
0 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Number Theory 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!