Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

MATHEMATICAL MODELLING OF SPATIAL-TEMPORAL STRUCTURE INTERFERENCE INFLUENCE FROM THE AERODYNAMIC TRANSMITTERS

Sagdeev K.M. 1 Sagdeevа E.K. 1
1 Federal State Autonomous Educational Institution Higher Professional Education «North-Caucasian Federal University»
The article deals with the scientific problem related to the problem of electronic protection channels of the trunked radio communication system in the conditions of the massive use of aerodynamic transmitters of interference. The purpose of the work is to obtain analytical dependence of the intensity distribution of the interference power at the antenna input suppressed a receiving device from elevation angle for a specific point in time. Solution of the scientific task is accomplished by statistical modeling of spatial-temporal distribution of angle-of-arrival interference in the vertical plane, with subsequent mathematical modeling of spatial-temporal structure of the interference generated by the aerodynamic group of transmitters of interference. To account for the influence of the medium of propagation imposed by its transfer function in the vertical plane. In the simulation result with a high confidence probability to get the average angular power spectrum of the interference in the vertical plane. The results of the work recommended for spatial selection interference from aerodynamic transmitters of interference.
aerodynamic transmitters of interference
spatial-temporal structure of the interference
angular power spectrum of the interference
mathematical model

Постановка научной задачи

Одной из актуальных задач информационной безопасности является радиоэлектронная защита (РЭЗ) радиоканалов, по которым осуществляется оперативное управление распределенными критически важными объектами. Анализ информационных источников [1–4] показывает, что подавление таких каналов, реализуемых на базе систем транкинговой радиосвязи, возможно за счет аэродинамических передатчиков помех (АПП), создающих широкополосные заградительные помехи. В настоящее время АПП могут быть легко реализованы с помощью беспилотных летательных аппаратов (БЛА) в виде дронов (квадрокоптеров) или парашютируемых передатчиков помех (ПП). Свое предназначение маломощные АПП осуществляют путем их барражирования и зависания в районе постановки помех с включенными на излучение ПП. Предполагается, что для эффективного РЭП транкинговой связи в определенной области пространства будет осуществляться массированное применение АПП, располагаемых в виде эквидистантной решетки относительно подавляемых радиоприемных устройств (РПУ), прежде всего базовых станций [1, 2].

Анализ особенностей в организации транкинговой радиосвязи и во взаимном расположении группировки АПП и подавляемого РПУ показывает, что для РЭЗ возможно использование пространственных различий между углами прихода сигнала и помех, но только по углу возвышения θ, т.е. в вертикальной плоскости (ВП). Следовательно, для решения задачи пространственной селекции помех, создаваемых группировкой АПП, необходимо наличие априорной информации об угловом распределении помехового поля (угловом спектре мощности помех) в вертикальной плоскости в раскрыве приемной антенны.

Научная задача, решаемая в статье, заключается в разработке математической модели пространственно-временной структуры (ПВС) помехового воздействия, создаваемого группировкой АПП, в вертикальной плоскости. Цель работы состоит в получении аналитической зависимости распределения интенсивности мощности помех sagd01.wmf на входе антенны подавляемого РПУ от угла возвышения q для конкретного момента времени t. Поскольку углы прихода помех являются случайными величинами, то модель должна обеспечивать получение статистических характеристик (математического ожидания М и СКО s), т.е.

sagd02.wmf (1)

где S – количество АПП в группировке, s – номер текущего АПП; sagd03.wmf – координаты s-го АПП в виде радиус-вектора и угла возвышения, sagd04.wmf – мощность помехи, создаваемая s-м АПП, на входе подавляемого РПУ; U – совокупность условий, влияющих на распространение помеховых радиоволн.

Отдельно следует остановиться на особенностях распространения помеховых радиоволн. Применяющиеся на практике транкинговые системы существенно различаются диапазоном рабочих длин волн (l от 0,345 до 2,174 м), поэтому на распространение помеховых радиоволн по-разному будет влиять земная поверхность, т.к. для разных l будет свой тангенс угла потерь в среде распространения [5]. Следовательно, эти обстоятельства должны быть учтены в разрабатываемой математической модели.

Решение научной задачи

Очевидно, что строгое решение волновых дифференциальных уравнений поля помех с учетом граничных эффектов раздела двух сред, а также различных свойств среды распространения и априорной неопределенности координат АПП чрезвычайно затруднено. В этом случае вероятностные характеристики создаваемого помехового воздействия можно получить на основе сочетания статистического подхода к пространственно-временному распределению (ПВР) углов прихода помех, создаваемых АПП в ВП [6, 7], с одной стороны, и феноменологического подхода к задаче определения ПВР помехового поля в точки подавления [8, 9], с другой стороны. Исходя из этого, требуемая модель ПВС разрабатывалась в два этапа:

1) статистическое моделирование ПВР углов прихода помех в ВП;

