Способ распознавания варианта наведения подвижного объекта на один из летательных аппаратов группы



Способ распознавания варианта наведения подвижного объекта на один из летательных аппаратов группы
Способ распознавания варианта наведения подвижного объекта на один из летательных аппаратов группы
Способ распознавания варианта наведения подвижного объекта на один из летательных аппаратов группы
Способ распознавания варианта наведения подвижного объекта на один из летательных аппаратов группы

Владельцы патента RU 2713212:

Федеральное государственное унитарное предприятие «Государственный научно-исследовательский институт авиационных систем» (ФГУП «ГосНИИАС») (RU)

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

 

Предполагаемое изобретение относится к области цифровой обработки радиолокационных сигналов и может быть использовано в радиолокационной станции (РЛС) для формирования достоверных оценок радиальных функционально-связанных координат (ФСК) взаимного перемещения группы летательных аппаратов (ЛА) (ЛА 1 - ведущий, носитель РЛС; ЛА 2, ЛА 3, ..., ЛА j - ведомые, где j - численный состав группы), и подвижного объекта, а также распознавания варианта наведения подвижного объекта на один из ЛА группы.

Известен способ распознавания направления самонаведения пущенной по группе самолетов ракеты с радиолокационной головкой самонаведения [1, 2], заключающийся в том, что сигнал, отраженный от подвижного объекта, подвергается узкополосной доплеровской фильтрации на основе процедуры быстрого преобразования Фурье (БПФ) и преобразуется в амплитудно-частотный спектр, определяется отсчет доплеровской частоты, соответствующий максимальной амплитуде спектральной составляющей спектра сигнала, который соответствует его отражениям от корпуса подвижного объекта, и поступает на вход фильтра сопровождения подвижного объекта, функционирующего в соответствии с процедурой (1)-(6) многомерной линейной дискретной калмановской фильтрации в соответствии с уравнениями [1, 2]

; (1)

; (2)

S( k+1 )= P ( k+1 ) H T ( k+1 ) Ψ 1 ( k+1 ) ; (3)

; (4)

X ( k+1 )=Ф( k ) X ( k )+S( k+1 )Z( k+1 ) ; (5)

P( k+1 )=[ IS( k+1 )H( k+1 ) ] P ( k+1 ) , (6)

где k = 0,1, …, К, …, - номер такта работы фильтра;

P (k+1) и P(k+1) - ковариационные матрицы (КМ) ошибок экстраполяции и фильтрации соответственно;

Ф(k) - переходная матрица состояния;

Q(k+1) и R(k+1) - КМ шумов возбуждения и наблюдения соответственно;

S(k+1) - матрица весовых коэффициентов;

I - единичная матрица;

X ^ (k) и X ^ (k+1) - вектор текущих и экстраполированных оценок радиальных ФСК взаимного перемещения группы ЛА и подвижного объекта;

Н(k+1) - матрица наблюдения;

Y(k) - вектор наблюдения;

Z(k+1) - матрица невязок измерения;

Ψ(k+1) - матрица априорных ошибок фильтрации;

"-1" - операция вычисления обратной матрицы;

"т" - операция транспонирования матрицы,

определяется оценка дальности D ^ (k) до подвижного объекта, определяется оценка угловой скорости ω ^ (k) вращения линии визирования «подвижный объект - ЛА 1» и определяется оценка скорости V ^ сбл (k) сближения ЛА 1 с подвижным объектом, на основе математической модели (ММ) полета подвижного объекта в непрерывном времени

D ˙ (t)= V сбл (t)=( V 0зур (t)+Δ V зур (t)+ V 0н (t)+Δ V н (t)) , D(0)= D 0 ; (7)

V ˙ 0зур (t)= a 0зур (t) , V 0зур (0)= V 00зур ; (8)

a ˙ 0зур (t)=0 , a 0зур (0)= a 00зур ; (9)

Δ V ˙ зур (t)= a зур (t) , Δ V зур (0)=Δ V 0зур ; (10)

a ˙ зур (t)= α Vзур a зур (t) β зур Δ V зур (t)+ n зур (t) 2 α Vзур σ aзур 2 , (11)

a зур (0)= a 0зур ; (12)

