Casestudie Breakdown prediction Contell PILOT - Transumo

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