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