V ˙ 0н (t)=0 , V 0н (0)= V 00н ; (13)

Δ V ˙ н (t)= a н (t) , Δ V н (0)=Δ V 0н ; (14)

a ˙ н (t)= α Vн a н (t) β н Δ V н (t)+ n н (t) 2 α Vн σ aн 2 , a н (0)= a 0н ; (15)

ϕ ˙ виз (t)=ω(t) , ϕ виз (0)= ϕ 0ур ; (16)

ω ˙ (t)= 2 V сбл (t) D(t) ω(t)+ ( j зур (t) j н (t)) D(t) , ω(0)= ω 0 ; (17)

j зур (t) dt = α j зур j зур (t)+ n j зур (t) 2 α j зур σ j зур 2 , j зур (0)= j 0ур ; (18)

j н (t) dt = n j н (t) 2 α j н σ j н 2 , j н (0)= j 0н , (19)

где D(t) - расстояние между подвижным объектом и ЛА 1 - носителем РЛС;

V сбл (t) - относительная скорость сближения подвижного объекта противника с ЛА 1;

V 0зур (t) и V 0н (t) - детерминированные составляющие радиальных скоростей соответственно подвижного объекта и ЛА 1;

Δ V зур (t) и Δ V н (t) - флюктуационные составляющие радиальных скоростей соответственно подвижного объекта и ЛА 1;

a 0зур (t) - детерминированная составляющая радиального ускорения подвижного объекта;

a зур (t) и a н (t) - флюктуационные составляющие радиальных ускорений соответственно подвижного объекта и ЛА 1;

ϕ виз (t) - угол визирования подвижного объекта с борта ЛА 1;

ω(t) - угловая скорость вращения линии визирования «ЛА 1 - подвижный объект»;

j зур (t) - нормальное ускорение подвижного объекта в горизонтальной плоскости;

j н (t) - нормальное ускорение ЛА 1 в горизонтальной плоскости;

α Vзур и α Vн - частоты скоростных флюктуаций соответственно подвижного объекта и ЛА 1; α j зур и α j н - параметры характеризующие маневренные свойства соответственно подвижного объекта и ЛА 1;

β зур и β н - квадраты частот собственных колебаний автокорреляционных функций соответственно подвижного объекта и ЛА 1, которые обусловлены скоростными флюктуациями;

σ aзур 2 и σ aн 2 - дисперсии флюктуаций радиальных ускорений соответственно подвижного объекта и ЛА 1;

σ j н 2 и σ j зур 2 - дисперсии флюктуаций нормальных ускорений соответственно подвижного объекта и ЛА 1;

n зур (t) , n н (t), n j зур (t) и n j н (t) - формирующие нормированные белые гауссовские шумы (БГШ) с нулевыми математическими ожиданиями (МО) и единичными интенсивностями;

V 00зур и V 00н , Δ V 0зур и Δ V 0н , j 0ур и j 0н , a 00зур , ϕ 0ур , ω 0 - начальные значения соответственно детерминированных составляющих радиальных скоростей подвижного объекта и ЛА 1, флюктуационных составляющих радиальных скоростей подвижного объекта и ЛА 1, нормальных ускорений подвижного объекта и ЛА 1, флюктуационной составляющей радиального ускорения подвижного объекта, угла визирования подвижного объекта с борта ЛА 1, угловой скорости вращения линии визирования «ЛА 1 - подвижный объект», представляемой в процедуре (1)-(6) матрицами , , размерностями 13×13,

на основе полученных оценок в соответствии с выражениями

h ^ t1г = ( D ^ cos(ϕ)) 2 ω ^ г V ^ сбл cos(ϕ) , h ^ t1в = ( D ^ cos(ε)) 2 ω ^ в V ^ сбл cos(ε) , h ^ t1 = h ^ t1г 2 + h ^ t1в 2 ,

где h ^ t1г и h ^ t1в - оценки проекций текущего промаха подвижного объекта относительно ЛА 1 соответственно на горизонтальную и вертикальную плоскости;

ω ^ г и ω ^ в - оценки проекций угловых скоростей вращения линии визирования «ЛА 1 - подвижный объект» соответственно на горизонтальную и вертикальную плоскости;

