exercise11.m
function exercise11()
% Lösung von Aufgabe 11 mit fzero
mdot = 5.668e-3; % in kg/s
Q12dot = 3225; % in W
t1 = 70; % in °C
M = 44.01e-3; % in kg/mol, für CO2
% Daten einlesen
data = load("cps.dat");
tD = data(:,1); % in °C
cpD = data(:,3)/M; % in J/(kg K)
% zu lösende Gleichung als Funktion von t2
f = @(t2) Q12dot - mdot*cp(t1, t2, tD, cpD)*(t2 - t1);
% Intervall für die Lösung finden und lösen
epsi = 1e-5;
cpmin = min(cpD);
t2 = fzero(f, [t1 + epsi, t1 + Q12dot/(mdot*cpmin)])
end
%---------------------------------------------------------------------
function y = cp(t1, t2, tD, cpD)
% berechnet die mittlere spezifische Wärmekapazität von t1 bis t2
% t1 Starttemperatur in °C
% t2 Endtemperatur in °C
% tD Tabellenwerte für t
% cpD Tabellenwerte für cp
% y Wärmekapazität in J/(kg K)
cp_t1 = interp1(tD, cpD, t1);
cp_t2 = interp1(tD, cpD, t2);
y = (t2*cp_t2 - t1*cp_t1)/(t2 - t1);
end