Casestudie Breakdown prediction Contell PILOT - Transumo

%Name and location of the source file
unit = 1;
filename = strcat('Channel', int2str(unit), '.csv');
path = 'C:\Dokumente und Einstellungen\Christian\Eigene
Dateien\Dokumente\Studium\Diplomarbeit\Monitoring Data Nijmegen
(Converted)\';
%Import Dataset
import = importdata(strcat(path, filename));
%Create Datevector (as serial date number):
date = datenum(import.textdata(:,1), 'dd-mm-yy HH:MM:SS');
%Import Interpolated Data from Disk
interpolation = importdata(strcat(path, filename, '- Interpolated.txt'));
%Definition of a Second
second = 1/(24*60*60);
%Definition of a Minute (For Performance Reasons)
minute = 1/(60*24);
%Definition of Day- and Nighttime
daybegin = (1/24)*6;
dayend = (1/24)*22;
%Start (Index of the imported data)
%1 = Begin of imported file, add 1440 per Day
start = 1 + 7*1440;
%-----Minima & Maxima-----
%(Use non interpolated data)
%Auxiliary variable
startposition = 1;
while floor(date(startposition) + second) < floor(interpolation(start,1) +
second);
startposition = startposition + 1;
end
jumplength = 0;
%Create Datevector for first column (per day)
dailydate = floor(date(startposition) + second);
%Create Minvector for second column (per day)
minvector = [];
%Create Mindaytimevector for third column (per day)
mindaytime = [];
%Create Minnighttimevector for forth column (per day)
minnighttime = [];
%The same for the maxima table
maxvector = [];
maxdaytime = [];
maxnighttime = [];
for i = startposition:length(date);
if isequal(floor(date(startposition)+ second), floor(date(i) + second));
jumplength = jumplength + 1;
else
A-116