ϕ и ε - проекции угла визирования подвижного объекта с борта ЛА 1 соответственно на горизонтальную и вертикальную плоскости,

определяется оценка h ^ t1 текущего промаха подвижного объекта относительно ЛА 1, с учетом информации от аппаратуры межсамолетной навигации, а именно оценок угла визирования α ^ виз на каждый ЛА группы, дальности d ^ до него и превышения Δ h ^ гр (принижения) определяется оценка h ^ t2 текущего промаха подвижного объекта относительно ЛА 2, в соответствии с выражением

Δ h ^ t = h ^ t2 h ^ t1 ,

определяется оценка разности текущих промахов Δ h ^ t , при этом положительное значение этой оценки соответствует наведению подвижного объекта на ЛА 1, а отрицательное - наведению подвижного объекта на ЛА 2.

Недостатками данного способа распознавания направления самонаведения пущенной по группе самолетов ракеты с радиолокационной головкой самонаведения являются:

1. Низкая достоверность оценки радиальных ФСК взаимного перемещения группы ЛА и подвижного объекта, и оценки варианта наведения этого объекта.

2. Неадекватность используемой ММ динамики ФСК фактическим взаимным перемещениям группы ЛА и подвижного объекта (модель не учитывает влияния варианта наведения подвижного объекта на динамику ФСК).

3. Отсутствие возможности комплексирования информации РЛС, измеряющей ФСК, и индикатора (обнаружителя) варианта наведения подвижного объекта, представленного бортовым комплексом обороны (БКО).

4. Допущение об отсутствии перенацеливания подвижного объекта (не учитываются априорные данные о смене варианта наведения подвижного объекта).

5. Отсутствие адаптации фильтра к варианту наведения подвижного объекта на один из ЛА группы.

6. Отсутствие определения варианта наведения подвижного объекта на один из ЛА для произвольного численного состава группы ЛА.

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

Для решения технической задачи в способе распознавания направления самонаведения пущенной по группе самолетов ракеты с радиолокационной головкой самонаведения [2], заключающемся в том, что сигнал, отраженный от подвижного объекта, подвергается узкополосной доплеровской фильтрации на основе БПФ и преобразуется в амплитудно-частотный спектр, определяется отсчет доплеровской частоты, соответствующий максимальной амплитуде спектральной составляющей спектра сигнала, который дополнительно поступает на вход многоканального фильтра совместного сопровождения подвижного объекта и распознавания варианта его наведения, функционирующего в соответствии с процедурой квазиоптимальной совместной фильтрации фазовых координат и распознавания состояния марковской структуры линейной стохастической динамической системы [3]

p ˜ k+1 ( s k+1 )= s k q k ( s k+1 | s k ) p ^ k ( s k ) ; (20)

x ˜ k+1 ( s k+1 )= p ˜ k+1 1 ( s k+1 ) s k q k ( s k+1 | s k ) p ^ k ( s k ) A k x ^ k ( s k ) ; (21)

R ˜ k+1 ( s k+1 )= p ˜ k+1 1 ( s k+1 ) s k q k ( s k+1 | s k ) p ^ k ( s k ) { A k R ^ k ( s k ) A k T +

+ F k G k F k T +[ x ˜ k ( s k+1 ) A k x ^ k ( s k )]×

× [ x ˜ k ( s k+1 ) A k x ^ k ( s k )] T } ; (22)

p ^ k+1 ( s k+1 )= α k+1 ( s k+1 ) s k+1 α k+1 ( s k+1 ) ; (23)

x ^ k+1 ( s k+1 )= x ˜ k+1 ( s k+1 )+ R ˜ k+1 ( s k+1 ) C k+1 T ( s k+1 ) Θ k+1 1 ( s k+1 ) Δ k+1 ( s k+1 ) ; (24)

R ^ k+1 ( s k+1 )= R ˜ k+1 ( s k+1 ) R ˜ k+1 ( s k+1 ) C k+1 T ( s k+1 ) Θ k+1 1 ( s k+1 ) C k+1 ( s k+1 ) R ˜ k+1 ( s k+1 ) ; (25)

α k+1 ( s k+1 )= p ˜ k+1 ( s k+1 ) π k ( r k+1 | r k , s k+1 ) [det Θ k+1 ( s k+1 )] 1/2 exp[ h k+1 ( s k+1 )] ; (26)

