dtlabvvx.m

来自「换热器模型的MATLAB实现」· M 代码 · 共 43 行

%% Returnerar dT f鰎 varm och kall str鰉 f鰎 
%% 
%% Motstr鰉
%% T(1) 鋜 kallvatten ut
%% T(2n) 鋜 varmvatten ut

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

function der = dTlabvvx(T,par)

n = length(T)/2;

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;

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

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

der = dT';