Untitled

>> Pilawski_Skomski_SAGN_skrypt_matlab_projekt1
-------------------------------------------WYNIKI--------------------------------------------------------
-------------------------------------------Obliczenie NINi oraz NDSP_INi----------------------------------
1. NINi przy zalozeniu dopuszczalnego prawdopodobienstwa straty B=0.002
 Brama(GWi)  |  Ai[Erl]  |  NINi
      1      |    150    |  180
      2      |    330    |  371
      3      |    410    |  454


 2. Przeliczone NINi na NPCM_30/32 dla kierunku IN
 Brama(GWi)  |  NPCM_30/32_INi
      1      |    6
      2      |    13
      3      |    16


 3. Rzeczywista liczba NR_INi laczy w kierunku In
 Brama(GWi)  |  NR_INi]
      1      |    180
      2      |    390
      3      |    480


 4. Liczba DSP dla kierunku IN przy zalozeniu n_dsp = 4
 Brama(GWi)  |  NDSP_INi
      1      |    45
      2      |    98
      3      |    120

-------------------------------------------Obliczenie NOUTi raz NDSP_OUTi----------------------------------

 1. Obliczenia Aout,PSTN_d
 Brama(GWi)  |  lambda_RT  |  lambda_G729  |  AOUT_PSTN_d
      1      |    11000    |    40        |  287.00
      2      |    11100    |    40        |  293.00
      3      |    11200    |    40        |  296.25


 2. Liczba NOUTi przy zalozeniu B = 0,002
 Brama(GWi)  |  AOUT_PSTN_d  |  NOUTi
      1      |    287.00     |   326
      2      |    293.00     |   332
      3      |    296.25     |   335


 3. Przeliczone NOUTi na NPCM_30/32 dla kierunku OUT
 Brama(GWi)  |  NPCM_30/32_OUTi
      1      |    11
      2      |    12
      3      |    12


 4. Rzeczywista liczba NR_OUTi laczy w kierunku OUT
 Brama(GWi)  |  NR_OUTi
      1      |    330
      2      |    360
      3      |    360


 5. Liczba DSP dla kierunku OUT przy zalozeniu n_dsp = 4
 Brama(GWi)  |  NDSP_OUTi
      1      |    83
      2      |    90
      3      |    90

-------------------------------------------Obliczenie NOUTi raz NDSP_OUTi----------------------------------
1. Obliczenia dla przypadku gdzie siec docelowa jest siecia PSTN/ISDN/GSM

 lambda_G711 = 100 [1/s]
 lb_pakiet = 80 [B]
 lBnag = 40 [B]
 bcal_G711 = 960 [b]

 cRT_11 = 0.00 [Mb/s]
 cRT_12 = 1.90 [Mb/s]
 cRT_13 = 6.74 [Mb/s]
 cRT_21 = 7.49 [Mb/s]
 cRT_22 = 0.00 [Mb/s]
 cRT_23 = 11.23 [Mb/s]
 cRT_31 = 9.22 [Mb/s]
 cRT_32 = 9.22 [Mb/s]
 cRT_33 = 0.00 [Mb/s]

 2. Obliczenia dla przypadku gdzie siec docelowa jest siecia IP

 lambda_G729 = 40 [1/s]
 lb_pakiet = 25 [B]
 lBnag = 40 [B]
 bcal_G729 = 520 [b]

 cRT_14 = 0.56 [Mb/s]
 cRT_15 = 0.56 [Mb/s]
 cRT_16 = 0.75 [Mb/s]
 cRT_24 = 1.62 [Mb/s]
 cRT_25 = 1.62 [Mb/s]
 cRT_26 = 0.81 [Mb/s]
 cRT_34 = 1.00 [Mb/s]
 cRT_35 = 2.50 [Mb/s]
 cRT_36 = 2.50 [Mb/s]

 3. Obliczenia dla przypadku gdzie siec zrodlowa jest siec PSTN/ISDN/GSM

 lambda_G711 = 100 [1/s]
 lb_pakiet = 80 [B]
 lBnag = 40 [B]
 bcal_G711 = 960 [b]

 cRT_11 = 0.00 [Mb/s]
 cRT_21 = 7.49 [Mb/s]
 cRT_31 = 9.22 [Mb/s]
 cRT_12 = 1.90 [Mb/s]
 cRT_22 = 0.00 [Mb/s]
 cRT_32 = 9.22 [Mb/s]
 cRT_13 = 6.74 [Mb/s]
 cRT_23 = 11.23 [Mb/s]
 cRT_33 = 0.00 [Mb/s]

 4. Obliczenia dla przypadku gdzie siec zrodlowa jest siec IP

 lambda_G729 = 40 [1/s]
 lb_pakiet = 25 [B]
 lBnag = 40 [B]
 bcal_G729 = 520 [b]

 cRT_41 = 0.57 [Mb/s]
 cRT_51 = 1.15 [Mb/s]
 cRT_61 = 1.16 [Mb/s]
 cRT_42 = 1.14 [Mb/s]
 cRT_52 = 1.15 [Mb/s]
 cRT_62 = 1.75 [Mb/s]
 cRT_43 = 1.14 [Mb/s]
 cRT_53 = 0.58 [Mb/s]
 cRT_63 = 1.16 [Mb/s]

 Obliczenia dla przypadku gdzie siec zrodlowa jest siec IP i docelowa jest IP

 cRT_44 = 0.00 [Mb/s]
 cRT_45 = 1.72 [Mb/s]
 cRT_46 = 1.14 [Mb/s]
 cRT_54 = 2.31 [Mb/s]
 cRT_55 = 0.00 [Mb/s]
 cRT_56 = 0.58 [Mb/s]
 cRT_64 = 1.16 [Mb/s]
 cRT_65 = 0.58 [Mb/s]
 cRT_66 = 0.00 [Mb/s]

 OSTATECZNIE
 Brama(GWi)  |  Kierunek A Ab. - B Ab.  |  Kierunek B Ab. - A Ab  |  cRT_IP
      1      |    10.51                 |        19.60            |   30.11
      2      |    22.78                 |        15.16            |   37.94
      3      |    24.42                 |        20.86            |   45.28

 ---------Obliczenia parametrów jakościowych na ustalonej drodze - IPLR, IPDT, IPDVmax-----------


