%圖解法確定精餾塔理論塔板數
function?tabanshu
clearclc
xd=0.747;
xw=0.016;
xf=0.063;%
R=3.91;
q=1.166;%
%繪制平衡線
vledata=[
????0.0?0.004?0.01?0.02?0.04?0.06?0.08?0.10?0.14?0.18?0.20?0.25?0.30?0.35?0.40?0.45?0.50?0.55?0.60?0.65?0.70?0.75?0.80?0.85?0.894?0.90?0.95?1.0
????0.0?0.053?0.11?0.175?0.273?0.34?0.392?0.43?0.482?0.513?0.525?0.551?0.575?0.595?0.614?0.635?0.657?0.678?0.698?0.725?0.755?0.785?0.82?0.855?0.894?0.898?0.942?1.0];
%plot（vledata（1:）vledata（2:）‘*‘）
hold?on
p=polyfit（vledata（1:）vledata（2:）9）;
x1=[0:0.001:1];
y1=polyval（px1）;
plot（x1y1‘LineWidth‘2）
hold?on;
plot（xwxw‘*‘）
hold?on;
plot（xdxd‘*‘）
hold?on;
line（[xdxd][0xd]‘LineWidth‘2‘Color‘‘k‘）
hold?on;
line（[xwxw][0xw]‘LineWidth‘2‘Color‘‘k‘）
hold?on;
%繪制對角線
line（[0.01.0][0.01.0]）
hol