Способ управления положением оси антенны бортовой радиолокационной станции при сопровождении меневрирующей воздушной цели

Изобретение относится к радиолокации и может быть использовано в радиолокационных станциях (РЛС) для непрерывной селекции, сопровождения интенсивно маневрирующей воздушной цели (МВЦ), в том числе информационного обеспечения процесса наведения летательного аппарата (ЛА). Техническим результатом заявленного изобретения является снижение ошибок управления приводом антенны РЛС ЛА по углу и угловой скорости при сопровождении МВЦ. В заявленном способе выполняют определение оптимальных, однозначно связанных с заданными параметрами привода антенны весовых коэффициентов ошибок управления приводом антенны по углу и угловой скорости рассчитанными через нахождение коэффициентов матрицы усиления сигналов управления, с учетом полной матрицы штрафов за точность в текущий момент времени управления. Расчет исключает трудоемкий, не оптимальный по результату поиск весов ошибок управления методом перебора (методом проб и ошибок). 12 ил.

 

Изобретение относится к радиолокации и может быть использовано в радиолокационных станциях (РЛС) для непрерывной селекции, сопровождения интенсивно маневрирующей воздушной цели (МВЦ), в том числе информационного обеспечения процесса наведения летательного аппарата (ЛА).

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

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

где:

Δϕ=ϕ_ц - φ_AB; Δω=ω_ц - ω_АВ;

ϕ_ц и ω_ц - пеленг цели и угловая скорость линии визирования цели;

φ_AB и ω_AB - угол поворота и угловая скорость поворота оси антенны;

- первая производная угловой скорости линии визирования;

- вторая производная угловой скорости линии визирования;

«∧» - символ оценок параметров;

K - штраф за величину сигнала управления и, обусловленный энергетическими ресурсами привода антенны;

b - коэффициент усиления привода антенны в плоскости управления;

Р₂₁ и Р₂₂ - элементы матрицы весовых коэффициентов текущего вектора состояния, принятые постоянными;

Т_пр - постоянная времени привода антенны.

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

Известен двухдиапазонный следящий угломер, приведенный в [2]. В данном устройстве по наблюдениям направления на цель ε, направления оси ЛА ψ и угла между направлением оси антенны и направлением оси ЛА ϕ находят текущий вектор состояния цели , текущий вектор состояния оси ЛА , где - обозначения угловой скорости и углового ускорения в векторах состояния цели, оси ЛА и оси антенны, соответственно. Сигнал управления приводом антенны u учитывает ошибки сопровождения по углу и угловой скорости , дополнительно учитывает ошибки по угловому ускорению , угловой скорости линии визирования и угловому ускорению линии визирования цели . Сигнал управления приводом антенны находится по выражению:

где - априорно известные постоянные коэффициенты усиления, записанные в запоминающем устройстве.

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

Известен способ управления инерционным приводом антенны в угломере, приведенный в [3, с. 284-285], используемый в качестве прототипа. В данном алгоритме сигнал управления и учитывает ошибки сопровождения по углу и угловой скорости в соответствии с формулой:

где:

- постоянный коэффициент усиления, определяющий вес ошибки управления по углу Δϕ;

- постоянный коэффициент усиления, определяющий вес сигнала ошибки сопровождения по угловой скорости Δω;

q₂₁ и q₂₂ - элементы матрицы штрафов в момент времени t_k окончания управления.

Недостатком способа-прототипа является относительно высокая ошибка управления приводом антенны по углу и угловой скорости при сопровождении МВЦ, связанная с не оптимальным выбором штрафов q₂₁ и q₂₂ за ошибки управления по углу и угловой скорости методом перебора (методом проб и ошибок).

Целью предлагаемого изобретения является снижение ошибок управления приводом антенны РЛС ЛА по углу и угловой скорости при сопровождении МВЦ.

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

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

где:

х_у - n-мерный вектор управляемых параметров привода антенны;

х_m - m-мерный вектор оцениваемых параметров отслеживаемой цели;

ƒ_y и ƒ_m - динамические матрицы состояния процессов (1) и (2);

u - r-мерный вектор сигналов управления (r меньше или равен n);

ξ_у и ξ_m - n-мерные векторы центрированных гауссовских возмущений процессов x_у и х_m;

- обобщенный вектор состояния;

z - m-мерный вектор наблюдений (m меньше или равен n);

H - матрица связей обобщенного вектора х с z;

ξ_u - m-мерные векторы центрированных гауссовских шумов измерений.

