dtlabvvx3der.m

换热器模型的MATLAB实现

%% Returnerar dT f鰎 varm och kall str鰉 f鰎 
%% 
%% Motstr鰉
%% F鰎 simulering ensam

%% T(1) 鋜 kallvatten ut
%% T(2n) 鋜 varmvatten ut
%% T(n) 鋜 kallvatten in
%% T(n+1) 鋜 varmvatten in
%%

%% 1 2 .. n
%% <-------- kv in
%% -------->
%%n+1 .. 2*n 

%In kommer 1:n-1, n+1:2n, dvs 2n-2 st

function der = dTlabvvx3der(t,T,flag,par)
n = length(T)/2; %9

rho = par(1);
cp = par(2);
Fk = par(3);
Fv = par(4);
k = par(5);
Tkin = par(6);
Tvin = par(7);
Atot = par(8);
Vtot = par(9);

V = Vtot/n;
A = Atot/n;

T = [Tkin; T(1:n); Tvin; T(n+1:2*n)];

%Kall
for i = 2:n+1,
   dT(i) = 1/(rho*cp*V)*(rho*cp*Fk*(T(i-1)-T(i))+k*A*(T(i+n+1)-T(i)));
end

%Varm
for i = (n+3):2*n+2,
   dT(i) = 1/(rho*cp*V)*(rho*cp*Fv*(T(i-1)-T(i))-k*A*(T(i)-T(i-n-1)));
end

der = [dT(2:n+1) dT(n+3:2*n+2)]';