<-- Klasa RT  -->
 Lacze    | Przepustowosc strumienia klasy RT[Mb/s]   | Natezenie ruchu A_RT klasy RT oferowane na lacze
 RB1_RR1  |           30.11                    |               0.201
 RR1_RR2  |           27.26                     |               0.182
 RR2_RR5  |           22.03                    |               0.000
 RR5_RB3  |           28.66                     |               0.000

 Lacze    |   IPLR         |    E(Tocz)[s]  |    E(Tnad)[s]  |  E(prop)[s]    |  IPDT[s]
 RB1_RR1  |  2.603909e-04  |  1.596310e-06  |  6.400000e-06  |  1.500000e-04  |  1.579963e-04
 RR1_RB2  |  1.621440e-04  |  1.414687e-06  |  6.400000e-06  |  3.500000e-04  |  3.578147e-04
 RR2_RB5  |  0  |  0  |  6.400000e-06  |  2.500000e-04  |  2.564000e-04
 RR5_RB3  |  0  |  0  |  6.400000e-06  |  3.500000e-04  |  3.564000e-04

 Lacze    |     IPDTmax[s]     |    IPDTmin[s]  |   IPDV
 RB1_RR1  |    1.141333e-04    |  6.400000e-06  |  1.077333e-04
 RR1_RR2  |    1.141333e-04    |  6.400000e-06  |  1.077333e-04
 RR2_RR5  |    1.141333e-04    |  6.400000e-06  |  1.077333e-04
 RR5_RB3  |    1.141333e-04    |  6.400000e-06  |  1.077333e-04

 Parametry jakosciowe klasy RT na zadanej drodze;

 IPLR = 4.225348e-04
 IPDT = 1.128611e-03
 IPDVmax = 4.309333e-04


 Przeplywsnosci strumieni NRT
 cNRT_44 = 0.00 [Mb/s]
 cNRT_45 = 67.76 [Mb/s]
 cNRT_46 = 67.76 [Mb/s]
 cNRT_54 = 70.84 [Mb/s]
 cNRT_55 = 0.00 [Mb/s]
 cNRT_56 = 70.84 [Mb/s]
 cNRT_64 = 110.39 [Mb/s]
 cNRT_65 = 47.31 [Mb/s]
 cNRT_66 = 0.00 [Mb/s]

<-- Klasa NRT  -->
 Lacze    | Przepustowosc strumienia klasy RT[Mb/s]  | Sumaryczna przeplywnosc dla strumienia klasy NRT[Mb/s]
 RB5_RR1  |           5.74                           |               67.76
 RR1_RB4  |           5.16                           |               67.76

 Lacze    | Natezenie ruchu oferowanego A_RT na lacze  | Natezenie ruchu oferowanego A_NRT na lacze
 RB5_RR1  |           0.04                             |               0.470
 RR1_RB4  |           0.03                             |               0.468

 Lacze    |     IPLR       |    IPDT[s]      |   IPDTmax[s]   |   IPDTmin[s]   |  IPDVmax[s]
 RB5_RR1  |  1.556289e-09  |  4.048838e-04   |  2.220448e-03  |  8.213333e-05  |  2.138314e-03
 RR1_RB4  |  1.406500e-09  |  3.543341e-04   |  2.211519e-03  |  8.213333e-05  |  2.129386e-03

 Parametry jakosciowe klasy NRT na zadanej drodze;
 IPLR = 2.962788e-09
 IPDT = 7.592180e-04
 IPDVmax = 4.267700e-03