Θ k+1 ( s k+1 )= С k+1 ( s k+1 ) R ˜ k+1 ( s k+1 ) С k+1 T ( s k+1 )+ E k+1 ( s k+1 ) Q k+1 E k+1 T ( s k+1 ) ; (27)

h( s k+1 )= 1 2 Δ k+1 T ( s k+1 ) Θ k+1 1 ( s k+1 ) Δ k+1 ( s k+1 ) ; (28)

Δ k+1 ( s k+1 )= z k+1 C k+1 ( s k+1 ) x ˜ ( s k+1 ) ; (29)

s ^ k =arg max s k p ^ k ( s k ) ; (30)

x ^ k = s k x ^ k ( s k ) p ^ ( s k ) ; (31)

R ^ k = s k [ R ^ ( s k )+ x ^ k ( s k ) x ^ k T ( s k )] p ^ ( s k ) x ^ k x ^ k T , (32)

основанной на априорных данных в виде ММ системы «подвижный объект - РЛС - БКО - группа ЛА» со случайной скачкообразной структурой (ССС), включающей модели линейной динамики радиальных ФСК взаимного перемещения группы ЛА и подвижного объекта

x k+1 = A k ( s k ) x k + F k ( s k ) ξ k ; (33)

их измерений в РЛС

z k = С k x k + E k ζ k ; (34)

марковской смены варианта наведения подвижного объекта

q k ( s k+1 | s k ) ; (35)

марковского индикатора варианта наведения подвижного объекта, представленного БКО,

π k+1 ( r k+1 | r k , s k+1 ) ; (36)

неуправляемых случайных возмущений и помех

G k , Q k ; (37)

при начальных условиях

p ^ 0 ( s 0 ), x ^ 0 ( s 0 ), R ^ 0 ( s 0 ) , (38)

где k - дискретный момент времени;

x k - вектор радиальных ФСК взаимного перемещения группы ЛА и подвижного объекта;

s k = 1,j ¯ - вариант наведения подвижного объекта ( s k =1 - объект наводится на ЛА 1 - носитель РЛС; s k =2 - объект наводится на ЛА 2; s k =3 - объект наводится на ЛА 3; s k =j - объект наводится на ЛА j);

z k - вектор измерений РЛС;

r k = 1,j ¯ - выходные показания индикатора варианта наведения подвижного объекта ( r k =1 - объект наводится на ЛА 1 - носитель РЛС; r k =2 - объект наводится на ЛА 2; r k =3 - объект наводится на ЛА 3; r k =j - объект наводится на ЛА j);

q k ( s k+1 | s k ) - условные вероятности перенацеливания подвижного объекта с ЛА s k на ЛА s k+1 ;

π k+1 ( r k+1 | r k , s k+1 ) - условные вероятности показаний бортовым комплексом обороны ЛА 1 наведения подвижного объекта на ЛА r k+1 на ( k+1)-м шаге при условии, что фактически наведение подвижного объекта осуществляется на ЛА s k+1 , а на предыдущем шаги было показано наведение на ЛА r k ;

p ˜ k+1 ( s k+1 ) , x ˜ k+1 ( s k+1 ), R ˜ k+1 ( s k+1 ) и p ^ k+1 ( s k+1 ), x ^ k+1 ( s k+1 ), R ^ k+1 ( s k+1 ) - прогнозируемые на один шаг дискретности вперед и апостериорные соответственно вероятности наведения подвижного объекта на ЛА s k+1 , условные математические ожидания ФСК при фиксированном варианте наведения подвижного объекта, условные КМ ошибок оценивания ФСК при фиксированном варианте наведения подвижного объекта;

s ^ k - квазиоптимальная по критерию максимума апостериорной вероятности оценка варианта наведения подвижного объекта;

x ^ k - апостериорное безусловное математическое ожидание ФСК;

R ^ k - апостериорная безусловная КМ ошибок оценивания ФСК;

G k , Q k - КМ соответственно векторов шумов возбуждения F k ξ k и помех E k ζ k ;

ξ k , ζ k - стандартные дискретные векторные белые шумы;

