kcpsr5

clear all

Ts = 6 * 10^(-8);   %szybko�� pr�bkowania odpowiedzi impulsowej

opoznienia_amplituda_TU50 = [0, -5.7;
                            21.7e-8, -7.6;
                            51.2e-8, -10.1;
                            51.4e-8, -10.2;
                            51.7e-8, -10.2;
                            67.4e-8, -11.5;
                            88.2e-8, -13.4;
                            123.0e-8, -16.3;
                            128.7e-8, -16.9;
                            131.1e-8, -17.1;
                            134.9e-8, -17.4
                            153.3e-8, -19;
                            153.5e-8, -19;
                            162.2e-8, -19.8;
                            181.8e-8, -21.5;
                            183.6e-8, -21.6;
                            188.4e-8, -22.1;
                            194.3e-8, -22.6;
                            204.8e-8, -23.5;
                            214.0e-8, -24.3];
 j = 1;
 for i=0:Ts:opoznienia_amplituda_TU50(20,1)+Ts
    probkowanie(j, 1) = i;
    probkowanie(j, 2) = 0;
    j = j + 1;
 end

 for i=1:1:length(opoznienia_amplituda_TU50)
    opoznienia_amplituda_TU50(i,2) = 10^(opoznienia_amplituda_TU50(i,2)/20);
 end

figure
stem(opoznienia_amplituda_TU50(:,1),opoznienia_amplituda_TU50(:,2),'filled');
title('Przed przepr�bkowaniem');

 probkowanie(1,2) = opoznienia_amplituda_TU50(1,2);
 for n=2:1:length(opoznienia_amplituda_TU50)
    for m=2:1:length(probkowanie)-1
        if(probkowanie(m,1) >= (opoznienia_amplituda_TU50(n,1) - Ts))
            roznica = opoznienia_amplituda_TU50(n,1) - probkowanie(m,1);
            probkowanie(m,2) = probkowanie(m,2) + opoznienia_amplituda_TU50(n,2)*(roznica/Ts);
            probkowanie(m+1,2) = probkowanie(m+1,2) + opoznienia_amplituda_TU50(n,2)*(Ts-roznica)/Ts;
            break;
        end
    end
 end

figure
stem(probkowanie(:,1),probkowanie(:,2),'filled');
title('Po przepr�bkowaniu');

rayleigh = zeros(1,40);
rayleigh_przemnozenie = zeros(1,20);
normalny = randn(1,80);
 j = 1;
 for i=1:2:80
    rayleigh(j) = sqrt(normalny(i)^2 + normalny(i+1)^2);
    j = j+1;
 end

 for i=1:1:length(probkowanie)
    rayleigh_przemnozenie(i) = rayleigh(i) * probkowanie(i,2);
 end

figure
stem(probkowanie(:,1),rayleigh_przemnozenie,'filled');
title('Przemno�enie');

suma_przed = 0;
for i=1:1:length(opoznienia_amplituda_TU50)
    suma_przed = suma_przed + opoznienia_amplituda_TU50(i,2);
end
suma_przed

suma_po = 0;
for i=1:1:length(probkowanie)
    suma_po = suma_po + probkowanie(i,2);
end
suma_po