wykresy

clc;
format compact
format short

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-3);
F=6;
IM=2;
h_BS= 6;
h_UT= 2;
h = 6;
W = 20;
SINR_dB = zeros(1, 140);

xp = 0;
yp = 85;

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

for i=1:7
    d2D(i) = (odleglosc_od_UT(i)*d_2D) / odleglosc_od_UT(1);
    d(i) = sqrt(d2D(i).^2 + (h_BS - h_UT)^2);


    % InH - Indoor Hotspot
    %PathL1(i) = 32.4 + 17.3*log10(d(i)) + 20*log10(f_c +0.1*i);
    %PathL2(i) = 17.3 + 38.3*log10(d(i)) + 24.9*log10(f_c +0.1*i);
    %PathL(i) = max(PathL1(i), PathL2(i));


    % UMa- Urban Macro, d_BP1 = 832
    %PathL1(i) = 28 + 40*log10(d(i)) + 20*log10(f_c +0.1*i) - 9*log10(832^2 + 23^2); % d_BP >832
    %PathL2(i) = 13.54 + 39.08*log10(d(i)) + 20*log10(f_c +0.1*i) - 6*(h_UT-1.5);
    %PathL(i) = max(PathL1(i), PathL2(i));

    % UMi- Urban Micro, d_BP1 = 832
    PathL1(i) = 32.4 + 21*log10(d(i)) +20*log10(f_c +0.1*i);  % d_BP < 832
    %PathL1(i) = 32.4 + 21*log10(d(i)) +20*log10(f_c +0.1*i)- 9.5*log10(832^2 + 23^2);  % d_BP < 832
    PathL2(i) = 22.4 + 35.3*log10(d(i)) + 21.3*log10(f_c +0.1*i) -0.3*(h_UT -1.5); % d_BP > 832
    PathL(i) = max(PathL1(i), PathL2(i));

    % RMa- Rural Macro, d_BP2 = 2722,7
    %PathL1(i) = 20*log10((40*pi*d(i)*(f_c +0.1*i))/3) + min((0.03*h^(1.72)), 10)*log10(d(i))- min((0.044*h^(1.72)), 14.77)+0.002*log10(h)*d(i) ;  % d_BP < 2722,7
    %PathL1(i) =  20*log10((40*pi*2722.7*(f_c +0.1*i))/3) + min((0.03*h^(1.72)), 10)*log10(2722.7)- min((0.044*h^(1.72)), 14.77)+0.002*log10(h)*2722.7 + 40*log10(d(i)/2722.7); % d_BP > 2722.7
    %PathL2(i) = 22.4 + 35.3*log10(d(i)) + 21.3*log10(f_c +0.1*i) -0.3*(h_UT -1.5); % d_BP > 2722,7
    %PathL(i) = max(PathL1(i), PathL2(i));


    S(i) = 43 + 17 - 2 + 0 - 0 - F - IM - PathL(i);
end
S_0 = 10^(S(1)/10 - 3); %% S z BS0


format short
for l=2:7
    I_W(l) = 10^(S(l)/10 - 3);
end

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

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

t=t+1;
end


%f1=figure    InH
%x=[1:140];
%plot(x, SINR_dB, 'k-');
%title("Wykres zależności poziomu SINR od odległości UT od stacji bazowej dla modelu InH");
%axis([0 100 0 85]);
%xlabel("d [m]");
%ylabel("SINR [dB]");


%f2=figure    UMa
%x=[1:10:1000];
%plot(x, SINR_dB, 'k-');
%title("Wykres zależności poziomu SINR od odległości UT od stacji bazowej dla modelu UMa");
%axis([0 1000 0 85]);
%xlabel("d [m]");
%ylabel("SINR [dB]");

%f3=figure    UMi
x=[1:5:700];
plot(x, SINR_dB, 'k-');
title("Wykres zależności poziomu SINR od odległości UT od stacji bazowej dla modelu UMi");
axis([0 450 0 85]);
xlabel("d [m]");
ylabel("SINR [dB]");

%f4=figure,   RMa
%x=[1:25:5000];
%plot(x, SINR_dB, 'k-');
%title("Wykres zależności poziomu SINR od odległości UT od stacji bazowej dla modelu RMa");
%axis([0 1000 0 70]);
%xlabel("d [m]");
%ylabel("SINR [dB]");
%subplot(); wiele wykresów w jednym oknie
%saveas(gcf, "RMa_wykres.png")


odleglosc_od_UT
d2D
d
PathL
S
S_0
I_W
I
SINR
SINR_dB