Радиотехника и электроника, 2020, T. 65, № 11, стр. 1091-1100

Метод наведения на воздушную цель в составе плотной группы

В. С. Верба a, И. Р. Загребельный ab*, В. И. Меркулов a

a АО Концерн радиостроения “Вега”
121170 Москва, Кутузовский просп., 34, Российская Федерация

b Московский государственный технический университет им. Н.Э. Баумана
105005 Москва, ул. 2-я Бауманская, 5, Российская Федерация

* E-mail: mail@vega.su

Поступила в редакцию 27.04.2020
После доработки 27.04.2020
Принята к публикации 11.06.2020

Полный текст (PDF)

Аннотация

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

ВВЕДЕНИЕ

В рамках все возрастающей роли группового противоборства в воздушно-космической сфере [1] как средств нападения, так и средств защиты [2], все более востребованной становится задача уничтожения отдельных летательных аппаратов (приоритетных целей [3]) в составе плотной группы.

Необходимо отметить, что к настоящему времени [4] подавляющее число средств поражения “воздух–воздух” ориентировано на поражение одиночных (разрешаемых) целей.

Если цели движутся в составе плотной группы и не разрешаются бортовой радиолокационной системой (РЛС) перехватчика ни по дальности, ни по скорости, ни по угловым координатам, то средства поражения будут наводиться на геометрический либо энергетический центр группы, не обеспечивая поражения интересующей цели [4].

В связи с этим весьма актуальной становится задача разработки методов индивидуального наведения на отдельные цели в составе плотной группы. В свою очередь эта задача разбивается на две подзадачи:

1) обеспечение разрешения бортовой РЛС до значений, позволяющих выделять отдельные цели на фоне других;

2) реализация высокоточного наведения на выбранную цель с минимизацией текущего промаха.

Среди различных видов разрешения наиболее сложно обеспечить разрешающую способность по углу ввиду естественных ограничений, определяемых шириной диаграммы направленности $\theta $ антенны бортовой РЛС носителя, которая зависит от соотношения размеров антенны $d$ и длины волны $\lambda $ измеряемого сигнала: $\theta \approx {\lambda \mathord{\left/ {\vphantom {\lambda d}} \right. \kern-0em} d}$ [5].

Перспективным направлением улучшения разрешающей способности по углу в однопозиционных бортовых РЛС является применение траекторного управления наблюдением [6], суть которого состоит в переходе от типового разрешения по углу, определяемого шириной диаграммы направленности, к разрешению по доплеровским частотам на основе эффекта доплеровского обострения луча (ДОЛ) [7], реализуемого в процессе полета носителя бортовой РЛС по специальной траектории.

Следует подчеркнуть, что для реализации ДОЛ полет летательного аппарата (ЛА) должен выполняться под некоторым вполне определенным углом относительно направления на цель, то время как для ее перехвата необходимо лететь либо на нее, либо в упрежденную точку встречи [8]. Таким образом, траектория должна удовлетворять двум противоречивым требованиям: обеспечивать разрешение цели в группе за счет ДОЛ и в то же время обеспечивать минимальный промах в процессе наведения на нее. То есть фактически необходимо решать двухкритериальную задачу.

Следует заметить, что особенностью этой задачи является возможность разнесения решаемых подзадач во времени, т.е. вначале обеспечить разрешение целей по угловым координатам, а на конечном участке – минимизировать промах. Однако для исполнения такого подхода необходимо иметь нестационарный закон управления [9], обеспечивающий автоматическое перераспределение функций управления от решения первой задачи ко второй.

В общем случае эта задача может быть решена различными способами [1012]. Однако их общим недостатком является использование как для разрешения, так и минимизации промаха в качестве управляемой координаты угловой скорости линии визирования, что с одной стороны ограничивает дальность их действия, а с другой – усложняет процедуру минимизации промаха.

Цель данной работы – синтез более совершенного метода наведения летательными аппаратами, обеспечивающего индивидуальное наведение на цель в составе плотной группы на значительно больших расстояниях, с автоматическим перераспределением приоритетов управления от решения задачи требуемого разрешения к решению задачи поражения цели.

При этом под методом наведения понимается закон формирования требуемой траектории, полет по которой позволит решить поставленную задачу [8].

1. ПОСТАНОВКА ЗАДАЧИ

Задача будет решаться на основе математического аппарата локального варианта статистической теории оптимального управления. Этот аппарат позволяет [13] для системы

(1)
${{{\mathbf{\dot {x}}}}_{{\text{у}}}}(t) = {{{\mathbf{F}}}_{{\text{у}}}}{{{\mathbf{x}}}_{{\text{у}}}}(t) + {{{\mathbf{B}}}_{{\text{у}}}}{\mathbf{u}}(t) + {{{\mathbf{\xi }}}_{{\text{у}}}}(t),$
предназначенной для отработки процесса

(2)
${{{\mathbf{\dot {x}}}}_{{\text{т}}}}(t) = {{{\mathbf{F}}}_{{\text{т}}}}{{{\mathbf{x}}}_{{\text{т}}}}(t)$

при наличии измерений

(3)
${\mathbf{z}}(t) = {\mathbf{Hx}}(t) + {{{\mathbf{\xi }}}_{{\text{и}}}}(t)$

сформировать сигнал управления

(4)
${\mathbf{u}}(t) = {{{\mathbf{K}}}^{{ - 1}}}{\mathbf{B}}_{{\text{у}}}^{T}{\mathbf{Q}}\left[ {{{{{\mathbf{\hat {x}}}}}_{{\text{т}}}}(t) - {{{{\mathbf{\hat {x}}}}}_{{\text{у}}}}(t)} \right],$

оптимальный по минимуму функционала качества

(5)
$\begin{gathered} I = {{M}_{{\text{у}}}}\left\{ {{{{\left[ {{{{\mathbf{x}}}_{{\text{т}}}}(t) - {{{\mathbf{x}}}_{{\text{у}}}}(t)} \right]}}^{T}}{\mathbf{Q}}\left[ {{{{\mathbf{x}}}_{{\text{т}}}}(t) - {{{\mathbf{x}}}_{{\text{у}}}}(t)} \right] + } \right. \\ \left. { + \,\,\int\limits_0^t {{\mathbf{u}}{{{(t)}}^{T}}{\mathbf{Ku}}(t)dt} } \right\}. \\ \end{gathered} $

