% Initialisation of variables:
D=2:1:8;
fminUplink=12.75*10^9;
fminDownlink=11.2*10^9;
c=3*10^8;
ntx=0.5;
nrx=0.6;
%a=0.091;
a=0.105;
b=[0.2 0.2 0.1]; %fix- step tracking - monopulse
ntrackingFixedUplink=[10.^(-1.2*(b(1)+(a.*D*fminUplink./(c*70))).^2)];
ntrackingStep=10^(-1.2*b(2)^2);
ntrackingMono=10^(-1.2*b(3)^2);
Tsky=10;%depending on frequency and elevation angle (here E=30° -> Tsky=10K)
Tground=50;
Tr=120;
Tf=290;
Loss=10^0.02; %Losses between antenna and receiver input.
Gfma=ntrackingFixedUplink.*ntx.*(pi*fminUplink/c*D).^2;
GfmaIndB=10*log10(Gfma);
Gstep=ntrackingStep*ntx*(pi*fminUplink/c*D).^2;
GstepIndB=10*log10(Gstep);
Gmono=ntrackingMono*ntx*(pi*fminUplink/c*D).^2;
GmonoIndB=10*log10(Gmono);

% Calcul of D=f(G):
Dfma=(c/(pi*fminUplink)).*(Gfma./(ntx.*ntrackingFixedUplink)).^(1/2);
Dstep=c/(pi*fminUplink).*(Gstep./(ntx*ntrackingStep)).^(1/2);
Dmono=c/(pi*fminUplink).*(Gmono./(ntx*ntrackingMono)).^(1/2);

% Calcul of cost:
Cfma=0.4.*Dfma.^2.6;
Cstep=0.4.*Dstep.^2.6+16.6.*Dstep.^0.17;
Cmono=0.4.*Dmono.^2.6+16.6.*Dmono.^0.17+40;
MatrixC=[Cfma;Cstep;Cmono]
MatrixGdb=[GfmaIndB;GstepIndB;GmonoIndB]

%Display the result:
figure(9)
subplot(2,1,1);
plot(MatrixGdb.',MatrixC.','-s','MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',5)
legend('fixed mounting','step tracking','monopulse',0);
text(MatrixGdb(1,1),MatrixC(1,1),'D=2m\rightarrow ','HorizontalAlignment','right')
text(MatrixGdb(2,1),MatrixC(2,1),'D=2m\rightarrow ','HorizontalAlignment','right')
text(MatrixGdb(3,1),MatrixC(3,1),'D=2m\rightarrow ','HorizontalAlignment','right')
text(MatrixGdb(1,7),MatrixC(1,7),'D=8m\rightarrow ','HorizontalAlignment','right')
text(MatrixGdb(2,7),MatrixC(2,7),'D=8m\rightarrow ','HorizontalAlignment','right')
text(MatrixGdb(3,7),MatrixC(3,7),'D=8m\rightarrow ','HorizontalAlignment','right')
title('Cost of the antenna in k$ depending on transmitting gain')
ylabel('Cost of the antenna in k$')
xlabel('Gtx in dB')

% Calcul of Tu at the receiver input:
Tu=(Tsky+Tground)/Loss+Tf*(1-1/Loss)+Tr;
% Calcul of ntrackingFixedDownlink
ntrackingFixedDownlink=[10.^(-1.2*(b(1)+(a.*D*fminDownlink./(c*70))).^2)];

% Calcul of Grx)ES
% Think to substract the loss between Antenna and LNA: L=0.2 dB
GrxFma=ntrackingFixedDownlink.*nrx.*(pi*fminDownlink/c*D).^2;
GrxFmaIndB=10*log10(GrxFma)
GrxStep=ntrackingStep*nrx*(pi*fminDownlink/c*D).^2;
GrxStepIndB=10*log10(GrxStep)
GrxMono=ntrackingMono*nrx*(pi*fminDownlink/c*D).^2;
GrxMonoIndB=10*log10(GrxMono)
% We need to take into account the feeder loss of 0.2 dB since Tu is calculated at the input of the LNA 
% G and T must be calculated at the same place
MatrixRxGdb=[GrxFmaIndB;GrxStepIndB;GrxMonoIndB]-0.2 ;% the loss between Antenna and LNA: L=0.2 dB
MatrixGoverTdb=MatrixRxGdb-10*log10(Tu)

%Display the result:
Matrix3RowsD=[D;D;D];
subplot(2,1,2);
plot(Matrix3RowsD.',MatrixGoverTdb.','-s','MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',5)
legend('fixed mounting','step tracking','monopulse',0);
title('Grx/T)ES in dB at the input of LNA depending on D')
ylabel('Grx/T)ES in dB at the input of LNA')
xlabel('D')