Θ k ( s k ) - условная КМ измерения при фиксированном варианте наведения подвижного объекта;

A k ( s k ) , F k ( s k ) - известные матрицы детерминированных функций от варианта s k наведения подвижного объекта;

C k , E k - известные матрицы коэффициентов;

Θ k 1 ( s k ) - обратная матрица по отношению к матрице Θ k ( s k );

T - операция транспонирования матрицы;

det Θ k ( s k ) - определитель матрицы Θ k ( s k );

exp[] - экспоненциальная функция,

определяется оценка s ^ k варианта наведения подвижного объекта, определяется оценка x ^ k безусловного математического ожидания ФСК, определяется оценка R ^ k безусловной КМ ошибок оценивания ФСК, на основе ММ (33) динамики радиальных ФСК взаимного перемещения группы ЛА и подвижного объекта, которая в непрерывном времени определяется усовершенствованной моделью (7)-(19), отличающейся зависимостью параметров α Vзур (s), α Vн (s), β зур (s), β н (s) от варианта наведения подвижного объекта s.

Новыми признаками, обладающими существенными отличиями, являются:

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

2. Учет влияния варианта наведения подвижного объекта на динамику ФСК через коэффициенты α Vзур (s), α Vн (s), β зур (s), β н (s) в модели (7)-(19).

3. Учет априорных данных о смене варианта наведения подвижного объекта в виде условных вероятностей переходов (35).

4. Комплексирование в (23) информации РЛС, измеряющей ФСК, и индикатора (обнаружителя) варианта наведения подвижного объекта, представленного БКО с моделью (36).

5. Коррекция оценок (21), (22), (31), (32) ФСК, полученных на основе модели (33) и измерений (34), по оцененным вероятностям (23) наведения подвижного объекта на соответствующий ЛА группы и априорным данным (35) о смене варианта наведения (адаптации фильтра к различным вариантам наведения подвижного объекта).

6. Определение варианта наведения подвижного объекта на один из ЛА для произвольного численного состава группы ЛА.

Данные признаки являются существенными и в известных технических решениях не обнаружены.

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

На фиг. 1. приведена блок-схема, поясняющая реализацию предлагаемого способа распознавания варианта наведения подвижного объекта на один из ЛА группы.

На фиг. 2. приведена структурная схема, поясняющая принцип функционирования многоканального фильтра 9.

Способ распознавания варианта наведения подвижного объекта на один из ЛА группы осуществляется следующим образом.

На вход известного блока 1 БПФ, используемого в [2], на промежуточной частоте с выхода приемника РЛС поступает сигнал S(t), отраженный от подвижного объекта, который подвергается узкополосной доплеровской фильтрации на основе процедуры БПФ и преобразуется в амплитудно-частотный спектр, одна из составляющих которого обусловлена отражениями сигнала от корпуса сопровождаемого подвижного объекта.

В известном формирователе 2 измерения, используемом в [2], во-первых, определяется отсчет доплеровской частоты F п (k+1), соответствующий максимальной амплитуде спектральной составляющей спектра сигнала, который соответствует его отражениям от корпуса подвижного объекта, во-вторых, данный отсчет доплеровской частоты преобразуется в значение скорости, как V п (k+1)=λ F п (k+1)/2 (где λ - рабочая длина волны РЛС), в-третьих, поступающее на вход измерение дальности в непрерывном времени D(t) преобразуется в дискретные отсчеты дальности D(k+1).

В результате на выходе блока 2 формируется измерение z k+1 = (D(k+1), V п (k+1)) T , которое поступает на вход нового многоканального фильтра 9 совместного сопровождения подвижного объекта и распознавания варианта его наведения, функционирующего в соответствии с известной процедурой квазиоптимальной совместной фильтрации фазовых координат и распознавания состояния марковской структуры линейной стохастической динамической системы (20)-(32), работающего на основе априорных данных в виде ММ системы «подвижный объект - РЛС - БКО - группа ЛА» со ССС, включающей (блок 10 памяти бортовой ЦВМ) модель (7)-(19) линейной динамики радиальных ФСК 3 взаимного перемещения группы ЛА и подвижного объекта, представленную матрицами ( A k ( s k ), F k ( s k )), модель измерений ФСК в РЛС 4, представленную матрицами ( С k , E k ), модель смены варианта наведения подвижного объекта 5, представленную переходными вероятностями q k (), модель индикатора варианта наведения подвижного объекта 6, представленную переходными вероятностями π k+1 (), модель неуправляемых случайных возмущений и помех 7 при начальных условиях 8, также поступающих на вход многоканального фильтра 9.