Движение ЛА, на котором размещен угломер РЛС, характеризуется относительно наблюдаемой МВЦ в плоскости XOY положением линии визирования (ε), бортовым углом пеленга цели (β); углом поворота антенны относительно продольной оси ЛА (φ_AB); угловым положением оси ЛА (ϑ) относительно оси X, совпадающей с горизонтальной проекцией вектора скорости ЛА в вертикальной (горизонтальной) плоскости (фиг. 1). Точки О и О_ц есть расположения ЛА и МВЦ соответственно в прямоугольной системе координат XOY, ОХ_оy - продольная ось ЛА, ОХ_а - равносигнальное направление, ОХ - горизонтальная ось, параллельная горизонтальной проекции вектора скорости ЛА, OY - предполагается вертикальной при рассмотрении движения ЛА и управления приводом антенны в вертикальной плоскости, при рассмотрении движения ЛА и управления приводом антенны в горизонтальной плоскости - горизонтальная.

Обобщенный вектор состояния х можно представить в виде:

где:

ω - угловая скорость линии визирования МВЦ;

Д и - дальность и скорость сближения с МВЦ;

а_Ц и а_ЛА - нормальные ускорения МВЦ и ЛА;

γ_Ц - коэффициент, учитывающий маневренные свойства цели (величина, обратная постоянной времени маневра цели);

ω_ϑ - угловая скорость положения оси ЛА (ϑ) в вертикальной (горизонтальной) плоскости;

γ_ϑ - коэффициент, характеризующий ширину спектра угловых колебаний оси ЛА;

Т_пр - постоянная времени привода антенны;

b - коэффициент усиления привода антенны;

u - сигнал управления приводом антенны;

ξ_Ц, ξ_ϑ, ξ_АВ - центрированные белые шумы с известными односторонними спектральными плотностями.

Ошибка управления в вертикальной плоскости определяется углом Δφ между равносигнальным направлением антенны и направлением на цель. Тогда для точного слежения за МВЦ необходимо обеспечить такое управление приводом антенны РЛС, чтобы ошибки сопровождения цели по углу Δφ=ε-ϑ-φ_АВ и угловой скорости Δω=ω-ω_ϑ-ω_АВ свести к нулю.

Уравнение (1) с учетом (9) и (10) в матричном виде запишется:

где:

B_y - матрица эффективности управления.

С учетом вектора требуемых параметров выражения (4), (5), (7) и (8) для формирования сигналов управления будем использовать линейно-квадратичный функционал качества [2]:

где:

α и γ - весовые коэффициенты ошибок управления;

K - штраф за величину сигналов управления;

Т - время управления, равное Т_пр - постоянной времени привода антенны; введение обусловлено приведением к единой размерности.

В соответствии с теорией оптимального управления [3, с. 234-251] получено, что минимизация линейно-квадратичного функционала качества (12) достигается при сигнале управления вида:

Выражение (13) - это закон управления, а выражение (14) - дифференциальное уравнение коэффициентов усиления сигналов управления. Требуемое управление приводом антенны РЛС в целях сопровождения МВЦ без срыва должно быть:

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

Разложение подинтегральной суммы первого интеграла (12) дает:

Исходя из выражений (16) и (17) полная матрица штрафов за точность управления примет вид:

С учетом (18) выражение (14) преобразуется к виду:

Определим элементы матрицы d₁₂ и d₂₂ в установившемся режиме, то есть после окончания переходных процессов, вызванных граничными условиями. Представим элементы матрицы D в виде:

Отсюда элементы матрицы будут равны:

Приравняв соответствующим элементам матрицы (19) получим систему уравнений:

Для решения системы уравнений были проведены следующие преобразования. Из (20) выражен с₁₁ и подставлен в (21), которое аналитически решено как квадратичное уравнение и найдены корни уравнения, решением является отрицательный корень, зависящий от с₂₂. Выражение (22) также решено как квадратичное уравнение, решением является отрицательный корень. Значения корней выражений (21) и (22) примут вид:

Совместным решением системы уравнений (23) и (24) находят с₁₂ и с₂₂.

Таким образом, поставляя с₁₂ и с₂₂ в (19), находят матрицу и далее, с учетом (13) и (14), находят закон управления приводом антенны следящего угломера РЛС:

Выражение (25) является оптимальным законом управления приводом антенны РЛС для сопровождения МВЦ, заданной дифференциальными уравнениями.

Предлагаемый способ работы РЛС на борту ЛА включает:

