m4c.m

来自「氧乐果控制过程」· M 代码 · 共 29 行

%*********************************************
%              Model for Control
%     function out=m4c(in1,in2,mat1,struct,W1,B1,W2,B2)
% 参数说明：out: 反应釜温度
%           in1: 投料阀开度
%           in2: 反应时刻
%           mat1: 液位矩阵
%           struct: 流量网络结构
%*********************************************

function out=m4c(in1,in2,mat1,struct,W1,B1,W2,B2)

if nargin<8 
    error('Not Enough Input Arguments!');
end

% 计算一甲胺流量
level=mat1;
pr1=minmax(level);
levelstd=(level-pr1(1,1))/(pr1(1,2)-pr1(1,1));      % 归一化处理
p1=[levelstd(in2+1);in1];
flux=sim(struct,p1);                                % 归一化流量值

% 计算反应釜温度
fluxsum=levelstd(in2+1)-levelstd(1);
time=(in2)/59;
p2=[flux;fluxsum;time];
t=purelin(W2*tansig(W1*p2,B1),B2);                  % 归一化温度值
out=-12/t;