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