Engee documentation

Wideband LOS Channel

Broadband distribution channel in line-of-sight conditions.

blockType: WidebandLOSChannel

Path in the library:

/Phased Array Systems/Environment and Target/Wideband LOS Channel

Description

Block Wideband LOS Channel simulates the propagation of signals from one point in space to several points or from several points back to one point via line-of-sight (LOS) channels. The block simulates propagation time, propagation losses in free space, Doppler shift, atmospheric and climatic losses. When modeling, it is assumed that the speed of signal propagation is much higher than the speed of the object, and in this case, the "stop-and-hop" model is valid.

When a signal propagates through the LOS channel to an object and back, 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 propagation delay in free space is not necessarily an integer multiple of the sampling interval, it may turn out that the total delay in both directions in calculations using a two-way propagation unit differs from the delay in calculations using two one-way propagation units. For this reason, it is recommended to use a single two-way distribution unit whenever possible.

Ports

Entrance

X — radiated pass signal:q[<br>] complex vector M by 1 | complex matrix M by N

The emitted signal is in the form of a complex column vector 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 array of reception points (or a subarray, if subarrays are supported).

Dimension The signal

Column vector M by 1

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

The M by N matrix

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 can be varied to simulate the changing duration of the signal. A change in size may occur, for example, in the case of a pulse 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 — position of the pass signal transmitter:q[<br>] real column vector 3 by 1 | real matrix 3 by N

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

If Pos1 is a column vector, then it has dimension 3 by 1. If Pos1 is a matrix, then each of its columns defines its position of the signal source and has dimension 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 pass signal receiver:q[<br>] real column vector 3 by 1 | real matrix 3 by N

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

If Pos2 is a column vector, then it has a dimension of 3 by 1. If Pos2 is a matrix, then each column defines its position of the signal receiver and has a dimension of 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 — pass signal transmitter speed:q[<br>] real column vector 3 by 1 | real matrix 3 by N

The speed of the transmitter, set as a real vector or matrix having the same dimension as Pos1.

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

Vel2 — pass signal receiver speed:q[<br>] real column vector 3 by 1 | real matrix 3 by N | real matrix N by 3

The speed of the signal receiver, specified as a real vector or matrix having the same dimension as Pos2.

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

Output

Port_1 — propagated pass signal:q[<br>] complex column vector M by 1 |complex matrix M by N

The propagated signal returned as a complex column vector M by 1 or a complex matrix M by N. Port_1 has the same dimension as the input port X: M is the duration of the signal, 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 propagation of the pass signal:q[<br>] 3e8 (default) | positive scalar

The propagation velocity of the signal in the form of a real positive scalar.

