RU2617110C1 - Способ сопровождения в радиолокационной станции групповой воздушной цели из класса "самолёты с турбореактивными двигателями" при воздействии уводящих по скорости помех - Google Patents
Способ сопровождения в радиолокационной станции групповой воздушной цели из класса "самолёты с турбореактивными двигателями" при воздействии уводящих по скорости помех Download PDFInfo
- Publication number
- RU2617110C1 RU2617110C1 RU2016108454A RU2016108454A RU2617110C1 RU 2617110 C1 RU2617110 C1 RU 2617110C1 RU 2016108454 A RU2016108454 A RU 2016108454A RU 2016108454 A RU2016108454 A RU 2016108454A RU 2617110 C1 RU2617110 C1 RU 2617110C1
- Authority
- RU
- Russia
- Prior art keywords
- group
- signal
- aircraft
- estimates
- doppler
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
- G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional [2D] tracking, e.g. combination of angle and range tracking, track-while-scan radar
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
Изобретение относится к области радиолокации и может быть использовано в радиолокационной станции (РЛС) для сопровождения групповой воздушной цели из класса «самолеты с турбореактивными двигателями» при воздействии уводящих по скорости помех. Достигаемый технический результат - повышение достоверности оценок доплеровских частот (ДЧ), обусловленных скоростью сближения носителя РЛС с каждым самолетом группы при воздействии уводящих по скорости помех. Способ заключается в параллельном сопровождении на основе калмановской фильтрации отсчетов ДЧ, обусловленных отражениями сигнала от планеров самолетов группы и центроида отсчетов ДЧ, обусловленных отражениями сигнала от лопаток рабочего колеса компрессора низкого давления двигателей самолетов; идентификации воздействия или отсутствия уводящих по скорости помех на основе вычисления модулей производных оценок разностей между оценками ДЧ, обусловленными отражениями сигнала от планера каждого самолета группы и центроидом ДЧ, обусловленных отражениями сигнала от лопаток рабочего колеса первых ступеней компрессора низкого давления двигателей самолетов группы; сравнении модулей производных оценок разностей ДЧ с порогом; при их непревышении установленного порога, что соответствует отсутствию воздействия уводящих по скорости помех, на выходе формируются оценки ДЧ, вычисляемые в соответствии с процедурой калмановской фильтрации на основе наблюдения, в противном случае принимается решение о воздействии уводящих по скорости помех и на выходе наряду с оценками ДЧ, которые не идентифицированы как уводящие по скорости помехи, формируются оценки ДЧ, вычисляемые на основе модели взаимного перемещения носителя РЛС и того самолета группы, отраженный от которого сигнал изначально еще не был идентифицирован как уводящая по скорости помеха. 2 ил.
Description
Изобретение относится к области радиолокации и может быть использовано в радиолокационной станции (РЛС) для сопровождения групповой воздушной цели (ГВЦ) из класса «самолеты с турбореактивными двигателями (ТРД)» при воздействии уводящих по скорости помех.
Известен способ сопровождения групповой воздушной цели, заключающийся в отслеживании ее центроида (среднего кинематического поведения группы) и боковых траекторий, распознавании на основе сравнения переменных состояния центральной и боковых траекторий отделяющихся целей от группы [1].
Недостатком данного способа сопровождения групповой воздушной цели является низкая достоверность оценок доплеровских частотах, обусловленных скоростью сближения носителя РЛС с каждым самолетом из состава группы при воздействии уводящих по скорости помех.
Известен способ сопровождения групповой воздушной цели из класса «самолеты с турбореактивными двигателями», заключающийся в том, что сигнал, отраженный от нее, на промежуточной частоте с выхода приемника РЛС подвергается узкополосной доплеровской фильтрации на основе процедуры быстрого преобразования Фурье (БПФ) и преобразуется в амплитудно-частотный спектр, составляющие которого обусловлены отражениями сигнала от планеров самолетов группы и вращающихся лопаток рабочего колеса компрессора низкого давления (КНД) их силовых установок, определяется отсчет доплеровской частоты, соответствующий максимальной амплитуде спектральной составляющей спектра сигнала, который обусловлен его отражением от планера одного из самолетов из их группы, слева и справа в полосе частот ±ΔF, где ΔF - априорно заданная полоса доплеровских частот, занимаемая спектральными составляющими, обусловленными отражениями сигнала от планеров самолетов группы, относительно частоты спектральной составляющей спектра сигнала, имеющей максимальную амплитуду, определяются частоты локальных максимумов спектра сигнала и их количество, которые превысили установленный порог, по выявленным отсчетам доплеровских частот локальных максимумов, находящимся в полосе частот ±ΔF относительно спектральной составляющей, имеющей максимальную амплитуду, вычисляется отсчет доплеровской частоты центроида, как среднее значение отсчетов доплеровских частот локальных максимумов, который поступает на вход оптимального фильтра сопровождения центроида групповой воздушной цели, работающего в соответствии с процедурой оптимальной многомерной линейной дискретной калмановской фильтрации, описываемой выражениями
где
K - общее количество тактов работы фильтра;
Р-(k+1) и P(k+1) - ковариационные матрицы ошибок экстраполяции и фильтрации соответственно;
Ф(k) - переходная матрица состояния;
Q(k+1) и R(k+1) - ковариационные матрицы шумов возбуждения и наблюдения соответственно;
S(k+1) - матрица весовых коэффициентов;
I - единичная матрица;
и - вектор текущих и экстраполированных оценок доплеровской частоты, обусловленной скоростью сближения центроида групповой воздушной цели с носителем станции ее сопровождения;
Н(k) - матрица наблюдения;
Y(k) - вектор наблюдения отсчетов доплеровских частот;
Z(k+1) - матрица невязок измерения;
Ψ(k+1) - матрица априорных ошибок фильтрации;
"-1" - операция вычисления обратной матрицы;
"т" - операция транспонирования матрицы, определяются отсчеты доплеровских частот локальных максимумов, имеющих наибольшие амплитуды, количество которых равно количеству локальных максимумов, расположенных в полосе доплеровских частот ±ΔF относительно спектральной составляющей с максимальной амплитудой, соответствующей отражению сигнала от планера одного из самолетов группы, и находящихся справа на частотах, превышающих значение ΔF, относительно частоты спектральной составляющей, имеющей максимальную амплитуду, которые поступают на вход оптимального фильтра сопровождения первых компрессорных составляющих спектра сигнала, обусловленных его отражениями от лопаток рабочего колеса первой ступени КНД двигателя каждого j-го, где ; J - количество самолетов в группе, самолета группы и работающего аналогично, как и оптимальный фильтр сопровождения центроида групповой воздушной цели, в соответствии с процедурой (1)-(6), при этом размерность матриц, входящих в процедуру (1)-(6), определяется количеством локальных максимумов спектральных составляющих сигнала, находящихся в полосе частот ±ΔF относительно спектральной составляющей с максимальной амплитудой, соответствующей отражению сигнала от планера одного из самолетов группы, на каждом k-ом такте работы обоих оптимальных фильтров сопровождения определяются оценки разности между оцененными значениями доплеровских частот, соответствующих центроиду групповой воздушной цели и отражениям сигнала от лопаток рабочего колеса первой ступени компрессора низкого давления двигателя каждого j-го самолета группы, которые соответствуют только одному из распознаваемых типу воздушной цели из класса «самолеты с турбореактивными двигателями», при этом весь диапазон возможных значений оценок разностей априорно разбивается на Q неперекрывающихся друг с другом поддиапазонов, нижняя FHq и верхняя FBq границы каждого q-го поддиапазона, , соответствующего i-му типу цели, где ; I - максимальное количество распознаваемых типов воздушных целей из класса «самолеты с турбореактивными двигателями», определяются выражениями
где
FPi - максимальная частота вращения ротора КНД силовой установки i-го типа самолета;
n1 и n2 - соответственно минимальное и максимальное значение величины относительных оборотов вращения ротора силовой установки, одинаковые для всех типов самолетов группы;
Nлi - количество лопаток рабочего колеса первой ступени КНД, вычисляется за К тактов работы обоих оптимальных фильтров вероятности Pqj попадания величин в каждый из априорно сформированный q-й поддиапазон, определяются номера i-х поддиапазонов, для которых величины вероятностей Pqj максимальны, максимальные значения величин Pqj max сравниваются с пороговым значением вероятности распознавания типа каждого самолета группы Pпор, при Pqj max≥Рпор принимается решение о том, что j-й самолет в группе имеет i-й тип из класса «самолеты с турбореактивными двигателями» с вероятностью Pqj max, не ниже заданной, в противном случае принимается решение о невозможности распознать тип самолета в группе с заданной вероятностью, количество величин разностей , попавших в q-й поддиапазон с вероятностями Pqj max, определяет количество самолетов i-го типа в группе [2].
Недостатком данного способа сопровождения групповой воздушной цели является низкая достоверность оценок доплеровских частот, обусловленных скоростью сближения носителя РЛС с каждым j-м, , самолетом группы при воздействии уводящих по скорости помех.
Цель изобретения - повышение достоверности оценок доплеровских частот, обусловленных скоростью сближения носителя РЛС с каждым самолетом группы при воздействии уводящих по скорости помех.
Указанная цель достигается тем, что в способе сопровождения в РЛС групповой воздушной цели из класса «самолеты с турбореактивными двигателями» заключающемся в том, что сигнал, отраженный от ГВЦ, на промежуточной частоте с выхода приемника РЛС подвергается узкополосной доплеровской фильтрации на основе процедуры БПФ и преобразуется в амплитудно-частотный спектр, составляющие которого обусловлены отражениями сигнала от планеров самолетов группы и вращающихся лопаток рабочего колеса КНД их силовых установок, определяется отсчет доплеровской частоты, соответствующей максимальной амплитуде спектральной составляющей спектра сигнала, который обусловлен его отражением от планера одного из самолетов группы, определяются слева и справа в полосе частот ±ΔF, где ΔF - априорно заданная полоса доплеровских частот, занимаемая спектральными составляющими, обусловленными отражениями сигнала от планеров самолетов группы, относительно частоты спектральной составляющей спектра сигнала, имеющей максимальную амплитуду, отсчеты доплеровских частот локальных максимумов спектра сигнала и их количество, которые превысили установленный порог, определяются отсчеты доплеровских частот локальных максимумов, имеющих наибольшие амплитуды, количество которых равно количеству локальных максимумов, расположенных в полосе доплеровских частот ±ΔF относительно спектральной составляющей с максимальной амплитудой, соответствующей отражению сигнала от планера одного из самолетов группы, и находящихся справа на частотах, превышающих значение ΔF, относительно частоты спектральной составляющей, имеющей максимальную амплитуду, которые обусловлены отражениями сигнала от лопаток рабочего колеса первых ступеней компрессора низкого давления двигателей самолетов группы, дополнительно вычисляется процедура оптимальной многомерной линейной дискретной калмановской фильтрации в соответствии с выражениями (1)-(6) при сопровождении по доплеровской частоте каждой j-й, , планерной составляющей спектра сигнала, обусловленной скоростью сближения каждого j-го, , самолета группы с носителем РЛС, при этом, размерность матриц, входящих в процедуру (1)-(6), определяется количеством локальных максимумов спектральных составляющих сигнала, находящихся в полосе частот ±ΔF относительно спектральной составляющей с максимальной амплитудой, соответствующей отражению сигнала от планера одного из самолетов группы, по выявленным отсчетам доплеровских частот локальных максимумов, находящихся справа на частотах, превышающих значение ΔF и соответствующих отражениями сигнала от лопаток рабочего колеса первых ступеней КНД двигателей самолетов группы, вычисляется отсчет доплеровской частоты центроида, как среднее значение отсчетов доплеровских частот локальных максимумов, который поступает на вход оптимального фильтра сопровождения центроида отсчетов доплеровских частот, обусловленных отражением сигнала от лопаток рабочего колеса первых ступеней КНД двигателей самолетов группы, работающего аналогично, как и оптимальный фильтр сопровождения планерных составляющих спектра сигнала в соответствии с процедурой (1)-(6), на каждом k-ом такте работы обоих оптимальных фильтров сопровождения вычисляются оценки разностей между оцененными значениями доплеровских частот, соответствующих центроиду доплеровской частоты , обусловленной отражениями сигнала от лопаток рабочего колеса первых ступеней КНД двигателей самолетов группы, и отражениям сигнала от планера каждого j-го самолета группы, , вычисляются модули производных оценок разностей между оцененными значениями доплеровских частот , обусловленными отражениями сигнала от планера каждого j-го самолета группы и центроидом доплеровской частоты ), обусловленной отражениями сигнала от лопаток рабочего колеса первых ступеней КНД двигателей самолетов группы, которые сравниваются с пороговым значением ε, близким к нулю, при выполнении условия для всех j-ых, , оценок планерных составляющих спектра сигнал, что соответствует отсутствию воздействия уводящих по скорости помех, в качестве выходной информации используются все оценки , , доплеровских частот, вычисляемые в соответствии с процедурой (1)-(6) на основе наблюдения Y(k), при не выполнении условия для всех , что соответствует воздействию уводящих по скорости помех, наряду с оценками доплеровских частот , где ; M - количество оценок доплеровских частот, которые не идентифицированы, как уводящие по скорости помехи, формируются оценки доплеровских частот , где ; N - количество оценок доплеровских частот, которые идентифицированы, как уводящие по скорости помехи; J=M+N, вычисляемые в соответствии с процедурой
где
- переходная матрица состояния, численные значения элементов которой определяются на основе корреляционного анализа текущих значений оценок доплеровских частот , , полученных до момента времени, когда выполнялось условие и идентифицированные затем, как уводящие по скорости помехи;
- оцененный вектор состояния, содержащий оценку доплеровской частоты , , полученную после ее идентификации, как уводящей по скорости помехи,
только на основе переходной матрицы состояния , учитывающей модель взаимного перемещения носителя РЛС и того самолета группы, отраженный от которого сигнал изначально еще не был идентифицирован, как уводящая по скорости помеха.
Новыми признаками, обладающими существенными отличиями, являются.
1. Идентификация отсутствия или воздействия уводящих по скорости помех на основе анализа совокупности модулей производных оценок разностей величин . При их не превышении установленного порога ε, близкого к нулю для всех оценок доплеровских частот , принимается решение об отсутствии воздействия уводящих по скорости помех, в противном случае, при выполнении условий либо для всех, либо части оценок доплеровских частот из их совокупности, принимается соответственно решение о том, что либо все, либо часть доплеровских частот из их совокупности, обусловленные скоростями сближения носителя РЛС с каждым j-м, , самолетом группы идентифицируются, как уводящие по скорости помехи.
2. Формирование достоверных оценок доплеровских частот , , обусловленных скоростью сближения носителя РЛС с каждым самолетом группы в соответствии с процедурой (1)-(6) при полном отсутствии уводящих по скорости помех, когда условия выполняются для всех оценок доплеровских частот или оценок , при идентификации M оценок без воздействия и , при идентификации N оценок при воздействии уводящих по скорости помех, когда условия соответственно выполняются только для М и не выполняются для N оценок доплеровских частот.
Данные признаки обладают существенными отличиями, так как в известных способах не обнаружены.
Применение новых признаков позволит идентифицировать воздействие или отсутствие уводящих по скорости помех при сопровождении групповой воздушной цели в РЛС и, в зависимости от результата идентификации, сформировать достоверные оценки доплеровских частот, обусловленные скоростью сближения носителя РЛС с каждым самолетом группы при воздействии уводящих по скорости помех.
На рисунке 1 приведена блок-схема, поясняющая предлагаемый способ сопровождения в РЛС групповой воздушной цели из класса «самолеты с турбореактивными двигателями» при воздействии уводящих по скорости помех, на рисунке 2 (а, б, в, г, д) - эпюры, поясняющие предлагаемый способ на примере сопровождения групповой воздушной цели, состоящей их двух самолетов, и воздействии одной уводящей по скорости помехи.
Способ сопровождения в РЛС групповой воздушной цели из класса «самолеты с турбореактивными двигателями» при воздействии уводящих по скорости помех осуществляется следующим образом.
На вход блока 1 БПФ (рисунок 1) на промежуточной частоте с выхода приемника РЛС поступает сигнал S(t) (рисунок 2а), отраженный от ГВЦ, который подвергается узкополосной доплеровской фильтрации на основе процедуры БПФ и преобразуется в амплитудно-частотный спектр S(f) (рисунок 2б), составляющие которого обусловлены отражениями сигнала от планеров самолетов группы и вращающихся лопаток рабочего колеса КНД их силовых установок.
В формирователе 2 (рисунок 1) отсчетов доплеровских частот планерных составляющих спектра сигнала, во-первых, определяется отсчет доплеровской частоты, соответствующий максимальной амплитуде спектральной составляющей спектра сигнала, который обусловлен его отражением от планера одного из самолетов из их группы (рисунок 2б, спектральная составляющая №2), во-вторых, слева и справа в полосе частот ±ΔF относительно частоты спектральной составляющей спектра сигнала, имеющей максимальную амплитуду, определяются частоты локальных максимумов спектра сигнала и их количество (рисунок 2б, Nлм=2), которые превысили установленный порог.
В формирователе 3 (рисунок 1) отсчета доплеровской частоты центроида первых компрессорных составляющих спектра сигнала, во-первых, определяются отсчеты доплеровских частот локальных максимумов, имеющие наибольшие амплитуды, количество которых равно количеству локальных максимумов, расположенных в полосе доплеровских частот ±ΔF относительно спектральной составляющей с максимальной амплитудой, соответствующей отражению сигнала от планера одного из самолетов группы, и находящихся справа на частотах, превышающих значение ΔF, относительно частоты спектральной составляющей, имеющей максимальную амплитуду, которые обусловлены отражениями сигнала от лопаток рабочих колеса первых ступеней КНД самолетов группы (рисунок 2б, локальные максимумы №1 и 2, находятся справа вне полосы частот ±ΔF), во-вторых, по выявленным отсчетам доплеровских частот локальных максимумов, находящихся справа на частотах, превышающих значение ΔF и соответствующих отражениями сигнала от лопаток рабочих колеса первых ступеней КНД двигателей самолетов группы, вычисляется отсчет доплеровской частоты центроида, как среднее значение отсчетов доплеровских частот локальных максимумов.
В результате на выходе формирователя 2 доплеровских отсчетов (рисунок 1) формируется вектор наблюдения Y(k) отсчетов доплеровских частот, обусловленных отражениями сигнала от планеров самолетов группы, которые поступают на вход оптимального фильтра 4 сопровождения групповой воздушной цели, работающего в соответствии с процедурой многомерной линейной дискретной калмановской фильтрации (1)-(6), а на выходе формирователя 3 отсчета доплеровской частоты центроида первых компрессорных составляющих спектра сигнала - вектор наблюдения Yц(k) отсчета центроида доплеровской частоты, который поступает на вход оптимального фильтра 5 сопровождения центроида первых компрессорных составляющей спектра сигнала, работающего аналогично, как и оптимальный фильтр 4, в соответствии с процедурой (1)-(6). При этом, размерность матриц, входящих в процедуру (1)-(6) для оптимального фильтра 4, определяется количеством локальных максимумов Nлм спектральных составляющих сигнала, находящихся в полосе частот ±ΔF.
На каждом k-ом такте работы обоих оптимальных фильтров 4 и 5 сопровождения в блоке 6 вычитания, состоящим из отдельных устройств вычитания, вычисляются оценки разностей между оцененными значениями доплеровских частот, соответствующих центроиду доплеровской частоты , обусловленной отражениями сигнала от лопаток рабочего колеса первых ступеней КНД двигателей самолетов группы, и отражениям сигнала от планера каждого j-го самолета группы, . Причем, данные оценки разностей при отсутствии воздействия уводящих по скорости помех являются постоянными величинами для каждого типа самолета в группе из класса «самолет с турбореактивным двигателем».
В блоке 7 дифференцирования, состоящим их отдельных устройств дифференцирования, вычисляются модули производных оценок разностей между оцененными значениями доплеровских частот , обусловленными отражениями сигнала от планера каждого j-го самолета группы и центроидом доплеровской частоты , которые поступают на соответствующие пороговые устройства блока 8 порогов, куда также, на все пороговые устройства поступает пороговое значение ε, близкое по величине к нулю.
При выполнении условия для всех j-ых, , оценок планерных составляющих спектра сигнал, что соответствует отсутствию воздействия уводящих по скорости помех, (поскольку (рисунки 2 в, г, д - временные участки [t0;t1] и [t2;t3]) величины - постоянны и их производные равны нулю), на выходе блока 8 порогов формируется цифровой код (на рисунке 1-«p(j=m)/з(j=m)»), который является разрешающим (индекс «p(j-m)») кодом для коммутаторов 9 и 10, запрещающим (индекс «з(j=m)») кодом для оперативного запоминающего устройства (ОЗУ) 11 и вычислителя 12 параметров автокорреляционных функций (АКФ) оценок доплеровских частот планерных составляющих спектра сигнала. В этом случае в качестве выходной информации используются все оценки , , J=M, доплеровских частот, вычисляемые в соответствии с процедурой (1)-(6) на основе наблюдения Y(k), которые также через первый коммутатор 9 (на его втором входе присутствует разрешающий код «p(j=m)»), поступают на вход ОЗУ 11 (построен на основе сдвиговых регистров), где осуществляется хранение текущих оценок доплеровских частот при отсутствии воздействия уводящей по скорости помехи. Поскольку для ОЗУ 11 в данном случае будет сформирован запрещающий код «зj=m)» на выходе блока 8 порогов, то хранящиеся в нем текущие оценки далее не поступают на вход вычислителя 12 параметров АКФ, на втором входе которого также будет присутствовать запрещающий код «з(j=m)» с выхода блока 8 порогов.
При не выполнении условия для всех или части оценок доплеровских частот, что соответствует воздействию уводящих по скорости помех, поскольку (рисунки 2 в, г, д - временной участок [t1;t2]) величины в этом случае будут постоянными, если изменение уводящей по скорости помехи осуществляется по линейному закону или изменяться по другому закону в соответствии с законом изменения уводящей по скорости помехи. В обоих случаях значения величин будут существенно отличаться от нуля, и на выходе блока 8 порогов сформируется результирующий код, который будет состоять из двух частей - запрещающего кода «з(n)» и разрешающего кода «р(m)» для коммутаторов 9 и 10, из разрешающего кода «p(n)» и запрещающего кода «з(n)» для ОЗУ 11 и вычислителя 12 параметров АКФ.
В результате на вход ОЗУ 11 продолжают поступать оценки доплеровских частот , , которые не были идентифицированы, как уводящие по скорости помехи (разрешающий код «р(m)» на входе первого коммутатора 9) и прекращается поступление текущих оценок , , которые были идентифицированы, как уводящие по скорости помехи (запрещающий код «з(n)» на входе первого коммутатора 9). Только те оценки доплеровских частот , которые изначально не были (рисунки 2 в, г, д - временные участки [t0;t1] и [t2;t3]), а затем были идентифицированы (временной участок [t1;t2]), как уводящие по скорости помехи, с выхода ОЗУ 11 (на его входе присутствует разрешающая «p(n)» и запрещающая «з(m)» составляющие кода) поступают на вход вычислителя 12 параметров АКФ (на его входе присутствует также разрешающая «p(n)» и запрещающая «з(m)» составляющие кода), в котором вычисляются АКФ оценок доплеровских частот в дискретном времени (которые были получены до их идентификации, как уводящие по скорости помехи) в соответствии с выражением [3]
где
Δt - интервал дискретизации;
z=1, 2, ..., p.
mF - тренд (изменяющееся в дискретном времени математическое ожидание) оценки доплеровской частоты.
Вычисленная в соответствии с выражением 9 автокорреляционная функция для каждой оценки доплеровской частоты, аппроксимируется спадающей по экспоненциальному закону косинусоидальной зависимостью вида [3]
где
σv, τv и fv - соответственно среднеквадратическое отклонение, время корреляции и собственная частота, которые являются параметрами АКФ.
Численные значения параметров АКФ (σv, τv, fv) каждой оценки , , доплеровской частоты, которая была получена до ее идентификации, как уводящей по скорости помехи, поступают на вход вычислителя 13 оценок доплеровских частот, которые вычисляются в соответствии с выражением (9).
При этом, в вычислителе 13, во-первых, на первом шаге рекуррентной процедуры вычисления оценок принимаются конечные значения оценок , , доплеровских частот до их идентификации, как уводящих по скорости помех, поступившие и хранящиеся в ОЗУ 11.
Во-вторых, матрицы , , размерностью 3×3, будут иметь следующие, отличные от нуля, элементы: ϕ11=1; ϕ12=Δt; ϕ21=-βΔt; ϕ22=1-αΔt; ϕ33=1+Δt, которые соответствуют динамической модели радиальных составляющих фазовых координат полета каждого самолета группы относительно носителя РЛС, описываемой следующей системой дифференциальных уравнений в непрерывном времени [3]
где
β=(2πfv)2 - квадрат частоты fv собственных колебаний АКФ, которые обусловлены скоростными флюктуациями полета каждого самолета группы и носителя РЛС;
n(t) - формирующий белый гауссовский шум с нулевыми средними значениями и единичными интенсивностями;
Vo - постоянная составляющая скорости полета каждого самолета группы;
Матрица наблюдения H(k+1), размерностью 1×3, будет иметь следующие, отличные от нуля, элементы: h11=h13=1.
В результате при воздействии уводящей по скорости помехи (рисунок 2в, временной участок [t1,t2]) на выходе вычислителя 13 будут формироваться оценки , , доплеровских частот, которые поступают на вход второго коммутатора 10. В результате на его выходе будут сформированы оценки доплеровских частот , на основе наблюдения Y(k) в соответствии с процедурой фильтрации (1)-(6) и которые не идентифицированы, как уводящие по скорости помехи (рисунок 2в, временные участки [t0;t1] и [t2;t3]), а также оценки доплеровских частот , , только на основе переходных матриц,, , которые были идентифицированы, как уводящие по скорости помехи (рисунок 2в, временной участок [t1;t2]). При этом общее количество оценок будет равно J=M+N.
Для оценки работоспособности предлагаемого способа было проведено его имитационное моделирование. При этом, были использованы радиолокационные сигналы, отраженные от группы из 4-х самолетов из класса «самолеты с турбореактивными двигателями», которые при проведении летно-экспериментальных исследований были зарегистрированы на промежуточной частоте с выхода линейной части приемника бортовой РЛС с фазированной антенной решеткой, построенной по импульсно-доплеровскому принципу обработки сигналов в сантиметровом диапазоне волн.
При узкополосном спектральном анализе зарегистрированных реальных радиолокационных сигналов в процедуре БПФ эквивалентная полоса пропускания его одного бина принималась равной 10 Гц.
Численные значения параметров динамических моделей, входящих в оптимальные фильтры 4 и 5, а также в вычислителе 13, принимались из примера, приведенного в [3].
В интервале времени (рисунок 2в) [t1;t2] имитировалась уводящая по доплеровской частоте помеха со скоростью 575 Гц/с.
Пороговое значение величины для всех пороговых устройств блока 8 порогов составляло ε=0,01.
В результате имитационного моделирования по реальным сигналам получены следующие обобщенные характеристики при отношениях сигнал/шум 14-24 дБ:
среднеквадратическая ошибка оценки доплеровской частоты:
без воздействия уводящей по скорости (доплеровской частоты)
помехи - 0,9-2,2 Гц;
при воздействии уводящей по скорости (доплеровской частоты)
помехи - 1,6-3,7 Гц,
что является приемлемым на практике.
Таким образом, применение предлагаемого изобретения позволит идентифицировать воздействие или отсутствие уводящих по скорости помех при сопровождении целей в РЛС и, в зависимости от результата идентификации, сформировать достоверные оценки доплеровских частот, обусловленных скоростью сближения каждого самолета группы с носителем РЛС.
Источники информации
1. Фарина А., Студер Ф. Цифровая обработка радиолокационной информации. Сопровождение целей. / Пер. с англ. - М.: Радио и связь, 1993, с. 246-248 (аналог).
2. Способ сопровождения групповой воздушной цели из класса «самолеты с турбореактивными двигателями». Патент на изобретение №2456633, 2011 (прототип).
3. Богданов А.В., Васильев О.В., Голубенко В.А., Маняшин С.М., Филонов А.А. Методика построения динамических моделей радиальных скоростей и ускорений пары воздушных целей, летящих в сомкнутом боевом порядке // Теория и системы управления, 2007 - №4 (страницы 139, 142, 145, 146 - формулы (2.1), (2.2), (3.2)-(3.5), (5.3), (5.4) и (5.11), пример п. 2, 3, 4, 7).
Claims (28)
- Способ сопровождения в радиолокационной станции групповой воздушной цели из класса «самолеты с турбореактивными двигателями» при воздействии уводящих по скорости помех, заключающийся в том, что сигнал, отраженный от групповой воздушной цели, на промежуточной частоте с выхода приемника радиолокационной станции подвергается узкополосной доплеровской фильтрации на основе процедуры быстрого преобразования Фурье и преобразуется в амплитудно-частотный спектр, составляющие которого обусловлены отражениями сигнала от планеров самолетов группы и вращающихся лопаток рабочего колеса компрессора низкого давления их силовых установок, определяется отсчет доплеровской частоты, соответствующий максимальной амплитуде спектральной составляющей спектра сигнала, который обусловлен его отражением от планера одного из самолетов группы, слева и справа в полосе частот ±ΔF, где ΔF - априорно заданная полоса доплеровских частот, занимаемая спектральными составляющими, обусловленными отражениями сигнала от планеров самолетов группы, относительно частоты спектральной составляющей спектра сигнала, имеющей максимальную амплитуду, определяются отсчеты доплеровских частот локальных максимумов спектра сигнала и их количество, которые превысили установленный порог, определяются отсчеты доплеровских частот локальных максимумов, имеющие наибольшие амплитуды, количество которых равно количеству локальных максимумов, расположенных в полосе доплеровских частот ±ΔF относительно спектральной составляющей с максимальной амплитудой, соответствующей отражению сигнала от планера одного из самолетов группы, и находящихся справа на частотах, превышающих значение ΔF, относительно частоты спектральной составляющей, имеющей максимальную амплитуду, которые обусловлены отражениями сигнала от лопаток рабочего колеса первых ступеней компрессора низкого давления двигателя самолетов группы, отличающийся тем, что вычисляется процедура оптимальной многомерной линейной дискретной калмановской фильтрации в соответствии с выражениями
- где
- К - общее количество тактов работы фильтра;
- P-(k+1) и P(k+1) - ковариационные матрицы ошибок экстраполяции и фильтрации соответственно;
- Ф(k) - переходная матрица состояния;
- Q(k+1) и R(k+1) - ковариационные матрицы шумов возбуждения и наблюдения соответственно;
- S(k+1) - матрица весовых коэффициентов;
- I - единичная матрица;
- J - количество самолетов в группе;
- Н(k) - матрица наблюдения;
- Y(k) - вектор наблюдения отсчетов доплеровских частот;
- Z(k+1) - матрица невязок измерения;
- Ψ(k+1) - матрица априорных ошибок фильтрации;
- "-1" - операция вычисления обратной матрицы;
- "т" - операция транспонирования матрицы,
- при сопровождении по доплеровской частоте каждой j-й, , планерной составляющей спектра сигнала, обусловленной скоростью сближения каждого j-го, , самолета группы с носителем радиолокационной станции, при этом размерность матриц, входящих в процедуру (1)-(6), определяется количеством локальных максимумов спектральных составляющих сигнала, находящихся в полосе частот ±ΔF относительно спектральной составляющей с максимальной амплитудой, соответствующей отражению сигнала от планера одного из самолетов группы, по выявленным отсчетам доплеровских частот локальных максимумов, находящихся справа на частотах, превышающих значение ΔF и соответствующих отражениями сигнала от лопаток рабочего колеса первых ступеней компрессора низкого давления двигателей самолетов группы, вычисляется отсчет доплеровской частоты центроида как среднее значение отсчетов доплеровских частот локальных максимумов, который поступает на вход оптимального фильтра сопровождения центроида отсчетов доплеровских частот, обусловленных отражением сигнала от лопаток рабочего колеса первых ступеней компрессора низкого давления двигателей самолетов группы, работающего аналогично, как и оптимальный фильтр сопровождения планерных составляющих спектра сигнала в соответствии с процедурой (1)-(6), на каждом k-м такте работы обоих оптимальных фильтров сопровождения вычисляются оценки разностей между оцененными значениями доплеровских частот, соответствующих центроиду доплеровской частоты , обусловленной отражениями сигнала от лопаток рабочего колеса первых ступеней компрессора низкого давления двигателей самолетов группы, и отражениям сигнала от планера каждого j-го самолета группы, , вычисляются модули производных оценок разностей между оцененными значениями доплеровских частот , обусловленными отражениями сигнала от планера каждого j-го самолета группы и центроидом доплеровской частоты , обусловленной отражениями сигнала от лопаток рабочего колеса первых ступеней компрессора низкого давления двигателей самолетов группы, которые сравниваются с пороговым значением ε, близким к нулю, при выполнении условия для всех j-х, , оценок планерных составляющих спектра сигнал, что соответствует отсутствию воздействия уводящих по скорости помех, в качестве выходной информации используются все оценки , , доплеровских частот, вычисляемые в соответствии с процедурой (1)-(6) на основе наблюдения Y(k), при невыполнении условия для всех , что соответствует воздействию уводящих по скорости помех, на выходе, наряду с оценками доплеровских частот , где ; М - количество оценок доплеровских частот, которые не идентифицированы как уводящие по скорости помехи, формируются оценки доплеровских частот , где ; N - количество оценок доплеровских частот, которые идентифицированы как уводящие по скорости помехи; J=M+N, вычисляемые в соответствии с процедурой
- где
- - оцененный вектор состояния, содержащий оценку доплеровской частоты , , полученную после ее идентификации как уводящая по скорости помеха, только на основе переходной матрицы состояния , учитывающей модель взаимного перемещения носителя РЛС и того самолета группы, отраженный от которого сигнал изначально еще не был идентифицирован как уводящая по скорости помеха.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2016108454A RU2617110C1 (ru) | 2016-03-09 | 2016-03-09 | Способ сопровождения в радиолокационной станции групповой воздушной цели из класса "самолёты с турбореактивными двигателями" при воздействии уводящих по скорости помех |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2016108454A RU2617110C1 (ru) | 2016-03-09 | 2016-03-09 | Способ сопровождения в радиолокационной станции групповой воздушной цели из класса "самолёты с турбореактивными двигателями" при воздействии уводящих по скорости помех |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| RU2617110C1 true RU2617110C1 (ru) | 2017-04-21 |
Family
ID=58643038
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2016108454A RU2617110C1 (ru) | 2016-03-09 | 2016-03-09 | Способ сопровождения в радиолокационной станции групповой воздушной цели из класса "самолёты с турбореактивными двигателями" при воздействии уводящих по скорости помех |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2617110C1 (ru) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2665031C1 (ru) * | 2018-01-31 | 2018-08-27 | Федеральное государственное казённое военное образовательное учреждение высшего образования "Военная академия воздушно-космической обороны имени Маршала Советского Союза Г.К. Жукова" Министерства обороны Российской Федерации | Способ сопровождения воздушной цели из класса "самолёт с турбореактивным двигателем" при воздействии уводящих по дальности и скорости помех |
| RU2705070C1 (ru) * | 2019-04-16 | 2019-11-05 | Федеральное государственное казённое военное образовательное учреждение высшего образования "Военная академия воздушно-космической обороны имени Маршала Советского Союза Г.К. Жукова" Министерства обороны Российской Федерации | Способ распознавания типа самолёта с турбореактивным двигателем в импульсно-доплеровской радиолокационной станции |
| RU2713635C1 (ru) * | 2019-05-27 | 2020-02-05 | Федеральное государственное унитарное предприятие "Государственный научно-исследовательский институт авиационных систем" (ФГУП "ГосНИИАС") | Способ сопровождения в радиолокационной станции воздушной цели из класса "самолёт с турбореактивным двигателем" при воздействии уводящих по дальности и скорости помех |
| RU2726869C1 (ru) * | 2019-08-01 | 2020-07-16 | Федеральное государственное унитарное предприятие "Государственный научно-исследовательский институт авиационных систем" (ФГУП "ГосНИИАС") | Способ распознавания функционального назначения летательных аппаратов пары по принципу "ведущий-ведомый" |
| RU2786518C1 (ru) * | 2022-02-03 | 2022-12-21 | Федеральное государственное казённое военное образовательное учреждение высшего образования "Военная академия воздушно-космической обороны имени Маршала Советского Союза Г.К. Жукова" Министерства обороны Российской Федерации | Способ распознавания типового состава групповой воздушной цели из класса "самолеты с турбореактивными двигателями" на основе калмановской фильтрации и нейронной сети |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5402131A (en) * | 1993-09-28 | 1995-03-28 | The United States Of America As Represented By The Secretary Of The Navy | Detection of radar targets using higher-order statistics |
| US5506817A (en) * | 1995-05-25 | 1996-04-09 | The United States Of America As Represented By The Secretary Of The Navy | Enhanced adaptive statistical filter providing sparse data stochastic mensuration for residual errors to improve performance for target motion analysis noise discrimination |
| WO1998036289A1 (de) * | 1997-02-14 | 1998-08-20 | Daimler-Benz Aerospace Ag | Verfahren zur zielklassifizierung |
| RU2324952C1 (ru) * | 2006-08-16 | 2008-05-20 | Общество с ограниченной ответственностью "ОКБ Траверз" | Способ сопровождения групповой воздушной цели |
| RU2392640C1 (ru) * | 2009-07-28 | 2010-06-20 | Дмитрий Геннадьевич Митрофанов | Способ выявления параметров траекторных нестабильностей малоразмерного воздушного объекта в виде радиального ускорения движения для режима сопровождения с использованием сигналов с поимпульсной перестройкой несущей частоты |
| WO2010121118A1 (en) * | 2009-04-17 | 2010-10-21 | Raytheon Company | Methods and apparatus for integration of distributed sensors and airport surveillance radar to mitigate blind spots |
| RU2456633C1 (ru) * | 2011-05-03 | 2012-07-20 | Федеральное государственное научное учреждение "Государственный научно-технологический центр "Наука" (ФГНУ "ГНТЦ "Наука") | Способ сопровождения групповой воздушной цели из класса "самолеты с турбореактивными двигателями" |
-
2016
- 2016-03-09 RU RU2016108454A patent/RU2617110C1/ru not_active IP Right Cessation
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5402131A (en) * | 1993-09-28 | 1995-03-28 | The United States Of America As Represented By The Secretary Of The Navy | Detection of radar targets using higher-order statistics |
| US5506817A (en) * | 1995-05-25 | 1996-04-09 | The United States Of America As Represented By The Secretary Of The Navy | Enhanced adaptive statistical filter providing sparse data stochastic mensuration for residual errors to improve performance for target motion analysis noise discrimination |
| WO1998036289A1 (de) * | 1997-02-14 | 1998-08-20 | Daimler-Benz Aerospace Ag | Verfahren zur zielklassifizierung |
| RU2324952C1 (ru) * | 2006-08-16 | 2008-05-20 | Общество с ограниченной ответственностью "ОКБ Траверз" | Способ сопровождения групповой воздушной цели |
| WO2010121118A1 (en) * | 2009-04-17 | 2010-10-21 | Raytheon Company | Methods and apparatus for integration of distributed sensors and airport surveillance radar to mitigate blind spots |
| RU2392640C1 (ru) * | 2009-07-28 | 2010-06-20 | Дмитрий Геннадьевич Митрофанов | Способ выявления параметров траекторных нестабильностей малоразмерного воздушного объекта в виде радиального ускорения движения для режима сопровождения с использованием сигналов с поимпульсной перестройкой несущей частоты |
| RU2456633C1 (ru) * | 2011-05-03 | 2012-07-20 | Федеральное государственное научное учреждение "Государственный научно-технологический центр "Наука" (ФГНУ "ГНТЦ "Наука") | Способ сопровождения групповой воздушной цели из класса "самолеты с турбореактивными двигателями" |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2665031C1 (ru) * | 2018-01-31 | 2018-08-27 | Федеральное государственное казённое военное образовательное учреждение высшего образования "Военная академия воздушно-космической обороны имени Маршала Советского Союза Г.К. Жукова" Министерства обороны Российской Федерации | Способ сопровождения воздушной цели из класса "самолёт с турбореактивным двигателем" при воздействии уводящих по дальности и скорости помех |
| RU2705070C1 (ru) * | 2019-04-16 | 2019-11-05 | Федеральное государственное казённое военное образовательное учреждение высшего образования "Военная академия воздушно-космической обороны имени Маршала Советского Союза Г.К. Жукова" Министерства обороны Российской Федерации | Способ распознавания типа самолёта с турбореактивным двигателем в импульсно-доплеровской радиолокационной станции |
| RU2713635C1 (ru) * | 2019-05-27 | 2020-02-05 | Федеральное государственное унитарное предприятие "Государственный научно-исследовательский институт авиационных систем" (ФГУП "ГосНИИАС") | Способ сопровождения в радиолокационной станции воздушной цели из класса "самолёт с турбореактивным двигателем" при воздействии уводящих по дальности и скорости помех |
| RU2726869C1 (ru) * | 2019-08-01 | 2020-07-16 | Федеральное государственное унитарное предприятие "Государственный научно-исследовательский институт авиационных систем" (ФГУП "ГосНИИАС") | Способ распознавания функционального назначения летательных аппаратов пары по принципу "ведущий-ведомый" |
| RU2786518C1 (ru) * | 2022-02-03 | 2022-12-21 | Федеральное государственное казённое военное образовательное учреждение высшего образования "Военная академия воздушно-космической обороны имени Маршала Советского Союза Г.К. Жукова" Министерства обороны Российской Федерации | Способ распознавания типового состава групповой воздушной цели из класса "самолеты с турбореактивными двигателями" на основе калмановской фильтрации и нейронной сети |
| RU2802653C1 (ru) * | 2022-10-04 | 2023-08-30 | Федеральное государственное казённое военное образовательное учреждение высшего образования "Военная академия воздушно-космической обороны имени Маршала Советского Союза Г.К. Жукова" Министерства обороны Российской Федерации | Способ распознавания типового состава групповой воздушной цели различных классов при различных условиях ее полета на основе калмановской фильтрации и нейронной сети |
| RU2842696C1 (ru) * | 2024-11-27 | 2025-07-01 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Ярославское высшее военное училище противовоздушной обороны" Министерства обороны Российской Федерации (RU) | Радиолокационная станция сопровождения воздушных целей с защитой от имитирующих уводящих по скорости помех |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2579353C1 (ru) | Способ сопровождения воздушной цели из класса "самолёт с турбореактивным двигателем" при воздействии уводящей по скорости помехи | |
| CN109975780B (zh) | 基于脉冲多普勒雷达时域回波的直升机型号识别算法 | |
| RU2419815C1 (ru) | Способ сопровождения воздушной цели из класса "самолет с турбореактивным двигателем" | |
| RU2665031C1 (ru) | Способ сопровождения воздушной цели из класса "самолёт с турбореактивным двигателем" при воздействии уводящих по дальности и скорости помех | |
| US6573861B1 (en) | Target classification method | |
| RU2456633C1 (ru) | Способ сопровождения групповой воздушной цели из класса "самолеты с турбореактивными двигателями" | |
| RU2617110C1 (ru) | Способ сопровождения в радиолокационной станции групповой воздушной цели из класса "самолёты с турбореактивными двигателями" при воздействии уводящих по скорости помех | |
| CN112859014B (zh) | 基于雷达信号分选的雷达干扰抑制方法、装置及介质 | |
| CN103176187B (zh) | 一种机载预警雷达地面高速公路目标过滤方法 | |
| CN112162273B (zh) | 一种基于奇异矢量的多旋翼无人机物理参数提取方法 | |
| CN103885043B (zh) | 基于广义匹配滤波的飞机目标杂噪稳健分类方法 | |
| CN111722188A (zh) | 基于stft预分选的pri变换雷达信号分选方法 | |
| CN108344982A (zh) | 基于长时间相参积累的小型无人机目标雷达检测方法 | |
| CN109870693A (zh) | 一种风力涡轮机杂波检测方法 | |
| CN117250594B (zh) | 一种雷达目标分类识别方法 | |
| RU2705070C1 (ru) | Способ распознавания типа самолёта с турбореактивным двигателем в импульсно-доплеровской радиолокационной станции | |
| RU2408031C2 (ru) | Способ сопровождения пилотируемой воздушной цели | |
| Molchanov et al. | Aerial target classification by micro-Doppler signatures and bicoherence-based features | |
| Ren et al. | Estimating physical parameters from multi-rotor drone spectrograms | |
| RU2735314C1 (ru) | Способ распознавания типа самолёта с турбореактивным двигателем в импульсно-доплеровской радиолокационной станции при воздействии имитирующих помех | |
| RU2732281C1 (ru) | Способ распознавания типа самолёта с турбореактивным двигателем в импульсно-доплеровской радиолокационной станции при воздействии уводящей по скорости помехи | |
| RU144505U1 (ru) | Устройство сопровождения воздушной цели из класса "самолёт с турбореактивным двигателем" | |
| CN109407069B (zh) | 一种基于选通加权频率分集的弹载雷达角闪烁频率抑制方法 | |
| RU2786518C1 (ru) | Способ распознавания типового состава групповой воздушной цели из класса "самолеты с турбореактивными двигателями" на основе калмановской фильтрации и нейронной сети | |
| CN106483513A (zh) | 一种飞机类目标微多普勒纹理特征提取方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20180310 |




