%% Pirmas
%Elektromagnetinės bangos dažnis f = 1, 02 GHz.
%Apskaičiuokite bangos ilgį ? [m].
f=1.02e09;
%sviesos greitis
c=2.998e08;
l=c/f;
disp(['Bangos ilgis ',num2str(l),' [m]']);

%% Antras
%Elektromagnetinės bangos ilgis ? = 15 cm. Apskaičiuokite jos dažnį f [MHz].
l=0.15;
%sviesos greitis
c=2.998e08;
f=c/l;
disp(['Daznis ',num2str(f),' [Hz]']);

%% Trecias
%Apskaičiuokite energijos srauto tankį WP atstumu d = 50 m nuo siųstuvo antenos, kai
%siųstuvo išėjime signalo stiprumas PS yra 10 W, siųstuvo antenos efektyvusis stiprinimas
%GS yra absoliučiais dydžiais 100. WP dydį išreikškite mW/cm2.

d=5000; %cm
P=10000; %mW
G=100;
W=(P*G)/(4*pi*d^2);
disp(['Energijos srauto tankis ',num2str(W),' [mW/cm2]']);

%% Ketvirtas

%L = 32,45 + 20*log10(d)+20*log10(f);

L = 138.6; %dB
f = 6350; %MHz
d = 10^((L-32.45-20*log10(f))/20);
disp(['Atstumas ',num2str(d),' [km]']);

%% Penktas

% a) Atstumas tarp siųstuvo ir imtuvo d = 27 km; elektromagnetinės bangos dažnis f = 15300 MHz;
d=27;
f=15300;
L1=32.45+ 20*log10(d) + 20*log10(f);
disp(['a) Nuostoliai L ',num2str(L1),' [dB]']);

% b) kaip pasikeis nuostoliai, jei elektromagnetinės bangos dažnį sumažinsime per pusę iki
% f = 7650 MHz, o atstumas d išliks nepakitęs?
f1=7650;
L2 = 32.45 + 20*log10(d) + 20*log10(f1);
L4=L2-L1;
disp(['b) Nuostoliai L ',num2str(L2),' [dB]']);
disp(['b) Pasikeite nuo pirmo ',num2str(L4),' [dB]']);


%c) kaip pasikeis nuostoliai, jei atstumą sumažinsime per pusę iki d = 13, 5 km, o dažnis
% išliks toks pat, kaip ir pirmuoju atveju, t.y. f = 15300 MHz?
d1=13.5;
L3 = 32.45 + 20*log10(d1)+20*log10(f);
L5=L3-L1;
disp(['c) Nuostoliai L ',num2str(L3),' [dB]']);
disp(['c) Pasikeite nuo pirmo ',num2str(L5),' [dB]']);


%% Sestas

d=27;
d1=13.5;
f=6500;
L = 32.45 + 20*log10(d) + 20*log10(f);
L1 = 32.45 + 20*log10(d1) + 20*log10(f);
L2=L1-L;
disp([' Nuostoliai L ',num2str(L),' [dB]']);
disp([' Sumazinus atstuma ',num2str(L1),' [dB]']);
disp([' Pasikeite ',num2str(L2),' [dB]']);

%% Septintas
d=2800;
d1=3400;
f=12.7;
L = 92.45 + 20*log10(d) + 20*log10(f);
L1 = 92.45 + 20*log10(d1) + 20*log10(f);
L2=L1-L;
disp([' Nuostoliai L ',num2str(L),' [dB]']);
disp([' Padidinus atstuma ',num2str(L1),' [dB]']);
disp([' Pasikeite ',num2str(L2),' [dB]']);

%% Astuntas
d1=10000;
d=20000;
f=15e09;
f1=20e09;
f2=10e09;
%sviesos greitis
c=2.998e08;
l=c/f;
l1=c/f1;
l2=c/f2;
disp(['Bangos ilgis 15GHz ',num2str(l),' [m]']);
disp(['Bangos ilgis 20GHz ',num2str(l1),' [m]']);
disp(['Bangos ilgis 10GHz ',num2str(l2),' [m]']);

R1 = sqrt(l*d1*(d-d1)/d);
disp(['Pirmas spindulys 15GHz ',num2str(R1),' [m]']);
R2 = sqrt(l1*d1*(d-d1)/d);
disp(['Antras spindulys 20GHz ',num2str(R2),' [m]']);
R3 = sqrt(l2*d1*(d-d1)/d);
disp(['Trecias spindulys 10GHz ',num2str(R3),' [m]']);