close all 
% folder = 'G:\BACHELORARBEIT\Auswertungen\3RC-Modell\Tabellen Messung zu Simulation';
% if ~exist(folder, 'dir')
%     mkdir(folder);
% end
% baseFileName = 'SpeicherGes_cvcd2.xlsx';
% fullFileName = fullfile(folder, baseFileName);
% % l=l+500;

e= 13300;

Y1=yout{1}.Values.Data;
Y2=yout{2}.Values.Data;
Y3=yout{3}.Values.Data;
t=(length(Y1));
Z1=theta1(1:t);
Z2=theta2(1:t);
Z3=theta3(1:t);
%% Plot
figure(1)
plot(1:t,Y1,1:t,Y2,1:t,Y3,1:t,Z1,1:t,Z2,1:t,Z3,'LineWidth',1.75);
set(gca,'FontSize',12);
set(gca,'FontName','Arial');
title('...','Color', 'none')
xlabel('Zeit [h]')
ylabel('Temperatur [°C]')
ylim([15 80]);
xlim([0 e]);
xticks(0:3600:e);
xticklabels({'0','1','2','3','4','5','6','7'});
yyaxis right
ylabel('...','Color', 'none');
yticks(0:0:1);
set(gca,'LooseInset',get(gca,'TightInset'))
% grid on
h = figure(1); %neues figure-fenster erzeugen
set(h,'Units','centimeters','Position',[1 1 16 10]); %Rand Rechts, Links Breite Höhe

legend('Simulation Speicher Oben','Simulation Speicher Mitte','Simulation Speicher Unten','Messung Speicher Oben','Messung Speicher Mitte','Messung Speicher Unten','Location','southeast');

%% Abweichung
M1 = rms(Y1-Z1); %Mittlerer Fehler
M2 = rms(Y2-Z2); 
M3 = rms(Y3-Z3); 
figure(2)
plot(1:t,Y1-Z1,1:t,Y2-Z2,1:t,Y3-Z3,'LineWidth',0.25);
set(gca,'FontSize',12);
set(gca,'FontName','Arial');
title('...','Color', 'none')
xlabel('Zeit [h]')
ylabel('Temperatur [K]')
xlim([0 l]);
xticks(0:3600:25200);
xticklabels({'0','1','2','3','4','5','6','7'});
set(gca,'LooseInset',get(gca,'TightInset'))
legend('Differenz Speicher Oben','Differenz Speicher Mitte','Differenz Speicher Unten');

%% Datenausgabe
beschriftung={'Zeit','ThetaObenMessung', 'ThetaMitteMessung','ThetaUntenMessung','ThetaObenSim','ThetaMitteSim','ThetaUntenSim'};
Dat=[(1:t)',Z1,Z2,Z3,Y1,Y2,Y3];
ausgabe=[beschriftung;num2cell(Dat)];
xlswrite(fullFileName ,ausgabe);
xlswrite(fullFileName ,ausgabe);