sparse/bVlogEp.m at master · egoolish/sparse · GitHub

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
epsilonrange_1 = 0.01:0.01:0.05;
epsilonrange_2 = 0.2:0.2:1;
epsilonrange = horzcat(epsilonrange_1, 0.1, epsilonrange_2);
e_values = log(1./epsilonrange);
b_values = [0.1393, -0.1345, -0.1060, -0.1773, -0.3847, -0.4060, -0.6948,...
    -0.9499, -1.1234, -1.2559, -1.4674];


f = figure;
poly = polyfit(b_values, e_values, 1);
a_value = poly(1);
b_value = poly(2);
plot(b_values, e_values, 'DisplayName', 'N^3');
xlabel(['b vs log(1/epsilon) ', 'a=', num2str(a_value), ' b=', num2str(b_value)]);
legend('show');

% N^3
% b_values = [-0.2918, -0.2512, -0.3364, -0.3416, -0.4932, -0.4391, -0.4858, -0.4031,...
%     -0.4513, -0.4993, -0.4661, -0.6927, -0.8072, -0.9505, -1.0359, -1.1057, -1.2115,...
%     -1.2345, -1.3123, -1.4582];


%N^3Rand
% b_values = [0.0126, -0.0615, -0.1006, -0.1832, -0.2561, -0.4351, -0.5817,...
%     -0.8467, -1.1029, -1.1915, -1.4589];

%N^2RandRun1
% b_values = [0.1253, 0.0415, -0.0126, -0.1270, -0.1431, -0.4428, -0.7620,...
%     -1.0066, -1.1989, -1.3515];

%N^2RandRun2
% b_values = [0.1175, -0.0283, -0.0931, -0.1411, -0.2316, -0.4064,...
%     -0.5421, -0.8580, -1.0679, -1.2944, -1.4332];

%N^2
% b_values = [-0.1162, 0.0696, -0.0024, -0.1308, -0.1480, -0.1498, -0.1718,...
%     -0.2993, -0.2478, -0.2121, -0.2964, -0.3644, -0.2800, -0.3818, -0.3842, -0.3906, -0.3788,...
%     -0.4235, -0.3159, -0.5044, -0.6101, -0.7691, -0.8642, -0.9211, -1.0702, -1.1906, -1.2148, -1.3634, -1.3547];