function CSTR_OptDes
% 連續(xù)攪拌槽式反應(yīng)器的最優(yōu)設(shè)計(jì)
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2002/05/06 $

clear all
clc

x0 = [0.5  0.5];    % 給定初值
options = [];       % 使用默認(rèn)參數(shù)
xL=[0  0];          % 設(shè)定下限
xU=[1  1];          % 設(shè)定上限

[x,f,flag] = fmincon(@Fee,x0,[],[],[],[],xL,xU)
CA0 = x(1);
xA = x(2);

fprintf('       優(yōu)化結(jié)果：\n\n')
fprintf('   轉(zhuǎn)化率為：%.3f\n',xA )
fprintf('   濃度CA0：%.3f 摩爾/升\n',CA0)
fprintf('   進(jìn)料量FA0為：%.3f 升/小時(shí)\n',100/(CA0*xA) )
fprintf('   反應(yīng)器體積V為：%.2f 升\n',50/CA0^2/(1-xA)^2 )
fprintf('   每小時(shí)總消耗費(fèi)用：%.3f 元\n',f )

% ------------------------------------------------------------------
function f = Fee(x)             % Objective function for Optimization
CA0 = x(1);
xA = x(2);
f = 400*CA0^0.4/xA + 4.183/CA0^1.2/(1-xA)^1.2;