-
Notifications
You must be signed in to change notification settings - Fork 1
/
figure_4_plot_params.m
68 lines (45 loc) · 1.29 KB
/
figure_4_plot_params.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
62
63
64
65
66
67
68
%% initialize
figdir = [pwd '\Figures\'];
if ~exist(figdir)
mkdir(figdir)
end
%% parameters - extero
tickparam = 0.5:0.1:1;
rtickparam = fliplr(tickparam);
zlab = 'Entropy (H)';
xlab = '\alpha';
ylab = '\beta';
param_mat_in = param_mat_extero;
%% figure
figure
cmap = cbrewer2('div', 'Spectral', 50);
s=surf(tickparam, tickparam, param_mat_in); zlabel(zlab), xlabel(xlab),ylabel(ylab); colormap(cmap)
set(gca, 'XTickLabel', tickparam, 'XTick',tickparam, 'YTickLabel', tickparam, 'YTick', tickparam)
az = 140; el = 20;
zlim([0 35]);
caxis([5 35]);
view(az,el);
pbaspect([1 1 1])
fname = [figdir 'fig4A_extero.pdf'];
title('Exteroceptive Uncertainty')
%print(fname, '-dpdf', '-r600');
%% parameters - intero
tickparam = 0.5:0.1:1;
rtickparam = fliplr(tickparam);
zlab = 'Entropy (H)';
xlab = '\alpha';
ylab = '\beta';
param_mat_in = param_mat_intero;
%% figure
figure
cmap = cbrewer2('div', 'Spectral', 50);
s=surf(tickparam, tickparam, param_mat_in); zlabel(zlab), xlabel(xlab),ylabel(ylab); colormap(cmap)
set(gca, 'XTickLabel', tickparam, 'XTick',tickparam, 'YTickLabel', tickparam, 'YTick', tickparam)
az = 140; el = 20;
zlim([0 35]);
caxis([5 35]);
view(az,el);
pbaspect([1 1 1])
title('Interoceptive Uncertainty')
fname = [figdir 'fig4C_intero.pdf'];
%print(fname, '-dpdf', '-r600');