>>


%------------------------------------DANE---------------------------------%
A_BramaGWi = [150 330 410];   %Natezenie w Erlangach
RT_IP = [11000 11100 11200];    %Intensywność pakietów 1/s, Klasy ruchu 'Real Time',
NRT_IP = [11000 11500 12800];   %Intensywność pakietów 1/s, Klasy ruchu 'No Real Time',

%IP-VBR1(RT - mowa), kodek G.729, RTP/UDP/IP/SDH
lRT = 20;   %średnia długość pakietu [B]
Tpak_G729 = 25; %[ms]
Tpak_G729_s = 0.025;    %[s]

%IP-VBR2(NRT), TCP/IP/SDH
lNRT = 1500;    %średnia długość pakietu [B]

%PSTN/ISDN/GSM na VBR1(RT), kodek G.711, RTP/UDP/IP/SDH
Tpak_G711 = 10; %[ms]
Tpak_G711_s = 0.01;  %[s]

WZ_RT = [0 0.11 0.39 0.15 0.15 0.2;
        0.2 0 0.3 0.2 0.2 0.1;
        0.2 0.2 0 0.1 0.25 0.25;
        0.1 0.2 0.2 0 0.3 0.2;
        0.2 0.2 0.1 0.4 0 0.1;
        0.2 0.3 0.2 0.2 0.1 0]; %macierz wspczynnikw zainteresowań dla klasy RT

WZ_NRT = [0 0 0 0 0 0;
          0 0 0 0 0 0;
          0 0 0 0 0 0;
          0 0 0 0 0.5 0.5;
          0 0 0 0.5 0 0.5;
          0 0 0 0.7 0.3 0]; %macierz wspczynnikw zainteresowań dla klasy NRT

K1 = 4; %pakietów
K2 = 25;  %pakietów

IPNAG_IP = 20;  %[B],   długość nagłówka IP
UDPNAG_UDP = 8; %[B],   długość nagłówka UDP
RTPNAG_RTP = 12;    %[B],   długość nagłówka RTP

%
%Obliczenia czas zacząć!!!!!!!!!!!!!!!!!!!!
%
%---------------------------Obliczenie NINi oraz NDSP_INi-----------------%
%1
%Dowiedzieć się czy prawdopodobieństwo straty
%jest równe 0.002
%Obliczanie N na podstawie E1,N
NINi = [0 0 0]; %Liczba łączy wejściowych
fprintf('-------------------------------------------WYNIKI--------------------------------------------------------  \n');
fprintf('-------------------------------------------Obliczenie NINi oraz NDSP_INi----------------------------------  \n');
fprintf('1. NINi przy zalozeniu dopuszczalnego prawdopodobienstwa straty B=0.002 \n Brama(GWi)  |  Ai[Erl]  |  NINi  \n');
for i = 1:1:3
    A = A_BramaGWi(i);
    E1 = 1;
    E1N1 = 1;
    while E1 > 0.002
        NINi(i) = NINi(i)+1;
        E1 = (A*E1N1)/(NINi(i) + A*E1N1);
        E1N1 = E1;
    end
    fprintf('      %d      |    %d    |  %d  \n', i,A_BramaGWi(i), NINi(i) );
end

%2
%Przeliczanie NIni na NPCM, dla kierunku IN
NPCM_30_32_INi = [0 0 0];   %Liczba łączy wejściowych w PCM, zaookrąglona w górę
fprintf('\n \n 2. Przeliczone NINi na NPCM_30/32 dla kierunku IN \n Brama(GWi)  |  NPCM_30/32_INi  \n');
for i = 1:1:3
    NPCM_30_32_INi(i) = ceil(NINi(i)/30);
    fprintf('      %d      |    %d      \n', i,NPCM_30_32_INi(i) );
end

%3
%Rzeczywista liczba NR_INi, dla kierunku IN
NR_INi = [0 0 0];   %Rzeczywista licza łączy wejściowych
fprintf('\n \n 3. Rzeczywista liczba NR_INi laczy w kierunku In \n Brama(GWi)  |  NR_INi]  \n');
for i = 1:1:3
    NR_INi(i) = NPCM_30_32_INi(i) * 30;
    fprintf('      %d      |    %d      \n', i,NR_INi(i) );
end

