cola_dv.m

国外经典书籍MULTIVARIABLE FEEDBACK CONTROL-多变量反馈控制 的源码

function xprime=cola_dv(t,X) 
% sample usage:   [t,x]=ode15s('cola_dv',[0 5000],0.5*ones(1,82));
%
% cola_dv - Subroutine for simulation with DV-configuration.
%
%
%            Inputs are distillate (D) and boilup (VB). Disturbances
%            are feedrate and feed composition. These are set by directly
%            altering 'cola_dv.m'. Outputs are liquid composition and
%            liquid hold up for stages 1 through NT, given in x. 

% Number of stages in the column
NT=41;

% Inputs and disturbances
D=0.5;                               % Distillate
VB=3.20629;                          % Boilup
F=1.0 + 0.00;                        % Feedrate
zF=0.5;                              % Feed composition
qF=1.0;                              % Feed liquid fraction

% P-Controllers for control of reboiler and condenser hold up.
KcB=10;  KcD=10;         % controller gains
MDs=0.5; MBs=0.5;        % Nominal holdups - these are rather small  
Ls=2.70629; Bs=0.5;      % Nominal flows
MB=X(NT+1);  MD=X(2*NT); % Actual reboiler and condenser holdup
LT=Ls+(MD-MDs)*KcD;       % Reflux
B=Bs+(MB-MBs)*KcB;       % Bottoms flow     

% Store all inputs and disturbances
U(1)=LT; U(2)=VB; U(3)=D; U(4)=B; U(5)=F; U(6)=zF; U(7)=qF;

xprime=colamod(t,X,U);