大师作文网matlab编程：图解法计算吸收塔理论塔板数
来源：学生作业帮 编辑：大师作文网作业帮 分类：综合作业 时间：2024/11/19 20:52:00
matlab编程：图解法计算吸收塔理论塔板数
图解法求吸收塔的理论塔板数的过程如下图所示.图中虚线称为操作线,由吸收塔的进出口组成确定,如图中的a 和b 点,两点坐标分别为（0,0.0101）和（0.00524,0.25）.图中实线表示汽液平衡线y=20.7x.求解过程从a 点出发,坐水平线与平衡线相交,交点处作垂直线与操作线相交,以此类推,直至垂直线与操作性的交点超过b 点.水平线的数目即为理论塔板数.试编写一个MATLAB程序重复以上求解过程,要求结果输出塔板数和图像.
clear
xa=0; %a点横坐标
ya=0.0101; %a点纵坐标
xb=0.00524; %b点横坐标
yb=0.25; %b点纵坐标
x=xa; %起点横坐标初值
y=ya; %起点纵坐标初值
k=0; %塔板数初值
while y<=yb
x1=y/20.7; %水平线与平衡线交点横坐标
y1=(yb-ya)/(xb-xa)*(x1-x)+y; %垂直线与操作线交点纵坐标
plot([x,x1],[y,y]) %绘制水平线段
hold on
plot([x1,x1],[y,y1]) %绘制垂直线段
k=k+1;
x=x1;
y=y1;
end
plot([xa,x],[ya,y],'--') %绘制操作线
plot([xa,xb],[ya,yb],'*') %绘制a,b两点
text(xa,ya+0.015,'a'); %标注a点
text(xb,yb+0.015,'b'); %标注b点
plot([0,x],[0,20.7*x],'-') %绘制平衡线
k %输出塔板数
k =
4
即理论塔板数为4
xa=0; %a点横坐标
ya=0.0101; %a点纵坐标
xb=0.00524; %b点横坐标
yb=0.25; %b点纵坐标
x=xa; %起点横坐标初值
y=ya; %起点纵坐标初值
k=0; %塔板数初值
while y<=yb
x1=y/20.7; %水平线与平衡线交点横坐标
y1=(yb-ya)/(xb-xa)*(x1-x)+y; %垂直线与操作线交点纵坐标
plot([x,x1],[y,y]) %绘制水平线段
hold on
plot([x1,x1],[y,y1]) %绘制垂直线段
k=k+1;
x=x1;
y=y1;
end
plot([xa,x],[ya,y],'--') %绘制操作线
plot([xa,xb],[ya,yb],'*') %绘制a,b两点
text(xa,ya+0.015,'a'); %标注a点
text(xb,yb+0.015,'b'); %标注b点
plot([0,x],[0,20.7*x],'-') %绘制平衡线
k %输出塔板数
k =
4
即理论塔板数为4