%4
%Obliczenie DSP dla kierunku IN, przy założeniu że ndsp = 4, ale rozumiem,
%że można wybrać 4,8 lub 16 ??????
nDSP = 4;   %Liczba strumieni obłusgiwanych przez jeden procesor sygnałowy
NDSP_INi = [0 0 0]; %Ilość procesorów sygnałowych dla kierunku IN
fprintf('\n \n 4. Liczba DSP dla kierunku IN przy zalozeniu n_dsp = 4 \n Brama(GWi)  |  NDSP_INi  \n');
for i = 1:1:3
    NDSP_INi(i) = ceil(NR_INi(i) / 4);
    fprintf('      %d      |    %d      \n', i,NDSP_INi(i) );
end


%
%------------------------Obliczenia NOUTi oraz NDSP_OUTi------------------%
%
%1
%Obliczenie liczby wyposażeń na wyjściu bramy GW do PSTN/ISDN/GSM
AOUT_PSTN_d = [0 0 0];
lambda_G729 = 1/Tpak_G729_s;
Ai_k = [0 0 0 0 0 0];

%Uzupełnienie Ai_k, pierwsze 3 to PSTN, a pozostałe 3 z IP
fprintf('\n-------------------------------------------Obliczenie NOUTi raz NDSP_OUTi----------------------------------  \n');
fprintf(' \n 1. Obliczenia Aout,PSTN_d \n Brama(GWi)  |  lambda_RT  |  lambda_G729  |  AOUT_PSTN_d    \n');
for i = 1:1:3
    Ai_k(i) = A_BramaGWi(i);
end
for i = 1:1:3
    Ai_k(i+3) = RT_IP(i)/lambda_G729;
end

%Oblicznie Aout
for i = 1:1:3   %Ilość PSTN
    for k = 1:1:6
    AOUT_PSTN_d(i) = AOUT_PSTN_d(i) + (WZ_RT(k,i)*Ai_k(k));
    end
     fprintf('      %d      |    %d    |    %d        |  %0.2f \n', i,RT_IP(i),lambda_G729, AOUT_PSTN_d(i) );
end

%2
%Obliczenia NOUTi dla AOUT
NOUTi=[0 0 0];  %%Liczba łączy wyjściowych
fprintf('\n \n 2. Liczba NOUTi przy zalozeniu B = 0,002 \n Brama(GWi)  |  AOUT_PSTN_d  |  NOUTi  \n');
for i = 1:1:3
    A = AOUT_PSTN_d(i);
    E1 = 1;
    E1N1 = 1;
    while E1 > 0.002
        NOUTi(i) = NOUTi(i)+1;
        E1 = (A*E1N1)/(NOUTi(i) + A*E1N1);
        E1N1 = E1;
    end
     fprintf('      %d      |    %0.2f     |   %d  \n', i,AOUT_PSTN_d(i), NOUTi(i) );
end

%3
%Przeliczonie NOUTi z Gw na PCM30/32
NPCM_30_32_OUTi = [0 0 0];   %Liczba łączy wyjściowych w PCM, zaookrąglona w górę
fprintf('\n \n 3. Przeliczone NOUTi na NPCM_30/32 dla kierunku OUT \n Brama(GWi)  |  NPCM_30/32_OUTi  \n');

for i = 1:1:3
    NPCM_30_32_OUTi(i) = ceil(NOUTi(i)/30);
    fprintf('      %d      |    %d      \n', i,NPCM_30_32_OUTi(i) );
end

%4
%Rzeczywista liczba NR łączy na kierunku OUT
NR_OUTi = [0 0 0];   %Rzeczywista licza łączy wejściowych
fprintf('\n \n 4. Rzeczywista liczba NR_OUTi laczy w kierunku OUT \n Brama(GWi)  |  NR_OUTi  \n');
for i = 1:1:3
    NR_OUTi(i) = NPCM_30_32_OUTi(i) * 30;
    fprintf('      %d      |    %d      \n', i,NR_OUTi(i) );
end

%5
%Obliczenie DSP dla kierunku OUT, przy założeniu że ndsp = 4
NDSP_OUTi = [0 0 0]; %Ilość procesorów sygnałowych dla kierunku IN
fprintf('\n \n 5. Liczba DSP dla kierunku OUT przy zalozeniu n_dsp = 4 \n Brama(GWi)  |  NDSP_OUTi  \n');
for i = 1:1:3
    NDSP_OUTi(i) = ceil(NR_OUTi(i) / 4);
    fprintf('      %d      |    %d      \n', i,NDSP_OUTi(i) );
end


%
%------------Obliczenia CRT_IP_i -----------------------------------------%
%
lB_pakiet = Tpak_G711_s/(125*10^(-6));  %[B], w mianowniku 125us,
lBnag = IPNAG_IP + UDPNAG_UDP + RTPNAG_RTP; %[B]
bcal_G711 = 8*(lB_pakiet + lBnag);  %[b]
lambda_G711 = 1/Tpak_G711_s;  %[1/s]

