medianvector = [medianvector; median(interpolation(startposition:startposition + jumplength - 1,2))]; modevector = [modevector; mode(interpolation(startposition:startposition + jumplength - 1,2))]; stdvector = [stdvector; std(interpolation(startposition:startposition + jumplength - 1,2))]; %Day- and Nighttime... nighttemp = []; daytemp = []; for j= startposition:startposition + jumplength - 1; if (mod(interpolation(j,1) + second, 1) >= daybegin) && (mod(interpolation(j,1) + second, 1)
xlswrite(strcat(path, 'Excel\', filename, '- Median'), [interpolateddailydate-693960, round(medianvector*10)/10, round(mediandaytime*10)/10, round(mediannighttime*10)/10], 'Median', 'A2'); xlswrite(strcat(path, 'Excel\', filename, '- Mode'), [interpolateddailydate-693960, round(modevector*10)/10, round(modedaytime*10)/10, round(modenighttime*10)/10], 'Mode', 'A2'); xlswrite(strcat(path, 'Excel\', filename, '- Standard Deviation'), [interpolateddailydate-693960, round(stdvector*10)/10, round(stddaytime*10)/10, round(stdnighttime*10)/10], 'Standard Deviation', 'A2'); xlswrite(strcat(path, 'Excel\', filename, '- Doordopenings'), [interpolateddailydate-693960, dailydooropenings, daytimedooropenings, nighttimedooropenings], 'Dooropenings', 'A2'); %Total Values totalmean = mean(interpolation(:,2)); totalmedian = median(interpolation(:,2)); totalmode = mode(interpolation(:,2)); totalstd = std(interpolation(:,2)); totaldooropenings = sum(interpolation(:,3)); %-----Temperature Distribution----- %Count total occurrences of single values %"Round" Command necessary in MATLAB. Otherwise some comparisons fail! %Contains[Temperature, Minutes of Occurence] totalOC = []; for i = min(interpolation(:,2)):0.1:max(interpolation(:,2)); totalOC = [totalOC; [i, sum(interpolation(:,2) == round(i*10)/10)]]; end %-----Plot Statements----- %Temperature Overview plot(date, import.data(:,2), 'k'); hold on; plot(interpolation(:,1), interpolation(:,4), '--b') plot(interpolation(:,1), interpolation(:,5), '--r') datetick('x',20, 'keeplimits'); title 'Temperature Overview'; xlabel 'Date'; ylabel 'Temperature (°C)'; axis([min(date) max(date) min(import.data(:,2)) max(import.data(:,2))]); hold off; print('-dtiff', strcat(path, 'Graphs\', filename, '- Temperature Overview.tif')); %Maximum Values per Day bar([dailydate(1):dailydate(length(maxvector))], maxvector, 'k'); hold on; plot(interpolation(:,1), interpolation(:,4), '--b') plot(interpolation(:,1), interpolation(:,5), '--r') datetick('x',20, 'keeplimits'); title 'Maximum Values per Day'; xlabel 'Date'; ylabel 'Temperature (°C)'; axis([dailydate(1) dailydate(length(maxvector)) min(maxvector) max(maxvector)]); A-121