Matlab 气分 精馏模拟,Matlab 模拟精馏塔

题目如图。

这是代码

%design parameters

xf=0.750;

xp=0.995;

xb=0.040;

p=[300 250 225];

a=[1.1486215 1.1630570 1.1730967];

b=[-0.0532519 -0.0608571 -0.0621885];

c=[-0.0174405 -0.0171428 -0.0217089];

%min reflux ration

for kk=1:3

alp=a(kk)+b(kk)*xf+c(kk)*xf^2;

yalp=alp*xf./(1+(alp-1)*xf);

slmn=(xp-yalp)./(xp-xf);

rmin(kk)=slmn./(1-slmn);

end

%stage count

n=[0 0 0 0 0 0 0 0 0 0 0 0];

for kk=1:3

rrst=rmin(kk)*1.05;

for jj=1:1600;       (我不明白这里的jj到底是什么,我就没有用)

%rectifying stage

rr=rrst*(1+jj/200);

nr=0;

slope=rr./(rr+1);

yr1=xp;

xr1=xp;

while xr1>xf

alr=a(kk)+b(kk)*xr1+c(kk)*xr1^2;

xr2