s171632_5

% 22.11.2020 r.
% 171632 Marta Trzaska
% Zad.5 Modelowanie odpowiedzi impulsowej kanału dla środowiska miejskiego

clc;
format compact;
%Szybkość próbkowania odpowiedzi impulsowej;
    N = 25000000;
    Ts = 1/N;

% 1. Generowanie tablicy odpowiedzi impulsowej kanału dla szybkości
% próbkowania N i ozdzielenie prążków przypadających między momentami próbkowania na
% dwa o określonych amplitudach

    opoznienie = [0 0.217 0.512 0.514 0.517 0.674 0.882 1.230 1.287 1.311 1.349 1.533 1.535 1.622 1.818 1.836 1.884 1.943 2.048 2.140];
    %opoznienie w mikrosekundach
    sr_ampl = [-5.7 -7.6 -10.1 -10.2 -10.2 -11.5 -13.4 -16.3 -16.9 -17.1 -17.4 -19 -19 -19.8 -21.5 -21.6 -22.1 -22.6 -23.5 -24.3];
    %sr_ampl w dB
    sr_ampl_lin = ones(1,20);
    odp_impul = ones(1,20);

    for j = 1:20
        sr_ampl_lin(j) = 10^(sr_ampl(j)/10);
    end

    i = 1;
    czas_max = opoznienie(20) * 10^(-6)

    for j= 0:Ts: czas_max + Ts
        moment_probk(i) = j;
        probki(i) = 0;
        i = i + 1;
    end

    probki(1) = sr_ampl_lin(1);
    for k = 2 : 20
        for l = 2 : (length(probki) - 1)
            if(moment_probk(l) >= ( opoznienie(k)*10^(-6) - Ts))
                roznica = opoznienie(k) - moment_probk(l);
                probki(l) = probki(l) + sr_ampl_lin(k) * (roznica/Ts);
                probki(l+1) = probki(l+1) + sr_ampl_lin(k) * ((Ts- roznica)/Ts);
                break;
            end
        end
    end

    rng default;
    x = randn(1, length(probki));
    y = randn(1, length(probki));

    for i = 1 : length(probki)
        rayleigh = sqrt( x(i)^2 + y(i)^2 );
        odp_impul(i) = probki(i) * rayleigh;
    end


% 3. %%%%%%%%%%%%%%%%%% WYKRESY %%%%%%%%%%%%%%%%%%

    f1=figure;
    stem(opoznienie, sr_ampl_lin); %sekwencja bitów wygenerowanych
    title("Względne śrendie amplitudy impulsów (skala liniowa)")
    xlabel("t [µs]");
    ylabel("Amplituda");
    axis([0 2.3 0 0.3]);
    %saveas(gcf, "sr_ampl_5.png")

    f2 = figure;
    stem(moment_probk, probki); %sekwencja bitów wygenerowanych
    title("Probki")
    xlabel("t [µs]");
    ylabel("Amplituda");
    %axis([0 2.3 0 0.3]);
    %saveas(gcf, "probki_5.png")

    f3 = figure;
    stem(moment_probk, odp_impul); %sekwencja bitów wygenerowanych
    title("Odpowiedź impulsowa kanału z uwzględnieniem rozkładu Rayleigha")
    xlabel("t [µs]");
    ylabel("Amplituda");
    %axis([0 2.3 0 0.3]);
    %saveas(gcf, "odp_impuls_5.png")