fprintf('\n-------------------------------------------Obliczenie NOUTi raz NDSP_OUTi----------------------------------  \n');

%1
%Obliczamy cRT_PSTN,i,j
fprintf('1. Obliczenia dla przypadku gdzie siec docelowa jest siecia PSTN/ISDN/GSM   \n');
fprintf('\n lambda_G711 = %d [1/s] \n lb_pakiet = %d [B] \n lBnag = %d [B] \n bcal_G711 = %d [b]   \n\n', lambda_G711, lB_pakiet, lBnag, bcal_G711  );

cRT_PSTN_z_d = [0 0 0;
                0 0 0;
                0 0 0];


for i = 1:1:3   %Sieć PSTN wejściowa
    for k = 1:1:3   %Sieć PSTN docelowa
    cRT_PSTN_z_d(i,k) = (WZ_RT(i,k)*NR_INi(i)*lambda_G711*bcal_G711)/10^(6);
    fprintf(' cRT_%d%d = %0.2f [Mb/s] \n', i, k, cRT_PSTN_z_d(i,k) );
    end
end

%2
%Obliczamy cRT_PSTN_i,j
lambda_G729 = 1/Tpak_G729_s;  %[s]
lB_pakiet_IP = 25;  %[B],
bcal_G729 = 8*(lB_pakiet_IP + lBnag);  %[b]

fprintf('\n 2. Obliczenia dla przypadku gdzie siec docelowa jest siecia IP   \n');
fprintf('\n lambda_G729 = %d [1/s] \n lb_pakiet = %d [B] \n lBnag = %d [B] \n bcal_G729 = %d [b]   \n\n', lambda_G729, lB_pakiet_IP, lBnag, bcal_G729  );
cRT_PSTN_z_IP_d = [0 0 0;
                   0 0 0;
                   0 0 0];


for i = 1:1:3   %Sieć PSTN docelowa
    for k = 4:1:6   %Sieć IP źródłowa
    cRT_PSTN_z_IP_d(i,k-3) = (WZ_RT(i,k)*NR_INi(i)*lambda_G729*bcal_G729)/10^(6); %'-3' żeby macierz byłą 3x3 a nie 3x6
    fprintf(' cRT_%d%d = %0.2f [Mb/s] \n', i, k, cRT_PSTN_z_IP_d(i,k-3) );
    end
end

%3
%Obliczamy cRT_jmPSTN_i
fprintf('\n 3. Obliczenia dla przypadku gdzie siec zrodlowa jest siec PSTN/ISDN/GSM   \n');
fprintf('\n lambda_G711 = %d [1/s] \n lb_pakiet = %d [B] \n lBnag = %d [B] \n bcal_G711 = %d [b]   \n\n', lambda_G711, lB_pakiet, lBnag, bcal_G711  );
cRT_PSTN_d_z = [0 0 0;
                0 0 0;
                0 0 0];
for i = 1:1:3   %Sieć PSTN docelowa
    for k = 1:1:3   %Sieć PSTN źródłowa
    cRT_PSTN_d_z(k,i) = (WZ_RT(k,i)*NR_INi(k)*lambda_G711*bcal_G711)/10^(6);
    fprintf(' cRT_%d%d = %0.2f [Mb/s] \n', k, i, cRT_PSTN_d_z(k,i) );
    end
end

%4
%Obliczamy cRTj,PSTNi
fprintf('\n 4. Obliczenia dla przypadku gdzie siec zrodlowa jest siec IP   \n');
fprintf('\n lambda_G729 = %d [1/s] \n lb_pakiet = %d [B] \n lBnag = %d [B] \n bcal_G729 = %d [b]   \n\n', lambda_G729, lB_pakiet_IP, lBnag, bcal_G729  );
cRT_PSTN_d_IP_z = [0 0 0;
                   0 0 0;
                   0 0 0];
for i = 1:1:3   %Sieć PSTN docelowa
   for k = 4:1:6   %Sieć IP źródłowa
    cRT_PSTN_d_IP_z(k-3,i) = (WZ_RT(k,i)*RT_IP(k-3)*bcal_G729)/10^(6); %'-3' żeby macierz byłą 3x3 a nie 3x6
    fprintf(' cRT_%d%d = %0.2f [Mb/s] \n', k, i, cRT_PSTN_d_IP_z(k-3,i) );
    end
end

cRT = [cRT_PSTN_z_d cRT_PSTN_z_IP_d; cRT_PSTN_d_IP_z zeros(3,3)];   %połączone macierze policzone w punktach od 1 do 4


