cstr_optdes.m

来自「这是一些关于MATLAB的小程序」· M 代码 · 共 33 行

function CSTR_OptDes
% 连续搅拌槽式反应器的最优设计
%
%   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 = [];       % 使用默认参数
xL=[0  0];          % 设定下限
xU=[1  1];          % 设定上限

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

fprintf('       优化结果：\n\n')
fprintf('   转化率为：%.3f\n',xA )
fprintf('   浓度CA0：%.3f 摩尔/升\n',CA0)
fprintf('   进料量FA0为：%.3f 升/小时\n',100/(CA0*xA) )
fprintf('   反应器体积V为：%.2f 升\n',50/CA0^2/(1-xA)^2 )
fprintf('   每小时总消耗费用：%.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;