Examples of link analysis from the reference

------------------------------------------------------------------------
Reference: Maral, G. and M. Bousquet. (1998.) Satellite Communications
           Systems, Third Edition. John Wiley.
------------------------------------------------------------------------
See also DBSignal, DBSignalToPower, DBWToPW, AntennaGain, Beamwidth3dB,
Gain3dB, LossAtmosphericGas, LossDepointing, LossFreeSpace,
LossPolarization, LossPrecipitation, TAntennaGround, TAttenuator,
TReceiver, TSky
------------------------------------------------------------------------

Contents

%-------------------------------------------------------------------------------
%	Copyright (c) 2001 Princeton Satellite Systems, Inc.
%   All rights reserved.
%-------------------------------------------------------------------------------

%-------------------------------------------------------------------------------

Example 1 pp. 23-24

%-------------------------------------------------------------------------------
fD            = 14; % GHz
pT            = DBSignal(100); % dBW
r             = 40000; % km
theta3DB      = 2;
effR          = 0.55;
effT          = 0.6;
dT            = 4;
fIllumination = 70;

gT            = AntennaGain(struct('type','circular aperture','area',pi*dT^2/4,'eff',effT),fD);
eIRPSL        = pT + gT;
lFS           = LossFreeSpace( fD, r );
gR            = Gain3dB( theta3DB, fIllumination, effR );
pR            = eIRPSL - lFS + gR;

fprintf('2.3.2 Example 1: The Uplink\n-----------------------------\n');
fprintf('Gain Transmit Antenna = %10.2f (dBW)\n',gT);
fprintf('EIRP SL               = %10.2f (dBW)\n',eIRPSL);
fprintf('Free Space Loss       = %10.2f (dB)\n',lFS);
fprintf('Gain Receive Antenna  = %10.2f (dB)\n',gR);
fprintf('Power Received        = %10.2f (dBW)\n',pR);
fprintf('Power Received        = %10.2f (pW)\n',DBWToPW(pR));


%-------------------------------------------------------------------------------
2.3.2 Example 1: The Uplink
-----------------------------
Gain Transmit Antenna =      53.15 (dBW)
EIRP SL               =      73.15 (dBW)
Free Space Loss       =     207.41 (dB)
Gain Receive Antenna  =      38.23 (dB)
Power Received        =     -96.03 (dBW)
Power Received        =     249.36 (pW)

Example 2 pp. 24-25

%-------------------------------------------------------------------------------
fD            = 12; % GHz
pT            = DBSignal(10); % dBW
r             = 40000; % km
theta3DB      = 2;
effT          = 0.55;
effR          = 0.6;
dR            = 4;
fIllumination = 70;

gT            = Gain3dB( theta3DB, fIllumination, effT );
eIRPSL        = pT + gT;
gR            = AntennaGain(struct('type','circular aperture','area',pi*dR^2/4,'eff',effR),fD);
lFS           = LossFreeSpace( fD, r );
pR            = eIRPSL - lFS + gR;

fprintf('\n\n2.3.2 Example 2: The Downlink\n-----------------------------\n');
fprintf('Gain Transmit Antenna = %10.2f (dBW)\n',gT);
fprintf('EIRP SL               = %10.2f (dBW)\n',eIRPSL);
fprintf('Free Space Loss       = %10.2f (dB)\n',lFS);
fprintf('Gain Receive Antenna  = %10.2f (dB)\n',gR);
fprintf('Power Received        = %10.2f (dBW)\n',pR);
fprintf('Power Received        = %10.2f (pW)\n',DBSignalToPower(pR));

%-------------------------------------------------------------------------------

2.3.2 Example 2: The Downlink
-----------------------------
Gain Transmit Antenna =      38.23 (dBW)
EIRP SL               =      48.23 (dBW)
Free Space Loss       =     206.07 (dB)
Gain Receive Antenna  =      51.81 (dB)
Power Received        =    -106.03 (dBW)
Power Received        =       0.00 (pW)

Example p. 32

%-------------------------------------------------------------------------------

tA   = 50;  % Antenna
tF   = 290; % Feed ambient
tERX = 50;  % Receiver
lFRX = 1;   % Feed loss in dB;
gFRX = 1/DBSignalToPower( 1 );

t1   = tA + TAttenuator( tF, lFRX) + tERX/gFRX;
t2   = t1*gFRX;

fprintf('\n\n2.4.6 Example\n-------------\n');
fprintf('Temperature at feeder input   = %10.2f (deg-K)\n',t1);
fprintf('Temperature at receiver input = %10.2f (deg-K)\n',t2);

%-------------------------------------------------------------------------------

2.4.6 Example
-------------
Temperature at feeder input   =     188.03 (deg-K)
Temperature at receiver input =     149.36 (deg-K)