fprintf('\n Obliczenia dla przypadku gdzie siec zrodlowa jest siec IP i docelowa jest IP  \n\n');
for i = 4:1:6   %Sieć IP docelowa
   for k = 4:1:6   %Sieć IP źródłowa
    cRT(i,k) = (WZ_RT(i,k)*RT_IP(i-3)*bcal_G729)/10^(6); %'-3' żeby macierz byłą 3x3 a nie 3x6
    fprintf(' cRT_%d%d = %0.2f [Mb/s] \n', i, k, cRT(i,k) );
    end
end

fprintf(' \n OSTATECZNIE \n Brama(GWi)  |  Kierunek A Ab. - B Ab.  |  Kierunek B Ab. - A Ab  |  cRT_IP    \n');
KierunekA_B = [0 0 0];  %Macierz dla kierunku A -> B
KierunekB_A = [0 0 0];  %MAcierz dla kierunku B -> A
cRT_IP = [0 0 0];   %suma przepływności

for i = 1:1:3
    for j = 1:1:6
        KierunekA_B(i) = KierunekA_B(i) + cRT(i,j);
        KierunekB_A(i) = KierunekB_A(i) + cRT(j,i);
    end
    cRT_IP(i) = KierunekA_B(i) + KierunekB_A(i);
    fprintf('      %d      |    %0.2f                 |        %0.2f            |   %0.2f \n', i,KierunekA_B(i), KierunekB_A(i), cRT_IP(i) );
end


%Zadanie 1.2
%Obliczanie parametrów jakościowych na ustalonej drodze -IPLR,IPDT,IPDVmax
%
%-------------------------------------------------Zadanie1.2---------------------------------------------%
%
fprintf('\n ---------Obliczenia parametrów jakościowych na ustalonej drodze - IPLR, IPDT, IPDVmax-----------  \n\n');
fprintf('\n<-- Klasa RT  -->\n');
%dpRB1,RB4 = {RB1, RR1, RR2, RR5, RB3} ---- RT -- /// 1 - RB1 - RR1 /// 2 - RR1 - RR2 /// 3 - RR2 - RR5 /// 4 - RR5 - RB3
%dpRB4,RB6 = {RB5, RR1, RB4} --------NRT

C = 150; %[Mb/s]
V = 200000;  %[Km/s]
RB1_RB3_S = [30 70 50 70];    %[km]
RB5_RB4_S = [50 40];
bcal_RT = 8*(lRT + lBnag);
bcal_NRT = 8*(lNRT + lBnag);    %bcal dla drogi NRT

RB1_RB3 = [0 0 0 0];

RB1_RB3(1) = cRT(1,2) + cRT(1,3) + cRT(1,4) + cRT(1,5) + cRT(1,6) + cRT(2,1) + cRT(3,1) + cRT(4,1) + cRT(5,1) + cRT(6,1);
RB1_RB3(2) = cRT(1,2) + cRT(1,3) + cRT(1,6) + cRT(2,1) + cRT(3,1) + cRT(6,1);
RB1_RB3(3) = cRT(1,3) + cRT(2,5) + cRT(4,3) + cRT(6,5) + cRT(3,1) + cRT(3,4) + cRT(5,2) + cRT(5,6);
RB1_RB3(4) = cRT(1,3) + cRT(4,3) + cRT(5,3) + cRT(2,1) + cRT(3,1) + cRT(3,4) + cRT(3,5);

RB1_RB3_ART = [0 0 0 0];

RB1_RB3_ART(1) = RB1_RB3(1)/C;
RB1_RB3_ART(2) = RB1_RB3(2)/C;

fprintf(' Lacze    | Przepustowosc strumienia klasy RT[Mb/s]   | Natezenie ruchu A_RT klasy RT oferowane na lacze  \n');
fprintf(' RB1_RR1  |           %0.2f                    |               %0.3f  \n', RB1_RB3(1),RB1_RB3_ART(1));
fprintf(' RR1_RR2  |           %0.2f                     |               %0.3f  \n', RB1_RB3(2), RB1_RB3_ART(2));
%dodane przez nas sciezki
fprintf(' RR2_RR5  |           %0.2f                    |               %0.3f  \n', RB1_RB3(3),RB1_RB3_ART(3));
fprintf(' RR5_RB3  |           %0.2f                     |               %0.3f  \n', RB1_RB3(4), RB1_RB3_ART(4));

RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT = [0 0 0 0 0;  %<----------------------- suma we wzorze
                                        0 0 0 0 0;
                                        0 0 0 0 0;
                                        0 0 0 0 0];
