Main Content

p618SiteDiversityOutage

Calculate outage probability due to rain attenuation with site diversity

Since R2021a

    Description

    Outage = p618SiteDiversityOutage(cfgsd) returns the outage probability due to rain attenuation with site diversity. The function calculates this value as per the ITU-R P.618 recommendation [1].

    This function requires MAT-files with digital maps from International Telecommunication Union (ITU) documents. If they are not available on the path, download and uncompress the data files from https://www.mathworks.com/supportfiles/spc/P618/ITURDigitalMaps.tar.gz to a location on the MATLAB® path.

    example

    Outage = p618SiteDiversityOutage(cfgsd,Name,Value) specifies additional options using one or more name-value pair arguments.

    example

    Examples

    collapse all

    This example requires MAT-files with digital maps from ITU documents. If they are not available on the path, execute the following commands to download and untar the MAT-files.

    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

    Create a P.618 site diversity configuration object with a signal frequency of 25 GHz.

    cfgsd = p618SiteDiversityConfig;
    cfgsd.Frequency = 25e9;

    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.

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

    Calculate the outage probability due to rain attenuation with site diversity.

    outage = p618SiteDiversityOutage(cfgsd)
    outage = 0.0338
    

    This example requires MAT-files with digital maps from ITU documents. If they are not available on the path, execute these commands to download and untar the MAT-files.

    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

    Create a default P.618 site diversity configuration object. Change the signal frequency to 25 GHz.

    cfgsd = p618SiteDiversityConfig;
    cfgsd.Frequency = 25e9;

    Specify the separation between two sites as 50 km and the attenuation threshold on the two links as [9 9] dB.

    cfgsd.SiteDistance = 50;
    cfgsd.AttenuationThreshold = [9 9];

    Calculate the outage probability for the specified site diversity configuration.

    outage = p618SiteDiversityOutage(cfgsd,'RainAnnualExceedances',[0.01 0.05 0.2],...
                                    'RainProbability1',0.3,...
                                    'RainProbability2',0.5) 
    outage = 0.0339
    

    Input Arguments

    collapse all

    P.618 site diversity configuration required for the calculation of the outage probability due to rain attenuation, specified as a p618SiteDiversityConfig object.

    Name-Value Arguments

    Specify optional pairs of arguments as Name1=Value1,...,NameN=ValueN, where Name is the argument name and Value is the corresponding value. Name-value arguments must appear after other arguments, but the order of the pairs does not matter.

    Before R2021a, use commas to separate each name and value, and enclose Name in quotes.

    Example: 'RainAnnualExceedances',[0.01 0.02 0.03 0.05] specifies the average annual time percentage of excess for the rain attenuation.

    Average annual time percentage of excess for the rain attenuation, specified as the comma-separated pair consisting of 'RainAnnualExceedances' and a nonnegative vector. The values in this vector must be less than the probability of rain at the two sites.

    If the local data is not available as an input, the function uses [0.01 0.02 0.03 0.05 0.1 0.2 0.3 0.5 1 2 3 5] as the default vector.

    Data Types: double | single

    Rain attenuations (in dB) at site 1, specified as the comma-separated pair consisting of 'RainAttenuations1' and a nonnegative vector. This value specifies the rain attenuation exceeded for the percentages given in the RainAnnualExceedances name-value pair argument. The dimension of this value must match that of the RainAnnualExceedances.

    If the local data is not available as an input, the function uses the method as defined in section 2.2.1.1 of the ITU-R P.618 [1] recommendation to calculate the rain attenuations for site 1.

    Note

    If you do not specify RainAttenuations1, then RainAnnualExceedances must be in the range from 0.01% to 5%.

    Data Types: double | single

    Rain attenuations (in dB) at site 2, specified as the comma-separated pair consisting of 'RainAttenuations2' and a nonnegative vector. This value specifies the rain attenuation exceeded for the percentages given in the RainAnnualExceedances name-value pair argument. The dimension of this value must match that of the RainAnnualExceedances.

    If the local data is not available as an input, the function uses the method as defined in section 2.2.1.1 of the ITU-R P.618 recommendation to calculate the rain attenuations for site 2.

    Note

    If you do not specify RainAttenuations2, then RainAnnualExceedances must be in the range from 0.01% to 5%.

    Data Types: double | single

    Probability of (in %) rain for site 1, specified as the comma-separated pair consisting of 'RainProbability1' and a nonnegative scalar.

    If the local measured rainfall rate data is not available as an input, the function uses the digital maps as defined in ITU-R P.837 Annex 1 [2] to calculate the rain probability for the sites.

    Data Types: double | single

    Probability of (in %) rain for site 2, specified as the comma-separated pair consisting of 'RainProbability2' and a nonnegative scalar.

    If the local measured rainfall rate data is not available as an input, the function uses the digital maps as defined in ITU-R P.837 Annex 1 [2] to calculate the rain probability for the sites.

    Data Types: double | single

    Output Arguments

    collapse all

    Outage probability due to rain attenuation with site diversity, returned as a nonnegative scalar. This argument predicts the joint probability (Pr(A1a1, A2a2)), where the attenuation on the path of the site 1 must exceed a1 and the attenuation on the path of the site 2 must exceed a2.

    References

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

    [2] International Telecommunication Union, ITU-R Recommendation P.837 (06/2017).

    [3] International Telecommunication Union, ITU-R Recommendation P.1511 (08/2019).

    [4] International Telecommunication Union, ITU-R Recommendation P.1510 (06/2017).

    [5] International Telecommunication Union, ITU-R Recommendation P.836 (12/2017).

    [6] International Telecommunication Union, ITU-R Recommendation P.840 (08/2019).

    [7] International Telecommunication Union, ITU-R Recommendation P.453 (08/2019).

    [8] International Telecommunication Union, ITU-R Recommendation P.839 (09/2013).

    [9] International Telecommunication Union, ITU-R Recommendation P.838 (03/2005).

    Extended Capabilities

    Version History

    Introduced in R2021a