СПОСОБ ОПРЕДЕЛЕНИЯ ДИНАМИЧЕСКОЙ ЭФФЕКТИВНОЙ ПЛОЩАДИ РАССЕЯНИЯ ГИПЕРЗВУКОВОЙ КРЫЛАТОЙ РАКЕТЫ
METHOD FOR DETERMINING DYNAMIC RADAR CROSS-SECTION OF HYPERSONIC CRUISE MISSILE

Изобретение относится к области моделирования и может быть использовано для определения зависимости эффективной площади рассеяния гиперзвуковой крылатой ракеты от времени.
В настоящее время для оценки эффективности зенитных ракетных комплексов (ЗРК) используются имитационные модели боевых действий, позволяющие исследовать процессы функционирования сложных систем с сохранением их логической структуры и последовательности чередования происходящих событий во времени. В свою очередь, эффективность ЗРК зависит от качества информации о противнике, получаемой от информационных средств, к которым прежде всего относятся радиолокационные станции и комплексы. При этом качество радиолокационной информации определяется вероятностью правильного обнаружения цели при фиксированном значении ложной тревоги, а также точностью определения ее координат и зависит от значения эффективной площади рассеяния (ЭПР) цели. В тоже время, в известных моделях боевых действий ЗРК процессы обнаружения воздушных целей, как правило, моделируются с использованием средних или медианных значений ЭПР в широком диапазоне ожидаемых ракурсов облучения цели. Такой подход позволяет решать задачи оценки эффективности группировок войск при ограничениях на вычислительные затраты, однако, снижает точность оценки дальности обнаружения целей, так как незначительное изменение ракурса облучения цели ведет к изменению ее ЭПР в десятки раз. Особенно это важно учитывать при моделировании процесса обнаружения перспективных гиперзвуковых целей, ракурс облучениях которых, в связи с высокой скоростью их полета (V=6-8 M), изменяется значительно быстрее, чем обычных средств воздушно-космического нападения, обуславливая непрерывное изменение ЭПР во времени.
В связи с вышеизложенным, целью настоящего изобретения является определение зависимости эффективной площади рассеяния гиперзвуковой крылатой ракеты от времени.
Известен способ моделирования радиолокационных сигналов, отраженных от космических объектов на внеатмосферном участке полета [Моделирование в радиолокации / А.И. Леонов, В.Н. Васенев, Ю.И. Гайдуков и др.; Под ред. А.И. Леонова. - М.: Сов. радио, 1979, рис. 4.4, с. 99], позволяющий определить эффективную площадь рассеяния только для космического объекта на баллистическом участке его полета.
Наиболее близким техническим решением, принятым за прототип, является способ моделирования динамических характеристик рассеяния с использованием статических радиолокационных характеристик [Радиолокационные характеристики летательных аппаратов / М.Е. Варганов, Ю.С. Зиновьев, Л.Ю. Астанин и др.; Под ред. Л.Т. Тучкова. - М.: Радио и связь, 1985, рис. 9.1, с. 174], который заключается в следующем.
ц
Экспериментально измеряется или аналитически рассчитывается диаграмма обратного отражения радиолокационного сигнала от цели σ(θ) и моделируется движение цели относительно РЛС обнаружения для получения ракурса ее облучения в различных точках траектории полета:
n
n
n
n
n
n
где ϕ, Δ, γ - углы ориентации цели относительно линии визирования РЛС (ϕ - угол нутации, λ - угол прецессии и γ- угол собственного вращения).
При этом углы ориентации рассчитываются с использованием элементов матрицы М направляющих косинусов между связанной и визирной системами координат [Радиолокационные характеристики летательных аппаратов / М.Е. Варганов, Ю.С. Зиновьев, Л.Ю. Астанин и др.; Под ред. Л.Т. Тучкова. - М.: Радио и связь, 1985, с. 189]:
p
где С - матрица перехода от топоцентрической к визирной системе координат, D - матрица перехода от топоцентрической системы координат к абсолютной, А - матрица перехода от абсолютной к орбитальной системе координат, В - матрица направляющих косинусов между связанной и орбитальной системами координат.
Элементы матриц А и В получают путем расчета составляющих кеплеровой орбиты и текущих параметров движения летательного аппарата (ЛА) в геоцентрической абсолютной системе координат. Однако, траекторию полета гиперзвуковой крылатой ракеты нельзя описать кеплеровыми зависимостями, в связи с чем применение известного способа не позволяет достичь поставленной цели.
Техническим результатом изобретения является создание способа определения динамической эффективной площади рассеяния гиперзвуковой крылатой ракеты.
Достижение заявленного технического результата и поставленной цели обеспечивается тем, что в способ определения динамической эффективной площади рассеяния гиперзвуковой крылатой ракеты, заключающийся в том, что в текущий момент времени рассчитываются значения углов нутации, прецессии и собственного вращения, характеризующие ориентацию цели относительно линии визирования радиолокационной станции обнаружения, в зависимости от которых определяется значение эффективной площади рассеяния гиперзвуковой крылатой ракеты, дополнительно включены следующие этапы: задание исходных данных для моделирования, расчет координат центра масс ГЗКР в геоцентрической сферической системе координат, расчет элементов матрицы направляющих косинусов между местной географической системой координат и геоцентрической гринвичской прямоугольной системой, расчет элементов матрицы направляющих косинусов между полускоростной системой координат и местной географической системой координат, расчет элементов матрицы направляющих косинусов между полускоростной системой координат и связанной системой координат, пересчет координат радиолокационной станции обнаружения из геоцентрической гринвичской прямоугольной в геодезическую систему координат, расчет элементов матрицы направляющих косинусов между геоцентрической гринвичской прямоугольной системой координат и топоцентрической (измерительной) системой, расчет координат центра масс гиперзвуковой крылатой ракеты в топоцентрической (измерительной) системе координат радиолокационной станции, вычисление элементов матрицы перехода от топоцентрической (измерительной) системы координат к визирной, расчет матрицы направляющих косинусов между связанной и визирной системами координат путем умножения транспонированной матрицы перехода от топоцентрической системы координат к визирной, матрицы направляющих косинусов между геоцентрической гринвичской прямоугольной системой координат и топоцентрической (измерительной) системой, транспонированной матрицы направляющих косинусов между местной географической системой координат и геоцентрической гринвичской, прямоугольной системой, транспонированной матрицы направляющих косинусов между полускоростной системой координат и местной географической системой координат, транспонированной матрицы направляющих косинусов между полускоростной системой координат и связанной системой координат.
Сущность изобретения поясняется рисунками, представленными на фиг. 1-3. На фиг. 1 представлена структурная схема способа определения динамической эффективной площади рассеяния гиперзвуковой крылатой ракеты; на фиг. 2 - устройство, реализующее способ определения динамической эффективной площади рассеяния гиперзвуковой крылатой ракеты; на фиг. 3 - зависимость ЭПР гиперзвуковой крылатой ракеты от времени.
Способ определения динамической эффективной площади рассеяния гиперзвуковой крылатой ракеты (фиг. 1) заключается в следующем:
1) на первом этапе способа задают исходные данные для моделирования:
i
- параметры траектории полета ГЗКР: множество координат, определяющих пространственную ориентацию ГЗКР в каждый момент времени ее полета t с заданной дискретностью
i
- параметры траектории полета ГЗКР: множество координат, определяющих пространственную ориентацию ГЗКР в каждый момент времени ее полета t с заданной дискретностью
ц
- диаграмма обратного отражения радиолокационного сигнала от цели σ(θ);
p
p
p
- координаты РЛС обнаружения x, y, z;
i
ц
ц
c
c
c
2) на втором этапе способа по полученным в момент времени tкоординатам центра масс ГЗКР х, у, z в геоцентрической гринвичской прямоугольной системе координат осуществляется расчет координат центра масс в геоцентрической сферической системе координат - радиуса r, углов λи ϕ (геоцентрические долгота и широта), Определяющих взаимную ориентацию осей геоцентрической гринвичской прямоугольной Oxyz и местной географической систем координат Cxyz:
Затем рассчитываются элементы матрицы Р направляющих косинусов между местной географической системой координат и геоцентрической гринвичской прямоугольной системой:
i
3) на третьем этапе способа по полученным в момент времени t значениям углов наклона траектории θ и поворота траектории ψ рассчитываются элементы матрицы L направляющих косинусов между полускоростной системой координат и местной географической системой координат:
i
4) на четвертом этапе способа по полученным в момент времени t значениям углов атаки α и скольжения β рассчитываются элементы матрицы S направляющих косинусов между полускоростной системой координат и связанной системой координат:
5) на пятом этапе способа для расчета координат центра масс ГЗКР относительно РЛС обнаружения осуществляется пересчет координат РЛС обнаружения из геоцентрической гринвичской прямоугольной в геодезическую систему координат. Для этого используются следующие выражения:
После этого рассчитываются элементы матрицы N направляющих косинусов между геоцентрической гринвичской прямоугольной системой координат и топоцентрической (измерительной) системой:
i
Далее в соответствии с заданными в исходных данных координатами РЛС и полученными в момент времени t координатами центра масс ГЗКР рассчитываются координаты центра масс ГЗКР в топоцентрической (измерительной) системе координат РЛС, сначала в прямоугольной:
а затем - в сферической:
гзкр
гзкр
гзкр
где d - наклонная дальность до цели, ε - угол места цели, β - азимут цели;
n
n
n
i
6) на шестом этапе способа производится расчет углов ориентации ГЗКР относительно линии визирования РЛС - углов нутации ϕ, прецессии λ и собственного вращения γ в момент времени t.
Для этого вычисляются элементы матрицы С перехода от топоцентрической (измерительной) системы координат к визирной:
Для этого вычисляются элементы матрицы С перехода от топоцентрической системы координат к визирной:
где
n
n
n
i
Далее рассчитываются значения углов ϕ, λ и γ на момент времени t. Для этого осуществляется расчет матрицы направляющих косинусов между связанной и визирной системой координат (матрицы М), при этом выражение (2) преобразуется к виду:
Элементы матрицы М используются для расчета искомых углов:
11
12
21
где m, m, m - соответствующие элементы матрицы М;
i
7) на седьмом этапе способа определяется значение динамической эффективной площади рассеяния ГЗКР на момент времени t в соответствии с рассчитанным на предыдущем этапе ракурсом облучения ГЗКР  и заданной в исходных данных диаграммой обратного отражения цели
Устройство (фиг. 2), реализующее способ определения динамической эффективной площади рассеяния гиперзвуковой крылатой ракеты, содержит соединенные между собой цифровыми линиями связи интерфейс 1 ввода исходных данных и анализа результатов моделирования, блок 2 расчета динамической эффективной площади рассеяния ГЗКР, блок 3 памяти исходных данных и результатов моделирования.
Интерфейс 1 ввода исходных данных и анализа результатов моделирования выполнен на базе ПЭВМ и содержит процессор, оперативное запоминающее устройство (ОЗУ), перепрограммируемое запоминающее устройство, вводно-выводное устройство (ВВУ) и дисплей. Интерфейс 1 ввода исходных данных и анализа результатов моделирования соединен по сигналам коммутации с блоком 2 расчета динамической эффективной площади рассеяния ГЗКР и блоком 3 памяти исходных данных и результатов модельных экспериментов.
Блок 2 расчета динамической эффективной площади рассеяния ГЗКР выполнен в виде программно-технического модуля на плате расширения ПЭВМ. Блок 2 расчета динамической эффективной площади рассеяния ГЗКР соединен по сигналам коммутации с интерфейсом 1 ввода исходных данных и анализа результатов моделирования и блоком 3 памяти.
Блок 3 памяти исходных данных и результатов модельных экспериментов содержит в цифровой форме наборы данных: по параметрам траектории полета ГЗКР, диаграмме обратного отражения цели, координатам РЛС обнаружения; и соединен по сигналам коммутации с интерфейсом 1 ввода исходных данных и анализа результатов моделирования, блоком 2 расчета динамической эффективной площади рассеяния ГЗКР.
Устройство функционирует следующим образом. С помощью интерфейса 1 ввода исходных данных и анализа результатов моделирования пользователем осуществляется задание исходных данных для моделирования, их запись в блок 3 памяти исходных данных и результатов моделирования и запуск моделирования. Далее в блоке 2 расчета динамической эффективной площади рассеяния ГЗКР осуществляются расчеты согласно этапов 2-7 способа. Затем результаты расчетов выводятся на дисплей интерфейса 1 ввода исходных данных и анализа результатов моделирования.
Таким образом, предложенный способ позволяет определить зависимость эффективной площади рассеяния гиперзвуковой крылатой ракеты от времени, что подтверждается расчетами, проведенными с помощью устройства. Результаты расчетов для конечного участка траектории полета (время полета 120 с) представлены на фиг. 3.