exercise12.m

function exercise12()
% Loesung von Aufgabe 12
R = 8.3145;  % in J/(mol K)

% cp-Daten einlesen
load 'cps.dat';
t = cps(:,1);    % in Grad Celsius
cp_O2 = cps(:,2);  % in J/(mol K)
cp_CO2 = cps(:,3); 
cp_NH3 = cps(:,4); 

% Berechnung unter der Annahme delta t = 100 = const.
anz = size(t,1) - 1;
ti = t(1:anz) + 50;
cpi_O2 = (cp_O2(2:anz+1).*t(2:anz+1) - cp_O2(1:anz).*t(1:anz))/100;
cpi_CO2 = (cp_CO2(2:anz+1).*t(2:anz+1) - cp_CO2(1:anz).*t(1:anz))/100;
cpi_NH3 = (cp_NH3(2:anz+1).*t(2:anz+1) - cp_NH3(1:anz).*t(1:anz))/100;

% Freiheitsgrade
f_O2 = 2*(cpi_O2/R - 1);
f_CO2 = 2*(cpi_CO2/R - 1);
f_NH3 = 2*(cpi_NH3/R - 1);

% und plotten
plot(ti, f_O2, ti, f_CO2, ti, f_NH3)
title('Zahl der Freiheitsgrade einiger Gase', 'FontSize', 16);
xlabel('t [\circ C]', 'FontSize', 14);
ylabel('f', 'FontSize', 14);
legend('O_2', 'CO_2', 'NH_3', 2); 
axis([0 1000 0  16]);

F = getframe(gcf);
imwrite(F.cdata, 'bild21.png');