hudongkang.m

直接带入数据..可以求解理论总塔板数以及提馏和精馏段塔板数

function hdk
clear all;clc
xi=[0.00 0.010 0.020 0.040 0.060 0.080 0.100 0.140 0.180 0.200 0.250 0.300 0.350 0.400 0.450 0.500 0.550 0.600 0.650 0.700 0.750 0.800 0.850 0.894 0.900 0.950 1];
yi=[0.00 0.110 0.170 0.270 0.340 0.392 0.430 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.820 0.855 0.894 0.898 0.942 1];
xf=0.1436;xd=0.8545;xw=0.0007836;
R=4;R1=4;
sp=spline(xi,yi);
x4plot=linspace(xi(1),xi(end),200);
y4plot=fnval(sp,x4plot);
plot(xi,yi,'*',x4plot,y4plot,'-')
N=quadl(@func1,xf,xd,[],[],sp,xd,R);
M=quadl(@func2,xw,xf,[],[],sp,xw,R1);
N=round(N+0.5);
M=round(M+0.5);
disp=('计算结果:');
fprintf('\n 精馏段理论塔板数N为:%d%s\n',N,'块')
fprintf('\n 提馏段理论塔板数M为:%d%s\n',M,'块')
fprintf('\n 共需理论塔板数N为:%d%s\n',N+M,'块')
function f=func1(x,sp,xd,R)
y=fnval(sp,x);
f=1./(y-x-(xd-y)/R);
function f=func2(x,sp,xw,R1)
y=fnval(sp,x);
f=1./(y-x-(y-xw)/R1);