% Initialisation of variables:
IBOmin=zeros(1,31);
Aes=79.4*ones(1,31);
EIRPmin=zeros(1,31);
IBO=-25:0.01:0;
% Coeff Ad for the zone A,B,C in dB:
ad=[13.5;13.7;13.8];
% Create an array with the C/Ni corresponding to each SourceZone using our downlink frequency plan
% in order to minimise downlink interferences:
CoverNiTotal(1)=MinCoverNiTotalA(2);
CoverNiTotal(2)=MinCoverNiTotalB(2);
CoverNiTotal(3)=MinCoverNiTotalC(2);

%Result of the precedent graphs:
%the minimum C/N is obtained with ad(1) and CoverNiTotal(1). So when sourceZone and destinationZone are A.
%For that specific case:
MatrixIBO_5Rows=[IBO;IBO;IBO;IBO;IBO];
MatrixIBO_32Rows=[MatrixIBO_5Rows;MatrixIBO_5Rows;MatrixIBO_5Rows;MatrixIBO_5Rows;MatrixIBO_5Rows;MatrixIBO_5Rows;IBO;IBO];

for GoverTes=6.5:1:36.5
    CoverNu=10.^((35.7+IBO)./10);
    CoverNd=10.^((ad(1)+GoverTes+IBO+6-6.4.*exp((IBO+6)./6.4))./10);
    CoverNim=10.^((-2.*IBO+23-6.4.*exp((IBO+6)./6.4)-6.25.*exp((IBO+11.75)./6.25))./10);
    CoverNi=10^(CoverNiTotal(1)/10);
    CoverNtot=1./((1./CoverNu)+(1./CoverNd)+(1./CoverNim)+(1./CoverNi));
    CoverNtotdb=10.*log10(CoverNtot);
    if GoverTes==6.5 %to initialize MatrixCoverNtotdb in the first loop
        Matrix5CoverNtotdb=CoverNtotdb;
	else %concatenation of each row corresponding to each G/T
        Matrix5CoverNtotdb=[Matrix5CoverNtotdb;CoverNtotdb];
	end
end

%Find the minimum IBO requested to fullfit the C/N required without margin for each G/T)ES:
CoverNrequired=12.3;
for i=1:31
	j=1;
	while Matrix5CoverNtotdb(i,j)<CoverNrequired
        j=j+1;   
	end
    Matrix5CoverNtotdb(i,j);
    IBOmin(i)=IBO(j);
end

for GoverTes=9:1:39
    CoverNu=10.^((35.7+IBO)./10);
    CoverNd=10.^((ad(1)+GoverTes+IBO+6-6.4.*exp((IBO+6)./6.4))./10);
    CoverNim=10.^((-2.*IBO+23-6.4.*exp((IBO+6)./6.4)-6.25.*exp((IBO+11.75)./6.25))./10);
    CoverNi=10^(CoverNiTotal(1)/10);
    CoverNtot=1./((1./CoverNu)+(1./CoverNd)+(1./CoverNim)+(1./CoverNi));
    CoverNtotdb=10.*log10(CoverNtot);
    if GoverTes==9 %to initialize MatrixCoverNtotdb in the first loop
        Matrix6CoverNtotdb=CoverNtotdb;
	else %concatenation of each row corresponding to each G/T
        Matrix6CoverNtotdb=[Matrix6CoverNtotdb;CoverNtotdb];
	end
end

%Find the minimum IBO requested to fullfit the C/N required with 2dB margin for each G/T)ES:
CoverNrequiredWithMargin=14.3;
for i=1:31
	j=1;
	while Matrix6CoverNtotdb(i,j)<CoverNrequiredWithMargin
        j=j+1;   
	end
    Matrix6CoverNtotdb(i,j);
    IBOminWithMargin(i)=IBO(j);
end

% Calcul of EIRPmin without margin:
fprintf('minimum EIRP to have the C/N required without margin for each G/T')
VectorGoverTes=6.5:1:36.5
EIRPmin=IBOmin+Aes
% Calcul of EIRPmin with 2dB margin:
fprintf('minimum EIRP to have the C/N required with 2dB margin for each G/T')
VectorGoverTesWithMargin=9:1:39
EIRPminWithMargin=IBOminWithMargin+Aes

% Add constant row with value 12.3
ConstRow=12.3*ones(1,2501);
Matrix5CoverNtotdb=[Matrix5CoverNtotdb;ConstRow];
%Display the result:
figure(8)
subplot(2,2,1);
plot(MatrixIBO_32Rows.',Matrix5CoverNtotdb.')
legend('G/T=6.5dB','G/T=7.5dB','G/T=8.5dB','G/T=9.5dB','G/T=10.5dB','G/T=11.5dB','G/T=12.5dB','G/T=13.5dB','G/T=14.5dB','G/T=15.5dB','G/T=16.5dB',0);
text(0,12.3,'\leftarrowC/N required without margin=12.3dB')
title('C/N)total depending on G/T)EarthStation')
ylabel('C/N)total')
xlabel('IBO')

subplot(2,2,3);
plot(VectorGoverTes,EIRPmin)
title('minimum EIRP)ES versus G/T)ES to have the C/N required without margin')
ylabel('minimum EIRP)ES')
xlabel('G/T)ES')

subplot(2,2,4);
plot(VectorGoverTesWithMargin,EIRPminWithMargin)
title('minimum EIRP)ES versus G/T)ES to have the C/N required with 2dB margin')
ylabel('minimum EIRP)ES')
xlabel('G/T)ES')