При функционировании многоканального фильтра 9 осуществляются следующие операции.

1. Классификатор экстраполяции 1: прогнозирование (20) вероятностей p ˜ k+1 ( s k+1 ) наведения подвижного объекта на каждый из ЛА группы на один шаг дискретности вперед на основе априорных данных о смене варианта наведения подвижного объекта, представленных соответственно начальными (38) и переходными (35) вероятностями цепи Маркова.

2. Условный экстраполятор 2: прогнозирование (21) на один шаг дискретности вперед условных математических ожиданий x ˜ k+1 ( s k+1 ) ФСК при фиксированном варианте наведения подвижного объекта с учетом найденных вероятностей на основе априорных данных (35) о смене варианта наведения подвижного объекта и альтернативных моделей (33) динамики ФСК взаимного перемещения группы ЛА и подвижного объекта.

3. Условный дисперсиометр экстраполяции 3: прогнозирование (22) на один шаг дискретности вперед условных КМ R ˜ k+1 ( s k+1 ) ошибок оценивания ФСК при фиксированном варианте наведения подвижного объекта с учетом найденных вероятностей (20) и МО (21) на основе априорных данных о смене варианта наведения подвижного объекта и альтернативных моделей динамики ФСК взаимного перемещения группы ЛА и подвижного объекта.

4. Классификатор фильтрации 4: оценка апостериорных вероятностей p ^ k+1 ( s k+1 ) наведения подвижного объекта на каждый из ЛА группы, по степени согласованности (26)-(29) спрогнозированных вероятностей (20), математических ожиданий ФСК (21) и КМ (22) ошибок их оценивания с результатами измерений в (29) и показаниями БКО в (26).

5. Условный фильтр 5: оценка (24) условных апостериорных математических ожиданий x ^ k+1 ( s k+1 ) ФСК взаимного перемещения группы ЛА и подвижного объекта при фиксированном варианте наведения подвижного объекта на основе спрогнозированных математических ожиданий (21) и КМ (22) ошибок прогноза с учетом результатов измерения в (29).

6. Условный дисперсиометр фильтрации 6: оценка (25) условных апостериорных КМ R ^ k+1 ( s k+1 ) ошибок оценивания ФСК при фиксированном варианте наведения подвижного объекта на основе спрогнозированных МО (21) и КМ (22) ошибок прогноза с учетом результатов измерения в (29).

7. Идентификатор 7: распознавание (30) такого s ^ k варианта наведения подвижного объекта, для которого найденная апостериорная вероятность (23) окажется больше.

8. Безусловный фильтр 8: нахождение (31) безусловной оценки x ^ k ФСК на основе апостериорных вероятностей (23) варианта подвижного объекта на каждый из ЛА группы и условных апостериорных оценок (24) ФСК, как безусловного математического ожидания.

9. Безусловный дисперсиометр фильтрации 9: нахождение (32) безусловной КМ R ^ k ошибок оценивания ФСК с учетом найденных апостериорных вероятностей (23) наведения подвижного объекта на каждый из ЛА группы, условных математических ожиданий (24) ФСК, условных КМ (25) ошибок их оценивания и безусловных оценок (31) ФСК.

Сущность распознавания варианта наведения подвижного объекта на один из ЛА группы наиболее полно отражают операции 1, 4, 7.

Сформированные на выходе многоканального фильтра 9 оценки s ^ k варианта наведения подвижного объекта, безусловного математического ожидания x ^ k ФСК взаимного перемещения группы ЛА и подвижного объекта, безусловной КМ R ^ k ошибок оценивания ФСК поступают на выход канала сопровождения подвижного объекта в РЛС.

