wyznaczanie_danych

clc;
close all;
format compact
format short

x = [470, 323, 323, 470, 617, 617, 470, 323, 176, 176, 176, 323, 470, 617, 764, 764, 764, 617, 470];
y = [400, 485, 315, 230, 315, 485, 570, 655, 570, 400, 230, 145, 60, 145, 230, 400, 570, 655, 740];
f_c = 4; % [GHz]
N_W = 1.38*10^(-23) * 300 * 10*10^6* 10^(5/10);
F=5;
IM=2;
h_BS= 3;
h_UT= 1.5;
h = 5;
W = 5;
xp = zeros(1, 85);
yp = zeros(1, 85);
SINR = zeros(1, 85);
%SINR_NLOS = zeros(1, 85);
%SINR_LOS = zeros(1, 85);
%SINR_dB_NLOS = zeros(1, 85);
%SINR_dB_LOS = zeros(1, 85);
SINR_dB = zeros(1, 85);

x_UT = 550;
y_UT = 250;
d_2D = 15; %1000, 650, 1000, 5000
t=1;

xp(1) = 470;
yp(1) = 400;

for  j = 2:85
xp(j) = xp(j-1) + 1;
yp(j) = (y_UT - y(1))./(x_UT - x(1))*(xp(j) - x(1)) + y(1);
end



for w=1:85
    for i=1:19
     odleglosc_od_UT(i) = sqrt((xp(w) - x(i)).^2 + (yp(w) - y(i)).^2 );
    end

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


    % InH - Indoor Hotspot
    PathL1(i) = 32.4 + 17.3*log10(d(i, w)) + 20*log10(f_c);
    PathL2(i) = 17.3 + 38.3*log10(d(i, w)) + 24.9*log10(f_c);
    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) = 3 + 7 - 2 + 0 - 0 - PathL(i);
    end
    S_0 = 10^(S(1)/10); %% S z BS0

format short
    I = 0;
    for l=2:19
        I_W(l) = 10^(S(l)/10);
        I = I + I_W(l);
    end

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

end

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


%f2=figure    UMa
%plot([1:85], 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