rainpl

RF signal attenuation due to rainfall using ITU model

Syntax

L = rainpl(range,freq,rainrate)

L = rainpl(range,freq,rainrate,elev)

L = rainpl(range,freq,rainrate,elev,tau)

L = rainpl(range,freq,rainrate,elev,tau,pct)

Description

L = rainpl(range,freq,rainrate) returns the signal attenuation, L, due to rain with a long-term statistical rain rate. In this syntax, attenuation is a function of signal path length, range, signal frequency, freq, and rain rate, rainrate. The path elevation angle and polarization tilt angles are assumed to be zero.

The rainpl function applies the International Telecommunication Union (ITU) rainfall attenuation model to calculate path loss of signals propagating in a region of rainfall [1]. The function applies when the signal path is contained entirely in a uniform rainfall environment. Rain rate does not vary along the signal path. The attenuation model applies only for frequencies at 1–1000 GHz.

L = rainpl(range,freq,rainrate,elev) also specifies the elevation angle, elev, of the propagation path.

L = rainpl(range,freq,rainrate,elev,tau) also specifies the polarization tilt angle, tau, of the signal.

L = rainpl(range,freq,rainrate,elev,tau,pct) also specifies the specified percentage of time, pct. pct is a scalar in the range of 0.001–1, inclusive. The attenuation, L, is computed from a power law using the long-term statistical 0.01% rain rate (in mm/h).

Examples

Signal Attenuation Due to Rainfall

Compute the signal attenuation due to rainfall for a 20 GHz signal over a distance of 10 km in light and heavy rain.

Propagate the signal in a light rainfall of 1 mm/hr.

rr = 1.0;
L = rainpl(10000,20.0e9,rr)

L = 
1.3009

Propagate the signal in a heavy rainfall of 10 mm/hr.

rr = 10.0;
L = rainpl(10000,20.0e9,rr)

L = 
8.1584

Signal Attenuation Due to Rainfall as Function of Frequency

Plot the signal attenuation due to a 20 mm/hr statistical rainfall for signals in the frequency range from 1 to 1000 GHz. The path distance is 10 km.

rr = 20.0;
freq = [1:1000]*1e9;
L = rainpl(10000,freq,rr);
semilogx(freq/1e9,L)
grid
xlabel('Frequency (GHz)')
ylabel('Attenuation (dB)')

Signal Attenuation Due to Rainfall as Function of Elevation Angle

Compute the signal attenuation due to heavy rain as a function of elevation angle. Elevation angles vary from 0 to 90 degrees. Assume a path distance of 100 km and a signal frequency of 100 GHz.

Set the rain rate to 10 mm/hr.

rr = 10.0;

Set the elevation angles, frequency, range.

elev = [0:1:90];
freq = 100.0e9;
rng = 100000.0*ones(size(elev));

Compute and plot the loss.

L = rainpl(rng,freq,rr,elev);
plot(elev,L)
grid
xlabel('Path Elevation (degrees)')
ylabel('Attenuation (dB)')

Signal Attenuation Due to Rainfall as Function of Polarization

Compute the signal attenuation due to heavy rainfall as a function of the polarization tilt angle. Assume a path distance of 100 km, a signal frequency of 100 GHz, and a path elevation angle of 0 degrees. Set the rainfall rate to 10 mm/hour. Plot the signal attenuation versus polarization tilt angle.

Set the polarization tilt angle to vary from -90 to 90 degrees.

tau = -90:90;

Set the elevation angle, frequency, path distance, and rain rate.

elev = 0;
freq = 100.0e9;
rng = 100e3*ones(size(tau));
rr = 10.0;

Compute and plot the attenuation.

L = rainpl(rng,freq,rr,elev,tau);
plot(tau,L)
grid
xlabel('Tilt Angle (degrees)')
ylabel('Attenuation (dB)')

Input Arguments

range — Signal path length in meters

nonnegative real-valued scalar | nonnegative real-valued M-by-1 column vector | nonnegative real-valued 1-by-M row vector

Signal path length in meters, specified as a nonnegative real-valued scalar, or as a M-by-1 or 1-by-M vector.

Example: [13000.0,14000.0]

freq — Signal frequency in Hz

positive real-valued scalar | nonnegative real-valued N-by-1 column vector | nonnegative real-valued 1-by-N row vector

Signal frequency in Hz, specified as a positive real-valued scalar, or as a nonnegative N-by-1 or 1-by-N vector. Frequencies must lie in the range 1–1000 GHz.

Example: [1400.0e6,2.0e9]

rainrate — Long-term statistical rain rate, in mm/h

nonnegative real-valued scalar

Long-term statistical rain rate, in mm/h, specified as a nonnegative real-valued scalar. The long-term statistical rain rate is the rain rate that is exceeded 0.01% of the time. You can adjust the percent of time using the pct argument.