В (1)–(5): ${{{\mathbf{x}}}_{{\text{т}}}}$ и ${{{\mathbf{x}}}_{{\text{у}}}}$$n$‑мерные векторы требуемых и управляемых координат состояния в текущие моменты времени $t;$ ${{{\mathbf{F}}}_{{\text{у}}}}$ и ${{{\mathbf{F}}}_{{\text{т}}}}$ – матрицы, определяющие динамические свойства подсистем (1) и (2); ${\mathbf{u}}$$r$‑мерный ($r \leqslant n$) вектор управления; ${{{\mathbf{B}}}_{{\text{у}}}}$ – матрица эффективности управления; ${\mathbf{z}}$$m$‑мерный ($m \leqslant 2n$) вектор измерений; ${\mathbf{H}}$ – матрица связи ${\mathbf{x}}$ и ${\mathbf{z}};$ ${{{\mathbf{\xi }}}_{{\text{у}}}},$ ${{{\mathbf{\xi }}}_{{\text{и}}}}$ – гауссовские центрированные векторы шумов состояния и измерений с известными матрицами спектральных плотностей; ${\mathbf{Q}}$ – неотрицательно определенная матрица штрафов за точность функционирования, ${\mathbf{K}}$ – положительно определенная матрица штрафов за величину сигналов управления; ${{{\mathbf{\hat {x}}}}_{{\text{т}}}}$ и ${{{\mathbf{\hat {x}}}}_{{\text{у}}}}$ – векторы оптимальных оценок процессов ${{{\mathbf{x}}}_{{\text{т}}}}$ и ${{{\mathbf{x}}}_{{\text{у}}}};$ My{·} –. операция математического ожидания при условии, что имеются измерения (3), Т – знак транспонирования

Отличительной особенностью закона управления (4), кроме простоты, является то, что он не требует априорного знания временного интервала функционирования, обеспечивая оптимизацию системы наведения на каждый момент времени. Кроме того, непосредственно для его реализации не требуется знания матриц ${{{\mathbf{F}}}_{{\text{у}}}}$ и ${{{\mathbf{F}}}_{{\text{т}}}}$ внутренних связей в (1) и (2), а достаточно знания оценок ${{{\mathbf{\hat {x}}}}_{{\text{т}}}}$ и ${{{\mathbf{\hat {x}}}}_{{\text{у}}}},$ полученных тем или иным способом.

Для реализации этой целевой установки необходимо найти:

а) закон формирования ${{{\mathbf{x}}}_{{\text{т}}}}$ в (4), (5), обеспечивающий разрешение целей по угловым координатам;

б) механизм перераспределения приоритетов управления от разрешения к устранению промаха;

в) закон управления, минимизирующий промахи на конечном участке.

Необходимо отметить, что, в общем случае, наведение необходимо осуществлять как в горизонтальной, так и в вертикальной плоскостях. Ниже для простоты рассмотрен синтез закона управления лишь для горизонтальной плоскости, в предположении, что закон управления для вертикальной плоскости будет аналогичным.

2. ОБОСНОВАНИЕ ТРЕБУЕМОЙ ТРАЕКТОРИИ ПОЛЕТА НОСИТЕЛЯ БОРТОВОЙ РЛС

Рассмотрим зависимость линейного разрешения бортовой РЛС по азимуту от взаимного расположения самолета-носителя и двух, близкорасположенных целей, не разрешаемых по углу и движущихся с одинаковыми скоростями и в одном направлении.

Пусть на носителе, находящемся в точке ${{O}_{{\text{н}}}}$ (рис. 1) на удалении $D$ от цели ${{O}_{{{\text{ц1}}}}}$ и близкорасположенной к ней цели ${{O}_{{{\text{ц2}}}}},$ движущемся со скоростью ${{V}_{{\text{н}}}}$ под углом ${{\varphi }_{{\text{н}}}}$ к ${{O}_{{{\text{ц1}}}}},$ формируется сигнал подсвета с длиной волны $\lambda $, который после отражения от цели принимается приемником бортовой РЛС. На рисунке: $\Delta \varphi < \theta $ – угол между целями ${{O}_{{{\text{ц1}}}}}$ и ${{O}_{{{\text{ц2}}}}},$ при котором они не разрешаются по углу; ${{\varphi }_{{\text{ц}}}}$ – угол между направлением скорости цели и линией визирования; $\omega $ – угловая скорость вращения линии визирования цели ${{O}_{{{\text{ц1}}}}};$ ${{j}_{{\text{н}}}}$ – поперечное ускорение носителя; $\theta $ – ширина диаграммы направленности бортовой РЛС. При этом минимальное расстояние $\Delta {{l}_{{\text{т}}}}$ между целями, при котором принимается решение о наличии двух целей, именуемое разрешающей способностью, существенно (более, чем на порядок) меньше дальности до цели:

(6)
$\Delta {{l}_{т}} \ll D.$
Рис. 1.

Геометрия взаимного расположения носителя бортовой РЛС ${{O}_{{\text{н}}}}$ и целей ${{O}_{{{\text{ц1}}}}},$ ${{O}_{{{\text{ц2}}}}}.$

Тогда доплеровская частота принимаемого от цели ${{O}_{{{\text{ц1}}}}}$ сигнала

(7)
${{F}_{{d1}}} = \frac{{2{{V}_{{\text{н}}}}\cos {{\varphi }_{{\text{н}}}}}}{\lambda } + \frac{{2{{V}_{{\text{ц}}}}\cos {{\varphi }_{{\text{ц}}}}}}{\lambda }.$

В свою очередь, доплеровская частота сигнала, принимаемого от цели ${{O}_{{{\text{ц2}}}}},$ удаленной от цели ${{O}_{{{\text{ц1}}}}}$ на расстояние $\Delta {{l}_{т}} \ll D,$ определяется как

