-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathMC_Plots.m
161 lines (138 loc) · 9.19 KB
/
MC_Plots.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
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
%% Generates box plots in Figures 5-6 for HB. %%
T1_F = 0.45; T1_S = 1.4; T2_F = 0.015; T2_S = 0.09; M0_F = 0.15; k_FS = 8; Delta = 0; PC1 = 0 + Delta; PC2 = pi + Delta;
cm = [lines(4);lines(1);lines(1)];
Solution_Bouhrara = squeeze(Solution_Bouhrara);
Solution_Deoni = squeeze(Solution_Deoni);
Solution_Wood = squeeze(Solution_Wood);
Solution_Zhang = squeeze(Solution_Zhang);
Solution_ASDeoni = squeeze(Solution_ASDeoni);
Solution_ASRui = squeeze(Solution_ASRui);
% T1S
figure(1); subplot(3,2,1);
x = [Solution_Bouhrara(:,1);Solution_Deoni(:,1);Solution_Wood(:,1);Solution_Zhang(:,1);Solution_ASDeoni(:,1);Solution_ASRui(:,1)];
g = [ones(1000,1); 2*ones(1000,1); 3*ones(1000,1); 4*ones(1000,1); 5*ones(1000,1); 6*ones(1000,1)];
BP1 = boxplot(x,g,'Colors',cm,'Symbol','','Labels',{'S1 B1','S1 B2','S1 B3','S1 B4','S2 B1','S3 B1'});
for ih=1:7
set(BP1(ih,:),'LineWidth',2);
end
ylabel('T_{1S} (s)','FontSize',20); hline(T1_S,'k--'); ylim([1 2.5]);
get(gca, 'XTick'); set(gca, 'FontSize', 18); get(gca, 'YTick'); set(gca, 'FontSize', 18);
grid on; grid minor;
% T1F
subplot(3,2,2);
x = [Solution_Bouhrara(:,2);Solution_Deoni(:,2);Solution_Wood(:,2);Solution_Zhang(:,2);Solution_ASDeoni(:,2);Solution_ASRui(:,2)];
g = [ones(1000,1); 2*ones(1000,1); 3*ones(1000,1); 4*ones(1000,1); 5*ones(1000,1); 6*ones(1000,1)];
BP1 = boxplot(x,g,'Colors',cm,'Symbol','','Labels',{'S1 B1','S1 B2','S1 B3','S1 B4','S2 B1','S3 B1'});
for ih=1:7
set(BP1(ih,:),'LineWidth',2);
end
ylabel('T_{1F} (s)','FontSize',20); hline(T1_F,'k--');
get(gca, 'XTick'); set(gca, 'FontSize', 18); get(gca, 'YTick'); set(gca, 'FontSize', 18);
grid on; grid minor;
% T2S
subplot(3,2,3);
x = [Solution_Bouhrara(:,3);Solution_Deoni(:,3);Solution_Wood(:,3);Solution_Zhang(:,3);Solution_ASDeoni(:,3);Solution_ASRui(:,3)];
g = [ones(1000,1); 2*ones(1000,1); 3*ones(1000,1); 4*ones(1000,1); 5*ones(1000,1); 6*ones(1000,1)];
BP1 = boxplot(x,g,'Colors',cm,'Symbol','','Labels',{'S1 B1','S1 B2','S1 B3','S1 B4','S2 B1','S3 B1'});
for ih=1:7
set(BP1(ih,:),'LineWidth',2);
end
ylabel('T_{2S} (s)','FontSize',20); hline(T2_S,'k--');
get(gca, 'XTick'); set(gca, 'FontSize', 18); get(gca, 'YTick'); set(gca, 'FontSize', 18);
grid on; grid minor;
% T2F
subplot(3,2,4);
x = [Solution_Bouhrara(:,4);Solution_Deoni(:,4);Solution_Wood(:,4);Solution_Zhang(:,4);Solution_ASDeoni(:,4);Solution_ASRui(:,4)];
g = [ones(1000,1); 2*ones(1000,1); 3*ones(1000,1); 4*ones(1000,1); 5*ones(1000,1); 6*ones(1000,1)];
BP1 = boxplot(x,g,'Colors',cm,'Symbol','','Labels',{'S1 B1','S1 B2','S1 B3','S1 B4','S2 B1','S3 B1'});
for ih=1:7
set(BP1(ih,:),'LineWidth',2);
end
ylabel('T_{2F} (s)','FontSize',20); hline(T2_F,'k--');
get(gca, 'XTick'); set(gca, 'FontSize', 18); get(gca, 'YTick'); set(gca, 'FontSize', 18);
grid on; grid minor;
% M0F
subplot(3,2,5);
x = [Solution_Bouhrara(:,5);Solution_Deoni(:,5);Solution_Wood(:,5);Solution_Zhang(:,5);Solution_ASDeoni(:,5);Solution_ASRui(:,5)];
g = [ones(1000,1); 2*ones(1000,1); 3*ones(1000,1); 4*ones(1000,1); 5*ones(1000,1); 6*ones(1000,1)];
BP1 = boxplot(x,g,'Colors',cm,'Symbol','','Labels',{'S1 B1','S1 B2','S1 B3','S1 B4','S2 B1','S3 B1'});
for ih=1:7
set(BP1(ih,:),'LineWidth',2);
end
ylabel('MWF','FontSize',20); hline(M0_F,'k--');
get(gca, 'XTick'); set(gca, 'FontSize', 18); get(gca, 'YTick'); set(gca, 'FontSize', 18);
grid on; grid minor;
% kFS
if size(Solution_Bouhrara,2) == 6
subplot(3,2,6);
x = [Solution_Bouhrara(:,6);Solution_Deoni(:,6);Solution_Wood(:,6);Solution_Zhang(:,6);Solution_ASDeoni(:,6);Solution_ASRui(:,6)];
g = [ones(1000,1); 2*ones(1000,1); 3*ones(1000,1); 4*ones(1000,1); 5*ones(1000,1); 6*ones(1000,1)];
BP1 = boxplot(x,g,'Colors',cm,'Symbol','','Labels',{'S1 B1','S1 B2','S1 B3','S1 B4','S2 B1','S3 B1'});
for ih=1:7
set(BP1(ih,:),'LineWidth',2);
end
ylabel('k_{FS} (s^{-1})','FontSize',20); hline(k_FS,'k--');
get(gca, 'XTick'); set(gca, 'FontSize', 18); get(gca, 'YTick'); set(gca, 'FontSize', 18);
grid on; grid minor;
text(-5.366990291262136,156.5184396214944,0,'Monte Carlo Simulations Including Exchange (Model and Data)','FontSize',22,'FontWeight','Bold');
else
text(2.626168224299065,1.685245897770786,0,'Monte Carlo Simulations Excluding Exchange (Model and Data)','FontSize',22,'FontWeight','Bold');
end
%% Generates MC search space plots in Figure 9 and SF 6. %%
BWidth = 0.01; MaxAxis = 250; cm = lines(4);
figure(1); subplot(3,1,1);
histogram(M0F_Picked_T1,'BinWidth',BWidth,'EdgeColor', 'k', 'FaceColor', 'k', 'FaceAlpha', 0.5,'EdgeAlpha',0.5); hold on
get(gca, 'XTick'); set(gca, 'FontSize', 20); get(gca, 'YTick'); set(gca, 'FontSize', 20);
[N1,E1] = histcounts(Solution_T1_Bouhrara(:,5),'BinWidth',BWidth);
plot(E1(2:end),N1,'o--','Color',cm(1,:),'LineWidth',2)
[N2,E2] = histcounts(Solution_T1_Deoni(:,5),'BinWidth',BWidth);
plot(E2(2:end),N2,'o--','Color',cm(2,:),'LineWidth',2)
[N3,E3] = histcounts(Solution_T1_Wood(:,5),'BinWidth',BWidth);
plot(E3(2:end),N3,'o--','Color',cm(3,:),'LineWidth',2)
[N4,E4] = histcounts(Solution_T1_Zhang(:,5),'BinWidth',BWidth);
plot(E4(2:end),N4,'o--','Color',cm(4,:),'LineWidth',2)
ll = legend('Top 1000 Solutions','SRC + B1','SRC + B2','SRC + B3','SRC + B4');
ll.FontSize = 18; ll.AutoUpdate = 'off'; ll.Location = 'northwest'; legend('boxoff')
line([0.25 0.25],[0 MaxAxis],'LineStyle','--','Color','k','LineWidth',2)
line([0 0],[0 MaxAxis],'LineStyle',':','Color',cm(1,:),'LineWidth',2); line([0.5 0.5],[0 MaxAxis],'LineStyle',':','Color',cm(1,:),'LineWidth',2)
line([0 0],[0 MaxAxis],'LineStyle',':','Color',cm(2,:),'LineWidth',2); line([0.35 0.35],[0 MaxAxis],'LineStyle',':','Color',cm(2,:),'LineWidth',2)
line([0.001 0.001],[0 MaxAxis],'LineStyle',':','Color',cm(3,:),'LineWidth',2); line([0.35 0.35],[0 MaxAxis],'LineStyle',':','Color',cm(3,:),'LineWidth',2)
line([1e-7 1e-7],[0 MaxAxis],'LineStyle',':','Color',cm(4,:),'LineWidth',2); line([0.3 0.3],[0 MaxAxis],'LineStyle',':','Color',cm(4,:),'LineWidth',2)
tt = title('WML'); tt.FontSize = 24; tt.Position = [0.530991661062542,251.194408952268,1.4210854715202e-14]; grid on; grid minor;
xlabel('MWF','FontSize',20); ylabel('Count','FontSize', 20); xlim([0 0.55]); ylim([0 MaxAxis])
subplot(3,1,2);
histogram(M0F_Picked_T2,'BinWidth',BWidth, 'EdgeColor', 'k', 'FaceColor', 'k', 'FaceAlpha', 0.5,'EdgeAlpha',0.5); hold on
get(gca, 'XTick'); set(gca, 'FontSize', 20); get(gca, 'YTick'); set(gca, 'FontSize', 20);
[N4,E4] = histcounts(Solution_T2_Bouhrara(:,5),'BinWidth',BWidth);
plot(E4(2:end),N4,'o--','Color',cm(1,:),'LineWidth',2)
[N5,E5] = histcounts(Solution_T2_Deoni(:,5),'BinWidth',BWidth);
plot(E5(2:end),N5,'o--','Color',cm(2,:),'LineWidth',2)
[N6,E6] = histcounts(Solution_T2_Wood(:,5),'BinWidth',BWidth);
plot(E6(2:end),N6,'o--','Color',cm(3,:),'LineWidth',2)
[N7,E7] = histcounts(Solution_T2_Zhang(:,5),'BinWidth',BWidth);
plot(E7(2:end),N7,'o--','Color',cm(4,:),'LineWidth',2)
line([0.175 0.175],[0 MaxAxis],'LineStyle','--','Color','k','LineWidth',2)
line([0 0],[0 MaxAxis],'LineStyle',':','Color',cm(1,:),'LineWidth',2); line([0.5 0.5],[0 MaxAxis],'LineStyle',':','Color',cm(1,:),'LineWidth',2)
line([0 0],[0 MaxAxis],'LineStyle',':','Color',cm(2,:),'LineWidth',2); line([0.35 0.35],[0 MaxAxis],'LineStyle',':','Color',cm(2,:),'LineWidth',2)
line([0.001 0.001],[0 MaxAxis],'LineStyle',':','Color',cm(3,:),'LineWidth',2); line([0.35 0.35],[0 MaxAxis],'LineStyle',':','Color',cm(3,:),'LineWidth',2)
line([1e-7 1e-7],[0 MaxAxis],'LineStyle',':','Color',cm(4,:),'LineWidth',2); line([0.3 0.3],[0 MaxAxis],'LineStyle',':','Color',cm(4,:),'LineWidth',2)
tt = title('INT'); tt.FontSize = 24; tt.Position = [0.5358,250.0366442100638,1.4210854715202e-14]; grid on; grid minor
xlabel('MWF','FontSize',20); ylabel('Count','FontSize', 20); xlim([0 0.55]); ylim([0 MaxAxis])
subplot(3,1,3);
histogram(M0F_Picked_T3,'BinWidth',BWidth, 'EdgeColor', 'k', 'FaceColor', 'k', 'FaceAlpha', 0.5,'EdgeAlpha',0.5); hold on
get(gca, 'XTick'); set(gca, 'FontSize', 20); get(gca, 'YTick'); set(gca, 'FontSize', 20);
[N8,E8] = histcounts(Solution_T3_Bouhrara(:,5),'BinWidth',BWidth);
plot(E8(2:end),N8,'o--','Color',cm(1,:),'LineWidth',2)
[N9,E9] = histcounts(Solution_T3_Deoni(:,5),'BinWidth',BWidth);
plot(E9(2:end),N9,'o--','Color',cm(2,:),'LineWidth',2)
[N10,E10] = histcounts(Solution_T3_Wood(:,5),'BinWidth',BWidth);
plot(E10(2:end),N10,'o--','Color',cm(3,:),'LineWidth',2)
[N11,E11] = histcounts(Solution_T3_Zhang(:,5),'BinWidth',BWidth);
plot(E11(2:end),N11,'o--','Color',cm(4,:),'LineWidth',2)
line([0.1 0.1],[0 MaxAxis],'LineStyle','--','Color','k','LineWidth',2)
line([0 0],[0 MaxAxis],'LineStyle',':','Color',cm(1,:),'LineWidth',2); line([0.5 0.5],[0 MaxAxis],'LineStyle',':','Color',cm(1,:),'LineWidth',2)
line([0 0],[0 MaxAxis],'LineStyle',':','Color',cm(2,:),'LineWidth',2); line([0.35 0.35],[0 MaxAxis],'LineStyle',':','Color',cm(2,:),'LineWidth',2)
line([0.001 0.001],[0 MaxAxis],'LineStyle',':','Color',cm(3,:),'LineWidth',2); line([0.35 0.35],[0 MaxAxis],'LineStyle',':','Color',cm(3,:),'LineWidth',2)
line([1e-7 1e-7],[0 MaxAxis],'LineStyle',':','Color',cm(4,:),'LineWidth',2); line([0.3 0.3],[0 MaxAxis],'LineStyle',':','Color',cm(4,:),'LineWidth',2)
tt = title('GML'); tt.FontSize = 24; tt.Position = [0.532769872225958,250.4403658536585,1.4210854715202e-14]; grid on; grid minor
xlabel('MWF','FontSize',20); ylabel('Count','FontSize', 20); xlim([0 0.55]); ylim([0 MaxAxis])