Example: 1.5

elev — Signal path elevation angle, in degrees

0.0 (default) | real-valued scalar | real-valued M-by-1 column vector | real-valued 1-by-M row vector

Signal path elevation angle, in degrees, specified as a real-valued scalar, an M-by-1 vector, or a 1-by- M vector. Units are in degrees between –90° and 90°. If elev is a scalar, all propagation paths have the same elevation angle. If elev is a vector, its length must match the dimension of range and each element in elev corresponds to a propagation range in range.

Example: [0,45]

tau — Tilt angle of polarization ellipse, in degrees

0.0 (default) | real-valued scalar | real-valued M-by-1 column vector | real-valued 1-by-M row vector

Tilt angle of the signal polarization ellipse, in degrees, specified as a real-valued scalar, or as an M-by-1 or 1-by- M vector. Units are in degrees between –90° and 90°. If tau is a scalar, all signals have the same tilt angle. If tau is a vector, its length must match the dimension of range. In that case, each element in tau corresponds to a propagation path in range.

The tilt angle is defined as the angle between the semi-major axis of the polarization ellipse and the x-axis. Because the ellipse is symmetrical, a tilt angle of 100° corresponds to the same polarization state as a tilt angle of -80°. Thus, the tilt angle need only be specified between ±90°.

Example: [45,30]

pct — Exceedance percentage of rainfall, in %

0.01 (default) | positive scalar between 0.001 and 1

Exceedance percentage of rainfall, in %, specified as a positive scalar between 0.001 and 1. The long-term statistical rain rate is the rain rate that is exceeded pct of the time.

Data Types: double

Output Arguments

L — Signal attenuation in dB

real-valued M-by-N matrix

Signal attenuation in dB, returned as a real-valued M-by-N matrix. Each matrix row represents a different path where M is the number of paths. Each column represents a different frequency where N is the number of frequencies.

More About

Rainfall Attenuation Model

This model calculates the attenuation of signals that propagate through regions of rainfall. Rain attenuation is a dominant fading mechanism and can vary from location-to-location and from year-to-year.

Electromagnetic signals are attenuated when propagating through a region of rainfall. Rainfall attenuation is computed according to the ITU rainfall model Recommendation ITU-R P.838-3: Specific attenuation model for rain for use in prediction methods. The model computes the specific attenuation (attenuation per kilometer) of a signal as a function of rainfall rate, signal frequency, polarization, and path elevation angle. The specific attenuation, ɣ_R, is modeled as a power law with respect to rain rate

$${\gamma }_R=k{R}^{\alpha },$$

where R is rain rate. Units are in mm/hr. The parameter k and exponent α depend on the frequency, the polarization state, and the elevation angle of the signal path. The specific attenuation model is valid for frequencies from 1–1000 GHz.

To compute the total attenuation for narrowband signals along a path, the function multiplies the specific attenuation by the an effective propagation distance, d_eff. Then, the total attenuation is L = d_effγ_R.

The effective distance is the geometric distance, d, multiplied by a scale factor

$$r=\frac{1}{0.477d^{0.633}{R}_{0.01}^{0.073\alpha }{f}^{0.123}-10.579\left(1-\exp \left(-0.024d\right)\right)}$$

where f is the frequency. The article Recommendation ITU-R P.530-17 (12/2017): Propagation data and prediction methods required for the design of terrestrial line-of-sight systems presents a complete discussion for computing attenuation.

The rain rate, R, used in these computations is the long-term statistical rain rate, R_0.01. This is the rain rate that is exceeded 0.01% of the time. The calculation of the statistical rain rate is discussed in Recommendation ITU-R P.837-7 (06/2017): Characteristics of precipitation for propagation modelling. This article also explains how to compute the attenuation for other percentages from the 0.01% value.

You can also apply the attenuation model to wideband signals. First, divide the wideband signal into frequency subbands and apply attenuation to each subband. Then, sum all attenuated subband signals into the total attenuated signal.

References

[1] Radiocommunication Sector of International Telecommunication Union. Recommendation ITU-R P.838-3: Specific attenuation model for rain for use in prediction methods. 2005.

[2] Radiocommunication Sector of International Telecommunication Union. Recommendation ITU-R P.530-17: Propagation data and prediction methods required for the design of terrestrial line-of-sight systems. 2017.

[3] Recommendation ITU-R P.837-7: Characteristics of precipitation for propagation modelling

[4] Seybold, J. Introduction to RF Propagation. New York: Wiley & Sons, 2005.

Extended Capabilities

C/C++ Code Generation

Generate C and C++ code using MATLAB® Coder™.

Usage notes and limitations:

Does not support variable-size inputs.

Version History

Introduced in R2016a