fprintf('\n Lacze    |   IPLR         |    E(Tocz)[s]  |    E(Tnad)[s]  |  E(prop)[s]    |  IPDT[s]   \n');
for i = 1:1:4   %Obliczanie elemenów tabeli z IPLR_ETocz_ETnad_Etprop_IPDT

    A = RB1_RB3_ART(i);
    RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,1) = ((1 - A)/(1 - (A^(K1+2))))*(A^(K1+1));
    RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,4) = RB1_RB3_S(i)/V;
    RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,3) = bcal_G711/(C*10^6);

    mi = 1/RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,3);
    X = (A/mi);
    Y = (1+(A^K1)*(K1*A - (K1+1)));
    Z = (1 - A)*(1 - (A^(K1+2)));

    RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,2) = (X*Y)/Z;

    RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,5) = RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,2) + RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,3) + RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,4);

end
fprintf(' RB1_RR1  |  %d  |  %d  |  %d  |  %d  |  %d    \n', RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(1,1),RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(1,2), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(1,3), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(1,4), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(1,5));
fprintf(' RR1_RB2  |  %d  |  %d  |  %d  |  %d  |  %d    \n', RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(2,1),RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(2,2), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(2,3), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(2,4), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(2,5));
fprintf(' RR2_RB5  |  %d  |  %d  |  %d  |  %d  |  %d    \n', RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(3,1),RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(3,2), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(3,3), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(3,4), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(3,5));
fprintf(' RR5_RB3  |  %d  |  %d  |  %d  |  %d  |  %d    \n', RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(4,1),RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(4,2), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(4,3), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(4,4), RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(4,5));


RB1_RB3_IPDTmax_IPDTmin_IPDVmax = [0 0 0;
                                   0 0 0;
                                   0 0 0;
                                   0 0 0];


fprintf('\n Lacze    |     IPDTmax[s]     |    IPDTmin[s]  |   IPDV   \n');
for i = 1:1:4

    Et_NRT = bcal_NRT/(C*10^6);

    RB1_RB3_IPDTmax_IPDTmin_IPDVmax(i,2) = RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,3);
    RB1_RB3_IPDTmax_IPDTmin_IPDVmax(i,1) = (K1 + 1)*RB1_RB3_IPDTmax_IPDTmin_IPDVmax(i,2) + Et_NRT;
    RB1_RB3_IPDTmax_IPDTmin_IPDVmax(i,3) = RB1_RB3_IPDTmax_IPDTmin_IPDVmax(i,1) - RB1_RB3_IPDTmax_IPDTmin_IPDVmax(i,2);

end
fprintf(' RB1_RR1  |    %d    |  %d  |  %d      \n', RB1_RB3_IPDTmax_IPDTmin_IPDVmax(1,1), RB1_RB3_IPDTmax_IPDTmin_IPDVmax(1,2), RB1_RB3_IPDTmax_IPDTmin_IPDVmax(1,3));
fprintf(' RR1_RR2  |    %d    |  %d  |  %d      \n', RB1_RB3_IPDTmax_IPDTmin_IPDVmax(2,1), RB1_RB3_IPDTmax_IPDTmin_IPDVmax(2,2), RB1_RB3_IPDTmax_IPDTmin_IPDVmax(2,3));
fprintf(' RR2_RR5  |    %d    |  %d  |  %d      \n', RB1_RB3_IPDTmax_IPDTmin_IPDVmax(3,1), RB1_RB3_IPDTmax_IPDTmin_IPDVmax(3,2), RB1_RB3_IPDTmax_IPDTmin_IPDVmax(3,3));
fprintf(' RR5_RB3  |    %d    |  %d  |  %d      \n', RB1_RB3_IPDTmax_IPDTmin_IPDVmax(4,1), RB1_RB3_IPDTmax_IPDTmin_IPDVmax(4,2), RB1_RB3_IPDTmax_IPDTmin_IPDVmax(4,3));


%Parametry jakościowe klasy RT !!!!!!!
fprintf('\n Parametry jakosciowe klasy RT na zadanej drodze;  \n');

IPLR = 0;
IPDT = 0;
IPDV_MAX = 0;

%---------------------------------------

for i = 1:1:4
    IPLR = IPLR + RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,1);
    IPDT = IPDT + RB1_RB3_IPLR_ETocz_ETnad_Etprop_IPDT(i,5);
    IPDV_MAX = IPDV_MAX + RB1_RB3_IPDTmax_IPDTmin_IPDVmax(i,3);
end
fprintf('\n IPLR = %d \n IPDT = %d \n IPDVmax = %d \n   \n', IPLR, IPDT, IPDV_MAX);


fprintf('\n Przeplywsnosci strumieni NRT  \n');
cNRT = [0 0 0 0 0 0;
        0 0 0 0 0 0;
        0 0 0 0 0 0;
        0 0 0 0 0 0;
        0 0 0 0 0 0;
        0 0 0 0 0 0];