Результаты сравнительного моделирования предлагаемого способа распознавания варианта самонаведения подвижного объекта на один из ЛА группы на основе нового многоканального фильтра сопровождения подвижного объекта, функционирующего в соответствии с процедурой квазиоптимальной совместной фильтрации фазовых координат и распознавания состояния марковской структуры линейной стохастической динамической системы, и известного способа распознавания направления самонаведения пущенной по группе самолетов ракеты с радиолокационной головкой самонаведения [1, 2] на основе процедуры линейной дискретной калмановской фильтрации свидетельствуют с доверительной вероятностью 0,95 о снижении среднеквадратического отклонения ошибки фильтрации на 16% и о повышении вероятности правильной распознавания варианта наведения подвижного объекта на 10%.

Таким образом, применение предлагаемого способа позволит повысить достоверность распознавания варианта наведения подвижного объекта на один из ЛА группы и оценку радиальных ФСК взаимного перемещения группы ЛА и подвижного объекта путем приближения получаемых оценок к их оптимальным значениям за счет учета влияния направления наведения подвижного объекта на динамику ФСК, комплексирования информации РЛС и БКО, учета априорных данных о смене направления наведения подвижного объекта и адаптации фильтра к этим сменам.

ИСТОЧНИКИ ИНФОРМАЦИИ

1. Обнаружение и распознавание пущенной противником ракеты в бортовой радиолокационной станции истребителя: монография / под. ред. В.О. Червакова. - Красноярск: Сиб. федер. ун-т, 2015.

2. Анцифиров А.А., Богданов А.В., Коротков С.С., Филонов А.А. Способ распознавания направления самонаведения пущенной по группе самолетов ракеты с радиолокационной головкой самонаведения // Патент России №26009530. 2017. Бюл. № 4. (прототип).

3. Бухалев, В.А. Оптимальное сглаживание в системах со случайной скачкообразной структурой / В.А. Бухалев. М.: ФИЗМАТЛИТ, 2013, страницы 117-120.

Способ распознавания варианта наведения подвижного объекта на один из летательных аппаратов группы, заключающийся в том, что сигнал, отражённый от подвижного объекта, подвергается узкополосной доплеровской фильтрации на основе процедуры быстрого преобразования Фурье и преобразуется в амплитудно-частотный спектр, определяется отсчёт доплеровской частоты, соответствующий максимальной амплитуде спектральной составляющей спектра сигнала, который соответствует его отражениям от корпуса подвижного объекта, отличающийся тем, что выделенные отсчеты доплеровских частот поступают на вход многоканального фильтра совместного сопровождения подвижного объекта и распознавания варианта его наведения, функционирующего в соответствии с процедурой квазиоптимальной совместной фильтрации фазовых координат и распознавания состояния марковской структуры линейной стохастической динамической системы, основанной на априорных данных в виде математической модели системы «подвижный объект – радиолокационная станция – бортовой комплекс обороны – группа летательных аппаратов» со случайной скачкообразной структурой, включающей модели линейной динамики радиальных функционально связанных координат взаимного перемещения группы летательных аппаратов и подвижного объекта, их измерений в радиолокационной станции ведущего летательного аппарата, марковской смены варианта наведения подвижного объекта, марковского индикатора варианта наведения подвижного объекта, представленного бортовым комплексом обороны ведущего летательного аппарата, неуправляемых случайных возмущений и помех, при начальных условиях, на основе априорных данных о смене варианта наведения подвижного объекта и альтернативных моделей динамики функционально связанных координат по результатам измерений и показаниям бортового комплекса обороны осуществляется совместное оценивание функционально связанных координат и варианта наведения подвижного объекта на основе метода двухмоментной параметрической аппроксимации в нескольких фильтрах, различающихся гипотезой о варианте наведения подвижного объекта, при этом в каждом фильтре на основе априорных данных о смене варианта наведения подвижного объекта, представленных соответственно начальными и переходными вероятностями цепи Маркова, прогнозируются вероятности наведения подвижного объекта на каждый из летательных аппаратов группы на один шаг дискретности вперед, с учетом найденных вероятностей на основе априорных данных о смене варианта наведения подвижного объекта и альтернативных моделей динамики функционально связанных координат взаимного перемещения группы летательных аппаратов и подвижного объекта прогнозируются на один шаг дискретности вперед условные математические ожидания функционально связанных координат при фиксированном варианте наведения подвижного объекта, с учетом найденных вероятностей и математических ожиданий на основе априорных данных о смене варианта наведения подвижного объекта и альтернативных моделей динамики функционально связанных координат взаимного перемещения группы летательных аппаратов и подвижного объекта прогнозируются на один шаг дискретности вперед условные ковариационные матрицы ошибок оценивания функционально связанных координат при фиксированном варианте наведения подвижного объекта, по степени согласованности спрогнозированных вероятностей, математических ожиданий функционально связанных координат и ковариационных матриц ошибок их оценивания с результатами измерений и показаниями бортового комплекса обороны, представленной функцией правдоподобия, осуществляется оценка апостериорных вероятностей наведения подвижного объекта на каждый из летательных аппаратов группы, на основе спрогнозированных математических ожиданий и ковариационных матриц ошибок прогноза с учетом результатов измерения находятся условные апостериорные математические ожидания функционально связанных координат взаимного перемещения группы летательных аппаратов и подвижного объекта при фиксированном варианте наведения подвижного объекта, на основе спрогнозированных ковариационных матриц ошибок прогноза с учетом результатов измерения находятся условные апостериорные ковариационные матрицы ошибок оценивания функциональносвязанных координат при фиксированном варианте наведения подвижного объекта, распознавание варианта наведения подвижного объекта на один из летательных аппаратов группы осуществляется выбором такого варианта наведения, для которого найденная апостериорная вероятность окажется больше, оценивание радиальных функционально связанных координат взаимного перемещения группы летательных аппаратов и наводящегося подвижного объекта осуществляется на основе апостериорных вероятностей наведения подвижного объекта на каждый из летательных аппаратов группы и условных апостериорных оценок функционально связанных координат, нахождением безусловного математического ожидания, с учетом найденных апостериорных вероятностей наведения подвижного объекта на каждый из летательных аппаратов группы, условных математических ожиданий функционально связанных координат, условных ковариационных матриц ошибок их оценивания и безусловных оценок функционально связанных координат находится безусловная ковариационная матрица ошибок оценивания функционально связанных координат.



 

