-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathPaperRevision_Code_Prep.m
287 lines (254 loc) · 14.8 KB
/
PaperRevision_Code_Prep.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
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
%%% Test stability of SRC estimates for different no. realisations. %%%
close all; clear all;
load('Fig5_Data.mat')
Bouhrara_Mean = mean(Solution_Bouhrara);
Bouhrara_Mean_US900 = mean(Solution_Bouhrara(1:900,:));
Bouhrara_Mean_US800 = mean(Solution_Bouhrara(1:800,:));
Bouhrara_Mean_US700 = mean(Solution_Bouhrara(1:700,:));
Bouhrara_Mean_US600 = mean(Solution_Bouhrara(1:600,:));
Bouhrara_Mean_US500 = mean(Solution_Bouhrara(1:500,:));
Bouhrara_Mean_US400 = mean(Solution_Bouhrara(1:400,:));
Bouhrara_Mean_US300 = mean(Solution_Bouhrara(1:300,:));
Bouhrara_Mean_US200 = mean(Solution_Bouhrara(1:200,:));
Bouhrara_Mean_US100 = mean(Solution_Bouhrara(1:100,:));
Bouhrara_SD = std(Solution_Bouhrara);
Bouhrara_SD_US900 = std(Solution_Bouhrara(1:900,:));
Bouhrara_SD_US800 = std(Solution_Bouhrara(1:800,:));
Bouhrara_SD_US700 = std(Solution_Bouhrara(1:700,:));
Bouhrara_SD_US600 = std(Solution_Bouhrara(1:600,:));
Bouhrara_SD_US500 = std(Solution_Bouhrara(1:500,:));
Bouhrara_SD_US400 = std(Solution_Bouhrara(1:400,:));
Bouhrara_SD_US300 = std(Solution_Bouhrara(1:300,:));
Bouhrara_SD_US200 = std(Solution_Bouhrara(1:200,:));
Bouhrara_SD_US100 = std(Solution_Bouhrara(1:100,:));
Bouhrara_MeanVector = [Bouhrara_Mean_US100;Bouhrara_Mean_US200;Bouhrara_Mean_US300;Bouhrara_Mean_US400;Bouhrara_Mean_US500;Bouhrara_Mean_US600;Bouhrara_Mean_US700;Bouhrara_Mean_US800;Bouhrara_Mean_US900;Bouhrara_Mean];
Bouhrara_SDVector = [Bouhrara_SD_US100;Bouhrara_SD_US200;Bouhrara_SD_US300;Bouhrara_SD_US400;Bouhrara_SD_US500;Bouhrara_SD_US600;Bouhrara_SD_US700;Bouhrara_SD_US800;Bouhrara_SD_US900;Bouhrara_SD];
Deoni_Mean = mean(Solution_Deoni);
Deoni_Mean_US900 = mean(Solution_Deoni(1:900,:));
Deoni_Mean_US800 = mean(Solution_Deoni(1:800,:));
Deoni_Mean_US700 = mean(Solution_Deoni(1:700,:));
Deoni_Mean_US600 = mean(Solution_Deoni(1:600,:));
Deoni_Mean_US500 = mean(Solution_Deoni(1:500,:));
Deoni_Mean_US400 = mean(Solution_Deoni(1:400,:));
Deoni_Mean_US300 = mean(Solution_Deoni(1:300,:));
Deoni_Mean_US200 = mean(Solution_Deoni(1:200,:));
Deoni_Mean_US100 = mean(Solution_Deoni(1:100,:));
Deoni_SD = std(Solution_Deoni);
Deoni_SD_US900 = std(Solution_Deoni(1:900,:));
Deoni_SD_US800 = std(Solution_Deoni(1:800,:));
Deoni_SD_US700 = std(Solution_Deoni(1:700,:));
Deoni_SD_US600 = std(Solution_Deoni(1:600,:));
Deoni_SD_US500 = std(Solution_Deoni(1:500,:));
Deoni_SD_US400 = std(Solution_Deoni(1:400,:));
Deoni_SD_US300 = std(Solution_Deoni(1:300,:));
Deoni_SD_US200 = std(Solution_Deoni(1:200,:));
Deoni_SD_US100 = std(Solution_Deoni(1:100,:));
Deoni_MeanVector = [Deoni_Mean_US100;Deoni_Mean_US200;Deoni_Mean_US300;Deoni_Mean_US400;Deoni_Mean_US500;Deoni_Mean_US600;Deoni_Mean_US700;Deoni_Mean_US800;Deoni_Mean_US900;Deoni_Mean];
Deoni_SDVector = [Deoni_SD_US100;Deoni_SD_US200;Deoni_SD_US300;Deoni_SD_US400;Deoni_SD_US500;Deoni_SD_US600;Deoni_SD_US700;Deoni_SD_US800;Deoni_SD_US900;Deoni_SD];
Wood_Mean = mean(Solution_Wood);
Wood_Mean_US900 = mean(Solution_Wood(1:900,:));
Wood_Mean_US800 = mean(Solution_Wood(1:800,:));
Wood_Mean_US700 = mean(Solution_Wood(1:700,:));
Wood_Mean_US600 = mean(Solution_Wood(1:600,:));
Wood_Mean_US500 = mean(Solution_Wood(1:500,:));
Wood_Mean_US400 = mean(Solution_Wood(1:400,:));
Wood_Mean_US300 = mean(Solution_Wood(1:300,:));
Wood_Mean_US200 = mean(Solution_Wood(1:200,:));
Wood_Mean_US100 = mean(Solution_Wood(1:100,:));
Wood_SD = std(Solution_Wood);
Wood_SD_US900 = std(Solution_Wood(1:900,:));
Wood_SD_US800 = std(Solution_Wood(1:800,:));
Wood_SD_US700 = std(Solution_Wood(1:700,:));
Wood_SD_US600 = std(Solution_Wood(1:600,:));
Wood_SD_US500 = std(Solution_Wood(1:500,:));
Wood_SD_US400 = std(Solution_Wood(1:400,:));
Wood_SD_US300 = std(Solution_Wood(1:300,:));
Wood_SD_US200 = std(Solution_Wood(1:200,:));
Wood_SD_US100 = std(Solution_Wood(1:100,:));
Wood_MeanVector = [Wood_Mean_US100;Wood_Mean_US200;Wood_Mean_US300;Wood_Mean_US400;Wood_Mean_US500;Wood_Mean_US600;Wood_Mean_US700;Wood_Mean_US800;Wood_Mean_US900;Wood_Mean];
Wood_SDVector = [Wood_SD_US100;Wood_SD_US200;Wood_SD_US300;Wood_SD_US400;Wood_SD_US500;Wood_SD_US600;Wood_SD_US700;Wood_SD_US800;Wood_SD_US900;Wood_SD];
Zhang_Mean = mean(Solution_Zhang);
Zhang_Mean_US900 = mean(Solution_Zhang(1:900,:));
Zhang_Mean_US800 = mean(Solution_Zhang(1:800,:));
Zhang_Mean_US700 = mean(Solution_Zhang(1:700,:));
Zhang_Mean_US600 = mean(Solution_Zhang(1:600,:));
Zhang_Mean_US500 = mean(Solution_Zhang(1:500,:));
Zhang_Mean_US400 = mean(Solution_Zhang(1:400,:));
Zhang_Mean_US300 = mean(Solution_Zhang(1:300,:));
Zhang_Mean_US200 = mean(Solution_Zhang(1:200,:));
Zhang_Mean_US100 = mean(Solution_Zhang(1:100,:));
Zhang_SD = std(Solution_Zhang);
Zhang_SD_US900 = std(Solution_Zhang(1:900,:));
Zhang_SD_US800 = std(Solution_Zhang(1:800,:));
Zhang_SD_US700 = std(Solution_Zhang(1:700,:));
Zhang_SD_US600 = std(Solution_Zhang(1:600,:));
Zhang_SD_US500 = std(Solution_Zhang(1:500,:));
Zhang_SD_US400 = std(Solution_Zhang(1:400,:));
Zhang_SD_US300 = std(Solution_Zhang(1:300,:));
Zhang_SD_US200 = std(Solution_Zhang(1:200,:));
Zhang_SD_US100 = std(Solution_Zhang(1:100,:));
Zhang_MeanVector = [Zhang_Mean_US100;Zhang_Mean_US200;Zhang_Mean_US300;Zhang_Mean_US400;Zhang_Mean_US500;Zhang_Mean_US600;Zhang_Mean_US700;Zhang_Mean_US800;Zhang_Mean_US900;Zhang_Mean];
Zhang_SDVector = [Zhang_SD_US100;Zhang_SD_US200;Zhang_SD_US300;Zhang_SD_US400;Zhang_SD_US500;Zhang_SD_US600;Zhang_SD_US700;Zhang_SD_US800;Zhang_SD_US900;Zhang_SD];
ASDeoni_Mean = mean(Solution_ASDeoni);
ASDeoni_Mean_US900 = mean(Solution_ASDeoni(1:900,:));
ASDeoni_Mean_US800 = mean(Solution_ASDeoni(1:800,:));
ASDeoni_Mean_US700 = mean(Solution_ASDeoni(1:700,:));
ASDeoni_Mean_US600 = mean(Solution_ASDeoni(1:600,:));
ASDeoni_Mean_US500 = mean(Solution_ASDeoni(1:500,:));
ASDeoni_Mean_US400 = mean(Solution_ASDeoni(1:400,:));
ASDeoni_Mean_US300 = mean(Solution_ASDeoni(1:300,:));
ASDeoni_Mean_US200 = mean(Solution_ASDeoni(1:200,:));
ASDeoni_Mean_US100 = mean(Solution_ASDeoni(1:100,:));
ASDeoni_SD = std(Solution_ASDeoni);
ASDeoni_SD_US900 = std(Solution_ASDeoni(1:900,:));
ASDeoni_SD_US800 = std(Solution_ASDeoni(1:800,:));
ASDeoni_SD_US700 = std(Solution_ASDeoni(1:700,:));
ASDeoni_SD_US600 = std(Solution_ASDeoni(1:600,:));
ASDeoni_SD_US500 = std(Solution_ASDeoni(1:500,:));
ASDeoni_SD_US400 = std(Solution_ASDeoni(1:400,:));
ASDeoni_SD_US300 = std(Solution_ASDeoni(1:300,:));
ASDeoni_SD_US200 = std(Solution_ASDeoni(1:200,:));
ASDeoni_SD_US100 = std(Solution_ASDeoni(1:100,:));
ASDeoni_MeanVector = [ASDeoni_Mean_US100;ASDeoni_Mean_US200;ASDeoni_Mean_US300;ASDeoni_Mean_US400;ASDeoni_Mean_US500;ASDeoni_Mean_US600;ASDeoni_Mean_US700;ASDeoni_Mean_US800;ASDeoni_Mean_US900;ASDeoni_Mean];
ASDeoni_SDVector = [ASDeoni_SD_US100;ASDeoni_SD_US200;ASDeoni_SD_US300;ASDeoni_SD_US400;ASDeoni_SD_US500;ASDeoni_SD_US600;ASDeoni_SD_US700;ASDeoni_SD_US800;ASDeoni_SD_US900;ASDeoni_SD];
ASRui_Mean = mean(Solution_ASRui);
ASRui_Mean_US900 = mean(Solution_ASRui(1:900,:));
ASRui_Mean_US800 = mean(Solution_ASRui(1:800,:));
ASRui_Mean_US700 = mean(Solution_ASRui(1:700,:));
ASRui_Mean_US600 = mean(Solution_ASRui(1:600,:));
ASRui_Mean_US500 = mean(Solution_ASRui(1:500,:));
ASRui_Mean_US400 = mean(Solution_ASRui(1:400,:));
ASRui_Mean_US300 = mean(Solution_ASRui(1:300,:));
ASRui_Mean_US200 = mean(Solution_ASRui(1:200,:));
ASRui_Mean_US100 = mean(Solution_ASRui(1:100,:));
ASRui_SD = std(Solution_ASRui);
ASRui_SD_US900 = std(Solution_ASRui(1:900,:));
ASRui_SD_US800 = std(Solution_ASRui(1:800,:));
ASRui_SD_US700 = std(Solution_ASRui(1:700,:));
ASRui_SD_US600 = std(Solution_ASRui(1:600,:));
ASRui_SD_US500 = std(Solution_ASRui(1:500,:));
ASRui_SD_US400 = std(Solution_ASRui(1:400,:));
ASRui_SD_US300 = std(Solution_ASRui(1:300,:));
ASRui_SD_US200 = std(Solution_ASRui(1:200,:));
ASRui_SD_US100 = std(Solution_ASRui(1:100,:));
ASRui_MeanVector = [ASRui_Mean_US100;ASRui_Mean_US200;ASRui_Mean_US300;ASRui_Mean_US400;ASRui_Mean_US500;ASRui_Mean_US600;ASRui_Mean_US700;ASRui_Mean_US800;ASRui_Mean_US900;ASRui_Mean];
ASRui_SDVector = [ASRui_SD_US100;ASRui_SD_US200;ASRui_SD_US300;ASRui_SD_US400;ASRui_SD_US500;ASRui_SD_US600;ASRui_SD_US700;ASRui_SD_US800;ASRui_SD_US900;ASRui_SD];
USFactors = 100:100:1000;
figure(1); subplot(3,2,1);
plot(USFactors,Bouhrara_MeanVector(:,1),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_MeanVector(:,1),'-o','LineWidth',2);
plot(USFactors,Wood_MeanVector(:,1),'-o','LineWidth',2);
plot(USFactors,Zhang_MeanVector(:,1),'-o','LineWidth',2);
plot(USFactors,ASDeoni_MeanVector(:,1),'-o','LineWidth',2);
plot(USFactors,ASRui_MeanVector(:,1),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16); ylabel('Mean T_{1S} (s)','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])
ll = legend('S1+B1','S1+B2','S1+B3','S1+B4','S2+B1','S3+B1'); ll.FontSize = 16; legend 'boxoff'; ll.Orientation = 'horizontal';
subplot(3,2,2);
plot(USFactors,Bouhrara_MeanVector(:,2),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_MeanVector(:,2),'-o','LineWidth',2);
plot(USFactors,Wood_MeanVector(:,2),'-o','LineWidth',2);
plot(USFactors,Zhang_MeanVector(:,2),'-o','LineWidth',2);
plot(USFactors,ASDeoni_MeanVector(:,2),'-o','LineWidth',2);
plot(USFactors,ASRui_MeanVector(:,2),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16); ylabel('Mean T_{1F} (s)','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])
subplot(3,2,3);
plot(USFactors,Bouhrara_MeanVector(:,3),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_MeanVector(:,3),'-o','LineWidth',2);
plot(USFactors,Wood_MeanVector(:,3),'-o','LineWidth',2);
plot(USFactors,Zhang_MeanVector(:,3),'-o','LineWidth',2);
plot(USFactors,ASDeoni_MeanVector(:,3),'-o','LineWidth',2);
plot(USFactors,ASRui_MeanVector(:,3),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16); ylabel('Mean T_{2S} (s)','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])
subplot(3,2,4);
plot(USFactors,Bouhrara_MeanVector(:,4),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_MeanVector(:,4),'-o','LineWidth',2);
plot(USFactors,Wood_MeanVector(:,4),'-o','LineWidth',2);
plot(USFactors,Zhang_MeanVector(:,4),'-o','LineWidth',2);
plot(USFactors,ASDeoni_MeanVector(:,4),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16); ylabel('Mean T_{2F} (s)','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])
subplot(3,2,5);
plot(USFactors,Bouhrara_MeanVector(:,5),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_MeanVector(:,5),'-o','LineWidth',2);
plot(USFactors,Wood_MeanVector(:,5),'-o','LineWidth',2);
plot(USFactors,Zhang_MeanVector(:,5),'-o','LineWidth',2);
plot(USFactors,ASDeoni_MeanVector(:,5),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16);ylabel('Mean MWF','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])
subplot(3,2,6);
plot(USFactors,Bouhrara_MeanVector(:,6),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_MeanVector(:,6),'-o','LineWidth',2);
plot(USFactors,Wood_MeanVector(:,6),'-o','LineWidth',2);
plot(USFactors,Zhang_MeanVector(:,6),'-o','LineWidth',2);
plot(USFactors,ASDeoni_MeanVector(:,6),'-o','LineWidth',2);
plot(USFactors,ASRui_MeanVector(:,6),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16); ylabel('Mean k_{FS} (s^{-1})','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])
figure(2); subplot(3,2,1);
plot(USFactors,Bouhrara_SDVector(:,1),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_SDVector(:,1),'-o','LineWidth',2);
plot(USFactors,Wood_SDVector(:,1),'-o','LineWidth',2);
plot(USFactors,Zhang_SDVector(:,1),'-o','LineWidth',2);
plot(USFactors,ASDeoni_SDVector(:,1),'-o','LineWidth',2);
plot(USFactors,ASRui_SDVector(:,1),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16); ylabel('SD T_{1S} (s)','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])
ll = legend('S1+B1','S1+B2','S1+B3','S1+B4','S2+B1','S3+B1'); ll.FontSize = 16; legend 'boxoff'; ll.Orientation = 'horizontal';
subplot(3,2,2);
plot(USFactors,Bouhrara_SDVector(:,2),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_SDVector(:,2),'-o','LineWidth',2);
plot(USFactors,Wood_SDVector(:,2),'-o','LineWidth',2);
plot(USFactors,Zhang_SDVector(:,2),'-o','LineWidth',2);
plot(USFactors,ASDeoni_SDVector(:,2),'-o','LineWidth',2);
plot(USFactors,ASRui_SDVector(:,2),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16); ylabel('SD T_{1F} (s)','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])
subplot(3,2,3);
plot(USFactors,Bouhrara_SDVector(:,3),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_SDVector(:,3),'-o','LineWidth',2);
plot(USFactors,Wood_SDVector(:,3),'-o','LineWidth',2);
plot(USFactors,Zhang_SDVector(:,3),'-o','LineWidth',2);
plot(USFactors,ASDeoni_SDVector(:,3),'-o','LineWidth',2);
plot(USFactors,ASRui_SDVector(:,3),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16); ylabel('SD T_{2S} (s)','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])
subplot(3,2,4);
plot(USFactors,Bouhrara_SDVector(:,4),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_SDVector(:,4),'-o','LineWidth',2);
plot(USFactors,Wood_SDVector(:,4),'-o','LineWidth',2);
plot(USFactors,Zhang_SDVector(:,4),'-o','LineWidth',2);
plot(USFactors,ASDeoni_SDVector(:,4),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16); ylabel('SD T_{2F} (s)','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])
subplot(3,2,5);
plot(USFactors,Bouhrara_SDVector(:,5),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_SDVector(:,5),'-o','LineWidth',2);
plot(USFactors,Wood_SDVector(:,5),'-o','LineWidth',2);
plot(USFactors,Zhang_SDVector(:,5),'-o','LineWidth',2);
plot(USFactors,ASDeoni_SDVector(:,5),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16);ylabel('SD MWF','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])
subplot(3,2,6);
plot(USFactors,Bouhrara_SDVector(:,6),'-o','LineWidth',2); hold on;
plot(USFactors,Deoni_SDVector(:,6),'-o','LineWidth',2);
plot(USFactors,Wood_SDVector(:,6),'-o','LineWidth',2);
plot(USFactors,Zhang_SDVector(:,6),'-o','LineWidth',2);
plot(USFactors,ASDeoni_SDVector(:,6),'-o','LineWidth',2);
plot(USFactors,ASRui_SDVector(:,6),'-o','LineWidth',2);
xlabel('No. Realisations','FontSize',16); ylabel('SD k_{FS} (s^{-1})','FontSize',16);
get(gca, 'XTick'); set(gca, 'FontSize', 16); get(gca, 'YTick'); set(gca, 'FontSize', 16);
grid on; grid minor; xlim([100 1000])