-
Notifications
You must be signed in to change notification settings - Fork 1
/
plotActUtil.m
61 lines (52 loc) · 1.6 KB
/
plotActUtil.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
function plotActUtil(act_util, times)
tt=times(2,1,:);
% plot activity utility
% 1: home, 2:work
figure
box on
% subplot(2,1,1)
grid on
hold on
x = genTimes(13, 12+size(act_util,1));
colors = lines(2);
u1 = flipud(cumsum(flipud(act_util(:,1))));
u2 = cumsum(act_util(:,2));
[AX,H1,H2] = plotyy(x,u2,x,u1,'plot');
set(H1,'LineStyle','--', ...
'Marker', '^', ...
'Color', colors(1,:), ...
'MarkerFaceColor', colors(1,:));
set(H2,'LineStyle',':', ...
'Marker', 'v', ...
'Color', colors(2,:), ...
'MarkerFaceColor', colors(2,:));
set(get(AX(1),'Ylabel'),'String','Cummulative utility (Work)')
set(get(AX(2),'Ylabel'),'String','Cummulative utility (In-Home)')
pbaspect(AX(1), [2 1 1])
pbaspect(AX(2), [2 1 1])
datetick(AX(1), 'x','HHPM')
xt = genTimes(12, 12+size(act_util,1)+1);
set(AX(1), 'XTick', xt(1:3:13))
set(AX(2), 'XTick', [])
set(AX(1), 'YTick', 0:20:150)
set(AX(2), 'YTick', 0:20:150)
ylim(AX(1), [0 150])
ylim(AX(2), [0 150])
legend('Work', 'Home', 'Location', 'NorthWest')
xlabel('Time of the day')
% ylim([0 max(cumsum(act_util(:)))*1.5])
% legend('Home', 'Work', 'Location', 'NorthWest')
% % plot travel time
% subplot(2,1,2)
% hold on
% line(x, -tt(:).*alpha, 'LineStyle', ':', ...
% 'Color', 'r', ...
% 'Marker','o', ...
% 'MarkerFaceColor', 'r');
% % datetick('x','HHPM');
% datetick('x','HHPM');
% xlabel('Time of the day')
% ylim([min(-tt(:).*alpha)*1.5 0])
% ylabel('Travel disutility')
% legend_h = legend('Travel from W to H', 'Location', 'NorthEast');
% set(legend_h, 'Color', 'none')