-
Notifications
You must be signed in to change notification settings - Fork 9
/
diff_rank_rho_lrr.m
55 lines (35 loc) · 1.24 KB
/
diff_rank_rho_lrr.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
clear;
config = diff_rank_rho_config();
method = 'lrr';
all_d = config.d;
all_rho = config.rho;
all_rep = config.repetitions;
K = config.K;
p = config.p;
n = config.n;
epochs = config.epoch;
%% add path
addpath('LRR/');
%% compute EV for LRR
for k=1:length(all_rep)
rep = all_rep(k);
for j=1:length(all_d)
d = all_d(j);
for i=1:length(all_rho)
rho = all_rho(i);
fprintf('LRR: rep = %d, d = %d, rho = %.2f\n', rep, d, rho);
result_file = sprintf(config.result_file_format, method, method, d, rho, rep);
data_file = sprintf(config.data_file_format, d, rho, rep);
load(data_file);
U_all = cell2mat(U);
UUt = U_all * U_all';
traceUUt = trace(UUt);
[X, ~] = solve_lrr(Z, Z, 1/sqrt(n), 1, 1);
[D, ~, ~] = svds(Z*X, d);
ev = js_compute_EV(D, UUt, traceUUt);
EV = ev * ones(1, n);
save(result_file, 'EV');
fprintf('EV = %g, save to %s\n', ev, result_file);
end
end
end