%% Mittelwert Tagestemperatur
% Importieren der Daten

data = readtable(['temperatur/data_1.csv'])
T = table2array(data(:,2));
t = hours(table2array(data(:,1)));

%
% Plot

figure
plot(t,T)
grid on
xlabel('Time [h]')
ylabel('Temperature [DegC]')

%
% Berechnung Mittelwert
% \overline{f} = \frac{1}{b-a}\int_a^b f(x)dx

fbar = trapz(t,T) / (t(end) - t(1))

% Plot mit Mittelwert

m1 = [t(1),t(end)]
m = [fbar,fbar]
plot(t,T,m1,m)
grid on
xlabel('Time [h]')
ylabel('Temperature [DegC]')
mean(T) % 05:10, 05:20,... fehlen!