Engee documentation

Wideband LOS Channel

Broadband propagation channel in line-of-sight conditions.

wideband los channel

Description

The Wideband LOS Channel block simulates the propagation of signals from one point in space to multiple points or from multiple points back to one point via line-of-sight (LOS) channels. The block simulates propagation time, free-space propagation loss, Doppler shift, atmospheric and climatic losses. The modelling assumes that the propagation velocity of the signal is much greater than the velocity of the object, in which case the stop-and-hop model is valid.

When the signal propagates through the LOS channel to and from the object, one block can be used to calculate the two-way propagation delay (back and forth) in the LOS channel or two blocks can be used to calculate the one-way propagation delay in each direction. Since the free-space propagation delay is not necessarily an integer multiple of the sampling interval, it may be that the total delay in both directions in the calculation when using a two-way propagation block is different from the delay in the calculation when using two one-way propagation blocks. For this reason, it is recommended that a single two-way propagation unit be used whenever possible.

Ports

Input

X - radiated signal
complex vector M by 1 | complex matrix M by N

The emitted signal in the form of a complex vector-column M by 1 or a complex matrix M by N. The value M is the duration of the signal, and N is the number of elements of the receive point array (or subarray, if subarrays are supported).

Dimension Signal

Vector-column M by 1

The same signal is emitted from all elements of the array (subarray).

Matrix M by N

Each column corresponds to the signal emitted by the corresponding element of the array (subarray).

The size of the first dimension of the input matrix may be varied to simulate the varying duration of the signal. The size variation may occur, for example, in the case of a pulsed signal with a variable pulse repetition rate.

Data types: Float16, Float32, Float64, Int8, Int16, Int32, Int64, UInt8, UInt16, UInt32, UInt64, Bool.

Support for complex numbers: Yes

Pos1 is the position of the signal transmitter.
real vector-column 3 by 1 | real matrix 3 by N

The signal position is specified as a real vector-column 3 by 1 or a real matrix 3 by N. The value N is the number of source positions. Examples of source positions are coordinates of transmitters, array elements or subarrays. The units of measurement of position coordinates are metres.

If Pos1 is a column vector, it has dimensionality 3 by 1. If Pos1 is a matrix, each column defines a different source position and has dimensionality 3 by N.

If Pos1 has more than one column, then Pos2 can have only one column. Pos1 and Pos2 cannot both be specified as matrices.

Data types: Float16, Float32, Float64, Int8, Int16, Int32, Int64, UInt8, UInt16, UInt32, UInt64.

Pos2 - position of the signal receiver
real vector-column 3 by 1 | real matrix 3 by N

The position of the signal receiver is specified as a real vector-column 3 by 1 or a real matrix 3 by N. The value N is the number of receivers, for example, the positions of array or subarray elements. The units of measurement of position coordinates are metres.

If Pos2 is a column vector, it has dimensionality 3 by 1. If Pos2 is a matrix, each column defines a different position of the signal receiver and has dimensionality 3 by N.

If Pos1 has more than one column, then Pos2 can have only one column. Pos1 and Pos2 cannot both be specified as matrices.

Data types: Float16, Float32, Float64, Int8, Int16, Int32, Int64, UInt8, UInt16, UInt32, UInt64.

Vel1 - speed of the signal transmitter
real vector-column 3 by 1 | real matrix 3 by N

Transmitter velocity given as a real vector or matrix having the same dimensionality as Pos1.

Data types: Float16, Float32, Float64, Int8, Int16, Int32, Int64, UInt8, UInt16, UInt32, UInt64.

Vel2 - speed of the signal receiver
real vector-column 3 by 1 | real matrix 3 by N | real matrix N by 3

The receiver velocity of a signal, given as a real vector or matrix having the same dimensionality as Pos2.

Data types: Float16, Float32, Float64, Int8, Int16, Int32, Int64, UInt8, UInt16, UInt32, UInt64.

Output

Port_1 - propagated signal
complex vector-column M by 1 | complex matrix M by N

The propagated signal returned as a complex vector-column M by 1 or a complex matrix M by N. Port_1 has the same dimensionality as the input port X: M is the signal duration and N is the number of signals.

Data types: Float16, Float32, Float64, Int8, Int16, Int32, Int64, UInt8, UInt16, UInt32, UInt64, Bool.

Support for complex numbers: Yes

Parameters

Propagation speed (m/s) - speed of signal propagation
3e8 (by default) | positive scalar

Signal propagation speed as a real positive scalar.

By default, the value of the speed of light is 3e8.

The unit of measurement is m/c.

Data types: Float16, Float32, Float64, Int8, Int16, Int32, Int64, UInt8, UInt16, UInt32, UInt64.

Signal carrier frequency (Hz) - carrier frequency of the signal
3e8 (By default) | positive scalar

Carrier frequency of the signal as a positive real scalar. The unit of measurement is Hz.

Data types: Float16, Float32, Float64, Int8, Int16, Int32, Int64, UInt8, UInt16, UInt32, UInt64.

