p618SiteDiversityConfig

Create P.618 site diversity configuration object

Description

The p618SiteDiversityConfig object sets P.618 site diversity configuration parameters required for the calculation of outage probability due to rain attenuation, as defined in the ITU-R P.618-14 recommendation [1].

Note

If you are using MATLAB^® in version R2025b or earlier, refer to the Tips section for instructions on how to run the examples by using ITU-R P.618-13.

Creation

Syntax

cfgSD = p618SiteDiversityConfig

cfgSD = p618SiteDiversityConfig(Name,Value)

Description

cfgSD = p618SiteDiversityConfig creates a P.618 site diversity configuration object with default property values.

cfgSD = p618SiteDiversityConfig(Name,Value) specifies Properties using one or more name-value pair arguments. Enclose each property name in quotes. For example, p618SiteDiversityConfig('Frequency',14.25e9,'ElevationAngle',[52.4099 52.4852]) configures a P.618 site diversity configuration object with a 14.25 GHz signal frequency and an elevation angle for two sites as [52.4099 52.4852].

Properties

Frequency — Signal frequency

14.25e9 (default) | scalar in the range [1e9, 55e9]

Signal frequency in Hz, specified as a scalar in the range [1e9, 55e9].

Data Types: double | single

ElevationAngle — Elevation angle of two sites

[52.4099 52.4852] (default) | two-element vector of values in the range [0, 90]

Elevation angle of the two sites in degrees, specified as a two-element vector of values in the range [0, 90].

Data Types: double | single

Latitude — Latitude of two sites

[25.768 25.463] (default) | two-element vector of values in the range [-90, 90]

Latitude of the two sites in degrees, specified as a two-element vector of values in the range [-90, 90]. A positive value corresponds to a North latitude, and a negative value corresponds to a South latitude.

Data Types: double | single

Longitude — Longitude of two sites

[-80.205 -80.486] (default) | two-element vector of values in the range [-180, 180]

Longitude of the two sites in degrees, specified as a two-element vector of values in the range [-180, 180]. A positive value corresponds to East longitude, and a negative value corresponds to West longitude.

Data Types: double | single

PolarizationTiltAngle — Polarization tilt angle for two sites

[0 0] (default) | two-element vector of values in the range [-90, 90]

Polarization tilt angle for the two sites in degrees, specified as a two-element vector of values in the range [-90, 90].

Data Types: double | single

SiteDistance — Separation between two sites

44.0256 (default) | positive scalar

Separation between the two sites in km, specified as a positive scalar.

Data Types: double | single

AttenuationThreshold — Attenuation threshold on two links

[9 3] (default) | two-element vector

Attenuation threshold on the two links in dB, specified as a two-element vector. The attenuation threshold on an earth space link is the maximum allowed attenuation on the path. Any attenuation value above this property value is considered an outage in the link.

Data Types: double | single

Object Functions

Specific to This Object

p618SiteDiversityOutage

Calculate outage probability due to rain attenuation with site diversity

Examples

Create P.618 Site Diversity Configuration Object

Create a default P.618 site diversity configuration object.

cfg = p618SiteDiversityConfig;

Specify the polarization tilt angles for two sites as [-90 90] degrees, separation between the two sites as 50 km, and attenuation threshold on the two links as [9 9] dB.

cfg.PolarizationTiltAngle = [-90 90];
cfg.SiteDistance = 50;
cfg.AttenuationThreshold = [9 9];

Set the direction of each earth station.

cfg.Latitude = [30 60];    % North direction
cfg.Longitude = [120 150]; % East direction

Display the properties of the configuration object.

disp(cfg);

  p618SiteDiversityConfig with properties:

                Frequency: 1.4500e+10
           ElevationAngle: [52.4099 52.4852]
                 Latitude: [30 60]
                Longitude: [120 150]
    PolarizationTiltAngle: [-90 90]
             SiteDistance: 50
     AttenuationThreshold: [9 9]

Outage Probability Due to Rain Attenuation with Site Diversity

This example requires a MAT file with digital maps from ITU documents. If the MAT file is not available on the path, execute the following commands to download and extract it.

matFile = exist('iturp618v14Maps.mat','file');
if ~matFile
    if ~exist('ITURP618v14DigitalMaps.tar.gz','file')
        url = 'https://www.mathworks.com/supportfiles/spc/P618/ITURP618v14DigitalMaps.tar.gz';
        websave('ITURP618v14DigitalMaps.tar.gz',url)
        untar('ITURP618v14DigitalMaps.tar.gz')
    else
        untar('ITURP618v14DigitalMaps.tar.gz')
    end
    addpath(cd)
end

This example shows how to calculate outage probability due to rain attenuation with site diversity.

The p618SiteDiversityOutage function applies to unbalanced and balanced systems and computes the joint probability of exceeding attenuation thresholds.

Configure P.618 Site Diversity Parameters

Create a default P.618 site diversity configuration object. Change property values, and then display the object properties.

cfgSD = p618SiteDiversityConfig;
cfgSD.Frequency = 25e9;                 % Signal frequency in Hz   
cfgSD.Latitude = [30 60];               % North direction
cfgSD.Longitude = [120 150];            % East direction
cfgSD.PolarizationTiltAngle = [-90 90];
cfgSD.AttenuationThreshold = [7 7];     % Attenuation threshold on the two links
cfgSD.SiteDistance = 50;                % Separation between the two sites
disp(cfgSD);

  p618SiteDiversityConfig with properties:

                Frequency: 2.5000e+10
           ElevationAngle: [52.4099 52.4852]
                 Latitude: [30 60]
                Longitude: [120 150]
    PolarizationTiltAngle: [-90 90]
             SiteDistance: 50
     AttenuationThreshold: [7 7]

Calculate Outage Probability

Calculate the outage probability due to rain attenuation for the specified site diversity configuration.

outage = p618SiteDiversityOutage(cfgSD, ...
                                RainAnnualExceedances=[0.01 0.01 0.03 0.05 0.1 0.2], ...
                                RainProbability1=0.3, ...
                                RainProbability2=0.4);
disp(outage);

    0.0030

Tips

For MATLAB version R2025b or earlier, the p618SiteDiversityConfig object uses ITU-R P.618-13 for its calculations instead of ITU-R P.618-14.

Use this code to download and unzip the MAT files containing digital maps from the ITU-R P.618-13 documents. The code also adds the extracted files to your MATLAB path.

maps = exist('maps.mat','file');
p836 = exist('p836.mat','file');
p837 = exist('p837.mat','file');
p840 = exist('p840.mat','file');
matFiles = [maps p836 p837 p840];
if ~all(matFiles)
    if ~exist('ITURDigitalMaps.tar.gz','file')
        url = 'https://www.mathworks.com/supportfiles/spc/P618/ITURDigitalMaps.tar.gz';
        websave('ITURDigitalMaps.tar.gz',url);
        untar('ITURDigitalMaps.tar.gz');
    else
        untar('ITURDigitalMaps.tar.gz');
    end
    addpath(cd);
end

To run the examples, replace the corresponding code fragment that downloads and unzips the latest digital maps from ITU-R P.618-14 in each example with this to ensure compatibility with ITU-R P.618-13.

References

[1] International Telecommunication Union, ITU-R Recommendation P.618 (08/2023).

Extended Capabilities

C/C++ Code Generation
Generate C and C++ code using MATLAB® Coder™.

Version History

Introduced in R2021a

R2026a: Support for ITU-R P.618-14

The p618SiteDiversityConfig feature now uses ITU-R P.618-14 for its calculations.