batchdist.m

来自「黄华江所著《实用化工计算机模拟—MATLAB在化学工程中的应用》的所有例子源程序」· M 代码 · 共 20 行

function BatchDist
% 简单蒸馏计算
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2003/02/26 $

clear all
clc

alpha = 2.5;    % 相对挥发度
x0 = 0.60;      % 开始时的苯含量    
var = 2;        % var = F0/F

[x,fval] = fzero(@DistEq,x0,[],alpha,x0,var)

% ------------------------------------------------------------------
function f = DistEq(x,alpha,x0,var)
f = log(var) - 1/(alpha-1)*(log(x0/x)+alpha*log((1-x)/(1-x0)));