MATLAB图解精馏塔理论塔板数
程序代码(总6页)
-CAL-FENGHAI.-(YICAI)-Company One1 -CAL-本页仅作为文档封面，使用请直接删除
MATLAB图解精馏塔理论塔板数程序代码
function distillation %文件名“distillation”可以更改 % 输入计算所需参数 q=1;%输入进料热状况参数 R=;%输入回流比
xD=;%输入塔顶轻组分摩尔分数 xW=;%输入塔底轻组分摩尔分数 xF=;%输入进料轻组分摩尔分数
%以下输入相平衡数据 x0=[0
1]; y0=[0
1];
Yr=@(x)R/(R+1).*x+xD/(R+1);%精馏段操作线
fun=@(x)(q-1)*(R/(R+1).*x+xD/(R+1))-(q*(x-xF)+(q-1)*xF); xQ=fzero(fun,;%求操作点 yQ=Yr(xQ); xOP=[xW,xQ,xD]; yOP=[xW,yQ,xD];
yfit=linspace(0,1,1001);
xfit=interp1(y0,x0,yfit,'pchip'); %%绘制图形 hold on box on
plot([0 1],[0 1],'k'); xlabel('x') ylabel('y') plot(x0,y0,'r')
plot(xfit,yfit,'r-') plot(xF,xF,'b*') plot(xQ,yQ,'bo') plot(xOP,yOP,'b-') k=1; yn(1)=xD;
xn(1)=interp1(y0,x0,yn(1),'pchip'); plot([xD,xn(1)],[yn(1),yn(1)],'b-') text(xn(1),yn(1),num2str(1),...
'HorizontalAlignment','center','VerticalAlignment','bottom') while xn(k)>xW
yn(k+1)=interp1(xOP,yOP,xn(k)); k=k+1;
xn(k)=interp1(y0,x0,yn(k),'pchip');
plot([xn(k-1),xn(k-1)],[yn(k-1),yn(k)],'b-') plot([xn(k-1),xn(k)],[yn(k),yn(k)],'b-') text(xn(k),yn(k),num2str(k),...
'HorizontalAlignment','center','VerticalAlignment','bottom' ) end N=k;
plot([xn(N),xn(N)],[yn(N),xn(N)],'b-') text(xn(N),yn(N),num2str(N),...
'HorizontalAlignment','center','VerticalAlignment','bottom' )
N_Feed=find(xntext,,{strcat('所需理论板：',num2str(N)),... strcat('进料板位置：',num2str(N_Feed))},...'HorizontalAlignment','left','VerticalAlignment','top')
%以下代码是为了去掉顶端和右边坐标轴的刻度 box off
ax2 = axes('Position',get(gca,'Position'),... 'XAxisLocation','top',... 'YAxisLocation','right',... 'Color','none',...
'XColor','k','YColor','k'); set(ax2,'YTick', []); set(ax2,'XTick', []); box on
运行结果如图1所示：
10.90.80.780.60.50.4110.30.2130.114000.10.20.30.40.5x0.60.70.80.9112910所需理论板：14进料板位置：675231y
图1 图解苯-甲苯精馏塔理论塔板数和进料位置
