exercise11.m

% Loesung von Aufgabe 11

% Daten einlesen
load -ascii 'kreisprozess.dat';
V = kreisprozess(:,1);
pu = kreisprozess(:,2);
po = kreisprozess(:,3);

% Kreisprozess zeichnen
plot(V, po, 'k', V, pu, 'k');
title('Kreisprozess eines Verbrennungsmotors');
xlabel('V [cm^3]');
ylabel('p [bar]');

% Umrechnung auf SI-Einheiten
V = 1e-6*V;     % in m^3
pu = 1e5*pu;    % in Pa
po = 1e5*po;    % in Pa

% Bestimmung der Kreisprozess-Arbeit
intu = trapz(V, pu);
into = trapz(V, po);
W = into - intu

% Bestimmung der Leistung
n = 3000/60;   % in 1/s
P = W * n      % in W