-
Notifications
You must be signed in to change notification settings - Fork 1
/
plotCC.m
29 lines (25 loc) · 731 Bytes
/
plotCC.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
function plotCC(Um, b, c, filename)
% plot marginal utility function
figure
hold on
grid on
box off
tt = 0:30:1440;
pdf = @(u0, b0, c0) u0./pi./c0./( ( (tt-b0)/c0 ).^2 + 1 );
lsty = {'-ro', '-rv', '-r^'};
for i = 1:length(Um)
if i == 1
plot(tt/1440.0, pdf(Um(i), b(i)+1440, c(i)), lsty{i}, 'LineWidth', 0.8)
end
plot(tt/1440.0, pdf(Um(i), b(i), c(i)), lsty{i}, 'LineWidth', 0.8)
end
set(gca, 'XTick', .0:1.0/6.0:1.0);
datetick('x', 'HH:MM', 'keepticks');
pbaspect([2 1 1])
xlabel('Time of the day')
ylabel('Marginal utility')
axis([.0 1.0 0.0 6])
pbaspect([2 1 1])
legend({'Home-AM', 'Home-PM', 'Work', 'Shopping'}, 'Location', 'NorthWest')
legend('boxoff')
export_fig(['FIGURES/' filename], '-pdf', '-jpg', '-r150')