Casestudie Breakdown prediction Contell PILOT - Transumo

%Contains [Duration, Number of Occurences, Percentage, Accumulated
%Percentage]
DurationDL = [];
OccurenceTemp = 0; %For Performance Reason
for i = min(AlarmsDL(:,3)):max(AlarmsDL(:,3));
OccurenceTemp = histc(AlarmsDL(:,3), round(i*10)/10);
if isempty(DurationDL);
DurationDL = [DurationDL; [i, OccurenceTemp,
(OccurenceTemp/length(AlarmsDL(:,1)))*100,
(OccurenceTemp/length(AlarmsDL(:,1)))*100]];
OccurenceTemp = 0;
else
if OccurenceTemp > 0;
DurationDL = [DurationDL; [i, OccurenceTemp,
(OccurenceTemp/length(AlarmsDL(:,1)))*100,
sum(DurationDL(1:length(DurationDL(:,1)),3)) +
(OccurenceTemp/length(AlarmsDL(:,1)))*100]];
end
OccurenceTemp = 0;
end
end
%-----Check Probability of Current Temperature (within Alarming Situations)
%(Calculated by using the maximum values per alarm)
%Contains [Maximum Temperature, Number of Occurences, Percentage,
Accumulated
%Percentage]
ProbabilityDL = [];
OccurenceTemp = 0; %For Performance Reason
for i = min(AlarmsDL(:,4)):0.1:max(AlarmsDL(:,4));
OccurenceTemp = histc(AlarmsDL(:,4), round(i*10)/10);
if isempty(ProbabilityDL);
ProbabilityDL = [ProbabilityDL; [i, OccurenceTemp,
(OccurenceTemp/length(AlarmsDL(:,1)))*100,
(OccurenceTemp/length(AlarmsDL(:,1)))*100]];
OccurenceTemp = 0;
else
if OccurenceTemp > 0;
ProbabilityDL = [ProbabilityDL; [i, OccurenceTemp,
(OccurenceTemp/length(AlarmsDL(:,1)))*100,
sum(ProbabilityDL(1:length(ProbabilityDL(:,1)),3)) +
(OccurenceTemp/length(AlarmsDL(:,1)))*100]];
end
OccurenceTemp = 0;
end
end
%-----Durations of Dooropenings-----
%Contains: [Date, Duration (in seconds)]
Dooropeningtime = [];
%Get Startingposition for non interpolated data
startposition = 1;
while floor(date(startposition) + second) < floor(interpolation(start,1) +
second);
startposition = startposition + 1;
end
% Get Duration of Dooropenings
A-131