matlab塔板数

%图解法确定精馏塔理论塔板数
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