${{F}_{{d2}}} = \frac{{2{{V}_{{\text{н}}}}\cos \left( {{{\varphi }_{{\text{н}}}} + \Delta \varphi } \right)}}{\lambda } + \frac{{2{{V}_{{\text{ц}}}}\cos \left( {{{\varphi }_{{\text{ц}}}} - \Delta \varphi } \right)}}{\lambda }.$

Разложив это выражение в ряд Тейлора с линейным приближением, получим:

(8)
$\begin{gathered} {{F}_{{d2}}} = \frac{2}{\lambda }\left[ {{{V}_{{\text{н}}}}\cos {{\varphi }_{{\text{н}}}} - {{V}_{{\text{н}}}}\left( {\sin {{\varphi }_{{\text{н}}}}} \right)\Delta \varphi + } \right. \\ \left. { + \,\,{{V}_{{\text{ц}}}}\cos {{\varphi }_{{\text{ц}}}} + {{V}_{{\text{ц}}}}\left( {\sin {{\varphi }_{{\text{ц}}}}} \right)\Delta \varphi } \right] = \\ = {{F}_{{d1}}} + \frac{2}{\lambda }\left( {{{V}_{{\text{ц}}}}\sin {{\varphi }_{{\text{ц}}}} - {{V}_{{\text{н}}}}\sin {{\varphi }_{{\text{н}}}}} \right)\Delta \varphi = \\ = {{F}_{{d1}}} + \frac{{2\Delta {{l}_{{\text{т}}}}}}{{\lambda D}}\left( {{{V}_{{\text{ц}}}}\sin {{\varphi }_{{\text{ц}}}} - {{V}_{{\text{н}}}}\sin {{\varphi }_{{\text{н}}}}} \right), \\ \end{gathered} $

где было учтено равенство (7) и соотношение $\Delta \varphi \approx {{\Delta {{l}_{{\text{т}}}}} \mathord{\left/ {\vphantom {{\Delta {{l}_{{\text{т}}}}} D}} \right. \kern-0em} D}$ для малых углов.

Чтобы разрешить цели ${{O}_{{{\text{ц1}}}}}$ и ${{O}_{{{\text{ц2}}}}}$ по доплеровской частоте, полоса пропускания доплеровского фильтра бортовой РЛС не должна превышать величины

(9)
$\Delta F = {{F}_{{d1}}} - {{F}_{{d2}}} = \frac{{2\Delta {{l}_{{\text{т}}}}}}{{\lambda D}}\left( {{{V}_{{\text{н}}}}\sin {{\varphi }_{{\text{н}}}} - {{V}_{{\text{ц}}}}\sin {{\varphi }_{{\text{ц}}}}} \right).$

Если полоса доплеровского фильтра $\Delta F$ и длина волны $\lambda $ известны, то можно предъявить требования к бортовому пеленгу цели ${{\varphi }_{{\text{н}}}},$ при котором цели будут разрешаться, находясь друг от друга на расстоянии $\Delta {{l}_{{\text{т}}}}.$ Из (9) можно получить:

(10)
${{V}_{{\text{н}}}}\sin {{\varphi }_{{\text{н}}}} - {{V}_{{\text{ц}}}}\sin {{\varphi }_{{\text{ц}}}} = \frac{{\Delta F\lambda D}}{{2\Delta {{l}_{{\text{т}}}}}},$
(11)
$\sin {{\varphi }_{{\text{н}}}} = \frac{{\Delta F\lambda D}}{{2\Delta {{l}_{{\text{т}}}}{{V}_{{\text{н}}}}}} + \frac{{{{V}_{{\text{ц}}}}}}{{{{V}_{{\text{н}}}}}}\sin {{\varphi }_{{\text{ц}}}}.$

Из (11) следует, что для реализации требуемого линейного разрешения $\Delta {{l}_{{\text{т}}}}$ траектория полета носителя должна быть такой, чтобы требуемый бортовой пеленг цели с носителя обладал значением не меньшим, чем

(12)
${{\varphi }_{{{\text{нт}}}}} = \arcsin \left[ {\frac{{\Delta F\lambda D}}{{2\Delta {{l}_{{\text{т}}}}{{V}_{{\text{н}}}}}} + \frac{{{{V}_{{\text{ц}}}}}}{{{{V}_{{\text{н}}}}}}\sin {{\varphi }_{{\text{ц}}}}} \right].$

При этом, учитывая, что ${{\left( {{{V}_{{\text{н}}}}\sin {{\varphi }_{{\text{н}}}} - {{V}_{{\text{ц}}}}\sin {{\varphi }_{{\text{ц}}}}} \right)} \mathord{\left/ {\vphantom {{\left( {{{V}_{{\text{н}}}}\sin {{\varphi }_{{\text{н}}}} - {{V}_{{\text{ц}}}}\sin {{\varphi }_{{\text{ц}}}}} \right)} D}} \right. \kern-0em} D} = $ ω, из (10) получаем требования к угловой скорости линии визирования:

(13)
${{\omega }_{{\text{т}}}} = \frac{{\Delta F\lambda }}{{2\Delta {{l}_{{\text{т}}}}}}.$

Анализ (12) и (13) позволяет заключить, что требуемые значения бортового пеленга ${{\varphi }_{{{\text{нт}}}}}$ и угловой скорости ${{\omega }_{{\text{т}}}},$ необходимые для достижения требуемого линейного разрешения $\Delta {{l}_{{\text{т}}}},$ зависят как от параметров системы обработки ($\lambda ,$ $\Delta F$), так и от условий применения ($D,$ ${{V}_{{\text{н}}}},$ ${{V}_{{\text{ц}}}}$ и ${{\varphi }_{{\text{ц}}}}$).

Необходимо подчеркнуть, что (12) и (13) имеют естественное ограничение по дальности применения, обусловленное, с одной стороны, вытекающим из (11) требованием $\sin {{\varphi }_{{{\text{нт}}}}} \leqslant 1,$ а с другой стороны сложностью реализации требуемых значений (13) на больших дальностях в реальном диапазоне скоростей движения целей и носителя. При этом последнее ограничение более жесткое.