for i = 4:1:6
    for j = 4:1:6
        cNRT (i,j) = (WZ_NRT(i,j)*NRT_IP(i-3)*bcal_NRT)/10^(6);
        fprintf(' cNRT_%d%d = %0.2f [Mb/s] \n', i, j, cNRT (i,j) );
    end
end

%--------------------------------------------------NRT-------------------------------
fprintf('\n<-- Klasa NRT  -->\n');
fprintf(' Lacze    | Przepustowosc strumienia klasy RT[Mb/s]  | Sumaryczna przeplywnosc dla strumienia klasy NRT[Mb/s]  \n');
RB5_RB4 = [0 0; %Sumaryczne przepływności dla sktrumienia klasy RT i NRT
           0 0];

RB5_RB4(1,1) = cRT(5,1) + cRT(5,4) + cRT(1,5) + cRT(4,5);
RB5_RB4(1,2) = cNRT(4,6);
RB5_RB4(2,1) = cRT(1,4) + cRT(5,4) + cRT(4,1) + cRT(4,5);
RB5_RB4(2,2) = cNRT(4,6);

fprintf(' RB5_RR1  |           %0.2f                           |               %0.2f  \n', RB5_RB4(1,1),RB5_RB4(1,2));
fprintf(' RR1_RB4  |           %0.2f                           |               %0.2f  \n', RB5_RB4(2,1),RB5_RB4(2,2));

fprintf('\n Lacze    | Natezenie ruchu oferowanego A_RT na lacze  | Natezenie ruchu oferowanego A_NRT na lacze  \n');
RB5_RB4_ART_ANRT = [0 0;    %Natezenie ruchu A_RT i A_NRT na łącze
                    0 0];

for i = 1:1:2
    RB5_RB4_ART_ANRT(i,1) = RB5_RB4(i,1)/C;
    RB5_RB4_ART_ANRT(i,2) = RB5_RB4(i,2)/(C-RB5_RB4(i,1));
end
fprintf(' RB5_RR1  |           %0.2f                             |               %0.3f  \n', RB5_RB4_ART_ANRT(1,1),RB5_RB4_ART_ANRT(1,2));
fprintf(' RR1_RB4  |           %0.2f                             |               %0.3f  \n', RB5_RB4_ART_ANRT(2,1),RB5_RB4_ART_ANRT(2,2));

fprintf('\n Lacze    |     IPLR       |    IPDT[s]      |   IPDTmax[s]   |   IPDTmin[s]   |  IPDVmax[s]  \n');
RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax = [0 0 0 0 0;    %Parametry jakościowe klasy NRT
                                              0 0 0 0 0];
%Parametry jakościowe klasy NRT !!!!!!!
IPLR_NRT = 0;
IPDT_NRT = 0;
IPDV_MAX_NRT = 0;
%---------------------------------------
for i = 1:1:2
    A = RB5_RB4_ART_ANRT(i,2);
    A1 = RB5_RB4_ART_ANRT(i,1);
    RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(i,1) = ((1 - A)/(1 - (A^(K2+2))))*(A^(K2+1));

    ET_nad =  bcal_NRT/(C*10^6);
    mi = 1/ET_nad;
    X = (A/mi);
    Y = (1+(A^K2)*(K2*A - (K2+1)));
    Z = (1 - A)*(1 - (A^(K2+2)));

    Et_ocz = (X*Y)/Z;
    Et_prop = RB5_RB4_S(i)/V;

    RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(i,2) = ET_nad + Et_ocz + Et_prop;

    RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(i,4) = ET_nad;
    RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(i,3) = (K2+1)*ET_nad/(1-A1);

    RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(i,5) = RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(i,3) - RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(i,4);

    % Sumowanie wartości
    IPLR_NRT = IPLR_NRT + RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(i,1);
    IPDT_NRT = IPDT_NRT + RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(i,2);
    IPDV_MAX_NRT = IPDV_MAX_NRT + RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(i,5);

end
fprintf(' RB5_RR1  |  %d  |  %d   |  %d  |  %d  |  %d  \n', RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(1,1),RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(1,2), RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(1,3), RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(1,4),RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(1,5));
fprintf(' RR1_RB4  |  %d  |  %d   |  %d  |  %d  |  %d  \n', RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(2,1),RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(2,2), RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(2,3), RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(2,4),RB5_RB4_IPLR_IPDT_IPDTmax_IPDT_min_IPDVmax(2,5));
fprintf('\n Parametry jakosciowe klasy NRT na zadanej drodze;  \n');
fprintf(' IPLR = %d \n IPDT = %d \n IPDVmax = %d \n   \n', IPLR_NRT, IPDT_NRT, IPDV_MAX_NRT);