Files
FAR_CS/Block RIP/bound/get_wl_ub.m
T
2024-02-21 16:44:41 +08:00

14 lines
392 B
Matlab

function result = get_wl_ub(N, M, eps)
k1 = sqrt((M - 1) / N);
k2 = log(M * N);
k3 = sqrt(k2 - log(eps));
k4 = sqrt(log(2 * M) - log(eps));
delta_1 = 24 * k1 * k2 * (2 * k3 + 1);
delta_2 = 1.5 * k1 * (2 * k4 + 1);
numerator = N * (1/8 - delta_1 - delta_2) ^ 2;
denominator = 81 * M * k2 * (1 + 2 * delta_2 / 3);
result = numerator / denominator;
end