promien

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% WERSJA 1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc;
x = [0, -170, -170, 0, 170, 170, 0];
y = [0, -80, 77, 170, 77, -80, -170];
f_c = 2.6; % [GHz]
N = -174;
N_W = 10^(N/10)
F=6;
IM=2;

xp = 0;
yp = 90;

d_2D = 100;

for i=1:7
    odleglosc_od_UT(i) = sqrt((xp - x(i))^2 + (yp - y(i))^2 );
end

for i=1:7
    d(i) = (odleglosc_od_UT(i)*d_2D) / odleglosc_od_UT(1);
    %PathL(i) = 28 + 22*log10(d(i)) + 20*log10(f_c); % -6*(h_UT-1.5);
    PathL(i) = 13.54 + 39.08*log10(d(i)) + 20*log10(f_c) - 9*log10(832^2 + 23^2);
    S(i) = 43 + 17 - 2 + 0 - 0 - F - IM - PathL(i);
end
S_0 = 10^(-0.50835); %% S z BS0
I_W(2) = 10^(-2.17598);
I_W(3) = 10^(-1.59271);
I_W(4) = 10^(-0.30845);
I_W(5) = 10^(-1.59271);
I_W(6) = 10^(-2.17598);
I_W(7) = 10^(-2.30889);

I = 0;
for k=2:7
    I = I + I_W(k);
end

%I = I_W(2) + I_W(3) + I_W(4) + I_W(5) + I_W(6) + I_W(7);

SINR =  S_0/(I + N_W);
SINR_dB = 10*log10(S_0/(I + N_W));


odleglosc_od_UT
d
PathL
S
S_0
I_W
I
SINR
SINR_dB

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% WERSJA 2%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clc;
x = [0, -170, -170, 0, 170, 170, 0];
y = [0, -80, 77, 170, 77, -80, -170];
f_c = 2.6; % [GHz]
N = -174;
N_W = 10^(N/10)
F=6;
IM=2;

xp = 0;
yp = 90;

d_2D = 100;

for i=1:7
    odleglosc_od_UT(i) = sqrt((xp - x(i))^2 + (yp - y(i))^2 );
end

for i=1:7
    d(i) = (odleglosc_od_UT(i)*d_2D) / odleglosc_od_UT(1);
    %PathL(i) = 28 + 22*log10(d(i)) + 20*log10(f_c); % -6*(h_UT-1.5);
    %PathL(i) = 13.54 + 39.08*log10(d(i)) + 20*log10(f_c +0.1*i) - 9*log10(832^2 + 23^2);
    PathL(i) = 13.54 + 39.08*log10(d(i)) + 20*log10(f_c) - 9*log10(312^2 + 23^2);
    S(i) = 43 + 17 - 2 + 0 - 0 - F - IM - PathL(i);
end
S_0 = 10^(S(1)/10); %% S z BS0

for l=2:7
    I_W(l) = 10^(S(l)/10);
end

I = 0;
for k=2:7
    I = I + I_W(k);
end

SINR =  S_0/(I + N_W);
SINR_dB = 10*log10(S_0/(I + N_W));


odleglosc_od_UT
d
PathL
S
S_0
I_W
I
SINR
SINR_dB