% ResultsSort
% Sort journey lengths by time of day to create energy used per time of day
% distribution. Uses mid-time of the journey for time stamp.
%
% Author George Hilton. email: g.hilton@soton.ac.uk
% Date 2/2/2017

load('Results.mat')

% Order journeys into time of day.

TimeofDayStart = regexp([Results{:,2}],'T','split');
TimeofDayStart = TimeofDayStart.';
TimeofDayEnd = regexp([Results{:,3}],'T','split');
TimeofDayEnd = TimeofDayEnd.';

TimeofDay = zeros(2287,2);

for i1 = 1:length(TimeofDay)
    
    hour_start = str2double(TimeofDayStart{i1,1}{2}(1:2));
    minute_start = str2double(TimeofDayStart{i1,1}{2}(4:5));
    
    hour_end = str2double(TimeofDayEnd{i1,1}{2}(1:2));
    minute_end = str2double(TimeofDayEnd{i1,1}{2}(4:5));
   
    minutes_start = hour_start*60 + minute_start;
    minutes_end = hour_end*60 + minute_end;
    
    minutes_middle(i1) = (minutes_end + minutes_start)/2;
    
    hour = floor(minutes_middle(i1)/60);
    minute = minutes_middle(i1) - hour*60;
    
    TimeofDay(i1,1) = hour;
    TimeofDay(i1,2) = minute; 
end

% B is time figures, I is order
[B,I] = sort(minutes_middle);

for i2 = 1:length(Results)
    SortedDistanceResults(i2) = Results(I(i2),1);
end

SortedDistanceResultsMat = cell2mat(SortedDistanceResults);

for i3 = 1:length(I)
TimeofDaySorted(i3) = TimeofDay(I(i3));
end

% ts = timeseries(SortedDistanceResultsMat,TimeofDaySorted);
% plot(ts)
% hold 
% 
% ts1 = resample(ts,[0:0.025:1]);
% 
% plot(ts1)

% times = datetime(TimeofDaySorted);

for i4 = 1:24
    
    i0 = i4 - 1;
    
    data = [];
    
    for i5 = 1 : length(TimeofDay)

        if TimeofDaySorted(i5) < i4 && TimeofDaySorted(i5) >= i0
            if SortedDistanceResultsMat(i5)>1
                data = vertcat(data,SortedDistanceResultsMat(i5));
            end
        end
    end
    
    str = sprintf('hourmidlong%d',i4);
    save(str,'data')
    
end