Похожие патенты:

Изобретение относится к радиолокации и может быть использовано для межпозиционного отождествления результатов измерений и определения координат воздушных целей в многопозиционной радиолокационной системе (МПРЛС) в условиях многоцелевой обстановки.

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

Изобретение относится к радиолокации и может использоваться в радиотехнических системах непрерывного излучения, установленных на подвижных объектах, для получения радиолокационного изображения (РЛИ) в процессе дистанционного зондирования земной (водной) поверхности.

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

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

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

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

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

Изобретение относится к радиолокационным системам и заключается в том, что по принятым от радиолокационного объекта (РЛО) радиосигналам оценивают значения расстояния от летательного аппарата (ЛА) - носителя РЛС до РЛО.

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

Изобретение относится к радиолокационным системам и заключается в том, что по принятым от радиолокационного объекта (РЛО) радиосигналам оценивают значения расстояния от летательного аппарата (ЛА) - носителя РЛС до РЛО.

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

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

Изобретение относится к области радиотехники и может использоваться в системах пассивной радиолокации и радиотехнического наблюдения для однопозиционного определения скоростей, координат и траекторий перемещающихся в пространстве радиоизлучающих объектов (РИО).

Изобретение относится к вычислительной технике. Технический результат – повышение точности классификации объектов на радиолокационном изображении.

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

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

Изобретение относится к области радиолокации с частотной манипуляцией непрерывного излучения (ЧМНИ) радиоволн и может быть использовано для обнаружения движущихся целей, измерения расстояния до объекта локации, скорости и направления движения.

Изобретение относится к аэродромным радиотехническим системам обеспечения захода самолета на посадку. Достигаемый технический результат - регулировка информационного параметра курсо-глиссадных радиомаяков, в частности совмещение с высокой точностью линии курса, формируемой курсовым радиомаяком (КРМ) системы посадки самолетов с направлением оси взлетно-посадочной полосы (ВПП), совмещение двух линий курса с направлением оси ВПП при работе КРМ на частотах в двух диапазонах одновременно, выставление с высокой точностью угла глиссады при неизменных высотах подвеса антенн глиссадного радиомаяка (ГРМ).

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

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

Наверх