[Jozgoo]

主程序ScriptEx5.m 调用的EquationEx5.m (ode15s函数中调用) 如下

function dstate = EquationEx5(t,state,pars)

CA = state(1); %
CB = state(2); %
T = state(3); %
CC = state(4); %
CD = state(5); %
CE = state(6); %

%%%%%%%%%%%%%%%%%%
Q = pars (1);
V = pars(2);
CA_0 = pars(3);
CB_0 = pars(4);

MyA = pars(5);
MyB = pars(6);
CA_s = pars(7);
CB_s = pars(8);
T_s = pars(9);
UA = pars(10);
Rho = pars(11);
Cp = pars(12);
T_a = pars(13);

DHR = pars(14);
A = pars(15);
B = pars(16);
AlphaA = pars(17);
AlphaB = pars(18);
T_0 = pars(19);
CC_0 = pars(20);
CD_0 = pars(21);
CE_0 = pars(22);
MyC = pars(23);
MyD = pars(24);
MyE = pars(25);
CC_s = pars(26);
CD_s = pars(27);
CE_s = pars(28);

%%%%%%%%%%%%%%%%%%


if((t >= 400) && (t <= 414))
Q = 0;
elseif ((t >= 900) && (t <= 910))
Q = 0;
else
Q = 5.3*10^(-6);
end

% Balances

dCAdt = Q/V*(CA_0 - CA) + MyA*(exp(A - B/T)*CA^AlphaA*CB^AlphaB);

dCBdt = Q/V*(CB_0 - CB) + MyB*(exp(A - B/T)*CA^AlphaA*CB^AlphaB);

dTdt = Q/V*(T_0 - T) + UA/(V*Rho*Cp)*(T_a - T) + (-DHR)/(Rho*Cp)*(exp(A - B/T)*CA^AlphaA*CB^AlphaB);

dCCdt = Q/V*(CC_0 - CC) + MyC*(exp(A - B/T)*CA^AlphaA*CB^AlphaB);

dCDdt = Q/V*(CD_0 - CD) + MyD*(exp(A - B/T)*CA^AlphaA*CB^AlphaB);

dCEdt = Q/V*(CE_0 - CE) + MyE*(exp(A - B/T)*CA^AlphaA*CB^AlphaB);


dstate = [dCAdt dCBdt dTdt dCCdt dCDdt dCEdt]';