精馏塔理论塔板数,第二讲 图解法确定精馏塔理论板层数(...

第二讲 图解法确定精馏塔理论板层数（matlab code）
（主讲人：北方星）
各位好！今天给大家列举了的是MATLAB软件的一个应用的实例,实例如下：



（matlab软件运行下的代码生成的图形）
其他的本人就不多说了，请大家练习，matlab功能强大，作为计算及绘图工具软件，掌握了后会给大家带来高效的实验数据处理能力，她的好处就不多讲了，本讲只是作为一个抛砖引玉，希望各位在上网看视频，玩游戏累了的时候，可以玩玩此软件，权当休息休息。
%图解法确定精馏塔理论板层数
function fdistillation
%图解法确定理论板层数
clear,clc
%1.已知数据
xD = 0.975;%塔顶组成
xW = 0.0235;%塔釜组成
xF = 0.44;%进料组成
R = 3.5;%回流比
q = 1.362;%进料热状况
vledata = [%汽液平衡数据
1.00 0.780 0.581 0.412 0.258 0.130 0
1.00 0.900 0.777 0.633 0.456 0.262 0];
%2 绘制平衡线
plot (vledata(1,:),vledata(2,:),'*')
hold _disibledevent=>
x1 = [0.0 :0.01:1.0];
y1 = polyval(p,x1);
plot(x1,y1)
hold _disibledevent=>
yc = recope(xc,R,xD);
% 7 绘制提馏段操作线
line([xc,xW],[yc,xW],'Color','c')
% 8 绘制梯级
x0 = xD;
y0 = xD;
n=0;
while 1
n = n+1;
[x0,y0] = tri(x0,y0,R ,xD,p ,xc ,yc,xW,n);
if x0<xW
break;
end
end
% 9修饰
s = sprintf('理论板数为%d',n);
text(0.1,0.9,s)
xlim([0,1]),ylim([0,1]),
title('图解法确定精馏塔理论板数'),
xlabel('x'),
ylabel('y')
%----精馏段操作线----
function y = recope(x,R,xD)
y = R*x/(R+1)+xD/(R+1);
%----进料线方程------
function y = qline(x,q,xF)
y = q*x/(q-1)-xF/(q-1);
%----交点------
function y = cross(x,R,xD,q,xF)
y = recope(x,R,xD)-qline(x,q,xF);
%-----平衡关系逆运算----
function y = vlei(x,p,yn)
y = polyval(p,x)-yn;
%-----绘制梯级------
function [x,y] = tri(x0,y0,R,xD,p,xc,yc,xW,n)
x = fzero(@vlei,x0,optimset('fzero'),p,y0);
if x>xc
y = recope(x,R,xD);
else
y = yc + (yc-xW)*(x-xc)/(xc-xW);
end
line ([x0,x,x],[y0,y0,y],'Color','m'),
text(x,y0,int2str(n),'verticalAlignment','bottom')
grid on

延伸阅读