From 8f965c14663176c0dfbebc3d784ae2dfe575e7ea Mon Sep 17 00:00:00 2001 From: Ksyer Date: Sat, 27 Jan 2024 20:57:40 +0800 Subject: [PATCH] =?UTF-8?q?Add=20=E2=80=9CRandomized=20Stepped=20FR=20Expl?= =?UTF-8?q?oiting=20BS=20of=20ET=E2=80=9D?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- .../test.m | 17 +++++++++++++++++ 1 file changed, 17 insertions(+) create mode 100644 Randomized Stepped FR Exploiting BS of ET/test.m diff --git a/Randomized Stepped FR Exploiting BS of ET/test.m b/Randomized Stepped FR Exploiting BS of ET/test.m new file mode 100644 index 0000000..6c541c1 --- /dev/null +++ b/Randomized Stepped FR Exploiting BS of ET/test.m @@ -0,0 +1,17 @@ +N = 32e9; +M = 3; +eps = 1e-4; + +delta_1 = 24 * sqrt((M-1)/N) * log(M*N) * (2*sqrt(log(M * N) - log(eps)) + 1); +delta_2 = 3/2 * sqrt((M-1)/N) * (2*sqrt(log(M * 2) - log(eps)) + 1); + +K = N * (1/8 - delta_1 - delta_2)^2 / (81 * M * log(M * N) * (1 + 2/3 * delta_2)); + +x = 0:25; +Ns = N * (1+randn(size(x))); +delta_1 = 24 * sqrt((M-1)./N) * log(M.*N) * (2*sqrt(log(M .* N) - log(eps)) + 1); +delta_2 = 3/2 * sqrt((M-1)./N) * (2*sqrt(log(M * 2) - log(eps)) + 1); +Ks = Ns * (1/8 - delta_1 - delta_2)^2 / (81 * M * log(M .* Ns) * (1 + 2/3 * delta_2)); + +rate = K .* M .* log(M .* Ns) ./ Ns; +plot(x, rate)