одновременное измерение углового положения линии визирования цели ε в вертикальной (горизонтальной) плоскости, скорости ее изменения ω, углового положения оси ЛА ϑ в вертикальной (горизонтальной) плоскости, скорости его изменения ω_ϑ в вертикальной (горизонтальной) плоскости, углового положения оси антенны относительно продольной оси ЛА (в дальнейшем называемом углом поворота) φ_АВ в вертикальной (горизонтальной) плоскости и скорости его изменения в вертикальной (горизонтальной) плоскости ω_AB, формирование сигналов ошибок сопровождения цели по углу Δφ и угловой скорости Δω, значение ошибки сопровождения цели по углу Δφ получают из выражения:

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

сигнал управления приводом антенны и получают весовой суммой ошибок сопровождения цели по углу Δφ и угловой скорости Δω с постоянными весами ω_Δφ и ω_Δω соответственно, зависящими от свойств привода антенны:

отличающийся тем, что в качестве w_Δφ и w_Δω, используются значения, определяемые расчетом в последовательности: находят значения коэффициентов с₁₂ и с₂₂ из системы уравнений и подставляют их в выражения для расчета весов w_Δφ и w_Δω по выражениям:

>

где:

b - коэффициент усиления привода антенны в вертикальной (горизонтальной) плоскости;

α и γ - весовые коэффициенты ошибок управления в вертикальной (горизонтальной) плоскости;

K - штраф за величину сигнала управления в вертикальной (горизонтальной) плоскости;

Т_пр - постоянная времени привода антенны.

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

Сущность предлагаемого способа управления приводом антенны угломера РЛС ЛА при сопровождении МВЦ поясняется дальнейшим описанием и чертежами.

На фиг. 1 представлена геометрия слежения за МВЦ.

На фиг. 2 представлены результаты моделирования отработки угла поворота антенны предлагаемым способом.

На фиг. 3 представлены результаты моделирования отработки угла поворота антенны способом-прототипом.

На фиг. 4 представлены результаты моделирования отработки угловой скорости поворота антенны предлагаемым способом.

На фиг. 5 представлены результаты моделирования отработки угловой скорости поворота антенны способом-прототипом.

На фиг. 6 представлена ошибка по углу поворота антенны предлагаемым способом и способом-прототипом.

На фиг. 7 представлена ошибка по угловой скорости поворота антенны предлагаемым способом и способом-прототипом.

На фиг. 8 представлены траектории движения ЛА и МВЦ, выполняющей маневр типа «змейка».

На фиг. 9 представлена ошибка по углу поворота антенны предлагаемым способом и способом-прототипом при сопровождении МВЦ, выполняющей маневр типа «змейка».

На фиг. 10 представлена ошибка по угловой скорости поворота антенны предлагаемым способом и способом-прототипом при сопровождении МВЦ, выполняющей маневр типа «змейка».

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

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

1) Δϕ(0)<0, Δω(0)<0;

2) Δϕ(0)>0, Δω(0)>0;

3) Δϕ(0)<0, Δω(0)>0;

4) Δϕ(0)>0, Δω(0)<0.

Предлагаемый способ формирования сигналов управления приводом антенны РЛС при сопровождении МВЦ в вертикальной (горизонтальной) плоскости реализуется следующим образом.

Одновременно в вертикальной (горизонтальной) плоскости измеряют значения угла линии визирования цели ε и скорость ее изменения ω, угловое положение оси ЛА ϑ в вертикальной (горизонтальной) плоскости, скорости его изменения ω_ϑ в вертикальной (горизонтальной) плоскости, угол поворота антенны РЛС φ_АВ и скорость ее изменения ω_AB.

Определяют требуемое значение угла поворота антенны РЛС Δφ, соответствующее ошибке сопровождения цели по углу по выражению:

Δφ=ε-ϑ-φ_АВ,

значение ошибки сопровождения цели по угловой скорости получают из выражения:

Δω=ω-ω_ϑ-ω_AB.

Из системы уравнений

находят значения коэффициентов с₁₂ и с₂₂.

Вычисляют весовые коэффициенты w_Δφв и w_Δωв по формулам:

Формируют сигнал управления приводом антенны РЛС в вертикальной (горизонтальной) плоскости по выражению:

Для определения эффективности предлагаемого способа было проведено математическое моделирование процесса управления приводом антенны РЛС при сопровождении МВЦ в вертикальной плоскости, а также способа-прототипа, в котором сигнал управления определен по выражению

где:

При этом рассматривалась типовые траектории движения ЛА и МВЦ на встречных курсах со следующими параметрами моделирования: скорость ЛА 300 м/с и МВЦ 250 м/с; начальные координаты ЛА (0;10⁴) и МВЦ (2×10⁴;1,1×10⁴); параметры привода антенны для обоих случаев: коэффициент усиления привода антенны b=1 и постоянная времени привода антенны Т_пр=0,5 с. Расчеты проводились при одинаковых весовых коэффициентах ошибок управления. На фиг. 2, 3, 4, 5 представлены зависимости изменения углов поворота для предлагаемого способа и способа прототипа, угловых скоростей поворота для предлагаемого способа и способа-прототипа с учетом требуемых значений по углу и угловой скорости.

В качестве показателей эффективности рассматривался модуль абсолютных ошибок сопровождения по углу и угловой скорости. Результаты моделирования в виде временных зависимостей показателей эффективности представлены на фиг. 6 и 7. Их анализ показывает, что при использовании предлагаемого способа ошибки управления снижаются. Так, в рассмотренных условиях ошибка сопровождения по углу при t больше 10 с снижается до 80%, а по угловой скорости до 50% по сравнению с прототипом. В интересах определения типа сопровождаемых целей проведено моделирование процесса сопровождения МВЦ, выполняющей маневр типа «змейка», характерная для гиперзвуковых ЛА. Траектории Л А и МВЦ показаны на фиг. 8. На фиг. 9 и 10 представлены ошибки сопровождения по углу и угловой скорости при слежении за МВЦ, выполняющей маневр типа «змейка» для предлагаемого способа и способа-прототипа. Анализ данных зависимостей показывает, что предлагаемый способ возможно использовать для управления приводом антенны РЛС в интересах сопровождения высокоскоростных, в том числе гиперзвуковых, целей.

На фиг. 11 и 12 показаны ошибки сопровождения по углу и угловой скорости при различных сочетаниях начальных условий:

1) Δϕ(0)<0, Δω(0)<0;

2) Δϕ(0)>0, Δω(0)>0; Δϕ(0)<0, Δω(0)>0; Δϕ(0)>0, Δω(0)<0.

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

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

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

Источники информации:

1. Патент России 2571363 от 20.12.2015, кл. H01Q 25/00, G01S 13/66. Способ управления инерционным приводом антенны, обеспечивающий устойчивое сопровождение интенсивно маневрирующих и высокоскоростных воздушных объектов.

2. Меркулов В.И., Лепин В.Н. Авиационные системы радиоуправления. - М.: Радио и связь. 1996. 396 с.

3. Бухалев В.А. Основы автоматики и теории управления. - М.: Изд. ВВИА им. Н.Е. Жуковского. 2006. 406 с.

Способ управления положением оси антенны бортовой РЛС, обеспечивающий устойчивое сопровождение интенсивно маневрирующих и высокоскоростных летательных аппаратов, включающий одновременное измерение углового положения линии визирования цели ε в вертикальной (горизонтальной) плоскости, скорости ее изменения ω, углового положения оси летательного аппарата (ЛA) ϑ в вертикальной (горизонтальной) плоскости, скорости его изменения ω_ϑ в вертикальной (горизонтальной) плоскости относительно горизонтальной проекции вектора скорости ЛА, углового положения оси антенны относительно продольной оси ЛА (в дальнейшем называемом углом поворота) φ_АВ вертикальной (горизонтальной) плоскости и скорости его изменения ω_АВ в вертикальной (горизонтальной) плоскости, формирование ошибок сопровождения цели по углу Δφ и угловой скорости Δω из выражений

Δφ=ε-ϑ-φ_АВ,

Δω=ω-ω_ϑ-ω_АВ;

сигнал управления приводом антенны получают весовой суммой ошибок сопровождения цели по углу Δφ и угловой скорости цели Δω с весами по углу w_Δφ и угловой скорости w_Δω по выражению

u₁=w_ΔφΔφ+w_ΔωΔω;

отличающийся тем, что в качестве весов по углу w_Δφ и угловой скорости w_Δω используются значения, определяемые расчетом в последовательности: находят из системы уравнений

значения коэффициентов с₁₂ и с₂₂; вычисляют веса ошибок по углу w_Δφ и угловой скорости цели w_Δω с учетом найденных коэффициентов c₁₂ и c₂₂ по выражениям

где b - коэффициент усиления привода антенны в вертикальной (горизонтальной) плоскости,

α и γ - весовые коэффициенты ошибок управления в вертикальной (горизонтальной) плоскости,

K - штраф за величину сигнала управления в вертикальной (горизонтальной) плоскости,

Т_пр - постоянная времени привода антенны.