2.5.6 and 2.5.7 Example 1 and 2 pp. 39-42.

%-------------------------------------------------------------------------------
fD            = 14; % GHz
pT            = DBSignal(100); % dBW
lFTX          = 0.5; % dB
thetaT        = 0.1; % deg
r             = 40000; % km
fIllumination = 70;
dT            = 4;
theta3DB      = Beamwidth3dB( dT, fD, fIllumination );
effR          = 0.55;
effT          = 0.6;
fIllumination = 70;
E             = 10;

Uplink antenna

%---------------
gT            = AntennaGain(struct('type','circular aperture','area',pi*dT^2/4,'eff',effT),fD);
lT            = LossDepointing( thetaT, theta3DB );
eIRPSL        = pT + gT - lT - lFTX;

lFS           = LossFreeSpace( fD, r );
lA            = LossAtmosphericGas( fD, E );
lU            = lA + lFS;

theta3DB      = 2;
gR            = Gain3dB( theta3DB, fIllumination, effR );
pR            = eIRPSL - lFS + gR;

lPol          = LossPolarization( 0 );
lFRX          = 1;
fR            = 3; % dB
tAmbient      = 290;
gOverT        = gR - fR - lFRX - DBSignal(TReceiver(290,1,290,290));
cOverN0       = eIRPSL - lU + gOverT - DBSignal(1.381e-23);
lRain         = LossPrecipitation( 14, 'k', 0.01, 20, 0, 0, [40000;0;0], 0 );
cOverN0Rain   = eIRPSL - lU - lRain + gOverT - DBSignal(1.381e-23);

fprintf('\n\n2.5.6 Example 1 and 2.5.7 Example 2: Uplink\n-------------------------------------------\n');
fprintf('Gain Transmit Antenna = %10.2f (dBW)\n',gT);
fprintf('Depointing Loss       = %10.2f (dB)\n',lFTX);
fprintf('Transmitter Loss      = %10.2f (dB)\n',lT);
fprintf('Theta 3 dB            = %10.2f (dB)\n',theta3DB);
fprintf('EIRP SL               = %10.2f (dBW)\n',eIRPSL);
fprintf('Free Space Loss       = %10.2f (dB)\n',lFS);
fprintf('Atmospheric Loss      = %10.2f (dB)\n',lA);
fprintf('Upward Path Loss      = %10.2f (dB)\n',lU);
fprintf('Rain Loss             = %10.2f (dB)\n',lRain);
fprintf('Gain Receive Antenna  = %10.2f (dB)\n',gR);
fprintf('Polarization Loss     = %10.2f (dB)\n',lPol);
fprintf('Receiver Loss         = %10.2f (dB)\n',lFRX);
fprintf('(G/T)SL               = %10.2f (dB)\n',gOverT);
fprintf('C/N0 Clear Sky        = %10.2f (dB)\n',cOverN0);
fprintf('C/N0 Rain             = %10.2f (dB)\n',cOverN0Rain);

%-------------------------------------------------------------------------------

2.5.6 Example 1 and 2.5.7 Example 2: Uplink
-------------------------------------------
Gain Transmit Antenna =      53.15 (dBW)
Depointing Loss       =       0.50 (dB)
Transmitter Loss      =       0.85 (dB)
Theta 3 dB            =       2.00 (dB)
EIRP SL               =      71.80 (dBW)
Free Space Loss       =     207.41 (dB)
Atmospheric Loss      =       0.29 (dB)
Upward Path Loss      =     207.70 (dB)
Rain Loss             =      11.80 (dB)
Gain Receive Antenna  =      38.23 (dB)
Polarization Loss     =      -0.00 (dB)
Receiver Loss         =       1.00 (dB)
(G/T)SL               =       6.59 (dB)
C/N0 Clear Sky        =      99.29 (dB)
C/N0 Rain             =      87.48 (dB)

2.5.8 and 2.5.9 Example 1 and 2 pp. 42-46.

%-------------------------------------------------------------------------------

Spacecraft

%-----------
fD            = 12; % GHz
pT            = DBSignal(10); % dBW
lFTX          = 1; % dB
theta3DB      = 2; % deg
effT          = 0.55;
r             = 40000; % km
E             = 10;
lA            = LossAtmosphericGas( fD, E );
fIllumination = 70;
lEdge         = 3; % dB

Ground station

%---------------
f             = 1;
lFRX          = 0.5;
tF            = 290;
dR            = 4;
effR          = 0.6;
thetaR        = 0.1;
tGround       = 45;

EIRP Satellite

%---------------
gT            = Gain3dB( theta3DB, fIllumination, effT );
eIRPSL        = pT + gT - lEdge - lFTX;

