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);