3. СИНТЕЗ ЗАКОНА УПРАВЛЕНИЯ, ОБЕСПЕЧИВАЮЩЕГО НАВЕДЕНИЕ НА ЦЕЛЬ В СОСТАВЕ ПЛОТНОЙ ГРУППЫ

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

Требуемое линейное разрешение $\Delta {{l}_{{\text{т}}}},$ определяемое интервалом между целями в группе, можно обеспечить соответствующим выбором угла ${{\varphi }_{{\text{н}}}} = {{\varphi }_{{{\text{нт}}}}}$ и угловой скорости линии визирования $\omega = {{\omega }_{{\text{т}}}}.$ В свою очередь, точность наведения на цель характеризуется величиной текущего промаха [8]

(14)
$h \approx \frac{{{{D}^{2}}\omega }}{{{{V}_{{{\text{сб}}}}}}},$

где ${{V}_{{{\text{сб}}}}} = - \dot {D}$ – скорость сближения объекта управления с целью, и для его уменьшения необходимо минимизировать величину $\omega .$

Спецификой модели состояния системы носитель – плотная группа является использование эквивалента носитель – одиночная цель, в котором роль одиночной цели выполняет элемент, координатами которого являются расстояние до центра группы и средний бортовой пеленг, соответствующий центру группы. Эти условия обычно выполняются на практике, если группа пространственно укладывается в элемент разрешения (рис. 2), а расстояния между участниками группы удовлетворяют соотношениям:

(15)
$\Delta {{l}_{{\text{д}}}} \ll {{D}_{{{\text{ср}}}}},\,\,\,\,\Delta {{l}_{{{\text{аз}}}}} \ll {{D}_{{{\text{ср}}}}},$
Рис. 2.

Определение дальности ${{D}_{{{\text{ср}}}}}$ и бортового пеленга ${{\varphi }_{{{\text{ср}}}}}$ центра группы.

где ${{D}_{{{\text{ср}}}}}$ – среднее расстояние до центра группы, $\Delta {{l}_{{\text{д}}}}$ – линейное разрешение по дальности, $\Delta {{l}_{{{\text{аз}}}}}$ – линейное разрешение по азимуту.

Ввиду выполнения условия (15) правомочно использование модели состояния [8]:

(16)
$\begin{gathered} {{{\dot {\varphi }}}_{{\text{н}}}} = \omega ,\,\,\,\,{{\varphi }_{{\text{н}}}}(0) = {{\varphi }_{{{\text{н0}}}}}; \\ \dot {\omega } = - \frac{{2\dot {D}}}{D}\omega - \frac{1}{D}{{j}_{{\text{н}}}} + {{\xi }_{{\omega }}},\,\,\,\,\omega (0) = {{\omega }_{{\text{0}}}}, \\ \end{gathered} $

которая получена с учетом того, что Δlт$ \ll $ D, $ \ll $ а цели в плотной группе не маневрируют. Достоинством модели (16), в которой ${{\xi }_{{\omega }}}$ – шумы состояния, является ее адаптируемость к условиям применения, обусловленная учетом маневра носителя за счет поперечного ускорения ${{j}_{{\text{н}}}},$ и зависимость $\omega $ от дальности и скорости сближения. Следует подчеркнуть, что средние значения ${{\varphi }_{{\text{н}}}},$ $\omega ,$ ${{j}_{{\text{н}}}},$ $D$ и $\dot {D}$ достаточно просто оцениваются существующими датчиками.

Синтез закона управления, совместно наилучшего по разрешающей способности, точности наведения самолета и экономичности будем выполнять при условии, что соблюдаются следующие допущения:

1) известны интервалы между целями $\Delta {{l}_{{\text{т}}}},$ определяющие требуемую линейную разрешающую способность по азимуту;

2) цели движутся в произвольном направлении с известной постоянной скоростью;

3) известен диапазон располагаемых перегрузок и допустимый промах, величина которого определяется согласно (14);

4) на борту наводимого объекта в рамках первичной обработки в РЛС проводится спектральный анализ принимаемых сигналов, а в рамках вторичной формируются измерения дальности до целей, скорости сближения с ними и бортовых пеленгов [14, 15];

5) канал наведения в вертикальной плоскости не влияет на канал наведения в горизонтальной плоскости.

Функционал качества (5), характеризующий эффективность закона управления с разрешением целей в плотной группе, должен учитывать ошибки по углу ${{\varphi }_{{{\text{нт}}}}} - {{\varphi }_{{\text{н}}}}$ и угловой скорости ${{\omega }_{{\text{т}}}} - \omega ,$ используемые для обеспечения требуемого линейного разрешения, ошибки по угловой скорости ${{\omega }_{h}} - \omega ,$ по которым можно минимизировать величину промаха (14), и расход энергии на сигналы управления ${{j}_{{\text{н}}}}$ в (16). Здесь ${{\omega }_{h}} = 0$ – требуемая для минимизации промаха величина угловой скорости линии визирования. С учетом этих особенностей минимизируемый функционал качества (5) может быть представлен в следующем виде:

(17)
$\begin{gathered} I = M\left\{ {{{{\left[ {\begin{array}{*{20}{c}} {{{\varphi }_{{{\text{нт}}}}} - {{\varphi }_{{\text{н}}}}} \\ {{{\omega }_{{\text{т}}}} - \omega } \\ {{{\omega }_{h}} - \omega } \end{array}} \right]}}^{Т}}\left[ {\begin{array}{*{20}{c}} {{{q}_{{11}}}}&{{{q}_{{12}}}}&{{{q}_{{13}}}} \\ {{{{{q}_{{21}}}D} \mathord{\left/ {\vphantom {{{{q}_{{21}}}D} {{{V}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{сб}}}}}}}}&{{{{{q}_{{22}}}D} \mathord{\left/ {\vphantom {{{{q}_{{22}}}D} {{{V}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{сб}}}}}}}}&{{{q}_{{23}}}} \\ {{{{{q}_{{31}}}D} \mathord{\left/ {\vphantom {{{{q}_{{31}}}D} {{{V}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{сб}}}}}}}}&{{{{{q}_{{32}}}D} \mathord{\left/ {\vphantom {{{{q}_{{32}}}D} {{{V}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{сб}}}}}}}}&{{{q}_{{33}}}} \end{array}} \right]} \right. \times \\ \left. { \times \,\,\left[ {\begin{array}{*{20}{c}} {{{\varphi }_{{{\text{нт}}}}} - {{\varphi }_{{\text{н}}}}} \\ {{{\omega }_{{\text{т}}}} - \omega } \\ {{{\omega }_{h}} - \omega } \end{array}} \right] + \int\limits_0^t {j_{{\text{н}}}^{2}{{k}_{j}}dt} } \right\}, \\ \end{gathered} $