The default value is the speed of light: `3e8'.

The units of measurement are m/s.

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

Signal carrier frequency (Hz) — carrier frequency of the pass signal:q[<br>] 3e8 (default) | positive scalar

The carrier frequency of the signal in the form of a positive real scalar. The units of measurement are Hz.

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

Number of subbands — number of pass subbands:q[<br>] 68 (default) | positive integer

The number of processed sub-ranges, set as a positive integer.

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

Specify atmospheric parameters — accounting for signal attenuation in the atmosphere
disabled (by default) | enabled

Select this option to enable recording of signal attenuation in the atmosphere.

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

Temperature (degrees Celsius) — ambient temperature pass:Q[<br>] 15 (default) | scalar

The ambient temperature, set 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) — atmospheric pressure of dry air
101325 (default) | positive scalar

Atmospheric pressure of dry air, given as a positive real scalar. The default value of this parameter 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 (default) | positive scalar

The density of water vapor 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) — density of liquid water
0.0 (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 for the density of liquid water are 0.05 for medium fog and 0.5 for thick fog.

Dependencies

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

Rain rate (mm/h) — precipitation intensity
0.0 (default) | non-negative scalar

Precipitation intensity, set 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
disabled (by default) | enabled

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

Inherit sample rate — inheritance of the sample rate
enabled (by default) | disabled

Check the box to inherit the sampling rate from higher-level blocks. Otherwise, set the sampling rate using the Sample rate (Hz) parameter.

Sample rate (Hz) — pass sampling rate:q[<br>] 1e6 (default) | positive scalar

The sampling frequency of the signal in the form of a positive scalar. The unit of measurement is Hz.

Dependencies

To use this option, 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 (default) | positive scalar

The maximum distance in meters between the transmitter and receiver of the 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 turn to zero.

Algorithms

Attenuation and loss factors

Signal attenuation or loss in a LOS broadband channel consists of four components:

,

where

  • — attenuation of the signal during propagation in free space.

  • — attenuation of the signal during propagation in the atmosphere.

  • — attenuation of the signal during propagation due to the presence of fog and clouds.

  • — signal attenuation during propagation due to precipitation.

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

Propagation delay, Doppler shift, and free space loss

When the signal source and receiver are stationary relative to each other, the output signal of the free space channel can be written as τ , where τ — signal delay, and — losses during propagation in free space. The signal delay is calculated as τ , where is the propagation distance, and — the speed of propagation. The losses on the way in free space are determined by the expression

,

where — 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 formula for losses along the propagation path in free space is invalid and can lead to loss values less than one, which is equivalent to signal amplification. For this reason, for the values losses are set at one.

If there is relative motion between the source and receiver, then the Doppler frequency shift is taken into account. The frequency shift is for unilateral distribution and for two-way distribution. Value — this is the relative speed of the receiver relative to the source.

A model of signal attenuation in the atmosphere

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

Electromagnetic signals weaken as they propagate through the atmosphere. This effect is mainly due to the lines of resonant absorption of oxygen and water vapor, 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 calculations. The model calculates the specific attenuation (attenuation per kilometer) as a function of temperature, pressure, water vapor density, and signal frequency. The atmospheric gas model is valid for frequencies of 1-1000 GHz and is applicable to polarized and unpolarized fields.

The formula for the specific attenuation at each frequency is:

.

Value It is an 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 a localized function of the frequency band and the intensity of the spectrum line . For atmospheric oxygen, the intensity of the spectral line is:

.

For atmospheric water vapor, the intensity of the spectral line is:

,

where:

  • — dry air pressure.

  • — partial pressure of water vapor.

  • — ambient temperature.

The units of measurement for pressure are hectopascals (gPa), and temperatures are degrees Kelvin.

Partial pressure of water vapor it is related to the density of water vapor , as follows:

.

The total atmospheric pressure is .

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

Localized frequency band functions they are complex functions of frequency. These functions depend on the empirical parameters of the model.

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

.

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

A model of signal attenuation in fog and clouds

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 calculations. The model calculates the specific attenuation (attenuation per kilometer) of a signal as a function of liquid water density, signal frequency, and temperature. The model is applicable to polarized and unpolarized fields. The formula for the specific attenuation at each frequency is:

,

where:

  • — the density of liquid water in gm/m ^3 ^.

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

The attenuation model in clouds and fog is valid for frequencies of 10-1000 GHz. The units of measurement 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:

.

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

A model of signal attenuation in the presence of precipitation

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

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

,

where:

  • — precipitation intensity. The units of measurement are mm/h.

  • parameter and the exponent They depend on the frequency, the state of polarization, and the elevation angle of the signal path.

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

The total attenuation of a narrowband signal is calculated by multiplying the specific attenuation by the effective propagation distance, . Then the total attenuation is:

.

The effective distance is the geometric distance multiplied by the scale factor:

,

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

Precipitation intensity , used in these calculations, represents the long-term statistical intensity of precipitation [5]. This is the intensity of precipitation, which is exceeded in 0.01% of cases.

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 the attenuated subband signals into a total attenuated signal.

Processing of sub-bands of frequencies

Subband processing divides the broadband signal into several subbands and applies narrowband processing to the signal in each subband. The signals of all sub-bands are combined to form an output signal.

When using broadband system objects or blocks, the number of sub-ranges is set. , into which the broadband signal needs to be decomposed. The center frequencies and subband widths are automatically calculated based on the total bandwidth and the number of subbands. The total frequency band is centered on the carrier or operating frequency . The total bandwidth is determined by the sampling rate . The width of the frequency subband is defined as . The central frequencies of the sub-bands are defined as

— if - even number,

— if "odd."

Some system objects allow you to get the central frequencies of the sub-bands as output data when starting the object. The returned subband frequencies are ordered according to the order of the discrete Fourier transform. The frequencies above the carrier are displayed first, then the frequencies below the carrier.

  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.