Downlink attenuation

%---------------------
lFS           = LossFreeSpace( fD, r );
lD            = lA + lFS;

Downlink system noise temperature

%----------------------------------
tSky          = TSky( 12, 10 );
tA            = TAntennaGround( tSky, 0, tGround );
tERX          = TAttenuator( tF, f );
tD            = TReceiver( tA, lFRX, tF, tERX );
gR            = AntennaGain(struct('type','circular aperture','area',pi*dT^2/4,'eff',effR),fD);

Pointing loss

%--------------
theta3DB      = Beamwidth3dB( dR, fD, fIllumination );
lP            = LossDepointing( thetaT, theta3DB );

gOverT        = gR - lP - lFRX - DBSignal(tD);
cOverN0       = eIRPSL - lD + gOverT - DBSignal(1.381e-23);

With rain

%----------

fprintf('\n\n2.5.8 Example 1 and 2.5.9 Example 2: Downlink\n-------------------------------------------\n');
fprintf('Gain Transmit Antenna             = %10.2f (dBW)\n',gT);
fprintf('Transmitter Loss                  = %10.2f (dB)\n',lFTX);
fprintf('Theta 3 dB                        = %10.2f (dB)\n',theta3DB);
fprintf('EIRP SL                           = %10.2f (dBW)\n',eIRPSL);
fprintf('Free Space Loss                   = %10.2f (dB)\n',lFS);
fprintf('Atmospheric Loss                  = %10.2f (dB)\n',lA);
fprintf('Downward Path Loss                = %10.2f (dB)\n',lD);
fprintf('Gain Receive Antenna              = %10.2f (dB)\n',gR);
fprintf('Receiver Loss                     = %10.2f (dB)\n',lFRX);
fprintf('Noise Temperature Receiver        = %10.2f (deg-K)\n',tERX);
fprintf('Pointing Loss                     = %10.2f (dB)\n',lP);
disp('---- Clear Sky ----')
fprintf('Noise Temperature Ground Antenna  = %10.2f (deg-K)\n',tA);
fprintf('Noise Temperature Downlink        = %10.2f (deg-K)\n',tD);
fprintf('(G/T)ES  Clear Sky                = %10.2f (dB)\n',gOverT);
fprintf('C/N0 Clear Sky                    = %10.2f (dB)\n',cOverN0);

disp('---- Rain ----')
lRain         = LossPrecipitation( fD, 'e', 0.01, 45, 0, 0, [40000;0;0], 0 );
tA            = TAntennaGround( tSky, lRain, tGround );
tD            = TReceiver( tA, lFRX, tF, tERX );
lD            = lA + lFS + lRain;

gOverT        = gR - lP - lFRX - DBSignal(tD);
cOverN0       = eIRPSL - lD + gOverT - DBSignal(1.381e-23);

fprintf('Rain Loss                         = %10.2f (dB)\n',lRain);
fprintf('Downward Path Loss                = %10.2f (dB)\n',lD);
fprintf('Noise Temperature Ground Antenna  = %10.2f (deg-K)\n',tA);
fprintf('Noise Temperature Downlink        = %10.2f (deg-K)\n',tD);
fprintf('(G/T)ES                           = %10.2f (dB)\n',gOverT);
fprintf('C/N0                              = %10.2f (dB)\n',cOverN0);


%--------------------------------------
% PSS internal file version information
%--------------------------------------

2.5.8 Example 1 and 2.5.9 Example 2: Downlink
-------------------------------------------
Gain Transmit Antenna             =      38.23 (dBW)
Transmitter Loss                  =       1.00 (dB)
Theta 3 dB                        =       0.44 (dB)
EIRP SL                           =      44.23 (dBW)
Free Space Loss                   =     206.07 (dB)
Atmospheric Loss                  =       0.29 (dB)
Downward Path Loss                =     206.36 (dB)
Gain Receive Antenna              =      51.81 (dB)
Receiver Loss                     =       0.50 (dB)
Noise Temperature Receiver        =      75.09 (deg-K)
Pointing Loss                     =       0.63 (dB)
---- Clear Sky ----
Noise Temperature Ground Antenna  =      65.02 (deg-K)
Noise Temperature Downlink        =     164.57 (deg-K)
(G/T)ES  Clear Sky                =      28.52 (dB)
C/N0 Clear Sky                    =      94.99 (dB)
---- Rain ----
Rain Loss                         =       3.63 (dB)
Downward Path Loss                =     209.99 (dB)
Noise Temperature Ground Antenna  =     209.44 (deg-K)
Noise Temperature Downlink        =     293.29 (deg-K)
(G/T)ES                           =      26.01 (dB)
C/N0                              =      88.85 (dB)