где ${{q}_{{11}}},$ …, ${{q}_{{33}}}$ – веса штрафов за точность выдерживания линейного разрешения по азимуту и величину промаха, ${{k}_{j}}$ – штраф за величину сигнала управления.

Сравнивая (16) с (1), а (17) с (5), получаем:

(18)
$\begin{gathered} {{{\mathbf{x}}}_{{\text{т}}}} = \left[ {\begin{array}{*{20}{c}} {{{\varphi }_{{{\text{нт}}}}}} \\ {{{\omega }_{{\text{т}}}}} \\ {{{\omega }_{h}}} \end{array}} \right],\,\,\,\,{{{\mathbf{x}}}_{{\text{у}}}} = \left[ {\begin{array}{*{20}{c}} {{{\varphi }_{{\text{н}}}}} \\ \omega \\ \omega \end{array}} \right],\,\,\,\,{{{\mathbf{F}}}_{{\text{у}}}} = \left[ {\begin{array}{*{20}{c}} 0&1&0 \\ 0&{{{ - 2\dot {D}} \mathord{\left/ {\vphantom {{ - 2\dot {D}} D}} \right. \kern-0em} D}}&0 \\ 0&0&{{{ - 2\dot {D}} \mathord{\left/ {\vphantom {{ - 2\dot {D}} D}} \right. \kern-0em} D}} \end{array}} \right], \\ {{{\mathbf{B}}}_{{\text{у}}}} = \left[ {\begin{array}{*{20}{c}} 0 \\ { - {1 \mathord{\left/ {\vphantom {1 D}} \right. \kern-0em} D}} \\ { - {1 \mathord{\left/ {\vphantom {1 D}} \right. \kern-0em} D}} \end{array}} \right],\,\,\,\,{\mathbf{Q}} = \left[ {\begin{array}{*{20}{c}} {{{q}_{{11}}}}&{{{q}_{{12}}}}&{{{q}_{{13}}}} \\ {{{{{q}_{{21}}}D} \mathord{\left/ {\vphantom {{{{q}_{{21}}}D} {{{V}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{сб}}}}}}}}&{{{{{q}_{{22}}}D} \mathord{\left/ {\vphantom {{{{q}_{{22}}}D} {{{V}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{сб}}}}}}}}&{{{q}_{{23}}}} \\ {{{{{q}_{{31}}}D} \mathord{\left/ {\vphantom {{{{q}_{{31}}}D} {{{V}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{сб}}}}}}}}&{{{{{q}_{{32}}}D} \mathord{\left/ {\vphantom {{{{q}_{{32}}}D} {{{V}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{сб}}}}}}}}&{{{q}_{{33}}}} \end{array}} \right], \\ {\mathbf{K}} = \left[ {{{k}_{j}}} \right],\,\,\,\,{\mathbf{u}} = \left[ {{{j}_{{\text{н}}}}} \right]. \\ \end{gathered} $

Спецификой матрицы штрафов ${\mathbf{Q}}$ в функционале (17) является ее нестационарный характер, обусловленный использованием сомножителей ${D \mathord{\left/ {\vphantom {D {{{V}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{сб}}}}}}},$ где ${{V}_{{{\text{сб}}}}} = - \dot {D}.$ Целесообразность использования нестационарных штрафов обусловлена требованием обеспечения автоматического перераспределения управления от формирования траектории, обеспечивающей разрешение целей в плотной группе к поражению цели на последнем участке.

Используя (18) в (4), имеем:

(19)
$\begin{gathered} {{j}_{{\text{н}}}} = \frac{1}{{{{k}_{j}}}}\left[ {\begin{array}{*{20}{c}} 0&{ - \frac{1}{{\hat {D}}}}&{ - \frac{1}{{\hat {D}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{{q}_{{11}}}}&{{{q}_{{12}}}}&{{{q}_{{13}}}} \\ {{{{{q}_{{21}}}\hat {D}} \mathord{\left/ {\vphantom {{{{q}_{{21}}}\hat {D}} {{{{\hat {V}}}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{{\hat {V}}}_{{{\text{сб}}}}}}}}&{{{{{q}_{{22}}}\hat {D}} \mathord{\left/ {\vphantom {{{{q}_{{22}}}\hat {D}} {{{{\hat {V}}}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{{\hat {V}}}_{{{\text{сб}}}}}}}}&{{{q}_{{23}}}} \\ {{{{{q}_{{31}}}\hat {D}} \mathord{\left/ {\vphantom {{{{q}_{{31}}}\hat {D}} {{{{\hat {V}}}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{{\hat {V}}}_{{{\text{сб}}}}}}}}&{{{{{q}_{{32}}}\hat {D}} \mathord{\left/ {\vphantom {{{{q}_{{32}}}\hat {D}} {{{{\hat {V}}}_{{{\text{сб}}}}}}}} \right. \kern-0em} {{{{\hat {V}}}_{{{\text{сб}}}}}}}}&{{{q}_{{33}}}} \end{array}} \right] \times \\ \times \,\,\left[ {\begin{array}{*{20}{c}} {{{{\hat {\varphi }}}_{{{\text{нт}}}}} - {{{\hat {\varphi }}}_{{\text{н}}}}} \\ {{{{\hat {\omega }}}_{{\text{т}}}} - \hat {\omega }} \\ {{{{\hat {\omega }}}_{h}} - \hat {\omega }} \end{array}} \right] = - \frac{{{{q}_{{21}}} + {{q}_{{31}}}}}{{{{k}_{j}}{{{\hat {V}}}_{{{\text{сб}}}}}}}\left( {{{{\hat {\varphi }}}_{{{\text{нт}}}}} - {{{\hat {\varphi }}}_{{\text{н}}}}} \right) - \\ - \,\,\frac{{{{q}_{{22}}} + {{q}_{{32}}}}}{{{{k}_{j}}{{{\hat {V}}}_{{{\text{сб}}}}}}}\left( {{{{\hat {\omega }}}_{{\text{т}}}} - \hat {\omega }} \right) - \frac{{{{q}_{{23}}} + {{q}_{{33}}}}}{{{{k}_{j}}\hat {D}}}\left( {{{{\hat {\omega }}}_{h}} - \hat {\omega }} \right). \\ \end{gathered} $

Принимая во внимание, что требуемое для минимизации промаха значение угловой скорости ${{\omega }_{h}} = 0,$ закон управления можно представить в виде:

(20)
$\begin{gathered} {{j}_{{\text{н}}}} = - \frac{{{{q}_{{21}}} + {{q}_{{31}}}}}{{{{k}_{j}}{{{\hat {V}}}_{{{\text{сб}}}}}}}\left( {{{{\hat {\varphi }}}_{{{\text{нт}}}}} - {{{\hat {\varphi }}}_{{\text{н}}}}} \right) - \\ - \,\,\frac{{{{q}_{{22}}} + {{q}_{{32}}}}}{{{{k}_{j}}{{{\hat {V}}}_{{{\text{сб}}}}}}}\left( {{{{\hat {\omega }}}_{{\text{т}}}} - \hat {\omega }} \right) + \frac{{{{q}_{{23}}} + {{q}_{{33}}}}}{{{{k}_{j}}\hat {D}}}\hat {\omega }. \\ \end{gathered} $

Анализ закона (20) позволяет прийти к следующим заключениям:

а) сигнал управления определяется ошибками, которые характеризуются несоответствием текущих значений ${{\varphi }_{{\text{н}}}},$ $\omega $ требуемым ${{\varphi }_{{{\text{нт}}}}},$ ${{\omega }_{{\text{т}}}}.$ Необходимо подчеркнуть, что на больших расстояниях, когда $\omega \approx 0$ и $\frac{{{{q}_{{23}}} + {{q}_{{33}}}}}{{{{k}_{j}}D}} \approx 0,$ основной вклад в управление вносят первые два слагаемые, обеспечивая требуемое линейное разрешение по азимуту. При этом первые два слагаемых образуют высокоточный пропорционально-дифференциальный метод последовательных упреждений [8], в котором второе слагаемое выполняет роль корректирующих сигналов, повышая точность отработки ${{\varphi }_{{{\text{нт}}}}}.$ По мере уменьшения $D$ увеличивается вклад третьего слагаемого и на небольших расстояниях превалирующей становится задача обеспечения минимального промаха;

б) веса ошибок управления зависят от соотношения штрафов ${{\left( {{{q}_{{21}}} + {{q}_{{31}}}} \right)} \mathord{\left/ {\vphantom {{\left( {{{q}_{{21}}} + {{q}_{{31}}}} \right)} {{{k}_{j}}}}} \right. \kern-0em} {{{k}_{j}}}},$ ${{\left( {{{q}_{{22}}} + {{q}_{{32}}}} \right)} \mathord{\left/ {\vphantom {{\left( {{{q}_{{22}}} + {{q}_{{32}}}} \right)} {{{k}_{j}}}}} \right. \kern-0em} {{{k}_{j}}}}$ и ${{\left( {{{q}_{{23}}} + {{q}_{{33}}}} \right)} \mathord{\left/ {\vphantom {{\left( {{{q}_{{23}}} + {{q}_{{33}}}} \right)} {{{k}_{j}}}}} \right. \kern-0em} {{{k}_{j}}}}.$ При этом, манипулируя их величинами можно управлять моментом перехода от управления, обеспечивающего стабилизацию линейного разрешения, к управлению, минимизирующему промах. Признаком разрешения сигналов по результатам спектрального анализа [14] является регистрация появления сигналов в других доплеровских фильтрах [11], позволяющих определить количественный состав группы и выделить цель, предназначенную для поражения;

в) для реализации метода наведения по закону (20) бортовая РЛС должна формировать типовой набор оценок дальности до цели, скорости сближения с ней, углов и угловой скорости линии визирования, варианты алгоритмов формирования которых рассмотрены в [14, 15].

4. ИССЛЕДОВАНИЕ ЗАКОНА УПРАВЛЕНИЯ, ОБЕСПЕЧИВАЮЩЕГО НАВЕДЕНИЕ НА ЦЕЛЬ В СОСТАВЕ ПЛОТНОЙ ГРУППЫ

Целью исследований являлась проверка возможности построения траектории на основе (20), полет по которой позволит разрешить цели в плотной группе с дальнейшим поражением одной из них. В связи с этим подробности спектрального анализа сигналов в бортовой РЛС далее не рассматривались.

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

1) скорость перехватчика составляет Vн = $ = 300\,\,{{\text{м}} \mathord{\left/ {\vphantom {{\text{м}} {\text{с}}}} \right. \kern-0em} {\text{с}}},$ а скорость целей, движущихся по прямолинейной траектории – ${{V}_{{\text{ц}}}} = 250\,\,{{\text{м}} \mathord{\left/ {\vphantom {{\text{м}} {\text{с}}}} \right. \kern-0em} {\text{с}}};$

2) шумы состояния в (16) отсутствуют и управление безынерционно;

3) в качестве показателей эффективности функционирования метода наведения на цель в составе плотной группы были приняты промах $h$ (14), текущее поперечное ускорение ${{j}_{{\text{н}}}}$ управляемого объекта и текущее линейное разрешение бортовой РЛС перехватчика по азимуту $\Delta l.$

При моделировании полагалось, что точка начала наведения перехватчика имеет координаты ${{z}_{{{\text{н0}}}}} = 20{\kern 1pt} {\kern 1pt} 000\,\,{\text{м,}}$ ${{x}_{{{\text{н0}}}}} = 10{\kern 1pt} {\kern 1pt} 000\,\,{\text{м,}}$ координаты целей ${{z}_{n}} = {{z}_{{n + 1}}}$ = ${{z}_{{n - 1}}} = 80{\kern 1pt} {\kern 1pt} 000\,\,{\text{м,}}$ ${{x}_{n}} = 60{\kern 1pt} {\kern 1pt} 000\,\,{\text{м,}}$ ${{x}_{{n - 1}}} = 59{\kern 1pt} {\kern 1pt} 750\,\,{\text{м,}}$ ${{x}_{{n + 1}}} = 60{\kern 1pt} {\kern 1pt} 250\,\,{\text{м,}}$ т.е. интервал между объектами в плотной группе составляет $\Delta {{l}_{{\text{т}}}} = 250\,\,{\text{м}}{\text{.}}$ При этом считалось, что на носителе используется типовая бортовая импульсно-доплеровская РЛС [14] с полосой пропускания доплеровских фильтров $\Delta F = 50\,\,{\text{Гц}}{\text{.}}$ Цель наведения – обеспечивая в процессе перехвата требуемое разрешение, поразить выбранную цель ${{{\text{Ц}}}_{n}}$ в плотной группе.

При этом проверялись три варианта управления – для обеспечения линейного разрешения $\Delta {{l}_{{{\text{т1}}}}} = 300\,\,{\text{м,}}$ $\Delta {{l}_{{{\text{т2}}}}} = 250\,\,{\text{м}}$ и $\Delta {{l}_{{{\text{т3}}}}} = 200\,\,{\text{м,}}$ которым на графиках соответствуют эпюры 1, 2 и 3. На рис. 3 приведены траектории участников наведения, полученные в результате эксперимента.

Рис. 3.

Траектории целей и носителя для трех вариантов требуемого линейного разрешения: кривая 1$\Delta {{l}_{{{\text{т1}}}}} = 300\,\,{\text{м,}}$ кривая 2$\Delta {{l}_{{{\text{т2}}}}} = 250\,\,{\text{м,}}$ кривая 3$\Delta {{l}_{{{\text{т3}}}}} = 200\,\,{\text{м}}{\text{.}}$

Результаты иллюстрируют успешное наведение на цель по плавной криволинейной траектории для всех трех вариантов $\Delta {{l}_{{\text{т}}}}.$ При этом, чем меньше требуемое линейное разрешение, тем большей кривизной обладает траектория. Маркерами выделены точки начала движения объектов и точки перехвата цели.

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

Рис. 4.

Зависимости текущих и требуемых (нумерация со штрихом) линейных разрешений по азимуту от времени для трех вариантов требуемого линейного разрешения: кривая 1$\Delta {{l}_{{{\text{т1}}}}} = 300\,\,{\text{м,}}$ кривая 2$\Delta {{l}_{{{\text{т2}}}}} = 250\,\,{\text{м,}}$ кривая 3$\Delta {{l}_{{{\text{т3}}}}} = 200\,\,{\text{м}}{\text{.}}$

На рис. 5 и 6 приведены графики дальностей до цели и промахов, рассчитываемых по формуле (14). Из рисунков видно, что дальности и промахи эффективно уменьшаются после начала наведения. При этом, чем более жесткие требования предъявляются к обеспечению заданного линейного разрешения (и, соответственно, чем больше кривизна траектории), тем большее время занимает процесс перехвата.

Рис. 5.

Зависимость дальностей до цели от времени для трех вариантов требуемого линейного разрешения: кривая 1$\Delta {{l}_{{{\text{т1}}}}} = 300\,\,{\text{м,}}$ кривая 2$\Delta {{l}_{{{\text{т2}}}}} = 250\,\,{\text{м,}}$ кривая 3$\Delta {{l}_{{{\text{т3}}}}} = 200\,\,{\text{м}}{\text{.}}$

Рис. 6.

Зависимость текущих промахов от времени для трех вариантов требуемого линейного разрешения: кривая 1$\Delta {{l}_{{{\text{т1}}}}} = 300\,\,{\text{м,}}$ кривая 2$\Delta {{l}_{{{\text{т2}}}}} = 250\,\,{\text{м,}}$ кривая 3$\Delta {{l}_{{{\text{т3}}}}} = 200\,\,{\text{м}}{\text{.}}$

Для оценки возможности практического применения закона управления (20) представляют интерес значения поперечных ускорений ${{j}_{{\text{н}}}}$ (перегрузки, испытываемые наводимым объектом) в процессе наведения. График зависимости поперечных ускорений от времени представлен на рис. 7.

Рис. 7.

Зависимость поперечного ускорения от времени для трех вариантов требуемого линейного разрешения: кривая 1$\Delta {{l}_{{{\text{т1}}}}} = 300\,\,{\text{м,}}$ кривая 2$\Delta {{l}_{{{\text{т2}}}}} = 250\,\,{\text{м,}}$ кривая 3$\Delta {{l}_{{{\text{т3}}}}} = 200\,\,{\text{м}}{\text{.}}$

Наибольшие перегрузки носитель испытывает в начальные моменты времени, что связано с выходом на требуемое направление движения. При этом минимальные значения разрешающей способности ограничиваются, с одной стороны, значением полосы пропускания $\Delta F$ доплеровских фильтров, а с другой стороны, допустимой перегрузкой летательного аппарата. Чем шире полоса пропускания фильтров, тем больше кривизна требуемой траектории и тем больше значения перегрузок, которые испытывает летательный аппарат.

После выхода на траекторию носитель обладает небольшими значениями поперечного ускорения вплоть до поражения цели.

Необходимо отметить, что в исследуемом законе (20) за стабилизацию линейного разрешения по азимуту отвечает два первых слагаемых, реализуя пропорционально-дифференциальный закон [8], в котором первое слагаемое играет основную роль, а второе – роль корректирующего сигнала. В связи с этим представляет интерес исследовать влияние каждого слагаемого на показатели эффективности наведения. Результаты исследования упрощенного варианта метода наведения, без второго корректирующего сигнала, показали, что несколько ухудшается точность выдерживания линейного разрешения (рис. 8), в то время как остальные показатели эффективности практически не изменяются.

Рис. 8.

Зависимости текущих и требуемых (нумерация со штрихом) линейных разрешений по азимуту от времени при законе управления без второго слагаемого для трех вариантов требуемого линейного разрешения: кривая 1$\Delta {{l}_{{{\text{т1}}}}} = 300\,\,{\text{м,}}$ кривая 2$\Delta {{l}_{{{\text{т2}}}}} = 250\,\,{\text{м,}}$ кривая 3$\Delta {{l}_{{{\text{т3}}}}} = 200\,\,{\text{м}}{\text{.}}$

В то же время исключение первого слагаемого в значительной мере влияет на все показатели эффективности. Графики результатов исследований этого варианта приведены на рис. 9–12, по которым видно, что особенно сильно, по сравнению с полным вариантом закона (20), изменяются разрешающие способности, промахи и перегрузки, представленные на рис. 10–12.

Рис. 9.

Траектории целей и носителя для трех вариантов требуемого линейного разрешения при законе управления без первого слагаемого: кривая 1$\Delta {{l}_{{{\text{т1}}}}} = 300\,\,{\text{м,}}$ кривая 2$\Delta {{l}_{{{\text{т2}}}}} = 250\,\,{\text{м,}}$ кривая 3$\Delta {{l}_{{{\text{т3}}}}} = 200\,\,{\text{м}}{\text{.}}$

Рис. 10.

Зависимости текущих и требуемых (нумерация со штрихом) линейных разрешений по азимуту от времени для трех вариантов требуемого линейного разрешения: кривая 1$\Delta {{l}_{{{\text{т1}}}}} = 300\,\,{\text{м,}}$ кривая 2$\Delta {{l}_{{{\text{т2}}}}} = 250\,\,{\text{м,}}$ кривая 3$\Delta {{l}_{{{\text{т3}}}}} = 200\,\,{\text{м}}{\text{.}}$

Рис. 11.

Зависимость текущих промахов от времени при законе управления без первого слагаемого для трех вариантов требуемого линейного разрешения: кривая 1$\Delta {{l}_{{{\text{т1}}}}} = 300\,\,{\text{м,}}$ кривая 2$\Delta {{l}_{{{\text{т2}}}}} = 250\,\,{\text{м,}}$ кривая 3$\Delta {{l}_{{{\text{т3}}}}} = 200\,\,{\text{м}}{\text{.}}$

Рис. 12.

Зависимость поперечного ускорения от времени при законе управления без первого слагаемого для трех вариантов требуемого линейного разрешения: кривая 1$\Delta {{l}_{{{\text{т1}}}}} = 300\,\,{\text{м,}}$ кривая 2$\Delta {{l}_{{{\text{т2}}}}} = 250\,\,{\text{м,}}$ кривая 3$\Delta {{l}_{{{\text{т3}}}}} = 200\,\,{\text{м}}{\text{.}}$

Объясняется это, прежде всего, тем, что на больших расстояниях угловые скорости линии визирования настолько малы, что закон управления на них не реагирует. Поэтому эта модификация закона управления начинает выполнять свои функции только после значительного сближения с целью (рис. 12).

ЗАКЛЮЧЕНИЕ

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

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

Для осуществления наведения по закону (20) от бортовой РЛС требуется формирование оценок дальности до цели, скорости сближения с ней, ее бортового пеленга и угловой скорости линии визирования.

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

Список литературы

  1. Попов И.М., Хамзатов М.И. Война будущего: концептуальные основы и практические выводы. М.: Кучково поле, 2017.

  2. Федосов Е.А. // Фазотрон. 2007. № 1, 2. С. 11.

  3. Верба В.С., Иевлев Д.И., Меркулов В.И. // Радиотехника. 2020. № 2(3). С. 5.

  4. Авиация ПВО России и научно-технический прогресс. Боевые комплексы и системы вчера, сегодня, завтра / Под ред. Е.А. Федосова. М.: Дрофа, 2005.

  5. Справочник по радиолокации. Т. 2. Радиолокационные антенные устройства / Под ред. М. Сколника. М.: Сов. радио, 1977.

  6. Григорьев Ф.Н., Кузнецов Н.А., Серебровский Л.П. Управление наблюдением в автоматических системах. М.: Наука, 1986.

  7. Авиационные радиолокационные комплексы и системы / Под ред. П.И. Дудника. М.: Изд. ВВИА им. проф. Н.Е. Жуковского, 2006.

  8. Авиационные системы радиоуправления / Под ред. В.И. Меркулова. М.: Изд-во ВВИА им. проф. Н.Е. Жуковского, 2008.

  9. Меркулов В.И. // Вестник ВКО. 2020. № 1(25). С. 25.

  10. Меркулов В.И., Самарин О.Ф., Францев В.В., Челей Г.С. Способ наведения летательного аппарата на отдельную воздушную цель в составе плотной группы. Патент РФ № 2253082. Опубл. офиц. бюл. “Изобретения. Полезные модели” № 15. 27.05.2005.

  11. Меркулов В.И., Забелин И.В. // Радиотехника. 2012. № 10. С. 107.

  12. Анциферов А.А., Богданов А.В., Бондарев В.Н. и др. // Журн. Сибирского федерального университета. Техника и технологии. 2017. Т. 10. № 2. С. 169.

  13. Информационно-измерительные и управляющие радиоэлектронные системы и комплексы / Под ред. В.С. Вербы. М.: Радиотехника, 2020.

  14. Многофункциональные радиолокационные комплексы истребителей / Под ред. В.Н. Лепина. М.: Радиотехника, 2014.

  15. Меркулов В.И., Верба В.С., Ильчук А.Р., Колтышев Е.Е. Автоматическое сопровождение целей в РЛС интегрированных авиационных комплексов. Т. 2. Сопровождение одиночных целей. М.: Радиотехника, 2018.

  16. Аврамов А.В., Иванов С.Л., Шевченко В.В., Янковский В.Т. Система управления вооружением модернизированного истребителя: учебное пособие. Воронеж: ВУНЦ ВВС “ВВА”, 2019.

Дополнительные материалы отсутствуют.