# How to use datetime to plot value vs. hour of day, i.e. daily cycle?

24 views (last 30 days)
K E on 21 Jul 2016
Commented: Peter Perkins on 3 Aug 2016
Let's say I have a time series of humidity with datetime time stamps, which are randomly sampled in time over several days though sample time is monotonic. I would like to plot humidity vs. the monotonic hour of day that the measurement was made, so I can see the daily cycle of humidity (high values around 3pm, regardless of the day the data was collected, for example). My first plot below with datetime fails at doing this. My second plot with datenum shows how I want it to look (monotonic hours on x-axis). I want to know if I can make the same plot using datetime - this will tell me whether to bother switching over to datetime, which has some nice features, or keep using datenum in my code.
%%Failed attempt to make plot using datetime
figure('Name', 'Failed Plot Using Datetime');
% Time vector is monotonic but randomly sampled in time
tDays = datetime(2014,6,28) + days(0:1/24:10) + ...
rand(size(days(0:1/24:10)))*1/24;
% Humidity measurement corresponding with time stamp
humidity = rand(1,length(tDays));
% Attempt to convert from days to hours of day
tHours = datetime(tDays, 'Format', 'HH');
plot(tHours, humidity, 'x');
xlabel('Not Monotonic Hours between 0-24');
ylabel('Humidity');
title('Failed Plot Using Datetime');
%%How I want the plot to look
% However this uses datenum and I want to use datetime instead
figure('Name', 'Plot Using Datenum');
tDays = datenum(2014,6,28) + (0:1/24:10) + rand(size(0:1/24:10));
humidity = rand(size(tDays));
% Conversion from days to hours
tHours = (tDays - floor(tDays))*24;
plot(tHours, humidity, 'x');
xlabel('Monotonic Hours between 0-24');
ylabel('Humidity');
title('Desired Plot, But Not Using Datetime');

Sean de Wolski on 25 Jul 2016
Edited: Sean de Wolski on 26 Jul 2016
EDIT
[tHours.Year, tHours.Month, tHours.Day] = deal(0);
plot(tHours,humidity,'x','DatetimeTickFormat','HH')
Peter Perkins on 3 Aug 2016
This is not a good way to do this. You want to plot vs. durations, created using the timeofday function:
plot(timeofday(tDays),humidity,'o')