clear all
clc

% Переменные по заданию
R = 2000;
P_tx = 320;
f_0 = 1800;
delta_f = 5*10^6;
k_n = 3;
% Окумура-Хата, большой город

% Константы
T = 300;
k = 1.38*(10^-23);

% Переменные
N = 6;
h_BS = 10 % БС всегда на одном и том же уровне
AB_d = randi([1 10], 1, N)
for i = 1:N
    AB_h_RX(i) = randi([1 h_BS]);
end
AB_h_RX
S = 3; % Т.к. город большой

% Рассчитываем мощьность шума
P_n = T*k*k_n*delta_f;

% Рассчитываем потери L
a = (1.1*log10(f_0) - 0.7).*AB_h_RX - (1.56*log10(f_0) - 0.8);
L = 46.3 + 33.9*log10(f_0) - 13.82*log10(h_BS) - a + (44.9 - 6.55.*log10(AB_h_RX)).*log10(AB_d) + S; % В децибеллах
L = 10.^(L./10); % Разы

% Рассчитываем мощьность, принятую АБ
P_RX = P_tx./L;

% Рассчтываем отношения сигнал/шум
SNR = P_RX./P_n;

% Рассчитываем пропускную способность канала связи
C = delta_f.*log2(1 + SNR);
disp('при вышеуканных условиях пропускная сопопсобность в Мбит/с равна ')
(C./1024^2)