Number of subbands - number of subbands
68 (By default) | Positive integer

The number of subbands to be processed, specified as a positive integer.

Data types: Int8, Int16, Int32, Int64, UInt8, UInt16, UInt32, UInt64.

Specify atmospheric parameters - take into account signal attenuation in the atmosphere
Off (By default) | On

Select this check box to enable atmospheric attenuation accounting.

When selected, the parameters Temperature (degrees Celsius), Dry air pressure (Pa), Water vapour density (g/m^3), Liquid water density (g/m^3), and Rain rate (mm/h) appear in the dialogue box.

Temperature (degrees Celsius) is the ambient temperature
15 (By default) | scalar.

Ambient temperature specified as a real scalar. The unit of measurement is degrees Celsius.

Dependencies

To use this parameter, select the Specify atmospheric parameters checkbox.

Dry air pressure (Pa) is the atmospheric pressure of dry air
101325 (By default) | positive scalar

Atmospheric pressure of dry air specified as a positive real scalar. The value of this parameter by default corresponds to one standard atmosphere. The unit of measurement is Pa.

Dependencies

To use this parameter, select the Specify atmospheric parameters checkbox.

Water vapour density (g/m^3) - the density of water vapour in the atmosphere
7.5 (by default) | `positive scalar'.

The density of water vapour in the atmosphere, given as a positive real scalar. The unit of measurement is g/m3.

Dependencies

To use this parameter, select the Specify atmospheric parameters checkbox.

Liquid water density (g/m^3) - liquid water density
0.0 (by default) | `non-negative scalar'.

The density of liquid water in fog or clouds, given as a non-negative real scalar. The unit of measurement is g/m3. Typical values of liquid water density are 0.05 for medium fog and 0.5 for dense fog.

Dependencies

To use this parameter, select the Specify atmospheric parameters checkbox.

Rain rate (mm/h) - rainfall intensity
0.0 (By default) | `non-negative scalar'.

Precipitation intensity specified as a non-negative real scalar. The unit of measurement is mm/h.

Dependencies

To use this parameter, select the Specify atmospheric parameters checkbox.

Perform two-way propagation - enable two-way propagation
Off (By default) | `on

Select this check box to perform two-way propagation between the signal transmitter and receiver. Otherwise, the unit performs one-way propagation from the transmitter to the signal receiver.

Inherit sample rate - sample rate inheritance
On (By default) | Off

Select the checkbox to inherit sample rate from upstream blocks. Otherwise, set the sample rate using the Sample rate (Hz) parameters.

Sample rate (Hz) - sampling rate
1e6 (By default) | Positive scalar

Sampling frequency of the signal as a positive scalar. The unit of measurement is Hz.

Dependencies

To use this parameters, uncheck the Inherit sample rate checkbox.

Data types: Float16, Float32, Float64, Int8, Int16, Int32, Int64, UInt8, UInt16, UInt32, UInt64.

Maximum one-way propagation distance (m) - maximum one-way propagation distance
10e3 (by default) | positive scalar

The maximum distance in metres between the transmitter and receiver of a signal in the form of a positive real scalar. The unit of measurement is m. The amplitudes of any signals that propagate beyond this distance are zero.

Algorithms

Attenuation and loss factors

The attenuation or signal loss in a wideband LOS channel consists of four components:

,

where

  • - is the attenuation of the signal when propagating in free space.

  • - signal attenuation in atmospheric propagation.

  • - signal attenuation for propagation due to fog and clouds.

  • - signal attenuation for propagation due to the presence of precipitation.

Each component is measured in units of magnitude, not dB.

Propagation delay, Doppler shift and free-space loss

When the signal source and receiver are stationary relative to each other, you can write the output of the free-space channel as τ , where τ is the signal delay and is the free-space propagation loss. The signal delay is calculated as τ , where is the propagation distance and is the propagation velocity. The free space path loss is defined by the expression

,

where is the wavelength of the signal.

This formula assumes that the target is in the far field of the transmitting element or array. In the near field, the free-space propagation path loss formula is invalid and can result in loss values less than unity, which is equivalent to signal amplification. For this reason, the loss is set to unity for the values of .

If there is relative motion between the source and receiver, the Doppler frequency shift is taken into account. The frequency shift is for one-way propagation and for two-way propagation. The value is the relative velocity of the receiver relative to the source.

Model of signal attenuation in the atmosphere

This model calculates the attenuation of signals propagating through atmospheric gases.

Electromagnetic signals are attenuated when propagating through the atmosphere. This effect is mainly due to the resonant absorption lines of oxygen and water vapour, with a smaller contribution from nitrogen gas. The model also includes a continuous absorption spectrum below 10 GHz. The ITU (International Telecommunication Union) model given in [1] is used for the calculations. The model calculates specific attenuation (attenuation per kilometre) as a function of temperature, pressure, water vapour density and signal frequency. The atmospheric gas model is valid for frequencies 1-1000 GHz and is applicable to polarised and unpolarised fields.

The formula for the specific attenuation at each frequency is as follows:

.

The value is the imaginary part of the complex refractive index of the atmosphere and consists of a spectral line component and a continuous component:

.

The spectral component consists of the sum of discrete spectral terms, which are the product of the localised frequency band function and the spectrum line intensity . For atmospheric oxygen, the line intensity of the spectrum is:

.

For atmospheric water vapour, the intensity of the spectrum line is:

,

where:

  • - dry air pressure.

  • - partial pressure of water vapour.

  • - ambient temperature.

The units of pressure are hectopascals (hPa) and the units of temperature are degrees Kelvin.

The partial pressure of water vapour is related to the density of water vapour , as follows:

.

The total atmospheric pressure is .

For each oxygen line, depends on two parameters: and . Similarly, each water vapour line depends on two parameters: and .

The localised bandwidth functions are complex functions of frequency. These functions depend on the empirical parameters of the model.

The calculation of the total attenuation of a narrowband signal is done by multiplying the specific attenuation by the path length . Then the total attenuation is equal to:

.

This model of signal attenuation can be applied to broadband signals. To do this, first divide the broadband signal into frequency sub-bands, calculate the attenuation for each sub-band, and then sum all the attenuated sub-band signals into a total attenuated signal.

Fog and Cloud Attenuation Model

This model calculates the attenuation of signals propagating through fog or clouds.

Fog or clouds are the same atmospheric phenomenon. The ITU model given in [2] is used for the calculations. The model calculates the specific attenuation (attenuation per kilometre) of a signal as a function of liquid water density, signal frequency and temperature. The model is applicable to polarised and unpolarised fields. The formula for the specific attenuation at each frequency is as follows:

,

where:

  • - density of liquid water in gm/m3.

  • - is the specific attenuation coefficient and depends on frequency.

The model of attenuation in clouds and fog is valid for frequencies 10-1000 GHz. The units of the specific attenuation coefficient are (dB/km)/(g/m3).

The total attenuation of a narrowband signal is calculated by multiplying the specific attenuation by the path length . Then the total attenuation is equal to:

.

This model of signal attenuation can be applied to broadband signals. To do this, first divide the broadband signal into frequency sub-bands, calculate the attenuation for each sub-band, and then sum all the attenuated sub-band signals into a total attenuated signal.

Signal attenuation model in the presence of precipitation

This model calculates the attenuation of signals propagating through areas of precipitation. Precipitation attenuation is the dominant attenuation mechanism and can vary from location to location and from year to year.

Electromagnetic signals are attenuated as they propagate through the precipitation region. The ITU model given in [3] is used for the calculations. The model calculates the specific attenuation (attenuation per kilometre) of a signal as a function of precipitation intensity, signal frequency, polarisation and elevation angle. The specific attenuation γ depends on the precipitation intensity according to a power law:

,

where:

  • - precipitation intensity. The units are mm/h.

  • parameters and degree index depend on frequency, polarisation state and signal path elevation angle.

This attenuation model is valid for frequencies 1-1000 GHz.

The calculation of the total attenuation of a narrowband signal is done by multiplying the specific attenuation by the effective propagation distance, . Then the total attenuation is equal to:

.

The effective distance is the geometric distance , multiplied by a scaling factor:

,

where is the frequency. A more detailed description of the attenuation calculation is given in [4].

The precipitation intensity , used in these calculations, is the long-term statistical precipitation intensity [5]. This is the precipitation intensity that is exceeded 0.01 % of the time.

The attenuation model can be applied to broadband signals. First divide the broadband signal into frequency sub-bands and apply the model to each sub-band. Then sum all attenuated sub-band signals into a total attenuated signal.

Processing the frequency sub-bands

Sub-band processing divides a wideband signal into several sub-bands and applies narrowband processing to the signal in each sub-band. The signals from all sub-bands are summed to form the output signal.

When using wideband system objects or blocks, the number of sub-bands , into which the wideband signal is to be decomposed, is specified. The centre frequencies and sub-band widths are automatically calculated based on the total bandwidth and the number of sub-bands. The total bandwidth is centred on the carrier or operating frequency . The total bandwidth is determined by the sampling frequency . The sub-band frequency bandwidth is defined as . The centre frequencies of the sub-bands are defined as

- if is odd,

- if is odd.

Some system objects allow to get the centre frequencies of sub-bands as output data when the object is started. The returned sub-band frequencies are ordered according to the order of the discrete Fourier transform. Frequencies above the carrier are displayed first, then frequencies below the carrier.

References

  1. Recommendation ITU-R P.676-10: Attenuation by atmospheric gases.

  2. Recommendation ITU-R P.840-6: Attenuation due to clouds and fog.

  3. Recommendation ITU-R P.838-3: Specific attenuation model for rain for use in prediction methods.

  4. Recommendation ITU-R P.530-17 (12/2017): Propagation data and prediction methods required for the design of terrestrial line-of-sight systems.

  5. Recommendation ITU-R P.837-7 (06/2017): Characteristics of precipitation for propagation modelling.