2) математическое моделирование ПВС помехового воздействия в ВП.

Исходными данными для моделирования выступают технические характеристики современных маломощных ПП, а также тактика группового применения АПП [1–4]. В основу моделирования помехового воздействия положена модель точечных источников влияния [7], которая конкретизирована применительно к АПП. Каждый АПП характеризуется функцией воздействия, имеющей физический смысл мощности помех, приведенной к полосе пропускания РПУ, с учетом характеристики направленности антенны ПП, т.е.

sagd05.wmf, (1')

где sagd06.wmf – мощность передатчика помех; sagd07.wmf – ширина спектра заградительной помехи; sagd08.wmf – полоса пропускания РПУ; sagd09.wmf – диаграмма направленности антенны ПП.

В результате каждый АПП создает в точке расположения подавляемого РПУ поле воздействия, определяемое как произведение функции воздействия ρ(θ) на передаточную характеристику sagd10.wmf среды распространения.

В первой статистической части произведено моделирование стохастической взаимосвязи подавляемого РПУ и группировки АПП в виде эквидистантной решетки. Для этого рассмотрим некоторую область U конечного объема в 4-мерном евклидовом пространстве R4, которая является областью физического пространства и времени. В этой области рассеяна в виде эквидистантной решетки счетная система из S АПП, расположенных в точках sagd11.wmf (s = 1, S). Очевидно, что максимальное РЭП наступит в момент времени sagd12.wmf, когда РПУ оказывается в объеме U0, который ограничен расположением четверки ближайших АПП, находящихся в центре решетки. При этом расстояние LT до траверса ближайшего АПП в данный момент времени является случайной величиной, распределенной по равномерному закону с плотностью sagd13.wmf, где Lи – интервал между ПП в ряду решетки, а Lд – дистанция между рядами.

Разбив область U на четыре непересекающиеся области (подрешетки) U1, U2, U3 и U4 относительно центра решетки, с учетом детерминированности взаимного расположения АПП, для конкретной z-й реализации случайных координат sagd14.wmf, sagd15.wmf и sagd16.wmf ближайшего АПП относительно местоположения РПУ, получим выражение (2), определяющее расстояние до траверса s-ro АПП в i-й подрешетке, т.е.

sagd17.wmf (2)

где сi и bi (i = 1,4) – соответственно номер ряда в i-й подрешетке и номер АПП в ряду; sagd18.wmf – величина погрешности горизонтальной координаты АПП, обусловленная неточностью их наведения и турбулентностью воздушных потоков; Z – объем выборки реализаций случайного взаимного расположения РПУ и решетки АПП. Величины Vи, Vд представляют собой проекции вектора скорости перемещения АПП в горизонтальной плоскости, они введены в (2) для обеспечения адекватности модели на случай, когда для подавления применяется группировка парашютируемых ПП, либо когда АПП барражируют в районе подавления.

В свою очередь угол прихода помех, определяемый как угол возвышения s-го АПП над горизонтом, равен

sagd19.wmf, (3)

где Нзад – заданная высота полета АПП, sagd20.wmf – величина неопределенности вертикальной координаты АПП, а Vв – проекция вектора скорости перемещения АПП на вертикальную ось.

Таким образом, для получения функции распределений sagd21.wmf углов прихода помех в конкретный момент времени необходимо методом статистических испытаний разыграть различные варианты взаимного расположении решетки АПП относительно подавляемого РПУ. По результатам накопленной статистики для диапазона углов возвышения sagd22.wmf дискретно с шагом Δθ определяется гистограмма распределения углов прихода помех от интервала углов возвышения. При этом среднее значение (математическое ожидание) и СКО функции распределения углов прихода помех определяются следующими выражениями:

sagd23.wmf (4)

sagd24.wmf (5)

В результате моделирования получена асимптотически не смещенная оценка среднего количества радиолучей помехового воздействия от углов прихода, попадающих в интервалы углов возвышения. Таким образом, полученные выражения (4) и (5) представляют собой статистическую модель ПВР углов прихода помех, создаваемых группировкой АПП.

На 2-м этапе моделирования произведено математическое описание структуры помехового поля в ВП с учетом явления интерференции помеховых радиоволн за счет влияния земной поверхности. Поскольку ПВР углов прихода помех в ВП в точке подавления имеет нестационарный и случайный характер, то среда распространения может быть представлена в виде линейного пространственного фильтра (ПФ), у которого передаточная характеристика (ПХ) имеет нестационарный характер, так как зависит от текущего времени воздействия.

Расстояния от АПП до точки подавления (раскрыва антенны РПУ) для прямой sagd25.wmf и отраженной sagd26.wmf радиоволны связаны с углом возвышения θ соответственно следующими соотношениями:

sagd27.wmf,

sagd28.wmf, (6)

где H – начальная высота полета АПП, h – высота подъема приемной антенны РПУ. Тогда для прямой волны на расстоянии sagd29.wmf мощность заградительной помехи на входе РПУ, в пределах его полосы пропускания (при условии, что ΔFЗП >> Δfпр) будет равна

sagd30.wmf, (7)

где sagd31.wmf – диаграмма направленности антенны РПУ; vП – коэффициент, учитывающий различия поляризации помехи и сигнала; l – рабочая длина волны радиосигнала.

Поскольку антенна подавляемого РПУ и АПП находятся высоко над земной поверхностью, то с учетом влияния подстилающей поверхности интерференционный множитель ослабления определяется формулой

sagd32.wmf, (8)

где sagd33.wmf и sagd34.wmf – соответственно модуль и фаза комплексного коэффициента отражения (Френеля) в ВП; sagd35.wmf – разность хода лучей прямых и отраженных помеховых радиоволн.

В результате ПХ среды распространения в ВП, как ПФ, будет определяться формулой (9), показывающей, как ослабляется энергия помехи в зависимости от угла ее прихода в ВП:

sagd36.wmf. (9)

Полученные в результате моделирования формализованные представления углов прихода помех в виде функции распределения sagd37.wmf и среды распространения в виде передаточной характеристики sagd38.wmf позволяют описать помеховое воздействие в точке приема (на входе антенны) в виде среднестатистического углового спектра мощности помех (УСМП) в ВП в конкретный момент времени t. При этом математическое выражение, описывающее среднее значение УСМП в ВП, в общем виде имеет вид

sagd39.wmf, [дБ], sagd40.wmf, (10)

где Т – время помехового воздействия АПП.

Выражение (10) представляет собой математическую модель, которая в абстрактной форме описывает ПВС помехового воздействия, создаваемого группировкой АПП в виде эквидистантной решетки. Физический смысл УСМП заключается в том, что sagd41.wmf демонстрирует, как в точке подавления распределена интенсивность мощности помех от угла возвышения (прихода помех) в конкретный момент помехового воздействия.

Результаты исследования и их обсуждение

С помощью разработанной модели с доверительной вероятностью 0,95 получено ПВР углов прихода помех в ВП, представляющее собой сглаженные статистические гистограммы. Моделирование производилось для случая массированного подавления РПУ базовой станции решеткой из 64-х АПП, равномерно распределенных в пространстве с интервалом и дистанцией в 500 м и включившихся на излучение на высоте 500 м и плавно снижающихся со скоростью 1–2 м/с. Полученные зависимости среднего значения sagd42.wmf и СКО sagd43.wmf представлены на рис. 1, а, б.

В результате исследований установлено, что ПВР углов прихода помех для конкретного момента времени сводится к устойчивому распределению. В фиксированный момент времени в ПВР имеется три характерных сектора углов возвышения:

1) сектор, в котором помехи отсутствуют, если sagd44.wmf;

2) детерминированный сектор, в котором помехи приходят всегда, если sagd45.wmf;

3) сектор, в котором помехи появляются спорадически, если sagd46.wmf.

Данные результаты следует использовать для пространственной селекции помех, создаваемых группировкой АПП.

Из анализа графиков, приведенных на рис. 2, а, б, следует, что ПХ среды распространения в ВП имеет сложный осциллирующий характер, который уменьшается при приближении угла возвышения к углу Брюстера (полного преломления радиоволн). При этом коэффициент передачи плавно возрастает с увеличением угла возвышения, поскольку при этом уменьшается расстояние от точки подавления до места расположения АПП.

Графики УСМП в точке подавления в начальный момент времени, а также через 120 секунд, приведены на рис. 3, а, б. Утолщенными линиями обозначены средние значения УСМП, а тонкими – его максимальные значения в пределах полученной дисперсии. Исследование УСМП показало, что его детерминированная часть вносит основной вклад в интегральную мощность помех и характеризуется наличием характерных экстремумов. Кроме того установлено, что неустойчивость УСМП течение времени t проявляется в сужении сектора, в котором помехи отсутствуют, в снижении углов характерных экстремумов детерминированного сектора и в увеличении спорадического сектора. Полученная в результате моделирования ПВС помехового воздействия в виде УСМП в ВП позволяет сформулировать требования к диаграмме направленности антенной системы подавляемого РПУ.

sagd1.tif

Рис. 1. Пространственно-временное распределение углов прихода помех в ВП

sagd2.tif

Рис. 2. Передаточная характеристика среды распространения в ВП

sagd3.tif

Рис. 3. Угловой спектр мощности помех в ВП

Заключение

В статье разработана математическая модель помехового воздействия в ВП, создаваемого группировкой АПП в виде эквидистантной решетки. Данная модель представляет собой математическое описание ПВС помехи на входе антенны подавляемого РПУ в виде среднестатистического УСМП в ВП для конкретного момента времени. Полученные в результате моделирования априорные данные об УСМП в ВП следует использовать для пространственной селекции помех, создаваемых АПП.